GMDR: Versatile Software for Detecting Gene-Gene and Gene-Environ- ment Interactions Underlying Complex Traits

Deciphering the Role of Calcium Signaling Pathway-Associated Single Nucleotide Variants in Susceptibility to Hypertension

BACKGROUND

As a complex disease, hypertension (HTN) is influenced by both genetic and environmental factors and their interaction. The calcium signaling pathway is known to be involved in the regulation of blood pressure, and dysfunction in this pathway may contribute to the development of hypertension. Genome-wide association studies (GWAS) have identified several genes in the calcium signaling pathway associated with susceptibility to HTN, including PLCB1, ATP2B1, and ADRB1. The aim of this study was to investigate the possible association between single nucleotide variants (SNVs) in the calcium signaling pathway and HTN.

METHODS

We genotyped three SNVs: rs1801253 (ADRB1), rs6108168 (PLCB1), and rs17249754 (ATP2B1) in a population of 131 patients with hypertension and 115 healthy controls from Kermanshah province, Iran.

RESULTS

Our results showed a strong and significant association between the G allele and the GG and CG genotypes of the rs1801253 variant in the ADRB1 gene and susceptibility to hypertension. However, no significant association was found for the other two SNVs. In addition, the presence of the GCG haplotype (alleles from left to right: rs1801253, rs6108168, and rs17249754) appeared to confer a protective role against HTN.

CONCLUSION

This study has made a significant contribution towards enhancing the comprehension of hypertension development, as well as the early identification of individuals who are at risk.

Genomics and transcriptomics identify quantitative trait loci affecting growth-related traits in silver pomfret (Pampus argenteus)

Pampus argenteus, a species distributed throughout the Indo-West Pacific, plays a significant role in the yield of aquaculture species. However, cultured P. argenteus has always been characterised by unbalanced growth synchronisation among individuals, slow growth rate, and lack of excellent germplasm resources. Therefore, we conducted mass selection for fast-growing strain P. argenteus for several consecutive years. Various genetic improvement programs have modified its genome sequence through selective pressure, leaving nucleotide signals that can be detected at the genomic level. In the present study, we combined bulked segregant analysis and transcriptome sequencing to identify candidate single nucleotide polymorphisms (SNPs) and key genes for growth-related traits in P. argenteus. A total of 7,280,936 SNPs and 2,212,379 insertions/deletions were identified in the extreme phenotypes of the fast-growing and slow-growing groups. Based on the examination of SNP frequency differences and sliding-window analysis, 42 SNPs were identified as candidate markers. Moreover, 14 of the 42 SNPs linked to growth-related traits were confirmed to be credible SNPs, and eight growth-related genes were screened, namely myb-binding protein 1 A, insulin A/B chains, α-1B adrenoceptor, engulfment and cell motility protein 3, myosin light chain kinase family member 4, insulin receptor located, unconventional myosin-9b, and matrilin-1. An optimal three-factor model (SNP4&SNP12&SNP14) was constructed using the generalized multifactor dimensionality reduction method, and its accuracy was verified as 67.72 %. These results may benefit genetic studies and accelerate genetic improvement of fast-growing strains of P. argenteus.

Interaction Between Hypoxia-Inducible Factor 1-alpha Gene Polymorphism and Helicobacter pylori Infection on Gastric Cancer in a Chinese Tibetan Population

To investigate the impact of four single nucleotide polymorphisms (SNPs) of the HIF1α gene and its interaction with Helicobacter pylori (H. pylori) infection on susceptibility to gastric cancer (GC).Logistic regression was used to test the relationship between four SNPs of HIF1α gene and the susceptibility of GC. A generalized multifactor dimensionality reduction (GMDR) model was used to assess the HIF1α gene-H. pylori infection interaction.Logistic regression analysis indicated that both the rs11549465-CT genotype and the T allele were associated with an increased risk of GC, adjusted OR (95% CI) were 1.63 (1.09-2.20) (CT vs. CC) and 1.70 (1.13-2.36) (T vs. C), respectively. We also found that both the rs11549467-A allele and rs11549467-GA genotype were associated with an increased risk of GC, and adjusted OR (95% CI) were 2.21 (1.61-2.86) (GA vs. GG), 2.13 (1.65-2.65) (A vs. G), respectively. However, no statistically significant impact of rs2057482 or rs1957757 on risk of GC was found. The GMDR model indicated a statistically significant two-dimensional model combination (including rs11549467 and H. pylori infection). The selected model had testing balanced accuracy of 0.60 and the best cross-validation consistencies of 10/10 (p = 0.0107). Compared with H. pylori infection negative participants with rs11549467-GG genotype, H. pylori positive participants with the rs11549467-GA genotype had the highest GC risk, the OR (95% CI) was 3.04 (1.98-4.12).The rs11549467-A allele and rs11549467-GA genotype was associated with increased GC risk. Additionally, the gene-environment interaction between HIF-1α-rs11549467 and H. pylori infection was also correlated with an increased risk of GC.

Gene‒environment interaction effect of hypothalamic‒pituitary‒adrenal axis gene polymorphisms and job stress on the risk of sleep disturbances

Background

Studies have shown that chronic exposure to job stress may increase the risk of sleep disturbances and that hypothalamic‒pituitary‒adrenal (HPA) axis gene polymorphisms may play an important role in the psychopathologic mechanisms of sleep disturbances. However, the interactions among job stress, gene polymorphisms and sleep disturbances have not been examined from the perspective of the HPA axis. This study aimed to know whether job stress is a risk factor for sleep disturbances and to further explore the effect of the HPA axis gene × job stress interaction on sleep disturbances among railway workers.

Methods

In this cross-sectional study, 671 participants (363 males and 308 females) from the China Railway Fuzhou Branch were included. Sleep disturbances were evaluated with the Pittsburgh Sleep Quality Index (PSQI), and job stress was measured with the Effort-Reward Imbalance scale (ERI). Generalized multivariate dimensionality reduction (GMDR) models were used to assess gene‒environment interactions.

Results

We found a significant positive correlation between job stress and sleep disturbances (P < 0.01). The FKBP5 rs1360780-T and rs4713916-A alleles and the CRHR1 rs110402-G allele were associated with increased sleep disturbance risk, with adjusted ORs (95% CIs) of 1.75 [1.38-2.22], 1.68 [1.30-2.18] and 1.43 [1.09-1.87], respectively. However, the FKBP5 rs9470080-T allele was a protective factor against sleep disturbances, with an OR (95% CI) of 0.65 [0.51-0.83]. GMDR analysis indicated that under job stress, individuals with the FKBP5 rs1368780-CT, rs4713916-GG, and rs9470080-CT genotypes and the CRHR1 rs110402-AA genotype had the greatest risk of sleep disturbances.

Conclusions

Individuals carrying risk alleles who experience job stress may be at increased risk of sleep disturbances. These findings may provide new insights into stress-related sleep disturbances in occupational populations.

Impact of interaction between interleukin-6 gene polymorphism and Helicobacter pylori infection on susceptibility to gastric cancer

OBJECTIVE

This study aimed to evaluate the association between four single nucleotide polymorphisms (SNPs) of the interleukin-6 (IL-6) gene and gastric cancer (GC), and impact of interaction between IL-6 SNPs and Helicobacter pylori (H. pylori ) infection on susceptibility to GC.

METHODS

Logistic regression was used to test the relationships between four SNPs of IL-6 gene and GC susceptibility. A generalized multifactor dimensionality reduction (GMDR) model was employed to assess the interaction effect between IL-6 gene and H. pylori infection on GC risk.

RESULTS

