Polymorphisms within Autophagy-Related Genes as Susceptibility Biomarkers for Multiple Myeloma: A Meta-Analysis of Three Large Cohorts and Functional Characterization.

Multiple myeloma (MM) arises following malignant proliferation of plasma cells in the bone marrow, that secrete high amounts of specific monoclonal immunoglobulins or light chains, resulting in the massive production of unfolded or misfolded proteins. Autophagy can have a dual role in tumorigenesis, by eliminating these abnormal proteins to avoid cancer development, but also ensuring MM cell survival and promoting resistance to treatments. To date no studies have determined the impact of genetic variation in autophagy-related genes on MM risk. We performed meta-analysis of germline genetic data on 234 autophagy-related genes from three independent study populations including 13,387 subjects of European ancestry (6863 MM patients and 6524 controls) and examined correlations of statistically significant single nucleotide polymorphisms (SNPs; p < 1 × 10-9) with immune responses in whole blood, peripheral blood mononuclear cells (PBMCs), and monocyte-derived macrophages (MDM) from a large population of healthy donors from the Human Functional Genomic Project (HFGP). We identified SNPs in six loci, CD46, IKBKE, PARK2, ULK4, ATG5, and CDKN2A associated with MM risk (p = 4.47 × 10-4-5.79 × 10-14). Mechanistically, we found that the ULK4rs6599175 SNP correlated with circulating concentrations of vitamin D3 (p = 4.0 × 10-4), whereas the IKBKErs17433804 SNP correlated with the number of transitional CD24+CD38+ B cells (p = 4.8 × 10-4) and circulating serum concentrations of Monocyte Chemoattractant Protein (MCP)-2 (p = 3.6 × 10-4). We also found that the CD46rs1142469 SNP correlated with numbers of CD19+ B cells, CD19+CD3- B cells, CD5+IgD- cells, IgM- cells, IgD-IgM- cells, and CD4-CD8- PBMCs (p = 4.9 × 10-4-8.6 × 10-4) and circulating concentrations of interleukin (IL)-20 (p = 0.00082). Finally, we observed that the CDKN2Ars2811710 SNP correlated with levels of CD4+EMCD45RO+CD27- cells (p = 9.3 × 10-4). These results suggest that genetic variants within these six loci influence MM risk through the modulation of specific subsets of immune cells, as well as vitamin D3-, MCP-2-, and IL20-dependent pathways.

Transcriptome-wide association study of multiple myeloma identifies candidate susceptibility genes.

BACKGROUND: While genome-wide association studies (GWAS) of multiple myeloma (MM) have identified variants at 23 regions influencing risk, the genes underlying these associations are largely unknown. To identify candidate causal genes at these regions and search for novel risk regions, we performed a multi-tissue transcriptome-wide association study (TWAS). RESULTS: GWAS data on 7319 MM cases and 234,385 controls was integrated with Genotype-Tissue Expression Project (GTEx) data assayed in 48 tissues (sample sizes, N = 80-491), including lymphocyte cell lines and whole blood, to predict gene expression. We identified 108 genes at 13 independent regions associated with MM risk, all of which were in 1 Mb of known MM GWAS risk variants. Of these, 94 genes, located in eight regions, had not previously been considered as a candidate gene for that locus. CONCLUSIONS: Our findings highlight the value of leveraging expression data from multiple tissues to identify candidate genes responsible for GWAS associations which provide insight into MM tumorigenesis. Among the genes identified, a number have plausible roles in MM biology, notably APOBEC3C, APOBEC3H, APOBEC3D, APOBEC3F, APOBEC3G, or have been previously implicated in other malignancies. The genes identified in this TWAS can be explored for follow-up and validation to further understand their role in MM biology.

Enrichment of B cell receptor signaling and epidermal growth factor receptor pathways in monoclonal gammopathy of undetermined significance: a genome-wide genetic interaction study.

