Genetic interaction in the association between oxidative stress and diabetes in the Spanish population

Oxidative stress (OS) is a relevant intermediate mechanism involved in Type 2 Diabetes Mellitus (T2D) development. To date, the interaction between OS parameters and variations in genes related to T2D has not been analyzed.

AIMS

To study the genetic interaction of genes potentially related to OS levels (redox homeostasis, renin-angiotensin-aldosterone system, endoplasmic stress response, dyslipidemia, obesity and metal transport) and OS and T2D risk in a general population from Spain (the Hortega Study) in relation to the risk of suffering from T2D.

MATERIALS AND METHODS

One thousand five hundred and two adults from the University Hospital Rio Hortega area were studied and 900 single nucleotide polymorphisms (SNPs) from 272 candidate genes were analyzed.

RESULTS

There were no differences in OS levels between cases and controls. Some polymorphisms were associated with T2D and with OS levels. Significant interactions were observed between OS levels and two polymorphisms in relation to T2D presence: rs196904 (ERN1 gene) and rs2410718 (COX7C gene); and between OS levels and haplotypes of the genes: SP2, HFF1A, ILI8R1, EIF2AK2, TXNRD2, PPARA, NDUFS2 and ERN1.

CONCLUSIONS

Our results indicate that genetic variations of the studied genes are associated with OS levels and that their interaction with OS parameters may contribute to the risk of developing T2D in the Spanish general population. These data support the importance of analyzing the influence of OS levels and their interaction with genetic variations in order to establish their real impact in T2D risk. Further studies are required to identify the real relevance of interactions between genetic variations and OS levels and the mechanisms involved in them.

Gene-environment interaction analysis of redox-related metals and genetic variants with plasma metabolic patterns in a general population from Spain: The Hortega Study

Background

Limited studies have evaluated the joint influence of redox-related metals and genetic variation on metabolic pathways. We analyzed the association of 11 metals with metabolic patterns, and the interacting role of candidate genetic variants, in 1145 participants from the Hortega Study, a population-based sample from Spain.

Methods

Urine antimony (Sb), arsenic, barium (Ba), cadmium (Cd), chromium (Cr), cobalt (Co), molybdenum (Mo) and vanadium (V), and plasma copper (Cu), selenium (Se) and zinc (Zn) were measured by ICP-MS and AAS, respectively. We summarized 54 plasma metabolites, measured with targeted NMR, by estimating metabolic principal components (mPC). Redox-related SNPs (N = 291) were measured by oligo-ligation assay.

Results

In our study, the association with metabolic principal component (mPC) 1 (reflecting non-essential and essential amino acids, including branched chain, and bacterial co-metabolism versus fatty acids and VLDL subclasses) was positive for Se and Zn, but inverse for Cu, arsenobetaine-corrected arsenic (As) and Sb. The association with mPC2 (reflecting essential amino acids, including aromatic, and bacterial co-metabolism) was inverse for Se, Zn and Cd. The association with mPC3 (reflecting LDL subclasses) was positive for Cu, Se and Zn, but inverse for Co. The association for mPC4 (reflecting HDL subclasses) was positive for Sb, but inverse for plasma Zn. These associations were mainly driven by Cu and Sb for mPC1; Se, Zn and Cd for mPC2; Co, Se and Zn for mPC3; and Zn for mPC4. The most SNP-metal interacting genes were NOX1, GSR, GCLC, AGT and REN. Co and Zn showed the highest number of interactions with genetic variants associated to enriched endocrine, cardiovascular and neurological pathways.

Conclusions

Exposures to Co, Cu, Se, Zn, As, Cd and Sb were associated with several metabolic patterns involved in chronic disease. Carriers of redox-related variants may have differential susceptibility to metabolic alterations associated to excessive exposure to metals.

•

In a population-based sample, cobalt, copper, selenium, zinc, arsenic, cadmium and antimony exposures were related to some metabolic patterns.

•

Carriers of redox-related variants displayed differential susceptibility to metabolic alterations associated to excessive metal exposures.

•

