A Genome-Wide Association Study of Upper Aerodigestive Tract Cancers Conducted within the INHANCE Consortium

Genome-wide association studies (GWAS) have been successful in identifying common genetic variation involved in susceptibility to etiologically complex disease. We conducted a GWAS to identify common genetic variation involved in susceptibility to upper aero-digestive tract (UADT) cancers. Genome-wide genotyping was carried out using the Illumina HumanHap300 beadchips in 2,091 UADT cancer cases and 3,513 controls from two large European multi-centre UADT cancer studies, as well as 4,821 generic controls. The 19 top-ranked variants were investigated further in an additional 6,514 UADT cancer cases and 7,892 controls of European descent from an additional 13 UADT cancer studies participating in the INHANCE consortium. Five common variants presented evidence for significant association in the combined analysis (p≤5×10−7). Two novel variants were identified, a 4q21 variant (rs1494961, p = 1×10−8) located near DNA repair related genes HEL308 and FAM175A (or Abraxas) and a 12q24 variant (rs4767364, p = 2×10−8) located in an extended linkage disequilibrium region that contains multiple genes including the aldehyde dehydrogenase 2 (ALDH2) gene. Three remaining variants are located in the ADH gene cluster and were identified previously in a candidate gene study involving some of these samples. The association between these three variants and UADT cancers was independently replicated in 5,092 UADT cancer cases and 6,794 controls non-overlapping samples presented here (rs1573496-ADH7, p = 5×10−8; rs1229984-ADH1B, p = 7×10−9; and rs698-ADH1C, p = 0.02). These results implicate two variants at 4q21 and 12q24 and further highlight three ADH variants in UADT cancer susceptibility.

Association study of genetic variants of pro-inflammatory chemokine and cytokine genes in systemic lupus erythematosus

BACKGROUND. Several lines of evidence suggest that chemokines and cytokines play an important role in the inflammatory development and progression of systemic lupus erythematosus. The aim of this study was to evaluate the relevance of functional genetic variations of RANTES, IL-8, IL-1alpha, and MCP-1 for systemic lupus erythematosus. METHODS. The study was conducted on 500 SLE patients and 481 ethnically matched healthy controls. Genotyping of polymorphisms in the RANTES, IL-8, IL-1alpha, and MCP-1 genes were performed using a real-time polymerase chain reaction (PCR) system with pre-developed TaqMan allelic discrimination assay. RESULTS. No significant differences between SLE patients and healthy controls were observed when comparing genotype, allele or haplotype frequencies of the RANTES, IL-8, IL-1alpha, and MCP-1 polymorphisms. In addition, no evidence for association with clinical sub-features of SLE was found. CONCLUSION. These results suggest that the tested functional variation of RANTES, IL-8, IL-1alpha, and MCP-1 genes do not confer a relevant role in the susceptibility or severity of SLE in the Spanish population.

Genome-wide association study of multiple sclerosis confirms a novel locus at 5p13.1.

Multiple Sclerosis (MS) is the most common progressive and disabling neurological condition affecting young adults in the world today. From a genetic point of view, MS is a complex disorder resulting from the combination of genetic and non-genetic factors. We aimed to identify previously unidentified loci conducting a new GWAS of Multiple Sclerosis (MS) in a sample of 296 MS cases and 801 controls from the Spanish population. Meta-analysis of our data in combination with previous GWAS was done. A total of 17 GWAS-significant SNPs, corresponding to three different loci were identified:HLA, IL2RA, and 5p13.1. All three have been previously reported as GWAS-significant. We confirmed our observation in 5p13.1 for rs9292777 using two additional independent Spanish samples to make a total of 4912 MS cases and 7498 controls (ORpooled = 0.84; 95%CI: 0.80-0.89; p = 1.36 × 10-9). This SNP differs from the one reported within this locus in a recent GWAS. Although it is unclear whether both signals are tapping the same genetic association, it seems clear that this locus plays an important role in the pathogenesis of MS.

IL2RA/CD25 gene polymorphisms: uneven association with multiple sclerosis (MS) and type 1 diabetes (T1D).