BACKGROUND: Recent identification of 10 germline variants predisposing to monoclonal gammopathy of undetermined significance (MGUS) explicates genetic dependency of this asymptomatic precursor condition with multiple myeloma (MM). Yet much of genetic burden as well as functional links remain unexplained. We propose a workflow to expand the search for susceptibility loci with genome-wide interaction and for subsequent identification of genetic clusters and pathways. METHODS: Polygenic interaction analysis on 243 cases/1285 controls identified 14 paired risk loci belonging to unique chromosomal bands which were then replicated in two independent sets (case only study, 82 individuals; case/control study 236 cases/ 2484 controls). Further investigation on gene-set enrichment, regulatory pathway and genetic network was carried out with stand-alone in silico tools separately for both interaction and genome-wide association study-detected risk loci. RESULTS: Intronic-PREX1 (20q13.13), a reported locus predisposing to MM was confirmed to have contribution to excess MGUS risk in interaction with SETBP1, a well-established candidate predisposing to myeloid malignancies. Pathway enrichment showed B cell receptor signaling pathway (P < 5.3 × 10- 3) downstream to allograft rejection pathway (P < 5.6 × 10- 4) and autoimmune thyroid disease pathway (P < 9.3 × 10- 4) as well as epidermal growth factor receptor regulation pathway (P < 2.4 × 10- 2) to be differentially regulated. Oncogene ALK and CDH2 were also identified to be moderately interacting with rs10251201 and rs16966921, two previously reported risk loci for MGUS. CONCLUSIONS: We described novel pathways and variants potentially causal for MGUS. The methodology thus proposed to facilitate our search streamlines risk locus-based interaction, genetic network and pathway enrichment analyses.

Chemotherapy-induced peripheral neuropathy: evidence from genome-wide association studies and replication within multiple myeloma patients.

BACKGROUND: Based on the possible shared mechanisms of chemotherapy-induced peripheral neuropathy (CIPN) for different drugs, we aimed to aggregate results of all previously published genome-wide association studies (GWAS) on CIPN, and to replicate them within a cohort of multiple myeloma (MM) patients. METHODS: Following a systematic literature search, data for CIPN associated single nucleotide polymorphisms (SNPs) with P-values< 10- 5 were extracted; these associations were investigated within a cohort of 983 German MM patients treated with bortezomib, thalidomide or vincristine. Cases were subjects that developed CIPN grade 2-4 while controls developed no or sub-clinical CIPN. Logistic regression with additive model was used. RESULTS: In total, 9 GWASs were identified from the literature on CIPN caused by different drugs (4 paclitaxel, 2 bortezomib, 1 vincristine, 1 docetaxel, and 1 oxaliplatin). Data were extracted for 526 SNPs in 109 loci. One hundred fourty-eight patients in our study population were CIPN cases (102/646 bortezomib, 17/63 thalidomide and 29/274 vincristine). In total, 13 SNPs in 9 loci were replicated in our population (p-value< 0.05). The four smallest P-values relevant to the nerve function were 0.0006 for rs8014839 (close to the FBXO33 gene), 0.004 for rs4618330 (close to the INTU gene), 0.006 for rs1903216 (close to the BCL6 gene) and 0.03 for rs4687753 (close to the IL17RB gene). CONCLUSIONS: Replicated SNPs provide clues of the molecular mechanism of CIPN and can be strong candidates for further research aiming to predict the risk of CIPN in clinical practice, particularly rs8014839, rs4618330, rs1903216, and rs4687753, which showed relevance to the function of nervous system.

Bortezomib-induced peripheral neuropathy: A genome-wide association study on multiple myeloma patients.

The proteasome-inhibitor bortezomib was introduced into the treatment of multiple myeloma more than a decade ago. It is clinically beneficial, but peripheral neuropathy (PNP) is a side effect that may limit its use in some patients. To examine the possible genetic predisposing factors to PNP, we performed a genome-wide association study on 646 bortezomib-treated German multiple myeloma patients. Our aim was to identify genetic risk variants associated with the development of PNP as a serious side effect of the treatment. We identified 4 new promising loci for bortezomib-induced PNP at 4q34.3 (rs6552496), 5q14.1 (rs12521798), 16q23.3 (rs8060632), and 18q21.2 (rs17748074). Even though the results did not reach genome-wide significance level, they support the idea of previous studies, suggesting a genetic basis for neurotoxicity. The identified single nucleotide polymorphisms map to genes or next to genes involved in the development and function of the nervous system (CDH13, DCC, and TENM3). As possible functional clues, 2 of the variants, rs12521798 and rs17748074, affect enhancer histone marks in the brain. The rs12521798 may also impact expression of THBS4, which affects specific signal trasduction pathways in the nervous system. Further research is needed to clarify the mechanism of action of the identified single nucleotide polymorphisms in the development of drug-induced PNP and to functionally validate our in silico predictions.