Cobalt and zinc showed a number of statistical interactions with variants from genes sharing biological pathways with a role in chronic diseases.

•

The metabolic impact of metals combined with variation in redox-related genes might be large in the population, given metals widespread exposure.

Association between genetic variants in oxidative stress-related genes and osteoporotic bone fracture. The Hortega follow-up study

The most widely accepted etiopathogenesis hypothesis of the origin of osteoporosis and its complications is that they are a consequence of bone aging and other environmental factors, together with a genetic predisposition. Evidence suggests that oxidative stress is crucial in bone pathologies associated with aging. The aim of this study was to determine whether genetic variants in oxidative stress-related genes modified the risk of osteoporotic fracture. We analysed 221 patients and 354 controls from the HORTEGA sample after 12-14 years of follow up. We studied the genotypic and allelic distribution of 53 SNPs in 24 genes involved in oxidative stress. The results showed that being a carrier of the variant allele of the SNP rs4077561 within TXNRD1 was the principal genetic risk factor associated with osteoporotic fracture and that variant allele of the rs1805754 M6PR, rs4964779 TXNRD1, rs406113 GPX6, rs2281082 TXN2 and rs974334 GPX6 polymorphisms are important genetic risk factors for fracture. This study provides information on the genetic factors associated with oxidative stress which are involved in the risk of osteoporotic fracture and reinforces the hypothesis that genetic factors are crucial in the etiopathogenesis of osteoporosis and its complications.

NF-kB2 Genetic Variations are Significantly Associated with Non-Small Cell Lung Cancer Risk and Overall Survival

During the last decade, a growing number of publications implicate NF-kB2 in NSCLC pathogenesis. Here, we investigated the clinical relevance of NF-kB2 single nucleotide polymorphisms (SNPs) rs7897947, rs11574852 and rs12769316 in NSCLC and their association with NF-kB2 protein and mRNA levels. Our data show that TT (rs7897947T >G) and AA (rs12769316G >A) genotypes were strongly associated with an increased risk for NSCLC (P = 0.019 and P = 0.003, respectively). Additionally, in multivariate analysis, TT (rs7897947T >G) homozygosity was associated with worse 2- and 3-year survival rates (P = 0.030 and P = 0.028, respectively), especially among patients with stages III/IV, who had worse 2, 3 and 5-year survival (P = 0.001, P = 0.022 and P = 0.035, respectively). In chemotherapy-treated patients, TT (rs12769316G >A) homozygosity was also associated with worse 2- and 3-year survival compared to G allele carriers (P = 0.006 and P = 0.014, respectively). Furthermore, rs12769316 was correlated with survival outcome of stage I and II patients (P = 0.031 and P = 0.006, respectively). Interestingly, amongst the patients who developed metastases, A allele carriers had better 5-year survival (P = 0.020). In addition, rs12769316 was associated with NF-kB2 protein (P = 0.001) and mRNA expression (P = 0.017) as well as with tumor maximum diameter (P = 0.025). Overall, this study suggests that rs7897947 and rs12769316 are involved in NSCLC susceptibility, in treatment response and in clinical outcome.

Are IL18RAP gene polymorphisms associated with body mass regulation? A cross-sectional study

OBJECTIVES

To investigate the association between IL18RAP and body mass index (BMI) and obesity and to verify the effect of a polymorphism in the microRNA136 (MIR136) IL18RAP binding region.

DESIGN

We analysed samples from two Spanish cross-sectional studies, VALCAR (Spanish Mediterranean coast) and Hortega (Spanish centre). These studies aimed at analysing cardiovascular risk and development of cardiovascular disease in the general population. Both populations correspond to regions with different characteristics.

SETTING

Five IL18RAP single nucleotide polymorphisms were selected using the SYSNPs web tool and analysed by oligonucleotide ligation assay (SNPlex). For the MIR136 functional study, cells were transfected with plasmids containing different rs7559479 polymorphism alleles and analysed by luciferase reporter assays.

PARTICIPANTS

1970 individuals (Caucasian, both genders): VALCAR (468) and Hortega (1502).

