Pancreatic cyst fluid harbors a unique microbiome

Background It is clear that specific intestinal bacteria are involved in the development of different premalignant conditions along the gastrointestinal tract. An analysis of the microbial constituents in the context of pancreatic cystic lesions has, however, as yet not been performed. This consideration prompted us to explore whether endoscopically obtained pancreatic cyst fluids (PCF) contain bacterial DNA and to determine the genera of bacteria present in such material. Methods Total DNA was isolated from 69 PCF samples. Bacterial 16S rRNA gene-specific PCR was performed followed by Sanger sequencing and de novo deep sequencing for the V3-V4 variable region of 16S rRNA gene. Results We observed that 98.2% of the samples were positive in conventional PCR, and that 100% of selected PCF samples (n = 33) were positive for bacterial microbiota as determined by next generation sequencing (NGS). Comprehensive NGS data analysis of PCF showed the presence of 408 genera of bacteria, of which 17 bacterial genera were uniquely abundant to PCF, when compared to the Human Microbiome Project (HMP) database and 15 bacterial microbiota were uniquely abundant in HMP only. Bacteroides spp., Escherichia/Shigella spp., and Acidaminococcus spp. which were predominant in PCF, while also a substantial Staphylococcus spp. and Fusobacterium spp. component was detected. Conclusion These results reveal and characterize an apparently specific bacterial ecosystem in pancreatic cyst fluid samples and may reflect the local microbiota in the pancreas. Some taxa with potential deleterious functions are present in the bacterial abundance profiles, suggesting that the unique microbiome in this specific niche may contribute to neoplastic processes in the pancreas. Further studies are needed to explore the intricate relationship between pathophysiological status in the host pancreas and its microbiota. Electronic supplementary material The online version of this article (10.1186/s40168-017-0363-6) contains supplementary material, which is available to authorized users.


Background
Pancreatic cysts are fluid-filled neoplasms that can be detected with a frequency of up to 2% in the general asymptomatic adult population and have a low, but not negligible risk for evolving into pancreatic ductal adenocarcinoma, while accounting for up to 5% of the total incidence of pancreatic cancerous lesions [1,2]. The vast majority of cysts are coincidentally found during cross-sectional imaging done for other reasons than cyst-related symptoms. Optimal clinical management of pancreatic cysts remains controversial, but there is consensus in the field that increased insight into the molecular pathogenesis of pancreatic cysts may guide development of rational strategies in this respect. Unfortunately, the etiology of pancreatic cysts remains largely obscure.
Progress with respect to understanding the nature and natural history of pancreatic cysts is compounded by the presence of different types of this lesion. Grutzmann et al. and Farrell et al. have attempted to classify the different types of pancreatic cyst lesions and distinguish among others as intra-ductal papillary mucinous neoplasms (IPMNs), mucinous cystic neoplasms (MCN), serous cystadenomas (SCA), and pseudocysts [3][4][5]. In general, it is assumed that IPMN and MCN pose a higher risk of developing into carcinoma, with IPMN being more prevalent compared with MCN [1,6]. IPMNs are further classified as main branch, side branch, or mixed types, based on the involvement of the duct in the pancreas [5]. Presently, there are no validated biomarkers to identify cystic lesions that require surgical resection and this constitutes a major challenge in this field. Although pancreatic lesions develop into malignancy in only up to 3 % of cases, 10 % of patients with such lesions undergo resection [5] suggesting the need for superior clinical tests and patients' stratification prior surgery. Currently, the decision for resection of pancreatic cyst lesions and/or continued monitoring is made according to the Sendai guidelines after evaluation of different clinical tests [7][8][9]. The available clinical tests include different biochemical analyses, cytology, pathological assessment of fine needle biopsy or aspiration material, endoscopic ultra-sonography (EUS), and radiological diagnosis such as endoscopic retrograde cholangiopancreatography (ERCP), magnetic resonance cholangiopancreatography (MRCP), and whole-body computerized tomography (CT). The inter-observer agreement, however, between different modalities remains moderate [10]. Therefore, a set of preoperative biochemical analyses have been increasingly used in clinical decision-making. This includes the study of cyst fluids and serum for the characteristic presence of carcinoembryonic antigen (CEA), cancer antigen 19.9 (CA19.9), cancer antigen 72.4 (CA72.4), cancer antigen 15.3 (CA15.3), pancreatic amylase, and mucin antigens, along with other cyst characteristics [11,12]. Other tests are based on specific analysis of different genetic modalities like K-RAS mutation and integrity, but efforts attempting to provide clinical validation for such tests have largely proven unsuccessful [1,6,13,14]. Increased insight into the factors that facilitate the development of cystic lesions would evidently benefit the identification of tests capable of providing guidance for clinical management of asymptomatic patients exhibiting pancreatic cysts.
Intriguingly, the human gut microbiome has emerged recently as an important environmental factor linked to the development of different intestinal and extra-intestinal malignancies [15][16][17]. In the stomach, Helicobacter pylori remains the archetypical example of a prokaryotic organism that can initiate a cascade of molecular events finally leading to fullblown cancer, whereas in the colon, various organisms and especially Fusobacteria have been linked to the appearance of dysplasia (reviewed in [17]). Whether the appearance of pancreatic cysts is linked to the presence of bacteria per se and if so whether specific types of bacteria are associated to the presence of cystic pancreatic neoplasms has remained unexplored. If, however, the presence of pancreatic cyst can be linked to the microbiome, this would entail a significant step forward with respect to our understanding of pancreatic cystogenesis.
The above-mentioned considerations prompted us to explore the potential microbial component of pancreatic cyst fluid. The results show that presence of bacterial DNA is common to such material and that especially Fusobacterium spp. and Bacteroides spp. are prominently present in such material. As some of these bacterial species have been linked to dysplastic processes elsewhere in the tracts, a causal link between the presence of such bacteria and pancreatic neoplasm may also exist and the results may indicate that bacterial colonization of pancreatic cyst fluid is a regular phenomenon.