BACKGROUND: IL-2 receptor (IL2R) alpha is the specific component of the high affinity IL2R system involved in the immune response and in the control of autoimmunity. METHODS AND RESULTS: Here we perform a replication and fine mapping of the IL2RA gene region analyzing 3 SNPs previously associated with multiple sclerosis (MS) and 5 SNPs associated with type 1 diabetes (T1D) in a collection of 798 MS patients and 927 matched Caucasian controls from the south of Spain. We observed association with MS in 6 of 8 SNPs. The rs1570538, at the 3'- UTR extreme of the gene, previously reported to have a weak association with MS, is replicated here (P = 0.032). The most associated T1D SNP (rs41295061) was not associated with MS in the present study. However, the rs35285258, belonging to another independent group of SNPs associated with T1D, showed the maximal association in this study but different risk allele. We replicated the association of only one (rs2104286) of the two IL2RA SNPs identified in the recently performed genome-wide association study of MS. CONCLUSIONS: These findings confirm and extend the association of this gene with MS and reveal a genetic heterogeneity of the associated polymorphisms and risk alleles between MS and T1D suggesting different immunopathological roles of IL2RA in these two diseases.

HLA and non-HLA genes in Behçet's disease: a multicentric study in the Spanish population.

INTRODUCTION According to genome wide association (GWA) studies as well as candidate gene approaches, Behçet's disease (BD) is associated with human leukocyte antigen (HLA)-A and HLA-B gene regions. The HLA-B51 has been consistently associated with the disease, but the role of other HLA class I molecules remains controversial. Recently, variants in non-HLA genes have also been associated with BD. The aims of this study were to further investigate the influence of the HLA region in BD and to explore the relationship with non-HLA genes recently described to be associated in other populations. METHODS This study included 304 BD patients and 313 ethnically matched controls. HLA-A and HLA-B low resolution typing was carried out by PCR-SSOP Luminex. Eleven tag single nucleotide polymorphisms (SNPs) located outside of the HLA-region, previously described associated with the disease in GWA studies and having a minor allele frequency in Caucasians greater than 0.15 were genotyped using TaqMan assays. Phenotypic and genotypic frequencies were estimated by direct counting and distributions were compared using the χ(2) test. RESULTS In addition to HLA-B*51, HLA-B*57 was found as a risk factor in BD, whereas, B*35 was found to be protective. Other HLA-A and B specificities were suggestive of association with the disease as risk (A*02 and A*24) or protective factors (A*03 and B*58). Regarding the non-HLA genes, the three SNPs located in IL23R and one of the SNPs in IL10 were found to be significantly associated with susceptibility to BD in our population. CONCLUSION Different HLA specificities are associated with Behçet's disease in addition to B*51. Other non-HLA genes, such as IL23R and IL-10, play a role in the susceptibility to the disease.

The INSIG2 rs7566605 polymorphism is not associated with body mass index and breast cancer risk

Background The single nucleotide polymorphism rs7566605, located in the promoter of the INSIG2 gene, has been the subject of a strong scientific effort aimed to elucidate its possible association with body mass index (BMI). The first report showing that rs7566605 could be associated with body fatness was a genome-wide association study (GWAS) which used BMI as the primary phenotype. Many follow-up studies sought to validate the association of rs7566605 with various markers of obesity, with several publications reporting inconsistent findings. BMI is considered to be one of the measures of choice to evaluate body fatness and there is evidence that body fatness is related with an increased risk of breast cancer (BC). Methods we tested in a large-scale association study (3,973 women, including 1,269 invasive BC cases and 2,194 controls), nested within the EPIC cohort, the involvement of rs7566605 as predictor of BMI and BC risk. Results and Conclusions In this study we were not able to find any statistically significant association between this SNP and BMI, nor did we find any significant association between the SNP and an increased risk of breast cancer overall and by subgroups of age, or menopausal status.

Multiple sclerosis risk variant HLA-DRB1*1501 associates with high expression of DRB1 gene in different human populations.

The human leukocyte antigen (HLA) DRB1*1501 has been consistently associated with multiple sclerosis (MS) in nearly all populations tested. This points to a specific antigen presentation as the pathogenic mechanism though this does not fully explain the disease association. The identification of expression quantitative trait loci (eQTL) for genes in the HLA locus poses the question of the role of gene expression in MS susceptibility. We analyzed the eQTLs in the HLA region with respect to MS-associated HLA-variants obtained from genome-wide association studies (GWAS). We found that the Tag of DRB1*1501, rs3135388 A allele, correlated with high expression of DRB1, DRB5 and DQB1 genes in a Caucasian population. In quantitative terms, the MS-risk AA genotype carriers of rs3135388 were associated with 15.7-, 5.2- and 8.3-fold higher expression of DQB1, DRB5 and DRB1, respectively, than the non-risk GG carriers. The haplotype analysis of expression-associated variants in a Spanish MS cohort revealed that high expression of DRB1 and DQB1 alone did not contribute to the disease. However, in Caucasian, Asian and African American populations, the DRB1*1501 allele was always highly expressed. In other immune related diseases such as type 1 diabetes, inflammatory bowel disease, ulcerative colitis, asthma and IgA deficiency, the best GWAS-associated HLA SNPs were also eQTLs for different HLA Class II genes. Our data suggest that the DR/DQ expression levels, together with specific structural properties of alleles, seem to be the causal effect in MS and in other immunopathologies rather than specific antigen presentation alone.

