Exploring the Associations of Inflammatory and Oxidative Stress Biomarkers with Pancreatic Diseases: An Observational and Mendelian Randomisation Study.

Identifying biomarkers linked to pancreatic ductal adenocarcinoma (PDAC) and chronic pancreatitis (CP) is crucial for early detection, treatment, and prevention. Methods: Association analyses of 10 serological biomarkers involved in cell signalling (IFN-γ, IL-6, IL-8, IL-10), oxidative stress (superoxide dismutase (SOD) and glutathione peroxidase (GPx) enzyme activities, total glutathione (GSH), malondialdehyde (MDA) levels), and intestinal permeability proteins (zonulin, I-FABP2) were conducted across PDAC (n = 12), CP (n = 21) and control subjects (n = 23). A Mendelian randomisation (MR) approach was used to assess causality of the identified significant associations in two large genetic cohorts (FinnGen and UK Biobank). Results: Observational results showed a downregulation of SOD and GPx antioxidant enzyme activities in PDAC and CP patients, respectively, and higher MDA levels in CP patients. Logistic regression models revealed significant associations between CP and SOD activity (OR = 0.21, 95% CI [0.05, 0.89], per SD), GPx activity (OR = 0.28, 95% CI [0.10, 0.79], per SD), and MDA levels (OR = 2.05, 95% CI [1.36, 3.08], per SD). MR analyses, however, did not support causality. Conclusions: These findings would not support oxidative stress-related biomarkers as potential targets for pancreatic diseases prevention. Yet, further research is encouraged to assess their viability as non-invasive tools for early diagnosis, particularly in pre-diagnostic CP populations.

Role of plant relatedness in plant-soil feedback dynamics of sympatric Asclepias species.

Plants affect associated biotic and abiotic edaphic factors, with reciprocal feedbacks from soil characteristics affecting plants. These two-way interactions between plants and soils are collectively known as plant-soil feedbacks (PSFs). The role of phylogenetic relatedness and evolutionary histories have recently emerged as a potential driver of PSFs, although the strength and direction of feedbacks among sympatric congeners are not well-understood. We examined plant-soil feedback responses of Asclepias syriaca, a common clonal milkweed species, with several sympatric congeners across a gradient of increasing phylogenetic distances (A. tuberosa, A. viridis, A. sullivantii, and A. verticillata, respectively). Plant-soil feedbacks were measured through productivity and colonization by arbuscular mycorrhizal (AM) fungi. Asclepias syriaca produced less biomass in soils conditioned by the most phylogenetically distant species (A. verticillata), relative to conspecific-conditioned soils. Similarly, arbuscular mycorrhizal (AM) fungal colonization of A. syriaca roots was reduced when grown in soils conditioned by A. verticillata, compared with colonization in plants grown in soil conditioned by any of the other three Asclepias species, indicating mycorrhizal associations are a potential mechanism of observed positive PSFs. This display of differences between the most phylogenetically distant, but not close or intermediate, paring(s) suggests a potential phylogenetic threshold, although other exogenous factors cannot be ruled out. Overall, these results highlight the potential role of phylogenetic distance in influencing positive PSFs through mutualists. The role of phylogenetic relatedness and evolutionary histories have recently emerged as a potential driver of plant-soil feedbacks (PSFs), although the strength and direction of feedbacks among sympatric congeners are not well-understood. Congeneric, sympatric milkweeds typically generated positive PSFs in terms of productivity and AM fungal colonization, suggesting the low likelihood of coexistence among tested pairs, with a strength of feedback increasing as the phylogenetic distance increases.

Meta-Analysis and Validation of a Colorectal Cancer Risk Prediction Model Using Deep Sequenced Fecal Metagenomes.

