Enhanced antitumour immunity following neoadjuvant chemoradiotherapy mediates a favourable prognosis in women with resected pancreatic cancer.

BACKGROUND: This study investigates sex disparities in clinical outcomes and tumour immune profiles in patients with pancreatic ductal adenocarcinoma (PDAC) who underwent upfront resection or resection preceded by gemcitabine-based neoadjuvant chemoradiotherapy (nCRT). METHODS: Patients originated from the PREOPANC randomised controlled trial. Upfront surgery was performed in 82 patients, and 66 received nCRT before resection. The impact of sex on overall survival (OS) was investigated using Cox proportional hazards models. The immunological landscape within the tumour microenvironment (TME) was mapped using transcriptomic and spatial proteomic profiling. RESULTS: The 5-year OS rate differed between the sexes following resection preceded by nCRT, with 43% for women compared with 22% for men. In multivariate analysis, the female sex was a favourable independent prognostic factor for OS only in the nCRT group (HR 0.19; 95% CI 0.07 to 0.52). Multivariate heterogeneous treatment effects analysis revealed a significant interaction between sex and treatment, implying increased nCRT efficacy among women with resected PDAC. The TME of women contained fewer protumoural CD163+MRC1+M2 macrophages than that of men after nCRT, as indicated by transcriptomic and validated using spatial proteomic profiling. CONCLUSION: PDAC tumours of women are more sensitive to gemcitabine-based nCRT, resulting in longer OS after resection compared with men. This may be due to enhanced immunity impeding the infiltration of protumoral M2 macrophages into the TME. Our findings highlight the importance of considering sex disparities and mitigating immunosuppressive macrophage polarisation for personalised PDAC treatment.

New proposal to address mediation analysis interrogations by using genetic variants as instrumental variables.

The application of causal mediation analysis (CMA) considering the mediation effect of a third variable is increasing in epidemiological studies; however, this requires fitting strong assumptions on confounding bias. To address this limitation, we propose an extension of CMA combining it with Mendelian randomization (MRinCMA). We applied the new approach to analyse the causal effect of obesity and diabetes on pancreatic cancer, considering each factor as potential mediator. To check the performance of MRinCMA under several conditions/scenarios, we used it in different simulated data sets and compared it with structural equation models. For continuous variables, MRinCMA and structural equation models performed similarly, suggesting that both approaches are valid to obtain unbiased estimates. When noncontinuous variables were considered, MRinCMA presented, overall, lower bias than structural equation models. By applying MRinCMA, we did not find any evidence of causality of obesity or diabetes on pancreatic cancer. With this new methodology, researchers would be able to address CMA hypotheses by appropriately accounting for the confounding bias assumption regardless of the conditions used in their studies in different settings.

Transcriptional regulation by NR5A2 links differentiation and inflammation in the pancreas.

Chronic inflammation increases the risk of developing one of several types of cancer. Inflammatory responses are currently thought to be controlled by mechanisms that rely on transcriptional networks that are distinct from those involved in cell differentiation. The orphan nuclear receptor NR5A2 participates in a wide variety of processes, including cholesterol and glucose metabolism in the liver, resolution of endoplasmic reticulum stress, intestinal glucocorticoid production, pancreatic development and acinar differentiation. In genome-wide association studies, single nucleotide polymorphisms in the vicinity of NR5A2 have previously been associated with the risk of pancreatic adenocarcinoma. In mice, Nr5a2 heterozygosity sensitizes the pancreas to damage, impairs regeneration and cooperates with mutant Kras in tumour progression. Here, using a global transcriptomic analysis, we describe an epithelial-cell-autonomous basal pre-inflammatory state in the pancreas of Nr5a2+/- mice that is reminiscent of the early stages of pancreatitis-induced inflammation and is conserved in histologically normal human pancreases with reduced expression of NR5A2 mRNA. In Nr5a2+/-mice, NR5A2 undergoes a marked transcriptional switch, relocating from differentiation-specific to inflammatory genes and thereby promoting gene transcription that is dependent on the AP-1 transcription factor. Pancreatic deletion of Jun rescues the pre-inflammatory phenotype, as well as binding of NR5A2 to inflammatory gene promoters and the defective regenerative response to damage. These findings support the notion that, in the pancreas, the transcriptional networks involved in differentiation-specific functions also suppress inflammatory programmes. Under conditions of genetic or environmental constraint, these networks can be subverted to foster inflammation.