Functional germline variants in driver genes of breast cancer.

PURPOSE: Germline mutations in tumour suppressor genes cause various cancers. These genes are also somatically mutated in sporadic tumours. We hypothesized that there may also be cancer-related germline variants in the genes commonly mutated in sporadic breast tumours. METHODS: After excluding the well-characterized breast cancer (BC) genes, we screened 15 novel genes consistently classified as BC driver genes in next-generation sequencing approaches for single nucleotide polymorphisms (SNPs). Altogether 40 SNPs located in the core promoter, 5'- and 3'-UTR or which were nonsynonymous SNPs were genotyped in 782 Swedish incident BC cases and 1,559 matched controls. After statistical analyses, further evaluations related to functional prediction and signatures of selection were performed. RESULTS: TBX3 was associated with BC risk (rs2242442: OR = 0.76, 95% CI 0.64-0.92, dominant model) and with less aggressive tumour characteristics. An association with BC survival and aggressive tumour characteristics was detected for the genes ATR (rs2227928: HR = 1.63; 95% CI 1.00-2.64, dominant model), RUNX1 (rs17227210: HR = 3.50, 95% CI 1.42-8.61, recessive model) and TTN (rs2303838: HR = 2.36; 95% CI 1.04-5.39; rs2042996: HR = 2.28; 95% CI 1.19-4.37, recessive model). According to the experimental ENCODE data all these SNPs themselves or SNPs in high linkage disequilibrium with them (r 2 ≥ 0.80) were located in regulatory regions. RUNX1 and TTN showed also several signatures of positive selection. CONCLUSION: The study gave evidence that germline variants in BC driver genes may have impact on BC risk and/or survival. Future studies could discover further germline variants in known or so far unknown driver genes which contribute to cancer development.

Genetic Susceptibility to Bortezomib-Induced Peripheral Neuroropathy: Replication of the Reported Candidate Susceptibility Loci.

The introduction of proteasome inhibitors in the treatment of multiple myeloma (MM) patients has been a therapeutic success. Peripheral neuropathy (PNP) remains one of the most frequent side-effects experienced by patients who receive these novel agents. Recent investigations on the mechanisms of PNP in patients treated with bortezomib have suggested genetic susceptibility to neurotoxicity. We used data from a genome-wide association study conducted on 646 bortezomib-treated German MM patients to replicate the previously reported associations between single-nucleotide polymorphisms (SNPs) in candidate genes and PNP in MM patients, including 298 SNPs with a nominal significance (p value <0.05). Twelve associations were confirmed at a significance level p value <0.05. The corresponding SNPs are located in genes involved in drug metabolism (ABCC1, ABCC6), development and function of the nervous system (POGZ, NFAT pathway, EDN1), modulation of immune responses (IL17RD, IL10RA) and the NF-κB signaling pathway (PSMB4, BTCR, F2). We systematically investigated functional consequences of those variants using several bioinformatics tools, such as HaploRegV4.1, RegulomeDB and UCSC Genome Browser. Expression quantitative trait loci (eQTL) data suggested that some of the identified SNPs might influence gene expression through a differential recruitment of transcription factors. In conclusion, we confirmed some of the recently reported associations between germline variation and PNP. Elucidating the mechanisms underlying these associations will contribute to the development of new strategies for the prevention or reduction of PNP.

Genomewide association study on monoclonal gammopathy of unknown significance (MGUS).