RESULTS

rs2293225, rs2272127 and rs7559479 showed the following associations: rs7559479 G allele correlated with a higher obesity risk (P=0.01; OR=1.82; 95% CI 1.15 to 2.87 for the VALCAR group; P=0.033; OR=1.35; 95% CI 1.03 to 1.79 for the Hortega population) and higher body mass index (BMI) values (P=0.0045; P=0.1 for VALCAR and Hortega, respectively); a significant association with obesity (P=0.0024, OR=1.44, 95% CI 1.14 to 1.82) and increased BMI values (P=0.008) was found when considering both populations together. rs2293225 T allele was associated with lower obesity risk (P=0.036; OR=0.60; 95% CI 0.35 to 0.96) and lower BMI values (P=0.0038; OR=1.41) while the rs2272127 G allele was associated with lower obesity risk (P=0.028; OR=0.66; 95% CI 0.44 to 0.97) only in the VALCAR population. A reporter assay showed that the presence of the A allele in rs7559479 was associated with increased MIR136 binding to IL18RAP.

CONCLUSIONS

Our results suggest that polymorphisms in IL18RAP influence susceptibility to obesity. We demonstrated that the A allele in rs7559479 increases MIR136 binding, which regulates IL-18 system activity.

Urine cadmium levels and albuminuria in a general population from Spain: A gene-environment interaction analysis

BACKGROUND

The interaction of cadmium with genes involved in oxidative stress, cadmium metabolism and transport pathways on albuminuria can provide biological insight on the relationship between cadmium and albuminuria at low exposure levels.

OBJECTIVES

We tested the hypothesis that specific genotypes in candidate genes may confer increased susceptibility to cadmium exposure.

METHODS

Cadmium exposure was estimated by inductively coupled plasma mass spectrometry (ICPMS) in urine from 1397 men and women aged 18-85years participating in the Hortega Study, a representative sample of a general population from Spain. Urine albumin was measured by automated nephelometric immunochemistry. Abnormal albuminuria was defined as urine albumin greater than or equal to 30mg/g.

RESULTS

The weighted prevalence of abnormal albuminuria was 6.3%. The median level of urine cadmium was 0.39 (IQR, 0.23-0.65) μg/g creatinine. Multivariable-adjusted geometric mean ratios of albuminuria comparing the two highest to the lowest tertile of urine cadmium were 1.62 (95% CI, 1.43-1.84) and 2.94 (95% CI, 2.58-3.35), respectively. The corresponding odds ratios of abnormal albuminuria were 1.58 (0.83, 3.02) and 4.54 (2.58, 8.00). The association between urine cadmium and albuminuria was observed across all participant subgroups evaluated including participants without hypertension, diabetes or chronic kidney disease. We observed Bonferroni-corrected statistically significant interactions between urine cadmium levels and polymorphisms in gene SLC30A7 and RAC1.

CONCLUSIONS

Increasing urine cadmium concentrations were cross-sectionally associated with increased albuminuria in a representative sample of a general population from Spain. Genetic variation in oxidative stress and cadmium metabolism and transport genes may confer differential susceptibility to potential cadmium effects.

A gene-environment interaction analysis of plasma selenium with prevalent and incident diabetes: The Hortega study

Background

Selenium and single-nucleotide-polymorphisms in selenoprotein genes have been associated to diabetes. However, the interaction of selenium with genetic variation in diabetes and oxidative stress-related genes has not been evaluated as a potential determinant of diabetes risk.

Methods

We evaluated the cross-sectional and prospective associations of plasma selenium concentrations with type 2 diabetes, and the interaction of selenium concentrations with genetic variation in candidate polymorphisms, in a representative sample of 1452 men and women aged 18–85 years from Spain.

Results