GATA4 and GATA6 loss-of-expression is associated with extinction of the classical programme and poor outcome in pancreatic ductal adenocarcinoma.

OBJECTIVE: GATA6 is a key regulator of the classical phenotype in pancreatic ductal adenocarcinoma (PDAC). Low GATA6 expression associates with poor patient outcome. GATA4 is the second most expressed GATA factor in the pancreas. We assessed whether, and how, GATA4 contributes to PDAC phenotype and analysed the association of expression with outcome and response to chemotherapy. DESIGN: We analysed PDAC transcriptomic data, stratifying cases according to GATA4 and GATA6 expression and identified differentially expressed genes and pathways. The genome-wide distribution of GATA4 was assessed, as well as the effects of GATA4 knockdown. A multicentre tissue microarray study to assess GATA4 and GATA6 expression in samples (n=745) from patients with resectable was performed. GATA4 and GATA6 levels were dichotomised into high/low categorical variables; association with outcome was assessed using univariable and multivariable Cox regression models. RESULTS: GATA4 messenger RNA is enriched in classical, compared with basal-like tumours. We classified samples in 4 groups as high/low for GATA4 and GATA6. Reduced expression of GATA4 had a minor transcriptional impact but low expression of GATA4 enhanced the effects of GATA6 low expression. GATA4 and GATA6 display a partially overlapping genome-wide distribution, mainly at promoters. Reduced expression of both proteins in tumours was associated with the worst patient survival. GATA4 and GATA6 expression significantly decreased in metastases and negatively correlated with basal markers. CONCLUSIONS: GATA4 and GATA6 cooperate to maintain the classical phenotype. Our findings provide compelling rationale to assess their expression as biomarkers of poor prognosis and therapeutic response.

Metagenomic Identification of Microbial Signatures Predicting Pancreatic Cancer From a Multinational Study.

BACKGROUND & AIMS: To identify gut and oral metagenomic signatures that accurately predict pancreatic ductal carcinoma (PDAC) and to validate these signatures in independent cohorts. METHODS: We conducted a multinational study and performed shotgun metagenomic analysis of fecal and salivary samples collected from patients with treatment-naïve PDAC and non-PDAC controls in Japan, Spain, and Germany. Taxonomic and functional profiles of the microbiomes were characterized, and metagenomic classifiers to predict PDAC were constructed and validated in external datasets. RESULTS: Comparative metagenomics revealed dysbiosis of both the gut and oral microbiomes and identified 30 gut and 18 oral species significantly associated with PDAC in the Japanese cohort. These microbial signatures achieved high area under the curve values of 0.78 to 0.82. The prediction model trained on the Japanese gut microbiome also had high predictive ability in Spanish and German cohorts, with respective area under the curve values of 0.74 and 0.83, validating its high confidence and versatility for PDAC prediction. Significant enrichments of Streptococcus and Veillonella spp and a depletion of Faecalibacterium prausnitzii were common gut signatures for PDAC in all the 3 cohorts. Prospective follow-up data revealed that patients with certain gut and oral microbial species were at higher risk of PDAC-related mortality. Finally, 58 bacteriophages that could infect microbial species consistently enriched in patients with PDAC across the 3 countries were identified. CONCLUSIONS: Metagenomics targeting the gut and oral microbiomes can provide a powerful source of biomarkers for identifying individuals with PDAC and their prognoses. The identification of shared gut microbial signatures for PDAC in Asian and European cohorts indicates the presence of robust and global gut microbial biomarkers.

A faecal microbiota signature with high specificity for pancreatic cancer.