Patient samples and pancreatic cyst fluid collection
A cohort of 69 patients with suspected cystic lesions was established between the period of 2008 and 2013 (Table 1 and Additional file 1).The pancreatic cyst fluids (PCF) were collected after a signed informed consent from these patients, who were undergoing endoscopic ultrasound fine needle aspiration (EUS-FNA) at the Department of Gastroenterology, Erasmus MC, The Netherlands. Following collection, pancreatic cyst fluids were transferred to the laboratory and stored at − 150°C until analysis.

DNA isolation
Approximately 300 μl from 69 pancreatic cyst fluid samples were used for total DNA isolation. After bead beating (Fast Prep®-24 Instrument), the supernatant and pellet were separated by centrifugation at 13000 rpm for 1 min and both pellet and supernatant were used for total DNA isolation, using the Wizard DNA isolation kit as specified by the manufacturer's protocol (catalog no. A1620, Promega BNL B.V, The Netherlands). Isolated DNA was equilibrated in the DNA rehydration solution from the kit and quantified on nanodrop-2000 spectrophotometer (Isogen Life Science BV, De Meern, The Netherlands). Pancreatic cyst fluid DNA was diluted to 1 ng/μl for the PCR analyses and subsequently stored at − 20°C.

16S rRNA gene sequencing analyses Sanger sequencing
Sanger sequencing was done in order to identify the bacterial 16S rRNA genes present in the PCF. The PCR products generated using universal primers of 16S rRNA were sequenced using F primer 5′-CTTAGGAATGAGACAGAGATG-3′ through LGC Genomics GMBH, Germany. The chromatograms were analyzed, sequences were curated and identified using web databases of integrated microbial genomes and 16S rRNA-specific nucleotide BLAST. The sequences have been deposited at NCBI under gene accession identifiers MF061964-MF061990.

Deep sequencing data analysis
Two types of deep sequencing analyses were performed, one with paired-end reads alone discarding the single pair reads leading to 408 genera identification, and one with paired-end reads and single pair reads leading to identification of 785 genera [18]. The data generated via paired-end reads alone was further used for downstream analysis for comparison of the PCL bacteria population vs human microbiome project samples of 13 different body sites and comparisons between the types of clinically defined cysts types and resection (Additional file 1).