Logistic regression analysis indicated that the rs1800795-C allele was associated with increased GC risk, adjusted ORs (95% CI) were 1.80 (1.21-2.41) (CC vs. GG) and 1.68 (1.09-2.30) (C vs. G), respectively. The rs10499563-C allele was associated with decreased risk of GC, and adjusted ORs (95% CI) were 0.62 (0.31-0.93) (TC vs. TT), 0.52 (0.18-0.89) (CC vs. TT) and 0.60 (0.29-0.92) (C vs. T), respectively. GMDR methods found a two-dimensional model combination (including rs1800795 and H. pylori infection) was statistically significant. The selected model had testing balanced accuracy of 59.85% and the best cross-validation consistencies of 10/10 ( P  = 0.0107). Compared with H. pylori -negative subjects with rs1800795- GG genotype, H. pylori -positive participants with GC or CC genotype had the highest risk of GC, the OR (95% CI) was 3.34 (1.78-4.97).

CONCLUSION

The rs1800795-C allele was associated with increased GC risk and the rs10499563-C allele was associated with decreased GC risk. The interaction between rs1800795 and H. pylori infection was also correlated with increased risk of GC.

Associations between Catechol-O-methyltransferase (COMT) polymorphisms and cognitive impairments, psychiatric symptoms and tardive dyskinesia in schizophrenia

INTRODUCTION

Catechol-O-methyltransferase (COMT) is a crucial enzyme involved in dopamine metabolism and has been implicated in the etiology of tardive dyskinesia (TD). We aimed to investigate the associations between COMT gene polymorphisms and the occurrence and severity of TD in a Chinese population, as well as the impact on the psychiatric symptoms and cognitive impairments observed in TD patients.

METHODS

A total of 216 chronic schizophrenia patients, including 59 TD patients and 157 NTD patients, were recruited for this study. Three SNPs of the COMT gene (rs4680, rs165599 and rs4818) were selected and genotyped using matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF-MS). TD severity, 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

In TD patients, total AIMs scores were higher in carriers of the rs4680 AA genotype than in carriers of the AG and GG genotypes (p = 0.01, 0.006), carriers of the rs4818 GC and CC genotypes had higher orofacial scores than in GG genotypes (p = 0.032, 0.002). In male TD patients, carriers of the rs165599 GA genotype scored lower in the extremities and trunk scores than AA genotype carriers (p = 0.015). Moreover, in male TD patients, COMT rs4818 was associated with cognition, since the C allele carriers had significantly higher immediate memory (p = 0.043) and verbal function (p = 0.040) scores than the G allele carriers. In addition, rs165599 genotype interacted with TD diagnosis on depressed factor (p = 0.031).

CONCLUSION

Within the Chinese population, COMT gene polymorphisms could potentially serve as biomarkers for the symptoms and prognosis of TD patients.

MTNR1B gene variations and high pre-pregnancy BMI increase gestational diabetes mellitus risk in Chinese women

AIM

To investigate the association of melatonin receptor 1B (MTNR1B) gene variations and pre-pregnancy body mass index (BMI) with gestational diabetes mellitus (GDM).

MATERIALS AND METHOD

In this study, 1566 Chinese Han pregnant women were enrolled and multiple genetic models were used to evaluate the association between MTNR1B gene polymorphisms and the risk of GDM. The clinical value of pre-pregnancy BMI in predicting GDM was analyzed and evaluated using receiver operating characteristic (ROC) curves. Several methods of analysis were used to examine the impact of gene-gene and gene-BMI interactions on the incidence of GDM influence.

RESULTS

For the MTNR1B gene, rs1387153 (C > T), rs10830962 (C > G), rs4753426 (T > C), and rs10830963 (C > G) are all risk mutations associated with the susceptibility of GDM. The ROC curve analysis indicated that the BMI demonstrated an area under the curve (AUC) of 0.595. Alongside, the sensitivity and specificity stood at 0.676 and 0.474 respectively. The maximum Joden index was found to be 0.150, with a corresponding critical BMI value of 20.5691 kg/m2. Interaction analysis revealed that gene-gene and gene-BMI interactions had no significant effect on GDM occurrence.

CONCLUSION

MTNR1B genetic variations confers the risk to GDM in Chinese women. Furthermore, the high pre-pregnancy BMI (≥20.5691 kg/m2) significantly increases the risk of GDM in Chinese women.

Effects of different iodine levels on the DNA methylation of intrinsic apoptosis-associated genes and analysis of gene-environment interactions in patients with autoimmune thyroiditis

Iodine is an essential nutrient that may change the occurrence of autoimmune thyroiditis (AIT). Apoptosis and DNA methylation participate in the pathogenesis and destructive mechanism of AIT. We detected the methylation and the expression of mRNA of intrinsic apoptosis-associated genes (YWHAG, ING4, BRSK2 and GJA1) to identify the potential interactions between the levels of methylation in these genes and different levels of iodine. 176 adult patients with AIT in Shandong Province, China, were included. The MethylTargetTM assay was used to verify the levels of methylation. We used PCR to detect the mRNA levels of the candidate genes. Interactions between methylation levels of the candidate genes and iodine levels were evaluated with multiplicative and addictive interaction models and GMDR. In the AIT group, YWHAG_1 and six CpG sites and BRSK2_1 and eight CpG sites were hypermethylated, whereas ING4_1 and one CpG site were hypomethylated. A negative correlation was found between methylation levels of YWHAG and mRNA expression. The combination of iodine fortification, YWHAG_1 hypermethylation and BRSK2_1 hypermethylation was significantly associated with elevated AIT risk. A four-locus model (YWHAG_1 × ING4_1 × BRSK2_1 × iodine level) was found to be the best model of the gene-environment interactions. We identified abnormal changes in the methylation status of YWHAG, ING4 and BRSK2 in patients with AIT in different iodine levels. Iodine fortification not only affected the methylation levels of YWHAG and BRSK2 but also interacted with the methylation levels of these genes and may ultimately increase the risk of AIT.

Polymorphisms in gene UGT1A1 modify the association of prenatal exposure to polycyclic aromatic hydrocarbons with congenital heart diseases risk

INTRODUCTION

Prenatal exposure to polycyclic aromatic hydrocarbons (PAHs) is a risk factor for the occurrence of congenital heart diseases (CHDs). Genetic susceptibility to PAHs metabolism may modify the exposure-risk relationship. The role of uridine diphosphoglucuronosyl transferase 1A1 (UGT1A1) genetic polymorphisms for modulating the impacts of prenatal PAHs exposure on the risk of CHDs remains to be discovered.

OBJECTIVE

The aim of this study was to investigate whether maternal UGT1A1 genetic polymorphisms are associated with fetal susceptibility to CHDs and to assess whether the risk is modified by maternal PAHs exposure.

METHODS

Maternal urinary biomarker of PAHs exposure was determined in 357 pregnant women with CHDs fetuses and 270 controls (pregnant women carrying fetuses without major congenital malformations). Urinary 1-hydroxypyrene-glucuronide (1-OHPG) concentration, a sensitive biomarker for PAHs exposure, was measured using ultra-high-performance liquid chromatography coupled with tandem mass spectrometry. Maternal single nucleotide polymorphisms (SNPs) in UGT1A1, including rs3755319, rs887829, rs4148323, rs6742078, and rs6717546, were genotyped using an improved multiplex ligation detection reaction (iMLDR) technique. Unconditional logistic regression was performed to determine the impacts of UGT1A1 polymorphisms on the risks of CHDs and their subtypes. Generalized multifactor dimensionality reduction (GMDR) was used to analyze the gene-gene and gene-PAHs exposure interactions.

RESULTS

None of the selected UGT1A1 polymorphisms was independently associated with the risk of CHDs. The interaction between SNP rs4148323 and PAHs exposure was observed to be associated with CHDs (p< .05). Pregnant women with high-level PAHs exposure and rs4148323 had an increased risk of carrying CHDs fetuses (GA-AA vs. GG: aOR = 2.00, 95% CI = 1.06-3.79). Moreover, the joint effect of rs4148323 and PAHs exposure was found to be significantly associated with risks of septal defects, conotruncal heart defects, and right-sided obstructive malformations.

