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.

Effects of environmental and genetic interactions on job burnout in coal miners: interactions between occupational stress, coping styles, and NR3C2 gene polymorphisms

Objectives

To investigate the current situation regarding occupational burnout among coal miners, explore the relationship between NR3C2 gene polymorphism and occupational burnout, and analyze the influence of the interaction between environment and gene on occupational burnout. This study provides a scientific basis for formulating health strategies to combat job burnout.

Methods

A total of 1,500 first-line coal mine workers were selected by cluster random sampling, and the job burnout scale, job content questionnaire (JCQ), and simplified coping style questionnaire (SCSQ) were used for the questionnaire survey. A total of 150 workers were randomly selected from the high burnout group and the low burnout group, and a total of 300 workers were selected as the research objects to examine the relationship between gene polymorphism, environment-gene interactions and burnout. This study employed iMLDRTM genotyping technology for NR3C2 gene (rs5522, rs2070950) polymorphism analysis. The relationship between the occurrence of job burnout, occupational stress, coping styles and the NR3C2 gene was analyzed.

Results

Finally, a total of 1,282 valid questionnaires were retrieved, with an effective recovery rate of 85.5%. The study included 128 participants (10%) with zero burnout, 400 (31.2%) with mild burnout, 649 (50.6%) with moderate burnout and 105 (8.2%) with severe burnout. There were significant differences in the rate of burnout among miners with respect to sex, age, working years, educational level, shifts, and marital status (P < 0.05). The difference in occupational stress between the different job burnout groups was statistically significant (P < 0.05). Compared with the GG genotype of rs2070950 of the NR3C2 gene, the CC genotype was identified as a susceptibility gene for occupational burnout (P < 0.05). In respect to rs5522, rs2070950, occupational stress, positive coping, and negative coping, the low-risk group was unlikely to suffer from job burnout compared with the high-risk group (OR = 0.103, 95%CI: 0.058-0.182).

Conclusion

In addition to demographic characteristics, occupational stress and negative coping styles were also identified as risk factors for job burnout. The interaction between locus rs5522, locus rs2070950, occupational stress, positive response, and negative response were found to affect the incidence of occupational burnout.

Effect of ACE, ACE2 and CYP11B2 gene polymorphisms and noise on essential hypertension among steelworkers in China: a case-control study

Background

Previous studies on the relationship between ACE I/D, ACE2 G8790A and CYP11B2-344T/C gene polymorphisms and essential hypertension (EH) were inconsistent. Moreover, few studies have reported the combined effect of these gene polymorphisms and noise exposure on EH. The purpose of this study was to explore the combined and separate effects of ACE I/D, ACE2 G8790A and CYP11B2-344T/C gene polymorphisms and noise on EH among steelworkers.

Methods

A case–control study was conducted on 725 male workers between March 2014 and July 2014 in the Tangsteel Company, China. The noise exposure of the workers were measured. Logistic regression and crossover analysis were used to analyse the effects of the interactions on the EH among steelworkers. GMDR was used to determine the best combination model of gene–noise interactions.

Results

Multivariate logistic regression showed that noise exposure increased the odds of EH, and the OR is 1.52 (95% CI 1.04–2.22). The risk of having EH for ACE I/D DD genotype carriers was 1.99 times that for II genotype carriers (95% CI 1.14–3.51). There was a negative additive interaction between ACE2 G8790A and CYP11B2-344T/C on EH (U3 =  − 2.221, P = 0.026, and S = 0.128) and a positive multiplicative interaction between ACE I/D and CYP11B2-344T/C on essential hypertension (P = 0.041). In addition, there was no significant gene–noise interaction model through the GMDR method after adjusting the confounders.

Conclusions

The ACE DD genotype may make men susceptible to EH. Simultaneously carrying the DD genotype of ACE I/D and the TC genotype of CYP11B2-344T/C increased the risk of EH.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01177-0.

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

Introduction

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

Methods

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

Results

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

Conclusions

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

Trial registration

Chictr.org.cn, identifier, ChiCTR2000033576.

Supplementary Information

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

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

OBJECTIVES

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

METHODS

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

RESULTS

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

CONCLUSION

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

Vitamin D levels and Vitamin D-related gene polymorphisms in Chinese children with type 1 diabetes

Low vitamin D levels may play a role in type 1 diabetes (T1D) susceptibility. Since 25(OH)D synthesis is genetically regulated, single nucleotide polymorphisms (SNPs) of important genes have also been shown to modulate the risk of T1D, so this study aimed to investigate the relationship between five SNPs in CYP2R1, DHCR7, CYP24A1, VDR genes, serum 25(OH)D levels and T1D in Chinese children. This case-control study included 141 T1D patients and 200 age-matched healthy children.25 (OH) D concentration was determined, genotyping was performed by High resolution melting (HRM). There was a significant difference in the prevalence of vitamin D deficiency, insufficiency, and sufficiency between T1D and healthy controls. (χ ² = 10.86, p = 0.004), however no evidence of the association between any group of SNPs and circulating 25(OH) D levels was observed. The allele distribution of CYP2R1(rs1993116) was significantly different between T1D and control group (p = 0.040), and the C allele carriers of rs1993116 had a higher risk of T1D than the T allele carriers, Carriers of the CC and CT genotypes of rs1993116 have higher T1D risk than those carrying the TT genotype. GMDR analysis revealed a significant interaction between CYP2R1(rs12794714) and CYP2R1(rs1993116) in the risk of T1D with a maximum testing balance accuracy of 60.39%.

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