OBJECTIVES: To identify germ line variants contributing to the development of monoclonal gammopathy of undetermined significance (MGUS), an asymptomatic premalignant precursor for multiple myeloma (MM). METHODS: We conducted the first genomewide association study (GWAS) on MGUS on 243 German cases with a replication on 294 Czech cases. Identified loci were further analyzed in 1508 German MM patients. New MM loci recently reported in a meta-analysis were also tested in the MGUS GWAS. RESULTS: In GWAS, we identified 10 loci contributing to development of MGUS at P-value threshold of 10-5 . The Czech cohort gave support for two associations (6q26, rs6933936; 7p21.3 rs10251201). In GWAS, rs974120 (8p23.2) reached genomewide significance (P=2.94×10-9 ), with a nominal significance in MM. The locus of rs974120 shows marks of transcriptional activity in leukemia according to ENCODE data. rs10251201 (7p21.3), rs9318227 (13q22.1), and rs10405859 (19q13.32) were associated with markers related to leukemogenesis and immune and inflammatory responses. Two newly identified candidate loci for MM, rs1948915 (8q24.21) and rs8058578 (16p11.2), were nominally associated with MGUS. CONCLUSIONS: These data allow a cautious first proposal for a germ line architecture of MGUS with links to leukemia and autoimmune conditions, the latter agreeing with a family study showing clustering of MGUS with autoimmune diseases.

Variation at 10p12.2 and 10p14 influences risk of childhood B-cell acute lymphoblastic leukemia and phenotype.

Acute lymphoblastic leukemia (ALL) is the major pediatric cancer diagnosed in economically developed countries with B-cell precursor (BCP)-ALL, accounting for approximately 70% of ALL. Recent genome-wide association studies (GWAS) have provided the first unambiguous evidence for common inherited susceptibility to BCP-ALL, identifying susceptibility loci at 7p12.2, 9p21.3, 10q21.2, and 14q11.2. To identify additional BCP-ALL susceptibility loci, we conducted a GWAS and performed a meta-analysis with a published GWAS totaling 1658 cases and 4723 controls, with validation in 1449 cases and 1488 controls. Combined analysis identified novel loci mapping to 10p12.2 (rs10828317, odds ratio [OR] = 1.23; P = 2.30 × 10(-9)) and 10p14 marked by rs3824662 (OR = 1.31; P = 8.62 × 10(-12)). The single nucleotide polymorphism rs10828317 is responsible for the N215S polymorphism in exon 7 of PIP4K2A, and rs3824662 localizes to intron 3 of the transcription factor and putative tumor suppressor gene GATA3. The rs10828317 association was shown to be specifically associated with hyperdiploid ALL, whereas the rs3824662-associated risk was confined to nonhyperdiploid non-TEL-AML1 + ALL. The risk allele of rs3824662 was correlated with older age at diagnosis (P < .001) and significantly worse event-free survivorship (P < .0001). These findings provide further insights into the genetic and biological basis of inherited genetic susceptibility to BCP-ALL and the influence of constitutional genotype on disease development.

Polymorphisms within Autophagy-Related Genes Influence the Risk of Developing Colorectal Cancer: A Meta-Analysis of Four Large Cohorts.

The role of genetic variation in autophagy-related genes in modulating autophagy and cancer is poorly understood. Here, we comprehensively investigated the association of autophagy-related variants with colorectal cancer (CRC) risk and provide new insights about the molecular mechanisms underlying the associations. After meta-analysis of the genome-wide association study (GWAS) data from four independent European cohorts (8006 CRC cases and 7070 controls), two loci, DAPK2 (p = 2.19 × 10-5) and ATG5 (p = 6.28 × 10-4) were associated with the risk of CRC. Mechanistically, the DAPK2rs11631973G allele was associated with IL1 ß levels after the stimulation of peripheral blood mononuclear cells (PBMCs) with Staphylococcus aureus (p = 0.002), CD24 + CD38 + CD27 + IgM + B cell levels in blood (p = 0.0038) and serum levels of en-RAGE (p = 0.0068). ATG5rs546456T allele was associated with TNF α and IL1 ß levels after the stimulation of PBMCs with LPS (p = 0.0088 and p = 0.0076, respectively), CD14+CD16- cell levels in blood (p = 0.0068) and serum levels of CCL19 and cortisol (p = 0.0052 and p = 0.0074, respectively). Interestingly, no association with autophagy flux was observed. These results suggested an effect of the DAPK2 and ATG5 loci in the pathogenesis of CRC, likely through the modulation of host immune responses.

Loci associated with genomic damage levels in chronic kidney disease patients and controls.