Pre-processing of amplicon reads and sequence data QC
Illumina Miseq sequence data was quality checked using FastQC and Rqc software's [19]. Base call quality distribution, percentage of bases with quality scores above Q20 and Q30, GC percentage, and sequencing adapter contamination were used for assessing the data quality. The details of the number of reads and base quality obtained for each sample are provided in Additional file 2.

Taxonomic classification and OTU clustering
For downstream analysis of the metagenome data, "mothur" software bundle was used [20].The qualityfiltered sequence reads were imported into mothur, and the read pairs were aligned with each other to form contigs based on sequence consensus. This results in longer contigs that span the targeted V3-V4 hypervariable region. These contigs were screened and only those between 470 and 500 bp were taken for further analysis. Contigs with ambiguous base calls were rejected; high-quality contigs were checked for identical sequences and duplicates merged. The filtered contigs were aligned to a known database for 16S rRNA [Silva V.119] [21]. These classified contigs were filtered for any undesired lineage from the taxonomy file. Any ambiguous contigs aligning to untargeted regions [other than V3-V4] were discarded. Using UCHIME algorithm, chimeric contigs were flagged and removed; a known reference for chimeric sequence was cross referenced [22]. These final set of contigs were then phylotype binned into operational taxonomic units (OTU) based on the Silva V.119 database and the abundance of each OTU in the population was estimated.

Comparison of the taxonomic classification of the pancreatic microbiomes with those from other body sites
In order to identify microbial signature of pancreatic cysts, a comparative analysis was carried out against organ-specific bacterial profiles obtained from the NIH Human Microbiome Project (HMP) database [23]. The pancreatic microbial profile of the current study was compared with microbial profiles of 13 different body sites downloaded from HMP (SRA numbers provided in Additional file 3). The notion that bacteria in the gut  Table 1) were exposed to bead beating and centrifugation followed by DNA extraction from both the pellet as well as the supernatant. Subsequently, the DNA was probed for the presence of sequence coding for bacterial 16S RNA through PCR. Reaction products were resolved by gel electrophoresis employing an agarose gel. The results show that bacterial DNA is commonly found in pancreatic cyst fluids may affect diseases outside the gut and vice versa is gradually accepted, and we therefore considered both anatomical sited in close proximity to the pancreas as well as more distant sites.
The comparative analysis was carried out between anterior nares, antecubital fossa, buccal mucosa, gingiva, hard palate, mid vagina, palatine tonsils, posterior fornix, retroauricular crease, saliva, stool, throat, and tongue. Although phylotype-based analysis is limited as compared to distance-based OTU classification, it allows for investigating the relationship of the PCF microbiota to previously characterized microbes in the HMP database.
Two analyses were performed for identifying the unique microbiota when PCF commensals were compared to the bacterial communities associated with 13 different body sites. First, a pairwise binomial test was carried out against each organ in order to identify significantly (FDR < 0.05, FC > 3) and absolute abundance difference of 10 abundant bacterial species in pancreas. Second, in order to identify the relatedness of the PCF microbiomes with those from other body sites, we compared bacterial communities of the PCF samples to those of taxonomic profiles from 13 other body sites from the HMP database by using principal component analysis (PCA). Third, statistical comparisons of these microbiomes were performed using the STAMP analysis package. The statistical tests used was ANOVA and Welch's t test, with Benjamini-Hochberg multiple test corrections. Those taxa having a higher abundance in the pancreatic cysts (with corrected P value < 0.01) were identified as those that are specifically abundant in the PCF.

Diversity analysis between the cyst types
Using STAMP statistical analysis package, those taxa having a higher abundance between the test groups IPMN vs pseudocyst vs MCN vs SCA were identified as specific for the groups to be classified based on the bacterial population. The cysts were grouped based on the pathological classifications, and they are tested using White's non parametric t test, two tailed with Benjamini-Hochberg multiple test corrections of P value (P < 0.05) [24]. The measured levels of CEA and CA19.9 ranges were used for classification of cysts as cysts, benign cysts, and malignant cysts according to study conducted by Talar-Wojnarowska et al. [25]. Then the classified data was used for the specific bacterial taxa identification between the resected types using the same parameters.