The gut microbiome is a potential modifiable risk factor for colorectal cancer (CRC). We re-analyzed all eight previously published stool sequencing data and conducted an MWAS meta-analysis. We used cross-validated LASSO predictive models to identify a microbiome signature for predicting the risk of CRC and precancerous lesions. These models were validated in a new study, Colorectal Cancer Screening (COLSCREEN), including 156 participants that were recruited in a CRC screening context. The MWAS meta-analysis identified 95 bacterial species that were statistically significantly associated with CRC (FDR < 0.05). The LASSO CRC predictive model obtained an area under the receiver operating characteristic curve (aROC) of 0.81 (95%CI: 0.78−0.83) and the validation in the COLSCREEN dataset was 0.75 (95%CI: 0.66−0.84). This model selected a total of 32 species. The aROC of this CRC-trained model to predict precancerous lesions was 0.52 (95%CI: 0.41−0.63). We have identified a signature of 32 bacterial species that have a good predictive accuracy to identify CRC but not precancerous lesions, suggesting that the identified microbes that were enriched or depleted in CRC are merely a consequence of the tumor. Further studies should focus on CRC as well as precancerous lesions with the intent to implement a microbiome signature in CRC screening programs.

Effects of Commercial Arbuscular Mycorrhizal Inoculants on Plant Productivity and Intra-Radical Colonization in Native Grassland: Unintentional De-Coupling of a Symbiosis?

There has been a surge in industries built on the production of arbuscular mycorrhizal (AM) fungal-based inoculants in the past few decades. This is not surprising, given the positive effects of AM fungi on plant growth and nutritional status. However, there is growing concern regarding the quality and efficacy of commercial inoculants. To assess the potential benefits and negative consequences of commercial AM fungal inoculants in grasslands, we conducted a controlled growth chamber study assessing the productivity and AM fungal root colonization of nine grassland plant species grown in grassland soil with or without one of six commercial AM fungal products. Our research showed no evidence of benefit; commercial inoculants never increased native plant biomass, although several inoculants decreased the growth of native species and increased the growth of invasive plant species. In addition, two commercial products contained excessive levels of phosphorus or nitrogen and consistently reduced AM fungal root colonization, indicating an unintentional de-coupling of the symbiosis. As there is little knowledge of the ecological consequences of inoculation with commercial AM fungal products, it is critical for restoration practitioners, scientists, and native plant growers to assess the presence of local AM fungal communities before investing in unnecessary, or possibly detrimental, AM fungal products.

Arbuscular Mycorrhizal Fungi Taxa Show Variable Patterns of Micro-Scale Dispersal in Prairie Restorations.

Human land use disturbance is a major contributor to the loss of natural plant communities, and this is particularly true in areas used for agriculture, such as the Midwestern tallgrass prairies of the United States. Previous work has shown that arbuscular mycorrhizal fungi (AMF) additions can increase native plant survival and success in plant community restorations, but the dispersal of AMF in these systems is poorly understood. In this study, we examined the dispersal of AMF taxa inoculated into four tallgrass prairie restorations. At each site, we inoculated native plant species with greenhouse-cultured native AMF taxa or whole soil collected from a nearby unplowed prairie. We monitored AMF dispersal, AMF biomass, plant growth, and plant community composition, at different distances from inoculation. In two sites, we assessed the role of plant hosts in dispersal, by placing known AMF hosts in a "bridge" and "island" pattern on either side of the inoculation points. We found that AMF taxa differ in their dispersal ability, with some taxa spreading to 2-m in the first year and others remaining closer to the inoculation point. We also found evidence that AMF spread altered non-inoculated neighboring plant growth and community composition in certain sites. These results represent the most comprehensive attempt to date to evaluate AMF spread.

The root-associated arbuscular mycorrhizal fungal assemblages of exotic alien plants are simplified in invaded distribution ranges, but dominant species are retained: A trans-continental perspective.

Arbuscular mycorrhizal fungi (AMF) provide crucial support for the establishment of plants in novel environments. We hypothesized that the OTU/genus richness and diversity of soil- and root-associated AMF associated with alien plant species in their exotic ranges are lower than those in their native ranges. We examined the root-associated and soil-dwelling AMF of 11 invasive plant species in their native and exotic ranges in the United States and Europe by DNA sequencing of the ITS2 locus. Examined root-associated AMF assemblages were simplified, which manifested as the loss of several AMF genera in the exotic ranges of the plants. These fungal assemblages were also characterized by greater dominance and simplification of the fungal assemblages. The dominant fungal genera were present regardless of whether their host plants were in their native or exotic ranges. Interestingly, both the native and invaded soils hosted diverse local AMF assemblages. Therefore, alien plant invasions were not limited to soils with low AMF diversity. Some AMF taxa could be context-dependent passengers rather than drivers of alien plant invasions. Further studies should identify functions of AMF missing or less abundant in roots of plants growing in exotic ranges.