Chronic kidney disease (CKD) is a multifactorial disorder with an important genetic component, and several studies have demonstrated potential associations with allelic variants. In addition, CKD patients are also characterized by high levels of genomic damage. Nevertheless, no studies have established relationships between DNA damage, or genomic instability present in CKD patients, and gene polymorphisms. To fill in this gap, the potential role of polymorphisms in genes involved in base excision repair (OGG1, rs1052133; MUTYH, rs3219489; XRCC1, rs25487), nucleotide excision repair (ERCC2/XPD, rs1799793, rs171140, rs13181; ERCC4, rs3136166); phase II metabolism (GSTP1, rs749174; GSTO1, rs2164624; GSTO2, rs156697), and antioxidant enzymes (SOD1, rs17880135, rs1041740, rs202446; SOD2, rs4880; CAT, rs1001179; GPX1, rs17080528; GPX3, rs870406: GPX4, rs713041) were inquired. In addition, some genes involved in CKD (AGT, rs5050; GLO1, rs386572987; SHROOM3, rs17319721) were also evaluated. The genomic damage, the genomic instability, and oxidative damage were evaluated by using the micronucleus and the comet assay in 589 donors (415 CKD patients and 174 controls). Our results showed significant associations between genomic damage and genes directly involved in DNA repair pathways (XRCC1, and ERCC2), and with genes encoding for antioxidant enzymes (SOD1 and GPX1). GSTO2, as a gene involved in phase II metabolism, and MUTYH showed also an association with genomic instability. Interestingly, the three genes associated with CKD (AGT, GLO1, and SHROOM3) showed associations with both the high levels of oxidatively damaged DNA and genomic instability. These results support our view that genomic instability can be considered a biomarker of the CKD status.

Epistatic effect of TLR3 and cGAS-STING-IKKε-TBK1-IFN signaling variants on colorectal cancer risk.

OBJECTIVE: The TLR3/cGAS-STING-IFN signaling has recently been reported to be disturbed in colorectal cancer due to deregulated expression of the genes involved. Our study aimed to investigate the influence of potential regulatory variants in these genes on the risk of sporadic colorectal cancer (CRC) in a Czech cohort of 1424 CRC patients and 1114 healthy controls. METHODS: The variants in the TLR3, CGAS, TMEM173, IKBKE, and TBK1 genes were selected using various online bioinformatic tools, such as UCSC browser, HaploReg, Regulome DB, Gtex Portal, SIFT, PolyPhen2, and miRNA prediction tools. RESULTS: Logistic regression analysis adjusted for age and sex detected a nominal association between CRC risk and three variants, CGAS rs72960018 (OR: 1.68, 95% CI: 1.11-2.53, P-value = .01), CGAS rs9352000 (OR: 2.02, 95% CI: 1.07-3.84, P-value = .03) and TMEM173 rs13153461 (OR: 1.53, 95% CI: 1.03-2.27, P-value = .03). Their cumulative effect revealed a threefold increased CRC risk in carriers of 5-6 risk alleles compared to those with 0-2 risk alleles. Epistatic interactions between these genes and the previously genotyped IFNAR1, IFNAR2, IFNA, IFNB, IFNK, IFNW, IRF3, and IRF7 genes, were computed to test their effect on CRC risk. Overall, we obtained nine pair-wise interactions within and between the CGAS, TMEM173, IKBKE, and TBK1 genes. Two of them remained statistically significant after Bonferroni correction. Additional 52 interactions were observed when IFN variants were added to the analysis. CONCLUSIONS: Our data suggest that epistatic interactions and a high number of risk alleles may play an important role in CRC carcinogenesis, offering novel biological understanding for the CRC management.

Genome-wide interaction and pathway-based identification of key regulators in multiple myeloma.

Inherited genetic susceptibility to multiple myeloma has been investigated in a number of studies. Although 23 individual risk loci have been identified, much of the genetic heritability remains unknown. Here we carried out genome-wide interaction analyses on two European cohorts accounting for 3,999 cases and 7,266 controls and characterized genetic susceptibility to multiple myeloma with subsequent meta-analysis that discovered 16 unique interacting loci. These risk loci along with previously known variants explain 17% of the heritability in liability scale. The genes associated with the interacting loci were found to be enriched in transforming growth factor beta signaling and circadian rhythm regulation pathways suggesting immunoglobulin trait modulation, TH17 cell differentiation and bone morphogenesis as mechanistic links between the predisposition markers and intrinsic multiple myeloma biology. Further tissue/cell-type enrichment analysis associated the discovered genes with hemic-immune system tissue types and immune-related cell types indicating overall involvement in immune response.