Analysis of TNFAIP3, a feedback inhibitor of nuclear factor-κB and the neighbor intergenic 6q23 region in rheumatoid arthritis susceptibility

INTRODUCTION Genome-wide association studies of rheumatoid arthritis (RA) have identified an association of the disease with a 6q23 region devoid of genes. TNFAIP3, an RA candidate gene, flanks this region, and polymorphisms in both the TNFAIP3 gene and the intergenic region are associated with systemic lupus erythematosus. We hypothesized that there is a similar association with RA, including polymorphisms in TNFAIP3 and the intergenic region. METHODS To test this hypothesis, we selected tag-single nucleotide polymorphisms (SNPs) in both loci. They were analyzed in 1,651 patients with RA and 1,619 control individuals of Spanish ancestry. RESULTS Weak evidence of association was found both in the 6q23 intergenic region and in the TNFAIP3 locus. The rs582757 SNP and a common haplotype in the TNFAIP3 locus exhibited association with RA. In the intergenic region, two SNPs were associated, namely rs609438 and rs13207033. The latter was only associated in patients with anti-citrullinated peptide antibodies. Overall, statistical association was best explained by the interdependent contribution of SNPs from the two loci TNFAIP3 and the 6q23 intergenic region. CONCLUSIONS Our data are consistent with the hypothesis that several RA genetic factors exist in the 6q23 region, including polymorphisms in the TNFAIP3 gene, like that previously described for systemic lupus erythematosus.

Genetic variations in drug-induced liver injury (DILI): resolving the puzzle.

Despite stringent requirements for drug development imposed by regulatory agencies, drug-induced liver injury (DILI) is an increasing health problem and a significant cause for failure to approve drugs, market withdrawal of commercialized medications, and adoption of regulatory measures. The pathogenesis is yet undefined, though the rare occurrence of idiosyncratic DILI (1/100,000–1/10,000) and the fact that hepatotoxicity often recurs after re-exposure to the culprit drug under different environmental conditions strongly points toward a major role for genetic variations in the underlying mechanism and susceptibility. Pharmacogenetic studies in DILI have to a large extent focused on genes involved in drug metabolism, as polymorphisms in these genes may generate increased plasma drug concentrations as well as lower clearance rates when treated with standard medication doses. A range of studies have identified a number of genetic variants in drug metabolism Phase I, II, and III genes, including cytochrome P450 (CYP) 2E1, N-acetyltransferase 2, UDP-glucuronosyltransferase 2B7, glutathione S-transferase M1/T1, ABCB11, and ABCC2, that enhance DILI susceptibility (Andrade et al., 2009; Agundez et al., 2011). Several metabolic gene variants, such as CYP2E1c1 and NAT2 slow, have been associated with DILI induced by specific drugs based on individual drug metabolism information. Others, such as GSTM1 and T1 null alleles have been associated with enhanced risk of DILI development induced by a large range of drugs. Hence, these variants appear to have a more general role in DILI susceptibility due to their role in reducing the cell's antioxidative capacity (Lucena et al., 2008). Mitochondrial superoxide dismutase (SOD2) and glutathione peroxidase 1 (GPX1) are two additional enzymes involved in combating oxidative stress, with specific genetic variants shown to enhance the risk of developing DILI

Fine mapping and functional analysis of the multiple sclerosis risk gene CD6.