Gene-gene interaction of AhRwith and within the Wntcascade affects susceptibility to lung cancer.

BACKGROUND: Aberrant Wnt signalling, regulating cell development and stemness, influences the development of many cancer types. The Aryl hydrocarbon receptor (AhR) mediates tumorigenesis of environmental pollutants. Complex interaction patterns of genes assigned to AhR/Wnt-signalling were recently associated with lung cancer susceptibility. AIM: To assess the association and predictive ability of AhR/Wnt-genes with lung cancer in cases and controls of European descent. METHODS: Odds ratios (OR) were estimated for genomic variants assigned to the Wnt agonist and the antagonistic genes DKK2, DKK3, DKK4, FRZB, SFRP4 and Axin2. Logistic regression models with variable selection were trained, validated and tested to predict lung cancer, at which other previously identified SNPs that have been robustly associated with lung cancer risk could also enter the model. Furthermore, decision trees were created to investigate variant × variant interaction. All analyses were performed for overall lung cancer and for subgroups. RESULTS: No genome-wide significant association of AhR/Wnt-genes with overall lung cancer was observed, but within the subgroups of ever smokers (e.g., maker rs2722278 SFRP4; OR = 1.20; 95% CI 1.13-1.27; p = 5.6 × 10-10) and never smokers (e.g., maker rs1133683 Axin2; OR = 1.27; 95% CI 1.19-1.35; p = 1.0 × 10-12). Although predictability is poor, AhR/Wnt-variants are unexpectedly overrepresented in optimized prediction scores for overall lung cancer and for small cell lung cancer. Remarkably, the score for never-smokers contained solely two AhR/Wnt-variants. The optimal decision tree for never smokers consists of 7 AhR/Wnt-variants and only two lung cancer variants. CONCLUSIONS: The role of variants belonging to Wnt/AhR-pathways in lung cancer susceptibility may be underrated in main-effects association analysis. Complex interaction patterns in individuals of European descent have moderate predictive capacity for lung cancer or subgroups thereof, especially in never smokers.

Accounting for EGFR Mutations in Epidemiologic Analyses of Non-Small Cell Lung Cancers: Examples Based on the International Lung Cancer Consortium Data.

BACKGROUND: Somatic EGFR mutations define a subset of non-small cell lung cancers (NSCLC) that have clinical impact on NSCLC risk and outcome. However, EGFR-mutation-status is often missing in epidemiologic datasets. We developed and tested pragmatic approaches to account for EGFR-mutation-status based on variables commonly included in epidemiologic datasets and evaluated the clinical utility of these approaches. METHODS: Through analysis of the International Lung Cancer Consortium (ILCCO) epidemiologic datasets, we developed a regression model for EGFR-status; we then applied a clinical-restriction approach using the optimal cut-point, and a second epidemiologic, multiple imputation approach to ILCCO survival analyses that did and did not account for EGFR-status. RESULTS: Of 35,356 ILCCO patients with NSCLC, EGFR-mutation-status was available in 4,231 patients. A model regressing known EGFR-mutation-status on clinical and demographic variables achieved a concordance index of 0.75 (95% CI, 0.74-0.77) in the training and 0.77 (95% CI, 0.74-0.79) in the testing dataset. At an optimal cut-point of probability-score = 0.335, sensitivity = 69% and specificity = 72.5% for determining EGFR-wildtype status. In both restriction-based and imputation-based regression analyses of the individual roles of BMI on overall survival of patients with NSCLC, similar results were observed between overall and EGFR-mutation-negative cohort analyses of patients of all ancestries. However, our approach identified some differences: EGFR-mutated Asian patients did not incur a survival benefit from being obese, as observed in EGFR-wildtype Asian patients. CONCLUSIONS: We introduce a pragmatic method to evaluate the potential impact of EGFR-status on epidemiological analyses of NSCLC. IMPACT: The proposed method is generalizable in the common occurrence in which EGFR-status data are missing.