Association of NLRP1 Coding Polymorphism with Lung Function and Serum IL-1ß Concentration in Patients Diagnosed with Chronic Obstructive Pulmonary Disease (COPD).

Chronic obstructive pulmonary disease (COPD) is a chronic disease characterized by a progressive decline in lung function due to airflow limitation, mainly related to IL-1ß-induced inflammation. We have hypothesized that single nucleotide polymorphisms (SNPs) in NLRP genes, coding for key regulators of IL-1ß, are associated with pathogenesis and clinical phenotypes of COPD. We recruited 704 COPD individuals and 1238 healthy controls for this study. Twenty non-synonymous SNPs in 10 different NLRP genes were genotyped. Genetic associations were estimated using logistic regression, adjusting for age, gender, and smoking history. The impact of genotypes on patients' overall survival was analyzed with the Kaplan-Meier method with the log-rank test. Serum IL-1ß concentration was determined by high sensitivity assay and expression analysis was done by RT-PCR. Decreased lung function, measured by a forced expiratory volume in 1 s (FEV1% predicted), was significantly associated with the minor allele genotypes (AT + TT) of NLRP1 rs12150220 (p = 0.0002). The same rs12150220 genotypes exhibited a higher level of serum IL-1ß compared to the AA genotype (p = 0.027) in COPD patients. NLRP8 rs306481 minor allele genotypes (AG + AA) were more common in the Global Initiative for Chronic Obstructive Lung Disease (GOLD) definition of group A (p = 0.0083). Polymorphisms in NLRP1 (rs12150220; OR = 0.55, p = 0.03) and NLRP4 (rs12462372; OR = 0.36, p = 0.03) were only nominally associated with COPD risk. In conclusion, coding polymorphisms in NLRP1 rs12150220 show an association with COPD disease severity, indicating that the fine-tuning of the NLRP1 inflammasome could be important in maintaining lung tissue integrity and treating the chronic inflammation of airways.

Author Correction: Identification of multiple risk loci and regulatory mechanisms influencing susceptibility to multiple myeloma.

The original version of this Article contained an error in the spelling of a member of the PRACTICAL Consortium, Manuela Gago-Dominguez, which was incorrectly given as Manuela Gago Dominguez. This has now been corrected in both the PDF and HTML versions of the Article. Furthermore, in the original HTML version of this Article, the order of authors within the author list was incorrect. The PRACTICAL consortium was incorrectly listed after Richard S. Houlston and should have been listed after Nora Pashayan. This error has been corrected in the HTML version of the Article; the PDF version was correct at the time of publication.

Author Correction: Genome-wide association study of classical Hodgkin lymphoma identifies key regulators of disease susceptibility.

The original version of this Article contained an error in the spelling of a member of the PRACTICAL Consortium, Manuela Gago-Dominguez, which was incorrectly given as Manuela Gago Dominguez. This has now been corrected in both the PDF and HTML versions of the Article.

Levels of DNA damage (Micronuclei) in patients suffering from chronic kidney disease. Role of GST polymorphisms.

Chronic kidney disease (CKD) patients are characterized by the presence of high levels of DNA damage, and a poor response to ionizing radiation. In this study, we proposed that variants in GST genes could explain this fact. One-hundred twenty seven CKD patients and one-hundred forty five controls constituted the studied groups. Micronuclei (MN) frequency was determined in peripheral blood lymphocytes at both basal level, and after challenging the cells with 0.5 Gy of ionizing radiation. The following polymorphisms: GSTP1 (rs749174), GSTO1 (rs2164624), and GSTO2 (rs156697) were evaluated in the two groups. Results indicate that gene variants were distributed differentially between CKD patients and controls. Although GSTO1 and GSTO2 variants were associated with lower levels of MN, this was observed in both CKD patients and controls. When net MN values were determined after irradiation, GSTO1 and GSTO2 variants were also associated with lower MN-frequencies. On the contrary, individuals with the GSTP1 variant showed higher values of induced MN. In conclusion, we have demonstrate that the selected GST polymorphism play a role in the incidence of CKD, and affects the levels of MN. Interestingly, the modulating effects observed on both, the basal and induced levels of DNA damage, are characteristic of the overall population, not only of the CKD patients.