The geometric mean of plasma selenium levels in the study sample was 84.2 µg/L. 120 participants had diabetes at baseline. Among diabetes-free participants who were not lost during the follow-up (N=1234), 75 developed diabetes over time. The multivariable adjusted odds ratios (95% confidence interval) for diabetes prevalence comparing the second and third to the first tertiles of plasma selenium levels were 1.80 (1.03, 3.14) and 1.97 (1.14, 3.41), respectively. The corresponding hazard ratios (95% CI) for diabetes incidence were 1.76 (0.96, 3.22) and 1.80 (0.98, 3.31), respectively. In addition, we observed significant interactions between selenium and polymorphisms in PPARGC1A, and in genes encoding mitochondrial proteins, such as BCS1L and SDHA, and suggestive interactions of selenium with other genes related to selenoproteins and redox metabolism.

Conclusions

Plasma selenium was positively associated with prevalent and incident diabetes. While the statistical interactions of selenium with polymorphisms involved in regulation of redox and insulin signaling pathways provide biological plausibility to the positive associations of selenium with diabetes, further research is needed to elucidate the causal pathways underlying these associations.

Inflammatory-Related Genetic Variants in Non-Muscle-Invasive Bladder Cancer Prognosis: A Multimarker Bayesian Assessment

BACKGROUND

Increasing evidence points to the role of tumor immunologic environment on urothelial bladder cancer prognosis. This effect might be partly dependent on the host genetic context. We evaluated the association of SNPs in inflammation-related genes with non-muscle-invasive bladder cancer (NMIBC) risk-of-recurrence and risk-of-progression.

METHODS

We considered 822 NMIBC included in the SBC/EPICURO Study followed-up >10 years. We selected 1,679 SNPs belonging to 251 inflammatory genes. The association of SNPs with risk-of-recurrence and risk-of-progression was assessed using Cox regression single-marker (SMM) and multimarker methods (MMM) Bayes A and Bayesian LASSO. Discriminative abilities of the models were calculated using the c index and validated with bootstrap cross-validation procedures.

RESULTS

While no SNP was found to be associated with risk-of-recurrence using SMM, three SNPs in TNIP1, CD5, and JAK3 showed very strong association with posterior probabilities >90% using MMM. Regarding risk-of-progression, one SNP in CD3G was significantly associated using SMM (HR, 2.69; P = 1.55 × 10(-5)) and two SNPs in MASP1 and AIRE, showed a posterior probability ≥80% with MMM. Validated discriminative abilities of the models without and with the SNPs were 58.4% versus 60.5% and 72.1% versus 72.8% for risk-of-recurrence and risk-of-progression, respectively.

CONCLUSIONS

Using innovative analytic approaches, we demonstrated that SNPs in inflammatory-related genes were associated with NMIBC prognosis and that they improve the discriminative ability of prognostic clinical models for NMIBC.

IMPACT

This study provides proof of concept for the joint effect of genetic variants in improving the discriminative ability of clinical prognostic models. The approach may be extended to other diseases. Cancer Epidemiol Biomarkers Prev; 25(7); 1144-50. ©2016 AACR.

Pharmacogenomic study in patients with multiple sclerosis: Responders and nonresponders to IFN-β

Objectives:

We aimed to investigate the association between polymorphisms located in type I interferon (IFN)-induced genes, genes belonging to the toll-like receptor (TLR) pathway, and genes encoding neurotransmitter receptors and the response to IFN-β treatment in patients with multiple sclerosis (MS).

Methods:

In a first or screening phase of the study, 384 polymorphisms were genotyped in 830 patients with MS classified into IFN-β responders (n = 416) and nonresponders (n = 414) according to clinical criteria. In a second or validation phase, the most significant polymorphisms associated with IFN-β response were genotyped in an independent validation cohort of 555 patients with MS (281 IFN-β responders and 274 nonresponders).

Results:

Seven single nucleotide polymorphisms (SNPs) were selected from the screening phase for further validation: rs832032 (GABRR3; p = 0.0006), rs6597 (STUB1; p = 0.019), rs3747517 (IFIH1; p = 0.010), rs2277302 (PELI3; p = 0.017), rs10958713 (IKBKB; p = 0.003), rs2834202 (IFNAR1; p = 0.030), and rs4422395 (CXCL1; p = 0.017). None of these SNPs were significantly associated with IFN-β response when genotyped in an independent cohort of patients. Combined analysis of these SNPs in all patients with MS (N = 1,385) revealed 2 polymorphisms associated with IFN-β response: rs2277302 (PELI3; p = 0.008) and rs832032 (GABRR3; p = 0.006).