BACKGROUND: Recent evidence suggests a role for the microbiome in pancreatic ductal adenocarcinoma (PDAC) aetiology and progression. OBJECTIVE: To explore the faecal and salivary microbiota as potential diagnostic biomarkers. METHODS: We applied shotgun metagenomic and 16S rRNA amplicon sequencing to samples from a Spanish case-control study (n=136), including 57 cases, 50 controls, and 29 patients with chronic pancreatitis in the discovery phase, and from a German case-control study (n=76), in the validation phase. RESULTS: Faecal metagenomic classifiers performed much better than saliva-based classifiers and identified patients with PDAC with an accuracy of up to 0.84 area under the receiver operating characteristic curve (AUROC) based on a set of 27 microbial species, with consistent accuracy across early and late disease stages. Performance further improved to up to 0.94 AUROC when we combined our microbiome-based predictions with serum levels of carbohydrate antigen (CA) 19-9, the only current non-invasive, Food and Drug Administration approved, low specificity PDAC diagnostic biomarker. Furthermore, a microbiota-based classification model confined to PDAC-enriched species was highly disease-specific when validated against 25 publicly available metagenomic study populations for various health conditions (n=5792). Both microbiome-based models had a high prediction accuracy on a German validation population (n=76). Several faecal PDAC marker species were detectable in pancreatic tumour and non-tumour tissue using 16S rRNA sequencing and fluorescence in situ hybridisation. CONCLUSION: Taken together, our results indicate that non-invasive, robust and specific faecal microbiota-based screening for the early detection of PDAC is feasible.

Avoiding collider bias in Mendelian randomization when performing stratified analyses.

Mendelian randomization (MR) uses genetic variants as instrumental variables to investigate the causal effect of a risk factor on an outcome. A collider is a variable influenced by two or more other variables. Naive calculation of MR estimates in strata of the population defined by a collider, such as a variable affected by the risk factor, can result in collider bias. We propose an approach that allows MR estimation in strata of the population while avoiding collider bias. This approach constructs a new variable, the residual collider, as the residual from regression of the collider on the genetic instrument, and then calculates causal estimates in strata defined by quantiles of the residual collider. Estimates stratified on the residual collider will typically have an equivalent interpretation to estimates stratified on the collider, but they are not subject to collider bias. We apply the approach in several simulation scenarios considering different characteristics of the collider variable and strengths of the instrument. We then apply the proposed approach to investigate the causal effect of smoking on bladder cancer in strata of the population defined by bodyweight. The new approach generated unbiased estimates in all the simulation settings. In the applied example, we observed a trend in the stratum-specific MR estimates at different bodyweight levels that suggested stronger effects of smoking on bladder cancer among individuals with lower bodyweight. The proposed approach can be used to perform MR studying heterogeneity among subgroups of the population while avoiding collider bias.

Deciphering the complex interplay between pancreatic cancer, diabetes mellitus subtypes and obesity/BMI through causal inference and mediation analyses.

OBJECTIVES: To characterise the association between type 2 diabetes mellitus (T2DM) subtypes (new-onset T2DM (NODM) or long-standing T2DM (LSDM)) and pancreatic cancer (PC) risk, to explore the direction of causation through Mendelian randomisation (MR) analysis and to assess the mediation role of body mass index (BMI). DESIGN: Information about T2DM and related factors was collected from 2018 PC cases and 1540 controls from the PanGenEU (European Study into Digestive Illnesses and Genetics) study. A subset of PC cases and controls had glycated haemoglobin, C-peptide and genotype data. Multivariate logistic regression models were applied to derive ORs and 95% CIs. T2DM and PC-related single nucleotide polymorphism (SNP) were used as instrumental variables (IVs) in bidirectional MR analysis to test for two-way causal associations between PC, NODM and LSDM. Indirect and direct effects of the BMI-T2DM-PC association were further explored using mediation analysis. RESULTS: T2DM was associated with an increased PC risk when compared with non-T2DM (OR=2.50; 95% CI: 2.05 to 3.05), the risk being greater for NODM (OR=6.39; 95% CI: 4.18 to 9.78) and insulin users (OR=3.69; 95% CI: 2.80 to 4.86). The causal association between T2DM (57-SNP IV) and PC was not statistically significant (ORLSDM=1.08, 95% CI: 0.86 to 1.29, ORNODM=1.06, 95% CI: 0.95 to 1.17). In contrast, there was a causal association between PC (40-SNP IV) and NODM (OR=2.85; 95% CI: 2.04 to 3.98), although genetic pleiotropy was present (MR-Egger: p value=0.03). Potential mediating effects of BMI (125-SNPs as IV), particularly in terms of weight loss, were evidenced on the NODM-PC association (indirect effect for BMI in previous years=0.55). CONCLUSION: Findings of this study do not support a causal effect of LSDM on PC, but suggest that PC causes NODM. The interplay between obesity, PC and T2DM is complex.