CD6 has recently been identified and validated as risk gene for multiple sclerosis (MS), based on the association of a single nucleotide polymorphism (SNP), rs17824933, located in intron 1. CD6 is a cell surface scavenger receptor involved in T-cell activation and proliferation, as well as in thymocyte differentiation. In this study, we performed a haptag SNP screen of the CD6 gene locus using a total of thirteen tagging SNPs, of which three were non-synonymous SNPs, and replicated the recently reported GWAS SNP rs650258 in a Spanish-Basque collection of 814 controls and 823 cases. Validation of the six most strongly associated SNPs was performed in an independent collection of 2265 MS patients and 2600 healthy controls. We identified association of haplotypes composed of two non-synonymous SNPs [rs11230563 (R225W) and rs2074225 (A257V)] in the 2(nd) SRCR domain with susceptibility to MS (P max(T) permutation = 1×10(-4)). The effect of these haplotypes on CD6 surface expression and cytokine secretion was also tested. The analysis showed significantly different CD6 expression patterns in the distinct cell subsets, i.e. - CD4(+) naïve cells, P = 0.0001; CD8(+) naïve cells, P<0.0001; CD4(+) and CD8(+) central memory cells, P = 0.01 and 0.05, respectively; and natural killer T (NKT) cells, P = 0.02; with the protective haplotype (RA) showing higher expression of CD6. However, no significant changes were observed in natural killer (NK) cells, effector memory and terminally differentiated effector memory T cells. Our findings reveal that this new MS-associated CD6 risk haplotype significantly modifies expression of CD6 on CD4(+) and CD8(+) T cells.

Impact of obesity-related genes in Spanish population

Background: The objective was to investigate the association between BMI and single nucleotide polymorphisms previously identified of obesity-related genes in two Spanish populations. Forty SNPs in 23 obesity-related genes were evaluated in a rural population characterized by a high prevalence of obesity (869 subjects, mean age 46 yr, 62% women, 36% obese) and in an urban population (1425 subjects, mean age 54 yr, 50% women, 19% obese). Genotyping was assessed by using SNPlex and PLINK for the association analysis. Results: Polymorphisms of the FTO were significantly associated with BMI, in the rural population (beta 0.87, p-value <0.001). None of the other SNPs showed significant association after Bonferroni correction in the two populations or in the pooled analysis. A weighted genetic risk score (wGRS) was constructed using the risk alleles of the Tag-SNPs with a positive Beta parameter in both populations. From the first to the fifth quintile of the score, the BMI increased 0.45 kg/m2 in Hortega and 2.0 kg/m2 in Pizarra. Overall, the obesity predictive value was low (less than 1%). Conclusion: The risk associated with polymorphisms is low and the overall effect on BMI or obesity prediction is minimal. A weighted genetic risk score based on genes mainly acting through central nervous system mechanisms was associated with BMI but it yields minimal clinical prediction for the obesity risk in the general population.

Evaluation of the IL2/IL21, IL2RA and IL2RB genetic variants influence on the endogenous non-anterior uveitis genetic predisposition.

BACKGROUND Recently, different genetic variants located within the IL2/IL21 genetic region as well as within both IL2RA and IL2RB loci have been associated to multiple autoimmune disorders. We aimed to investigate for the first time the potential influence of the IL2/IL21, IL2RA and IL2RB most associated polymorphisms with autoimmunity on the endogenous non-anterior uveitis genetic predisposition. METHODS A total of 196 patients with endogenous non-anterior uveitis and 760 healthy controls, all of them from Caucasian population, were included in the current study. The IL2/IL21 (rs2069762, rs6822844 and rs907715), IL2RA (2104286, rs11594656 and rs12722495) and IL2RB (rs743777) genetic variants were genotyped using TaqMan® allelic discrimination assays. RESULTS A statistically significant difference was found for the rs6822844 (IL2/IL21 region) minor allele frequency in the group of uveitis patients compared with controls (P(-value)=0.02, OR=0.64 CI 95%=0.43-0.94) although the significance was lost after multiple testing correction. Furthermore, no evidence of association with uveitis was detected for the analyzed genetic variants of the IL2RA or IL2RB loci. CONCLUSION Our results indicate that analyzed IL2/IL21, IL2RA and IL2RB polymorphisms do not seem to play a significant role on the non-anterior uveitis genetic predisposition although further studies are needed in order to clear up the influence of these loci on the non-anterior uveitis susceptibility.

Influence of the LILRA3 Deletion on Multiple Sclerosis Risk: Original Data and Meta-Analysis.