Statistics and calculations
All the statistical analysis were done using excel, Graphpad Prism 5.0., and STAMP statistical analysis software.
When appropriate, results were Bonferroni corrected in graphpad analysis and Benjamini-Hochberg multiple test corrections in STAMP. P values < 0.05 (between cysts and resection types) and < 0.01 [between body sites (NGS)] were considered significant.

Study sample characteristics
For analysis of the potential microbial component of pancreatic cyst fluid (PCF), material collected from 69 patients was used. The characteristics of these patients are listed in Table 1 and Additional file 1. Of the PCF analyzed, 27 were obtained from patients harboring an IPMN (39.1%), including two patients with main branch IPMN (2.9%), four patients with mixed type IPMN (5.8%), 12 patients with a unifocal side-branch IPMN (17.4%), and four patients presenting with a multifocal side-branch IPMN (5.8%). In five patients, the type of IPMN was unclassified (7.2%). In addition, our cohort contained 13 patients with an MCN (18.8%), nine patients with pseudocysts (13.0%), and nine patients with serous cystadenomas (13.0%). Finally, the cohort contained 11 patients with apparently multiple forms of cystic lesions, gastrointestinal stromal tumor (GIST), neuroendocrine tumors (NET), or having no definitive clinical diagnosis and for this study, these patient were classified as "others" (15.9%). We concluded that this cohort would allow the study of potential microbiological constituents of PCF and to relate results to the clinical phenotype of the patient from which the fluid was obtained.
Bacterial DNA is commonly present in EUS-FNA-collected pancreatic cyst fluids As it is as yet unknown whether pancreatic cystic fluid hosts a microbiological component, analyzed the EUS-FNA-collected pancreatic cyst fluid obtained from our cohort for the presence of significant amounts of bacterial DNA. Importantly, we found that the majority of these fluid samples were rich for bacterial DNA, with 16S rRNA PCR demonstrating the presence of bacterial DNA in 64 (92.8%) out of the 69 samples (Table 2 and Fig. 1). The presence of bacterial DNA in cyst fluid did not statistically relate to the type of lesion from which it was obtained (P value > 0.99; χ 2 test): in mucinous cystic neoplasms, 100% of samples contained significant amounts of bacterial DNA; in IPMNs, this number was 92.6%; in pseudocysts, 88.9% of PCF were positive for bacterial DNA; in serous cystadenomas, 88.9%; whereas in the group of others which included GIST, NET, and clinically undefined samples, 90.9% displayed significant amounts of bacterial DNA (Additional file 1 and Table 2).

Microbial composition of pancreatic cyst fluid with a differential clinical aspect
The presence of a microbial component in PCF raises obvious questions as to the identity of the organisms apparently present in such fluid. Sanger sequencing of PCR products generated using universal 16S rRNA demonstrated that Fusobacterium spp. is present in 13 out of all 69 PCF samples (18.84%). Another predominating bacterium in PCF was Bacillus spp. which was present in 16 out of 69 (23.19%) samples. The presence of other bacteria was also noted which included Ruminococcus spp., Staphylococcus spp., Caldimonas spp., Arthrobacter spp., Acinetobacter spp., Bacteroides spp., Orpinomyces spp., and Anaerococcus spp. (Table 2).
To confirm the presence of these bacteria and gain more insight into the bacterial composition of the PCF, DNA obtained from the 23 fluids containing the highest apparent concentration of 16S rRNA copies were sent out for 16S rRNA de novo sequencing of V3-V4 region, using universal 16S rRNA primers which should allow identification of the bacteria present in such fluid at least on genus level. The results (shown in Fig. 2 and Additional file 2) are consistent with the presence of diverse bacterial ecosystems in such fluids, with as most predominant genera present Bacteroides, Escherichia/ Shigella, and Acidaminococcus, but in total, 408 different genera were detected in the 33 samples analyzed, of which 93 genera were found in at least 50% of PCF samples analyzed. Different types of cysts were not statistically different with respect to microbial composition: for none of the 408 genera, a Bonferroni-corrected statistically significant difference in abundance was detected when serous cystadenoma, pseudocysts, IPMN, or mucous cystic neoplasm-derived fluids were compared. Furthermore, when low P values (Bonferroni-uncorrected P value < 0.05 but > 0.0001) were considered, it appeared that such findings were limited to very low abundant organisms (< 0.05% of all bacteria) that were relative low in one of the groups compared, more indicative of technical detection problems rather than a reflection of true biological differences between the groups (Additional file 2). In apparent agreement with the notion that the ecological niche provided by cyst fluid is relatively similar between different cyst manifestations is also the observation that the Shannon index for ecological diversity (H′) (Fig. 3) is not different between the different groups (P value = 0.99; one-way ANOVA).