Plasma protein biomarkers for early detection of pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) has a poor prognosis, mainly due to late diagnosis at advanced tumor stages. In this study, we aimed to identify plasma protein biomarkers for early detection of PDAC. Totally, 135 PDAC patients (early PDAC, Stage I/II, n = 71; advanced PDAC, Stage III/IV, n = 64), 13 benign lesions/chronic pancreatitis patients and 72 healthy individuals, with corresponding plasma samples from a case-control study in Sweden were included. A proximity extension assay was used to detect 92 cancer-related proteins, and an enzyme-linked immunosorbent assay/electrochemiluminescence immunoassay was used to detect CA19-9. Predictive features were selected from these 93 candidate proteins and three covariates in the Swedish participants, and then validated in Spanish participants, including 37 early PDAC patients, 38 advanced PDAC patients, 19 chronic pancreatitis patients and 36 healthy controls. A panel of eight proteins discriminating early PDAC from healthy individuals was identified, and the cross-validated area under the curves (AUCs) were 0.85 (95% confidence interval, 95% CI, 0.78-0.91) and 0.81 (95% CI, 0.70-0.92) in the Swedish and Spanish participants, respectively. Another eight-protein panel was predictive for classifying advanced PDAC from healthy controls in two populations, with cross-validated AUCs of 0.89 (95% CI, 0.83-0.95) and 0.90 (95% CI, 0.83-0.98), respectively. In conclusion, eight protein biomarkers were identified and externally validated, potentially allowing early detection of PDAC patients if validated in additional prospective studies.

Association of patients' sex with treatment outcomes after intravesical bacillus Calmette-Guérin immunotherapy for T1G3/HG bladder cancer.

PURPOSE: To investigate the association of patients' sex with recurrence and disease progression in patients treated with intravesical bacillus Calmette-Guérin (BCG) for T1G3/HG urinary bladder cancer (UBC). MATERIALS AND METHODS: We analyzed the data of 2635 patients treated with adjuvant intravesical BCG for T1 UBC between 1984 and 2019. We accounted for missing data using multiple imputations and adjusted for covariate imbalance between males and females using inverse probability weighting (IPW). Crude and IPW-adjusted Cox regression analyses were used to estimate the hazard ratios (HR) with their 95% confidence intervals (CI) for the association of patients' sex with HG-recurrence and disease progression. RESULTS: A total of 2170 (82%) males and 465 (18%) females were available for analysis. Overall, 1090 (50%) males and 244 (52%) females experienced recurrence, and 391 (18%) males and 104 (22%) females experienced disease progression. On IPW-adjusted Cox regression analyses, female sex was associated with disease progression (HR 1.25, 95%CI 1.01-1.56, p = 0.04) but not with recurrence (HR 1.06, 95%CI 0.92-1.22, p = 0.41). A total of 1056 patients were treated with adequate BCG. In these patients, on IPW-adjusted Cox regression analyses, patients' sex was not associated with recurrence (HR 0.99, 95%CI 0.80-1.24, p = 0.96), HG-recurrence (HR 1.00, 95%CI 0.78-1.29, p = 0.99) or disease progression (HR 1.12, 95%CI 0.78-1.60, p = 0.55). CONCLUSION: Our analysis generates the hypothesis of a differential response to BCG between males and females if not adequately treated. Further studies should focus on sex-based differences in innate and adaptive immune system and their association with BCG response.

A combination of urinary biomarker panel and PancRISK score for earlier detection of pancreatic cancer: A case-control study.