Hepcidin-regulating iron metabolism genes and pancreatic ductal adenocarcinoma: a pathway analysis of genome-wide association studies.

BACKGROUND: Epidemiological studies have suggested positive associations for iron and red meat intake with risk of pancreatic ductal adenocarcinoma (PDAC). Inherited pathogenic variants in genes involved in the hepcidin-regulating iron metabolism pathway are known to cause iron overload and hemochromatosis. OBJECTIVES: The objective of this study was to determine whether common genetic variation in the hepcidin-regulating iron metabolism pathway is associated with PDAC. METHODS: We conducted a pathway analysis of the hepcidin-regulating genes using single nucleotide polymorphism (SNP) summary statistics generated from 4 genome-wide association studies in 2 large consortium studies using the summary data-based adaptive rank truncated product method. Our population consisted of 9253 PDAC cases and 12,525 controls of European descent. Our analysis included 11 hepcidin-regulating genes [bone morphogenetic protein 2 (BMP2), bone morphogenetic protein 6 (BMP6), ferritin heavy chain 1 (FTH1), ferritin light chain (FTL), hepcidin (HAMP), homeostatic iron regulator (HFE), hemojuvelin (HJV), nuclear factor erythroid 2-related factor 2 (NRF2), ferroportin 1 (SLC40A1), transferrin receptor 1 (TFR1), and transferrin receptor 2 (TFR2)] and their surrounding genomic regions (±20 kb) for a total of 412 SNPs. RESULTS: The hepcidin-regulating gene pathway was significantly associated with PDAC (P = 0.002), with the HJV, TFR2, TFR1, BMP6, and HAMP genes contributing the most to the association. CONCLUSIONS: Our results support that genetic susceptibility related to the hepcidin-regulating gene pathway is associated with PDAC risk and suggest a potential role of iron metabolism in pancreatic carcinogenesis. Further studies are needed to evaluate effect modification by intake of iron-rich foods on this association.

Tumor-Associated Microbiome: Where Do We Stand?

The study of the human microbiome in oncology is a growing and rapidly evolving field. In the past few years, there has been an exponential increase in the number of studies investigating associations of microbiome and cancer, from oncogenesis and cancer progression to resistance or sensitivity to specific anticancer therapies. The gut microbiome is now known to play a significant role in antitumor immune responses and in predicting the efficacy of immune-checkpoint inhibitors in cancer patients. Beyond the gut, the tumor-associated microbiome-microbe communities located either in the tumor or within its body compartment-seems to interact with the local microenvironment and the tumor immune contexture, ultimately impacting cancer progression and treatment outcome. However, pre-clinical research focusing on causality and mechanistic pathways as well as proof-of-concept studies are still needed to fully understand the potential clinical utility of microbiome in cancer patients. Moreover, there is a need for the standardization of methodology and the implementation of quality control across microbiome studies to allow for a better interpretation and greater comparability of the results reported between them. This review summarizes the accumulating evidence in the field and discusses the current and upcoming challenges of microbiome studies.

Smoking Modifies Pancreatic Cancer Risk Loci on 2q21.3.