Conclusions:

These findings do not support an association between polymorphisms located in genes related to the type I IFN or TLR pathways or genes encoding neurotransmitter receptors and the clinical response to IFN-β. Nevertheless, additional genetic and functional studies of PELI3 and GABRR3 are warranted.

Variant of BCL3 gene is strongly associated with five-year survival of non-small-cell lung cancer patients

OBJECTIVES

BCL3, a known atypical IκB family member, has been documented to be upregulated in hematological malignancies and in some solid tumors, functioning as a crucial player in tumor development. Recently, rs8100239, a tag-Single Nucleotide Polymorphism (SNP) in BCL3 (T>A) has been identified, but there are no data regarding its involvement in non-small-cell lung cancer (NSCLC) initiation and progression.

MATERIALS AND METHODS

To study the possible association of BCL3 with NSCLC, 268 patients and 279 healthy controls were genotyped for rs8100239. Moreover, BCL3 protein expression was also investigated in 112 NSCLC cases through an immunohistochemical analysis.

RESULTS

NSCLC patients with AA genotype displayed significantly worse prognosis compared to T allele carriers (P<0.001), who had less frequent intermediate nuclear BCL3 expression (P=0.042). In addition, overexpression of BCL3 was detected in tumor specimens, compared to normal tissue (P<0.001). Furthermore, BCL3 protein levels were associated with five-year survival (P=0.039), maximum diameter of lesion (P=0.012), grade (P=0.002) and relapse frequency (P=0.041).

CONCLUSIONS

The present study is the first to show a relationship between the genetic variation rs8100239 of BCL3 and cancer patients' survival. It also represents the first quantitative evaluation of BCL3 expression in NSCLC. Our findings indicate that rs8100239 may be considered as a novel prognostic indicator, demonstrating also the overexpression of BCL3 protein in NSCLC and implicating this pivotal molecule in the pathogenesis of NSCLC.

Polymorphisms in endothelin system genes, arsenic levels and obesity risk

BACKGROUND/OBJECTIVES

Obesity has been linked to morbidity and mortality through increased risk for many chronic diseases. Endothelin (EDN) system has been related to endothelial function but it can be involved in lipid metabolism regulation: Receptor type A (EDNRA) activates lipolysis in adipocytes, the two endothelin receptors mediate arsenic-stimulated adipocyte dysfunction, and endothelin system can regulate adiposity by modulating adiponectin activity in different situations and, therefore, influence obesity development. The aim of the present study was to analyze if single nucleotide polymorphisms (SNPs) in the EDN system could be associated with human obesity.

SUBJECTS/METHODS

We analyzed two samples of general-population-based studies from two different regions of Spain: the VALCAR Study, 468 subjects from the area of Valencia, and the Hortega Study, 1502 subjects from the area of Valladolid. Eighteen SNPs throughout five genes were analyzed using SNPlex.

RESULTS

We found associations for two polymorphisms of the EDNRB gene which codifies for EDN receptor type B. Genotypes AG and AA of the rs5351 were associated with a lower risk for obesity in the VALCAR sample (p=0.048, OR=0.63) and in the Hortega sample (p=0.001, OR=0.62). Moreover, in the rs3759475 polymorphism, genotypes CT and TT were also associated with lower risk for obesity in the Hortega sample (p=0.0037, OR=0.66) and in the VALCAR sample we found the same tendency (p=0.12, OR=0.70). Furthermore, upon studying the pooled population, we found a stronger association with obesity (p=0.0001, OR=0.61 and p=0.0008, OR=0.66 for rs5351 and rs3759475, respectively). Regarding plasma arsenic levels, we have found a positive association for the two SNPs studied with obesity risk in individuals with higher arsenic levels in plasma: rs5351 (p=0.0054, OR=0.51) and rs3759475 (p=0.009, OR=0.53).