BACKGROUND: Pancreatic ductal adenocarcinoma (PDAC) is one of the deadliest cancers, with around 9% of patients surviving >5 years. Asymptomatic in its initial stages, PDAC is mostly diagnosed late, when already a locally advanced or metastatic disease, as there are no useful biomarkers for detection in its early stages, when surgery can be curative. We have previously described a promising biomarker panel (LYVE1, REG1A, and TFF1) for earlier detection of PDAC in urine. Here, we aimed to establish the accuracy of an improved panel, including REG1B instead of REG1A, and an algorithm for data interpretation, the PancRISK score, in additional retrospectively collected urine specimens. We also assessed the complementarity of this panel with CA19-9 and explored the daily variation and stability of the biomarkers and their performance in common urinary tract cancers. METHODS AND FINDINGS: Clinical specimens were obtained from multiple centres: Barts Pancreas Tissue Bank, University College London, University of Liverpool, Spanish National Cancer Research Center, Cambridge University Hospital, and University of Belgrade. The biomarker panel was assayed on 590 urine specimens: 183 control samples, 208 benign hepatobiliary disease samples (of which 119 were chronic pancreatitis), and 199 PDAC samples (102 stage I-II and 97 stage III-IV); 50.7% were from female individuals. PDAC samples were collected from patients before treatment. The samples were assayed using commercially available ELISAs. Statistical analyses were performed using non-parametric Kruskal-Wallis tests adjusted for multiple comparisons, and multiple logistic regression. Training and validation datasets for controls and PDAC samples were obtained after random division of the whole available dataset in a 1:1 ratio. The substitution of REG1A with REG1B enhanced the performance of the panel to detect resectable PDAC. In a comparison of controls and PDAC stage I-II samples, the areas under the receiver operating characteristic curve (AUCs) increased from 0.900 (95% CI 0.843-0.957) and 0.926 (95% CI 0.843-1.000) in the training (50% of the dataset) and validation sets, respectively, to 0.936 in both the training (95% CI 0.903-0.969) and the validation (95% CI 0.888-0.984) datasets for the new panel including REG1B. This improved panel showed both sensitivity (SN) and specificity (SP) to be >85%. Plasma CA19-9 enhanced the performance of this panel in discriminating PDAC I-II patients from controls, with AUC = 0.992 (95% CI 0.983-1.000), SN = 0.963 (95% CI 0.913-1.000), and SP = 0.967 (95% CI 0.924-1.000). We demonstrate that the biomarkers do not show significant daily variation, and that they are stable for up to 5 days at room temperature. The main limitation of our study is the low number of stage I-IIA PDAC samples (n = 27) and lack of samples from individuals with hereditary predisposition to PDAC, for which specimens collected from control individuals were used as a proxy. CONCLUSIONS: We have successfully validated our urinary biomarker panel, which was improved by substituting REG1A with REG1B. At a pre-selected cutoff of >80% SN and SP for the affiliated PancRISK score, we demonstrate a clinically applicable risk stratification tool with a binary output for risk of developing PDAC ('elevated' or 'normal'). PancRISK provides a step towards precision surveillance for PDAC patients, which we will test in a prospective clinical study, UroPanc.

Resection of pancreatic cancer in Europe and USA: an international large-scale study highlighting large variations.

OBJECTIVE: Resection can potentially cure resectable pancreatic cancer (PaC) and significantly prolong survival in some patients. This large-scale international study aimed to investigate variations in resection for PaC in Europe and USA and determinants for its utilisation. DESIGN: Data from six European population-based cancer registries and the US Surveillance, Epidemiology, and End Results Program database during 2003-2016 were analysed. Age-standardised resection rates for overall and stage I-II PaCs were computed. Associations between resection and demographic and clinical parameters were assessed using multivariable logistic regression models. RESULTS: A total of 153 698 records were analysed. In population-based registries in 2012-2014, resection rates ranged from 13.2% (Estonia) to 21.2% (Slovenia) overall and from 34.8% (Norway) to 68.7% (Denmark) for stage I-II tumours, with great international variations. During 2003-2014, resection rates only increased in USA, the Netherlands and Denmark. Resection was significantly less frequently performed with more advanced tumour stage (ORs for stage III and IV versus stage I-II tumours: 0.05-0.18 and 0.01-0.06 across countries) and increasing age (ORs for patients 70-79 and ≥80 versus those <60 years: 0.37-0.63 and 0.03-0.16 across countries). Patients with advanced-stage tumours (stage III-IV: 63.8%-81.2%) and at older ages (≥70 years: 52.6%-59.5%) receiving less frequently resection comprised the majority of diagnosed cases. Patient performance status, tumour location and size were also associated with resection application. CONCLUSION: Rates of PaC resection remain low in Europe and USA with great international variations. Further studies are warranted to explore reasons for these variations.