Germline variation and smoking are independently associated with pancreatic ductal adenocarcinoma (PDAC). We conducted genome-wide smoking interaction analysis of PDAC using genotype data from four previous genome-wide association studies in individuals of European ancestry (7,937 cases and 11,774 controls). Examination of expression quantitative trait loci data from the Genotype-Tissue Expression Project followed by colocalization analysis was conducted to determine whether there was support for common SNP(s) underlying the observed associations. Statistical tests were two sided and P < 5 × 10-8 was considered statistically significant. Genome-wide significant evidence of qualitative interaction was identified on chr2q21.3 in intron 5 of the transmembrane protein 163 (TMEM163) and upstream of the cyclin T2 (CCNT2). The most significant SNP using the Empirical Bayes method, in this region that included 45 significantly associated SNPs, was rs1818613 [per allele OR in never smokers 0.87, 95% confidence interval (CI), 0.82-0.93; former smokers 1.00, 95% CI, 0.91-1.07; current smokers 1.25, 95% CI 1.12-1.40, P interaction = 3.08 × 10-9). Examination of the Genotype-Tissue Expression Project data demonstrated an expression quantitative trait locus in this region for TMEM163 and CCNT2 in several tissue types. Colocalization analysis supported a shared SNP, rs842357, in high linkage disequilibrium with rs1818613 (r 2 = 0. 94) driving both the observed interaction and the expression quantitative trait loci signals. Future studies are needed to confirm and understand the differential biologic mechanisms by smoking status that contribute to our PDAC findings. SIGNIFICANCE: This large genome-wide interaction study identifies a susceptibility locus on 2q21.3 that significantly modified PDAC risk by smoking status, providing insight into smoking-associated PDAC, with implications for prevention.

Integration of multiomic annotation data to prioritize and characterize inflammation and immune-related risk variants in squamous cell lung cancer.

Clinical trial results have recently demonstrated that inhibiting inflammation by targeting the interleukin-1ß pathway can offer a significant reduction in lung cancer incidence and mortality, highlighting a pressing and unmet need to understand the benefits of inflammation-focused lung cancer therapies at the genetic level. While numerous genome-wide association studies (GWAS) have explored the genetic etiology of lung cancer, there remains a large gap between the type of information that may be gleaned from an association study and the depth of understanding necessary to explain and drive translational findings. Thus, in this study we jointly model and integrate extensive multiomics data sources, utilizing a total of 40 genome-wide functional annotations that augment previously published results from the International Lung Cancer Consortium (ILCCO) GWAS, to prioritize and characterize single nucleotide polymorphisms (SNPs) that increase risk of squamous cell lung cancer through the inflammatory and immune responses. Our work bridges the gap between correlative analysis and translational follow-up research, refining GWAS association measures in an interpretable and systematic manner. In particular, reanalysis of the ILCCO data highlights the impact of highly associated SNPs from nuclear factor-κB signaling pathway genes as well as major histocompatibility complex mediated variation in immune responses. One consequence of prioritizing likely functional SNPs is the pruning of variants that might be selected for follow-up work by over an order of magnitude, from potentially tens of thousands to hundreds. The strategies we introduce provide informative and interpretable approaches for incorporating extensive genome-wide annotation data in analysis of genetic association studies.

Incidence and Survival Trends of Pancreatic Cancer in Girona: Impact of the Change in Patient Care in the Last 25 Years.

(1) Background: We investigated the incidence and survival trends for pancreatic cancer (PC) over the last 25 years in the Girona region, Catalonia, Spain; (2) Methods: Data were extracted from the population-based Girona Cancer Registry. Incident PC cases during 1994-2015 were classified using the International Classification of Diseases for Oncology Third Edition (ICD-O-3). Incidence rates age-adjusted to the European standard population (ASRE) and world standard population (ASRW) were obtained. Trends were assessed using the estimated annual percentage of change (EAPC) of the ASRE13. Observed and relative survivals (RS) were estimated with the Kaplan-Meier and Pohar Perme methods, respectively; (3) Results: We identified 1602 PC incident cases. According to histology, 44.4% of cases were exocrine PC, 4.1% neuroendocrine, and 51.1% malignant-non-specified. The crude incidence rate (CR) for PC was 11.43 cases-per-100,000 inhabitants/year. A significant increase of incidence with age and over the study period was observed. PC overall 5-year RS was 7.05% (95% confidence interval (CI) 5.63; 8.84). Longer overall survival was observed in patients with neuroendocrine tumours (5-year RS 61.45%; 95% CI 47.47; 79.55). Trends in 5-year RS for the whole cohort rose from 3.27% (95% CI 1.69-6.35) in 1994-1998 to 13.1% (95% CI 9.98; 17.2) in 2010-2015; (4) Conclusions: Incidence rates of PC in Girona have increased in the last two decades. There is a moderate but encouraging increase in survival thorough the study period. These results can be used as baseline for future research.