CONCLUSIONS

Our results support the hypothesis that polymorphisms of the EDNRB gene may influence the susceptibility to obesity and can interact with plasma arsenic levels.

Plasma selenium levels and oxidative stress biomarkers: a gene-environment interaction population-based study

The role of selenium exposure in preventing chronic disease is controversial, especially in selenium-repleted populations. At high concentrations, selenium exposure may increase oxidative stress. Studies evaluating the interaction of genetic variation in genes involved in oxidative stress pathways and selenium are scarce. We evaluated the cross-sectional association of plasma selenium concentrations with oxidative stress levels, measured as oxidized to reduced glutathione ratio (GSSG/GSH), malondialdehyde (MDA), and 8-oxo-7,8-dihydroguanine (8-oxo-dG) in urine, and the interacting role of genetic variation in oxidative stress candidate genes, in a representative sample of 1445 men and women aged 18-85 years from Spain. The geometric mean of plasma selenium levels in the study sample was 84.76 µg/L. In fully adjusted models the geometric mean ratios for oxidative stress biomarker levels comparing the highest to the lowest quintiles of plasma selenium levels were 0.61 (0.50-0.76) for GSSG/GSH, 0.89 (0.79-1.00) for MDA, and 1.06 (0.96-1.18) for 8-oxo-dG. We observed nonlinear dose-responses of selenium exposure and oxidative stress biomarkers, with plasma selenium concentrations above ~110 μg/L being positively associated with 8-oxo-dG, but inversely associated with GSSG/GSH and MDA. In addition, we identified potential risk genotypes associated with increased levels of oxidative stress markers with high selenium levels. Our findings support that high selenium levels increase oxidative stress in some biological processes. More studies are needed to disentangle the complexity of selenium biology and the relevance of potential gene-selenium interactions in relation to health outcomes in human populations.

Genomic and metabolomic profile associated to microalbuminuria

To identify factors related with the risk to develop microalbuminuria using combined genomic and metabolomic values from a general population study. One thousand five hundred and two subjects, Caucasian, more than 18 years, representative of the general population, were included. Blood pressure measurement and albumin/creatinine ratio were measured in a urine sample. Using SNPlex, 1251 SNPs potentially associated to urinary albumin excretion (UAE) were analyzed. Serum metabolomic profile was assessed by 1H NMR spectra using a Brucker Advance DRX 600 spectrometer. From the total population, 1217 (mean age 54 ± 19, 50.6% men, ACR>30 mg/g in 81 subjects) with high genotyping call rate were analysed. A characteristic metabolomic profile, which included products from mitochondrial and extra mitochondrial metabolism as well as branched amino acids and their derivative signals, were observed in microalbuminuric as compare to normoalbuminuric subjects. The comparison of the metabolomic profile between subjects with different UAE status for each of the genotypes associated to microalbuminuria revealed two SNPs, the rs10492025_TT of RPH3A gene and the rs4359_CC of ACE gene, with minimal or no statistically significant differences. Subjects with and without microalbuminuria, who shared the same genotype and metabolomic profile, differed in age. Microalbuminurics with the CC genotype of the rs4359 polymorphism and with the TT genotype of the rs10492025 polymorphism were seven years older and seventeen years younger, respectively as compared to the whole microalbuminuric subjects. With the same metabolomic environment, characteristic of subjects with microalbuminuria, the TT genotype of the rs10492025 polymorphism seems to increase and the CC genotype of the rs4359 polymorphism seems to reduce risk to develop microalbuminuria.

Genetic variation in the TP53 pathway and bladder cancer risk. a comprehensive analysis

INTRODUCTION

Germline variants in TP63 have been consistently associated with several tumors, including bladder cancer, indicating the importance of TP53 pathway in cancer genetic susceptibility. However, variants in other related genes, including TP53 rs1042522 (Arg72Pro), still present controversial results. We carried out an in depth assessment of associations between common germline variants in the TP53 pathway and bladder cancer risk.

MATERIAL AND METHODS