Ambient air pollution and incident bladder cancer risk: Updated analysis of the Spanish Bladder Cancer Study.

Although outdoor air pollution and particulate matter in outdoor air have been consistently linked with increased lung cancer risk, the evidence for associations at other cancer sites is limited. Bladder cancer shares several risk factors with lung cancer and some positive associations of ambient air pollution and bladder cancer risk have been observed. This study examined associations of ambient air pollution and bladder cancer risk in the large-scale Spanish Bladder Cancer Study. Estimates of ambient fine particulate matter (PM2.5 ) and nitrogen dioxide (NO2 ) concentrations were assigned to the geocoded participant residence of 938 incident bladder cancer cases and 973 hospital controls based on European multicity land-use regression models. Adjusted odds ratios (ORs) and 95% confidence intervals (CI) for associations of ambient air pollution and bladder cancer risk were estimated using unconditional logistic regression models. Overall, there was no clear association between either ambient PM2.5 (OR per 5.9 µg/m3 = 1.06, 95% CI 0.71-1.60) or NO2 (OR per 14.2 µg/m3 = 0.97, 95% CI 0.84-1.13) concentrations and incident bladder cancer risk. There was no clear evidence for effect modification according to age group, sex, region, education, cigarette smoking status, or pack-years. Results were also similar among more residentially stable participants and in two-pollutant models. Overall, there was no clear evidence for associations of ambient PM2.5 and NO2 concentrations and incident bladder cancer risk. Further research in other large-scale population studies is needed with detailed information on measured or modeled estimates of ambient air pollution concentrations and individual level risk factors.

Pancreatic cancer and autoimmune diseases: An association sustained by computational and epidemiological case-control approaches.

Deciphering the underlying genetic basis behind pancreatic cancer (PC) and its associated multimorbidities will enhance our knowledge toward PC control. The study investigated the common genetic background of PC and different morbidities through a computational approach and further evaluated the less explored association between PC and autoimmune diseases (AIDs) through an epidemiological analysis. Gene-disease associations (GDAs) of 26 morbidities of interest and PC were obtained using the DisGeNET public discovery platform. The association between AIDs and PC pointed by the computational analysis was confirmed through multivariable logistic regression models in the PanGen European case-control study population of 1,705 PC cases and 1,084 controls. Fifteen morbidities shared at least one gene with PC in the DisGeNET database. Based on common genes, several AIDs were genetically associated with PC pointing to a potential link between them. An epidemiologic analysis confirmed that having any of the nine AIDs studied was significantly associated with a reduced risk of PC (Odds Ratio (OR) = 0.74, 95% confidence interval (CI) 0.58-0.93) which decreased in subjects having ≥2 AIDs (OR = 0.39, 95%CI 0.21-0.73). In independent analyses, polymyalgia rheumatica, and rheumatoid arthritis were significantly associated with low PC risk (OR = 0.40, 95%CI 0.19-0.89, and OR = 0.73, 95%CI 0.53-1.00, respectively). Several inflammatory-related morbidities shared a common genetic component with PC based on public databases. These molecular links could shed light into the molecular mechanisms underlying PC development and simultaneously generate novel hypotheses. In our study, we report sound findings pointing to an association between AIDs and a reduced risk of PC.

Winner's Curse Correction and Variable Thresholding Improve Performance of Polygenic Risk Modeling Based on Genome-Wide Association Study Summary-Level Data.