CONCLUSIONS

Maternal genetic variations of UGT1A1 rs4148323 may modify the association between prenatal PAHs exposure and CHDs risk. This finding needs to be further confirmed in a larger-scale study.

Outdoor time influences VIPR2 polymorphism rs2071623 to regulate axial length in Han Chinese children

Clinical relevance

Identification of individuals with a higher risk of developing refractive error under specific gene and environmental backgrounds, especially myopia, could enable more personalized myopic control advice for patients.

Background

Refractive error is a common disease that affects visual quality and ocular health worldwide. Its mechanisms have not been elaborated, although both genes and the environment are known to contribute to the process. Interactions between genes and the environment have been shown to exert effects on the onset of refractive error, especially myopia. Axial length elongation is the main characteristic of myopia development and could indicate the severity of myopia. Thus, the purpose of the study was to investigate the interaction between environmental factors and genetic markers of VIPR2 and their impact on spherical equivalence and axial length in a population of Han Chinese children.

Methods

A total of 1825 children aged 13~15 years in the Anyang Childhood Eye Study (ACES) were measured for cycloplegic autorefraction, axial length, and height. Saliva DNA was extracted for genotyping three single-nucleotide polymorphisms (SNPs) in the candidate gene (VIPR2). The median outdoor time (2 h/day) was used to categorize children into high and low exposure groups, respectively. Genetic quality control and linear and logistic regressions were performed. Generalized multifactor dimensional reduction (GMDR) was used to investigate gene-environment interactions.

Results

There were 1391 children who passed genetic quality control. Rs2071623 of VIPR2 was associated with axial length (T allele, β=-0.11 se=0.04 p=0.006), while SNP nominally interacted with outdoor time (T allele, β=-0.17 se=0.08 p=0.029). Rs2071623 in children with high outdoor exposure had a significant interaction effect on axial length (p=0.0007, β=-0.19 se=0.056) compared to children with low outdoor exposure. GMDR further suggested the existence of an interaction effect between outdoor time and rs2071623.

Conclusions

Rs2071623 within VIPR2 could interact with outdoor time in Han Chinese children. More outdoor exposure could enhance the protective effect of the T allele on axial elongation.

Gene-environment interaction in myopia

Myopia is a health issue that has attracted global attention due to its high prevalence and vision-threatening complications. It is well known that the onset and progression of myopia are related to both genetic and environmental factors: more than 450 common genetic loci have been found to be associated with myopia, while near work and outdoor time are the main environmental risk factors. As for many complex traits, gene-environment interactions are implicated in myopia development. To date, several genetic loci have been found to interact with near work or educational level. Gene-environment interaction research on myopia could yield models that provide more accurate risk predictions, thus improving targeted treatments and preventive strategies. Additionally, such investigations might have the potential to reveal novel genetic information. In this review, we summarised the findings in this field and proposed some topics for future investigations.

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.

Association of PI3K/AKT/mTOR pathway autophagy-related gene polymorphisms with pulmonary tuberculosis susceptibility in a Chinese population

BACKGROUND

Autophagy can inhibit the survival of intracellular microorganisms including Mycobacterium tuberculosis (Mtb), and the PI3K/AKT/mTOR pathway plays a crucial role. This study investigated the association between PI3K/AKT/mTOR pathway autophagy-related gene polymorphisms and pulmonary tuberculosis (PTB) susceptibility.

METHODS

KEGG pathway and gene ontology (GO) databases were searched for genes belonging to the PI3K/AKT/mTOR and autophagy pathways. Thirty SNPs in nine genes were identified and tested for their associations with tuberculosis in 130 patients with PTB and 271 controls. We constructed genetic risk scores (GRSs) and divided the participants into 3 subgroups based on their GRSs:0-5, 6-10, and 11-16.

RESULTS

This analysis revealed that the AKT1 (rs12432802), RPTOR (rs11654508, rs12602885, rs2090204, rs2589144, and rs2672897), and TSC2 (rs2074969) polymorphisms were significantly associated with PTB risk. A decreasing trend was observed (P trend 0.020), in which a lower GRS was associated with a higher risk of PTB ([6-10] vs. [0-5]: OR (95%CI) 0.590 (0.374-0.931); [11-16] vs. [0-5]: OR (95%CI) 0.381 (0.160-0.906)).

CONCLUSIONS

Polymorphisms in AKT1, RPTOR, and TSC2 may influence susceptibility to PTB.

Relationship between APOE, PER2, PER3 and OX2R Genetic Variants and Neuropsychiatric Symptoms in Patients with Alzheimer's Disease

Alzheimer's disease (AD) is characterized by the presence of neuropsychiatric or behavioral and psychological symptoms of dementia (BPSD). BPSD have been associated with the APOE_ε4 allele, which is also the major genetic AD risk factor. Although the involvement of some circadian genes and orexin receptors in sleep and behavioral disorders has been studied in some psychiatric pathologies, including AD, there are no studies considering gene-gene interactions. The associations of one variant in PER2, two in PER3, two in OX2R and two in APOE were evaluated in 31 AD patients and 31 cognitively healthy subjects. Genotyping was performed using real-time PCR and capillary electrophoresis from blood samples. The allelic-genotypic frequencies of variants were calculated for the sample study. We explored associations between allelic variants with BPSD in AD patients based on the NPI, PHQ-9 and sleeping disorders questionnaires. Our results showed that the APOE_ε4 allele is an AD risk variant (p = 0.03). The remaining genetic variants did not reveal significant differences between patients and controls. The PER3_rs228697 variant showed a nine-fold increased risk for circadian rhythm sleep-wake disorders in Mexican AD patients, and our gene-gene interaction analysis identified a novel interaction between PERIOD and APOE gene variants. These findings need to be further confirmed in larger samples.

Genetic loci, rs17817449 and rs6567160, known for obesity and the risk of stroke events among middle-aged and older Chinese people

BACKGROUND

Fat Mass and Obesity-Associated (FTO) and the Melanocortin-4 Receptor (MC4R) genes are strongly associated with obesity, an established risk factor for stroke. We aimed to assess the associations between rs17817449 at the FTO and rs6567160 at the MC4R and the risk of stroke events in middle-aged and older Chinese people.

MATERIALS AND METHODS

Study data were obtained from the Guangzhou Biobank Cohort Study; a total of 148 participants with a self-reported history of stroke and an equal volume of age- and sex-matched participants were selected as the cases and the controls in a case-control study; a total of 13,967 participants at the first follow-up and all participants with fatal stroke (up to April 2021) were included in a retrospective cohort study. Conditional logistic regression and the Cox proportional hazards regression analyses were used to assess the associations of the two genetic loci with the risk of stroke events.

RESULTS

After adjusting for age, sex, education, job, smoking, alcohol consumption, body mass index, physical activity, hypertension, diabetes, and dyslipidemia, rs17817449 and rs6567160 shared minor alleles G and C, respectively, in the case-control analyses. The genotypes GG+GT of rs17817449 at the FTO were significantly associated with a decreased risk of fatal stroke occurrence, with fatal all strokes having an adjusted hazard ratio (aHR) of 0.71 (95% confidence intervals (CI) 0.52-0.97, P = 0.04) and fatal ischemic stroke having an aHR of 0.64 (95% CI 0.41-1.00, P = 0.05), when the genotype TT was taken as a reference and a series of multiplicities were adjusted; the risk of fatal all strokes was lowered by dyslipidemia (aHR = 0.63, 95% CI 0.39-1.00, P = 0.05) and non-diabetes (aHR = 0.68, 95% CI 0.46-0.99, P = 0.049) in the retrospective cohort analyses. Significances were observed neither in the associations between rs6567160 and the risk of stroke events nor in an interaction between rs17817449 and rs6567160 in the two-stage analyses.

CONCLUSION