Diversity within the cyst type and resections
There were no significant difference in the diversity of bacterial microbiota seen in the 33 sample group of deep sequenced data, between the cyst types, or between the resected vs non resected groups or between the CA19.9 or CEA range-based classification groups. The results were in accordance with the bacterial presence in the  Fig. 4 and Additional file 4). ANOVA analysis, which is performed to identify the effect size of the bacterial population, reveals that 314 bacterial microbiota show an effect size variation between the bacterial commensals in PCF versus selected body sites, high effect size (0.1 to 0.9 and P < 0.01) bacterial populations are shown in Fig. 5, and total effect size variation is shown in Additional file 5. The PCA of different body sites and PCF are shown in Additional file 6: Figure S1. The Welch's t test was performed to identify the specific bacterial populations having high and low abundance in the PCF when compared to HMP selected body sites. A total of 136 bacteria (Additional file 6: Figure S2 and Additional file 7) are identified, of which 17 abundant bacteria in PCF were highly unique and potentially pro-cancerous. They are, in ascending order of abundance in PCF, Coprococcus spp., Collinsella spp., Butyricicoccus spp., Ruminococcus spp., Parabacteroides spp., Alistipes spp., Clostridium XI spp., Gemmiger spp., Dorea spp., Lachnospiracea incertae sedis spp., Blautia spp., Bifidobacterium spp., Sphingomonas spp., Acidaminococcus spp., Faecalibacterium spp., Escherichia/Shigella, and Bacteroides spp. (P < 0.0001). A role for many of these bacteria in initiation and progression of colon, lung, and liver cancer has been suggested [26]. The following 15 bacterial genera were higher in all the 13 body sites used in this analysis: Streptococcus spp., Propionibacterium spp., Lactobacillus spp., Fusobacterium spp., Corynebacterium spp., Veillonella spp., Neisseria spp., Staphylococcus spp., Porphyromonas spp., Prevotella spp., Leptotrichia spp., Actinomyces spp., Capnocytophaga spp., Gemella spp., and Selenomonas spp. (P < 0.0001). These were also the least abundant bacteria present in the PCF. This comparison suggests that the PCF fluid bacterial colonization is unique and characterized by high abundant genera which have a tendency to propagate in the cancerous microenvironment, feeding the tumor (Fig. 6).