We investigated 184 tagSNPs from 18 genes in 1,058 cases and 1,138 controls from the Spanish Bladder Cancer/EPICURO Study. Cases were newly-diagnosed bladder cancer patients during 1998-2001. Hospital controls were age-gender, and area matched to cases. SNPs were genotyped in blood DNA using Illumina Golden Gate and TaqMan assays. Cases were subphenotyped according to stage/grade and tumor p53 expression. We applied classical tests to assess individual SNP associations and the Least Absolute Shrinkage and Selection Operator (LASSO)-penalized logistic regression analysis to assess multiple SNPs simultaneously.

RESULTS

Based on classical analyses, SNPs in BAK1 (1), IGF1R (5), P53AIP1 (1), PMAIP1 (2), SERINPB5 (3), TP63 (3), and TP73 (1) showed significant associations at p-value≤0.05. However, no evidence of association, either with overall risk or with specific disease subtypes, was observed after correction for multiple testing (p-value≥0.8). LASSO selected the SNP rs6567355 in SERPINB5 with 83% of reproducibility. This SNP provided an OR = 1.21, 95%CI 1.05-1.38, p-value = 0.006, and a corrected p-value = 0.5 when controlling for over-estimation.

DISCUSSION

We found no strong evidence that common variants in the TP53 pathway are associated with bladder cancer susceptibility. Our study suggests that it is unlikely that TP53 Arg72Pro is implicated in the UCB in white Europeans. SERPINB5 and TP63 variation deserve further exploration in extended studies.

Application of multi-SNP approaches Bayesian LASSO and AUC-RF to detect main effects of inflammatory-gene variants associated with bladder cancer risk

The relationship between inflammation and cancer is well established in several tumor types, including bladder cancer. We performed an association study between 886 inflammatory-gene variants and bladder cancer risk in 1,047 cases and 988 controls from the Spanish Bladder Cancer (SBC)/EPICURO Study. A preliminary exploration with the widely used univariate logistic regression approach did not identify any significant SNP after correcting for multiple testing. We further applied two more comprehensive methods to capture the complexity of bladder cancer genetic susceptibility: Bayesian Threshold LASSO (BTL), a regularized regression method, and AUC-Random Forest, a machine-learning algorithm. Both approaches explore the joint effect of markers. BTL analysis identified a signature of 37 SNPs in 34 genes showing an association with bladder cancer. AUC-RF detected an optimal predictive subset of 56 SNPs. 13 SNPs were identified by both methods in the total population. Using resources from the Texas Bladder Cancer study we were able to replicate 30% of the SNPs assessed. The associations between inflammatory SNPs and bladder cancer were reexamined among non-smokers to eliminate the effect of tobacco, one of the strongest and most prevalent environmental risk factor for this tumor. A 9 SNP-signature was detected by BTL. Here we report, for the first time, a set of SNP in inflammatory genes jointly associated with bladder cancer risk. These results highlight the importance of the complex structure of genetic susceptibility associated with cancer risk.

Recent human evolution has shaped geographical differences in susceptibility to disease

Background

Searching for associations between genetic variants and complex diseases has been a very active area of research for over two decades. More than 51,000 potential associations have been studied and published, a figure that keeps increasing, especially with the recent explosion of array-based Genome-Wide Association Studies. Even if the number of true associations described so far is high, many of the putative risk variants detected so far have failed to be consistently replicated and are widely considered false positives. Here, we focus on the world-wide patterns of replicability of published association studies.

Results

We report three main findings. First, contrary to previous results, genes associated to complex diseases present lower degrees of genetic differentiation among human populations than average genome-wide levels. Second, also contrary to previous results, the differences in replicability of disease associated-loci between Europeans and East Asians are highly correlated with genetic differentiation between these populations. Finally, highly replicated genes present increased levels of high-frequency derived alleles in European and Asian populations when compared to African populations.

Conclusions

Our findings highlight the heterogeneous nature of the genetic etiology of complex disease, confirm the importance of the recent evolutionary history of our species in current patterns of disease susceptibility and could cast doubts on the status as false positives of some associations that have failed to replicate across populations.