The G allele of rs17817449 at the FTO, not rs6567160 at the MC4R, was associated with a decreased risk of fatal stroke occurrence; its functional role in stroke should be explored in relatively healthy middle-aged to older Chinese people.

Association of Polymorphisms in NHEJ Pathway Genes with HIV-1 Infection and AIDS Progression in a Northern Chinese MSM Population

Background and Aims

Men who have sex with men (MSM) are at high risk of HIV infection. The nonhomologous end joining (NHEJ) pathway is the main way of double-stranded DNA break (DSB) repair in the higher eukaryotes and can repair the DSB timely at any time in cell cycle. It is also indicated that the NHEJ pathway is associated with HIV-1 infection since the DSB in host genome DNA occurs in the process of HIV-1 integration. The aim of the present investigation was to evaluate associations of single-nucleotide polymorphisms (SNPs) in NHEJ pathway genes with susceptibility to HIV-1 infection and AIDS progression among MSM residing in northern China.

Methods

A total of 481 HIV-1 seropositive men and 493 HIV-1 seronegative men were included in this case-control study. Genotyping of 22 SNPs in NHEJ pathway genes was performed using the SNPscan™ Kit.

Results

Positive associations were observed between XRCC6 rs132770 and XRCC4 rs1056503 genotypes and the susceptibility to HIV-1 infection. In gene-gene interaction analysis, significant SNP-SNP interactions of XRCC6 and XRCC4 genetic variations were found to play a potential role in the risk of HIV-1 infection. In stratified analysis, XRCC5 rs16855458 was significantly associated with CD4+ T cell counts in AIDS patients, whereas LIG4 rs1805388 was linked to the clinical phases of AIDS patients.

Conclusions

NHEJ gene polymorphisms can be considered to be risk factors of HIV-1 infection and AIDS progression in the northern Chinese MSM population.

Interaction Between Vascular Endothelial Growth Factor Gene Polymorphism and Smoking on Gastric Cancer Risk in Chinese Han Population

Aim: In this study, we aimed to evaluate the associations of vascular endothelial growth factor (VEGF) gene single nucleotide polymorphisms (SNPs) and its interaction with current smoking with gastric cancer (GC) risk in the Chinese Han population.

Methods: We used logistic regression model to test the association between VEGF gene polymorphism and the risk of GC. The association strength was evaluated by odds ratio (OR) and 95% confidence interval (CI) calculated using logistic regression. Generalized multifactor dimensionality reduction (GMDR) was used to analyze the effect of the interaction between VEGF gene and current smoking on GC risk.

Results: Logistic regression analysis showed that the risk of GC was significantly higher in rs10434 -G allele carriers than that in AA genotype carriers (AG + GG and AA), and the adjusted OR (95% CI) = 1.64 (1.24–2.08). In addition, we found a significantly higher GC risk in subjects with rs833061-T allele than those with CC allele (CT + TT and CC), adjusted or (95% CI) = 1.43 (1.10–1.87). We also found a statistically significant two- locus model (p = 0.018), including rs10434 and current smoking, indicating a significant interaction between rs10434 and current smoking on the risk of GC. Hierarchical analysis found that current smokers with AG or GG genotype have the highest GC risk, compared to never- smokers with AA genotype, OR (95% CI) = 2.43 (1.64–3.28).

Conclusion: We found that rs10434 -G and rs833061-T alleles, gene- environment interaction between rs10434, and current smoking were all related to increased GC risk.

Fibrillin-1 Gene Polymorphisms (rs145233125, rs11070646, rs201170905) Are Associated With the Susceptibility and Clinical Prognosis of DeBakey Type III Aortic Dissection in Chinese Han Population

ABSTRACT

We aim to investigate whether genetic variants of the Fibrillin-1 (FBN1) gene were associated with DeBakey type III aortic dissection (AD) and its clinical prognosis in Chinese Han population. Three single-nucleotide polymorphisms (SNPs) (rs145233125, rs11070646, rs201170905) in FBN1 were analyzed in patients with DeBakey type III AD (159) and healthy subjects (216). Gene-environment interactions were evaluated to use generalized multifactor dimensionality reduction. Haplotype analysis of the 3 SNPs in the FBN1 gene was performed by Haploview software. Patients were followed up for average 4 years. G carriers of rs11070646 and rs201170905 in FBN1 have an increased risk of DeBakey type III AD. The interaction of FBN1 and environmental factors facilitated to the increased risk of DeBakey type III AD (cross-validation consistency = 10/10, P = 0.001). One of the most common haplotypes revealed an increased risk of DeBakey type III AD (CGG, P = 0.009). Recessive models of rs145233125 CC genotype ( P < 0.05) and rs201170905 GG genotype ( P < 0.001) were associated with an increased risk of death and recurrent chest pain of DeBakey type III AD. In conclusions, FBN1 gene polymorphisms contribute to DeBakey type III AD susceptibility. The interactions of gene and environment are related with the risk of DeBakey type III AD. C carriers of rs145233125 and G carriers of rs201170905 may be the adverse prognostic indicators of death and recurrent chest pain in DeBakey type III AD.

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.

Folate metabolism abnormalities in infertile patients with endometriosis

Background: Homocysteine levels can be impacted by enzymes variations. Aim: To correlate MTHFR, MTR and MTRR variants with homocysteine levels in the blood and follicular fluid and assisted reproduction results. Material & methods: MTHFR (rs2274976, rs1801131, rs1801133), MTR (rs1805087) and MTRR (rs1801394) genotyping was performed by TaqMan assays and compared with homocysteine levels, measured by ELISA, to oocytes retrieved and to the pregnancy status of women with endometriosis and controls. Results: The MTR G allele and GG genotype were more common in patients with endometriosis. They also showed lower levels of homocysteine and more clinical gestations. Epistasis analysis showed a model associated with gestational results, composed of MTHFR+MTR variants (CC+AG). Conclusion: The summation effect of variants in genes participating in folate metabolism was associated with pregnancy status in Brazilian women. MTR variants were more observed in endometriosis patients, as well as lower follicular Hcy levels and increased clinical pregnancy results.

Association of gene polymorphisms in FBN1 and TGF-β signaling with the susceptibility and prognostic outcomes of Stanford type B aortic dissection

Background

This study is aimed at investigating the association of Fibrillin-1 (FBN1) and transforming growth factor β (TGF-β) signaling-related gene polymorphisms with the susceptibility of Stanford type B aortic dissection (AD) and its clinical prognostic outcomes.

Methods

Five single-nucleotide polymorphism (SNPs) (FBN1rs 145233125, rs201170905, rs11070646, TGFB1rs1800469, and TGFB2rs900) were analyzed in patients with Stanford type B AD (164) and healthy controls (317). Gene–gene and gene–environment interactions were assessed by generalized multifactor dimensionality reduction. A 4-year follow-up was performed for all AD patients.

Results

G carriers of FBN1 rs201170905 and TGFB1 rs1800469 have an increased risk of Stanford type B AD. The interaction of FBN1, TGFB1, TGFB2 and environmental promoted to the increased risk of type B AD (cross-validation consistency = 10/10, P = 0.001). Dominant models of FBN1rs145233125 TC + CC genotype (P = 0.028), FBN1 rs201170905 AG + GG (P = 0.047) and TGFB1 rs1800469 AG + GG (P = 0.052) were associated with an increased risk of death of Stanford type B AD. The recessive model of FBN1 rs145233125 CC genotype (P < 0.001), FBN1rs201170905 GG (P < 0.001), TGFB1 rs1800469 AG + GG genotype (P = 0.011) was associated with an increased risk of recurrence of chest pain in Stanford type B AD.

Conclusions

The interactions of gene–gene and gene–environment are related with the risk of Stanford type B AD. C carriers of rs145233125, G carriers of rs201170905 and G carriers of rs1800469 may be the poor clinical outcome indicators of mortality and recurrent chest pain in Stanford type B AD.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01213-z.