Identifying environmental risk factors for inflammatory bowel diseases: a Mendelian randomization study.

Several studies have examined environmental factors and inflammatory bowel diseases (IBD) using traditional approaches; however, provided results are still conflicting. Our aim was to determine whether lifestyle and nutrient exposures, related to IBD in observational meta-analyses, influence IBD risk using a Mendelian randomization (MR) approach. A two-sample MR approach was applied on summary-level genome-wide association results. Genetic variants strongly associated with measures of tobacco smoking, obesity and fat distribution, physical activity, and blood levels of vitamins and fatty acids were evaluated on genetic data from international IBD consortia including a total of 25,042 IBD cases (12,194 cases of Crohn's disease (CD) and 12,366 cases of ulcerative colitis (UC)) and 34,915 controls. Our results indicated that, among lifestyle exposures, being a smoker was positively associated with CD (OR 1.13, P = 0.02), but it was not associated with UC risk (OR 0.99, P = 0.88). Body-mass index (BMI) and body fat percentage were positively associated with CD (OR 1.11, P = 0.02, per standard deviation (SD) of 4.6 kg/m2; and OR 1.50, P = 3 × 10-10, per SD of 6.6%; respectively); while for UC, BMI was inversely associated (OR 0.85, P = 5 × 10-5; per SD) and body fat percentage showed a OR of 1.11 (P = 0.11; per SD). Additionally, among nutrient exposures, omega-3 fatty acids levels were inversely associated with CD (OR 0.67, P = 2 × 10-6). Our MR results did not support a protective effect for being a smoker on UC risk; however, they are compatible with a risk effect for higher body fat proportion and a protective role for higher levels of omega-3 fatty acids on CD etiology.

A new efficient method to detect genetic interactions for lung cancer GWAS.

BACKGROUND: Genome-wide association studies (GWAS) have proven successful in predicting genetic risk of disease using single-locus models; however, identifying single nucleotide polymorphism (SNP) interactions at the genome-wide scale is limited due to computational and statistical challenges. We addressed the computational burden encountered when detecting SNP interactions for survival analysis, such as age of disease-onset. To confront this problem, we developed a novel algorithm, called the Efficient Survival Multifactor Dimensionality Reduction (ES-MDR) method, which used Martingale Residuals as the outcome parameter to estimate survival outcomes, and implemented the Quantitative Multifactor Dimensionality Reduction method to identify significant interactions associated with age of disease-onset. METHODS: To demonstrate efficacy, we evaluated this method on two simulation data sets to estimate the type I error rate and power. Simulations showed that ES-MDR identified interactions using less computational workload and allowed for adjustment of covariates. We applied ES-MDR on the OncoArray-TRICL Consortium data with 14,935 cases and 12,787 controls for lung cancer (SNPs = 108,254) to search over all two-way interactions to identify genetic interactions associated with lung cancer age-of-onset. We tested the best model in an independent data set from the OncoArray-TRICL data. RESULTS: Our experiment on the OncoArray-TRICL data identified many one-way and two-way models with a single-base deletion in the noncoding region of BRCA1 (HR 1.24, P = 3.15 × 10-15), as the top marker to predict age of lung cancer onset. CONCLUSIONS: From the results of our extensive simulations and analysis of a large GWAS study, we demonstrated that our method is an efficient algorithm that identified genetic interactions to include in our models to predict survival outcomes.

Incorporating multiple sets of eQTL weights into gene-by-environment interaction analysis identifies novel susceptibility loci for pancreatic cancer.

It is of great scientific interest to identify interactions between genetic variants and environmental exposures that may modify the risk of complex diseases. However, larger sample sizes are usually required to detect gene-by-environment interaction (G × E) than required to detect genetic main association effects. To boost the statistical power and improve the understanding of the underlying molecular mechanisms, we incorporate functional genomics information, specifically, expression quantitative trait loci (eQTLs), into a data-adaptive G × E test, called aGEw. This test adaptively chooses the best eQTL weights from multiple tissues and provides an extra layer of weighting at the genetic variant level. Extensive simulations show that the aGEw test can control the Type 1 error rate, and the power is resilient to the inclusion of neutral variants and noninformative external weights. We applied the proposed aGEw test to the Pancreatic Cancer Case-Control Consortium (discovery cohort of 3,585 cases and 3,482 controls) and the PanScan II genome-wide association study data (replication cohort of 2,021 cases and 2,105 controls) with smoking as the exposure of interest. Two novel putative smoking-related pancreatic cancer susceptibility genes, TRIP10 and KDM3A, were identified. The aGEw test is implemented in an R package aGE.

Gallbladder disease, cholecystectomy, and pancreatic cancer risk in the International Pancreatic Cancer Case-Control Consortium (PanC4).

BACKGROUND: The association among gallbladder disease, cholecystectomy, and pancreatic cancer is unclear. Moreover, time interval between gallbladder disease or cholecystectomy and pancreatic cancer diagnosis is not considered in most previous studies. AIM: To quantify the association among gallbladder disease, cholecystectomy, and pancreatic cancer, considering time since first diagnosis of gallbladder disease or cholecystectomy. METHODS: We used data from nine case-control studies within the Pancreatic Cancer Case-Control Consortium, including 5760 cases of adenocarcinoma of the exocrine pancreas and 8437 controls. We estimated pooled odds ratios and the corresponding 95% confidence intervals by estimating study-specific odds ratios through multivariable unconditional logistic regression models, and then pooling the obtained estimates using fixed-effects models. RESULTS: Compared with patients with no history of gallbladder disease, the pooled odds ratio of pancreatic cancer was 1.69 (95% confidence interval, 1.51-1.88) for patients reporting a history of gallbladder disease. The odds ratio was 4.90 (95% confidence interval, 3.45-6.97) for gallbladder disease diagnosed <2 years before pancreatic cancer diagnosis and 1.11 (95% confidence interval, 0.96-1.29) when ≥2 years elapsed. The pooled odds ratio was 1.64 (95% confidence interval, 1.43-1.89) for patients who underwent cholecystectomy, as compared to those without cholecystectomy. The odds ratio was 7.00 (95% confidence interval, 4.13-11.86) for a surgery <2 years before pancreatic cancer diagnosis and 1.28 (95% confidence interval, 1.08-1.53) for a surgery ≥2 years before. CONCLUSIONS: There appears to be no long-term effect of gallbladder disease on pancreatic cancer risk, and at most a modest one for cholecystectomy. The strong short-term association can be explained by diagnostic bias and reverse causation.

Genome-Wide Gene-Diabetes and Gene-Obesity Interaction Scan in 8,255 Cases and 11,900 Controls from PanScan and PanC4 Consortia.

BACKGROUND: Obesity and diabetes are major modifiable risk factors for pancreatic cancer. Interactions between genetic variants and diabetes/obesity have not previously been comprehensively investigated in pancreatic cancer at the genome-wide level. METHODS: We conducted a gene-environment interaction (GxE) analysis including 8,255 cases and 11,900 controls from four pancreatic cancer genome-wide association study (GWAS) datasets (Pancreatic Cancer Cohort Consortium I-III and Pancreatic Cancer Case Control Consortium). Obesity (body mass index ≥30 kg/m2) and diabetes (duration ≥3 years) were the environmental variables of interest. Approximately 870,000 SNPs (minor allele frequency ≥0.005, genotyped in at least one dataset) were analyzed. Case-control (CC), case-only (CO), and joint-effect test methods were used for SNP-level GxE analysis. As a complementary approach, gene-based GxE analysis was also performed. Age, sex, study site, and principal components accounting for population substructure were included as covariates. Meta-analysis was applied to combine individual GWAS summary statistics. RESULTS: No genome-wide significant interactions (departures from a log-additive odds model) with diabetes or obesity were detected at the SNP level by the CC or CO approaches. The joint-effect test detected numerous genome-wide significant GxE signals in the GWAS main effects top hit regions, but the significance diminished after adjusting for the GWAS top hits. In the gene-based analysis, a significant interaction of diabetes with variants in the FAM63A (family with sequence similarity 63 member A) gene (significance threshold P < 1.25 × 10-6) was observed in the meta-analysis (P GxE = 1.2 ×10-6, P Joint = 4.2 ×10-7). CONCLUSIONS: This analysis did not find significant GxE interactions at the SNP level but found one significant interaction with diabetes at the gene level. A larger sample size might unveil additional genetic factors via GxE scans. IMPACT: This study may contribute to discovering the mechanism of diabetes-associated pancreatic cancer.

Associations between Genetically Predicted Blood Protein Biomarkers and Pancreatic Cancer Risk.

BACKGROUND: Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal malignancies, with few known risk factors and biomarkers. Several blood protein biomarkers have been linked to PDAC in previous studies, but these studies have assessed only a limited number of biomarkers, usually in small samples. In this study, we evaluated associations of circulating protein levels and PDAC risk using genetic instruments. METHODS: To identify novel circulating protein biomarkers of PDAC, we studied 8,280 cases and 6,728 controls of European descent from the Pancreatic Cancer Cohort Consortium and the Pancreatic Cancer Case-Control Consortium, using genetic instruments of protein quantitative trait loci. RESULTS: We observed associations between predicted concentrations of 38 proteins and PDAC risk at an FDR of < 0.05, including 23 of those proteins that showed an association even after Bonferroni correction. These include the protein encoded by ABO, which has been implicated as a potential target gene of PDAC risk variant. Eight of the identified proteins (LMA2L, TM11D, IP-10, ADH1B, STOM, TENC1, DOCK9, and CRBB2) were associated with PDAC risk after adjusting for previously reported PDAC risk variants (OR ranged from 0.79 to 1.52). Pathway enrichment analysis showed that the encoding genes for implicated proteins were significantly enriched in cancer-related pathways, such as STAT3 and IL15 production. CONCLUSIONS: We identified 38 candidates of protein biomarkers for PDAC risk. IMPACT: This study identifies novel protein biomarker candidates for PDAC, which if validated by additional studies, may contribute to the etiologic understanding of PDAC development.

Menstrual Factors, Reproductive History, Hormone Use, and Urothelial Carcinoma Risk: A Prospective Study in the EPIC Cohort.

BACKGROUND: Urothelial carcinoma is the predominant (95%) bladder cancer subtype in industrialized nations. Animal and epidemiologic human studies suggest that hormonal factors may influence urothelial carcinoma risk. METHODS: We used an analytic cohort of 333,919 women from the European Prospective Investigation into Cancer and Nutrition Cohort. Associations between hormonal factors and incident urothelial carcinoma (overall and by tumor grade, tumor aggressiveness, and non-muscle-invasive urothelial carcinoma) risk were evaluated using Cox proportional hazards models. RESULTS: During a mean of 15 years of follow-up, 529 women developed urothelial carcinoma. In a model including number of full-term pregnancies (FTP), menopausal status, and menopausal hormone therapy (MHT), number of FTP was inversely associated with urothelial carcinoma risk (HR≥5vs1 = 0.48; 0.25-0.90; P trend in parous women = 0.010) and MHT use (compared with nonuse) was positively associated with urothelial carcinoma risk (HR = 1.27; 1.03-1.57), but no dose response by years of MHT use was observed. No modification of HRs by smoking status was observed. Finally, sensitivity analyses in never smokers showed similar HR patterns for the number of FTP, while no association between MHT use and urothelial carcinoma risk was observed. Association between MHT use and urothelial carcinoma risk remained significant only in current smokers. No heterogeneity of the risk estimations in the final model was observed by tumor aggressiveness or by tumor grade. A positive association between MTH use and non-muscle-invasive urothelial carcinoma risk was observed. CONCLUSIONS: Our results support that increasing the number of FTP may reduce urothelial carcinoma risk. IMPACT: More detailed studies on parity are needed to understand the possible effects of perinatal hormone changes in urothelial cells.