Objectives

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

Methods

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

Results

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

Conclusions

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

Association of Variant Interactions in RANK, RANKL, OPG, TRAF6, and NFATC1 Genes with the Development of Osteonecrosis of the Femoral Head

Multiple gene polymorphisms have been demonstrated to correlate with the susceptibility to osteonecrosis of the femoral head (ONFH). However, as a complex disease induced by multiple genes, the development of ONFH has rarely been reported to involve in gene interaction. In this study, we first explored the association of 10 variants interactions in receptor activator of nuclear factor-kappa B (RANK), RANK ligand (RANKL), osteoprotegerin (OPG), tumor necrosis factor receptor-associated factor 6 (TRAF6), and nuclear factor of activated T cells cytoplasmic 1 (NFATC1) genes with the development and clinical phenotypes of ONFH in a 377 ONFH case-control study with using Mass ARRAY^® platform. Our results showed that not only a total of 6 interactional variants in the paired 10 variants interactions were significantly associated with the development of ONFH (OPG rs2073617 and NFATC1 rs754093, p < 0.019; OPG rs2073618 and NFATC1 rs754093, p < 0.008; OPG rs2073617 and RANKL rs1054016, p < 0.039, respectively) but also a total of 4 paired interactional variants were found to involve significantly in the increased risk of bilateral hip lesions in ONFH (OPG rs2073617 and TRAF6 rs5030411, p = 0.044; RANK rs884205 and TRAF6 rs5030411, p = 0.045, respectively). Moreover, the results from generalized multifactor dimensionality reduction also showed that the five best models were identified and associated significantly with ONFH risk, p = 0.001, 0.01, 0.01, 0.01, and 0.01, respectively. Our results first suggest that the variants in RANK/RANKL/OPG pathway genes affected the development of ONFH in gene interaction manner through the interaction of the paired variants and multiple variants.

The TXA2R rs1131882, P2Y1 rs1371097 and GPIIIa rs2317676 three-loci interactions may increase the risk of carotid stenosis in patients with ischemic stroke

Background

The genetic risk factors for carotid stenosis are not fully understood. The aim of this study is to investigate the relationship between variants in platelet activation-relevant genes and carotid stenosis in patients with ischemic stroke (IS).

Methods

Eleven variants of platelet activation-relevant genes, aggregates of platelet-leukocyte, and platelet aggregation were examined in 236 IS patients with carotid stenosis and 378 patients without carotid stenosis. High-resolution B-mode ultrasound was used to assess carotid stenosis. Generalized multifactor dimensionality reduction (GMDR) methods were applied in analyzing gene-gene interactions to determine whether there was any interactive role of assessed variants in affecting risk of carotid stenosis.

Results

Platelet aggregation and aggregates of platelet-leukocyte showed higher value in patients with carotid stenosis, compared with patients without carotid stenosis. Excluding potential disturbance variables, these 11 variants were not associated with carotid stenosis. However, according to the GMDR analysis, gene-gene interactions among TXA2R rs1131882, P2Y1 rs1371097 and GPIIIa rs2317676 had a synergistic influence on carotid stenosis. The high-risk interactions between the three variants showed a relationship with higher platelet activation, and have independent associations with risk of carotid stenosis (OR = 2.72, 95% CI: 1.28–7.82, P = 0.001).

Conclusion

The interactions among rs1131882, rs1371097 and rs2317676 perhaps increase the risk of symptomatic carotid stenosis, and maybe a potential marker for carotid stenosis. In this study, the combinatorial analysis made good use in elucidating complex risk factors in the heredity of carotid stenosis.

High genetic risk scores of SLIT3, PLEKHA5 and PPP2R2C variants increased insulin resistance and interacted with coffee and caffeine consumption in middle-aged adults

BACKGROUNDS AND AIMS

Insulin resistance is a common feature of metabolic syndrome that may be influenced by genetic risk factors. We hypothesized that genetic risk scores (GRS) of SNPs that influence insulin resistance and signaling interact with lifestyles to modulate insulin resistance in Korean adults.

METHODS AND RESULTS