Investigation of single and synergic effects of NLRC5 and PD-L1 variants on the risk of colorectal cancer.

Constitutive activation of interferon signaling pathways has been reported in colorectal cancer (CRC), leading to a strong CD8+ T cell response through stimulation of NLRC5 expression. Primed CD8+ T cell expansion, however, may be negatively regulated by PD-L1 expression. Additionally, aberrant PD-L1 expression enables cancer cells to escape the immune attack. Our study aimed to select potential regulatory variants in the NLRC5 and PD-L1 genes by using several online in silico tools, such as UCSC browser, HaploReg, Regulome DB, Gtex Portal, microRNA and transcription factor binding site prediction tools and to investigate their influence on CRC risk in a Czech cohort of 1424 CRC patients and 1114 healthy controls. Logistic regression analysis adjusted for age and gender reported a moderate association between rectal cancer risk and two NLRC5 SNPs, rs1684575 T>G (OR: 1.60, 95% CI: 1.13-2.27, recessive model) and rs3751710 (OR: 0.70, 95% CI: 0.51-0.96, dominant model). Given that a combination of genetic variants, rather than a single polymorphism, may explain better the genetic etiology of CRC, we studied the interplay between the variants within NLRC5, PD-L1 and the previously genotyped IFNGR1 and IFNGR2 variants, to evaluate their involvement in the risk of CRC development. Overall we obtained 18 pair-wise interactions within and between the NLRC5 ad PD-L1 genes and 6 more when IFNGR variants were added. Thirteen out of the 24 interactions were below the threshold for the FDR calculated and controlled at an arbitrary level q*<0.10. Furthermore, the interaction IFNGR2 rs1059293 C>T-NLRC5 rs289747 G>A (P<0.0001) remained statistically significant even after Bonferroni correction. Our data suggest that not only a single genetic variant but also an interaction between two or more variants within genes involved in immune regulation may play important roles in the onset of CRC, providing therefore novel biological information, which could eventually improve CRC risk management but also PD-1-based immunotherapy in CRC.

Identification of multiple risk loci and regulatory mechanisms influencing susceptibility to multiple myeloma.

Genome-wide association studies (GWAS) have transformed our understanding of susceptibility to multiple myeloma (MM), but much of the heritability remains unexplained. We report a new GWAS, a meta-analysis with previous GWAS and a replication series, totalling 9974 MM cases and 247,556 controls of European ancestry. Collectively, these data provide evidence for six new MM risk loci, bringing the total number to 23. Integration of information from gene expression, epigenetic profiling and in situ Hi-C data for the 23 risk loci implicate disruption of developmental transcriptional regulators as a basis of MM susceptibility, compatible with altered B-cell differentiation as a key mechanism. Dysregulation of autophagy/apoptosis and cell cycle signalling feature as recurrently perturbed pathways. Our findings provide further insight into the biological basis of MM.

Genetic correlation between multiple myeloma and chronic lymphocytic leukaemia provides evidence for shared aetiology.

The clustering of different types of B-cell malignancies in families raises the possibility of shared aetiology. To examine this, we performed cross-trait linkage disequilibrium (LD)-score regression of multiple myeloma (MM) and chronic lymphocytic leukaemia (CLL) genome-wide association study (GWAS) data sets, totalling 11,734 cases and 29,468 controls. A significant genetic correlation between these two B-cell malignancies was shown (Rg = 0.4, P = 0.0046). Furthermore, four of the 45 known CLL risk loci were shown to associate with MM risk and five of the 23 known MM risk loci associate with CLL risk. By integrating eQTL, Hi-C and ChIP-seq data, we show that these pleiotropic risk loci are enriched for B-cell regulatory elements and implicate B-cell developmental genes. These data identify shared biological pathways influencing the development of CLL and, MM and further our understanding of the aetiological basis of these B-cell malignancies.