Polymorphisms and Gene-Gene Interaction in AGER/IL6 Pathway Might Be Associated with Diabetic Ischemic Heart Disease

Background: Although the genetic susceptibility to diabetes and ischemic heart disease (IHD) has been well demonstrated, studies aimed at exploring gene variations associated with diabetic IHD are still limited; Methods: Our study included 204 IHD cases who had been diagnosed with diabetes before the diagnosis of IHD and 882 healthy controls. Logistic regression was used to find the association of candidate SNPs and polygenic risk score (PRS) with diabetic IHD. The diagnostic accuracy was represented with AUC. Generalized multifactor dimensionality reduction (GMDR) was used to illustrate gene-gene interactions; Results: For IL6R rs4845625, the CT and TT genotypes were associated with a lower risk of diabetic IHD than the CC genotype (OR = 0.619, p = 0.033; OR = 0.542, p = 0.025, respectively). Haplotypes in the AGER gene (rs184003-rs1035798-rs2070600-rs1800624) and IL6R gene (rs7529229-rs4845625-rs4129267-rs7514452-rs4072391) were both significantly associated with diabetic IHD. PRS was associated with the disease (OR = 1.100, p = 0.005) after adjusting for covariates, and the AUC were 0.763 (p < 0.001). The GMDR analysis suggested that rs184003 and rs4845625 were the best interaction model after permutation testing (p = 0.001) with a cross-validation consistency of 10/10; Conclusions: SNPs and haplotypes in the AGER and IL6R genes and the interaction of rs184003 and rs4845625 were significantly associated with diabetic IHD.

Gene-gene interaction and new onset of major depressive disorder: Findings from a Chinese freshmen nested case-control study

BACKGROUND

The gene-gene interaction is known to be the genetic cause of major depressive disorder (MDD). Several genes have been found to be related to MDD. The objectives of this study were to verify the susceptibility genes of MDD in a sample of university students in China, and to investigate possible gene-gene interactions in relation to the risk of MDD.

METHODS

7,627 Chinese Han freshmen were enrolled at baseline survey in 2018. After a 2-year follow-up, 170 new onset MDD cases and 680 controls with DNA samples reserved were sequenced and genotyped for 4 selected Single Nucleotide Polymorphisms (SNPs) in a nested case-control study (ratio of 1:4). Chi-square test was used to identify the relationships between SNPs and risk of MDD. Generalized multifactor dimensionality reduction (GMDR) was used to analyze the gene-gene interactions.

RESULTS

The 2-year incidence of MDD in Chinese college students was 3.75% (95% CI: 3.24%, 4.34%). There was no statistical difference in MDD incidences between males (3.74%, 95% CI: 3.12%, 4.49%) and females (3.77%, 95% CI: 2.97%, 4.78%) (p>0.05). TMEM161B (rs768705) was positively associated with new onset MDD (χ² = 0.75, p = 0.023). The AG genotype of rs768705 was significant (OR=1.640, 95%CI:1.414-2.358). The gene-gene interaction between TMEM161B (rs768705) and LHPP (rs35936514) was statistically significant in this nested case-control study (p = 0.011). The CV consistency was 9/10 and the testing accuracy was 0.5274.

LIMITATIONS

The results could not be inferred to other ethnics.

CONCLUSIONS

This study provided evidence that combined rs768705 (TMEM161B) and rs35936514 (LHPP) may modulate the risk of MDD.

Coronary Heart Disease in Type 2 Diabetes Mellitus: Genetic Factors and Their Mechanisms, Gene-Gene, and Gene-Environment Interactions in the Asian Populations

Asians are more susceptible to type 2 diabetes mellitus (T2D) and its coronary heart disease (CHD) complications than the Western populations, possibly due to genetic factors, higher degrees of obesity, insulin resistance, and endothelial dysfunction that could occur even in healthy individuals. The genetic factors and their mechanisms, along with gene-gene and gene-environment interactions associated with CHD in T2D Asians, are yet to be explored. Therefore, the objectives of this paper were to review the current evidence of genetic factors for CHD, summarize the proposed mechanisms of these genes and how they may associate with CHD risk, and review the gene-gene and gene-environment interactions in T2D Asians with CHD. The genetic factors can be grouped according to their involvement in the energy and lipoprotein metabolism, vascular and endothelial pathology, antioxidation, cell cycle regulation, DNA damage repair, hormonal regulation of glucose metabolism, as well as cytoskeletal function and intracellular transport. Meanwhile, interactions between single nucleotide polymorphisms (SNPs) from different genes, SNPs within a single gene, and genetic interaction with environmental factors including obesity, smoking habit, and hyperlipidemia could modify the gene's effect on the disease risk. Collectively, these factors illustrate the complexities of CHD in T2D, specifically among Asians.

Using Machine Learning to Predict Obesity Based on Genome-Wide and Epigenome-Wide Gene-Gene and Gene-Diet Interactions

Obesity is associated with many chronic diseases that impair healthy aging and is governed by genetic, epigenetic, and environmental factors and their complex interactions. This study aimed to develop a model that predicts an individual's risk of obesity by better characterizing these complex relations and interactions focusing on dietary factors. For this purpose, we conducted a combined genome-wide and epigenome-wide scan for body mass index (BMI) and up to three-way interactions among 402,793 single nucleotide polymorphisms (SNPs), 415,202 DNA methylation sites (DMSs), and 397 dietary and lifestyle factors using the generalized multifactor dimensionality reduction (GMDR) method. The training set consisted of 1,573 participants in exam 8 of the Framingham Offspring Study (FOS) cohort. After identifying genetic, epigenetic, and dietary factors that passed statistical significance, we applied machine learning (ML) algorithms to predict participants' obesity status in the test set, taken as a subset of independent samples (n = 394) from the same cohort. The quality and accuracy of prediction models were evaluated using the area under the receiver operating characteristic curve (ROC-AUC). GMDR identified 213 SNPs, 530 DMSs, and 49 dietary and lifestyle factors as significant predictors of obesity. Comparing several ML algorithms, we found that the stochastic gradient boosting model provided the best prediction accuracy for obesity with an overall accuracy of 70%, with ROC-AUC of 0.72 in test set samples. Top predictors of the best-fit model were 21 SNPs, 230 DMSs in genes such as CPT1A, ABCG1, SLC7A11, RNF145, and SREBF1, and 26 dietary factors, including processed meat, diet soda, French fries, high-fat dairy, artificial sweeteners, alcohol intake, and specific nutrients and food components, such as calcium and flavonols. In conclusion, we developed an integrated approach with ML to predict obesity using omics and dietary data. This extends our knowledge of the drivers of obesity, which can inform precision nutrition strategies for the prevention and treatment of obesity. Clinical Trial Registration: [www.ClinicalTrials.gov], the Framingham Heart Study (FHS), [NCT00005121].

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.

Genetic variants of PON1, GSTM1, GSTT1, and locus 9p21.3, and the risk for premature coronary artery disease in Yucatan, Mexico

OBJECTIVE

Genetic variants of PON1, rs70587, rs662, rs854560, GSTM1, and GSTT1 and two single nucleotide polymorphisms (SNP) at 9p21.3 locus, rs1333049, and rs2383207; were evaluated in association with the risk for premature coronary artery disease (CAD) in a population of Yucatan, Mexico. These genes are involved in the inactivation of pro-oxidants and pro-inflammatory mediators, lipid and xenobiotic metabolism, detoxification of reactive oxygen species, and regulation of cellular proliferation playing key roles in the pathogenesis of atherosclerosis.

METHODS

We conducted a matched case-control study with 98 CAD cases and 101 healthy controls. Genotyping of PON1 and 9p21.2 SNP was performed by real time-PCR and for GSTM1 and GSTT1 with multiplex-PCR. Odds ratios (OR) were calculated to estimate association and generalized multifactor dimensionality reduction (GMDR) algorithm to identify gene-gene and gene-environment interactions.