BACKGROUND Multiple sclerosis (MS) is a neurodegenerative, autoimmune disease of the central nervous system. Genome-wide association studies (GWAS) have identified over hundred polymorphisms with modest individual effects in MS susceptibility and they have confirmed the main individual effect of the Major Histocompatibility Complex. Additional risk loci with immunologically relevant genes were found significantly overrepresented. Nonetheless, it is accepted that most of the genetic architecture underlying susceptibility to the disease remains to be defined. Candidate association studies of the leukocyte immunoglobulin-like receptor LILRA3 gene in MS have been repeatedly reported with inconsistent results. OBJECTIVES In an attempt to shed some light on these controversial findings, a combined analysis was performed including the previously published datasets and three newly genotyped cohorts. Both wild-type and deleted LILRA3 alleles were discriminated in a single-tube PCR amplification and the resulting products were visualized by their different electrophoretic mobilities. RESULTS AND CONCLUSION Overall, this meta-analysis involved 3200 MS patients and 3069 matched healthy controls and it did not evidence significant association of the LILRA3 deletion [carriers of LILRA3 deletion: p = 0.25, OR (95% CI) = 1.07 (0.95-1.19)], even after stratification by gender and the HLA-DRB1*15:01 risk allele.

The association of education with body mass index and waist circumference in the EPIC-PANACEA study

Background To examine the association of education with body mass index (BMI) and waist circumference (WC) in the European Prospective Investigation into Cancer and Nutrition (EPIC). Method This study included 141,230 male and 336,637 female EPIC-participants, who were recruited between 1992 and 2000. Education, which was assessed by questionnaire, was classified into four categories; BMI and WC, measured by trained personnel in most participating centers, were modeled as continuous dependent variables. Associations were estimated using multilevel mixed effects linear regression models. Results Compared with the lowest education level, BMI and WC were significantly lower for all three higher education categories, which was consistent for all countries. Women with university degree had a 2.1 kg/m2 lower BMI compared with women with lowest education level. For men, a statistically significant, but less pronounced difference was observed (1.3 kg/m2). The association between WC and education level was also of greater magnitude for women: compared with the lowest education level, average WC of women was lower by 5.2 cm for women in the highest category. For men the difference was 2.9 cm. Conclusion In this European cohort, there is an inverse association between higher BMI as well as higher WC and lower education level. Public Health Programs that aim to reduce overweight and obesity should primarily focus on the lower educated population.

Association between pre-diagnostic circulating vitamin D concentration and risk of colorectal cancer in European populations:a nested case-control study

Objective. To examine the association between pre-diagnostic circulating vitamin D concentration, dietary intake of vitamin D and calcium, and the risk of colorectal cancer in European populations. Design Nested case-control study. Setting. The study was conducted within the EPIC study, a cohort of more than 520 000 participants from 10 western European countries. Participants: 1248 cases of incident colorectal cancer, which developed after enrolment into the cohort, were matched to 1248 controls. Main outcome measures. Circulating vitamin D concentration (25-hydroxy-vitamin-D, 25-(OH)D) was measured by enzyme immunoassay. Dietary and lifestyle data were obtained from questionnaires. Incidence rate ratios and 95% confidence intervals for the risk of colorectal cancer by 25-(OH)D concentration and levels of dietary calcium and vitamin D intake were estimated from multivariate conditional logistic regression models, with adjustment for potential dietary and other confounders. Results. 25-(OH)D concentration showed a strong inverse linear dose-response association with risk of colorectal cancer (P for trend <0.001). Compared with a pre-defined mid-level concentration of 25-(OH)D (50.0-75.0 nmol/l), lower levels were associated with higher colorectal cancer risk (<25.0 nmol/l: incidence rate ratio 1.32 (95% confidence interval 0.87 to 2.01); 25.0-49.9 nmol/l: 1.28 (1.05 to 1.56), and higher concentrations associated with lower risk (75.0-99.9 nmol/l: 0.88 (0.68 to 1.13); ≥100.0 nmol/l: 0.77 (0.56 to 1.06)). In analyses by quintile of 25-(OH)D concentration, patients in the highest quintile had a 40% lower risk of colorectal cancer than did those in the lowest quintile (P<0.001). Subgroup analyses showed a strong association for colon but not rectal cancer (P for heterogeneity=0.048). Greater dietary intake of calcium was associated with a lower colorectal cancer risk. Dietary vitamin D was not associated with disease risk. Findings did not vary by sex and were not altered by corrections for season or month of blood donation. Conclusions The results of this large observational study indicate a strong inverse association between levels of pre-diagnostic 25-(OH)D concentration and risk of colorectal cancer in western European populations. Further randomised trials are needed to assess whether increases in circulating 25-(OH)D concentration can effectively decrease the risk of colorectal cancer.