Genome-wide association studies (GWAS) of subjects aged 40-65 years who participated in the Ansung/Ansan cohorts (8842 adults) in Korea revealed 52 genetic variants that influence insulin resistance. The best gene-gene interaction model was explored using the generalized multifactor dimensionality reduction (GMDR) method. GRS from the best model were calculated and the GRS were divided into low, medium and high groups. The best model for representing insulin resistance included SLIT3_rs2974430, PLEKHA5_rs1077044, and PPP2R2C_rs16838853. The odds ratios for insulin resistance were increased by 150% in the High-GRS group compared to the Low-GRS group. However, ORs for insulin secretion capacity, measured by HOMA-B, were not associated with GRS. Coffee and caffeine intake and GRS had an interaction with insulin resistance: In subjects with high coffee (≥10 cups/week) or caffeine intake (≥220 mg caffeine/day), insulin resistance was significantly elevated in the High-GRS group, but not in the Low-GRS. However, alcohol intake, smoking and physical activity did not have an interaction with GRS. Insulin secretion capacity was not significantly influenced by GRS when evaluating the adjusted odds ratios.

CONCLUSIONS

Subjects with High-GRS may be susceptible to increased insulin resistance by 50% and its risk may be exacerbated by consuming more than 10 cups coffee/week or 220 mg caffeine/day.

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.

Potential Interactions Between Genetic Polymorphisms of the Transforming Growth Factor-β Pathway and Environmental Factors in Abdominal Aortic Aneurysms

OBJECTIVE/BACKGROUND

Evidence has accumulated that multiple polymorphisms in the transforming growth factor (TGF)-β pathway and renin-angiotensin system play important roles in determining susceptibility to abdominal aortic aneurysm (AAA). Few studies have considered interactions between these gene polymorphisms and environmental factors. The aim of this study was to evaluate the contribution of single nucleotide polymorphisms (SNPs) and complex gene-environment interactions in AAA.

METHODS

Six SNPs located in TGFB, TGFBR1, TGFBR2 and AGTR1 were selected. Genotyping of blood samples and collection of lifestyle factors were performed in 155 unrelated participants with AAAs and 310 non-AAA controls. Unconditional logistic regression was performed to assess the effects of SNPs on the risk of AAA. Generalized multifactor dimensionality reduction (GMDR) was used to evaluate gene-gene and gene-environment interactions.

RESULTS

Participants carrying TGFB1 rs1800469 TT (odds ratio [OR] 1.83, 95% confidence interval [CI] 1.18-2.85) or AGTR1 rs12695895 TT (OR 4.21, 95% CI 1.41-12.53) genotypes had a higher risk of AAA than those with the common CC genotype. The gene-gene interaction of AGTR1 rs5182, TGFBR1 rs1626340, and TGFB1 rs1800469 was found to be the best model according to the results of the GMDR analysis (cross validation consistency [CVC]) 10/10; p = .010). Smoking, dyslipidemia, and rs1800469 together contributed to the risk of AAA, which demonstrated a potential and complex gene-environment interaction among the three variants that might affect AAA risk (CVC 6/10; p = .001).

CONCLUSION

In this study of the Chinese population, homozygosity of TGFB1 rs1800469-T and AGTR1 rs12695895-T might be associated with increased risk of AAA. The complex gene-gene and gene-environment interactions might contribute to the risk of AAA. As a small study, the preliminary results need extensive validation and replication in larger populations.

Genetic polymorphisms of ALOX5AP and CYP3A5 increase susceptibility to ischemic stroke and are associated with atherothrombotic events in stroke patients

BACKGROUND

The contributions of gene-gene interactions to pathogenesis of stroke remain largely elusive. The present study was designed to investigate the associations between genetic variations and ischemic stroke risk, the roles of gene-gene interactions in ischemic stroke, and their associations with atherothrombotic events.

METHODS

Among 396 patients with ischemic stroke and 378 controls, we examined 8 variants from 5 genes, including ALOX5AP-SG13S32 (rs9551963), SG13S42 (rs4769060), SG13S89 (rs4769874), SG13S114 (rs10507391), EPHX2 G860A (rs751141), CYP2C9*2 C430T (rs1799853), CYP2C9*3 A1075C (rs1057910), and CYP3A5 A6986G (rs776746), using matrix-assisted laser desorption/ionization time of flight mass spectrometry. Gene-gene interactions were determined by the generalized multifactor dimensionality reduction (GMDR) method. All ischemic stroke patients were followed up 12 months for atherothrombotic events, including recurrent ischemic stroke and other vascular events.

RESULTS

Single-gene variant analysis showed no significant differences in the genotype distributions of the 8 variants between the 2 groups. However, the GMDR analysis showed a significant interaction between rs10507391 and rs776746, in those cases carrying rs10507391 AA and rs776746 GG, the risk of ischemic stroke increased by 2.014 times (95% confidence interval [CI], 1.896-6.299; P = .006), and the atherothrombotic events occurred more frequently in those patients during follow-up period (P < .001). Multiple Cox regression analysis showed that the interaction between rs10507391 AA and rs776746 GG was an independent risk factor for atherothrombotic events (relative risk = 2.921; 95% CI, 1.118-7.012; P = .008).

CONCLUSIONS

The interaction between rs10507391 and rs776746 increases the susceptibility to ischemic stroke and is associated with atherothrombotic events in stroke patients.