RESULTS

The distribution of all allele/genotype frequencies in controls was within Hardy-Weinberg expectations (p > .05) except for GSTM1. The allele/genotype frequencies of the GSTT1 null were significantly higher in CAD cases than in controls, suggesting association with higher risk for developing CAD. The other SNPs did not show any significant independent association with premature CAD. GMDR revealed a significant interaction between GSTT1 and LL55 genotype. Likewise, the body mass index (BMI) and smoking also showed an interaction with GSTT1.

CONCLUSION

The GSTT1 null allele/genotype is associated with an increased risk of developing premature CAD, the effect of which is not modified by cardiovascular risk factors in the population of Yucatan.

Association of gene polymorphisms in MYH11 and TGF-β signaling with the susceptibility and clinical outcomes of DeBakey type III aortic dissection

To investigate the association of myosin heavy chain protein 11 (MYH11) and transforming growth factor β signaling-related gene polymorphisms with the susceptibility of DeBakey type III aortic dissection (AD) and its clinical outcomes. Four single-nucleotide polymorphism (SNPs) (MYH11 rs115364997, rs117593370, TGFB1 rs1800469, and TGFBR1 rs1626340) were analyzed in patients with DeBakey III AD (173) and healthy participants (335). Gene-gene and gene-environment interactions were evaluated using generalized multifactor dimensionality reduction. The patients were followed up for a median of 55.7 months. MYH11 rs115364997 G or TGFBR1 rs1626340 A carriers had an increased risk of DeBakey type III AD. MYH11, TGFB1, TGFBR1, and environment interactions contributed to the risk of DeBakey type III AD (cross-validation consistency = 10/10, P = 0.001). Dominant models of MYH11 rs115364997 AG + GG genotype (HR = 2.443; 95%CI: 1.096-5.445, P = 0.029), TGFB1 rs1800469 AG + GG (HR = 2.303; 95%CI: 1.069-4.96, P = 0.033) were associated with an increased risk of mortality in DeBakey type III AD. The dominant model of TGFB1 rs1800469 AG + GG genotype was associated with an increased risk of recurrence of chest pain in DeBakey type III AD (HR = 1.566; 95%CI: 1.018-2.378, P = 0.041). In conclusions, G carriers of MYH11 rs115364997 or TGFB1 rs1800469 may be the poor prognostic indicators of mortality and recurrent chest pain in DeBakey type III AD. The interactions of gene-gene and gene-environment are associated with the risk of DeBakey type III AD.

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.

p53 p.Pro72Arg (rs1042522) and Mouse Double Minute 2 (MDM2) Single-Nucleotide Polymorphism (SNP) 309 Variants and Their Interaction in Chronic Lymphocytic Leukemia(CLL): A Survey in CLL Patients from Western Iran

Background: Chronic lymphocytic leukemia (CLL) is the most common leukemia in adults. The MDM2 and p53 are interacting proteins that play crucial roles in cell biology. Genetic variations of p53 and MDM2 have been identified in many cancers including CLL; among which are SNP309 in the promoter of MDM2 and SNP codon72 in p53. Materials and Methods: In this study, we sought to find the impact of two SNPs of p53 and MDM2 in the pathogenesis of CLL. A total of 100 CLL patients and 102 healthy controls were recruited. Genomic DNA was extracted, and genotyping was performed using the PCR-RFLP method. The allele and genotype associations were analyzed using the χ2 test. The gene-gene interaction analysis was studied using GMDR v0.9. Results: Our study found the absence of a significant difference between CLL patients and controls related to the allelic frequencies or genotypic distributions for both MDM2 SNP309 and p53 codon72. A significantly higher frequency of p53 C allele was found in patients with disease duration of more than 36 compared to those less than 36 months. However, GMDR analysis suggests genetic interaction between the genes under study. Conclusion: Our findings indicated each polymorphism of p53 codon72 and MDM2 (SNP309) was not a risk factor for CLL but the p53 C allele could be associated with the disease duration. Besides, the interaction between p53/MDM2 genotypes may confer susceptibility to CLL. Our study could be useful in genetic association studies of CLL and the role of gene-gene interactions in the susceptibility to the disease.

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.

Associations Between Cytochrome P450 (CYP) Gene Single-Nucleotide Polymorphisms and Second-to-Fourth Digit Ratio in Chinese University Students

BACKGROUND Cytochrome P450 (CYP) genes are necessary for the production or metabolism of fetal sex hormones during pregnancy. The second-to-fourth digit ratio (2D: 4D) is formed in the early stage of human fetal development and considered an indicator reflecting prenatal sex steroids levels. We explored the association between 2D: 4D and single-nucleotide polymorphisms (SNPs) of CYP. MATERIAL AND METHODS Correlation analysis between 2D: 4D and 8 SNPs, rs2687133 (CPY3A7), rs7173655 (CYP11A1), rs1004467, rs17115149, and rs2486758 (CYP17A1), and rs4646, rs2255192, rs4275794 (CYP19A1), was performed using data from 426 female and 412 male Chinese university students. SNP genotyping was conducted using PCR. Digit lengths were photographed and measured by image processing software. RESULTS rs2486758 (CYP17A1) correlated with left hand 2D: 4D in men (P=0.026), and rs1004467 (CYP17A1) correlated with right hand 2D: 4D in men (P=0.008) and the whole population (P=0.032). In men, allele G rs1004467 decreased right hand 2D: 4D, while allele C of rs2486758 increased left hand 2D: 4D. In women, left hand 2D: 4D was higher in genotypes with allele A of SNP rs4646 (CYP19A1) under the dominant genetic model; female DR-L was higher in genotypes with allele T of rs17115149 (CYP11A1). SNPs rs2687133 (CYP3A7) and rs1004467 (CYP17A1) were significantly correlated with right hand 2D: 4D (P=0.0107). CONCLUSIONS SNPs rs1004467 and rs2486758 of CYP17A1 are significant in the relationship between 2D: 4D and CYP gene polymorphisms under different conditions. SNP interactions between CYP genes probably impact 2D: 4D. The correlation between 2D: 4D and some sex hormone-related diseases may be due to the effect of CYP variants on the 2 phenotypes.

Fluoride exposure, dopamine relative gene polymorphism and intelligence: A cross-sectional study in China

BACKGROUND

Excessive fluoride exposure is related to adverse health outcomes, but whether dopamine (DA) relative genes are involved in the health effect of low-moderate fluoride exposure on children's intelligence remain unclear.

OBJECTIVES

We conducted a cross-sectional study to explore the role of DA relative genes in the health effect of low-moderate fluoride exposure in drinking water.

METHODS

