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.

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.

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.

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.

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.

CD40: Novel Association with Crohn's Disease and Replication in Multiple Sclerosis Susceptibility

Background: A functional polymorphism located at 21 from the start codon of the CD40 gene, rs1883832, was previously reported to disrupt a Kozak sequence essential for translation. It has been consistently associated with Graves’ disease risk in populations of different ethnicity and genetic proxies of this variant evaluated in genome-wide association studies have shown evidence of an effect in rheumatoid arthritis and multiple sclerosis (MS) susceptibility. However, the protective allele associated with Graves’ disease or rheumatoid arthritis has shown a risk role in MS, an effect that we aimed to replicate in the present work. We hypothesized that this functional polymorphism might also show an association with other complex autoimmune condition such as inflammatory bowel disease, given the CD40 overexpression previously observed in Crohn’s disease (CD) lesions. Methodology: Genotyping of rs1883832C.T was performed in 1564 MS, 1102 CD and 969 ulcerative colitis (UC) Spanish patients and in 2948 ethnically matched controls by TaqMan chemistry. Principal Findings: The observed effect of the minor allele rs1883832T was replicated in our independent Spanish MS cohort [p= 0.025; OR (95% CI)= 1.12 (1.01–1.23)]. The frequency of the minor allele was also significantly higher in CD patients than in controls [p= 0.002; OR (95% CI)= 1.19 (1.06–1.33)]. This increased predisposition was not detected in UC patients [p= 0.5; OR (95% CI)= 1.04 (0.93–1.17)]. Conclusion: The impact of CD40 rs1883832 on MS and CD risk points to a common signaling shared by these autoimmune conditions

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.

Breast cancer prognostic classification in the molecular era: the role of histological grade.

Breast cancer is a heterogeneous disease with varied morphological appearances, molecular features, behavior, and response to therapy. Current routine clinical management of breast cancer relies on the availability of robust clinical and pathological prognostic and predictive factors to support clinical and patient decision making in which potentially suitable treatment options are increasingly available. One of the best-established prognostic factors in breast cancer is histological grade, which represents the morphological assessment of tumor biological characteristics and has been shown to be able to generate important information related to the clinical behavior of breast cancers. Genome-wide microarray-based expression profiling studies have unraveled several characteristics of breast cancer biology and have provided further evidence that the biological features captured by histological grade are important in determining tumor behavior. Also, expression profiling studies have generated clinically useful data that have significantly improved our understanding of the biology of breast cancer, and these studies are undergoing evaluation as improved prognostic and predictive tools in clinical practice. Clinical acceptance of these molecular assays will require them to be more than expensive surrogates of established traditional factors such as histological grade. It is essential that they provide additional prognostic or predictive information above and beyond that offered by current parameters. Here, we present an analysis of the validity of histological grade as a prognostic factor and a consensus view on the significance of histological grade and its role in breast cancer classification and staging systems in this era of emerging clinical use of molecular classifiers.

Hepatic safety of antibiotics used in primary care

Antibiotics used by general practitioners frequently appear in adverse-event reports of drug-induced hepatotoxicity. Most cases are idiosyncratic (the adverse reaction cannot be predicted from the drug's pharmacological profile or from pre-clinical toxicology tests) and occur via an immunological reaction or in response to the presence of hepatotoxic metabolites. With the exception of trovafloxacin and telithromycin (now severely restricted), hepatotoxicity crude incidence remains globally low but variable. Thus, amoxicillin/clavulanate and co-trimoxazole, as well as flucloxacillin, cause hepatotoxic reactions at rates that make them visible in general practice (cases are often isolated, may have a delayed onset, sometimes appear only after cessation of therapy and can produce an array of hepatic lesions that mirror hepatobiliary disease, making causality often difficult to establish). Conversely, hepatotoxic reactions related to macrolides, tetracyclines and fluoroquinolones (in that order, from high to low) are much rarer, and are identifiable only through large-scale studies or worldwide pharmacovigilance reporting. For antibiotics specifically used for tuberculosis, adverse effects range from asymptomatic increases in liver enzymes to acute hepatitis and fulminant hepatic failure. Yet, it is difficult to single out individual drugs, as treatment always entails associations. Patients at risk are mainly those with previous experience of hepatotoxic reaction to antibiotics, the aged or those with impaired hepatic function in the absence of close monitoring, making it important to carefully balance potential risks with expected benefits in primary care. Pharmacogenetic testing using the new genome-wide association studies approach holds promise for better understanding the mechanism(s) underlying hepatotoxicity.

Insulin receptor substrate-1 expression is increased in circulating leukocytes of patients with acute coronary syndrome.

The mechanisms underlying the increased risk of cardiovascular disease associated with diabetes mellitus (DM) are not fully defined. Insulin resistance in human metabolic syndrome patients is associated with decreased expression of the insulin receptor substrate-2- (Irs2-) AKT2 axis in mononuclear leukocytes (MLs). Moreover, acute coronary syndrome (ACS) has been linked through genome-wide association studies to the 2q36-q37.3 locus, which contains the Irs1 gene. Here, we investigated the expression of insulin-signaling pathway genes in MLs from patients with DM, ACS, and ACS plus DM. Quantitative real-time PCR expression studies showed no differences in the mRNA levels of Irs2, Akt2, and Akt1 among all patients. However, Irs1 mRNA expression was significantly increased in patients with ACS-diabetics and nondiabetics-compared with diabetic patients without ACS (P < .02 and P < .005, resp.). The present study reveals for the first time an association between increased Irs1 mRNA levels in MLs of patients with ACS which is not related to DM.

Sources of pre-analytical variations in yield of DNA extracted from blood samples: analysis of 50,000 DNA samples in EPIC.

The European Prospective Investigation into Cancer and nutrition (EPIC) is a long-term, multi-centric prospective study in Europe investigating the relationships between cancer and nutrition. This study has served as a basis for a number of Genome-Wide Association Studies (GWAS) and other types of genetic analyses. Over a period of 5 years, 52,256 EPIC DNA samples have been extracted using an automated DNA extraction platform. Here we have evaluated the pre-analytical factors affecting DNA yield, including anthropometric, epidemiological and technical factors such as center of subject recruitment, age, gender, body-mass index, disease case or control status, tobacco consumption, number of aliquots of buffy coat used for DNA extraction, extraction machine or procedure, DNA quantification method, degree of haemolysis and variations in the timing of sample processing. We show that the largest significant variations in DNA yield were observed with degree of haemolysis and with center of subject recruitment. Age, gender, body-mass index, cancer case or control status and tobacco consumption also significantly impacted DNA yield. Feedback from laboratories which have analyzed DNA with different SNP genotyping technologies demonstrate that the vast majority of samples (approximately 88%) performed adequately in different types of assays. To our knowledge this study is the largest to date to evaluate the sources of pre-analytical variations in DNA extracted from peripheral leucocytes. The results provide a strong evidence-based rationale for standardized recommendations on blood collection and processing protocols for large-scale genetic studies.

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

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.

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.