Recent heritability analyses have indicated that genome-wide association studies (GWAS) have the potential to improve genetic risk prediction for complex diseases based on polygenic risk score (PRS), a simple modelling technique that can be implemented using summary-level data from the discovery samples. We herein propose modifications to improve the performance of PRS. We introduce threshold-dependent winner's-curse adjustments for marginal association coefficients that are used to weight the single-nucleotide polymorphisms (SNPs) in PRS. Further, as a way to incorporate external functional/annotation knowledge that could identify subsets of SNPs highly enriched for associations, we propose variable thresholds for SNPs selection. We applied our methods to GWAS summary-level data of 14 complex diseases. Across all diseases, a simple winner's curse correction uniformly led to enhancement of performance of the models, whereas incorporation of functional SNPs was beneficial only for selected diseases. Compared to the standard PRS algorithm, the proposed methods in combination led to notable gain in efficiency (25-50% increase in the prediction R2) for 5 of 14 diseases. As an example, for GWAS of type 2 diabetes, winner's curse correction improved prediction R2 from 2.29% based on the standard PRS to 3.10% (P = 0.0017) and incorporating functional annotation data further improved R2 to 3.53% (P = 2×10-5). Our simulation studies illustrate why differential treatment of certain categories of functional SNPs, even when shown to be highly enriched for GWAS-heritability, does not lead to proportionate improvement in genetic risk-prediction because of non-uniform linkage disequilibrium structure.

Genome-wide association study identifies inversion in the CTRB1-CTRB2 locus to modify risk for alcoholic and non-alcoholic chronic pancreatitis.

OBJECTIVE: Alcohol-related pancreatitis is associated with a disproportionately large number of hospitalisations among GI disorders. Despite its clinical importance, genetic susceptibility to alcoholic chronic pancreatitis (CP) is poorly characterised. To identify risk genes for alcoholic CP and to evaluate their relevance in non-alcoholic CP, we performed a genome-wide association study and functional characterisation of a new pancreatitis locus. DESIGN: 1959 European alcoholic CP patients and population-based controls from the KORA, LIFE and INCIPE studies (n=4708) as well as chronic alcoholics from the GESGA consortium (n=1332) were screened with Illumina technology. For replication, three European cohorts comprising 1650 patients with non-alcoholic CP and 6695 controls originating from the same countries were used. RESULTS: We replicated previously reported risk loci CLDN2-MORC4, CTRC, PRSS1-PRSS2 and SPINK1 in alcoholic CP patients. We identified CTRB1-CTRB2 (chymotrypsin B1 and B2) as a new risk locus with lead single-nucleotide polymorphism (SNP) rs8055167 (OR 1.35, 95% CI 1.23 to 1.6). We found that a 16.6 kb inversion in the CTRB1-CTRB2 locus was in linkage disequilibrium with the CP-associated SNPs and was best tagged by rs8048956. The association was replicated in three independent European non-alcoholic CP cohorts of 1650 patients and 6695 controls (OR 1.62, 95% CI 1.42 to 1.86). The inversion changes the expression ratio of the CTRB1 and CTRB2 isoforms and thereby affects protective trypsinogen degradation and ultimately pancreatitis risk. CONCLUSION: An inversion in the CTRB1-CTRB2 locus modifies risk for alcoholic and non-alcoholic CP indicating that common pathomechanisms are involved in these inflammatory disorders.

Integrative eQTL analysis of tumor and host omics data in individuals with bladder cancer.

Integrative analyses of several omics data are emerging. The data are usually generated from the same source material (i.e., tumor sample) representing one level of regulation. However, integrating different regulatory levels (i.e., blood) with those from tumor may also reveal important knowledge about the human genetic architecture. To model this multilevel structure, an integrative-expression quantitative trait loci (eQTL) analysis applying two-stage regression (2SR) was proposed. This approach first regressed tumor gene expression levels with tumor markers and the adjusted residuals from the previous model were then regressed with the germline genotypes measured in blood. Previously, we demonstrated that penalized regression methods in combination with a permutation-based MaxT method (Global-LASSO) is a promising tool to fix some of the challenges that high-throughput omics data analysis imposes. Here, we assessed whether Global-LASSO can also be applied when tumor and blood omics data are integrated. We further compared our strategy with two 2SR-approaches, one using multiple linear regression (2SR-MLR) and other using LASSO (2SR-LASSO). We applied the three models to integrate genomic, epigenomic, and transcriptomic data from tumor tissue with blood germline genotypes from 181 individuals with bladder cancer included in the TCGA Consortium. Global-LASSO provided a larger list of eQTLs than the 2SR methods, identified a previously reported eQTLs in prostate stem cell antigen (PSCA), and provided further clues on the complexity of APBEC3B loci, with a minimal false-positive rate not achieved by 2SR-MLR. It also represents an important contribution for omics integrative analysis because it is easy to apply and adaptable to any type of data.