We recruited 567 resident children, aged 6-11 years old, randomly from endemic and non-endemic fluorosis areas in Tianjin, China. Spot urine samples were tested for urinary fluoride concentration, combined Raven`s test was used for intelligence quotient test. Fasting venous blood were collected to analyze ANKK1 Taq1A (rs1800497), COMT Val158Met (rs4680), DAT1 40 bp VNTR and MAOA uVNTR. Multivariable linear regression models were used to assess associations between fluoride exposure and IQ scores. We applied multiplicative and additive models to appraise single gene-environment interaction. Generalized multifactor dimensionality reduction (GMDR) was used to evaluate high-dimensional interactions of gene-gene and gene-environment.

RESULTS

In adjusted model, fluoride exposure was inversely associated with IQ scores (β = -5.957, 95% CI: -9.712, -2.202). The mean IQ scores of children with high-activity MAOA genotype was significantly lower than IQ scores of those with low-activity (P = 0.006) or female heterozygote (P = 0.016) genotype. We detected effect modification by four DA relative genes (ANKK1, COMT, DAT1 and MAOA) on the association between UF and IQ scores. We also found a high-dimensional gene-environment interaction among UF, ANKK1, COMT and MAOA on the effect of IQ (testing balanced accuracy = 0.5302, CV consistency: 10/10, P = 0.0107).

CONCLUSIONS

Our study suggests DA relative genes may modify the association between fluoride and intelligence, and a potential interaction among fluoride exposure and DA relative genes on IQ.

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.

Association analysis between the interaction of RAS family genes mutations and papillary thyroid carcinoma in the Han Chinese population

Papillary thyroid carcinoma (PTC) is the major subtype of thyroid cancer, accounting for 75%-85% of all thyroid malignancies. This study aimed to identify the association between the interactions of single nucleotide polymorphisms (SNPs) in RAS family genes and PTC in the Han Chinese population, to provide clues to the pathogenesis and potential therapeutic targets for PTC. Hap Map and NCBI-db SNP databases were used to retrieve SNPs. Haploview 4.2 software was used to filter SNPs based on specific parameters, six SNPs of RAS gene (KRAS-rs12427141, KRAS-rs712, KRAS-rs7315339, HRAS-rs12628, NRAS-rs14804 and NRAS-rs2273267) were genotyped by matrix-assisted laser desorption/ionization time of flight mass spectrometry (MALDI-TOF-MS) in 673 PTC patients and 657 healthy controls, the interactive effect was evaluated by crossover analysis, logistic regression and GMDR software. We found that genetic mutation in rs712 have significant associations with PTC risk after Bonferroni correction (p<0.001). The interaction between KRAS-rs12427141 and HRAS-rs12628 increased the risk of PTC (U=-2.119, p<0.05), the interaction between KRAS-rs2273267 and HRAS-rs7315339 reduced the risk of PTC (U=2.195, p<0.05). GMDR analysis showed that the two-factor model (KRAS-rs712, NRAS-rs2273267) was the best (p=0.0107). Summarily, there are PTC-related interactions between RAS family genes polymorphisms in the Han Chinese population.

Pro-inflammatory cytokine polymorphisms in ONECUT2 and HNF4A and primary colorectal carcinoma: a post genome-wide gene-lifestyle interaction study

Immune-related molecular and genetic pathways that are connected to colorectal cancer (CRC) and lifestyles in postmenopausal women are incompletely characterized. In this study, we examined the role of pro-inflammatory biomarkers such as C-reactive protein (CRP) and interleukin-6 (IL-6) in those pathways. Through selection of the best predictive single-nucleotide polymorphisms (SNPs) and lifestyles, our goal was to improve the prediction accuracy and ability for CRC risk. Using large cohort data of postmenopausal women from the Women's Health Initiative Database for Genotypes and Phenotypes Study, we previously conducted a genome-wide association (GWA) for a CRP and IL-6 gene-behavioral interaction study. For the present study, we added GWA-SNPs from outside GWA studies, resulting in a total of 152 SNPs. Together with 41 selected lifestyles, we performed a 2-stage multimodal random survival forest analysis with generalized multifactor dimensionality reduction approach to construct CRC risk profiles. Overall and in obesity strata (by body mass index, waist circumference, waist-to-hip ratio, exercise, and dietary fat intake), we identified the best predictive genetic markers in inflammatory cytokines and lifestyles. Across the strata, 2 SNPs (ONECUT2 rs4092465 and HNF4A rs1800961) and 1 lifestyle factor (relatively short-term past use of oral contraceptives) were the most common and strongest predictive markers for CRC risk. The risk profile that combined those variables exhibited synergistically increased risk for CRC; this pattern appeared more strongly in obese and inactive subgroups. Our results may contribute to improved predictability for CRC and suggest genetically targeted lifestyle interventions for women carrying the inflammatory-risk genotypes, reducing CRC risk.

The MC4R SNPs, their haplotypes and gene-environment interactions on the risk of obesity

Background

Little is known about the correlation between the melanocortin 4 receptor gene (MC4R) single nucleotide polymorphisms (SNPs) and the risk of obesity. This research sought to test the MC4R rs17782313, rs476828 and rs12970134 SNPs, their haplotypes and gene-environment interactions on the risk of obesity in the Maonan ethnic group, an isolated minority in China.

Methods

A case-control study comprised of 1836 participants (obesity group, 858; and control group, 978) was conducted. Genotypes of the three SNPs were determined by the next-generation sequencing (NGS) technology.

Results

The genotypic frequencies of the three SNPs were different between the obesity and control groups (P <  0.05 for all). The minor allelic frequency of the MC4R rs17782313C, rs476828C and rs12970134A was higher in obesity than in control groups (13.8% vs. 8.3%, P <  0.001, 17.1% vs. 10.9%, P <  0.001; and 15.5% vs. 11.5%, P <  0.001; respectively). Additionally, the dominant model of rs17782313 and rs476828 SNPs revealed an increased morbidity function on the risk of obesity (P <  0.05). A correlation between SNP-environment and the risk of obesity was also observed. The rs17782313C-rs476828C-rs12970134A haplotype was associated with high risk of obesity (OR = 1.796, 95% CI = 1.447–2.229), whereas the rs17782313T-rs476828T-rs12970134G and rs17782313T-rs476828T-rs12970134A haplotypes were associated with low risk of obesity (OR = 0.699, 95% CI = 0.586–0.834 and OR = 0.620, 95% CI = 0.416–0.925; respectively). The interactions between haplotype and waist circumference on the risk of obesity were also noted.

Conclusions

We discovered that the MC4R rs17782313, rs476828 and rs12970134 SNPs and their haplotypes were associated with the risk of obesity in the Chinese Maonan population.

Association studies of dopamine synthesis and metabolism genes with multiple phenotypes of heroin dependence

BACKGROUND

Heroin dependence is a complex disease with multiple phenotypes. Classification of heroin users into more homogeneous subgroups on the basis of these phenotypes could help to identify the involved genetic factors and precise treatments. This study aimed to identify the association between genetic polymorphisms of DA synthesis and metabolism genes, including tyrosine hydroxylase (TH), DOPA decarboxylase (DDC), solute carrier family 6 member 3 (SLC6A3) and DA beta-hydroxylase (DBH), with six important phenotypes of heroin dependence.

METHODS

A total of 801 heroin dependent patients were recruited and fourteen potential functional single nucleotide polymorphisms (SNPs) were genotyped by SNaPshot. Associations between SNPs with six phenotypes were mainly assessed by binary logistic regression. Generalized multifactor dimensionality reduction was used to analyze the gene-by-gene and gene-by-environment interactions.

RESULTS

We found that DBH rs1611114 TT genotype had a protective effect on memory impairment after heroin dependence (P = 0.002, OR = 0.610). We also found that the income-rs12666409-rs129915-rs1611114 interaction yielded the highest testing balance accuracy and cross-validation consistency for memory change after heroin dependence.

CONCLUSIONS

Our results suggest that the memory change after heroin dependence was a result of a combination of genetics and environment. This finding could lead to a better understanding of heroin dependence and further improve personalized treatment.

Gene-gene and gene-environment interaction data for platinum-based chemotherapy in non-small cell lung cancer

Gene-gene (GXG) and gene-environment (GXE) interactions play important roles in pharmacogenetics study. Simultaneously incorporating multiple single nucleotide polymorphisms (SNPs) and clinical factors is needed to explore the association of their interactions with drug response and toxicity phenotypes. We genotyped 504 SNPs in a total of 490 Chinese non-small cell lung cancer (NSCLC) patients, and the correlation of GXG and GXE interactions with platinum-based chemotherapeutic efficacy and safety were analyzed. In this data descriptor, we shared our data set which could help others to reuse them. All kinds of file types needed for GXG and GXE analysis were supplied. The process of genotyping and data analysis was also introduced step by step.

Genetic polymorphisms in CDH1 are associated with endometrial carcinoma susceptibility among Chinese Han women

The cadherin 1 (CDH1) gene plays critical roles in the epithelial-mesenchymal transition process, potentially offering us a glimpse into the development of endometrial carcinoma (EC). The present study aimed to identify whether genetic variants in CDH1 affect EC susceptibility in Chinese Han women, using a strategy combining haplotype-tagging single nucleotide polymorphisms (htSNPs) association analysis with fine-scale mapping. A total of 9 htSNPs in CDH1 were genotyped among 516 cases and 706 age-matched cancer-free controls. Logistic regression analyses revealed 3 htSNPs (rs17715799, rs6499199 and rs13689) to be associated with increased EC risk and 3 htSNPs (rs12185157, rs10431923 and rs4783689) with decreased EC risk. Furthermore, 14 newly imputed SNPs of CDH1 were identified to be associated with EC risk (P<0.05) using genotype imputation analysis. Notably, multivariate logistic analysis demonstrated that rs13689, rs10431923 and rs10431924 could affect EC susceptibility independently (P≤0.001). Subsequent Generalized Multifactor Dimensionality Reduction analysis revealed several best fitting models for predicting EC risk, including SNP-SNP interactions among rs7100190, rs12185157, rs10431923, rs7186053, rs6499199, rs4783689, rs13689, rs6499197 and rs10431924, and SNP-environment interactions between related SNPs and number of childbirth. Moreover, functional annotations suggest that the majority of these susceptible variants may carry potential biological functions that affect certain gene regulatory elements. In summary, this study suggested that the genetic polymorphisms of CDH1 were indeed associated with EC susceptibility on several levels. If further additional functional studies could verify these findings, these genetic variants may serve as future personalized markers for the early prediction of endometrial cancer in Chinese Han women.

Triangular relationship between CYP2R1 gene polymorphism, serum 25(OH)D3 levels and T2DM in a Chinese rural population

BACKGROUND

A low serum vitamin D concentration is associated with an increased risk of type 2 diabetes mellitus (T2DM). Recently, several single nucleotid polymorphisms (SNPs) have been identified which influence vitamin D levels. If a causal relationship exists between vitamin D concentrations and T2DM, one would expect a similar association between the newly identified SNPs and T2DM risk. Therefore, this study investigated the association between four SNPs of cytochrome P450 family 2, subfamily R, peptide 1 (CYP2R1) gene, serum 25(OH)D₃ levels and T2DM.

METHODS

Three hundred and ninety-seven patients with confirmed T2DM, as well as 397 age- and gender-matched controls were enrolled in this case-control study. Genotyping was performed by TaqMan probe assays. Kruskal-Wallis one-way analysis and muitiple logistic regression analysis were performed to identify the possible risk genotype for vitamin D levels and T2DM, respectively. Generalized multifactor dimensionality reduction (GMDR) was used to analyze the gene-gene and gene-environment interactions.

RESULTS

The serum 25(OH)D₃ levels were significant lower in the T2DM group. Significant differences were observed between patients and controls in terms of the genotype distributions of rs1993116 (P = 0.048) and rs10766197 (P = 0.024). Similarly, rs1993116 and rs10766197 polymorphisms were found to be significantly associated with T2DM risk. AG + GG genotype carriers of the rs1993116 and rs10766197 polymorphisms could have an increased risk of developing T2DM compared with AA carriers, the OR and 95% CI were 1.64 (1.09-2.46) and 1.76 (1.18-2.65), respectively. However, none of the tested SNPs were independently associated with serum 25(OH)D₃ levels (P > 0.059). Gene-gene and gene-environment interaction analyses indicated that rs12794714-rs10766197 and rs12794714-vitamin D deficiency (VDD) models successfully predicted T2DM risk (P < 0.001).

CONCLUSIONS

Rs1993116 and rs10766197 polymorphisms of CYP2R1 gene may be novel genetic markers for T2DM in China. Given the lack of association between SNPs and serum 25(OH)D₃ levels, well-designed future studies should be conducted with larger sample sizes in rural areas of China.

Interactions among genes involved in reverse cholesterol transport and in the response to environmental factors in dyslipidemia in subjects from the Xinjiang rural area

Gene-gene and gene-environment interactions may be partially responsible for dyslipidemia, but studies investigating interactions in the reverse cholesterol transport system (RCT) are limited. We explored these interactions in a Xinjiang rural population by genotyping five SNPs using SNPShot technique in APOA1, ABCA1, and LCAT, which are involved in the RCT (690 patients, 743 controls). We conducted unconditional logistical regression analysis to evaluate associations and generalized multifactor dimensionality reduction to evaluate interactions. Results revealed significant differences in rs670 and rs2292318 allele frequencies between cases and controls (P<0.025). rs670 G allele carriers were more likely to develop dyslipidemia than A allele carriers (OR = 1.315, OR 95% CI: 1.067-2.620; P = 0.010). rs2292318 T allele carriers were more likely to develop dyslipidemia than A allele carriers (OR = 1.264, OR 95% CI: 1.037-1.541; P = 0.020). Gene-gene interaction model APOA1rs670-ABCA1rs1800976-ABCA1rs4149313-LCATrs1109166 (P = 0.0107) and gene-environment interaction model ABCA1rs1800976-ABCA1rs4149313-LCATrs1109166-obesity-smoking were optimal dyslipidemia predictors (P = 0.0107) and can interact (4). Differences in A-C-A-C-A and G-G-G-T-G haplotype frequencies were observed (P<0.05). Serum lipid profiles could be partly attributed to RCT gene polymorphisms. Thus, dyslipidemia is influenced by APOA1, ABCA1, LCAT, environmental factors, and their interactions.

UCP2 and UCP3 variants and gene-environment interaction associated with prediabetes and T2DM in a rural population: a case control study in China

Background

There are disparities for the association between uncoupling proteins (UCP) and type 2 diabetes (T2DM). The study was to examine the associations of genetic variants of UCP2 and UCP3 with prediabetes and T2DM in a rural Chinese population.

Methods

A population-based case-control study of 397 adults with T2DM, 394 with prediabetes and 409 with normal glucose tolerance (NGT) was carried out in 2014 in a rural community in eastern China. Three groups were identified through a community survey and the prediabetes and NGT groups were frequently matched by age and gender with the T2DM group and they were not relatives of T2DM subjects. With r² ≥ 0.8 and minor allele frequency (MAF) ≥0.05 for tag single nucleotide polymorphisms (SNPs) with potential function, three (rs660339, rs45560234 and rs643064) and six (rs7930460, rs15763, rs647126, rs1800849, rs3781907 and rs1685356) SNPs were selected respectively for UCP2 and UCP3 and genotyped in real time using the MassARRAY system (Sequenom; USA). The haplotypes, gene-environmental interaction and association between genetic variants of UCP2 and UCP3 and prediabetes or T2DM were explored.

Results

There were no significant differences in age and sex among three study groups. After the adjustment for possible covariates, the A allele of rs1800849 in UCP3 was significantly associated with prediabetes (aOR_{AA vs GG} = 1.68, 95% CI: 1.02–2.78), and the association was also significant under the recessive model (aOR _{AA vs GA + GG} = 1.64, 95% CI: 1.02–2.66). Also, rs15763 was found to be marginally significantly associated with T2DM under dominant model (OR_{GA + AA vs GG} = 0.73, 95% CI: 0.52–1.03, P = 0.072). No haplotype was significantly associated with prediabetes or T2DM. Multiplicative interactions for rs660339-overweight on T2DM were observed. In addition, the AA genotype of rs660339 was associated with an increased risk of T2DM in overweight subjects (OR = 1.48, 95%CI: 0.87–2.52) but with a decreased risk in those with normal weight (OR = 0.54, 95%CI: 0.28–1.05).

Conclusions

Rs1800849 in UCP3 was significantly associated with prediabetes. Overweight might modify the effects of rs660339 of UCP2 on T2DM.

Electronic supplementary material

The online version of this article (10.1186/s12881-018-0554-4) contains supplementary material, which is available to authorized users.