- Research article
- Open Access
Role of biliary stent and neoadjuvant chemotherapy in the pancreatic tumor microbiome
BMC Microbiology volume 21, Article number: 280 (2021)
Intra-tumor microbiota have been implicated in pancreatic ductal adenocarcinoma (PDAC) development, treatment response and post-treatment survivorship. Moreover, therapeutic interventions targeting microbiota may improve the response to chemotherapy and immunotherapy, further emphasizing the critical need to understand the origins of and growth of bacteria within the pancreatic tumor microenvironment. Here, we studied the role of several clinical factors on the bacterial colonization of PDAC.
We obtained matched tumor and normal pancreatic tissue specimens from 27 patients who had undergone surgical resection for PDAC between 2011 and 2015 from the University of Minnesota Biological Materials Procurement Network (BioNet). We found that 26 (48%) out of 54 pancreatic tissue samples harbored detectable bacterial communities using real-time PCR targeting the 16S rRNA gene. Bacterial colonization was detected significantly more frequently in samples from patients who had pancreatic head tumors, underwent Whipple procedure, or had preoperative biliary stent placement. There was also a significantly greater relative abundance of microbiota from the family Enterobacteriaceae among samples from patients who underwent biliary stent placement or neoadjuvant treatment with a combination of Gemcitabine and Paclitaxel.
These findings suggest that biliary stent placement and neoadjuvant chemotherapy are associated with specific alterations that promote the infiltration and growth of intra-tumor bacteria in the setting of PDAC. Further studies exploring whether specific bacterial communities could contribute to increased chemoresistance will be essential for optimizing medical therapies in the future.
Pancreatic ductal adenocarcinoma (PDAC) has an average five-year survival of 9% across all stages and is the fourth leading cause of cancer death in men and women in the United States . Fewer than 20% of patients present with resectable tumors, and most patients develop disease recurrence despite curative-intent surgery and adjuvant chemotherapy, with or without chemoradiation . A lack of early diagnostic markers, rapid disease progression, and poor sensitivity to adjuvant therapies contribute to this dismal prognosis. Therefore, novel strategies to advance current therapeutic options are needed.
Once considered a sterile organ, recent evidence demonstrates that the pancreas harbors a microbiota, which is even more abundant in pancreatic cancer tissue . The presence of intra-tumor microbiota in pancreatic cancer has many implications. The tumor microbiome has been linked to PDAC oncogenesis and cancer-associated inflammation . Tumor microbial signatures may be indicative of the likelihood of PDAC survival, as a distinct bacterial composition with greater diversity has been found in tumors from long-term survivors (LTS) compared to short-term survivors (STS) . In fact, transfer of gut microbiota from LTS induced an antitumor response in tumor-bearing mice that was not observed after transfer of microbiota from STS. In an orthotopic murine model of PDAC, gut bacterial ablation with antibiotics led to immunogenic reprogramming of the tumor microenvironment and improved the efficacy of immunotherapy . Furthermore, a landmark study by Geller et al. revealed that intra-tumor Gammaproteobacteria can metabolize the chemotherapy agent Gemcitabine, suggesting a potential role in mediating chemoresistance in PDAC .
Bacterial colonization of pancreatic cancer tissue is thought to be due to the migration of microbiota from oral, gastrointestinal, and hepatobiliary sources [7,8,9,10,11]. Routes of colonization of intestinal bacteria are suspected to include the bloodstream, lymphatic system, and direct reflux from the duodenum . Interventions such as preoperative biliary drainage and neoadjuvant chemotherapy have been shown to alter biliary and duodenal microbiota, and pancreatic head carcinomas have been strongly correlated with bactibilia [13,14,15,16]. Yet, the influence of these clinical factors on the infiltration and growth of bacteria within the pancreatic tumor microenvironment is poorly characterized. Geller et al. began to investigate this and found that patients who underwent pancreatic duct instrumentation had more bacteria in their pancreas tumors than those who had not undergone instrumentation . However, further research to define the impact of these clinical factors on the bacterial colonization of pancreatic tumor tissue is imperative. In this study, we performed a retrospective, exploratory analysis of pancreatic tissue samples from patients who underwent surgical resection for PDAC. We hypothesized that clinical parameters such as biliary stent placement and neoadjuvant chemotherapy likely contribute to intra-tumor bacterial colonization.
Matched malignant pancreatic tissue and normal adjacent pancreatic tissue samples were obtained from 27 patients (n = 54 samples) who underwent surgical resection for PDAC at the University of Minnesota between 2011 and 2015. Patient demographic and clinical characteristics are presented in Table 1. This cohort consisted of 18 men (67%) and 9 women (33%) with a mean age of 64 years and mean body mass index (BMI) of 29 kg/m2. Tumors were most commonly located in the head of the pancreas (n = 18, 67%) and classified as stage IIB (n = 14, 52%). A preoperative biliary stent was placed in 16 (59%) patients, with similar frequency in choice of plastic or metal stent material. Neoadjuvant chemotherapy was administered in 10 (37%) patients, and choice of specific chemotherapy regimen was left to the discretion of the treating oncologist. Few patients had antibiotic exposures within 6 weeks pre-operatively (n = 3, 11%); one treated for cholangitis until 10 days prior to surgery, one treated for acute cholecystitis until 14 days prior to surgery, and one treated for acute diverticulitis until 5 weeks prior to surgery.
A majority of patients underwent Whipple procedure (n = 18, 67%). Choice of perioperative antibiotic prophylaxis was based on surgeon preference, and the majority of patients received either Cefazolin (n = 12, 44%), Cefoxitin (n = 7, 26%), or Cefotetan (n = 5, 18%). Levofloxacin (n = 1, 4%) and Clindamycin (n = 1, 4%) were given to patients with penicillin allergy. Infectious postoperative complications occurred in 8 (30%) patients. Most patients were short-term survivors who lived less than 5 years post-operatively (n = 17, 63%), and few were long-term survivors who lived greater than 5 years post-operatively (n = 3, 11%).
Bacterial colonization associated with clinical features
We investigated whether pancreatic tissue samples harbored detectable bacterial DNA via real-time PCR targeting the bacterial 16S rRNA gene. Out of the 27 patients examined, bacteria were detected in both matched samples from 8 (29.6%) patients, only malignant tissue sample from 5 (18.5%) patients, only normal adjacent tissue sample from 5 (18.5%) patients, and neither matched sample in 9 (33.3%) patients. To understand which clinical parameters may contribute to the presence of microbiota on PDAC tumor tissue, we compared malignant pancreatic tissue samples with (n = 13, 48%) and without (n = 14, 52%) bacterial colonization at the time of surgery (Table 2).
No significant differences in bacterial colonization were observed based on sex, age, BMI, tumor stage, or receipt of neoadjuvant chemotherapy. Bacterial colonization occurred more frequently in samples from patients with tumors located in the head of the pancreas (χ2 = 6.0, P = 0.002) and those that required surgical resection with Whipple procedure (χ2 = 3.8, P < 0.001) compared to samples without bacterial colonization. Samples with bacterial colonization were also significantly more likely to be from patients who underwent preoperative biliary stent placement (χ2 = 3.8, P = 0.001) in contrast to samples without bacterial colonization. Yet, no difference based on type of biliary stent (χ2 = 6.0, P = 0.565) was observed. Nor were there significant differences in bacterial colonization associated with type of perioperative antibiotic prophylaxis, postoperative infectious complications, or short- versus long-term survival. This analysis was also performed with normal adjacent tissue and no significant differences in bacterial colonization were identified based on these clinical factors.
Enterobacteriaceae linked to stent placement and neoadjuvant chemotherapy
Next, we performed 16S rRNA gene-based amplicon sequencing of the 26 samples that harbored bacterial signatures, as determined by real-time PCR, to characterize taxonomic profiles. Among both malignant and normal adjacent tissue samples, the most abundant families of bacteria included Ruminococcaceae, Staphylococcaceae, Bacillaceae, Enterobacteriaceae, and Pseudomonadaceae (Fig. 1A-B). There were no significant differences in overall community composition between malignant and normal adjacent pancreatic tissue types (Fig. 1C, ANOSIM R = − 0.01, P = 0.60). There also were no significant differences in community composition based on receipt of preoperative biliary stent or neoadjuvant chemotherapy (R = 0.12 and 0.14, P = 0.23 and 0.13, respectively).
Despite similar overall community composition, differences in the relative abundances of members of Enterobacteriaceae, a predominant family of microbiota among samples, were observed based on receipt of preoperative interventions (Fig. 2). There was significantly greater relative abundance of Enterobacteriaceae in samples from patients who underwent biliary stent placement (Dunn’s post-hoc P = 0.004). Additionally, samples from patients who underwent neoadjuvant treatment with a combination of Gemcitabine and Paclitaxel had significantly greater relative abundances of Enterobacteriaceae compared with samples from patients who underwent Gemcitabine therapy alone (P = 0.04), or patients who did not receive neoadjuvant chemotherapy at all (P < 0.001).
Resistance to antitumor drug therapies remains a significant challenge in the management of pancreatic cancer . Tumor-associated microbiota have been implicated in reducing the efficacy of chemotherapeutics and could potentially serve as a novel target to optimize current medical therapies [7, 18]. In this study, we aimed to further understand the complex relationship between microbiota and PDAC by exploring clinical factors that may contribute to intra-tumor bacterial colonization.
Bacterial colonization of pancreatic tumor tissue was more likely in patients who had pancreatic head tumors and those that underwent Whipple procedure. This could be due to the fact that pancreatic head tumors, which often require Whipple procedure for resection, are closer in proximity to the gastrointestinal and hepatobiliary systems. Tumors at this location are also more likely to result in malignant biliary obstruction than distal tumors.
Intra-tumor microbial signatures were also more detectable in patients who had undergone preoperative biliary drainage with stent placement. Importantly, we found that biliary stent placement was correlated with greater relative abundances of Enterobacteriaceae, which include bacteria that have been postulated to confer resistance to the chemotherapeutic agent Gemcitabine . As stated previously, Scheufele, et al. reported that stent placement affected the composition of biliary microbiota . Migration of microbes from the biliary tract, adjacent to the pancreas, may influence the growth of intra-tumor microbiota. These investigators also postulated that plastic stents could lead to greater alterations in biliary microbiota because plastic stents have greater risk of recurrent biliary obstruction and poorer duration of patency compared to metal stents. Interestingly, we found that the type of biliary stent used did not impact the growth of intra-tumor bacteria.
We found that the multi-agent chemotherapy regimen comprised of Gemcitabine and Paclitaxel was also associated with greater relative abundances of Enterobacteriaceae compared to other chemotherapy regimens and no chemotherapy at all among our samples. This could be due to immunosuppression from the chemotherapy agent itself. Patients who undergo neoadjuvant therapy may additionally experience prolonged periods of stenting and increased exposure to antibiotics for stent-related cholangitis . Future studies investigating the effects of neoadjuvant chemotherapy on PDAC microbiota should include multi-agent regimens as they may have variable outcomes.
Our findings limit definitive conclusions due to the retrospective fashion of data retrieval. Pancreatic tissue samples were difficult to obtain, and the small sample size of this study restricted the ability to perform well powered multivariate statistics. This may introduce bias from the lack of adjustment for potentially confounding factors, such as choice of preoperative antibiotic. We also did not have access to pancreatic ductal fluid samples or associated culture data, which may have provided additional insight into microbial changes in the tumor microenvironment. We chose to focus on family-level taxonomic assignment due to the uncertainty in classification at greater resolution due to low bacterial biomass. We did not explore the role of viral or fungal microbiota, which have also been implicated in PDAC oncogenesis [7, 19]. Additionally, the loss to follow-up of 26% of patients limits interpretation of findings related to survival. Nevertheless, the impact of several clinical factors in the growth of intra-tumor bacteria in PDAC observed in this study is novel. A better understanding of the role of intra-tumor microbiota in PDAC may aid in identifying new approaches in microbial modulation to optimize current chemotherapy options and improve patient outcomes.
Preoperative biliary stent placement and neoadjuvant chemotherapy are associated with specific alterations in the pancreatic tumor microenvironment that may encourage bacterial colonization of PDAC tissue. Particularly, biliary stent placement and chemotherapy regimen of Gemcitabine and Paclitaxel may selectively promote the growth of Enterobacteriaceae, which have been postulated to contribute to chemoresistance. Future studies addressing whether these compositional alterations could contribute to increased chemoresistance may provide an opportunity to enhance the efficacy of medical therapies for pancreatic cancer.
Human sample collection
Matched tumor and normal adjacent pancreatic tissue samples were collected from 27 patients who had undergone surgical resection for PDAC between 2011 and 2015 from the University of Minnesota Biological Materials Procurement Network (BioNet). Surgical specimens were initially removed from patients intraoperatively by the surgeon and sent fresh to the pathologist. After gross examination, malignant and normal adjacent pancreatic tissue samples were sliced from the single surgical specimen. They were placed in clean petri dishes and given to BioNet, where they were frozen in liquid nitrogen in separate cryovials. Tissue handling by BioNet was performed in a clean environment, where there is possibility of contamination without complete sterility.
Normal adjacent pancreatic tissue was obtained from the site most distant from the tumor, with exact distance from the tumor depending on the size of the tumor and the entire specimen. Histologic quality control assessment was performed, and all normal adjacent tissue samples had 0% neoplastic tissue. Samples without histologic quality control assessment were excluded. Human specimens were obtained under approval by the Institutional Review Board (IRB). All specimens were stored at − 80 °C until further use.
DNA extraction and real-time PCR
DNA was extracted from pancreatic tissue samples using the DNeasy PowerSoil kit (QIAGEN, Hilden, Germany) in accordance with manufacturer instructions. Total bacterial DNA in pancreatic tissue samples was determined by non-quantitative, real-time PCR. The V4 hypervariable region of the 16S rRNA gene was amplified using forward primer 515F (5′- GTGYCAGCMGCCGCGGTAA − 3′) and reverse primer 806R (5′- GGACTACNVGGGTWTCTAAT − 3′). Briefly, the reaction mix contained QuantiTect SYBR Green PCR Master Mix (QIAGEN), 0.5 μM forward and reverse primers, nuclease-free water, and sample DNA. PCR cycling conditions were: initial denaturation at 95 °C for 15 min, followed by 40 cycles of 30 s at 95 °C, 30 s at 50 °C, and 30 s at 72 °C. A standard curve using the V4 gene fragment cloned into a GenBlocks vector (IDT DNA, Coralville, IA) was run at concentrations from 5 × 106 to 50 gene copies/reaction. Negative control (sterile water blank) reactions were also included. Reactions were run in triplicate on the LightCycler 480 Instrument II (Roche Diagnostics, Ltd., Basel, Switzerland). Due to low bacterial biomass and samples amplifying at the lower end of the standard curve, quantitative results were deemed unreliable, and samples were considered positive if they amplified before cycle 30 and had a melting temperature within the range of the standards (85.4 ± 0.7 °C).
16S rRNA gene-based amplicon sequencing
The V4 hypervariable region of the 16S rRNA gene was amplified using the 515F/806R primer set  by the University of Minnesota Genomics Center (UMGC, Minneapolis, MN, USA), as previously described . Paired-end, dual indexed sequencing at read length of 300 nucleotides was performed on the Illumina MiSeq platform (Illumina, Inc., San Diego, CA, USA) by UMGC . Negative sterile water controls were included in all sequencing runs and did not produce amplicons. Raw data are stored in the Sequence Read Archive under BioProject accession number SRP197553.
Amplicon sequence data were processed and analyzed using mothur software ver 1.41.1  with our previously published pipeline for quality screening and taxonomic annotation . Reads were paired-end joined, quality trimmed, and aligned against the SILVA database ver. 132 . Operational taxonomic units (OTUs) were binned at 97% similarity using the furthest-neighbor algorithm and taxonomic assignment was done against the Ribosomal Database Project ver. 16 . A mean (± standard deviation) of 6866 ± 1677 reads per sample were obtained and samples were rarefied to 1800 reads for statistical comparisons, resulting in a mean estimated Good’s coverage of 97.1 ± 0.7%. After rarefaction, a mean of 162.6 ± 28.0 OTUs were observed among all samples.
Characteristics of samples with and without bacterial colonization were compared using Chi-squared analysis for categorical variables and Mann-Whitney-U test for continuous variables. Differences in overall community composition (beta diversity) were done using Bray-Curtis dissimilarities  and evaluated using analysis of similarity (ANOSIM)  with Bonferroni correction when appropriate using mothur. Samples were visualized by ordination using principal coordinate analysis (PCoA). Mean relative abundances of predominant families (found at mean abundances > 2.0%) were compared between groups with non-parametric Kruskal-Wallis rank sum test, followed by Bonferroni-corrected Dunn’s post-hoc test for pairwise comparisons. Analyses were conducted with XLSTAT (version 2020.2.3; Addinsoft, Belmont, MA, USA). All statistics were evaluated at α = 0.05.
Availability of data and materials
Raw data are stored in the Sequence Read Archive under BioProject accession number SRP197553.
Siegel RL, Miller KD, Fuchs HE, Jemal A. Cancer statistics, 2021. CA Cancer J Clin. 2021;71(1):7–33.
Huang L, Jansen L, Balavarca Y, Molina-Montes E, Babaei M, Van Der Geest L, et al. Resection of pancreatic cancer in Europe and USA: an international large-scale study highlighting large variations. Gut. 2019;68(1):130–9.
Pushalkar S, Hundeyin M, Daley D, Zambirinis CP, Kurz E, Mishra A, et al. The pancreatic cancer microbiome promotes oncogenesis by induction of innate and adaptive immune suppression. Cancer Discov. 2018;8(4):403–16.
Guo W, Zhang Y, Guo S, Mei Z, Liao H, Dong H, et al. Tumor microbiome contributes to an aggressive phenotype in the basal-like subtype of pancreatic cancer. Commun Biol. 2021;4(1):1019.
Riquelme E, Zhang Y, Zhang L, Montiel M, Zoltan M, Dong W, et al. Tumor microbiome diversity and composition influence pancreatic Cancer outcomes. Cell. 2019;178(4):795–806.
Geller LT, Barzily-Rokni M, Danino T, Jonas OH, Shental N, Nejman D, et al. Potential role of intratumor bacteria in mediating tumor resistance to the chemotherapeutic drug gemcitabine. Science. 2017;357(6356):1156–60.
Wei MY, Shi S, Liang C, Meng QC, Hua J, Zhang YY, et al. The microbiota and microbiome in pancreatic cancer: more influential than expected. Mol Cancer. 2019;18(97):1–15.
Fan X, Alekseyenko AV, Wu J, Peters BA, Jacobs EJ, Gapstur SM, et al. Human oral microbiome and prospective risk for pancreatic cancer: a population-based nested case-control study. Gut. 2018;67(1):120–7.
Castillo E, Meier R, Chung M, Koestler DC, Chen T, Paster BJ, et al. The microbiomes of pancreatic and duodenum tissue overlap and are highly subject Speci fi c but differ between pancreatic Cancer and noncancer subjects. Cancer Epidemiol Biomark Prev. 2019;28(2):370–84.
Kohi S, Macgregor-das A, Dbouk M, Yoshida T, Chuidian M, Abe T, et al. Alterations in the duodenal fluid microbiome of patients with pancreatic cancer. Clin Gastroenterol Hepatol. 2021. in press.
Jin Y, Gao H, Chen H, Wang J, Chen M, Li G, et al. Identification and impact of hepatitis B virus DNA and antigens in pancreatic cancer tissues and adjacent non-cancerous tissues. Cancer Lett. 2013;335(2):447–54.
Thomas RM, Jobin C. Microbiota in pancreatic health and disease: the next frontier in microbiome research. Nat Rev Gastroenterol Hepatol. 2020;17(1):53–64.
Goel N, Nadler A, Reddy S, Hoffman JP, Pitt HA. Biliary microbiome in pancreatic cancer: alterations with neoadjuvant therapy. Hpb. 2019;21(12):1753–60.
Meijer LL, Zwart ES, Brandt BW, Mebius R, Deng D, Giovannetti E, et al. Tumor microbiome: pancreatic cancer and duodenal fluids contain multitudes, … but do they contradict themselves? Crit Rev Oncol Hematol. 2019;144:2016–8.
Scheufele F, Aichinger L, Jäger C, Demir IE, Schorn S, Sargut M, et al. Effect of preoperative biliary drainage on bacterial flora in bile of patients with periampullary cancer. Br J Surg. 2017;104(2):e182–8.
Serra N, Di CP, Gulotta G, Arpa F, Giammanco A, Colomba C, et al. Bactibilia in women affected with diseases of the biliary tract and pancreas. A STROBE guidelines-adherent cross-sectional study in southern Italy. J Med Microbiol. 2018;67:1090–5.
Zeng S, Pöttler M, Lan B, Grützmann R, Pilarsky C, Yang H. Chemoresistance in pancreatic cancer. Int J Mol Sci. 2019;20(18):1–19.
Geller LT, Straussman R. Intratumoral bacteria may elicit chemoresistance by metabolizing anticancer agents. Mol Cell Oncol. 2018;5(1):3–5.
Aykut B, Pushalkar S, Chen R, Li Q, Abengozar R, Kim JI, et al. The fungal mycobiome promotes pancreatic oncogenesis via activation of MBL. Nature. 2019;574(7777):264–7.
Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Huntley J, Fierer N, et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012 Aug;6(8):1621–4.
Gohl DM, Vangay P, Garbe J, MacLean A, Hauge A, Becker A, et al. Systematic improvement of amplicon marker gene methods for increased accuracy in microbiome studies. Nat Biotechnol. 2016;34(9):942–9.
Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75(23):7537–41.
Staley C, Kaiser T, Vaugh BP, Graiziger CT, Hamilton MJ, Khoruts A, et al. Predicting recurrence of Clostridium difficile infection following encapsulated fecal microbiota transplantation. Microbiome. 2018;6:166.
Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig WG, Peplies J, et al. SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 2007;35(21):7188–96.
Cole JR, Wang Q, Cardenas E, Fish J, Chai B, Farris RJ, et al. The ribosomal database project: improved alignments and new tools for rRNA analysis. Nucleic Acids Res. 2009;37:D141–5.
Bray JR, Curtis JT. An ordination of the upland forest communities of southern Wisconsin. Ecol Monogr. 1957;27:325–49.
Clarke KR. Non-parametric multivariate analyses of changes in community structure. Aust J Ecol. 1993;18(1):117–43.
We would like to thank the University of Minnesota Biological Materials Procurement Network (BioNet) for providing archived samples. In addition, we thank the University of Minnesota Genomics Center for assistance with amplification and Illumina sequencing. Computational processing and analyses were carried out using the resources of the Minnesota Supercomputing Institute.
This study was funded by SPORE NCI grant (5P50CA101955–09), which funded archiving of the tissue samples. HN received partial support from the Hubbard Broadcasting Foundation. Neither funding body contributed to the study design, analysis, or interpretation of data.
Ethics approval and consent to participate
Patient samples were collected at the time of procedure under an Institutional Review Board (IRB)-approved protocol and provided informed, written inform consent. Retrieval of de-identified tissues was exempted from human research concerns (IRB STUDY00002967). IRB-approval was obtained to identify patient samples for further investigation of clinical information without further contact with the patient (IRB STUDY00007496).
Consent for publication
The authors declare they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Nalluri, H., Jensen, E. & Staley, C. Role of biliary stent and neoadjuvant chemotherapy in the pancreatic tumor microbiome. BMC Microbiol 21, 280 (2021). https://doi.org/10.1186/s12866-021-02339-3
- Pancreatic tumor microbiota
- Pancreatic adenocarcinoma
- Biliary stent
- Neoadjuvant chemotherapy