Use of big data for drug development and for public and personal health and care.

The use of data analytics across the entire healthcare value chain, from drug discovery and development through epidemiology to informed clinical decision for patients or policy making for public health, has seen an explosion in the recent years. The increase in quantity and variety of data available together with the improvement of storing capabilities and analytical tools offer numerous possibilities to all stakeholders (manufacturers, regulators, payers, healthcare providers, decision makers, researchers) but most importantly, it has the potential to improve general health outcomes if we learn how to exploit it in the right way. This article looks at the different sources of data and the importance of unstructured data. It goes on to summarize current and potential future uses in drug discovery, development, and monitoring as well as in public and personal healthcare; including examples of good practice and recent developments. Finally, we discuss the main practical and ethical challenges to unravel the full potential of big data in healthcare and conclude that all stakeholders need to work together towards the common goal of making sense of the available data for the common good.

Asthma status is associated with decreased risk of aggressive urothelial bladder cancer.

Previous studies suggested an association between atopic conditions and specific cancers. The results on the association with urothelial bladder cancer (UBC) are scarce and inconsistent. To evaluate the association between asthma and risk of UBC, we considered 936 cases and 1,022 controls from the Spanish Bladder Cancer/EPICURO Study (86% males, mean age 65.4 years), a multicenter and hospital-based case-control study conducted during 1998-2001. Participants were asked whether they had asthma and detailed information about occupational exposures, smoking habits, dietary factors, medical conditions and history of medication was collected through face-to-face questionnaires performed by trained interviewers. Since asthma and UBC might share risk factors, association between patients' characteristics and asthma was studied in UBC controls. Association between UBC and asthma was assessed using logistic regression unadjusted and adjusted for potential confounders. The complex interrelationships, direct and mediating effect of asthma on UBC, were appraised using counterfactual mediation models. Asthma was associated with a reduced risk of UBC (odds ratio (OR) = 0.54, 95% confidence interval (CI) 0.37, 0.79) after adjusting for a wide range of confounders. No mediating effect was identified. The reduced risk associated with asthma was restricted to patients with high-risk non-muscle invasive (OR = 0.25, 95%CI 0.10, 0.62) and muscle invasive UBC (OR = 0.32, 95%CI 0.15, 0.69). Our results support that asthma is associated with a decreased risk of UBC, especially among aggressive tumors. Further work on the relationship between asthma and other atopic conditions and cancer risk should shed light on the relationship between immune response mechanisms and bladder carcinogenesis.

Identification of a novel susceptibility locus at 13q34 and refinement of the 20p12.2 region as a multi-signal locus associated with bladder cancer risk in individuals of European ancestry.

Candidate gene and genome-wide association studies (GWAS) have identified 15 independent genomic regions associated with bladder cancer risk. In search for additional susceptibility variants, we followed up on four promising single-nucleotide polymorphisms (SNPs) that had not achieved genome-wide significance in 6911 cases and 11 814 controls (rs6104690, rs4510656, rs5003154 and rs4907479, P < 1 × 10(-6)), using additional data from existing GWAS datasets and targeted genotyping for studies that did not have GWAS data. In a combined analysis, which included data on up to 15 058 cases and 286 270 controls, two SNPs achieved genome-wide statistical significance: rs6104690 in a gene desert at 20p12.2 (P = 2.19 × 10(-11)) and rs4907479 within the MCF2L gene at 13q34 (P = 3.3 × 10(-10)). Imputation and fine-mapping analyses were performed in these two regions for a subset of 5551 bladder cancer cases and 10 242 controls. Analyses at the 13q34 region suggest a single signal marked by rs4907479. In contrast, we detected two signals in the 20p12.2 region-the first signal is marked by rs6104690, and the second signal is marked by two moderately correlated SNPs (r(2) = 0.53), rs6108803 and the previously reported rs62185668. The second 20p12.2 signal is more strongly associated with the risk of muscle-invasive (T2-T4 stage) compared with non-muscle-invasive (Ta, T1 stage) bladder cancer (case-case P ≤ 0.02 for both rs62185668 and rs6108803). Functional analyses are needed to explore the biological mechanisms underlying these novel genetic associations with risk for bladder cancer.