Discussion
In the current study, we have demonstrated that 92.8% of a large and clinically representative collection of pancreatic cyst fluids (PCF) is host to a complicated ecosystem of bacteria. The distinct bacterial community had a rich species diversity which is very different from that observed elsewhere in the proximal human digestive tract. Although the EUS-FNA is not a sterile procedure, 7.8% of the PCF samples were negative for bacterial DNA arguing that the procedure did not have a major impact on our cultivation-independent bacteriological findings. Furthermore, the ecosystem observed by deep sequencing of pancreatic cyst material is inconsistent with that expected from contamination by the oral, nasal, pharyngeal, oesophageal, gastric, or duodenal flora. Members of the genus Streptococcus constitute by far the major component of the bacterial ecosystem in the esophagus [27], but while streptococcus was detected in all samples sequenced, it constitutes less than 1% of the total microbiota in any sample in our series, with the exception of one sample obtained from an IPMN. Likewise, the stomach is relatively sterile and dominated by H. pylori. The genus Helicobacter, however, was rare in our analysis, this genus the 67th genus in our ranking of our samples. Other bacteria found regularly in the stomach include Streptococcus, Neisseria, and Lactobacillus, and these species are not dominant in our analysis (being ranked as the 18th, 129th, and 54th most prevalent genus) [20]. Thus, the esophagus and stomach are an unlikely source of major contamination in our results. The duodenum can be home to substantial bacterial content, as also evident from conditions such as small intestinal bacterial overgrowth or SIBO. Valeria D' Argenio et al. recently reported on the bacterial composition of the duodenum in a cohort including 15 healthy individuals [28]. The most prevalent genus emerging from this analysis was Propionibacterium, which, although fairly often detected, was relatively rare in our analysis of pancreatic cyst fluids, with its abundance not even approaching 1% of the entire flora in any of the samples tested. Other genera dominating the duodenal microbiota include Porphyromonas (ranked 76th in cyst fluid), Streptococcus, Neisseria, and Heamophillus (ranked 21st, 58th, and 57th in cyst fluid). Hence, the bacteria we observed in cyst fluid widely diverges from that expected if duodenal contamination was a major factor. Indeed in an analysis of microbiome of the healthy proximal tractus in our own institution (manuscript in preparation), we observed that Streptococcus, Veillonella, Prevotella, and Pseudomonas have the highest relative abundance from distal esophagus to the jejunum and proximal ileum. This situation is markedly different from that observed in pancreatic cyst fluids (although Prevotella was the 7th most abundant genus in pancreatic cyst fluid, this is still markedly lower as that observed in the proximal gastrointestinal tract of volunteers). In apparent agreement, in a preliminary series of experiments employing fluorescent in situ hybridization on surgically obtained pancreatic cyst material, bacteria were apparent (not shown). Thus, the most straightforward interpretation of results is that the pancreatic cyst is home to a previously unsuspected and also unique bacterial ecosystem.
The presence of such a bacterial ecosystem in pancreatic cysts raises important questions as to the role of the bacteria present in such cysts in the development of such structures. In general bacteria have been linked to transformation and trans-differentiation of endodermal epithelia, with H. pylori being the most important example. The genus Helicobacter was only marginally detected in pancreatic cyst fluid; other bacterial species linked to transformation in the intestine were, however, more prevalent. In particular, F. nucleatum, which is relatively predominant in a fraction of the samples and has been convincingly linked to neoplasm formation in the colon, excites interest in this respect. Emerging evidence suggests that Fusobacterium species detected in pancreatic cancer tissues is associated with a worse clinical outcome in pancreatic cancer patients [29]. But also other bacteria seen in pancreatic cyst fluids, even if it is at low abundance, may have a role here. Using NGS technology, multiple studies have compared oral microbiota between healthy individuals and those with pancreatic cancer [4,30]. Although no correlation between known pro-oncogenic oral pathogens and pancreatic cancer was detected, Lin et al. showed that pancreatic cancer patients had significantly higher levels of Bacteroides genus compared to control subjects [30], which is in agreement with our finding that Bacteroides spp. were the most predominant genus in PCF. In addition, a recent pilot study using 16S rRNA gene sequencing on saliva specimens showed that while there was no difference in diversity of oral microbiota between patients with pancreatic ductal adenocarcinoma (PDAC), IMPN, or healthy controls, pancreatic cancer patients had higher levels of members of the phylum Firmicutes while healthy individuals were associated with higher relative abundances of Proteobacteria [31]. Evidence gathered from those studies supposes that pancreatic malignancy may be associated with changes in abundances of some groups of bacteria in the human digestive tract. Therefore, dysbiotic microbiota in the upper digestive tract including oral cavity may interrupt the unique ecosystem in pancreatic cyst fluids along the neoplastic process in pancreas. At bay with the notion of a causative role for the bacterial flora in pancreatic cyst formation, however, is the observation that both detection of bacteria per se or the composition of the cyst flora did not show correlation to clinical parameters and thus the role of the flora in cystogenesis, if any, awaits further study.
Pancreatic infections mainly arise from translocation of bacteria from the small bowel, and rarely from the colon and oropharyngeal route as demonstrated by study Fig. 4 Pancreatic cyst microbiome is unique among other human body microbiomes. Publicly accessible databases were mined for composition of microbiomes at different human body sites and compared to those observed in PCF. It appears that there are 27 to 314 bacterial genotypes differently present in the PCF when compared to the selected body sites when analyzed via pairwise binomial test with high abundance PCF bacteria(P < 0.0001) and ANOVA test, respectively on Veillonella and Bifidobacterium spp. which were identified in pancreatic abscesses [32]. A study of Brook et al. has identified 158 bacterial species from pancreatic abscesses, of which 77 isolates were aerobic and the remaining 81 were anaerobic bacteria [33]. The most commonly detected microorganisms in infected pancreatic pseudocysts include often not only opportunistic bacteria like E. coli, Enterobacter spp., Klebsiella spp., and Staphylococcus spp. but also fungal isolates including Candida albicans (15 case studies) [34]. Importantly, EUS FNAB procedure caused serious Clostridium perfringens infections in five patients leading to pancreatitis and pancreatic cyst formation, which required surgical interventions [35]. Such studies illustrate the nature of the bacterial transfer from the early to mid-gut commensal bacteria to the pancreas. Yet none of these earlier studies have directly proved the presence of bacteria in the pancreatic cyst and its fluid. Its apparent divergence from the flora in the duodenum is probably a reflection of the specific conditions in pancreatic cysts which include an absence of contact with the digestive nutrient stream, an exposure to high levels of pancreatic secretes and the presence of abundant mucus. At present, clinical behavior of pancreatic cysts is very difficult to predict and adequate management of cysts represents one of the largest challenges in clinical gastroenterology. It is tempting to speculate, however, that changing conditions in cysts related to transformation of the structure to full-blown malignant cancer would also influence the milieu interieur formed by cyst fluid and hence have a profound effect on the bacterial composition. As such, changes in the bacterial communities may serve to detect such transformation in cystforming structures and may become useful for guiding clinical management.
Irrespective, however, of its potential as future diagnostic and prognostic marker, the present study shows an as yet unknown bacterial ecosystem in pancreatic cyst fluid. As it is evident that bacteria influence physiology and pathophysiology of their interacting epithelia Fig. 5 a, b Contribution of individual bacterial genera to the unique aspect of PCF fluid. For 82 different genera the relative in abundance in PCF was compared to that 13 other body sites and the relative contribution to the PCF-specific nature of the microbiome spectrum was calculated. In additional files information on a further 232 genera can be found. The results show that PCF contains a microbiome that is characterized by an overall uniqueness that cannot be attributed to a single genus everywhere in the gut, it is well possible that such interactions also exist in the pancreas and its cysts, and that the biology of the pancreatic cysts is in important ways shaped by this ecosystem. Studies addressing this possibility are currently in progress.

Conclusions
The study reveals previously undescribed bacterial diversity present in human pancreas and its cyst fluids. As specific bacteria are associated with this body site, we propose that such bacteria may carry the potential to influence the development of pathophysiological processes in the pancreas. The study points out to the need to further explore the microbiome in this specific niche for diagnostic and therapeutic purposes.  Figure S2. Difference of mean of selected high and low abundance bacterial microbiome in PCF and 13 body sites together, respectively. Difference of mean between the Fig. 6 Bacterial commensals of PCF and selected different body sites comparison (HMP database) shows the difference in the distribution uniqueness in the PCF than the body sites. One hundred thirty-six bacterial genus of high and low abundance with P < 0.0001, and their difference of mean is plotted with 99% confidence interval, obtained via the Welch's t test comparison between the groups bacterial genus distribution of pancreatic cyst fluid and 13 different body site microbiome selected from Human Microbiome Project database. When compared 17 bacterial genus (with p < 0.01) showing high abundance in PCF and 15 bacterial genus (with p < 0.01) showing high abundance in 13 body sites selected. (PPTX 7223 kb) Additional file 7: Welch's t test statistics run for the bacterial genotypes comparison calculation of the PCF and 13 different body sites microbiota. (XLSX 9711 kb)