<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Oncol.</journal-id>
<journal-title>Frontiers in Oncology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Oncol.</abbrev-journal-title>
<issn pub-type="epub">2234-943X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fonc.2024.1374769</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Oncology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Detection and comparison of tumor cell-associated microbiota from different compartments of colorectal cancer</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" equal-contrib="yes">
<name>
<surname>Zuo</surname>
<given-names>Yanzhen</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="author-notes" rid="fn003">
<sup>&#x2020;</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author" equal-contrib="yes">
<name>
<surname>Lu</surname>
<given-names>Yanjie</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="author-notes" rid="fn003">
<sup>&#x2020;</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1312386"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Pang</surname>
<given-names>Jiayu</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Jin</surname>
<given-names>Shunkang</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhang</surname>
<given-names>Xinyu</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Zhao</surname>
<given-names>Enhong</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Li</surname>
<given-names>Yuhong</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Cancer Research Laboratory, Chengde Medical College</institution>, <addr-line>Chengde, Hebei</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Department of Gastrointestinal Surgery, Affiliated Hospital of Chengde Medical College</institution>, <addr-line>Chengde, Hebei</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Ales Lapanje, Institut Jo&#x17e;ef Stefan (IJS), Slovenia</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Ye Yang, University of Florida, United States</p>
<p>Zhengrui Li, Shanghai Jiao Tong University, China</p>
<p>Yefei Zhu, Tongji University, China</p>
<p>Valeriy Poroyko, Laboratory Corporation of America Holdings (LabCorp), United States</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Yuhong Li, <email xlink:href="mailto:youngcheer2003@foxmail.com">youngcheer2003@foxmail.com</email>; Enhong Zhao, <email xlink:href="mailto:Zheng13253189639@163.com">Zheng13253189639@163.com</email>
</p>
</fn>
<fn fn-type="equal" id="fn003">
<p>&#x2020;These authors have contributed equally to this work and share first authorship</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>21</day>
<month>05</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>14</volume>
<elocation-id>1374769</elocation-id>
<history>
<date date-type="received">
<day>22</day>
<month>01</month>
<year>2024</year>
</date>
<date date-type="accepted">
<day>30</day>
<month>04</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 Zuo, Lu, Pang, Jin, Zhang, Zhao and Li</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Zuo, Lu, Pang, Jin, Zhang, Zhao and Li</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<sec>
<title>Introduction</title>
<p>Intratumoral microbes play an important role in the development of colorectal cancer (CRC). However, studying intratumoral microbes in CRC faces technical challenges, as tumor microbe communities are often contaminated by fecal microbes due to the structure of the gut folds and villi. The present study aimed to develop a new method for isolating tumor cell-associated microbiota and comparing microbial populations from different compartments.</p>
</sec>
<sec>
<title>Materials and methods</title>
<p>The distribution of intestinal bacteria was detected using immunohistochemistry combined with 5R-16s rRNA gene sequencing to explore the effects of the sampling site and number of washes on the detection of microbiota. The 5R-16s rRNA gene sequencing was performed using 44 samples from 11 patients with CRC, including CRC tumor tissues (TT), normal tissues adjacent to CRC (NT), tumor cells (TC), and normal cells (NC). TC and NC were obtained from the TT and NT using an enzymatic digestion method. The microbiota and their potential functions in the four groups were analyzed and compared to determine the differential microbiota related to CRC.</p>
</sec>
<sec>
<title>Results</title>
<p>Bacteria were mainly distributed in the feces covering intestinal tissues and in the epithelial cells and macrophages within the tissues. Different sampling sites and number of washes led to detection of different microbiota distributions. Although the cleaning method could be controlled, sampling sites varied and led to different microbiota distributions. The phyla of Firmicutes and Bacteroidetes were highly abundant in the conventionally used tissue samples, whereas Proteobacteria was the most abundant phyla in the cell samples isolated with the new method (i.e., after cell enzymatic hydrolysis). Detection of CRC cell-associated microbiota using a cell enzymatic digestion method showed that some bacteria, such as <italic>Fusobacterium, Eikenella, Shewanella</italic>, and <italic>Listeria</italic>, were more abundant in TT than NT, whereas the abundance of <italic>Akkermansia</italic> was lower in TT than NT. The tumor/normal ratios of some bacteria, such as <italic>Gemella, Escherichia, Shigella</italic>, and <italic>Blautia</italic>, were different between the cell and tissue samples.</p>
</sec>
<sec>
<title>Conclusion</title>
<p>The cell enzymatic digestion method reduced fecal bacterial contamination, enabling low biomass intratumoral microbiota to be detected and allowing prediction of bacterial distributions.</p>
</sec>
</abstract>
<kwd-group>
<kwd>colorectal cancer</kwd>
<kwd>tumor cell-associated microbiota</kwd>
<kwd>microbial populations</kwd>
<kwd>contamination</kwd>
<kwd>enzymatic digestion</kwd>
</kwd-group>
<counts>
<fig-count count="7"/>
<table-count count="2"/>
<equation-count count="0"/>
<ref-count count="29"/>
<page-count count="14"/>
<word-count count="6073"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-in-acceptance</meta-name>
<meta-value>Gastrointestinal Cancers: Colorectal Cancer</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<label>1</label>
<title>Introduction</title>
<p>Colorectal cancer (CRC) is one of the most common malignancies of the digestive tract. According to the 2020 Global Cancer Statistics report, CRC has the third highest incidence and the second highest mortality in the world (<xref ref-type="bibr" rid="B1">1</xref>). The largest number of microbes in the human body exist in the gastrointestinal tract, which comprise a complex ecosystem known as the gut microbiome. With the development of new sequencing technology and bioinformatics, it is possible to further study the composition, metabolism, and metagenome of intestinal microecology (<xref ref-type="bibr" rid="B2">2</xref>). While many studies have revealed a strong relationship between gut microbiota and CRC, most have focused on fecal microbiota and their metabolites (<xref ref-type="bibr" rid="B3">3</xref>). Due to the complexity of fecal microbiota and the multi-factor regulation of their interactions with the intestinal mucosal epithelium, as well as the fact that most of the bacteria are &#x201c;passing bacteria&#x201d; rather than &#x201c;colonizing bacteria,&#x201d; exploring the direct evidence and biological behavior of bacterial enrichment and their invasion of tumor tissues and cells has become a new research direction in this field. In recent years, some studies have focused on tumor-associated microbiota. Tumor-associated microbiota are intrinsic components of the tumor microenvironment in human cancer types (<xref ref-type="bibr" rid="B4">4</xref>&#x2013;<xref ref-type="bibr" rid="B6">6</xref>). Intratumoral microbiota are integral tumor components that play critical roles in shaping the tumor microenvironment (<xref ref-type="bibr" rid="B7">7</xref>, <xref ref-type="bibr" rid="B8">8</xref>). Intratumoral bacteria are mostly intracellular and are present in both cancer and immune cells; each tumor type has a distinct microbiome composition (<xref ref-type="bibr" rid="B5">5</xref>). The functions of intratumoral bacteria are associated with the tumor type/subtype and their responses to immunotherapy (<xref ref-type="bibr" rid="B5">5</xref>). Cancer cells infected with bacteria invade the surrounding environment as single cells and recruit myeloid cells to the bacterial regions. Intratumoral microbiota are mainly distributed and enriched in microniches, influencing immune and epithelial cell functions that promote cancer progression (<xref ref-type="bibr" rid="B4">4</xref>). Tumor-associated bacteria play a crucial and direct role in inducing tumorigenesis, changing cell metabolism, promoting tumor metastasis, and preventing immune surveillance (<xref ref-type="bibr" rid="B9">9</xref>). Investigating tumor-associated microbiota may aid the discovery of novel treatment options for patients with cancer.</p>
<p>However, unlike other solid tumors, such as breast cancer, which have no direct connection with the external environment, CRC tumors are inherently contaminated by fecal microbes owing to the structure of the gut folds and villi (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1A</bold>
</xref>). Preventing contamination with and detection of tumor-associated microbiota are challenging due to technical difficulties. In the present study, we explored different methods for detecting the distribution of gut microbes and established a new method for isolating the tumor cell-associated microbiota of CRC. The microbiota and their potential functions were analyzed and compared using different methods, and the tumor cell-associated microbiota of CRC was obtained.</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>
<bold>(A)</bold> Schematic diagram of the distribution of colorectal and intestinal microbes. CRC tumors are inherently contaminated by fecal microbes owing to the structure of the gut folds and villi. &#x2460;, &#x2461;, &#x2462;, and &#x2463; are schematic diagrams of the sampling sites. CRC tumor tissues and the adjacent normal tissues from 11 patients were taken as shown in &#x2460; through &#x2463;. <bold>(B)</bold> a schematic diagram of isolation process of cell-associated microbiota. <bold>(C)</bold> Immunohistochemical analysis of LPS. LPS staining was mainly observed in the feces from the surface of intestinal tissues and in epithelial cells and macrophages within the tissues. LPS, lipopolysaccharide.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fonc-14-1374769-g001.tif"/>
</fig>
</sec>
<sec id="s2" sec-type="materials|methods">
<label>2</label>
<title>Materials and methods</title>
<p>This study was performed in accordance with the guidelines of the 1975 Declaration of Helsinki, revised in 1993. All participants provided informed consent and the study was approved by the Medical Ethics Committee of Chengde Medical University (No. 202215).</p>
<sec id="s2_1">
<label>2.1</label>
<title>Patients</title>
<p>CRC tumor samples were obtained from patients who underwent surgery at the Affiliated Hospital of Chengde Medical College, Chengde City, Hebei Province, China, between September and December of 2022. The inclusion criteria were: (1) patients who had been diagnosed with CRC; (2) no previous history of colorectal surgery; (3) no use of antibiotics or corticosteroids or probiotics one month prior to the surgery; (4) no familial adenomatous polyposis or hereditary nonpolyposis CRC; no inflammatory bowel disease, metabolic diseases (such as diabetes, obesity, hyperlipidemia), infectious diseases, severe liver disease, kidney diseases, or immunodeficiency; (5) no special dietary habits; and (6) patients who volunteered to participate in the study. The exclusion criteria were as follows: (1) pregnancy or lactation, (2) previous cancer diagnosis, (3) patients who had received neoadjuvant chemoradiotherapy, and (4) unwillingness to participate in the study. Based on the aforementioned inclusion and exclusion criteria, 12 patients with CRC were included in this study.</p>
</sec>
<sec id="s2_2">
<label>2.2</label>
<title>Colorectal cancer and control samples</title>
<p>Samples were collected from the 12 patients through surgery. Following radical resection of CRC, the tumor tissue (TT) and adjacent normal intestinal tissue (NT, more than 5 cm distant from the cancer) were separately cut with a sterile knife. The samples were placed in sterile tubes, immediately transported to the laboratory in an icebox, and washed three times with sterile saline. A portion of each sample was taken and fixed in 4% paraformaldehyde separately for histological and immunohistochemical (IHC) staining, and the remaining tissue was used in the following steps. TT and NT from 11 patients were taken as shown in <xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1A</bold>
</xref> &#x2460; through &#x2463;. They all contained areas of contact with feces and duplicate samples were taken, one being frozen in liquid nitrogen for subsequent 5R-16s rRNA gene sequencing and another being used to isolate tumor cells (TC) or normal cells (NC) as described below. The sampling of the twelfth patient was slightly different from the 11 patients mentioned above. As shown in <xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1A</bold>
</xref>, TT from the twelfth patient was sampled from three areas (&#x2460;, &#x2461;, and &#x2462;), part &#x2460; containing fecal bacteria adhering to the surface of the tumor tissue while part &#x2461; and &#x2462; containing no areas in contact with fecal bacteria. During the surgery, a pathologist identified the mucosal part and used a surgical knife to remove it. Part &#x2462; was washed three more times with sterile saline. Three NT samples were taken as shown in part &#x2463;, one being frozen in liquid nitrogen for subsequent 5R-16s rRNA gene sequencing, the second being washed three more times, while the third was washed six times with sterile saline. Six samples from the twelfth patient were frozen in liquid nitrogen for subsequent 5R-16s rRNA gene sequencing. All tissues were sampled in a sterile biosafety cabinet. Each sample had a wet weight of at least 200 mg. An empty tube was set up as an environmental control tube in parallel to the experimental tubes during sample collection procedures, with a total of 3 environmental control tubes used.</p>
</sec>
<sec id="s2_3">
<label>2.3</label>
<title>Histological and IHC staining</title>
<p>The paraffin-embedded tissues were cut into 4-&#x3bc;m thick sections. After deparaffinization and hydration, a representative section of each sample was stained with hematoxylin and eosin to confirm the pathological diagnosis of the tumor and normal tissues. Consecutive tissue sections were subjected to IHC staining. Antigen retrieval was performed by boiling the slides in sodium citrate solution. Then, 3% H<sub>2</sub>O<sub>2</sub> was used to eliminate the endogenous peroxidase activity. Next, tissue sections were incubated with mouse anti-<italic>E. coli</italic> lipopolysaccharide (LPS) antibodies (1:200 dilution; Abcam, cat# ab35654) overnight at 4&#x2103; followed by incubation with horseradish peroxidase-conjugated anti-mouse secondary antibodies. Color development was performed using diaminobenzidine. The tissue sections were counterstained with hematoxylin, dehydrated, mounted, and observed under a light microscope.</p>
</sec>
<sec id="s2_4">
<label>2.4</label>
<title>Isolation of tumor cell-associated microbiota of CRC</title>
<p>Tissues from 11 patients (see section 2.2) was enzymatically processed to obtain single cells (TC and NC), respectively. As shown in <xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1B</bold>
</xref>, the tissue was cut into ~1 mm<sup>3</sup> pieces and digested in a 4 mL solution with 1 mg/mL dispase (CAS: 42613-33-2, D6430, Solarbio) and 1 mg/mL collagenase (CAS: 9001-12-1, C8160, Solarbio) for 30 min at 37&#xb0;C. After enzymatic digestion, the samples were diluted in Dulbecco&#x2019;s modified Eagle&#x2019;s medium (DMEM) containing 20% fetal bovine serum. The cell suspension was filtered through a 40-&#x3bc;m cell strainer (BD 352340) to collect the single-cell suspension. After centrifugation at 300 x g for 5 min, the supernatant was discarded. The cell pellet was washed three times with 5 mL phosphate-buffered saline. Finally, the cells were resuspended in 500 &#xb5;L serum-free DMEM and frozen in liquid nitrogen for subsequent 5R-16s rRNA gene sequencing.</p>
</sec>
<sec id="s2_5">
<label>2.5</label>
<title>5R-16s rRNA gene sequencing</title>
<p>5R-16s rRNA gene sequencing was performed as described previously (<xref ref-type="bibr" rid="B10">10</xref>) (<xref ref-type="bibr" rid="B5">5</xref>). DNA from frozen samples were extracted with the CTAB method using kit DP302-02 (TianGen, Beijing, China) and were quantified using Qubit (Invitrogen, USA). 16S rRNA gene amplification and sequencing were performed by amplifying five (V2, V3, V5, V6, and V8) regions of the 16S rRNA gene in a multiplex (Primer sequences are described in <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Material S1</bold>
</xref>). PCR products were purified with AMPure XT beads (Beckman Coulter Genomics, Danvers, MA, USA) and quantified with Qubit (Invitrogen, USA). Purified PCR products were evaluated using an Agilent 2100 Bioanalyzer (Agilent, USA) and Illumina (Kapa Biosciences, Woburn, MA, USA) library quantification kit. The qualified library concentration should be above 2 nM. Sequencing was performed using NovaSeq 6000 sequencer with NovaSeq 6000 SP Reagent Kit (500 cycles). The Short Multiple Regions Framework (<xref ref-type="bibr" rid="B11">11</xref>) was used to analyze the sequences of the five amplified regions for bacterial taxonomic identification. We used the Greengenes database for this project (ver. May 2013). To reduce the effect of low-abundance noise on subsequent analyses, the number of sequence reads per sample was normalized to remove samples with less than 1000 total reads (including negative controls) and bacterial data with a relative abundance of less than 10<sup>-4</sup>. A total of 19 empty tubes were set up as negative control tubes during the DNA isolation and sequencing procedures, which were used to exclude contaminated bacteria together with the 3 environmental control tubes. Considering the low load biomass of the samples, bacteria introduced from the sampling environment and sequencing experiments will cause big interference to the microbiota of the samples, thus environmental controls (an empty tube was set up in parallel to the experimental tubes during sample collection procedures) and negative controls (an empty tube set up in parallel to the experimental tubes during DNA isolation and sequencing procedures) were set up. Based on the prevalence of bacteria in the control tubes, the contaminating bacteria at the sampling end and the sequencing end were determined. In this project, more than 50% prevalence was set as the threshold for identification of contaminating bacteria, that is, when &#x2265; 2 of the 3 environmental control tubes or &#x2265; 9 of the 19 negative controls showed the species, the species was considered as contaminated bacteria and was excluded from analyses. After removing contaminating bacteria, diversity analysis and differential flora identification within and between groups were performed based on the filtered flora data.</p>
</sec>
<sec id="s2_6">
<label>2.6</label>
<title>Data analysis</title>
<p>In the alpha diversity analysis, species richness and evenness were measured, as well as sequencing depth. The Chao1 and Observed_species indices reflected the species richness of the samples. Good coverage reflected the flora coverage of the sample. The Simpson and Shannon indices were mainly used to reflect the richness and evenness of species. Beta diversity analysis was performed to analyze species differences between groups. Principal component analysis (PCA), principal coordinate analysis (PCoA), non-metric multidimensional scaling (NMDS), analysis of similarities (ANOSIM) and other methods were utilized to observe differences between samples. The ratios of the abundance of each bacterial species to the abundance of all species at each taxonomic level (relative abundance) were calculated, and stacked bar charts of bacterial species distribution in different samples (or groups) at each taxonomic level were drawn to analyze and compare the changes in bacterial composition between different groups. A Venn diagram was used to analyze the number of common and unique species across multiple groups at each taxonomic level to visually show the similarity and specificity of species composition in different groups. Linear discriminant analysis effect size (LEfSe) (<xref ref-type="bibr" rid="B12">12</xref>) was used to analyze differences in the microbiota compositions. In this study, the threshold of the LEfSe analysis was set at a Linear Discriminate Analysis (LDA) value &gt; 3, <italic>p</italic> &lt; 0.05. The results are presented as a table of potential bacterial markers, a bar chart of the LDA value distribution, and an evolutionary branch map. The Mann Whitney U test was used to analyze significant differences between groups (<xref ref-type="bibr" rid="B13">13</xref>). Microbiota function was predicted using the latest version of Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) 2 (<xref ref-type="bibr" rid="B13">13</xref>). The results of gene functional annotation of the microbiota from databases, including COG, EC, KO, PFAM, and TIGRFAM profiling were obtained, and then statistical analysis of metagenomic profiles (STAMP) (<xref ref-type="bibr" rid="B14">14</xref>) was used for differential analysis to obtain significantly different gene functions between groups. In all these analyses, Fisher_test was used for between-sample statistics, with a <italic>p</italic> &lt; 0.05 considered statistically significant. Statistics of multiple replicate samples were performed using the rank sum test, with Wilcox_test for two-group comparisons, Kruskal_test for multi-group comparisons, and <italic>p</italic> &lt; 0.05 considered statistically significant.</p>
</sec>
</sec>
<sec id="s3" sec-type="results">
<label>3</label>
<title>Results</title>
<sec id="s3_1">
<label>3.1</label>
<title>Distribution of gut microbes</title>
<p>IHC analysis showed that LPS was mainly distributed in the feces found on the surface of intestinal tissues and epithelial cells and macrophages within the tissues (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1C</bold>
</xref>).</p>
</sec>
<sec id="s3_2">
<label>3.2</label>
<title>Influence of sampling site and cleaning on microbial detection</title>
<p>The compositions of microbiota were different among the six samples from the 12<sup>th</sup> patient (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2</bold>
</xref>). Samples E, F, and G were adjacent normal tissue samples (equivalent to part &#x2463; in <xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1A</bold>
</xref>) containing the contaminated intestinal mucosa, which were washed different times as shown in <xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2E</bold>
</xref>. Their respective species compositions are similar. Sample A was the tumor tissue sample (equivalent to part &#x2460; in <xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1A</bold>
</xref>) containing the contaminated intestinal mucosa. Sample B (equivalent to part &#x2461; in <xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1A</bold>
</xref>) and sample C (equivalent to part &#x2462; in <xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1A</bold>
</xref>) had the mucosa removed and different washing numbers (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2E</bold>
</xref>). There was little difference between the bacterial composition of A and that of E, F, and G (<xref ref-type="fig" rid="f2">
<bold>Figures&#xa0;2A, B</bold>
</xref>), indicating that their microbiomes were similar. However, there were significant differences between the species compositions of A and B or C (<italic>p</italic> &lt; 0.05). Firmicutes and Bacteroidetes were the dominant species in E, F, and G, whereas the relative abundance of Cyanobacteria was dominant in B and C (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2C</bold>
</xref>). At the genus level, the abundance of <italic>Prevotella, Faecalibacterium</italic>, and <italic>Bacteroides</italic> were reduced in B and C compared to that in the other samples, whereas <italic>Unknown_family_Unknown_genus39</italic> and <italic>Vibrio</italic> were present, and the relative abundance of <italic>Leptotrichia</italic> was increased in TT compared to that in NT (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2D</bold>
</xref>).</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Microbiota abundance with different sampling sites and cleaning conditions. These results were obtained from the analysis of six samples from the 12th patient with different treatment conditions. <bold>(A)</bold> PCA, the more similar the species composition of the samples, the closer they are in the PCA plot. <bold>(B)</bold> UPGMA analysis, the shorter the branch length between samples, the higher the similarity of the two samples. <bold>(C)</bold> Microflora distribution at the phylum level. <bold>(D)</bold> Microflora distribution at the genus level. <bold>(E)</bold> Profiles of the samples. PCA, principal component analysis; UPGMA, Unweighted Pair-Group Method with Arithmetic Mean.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fonc-14-1374769-g002.tif"/>
</fig>
</sec>
<sec id="s3_3">
<label>3.3</label>
<title>Tumor cell-associated microbiota obtained with cell enzymatic digestion combined with microbial detection</title>
<sec id="s3_3_1">
<label>3.3.1</label>
<title>Quality control of sample sequencing data</title>
<p>5R-16s RNA gene analysis was performed on 44 samples from 11 patients with CRC, including CRC tumor tissues (TT), normal tissues adjacent to CRC (NT), tumor cells (TC), normal cells (NC). TC and NC were obtained through separation and extraction of TT and NT tissues using enzymatic digestion. A total of 4,606,353 raw reads (&gt; 100,000 reads per sample) and 4,589,247 clean reads (&#x2265; 99588 reads per sample) were obtained from 44 samples, with an average Q20 of 98.71% and Q30 of 96.19% (<xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Exl S2</bold>
</xref>).</p>
</sec>
<sec id="s3_3_2">
<label>3.3.2</label>
<title>alpha diversity analysis</title>
<p>In the alpha diversity analysis, the alpha rarefaction curve showed a reasonable sequencing depth and a good value of 1 (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3A</bold>
</xref>), indicating complete coverage of microbiota sequencing. Chao1 and observed_species were used to estimate the number of species in a community. The results showed that the number of species in the cell groups after using the enzymatic digestion method was significantly lower than that in the tissue groups (<italic>p</italic> &lt; 0.05). The ratio of the microbe number in cells to that of tissues was higher in tumors than that in normal tissues (<xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>; <xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3B</bold>
</xref>) (Additional rarefaction curves are shown in <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Figure S3</bold>
</xref>).</p>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>Alpha diversity analysis. Alpha diversity reflects species richness and evenness, as well as sequencing depth. <bold>(A)</bold> Rarefaction curves of good index. When the curve tends to be flat, it indicates that the amount of sequencing data is rich enough, and further increasing the amount of data will contribute little to the discovery of more species. Good coverage reflected the flora coverage of the sample. <bold>(B)</bold> violin plot: the Chao1 and Observed analysis reflected the species richness of the samples, the Simpson and Shannon analyses are mainly reflected the richness and evenness of species. *: <italic>p&lt; 0.05</italic>, **: <italic>p&lt; 0.01</italic>, ***: <italic>p&lt; 0.001</italic>.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fonc-14-1374769-g003.tif"/>
</fig>
<table-wrap id="T1" position="float">
<label>Table&#xa0;1</label>
<caption>
<p>alpha_diversity_statistics.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="center">alpha_diversity</th>
<th valign="middle" align="center">TC</th>
<th valign="middle" align="center">NC</th>
<th valign="middle" align="center">TT</th>
<th valign="middle" align="center">NT</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="center">goods_coverage</td>
<td valign="middle" align="center">1.00&#xb1;0.00</td>
<td valign="middle" align="center">1.00&#xb1;0.00</td>
<td valign="middle" align="center">1.00&#xb1;0.00</td>
<td valign="middle" align="center">1.00&#xb1;0.00</td>
</tr>
<tr>
<td valign="middle" align="center">chao1</td>
<td valign="middle" align="center">93.06&#xb1;70.65</td>
<td valign="middle" align="center">86.09&#xb1;52.34</td>
<td valign="middle" align="center">210.29&#xb1;75.97</td>
<td valign="middle" align="center">236.37&#xb1;85.83</td>
</tr>
<tr>
<td valign="middle" align="center">observed_species</td>
<td valign="middle" align="center">92.91&#xb1;70.41</td>
<td valign="middle" align="center">86.09&#xb1;52.34</td>
<td valign="middle" align="center">209.82&#xb1;75.67</td>
<td valign="middle" align="center">235.91&#xb1;85.68</td>
</tr>
<tr>
<td valign="middle" align="center">shannon</td>
<td valign="middle" align="center">3.02&#xb1;2.03</td>
<td valign="middle" align="center">2.94&#xb1;2.15</td>
<td valign="middle" align="center">4.95&#xb1;1.05</td>
<td valign="middle" align="center">5.52&#xb1;0.75</td>
</tr>
<tr>
<td valign="middle" align="center">simpson</td>
<td valign="middle" align="center">0.63&#xb1;0.32</td>
<td valign="middle" align="center">0.59&#xb1;0.33</td>
<td valign="middle" align="center">0.89&#xb1;0.10</td>
<td valign="middle" align="center">0.94&#xb1;0.02</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s3_3_3">
<label>3.3.3</label>
<title>Beta diversity analysis</title>
<p>Beta diversity reflects species differences between environmental communities. Beta and alpha diversities constitute the overall diversity or biological heterogeneity of an environmental community. Beta diversity analysis mainly measures the differences between samples using PCA, PCoA, the Unweighted Pair-Group Method with Arithmetic mean (UPGMA), ANOSIM, and other methods. The results of PCA and PCoA showed no obvious differences between the TT and NT groups or between the TC and NC groups, but there were differences between the cells and tissues (<xref ref-type="fig" rid="f4">
<bold>Figures&#xa0;4A, B</bold>
</xref>). The UPGMA analysis results show that the microbiomes of eight out of 11 patients had few differences between tumor and normal tissue samples. However, after enzymatic hydrolysis, only the microbiomes of the tumor and normal cell samples of three patients were similar, whilst the microbiomes of the cell samples from most patients showed obvious differences (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4C</bold>
</xref>). This indicates that the data for each group were consistent and that the enzymatic hydrolysis of cells was effective. The ANOSIM results showed a significant difference in microbiota composition between the cell and tissue groups (<italic>R &gt; 0</italic>, <italic>p &lt; 0.05</italic>) (<xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref>).</p>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>Beta diversity analysis. Beta diversity reflects species differences between groups. <bold>(A)</bold> The more similar the species composition of the samples, the closer they are in the PCA plot. <bold>(B)</bold> A closer distance between two points indicates a smaller difference in community composition in the PCoA plot. <bold>(C)</bold> The shorter the branch length between samples, the more similar two samples are. PCA, principal component analysis; PCoA, principal coordinate analysis.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fonc-14-1374769-g004.tif"/>
</fig>
<table-wrap id="T2" position="float">
<label>Table&#xa0;2</label>
<caption>
<p>anosim_result.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="center">Method name</th>
<th valign="middle" align="center">
<italic>R</italic> statistic</th>
<th valign="middle" align="center">
<italic>p</italic>_value</th>
<th valign="middle" align="center">Groups</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="center">bray_curtis</td>
<td valign="middle" align="center">0.56</td>
<td valign="middle" align="center">0.00</td>
<td valign="middle" align="center">NC_vs_NT</td>
</tr>
<tr>
<td valign="middle" align="center">bray_curtis</td>
<td valign="middle" align="center">-0.03</td>
<td valign="middle" align="center">0.72</td>
<td valign="middle" align="center">TC_vs_NC</td>
</tr>
<tr>
<td valign="middle" align="center">bray_curtis</td>
<td valign="middle" align="center">0.47</td>
<td valign="middle" align="center">0.00</td>
<td valign="middle" align="center">TC_vs_TT</td>
</tr>
<tr>
<td valign="middle" align="center">bray_curtis</td>
<td valign="middle" align="center">-0.07</td>
<td valign="middle" align="center">0.93</td>
<td valign="middle" align="center">TT_vs_NT</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s3_3_4">
<label>3.3.4</label>
<title>Analysis of microbiota</title>
<p>The distribution of the microbiota differed among the groups at the phylum and genus levels. A total of 22 phyla and 518 genera were detected in the TT samples, 24 phyla and 594 genera in the NT samples, 21 phyla and 367 genera in the TC samples, and 24 phyla and 357 genera in the NC samples (<xref ref-type="fig" rid="f5">
<bold>Figures&#xa0;5A, B</bold>
</xref>). Microbial species that were present in more than half of the 11 patient samples were further analyzed and found that eleven phyla and 59 genera were detected in TT samples, 12 phyla and 79 genera in NT samples, 8 phyla and 37 genera in TC samples, and 9 phyla and 33 genera in NC samples (<xref ref-type="fig" rid="f5">
<bold>Figures&#xa0;5C, D</bold>
</xref>). At the phylum level (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5E</bold>
</xref>), Firmicutes, Bacteroidetes, and Proteobacteria were the most abundant in the tissue samples, followed by Actinobacteria and Fusobacteria. Notably, Firmicutes and Bacteroidetes were highly abundant in tissues, whereas Proteobacteria were the most abundant in the treated cell samples. The abundance of Firmicutes was lower in the tumor tissues than that in the normal tissue groups but was higher in the tumor cell groups. There was no significant difference in Actinobacteria abundance between the tumor tissue and normal tissue groups (<italic>p &gt; 0.05)</italic>. At the genus level, <italic>Pseudomonas, Bacteroides, Brevundimonas</italic>, and <italic>Lactobacillus</italic> were predominant in the TC group. <italic>Pseudomonas, Brevundimonas, Bacteroides</italic>, and <italic>Streptococcus</italic> were the dominant bacteria in the NC group, whereas <italic>Bacteroides, Prevotella, Fusobacterium</italic>, and others were the dominant bacteria in the TT group. The dominant bacteria in the NT group were <italic>Bacteroides, Prevotella</italic>, and <italic>Faecalibacterium</italic> (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5F</bold>
</xref>). The percentages of the top 10 species are shown in the pie chart in <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Figure S4</bold>
</xref>.</p>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>The distribution of microbiota among the groups at the phylum and genus levels. <bold>(A)</bold> Venn analysis of all species at the phylum level. <bold>(B)</bold> Venn analysis of all species at the genus level. <bold>(C, D)</bold> Venn analysis of species present in more than half of the samples at the phylum and genus levels. <bold>(E, F)</bold> Column stacking charts of species abundance at the phylum and genus levels. TC, TT, NC, and NT are the four groups. NC, normal cell; NT, normal tissue; TC, tumor cell; TT, tumor tissue.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fonc-14-1374769-g005.tif"/>
</fig>
</sec>
<sec id="s3_3_5">
<label>3.3.5</label>
<title>Analysis of biomarkers</title>
<p>To identify differences in the compositions of the microbiota among groups, the LEfSe algorithm was used. Using LDA&gt;3.0 as a threshold for discriminative features, as shown in the cladogram and LDA score bar graph (<xref ref-type="fig" rid="f6">
<bold>Figures&#xa0;6A, B</bold>
</xref>), we identified significant differences in the abundances of Veillonellaceae and Tissierellaceae at the family level and <italic>Eikenella</italic> at the genus level (<italic>p &lt; 0.05</italic>), while these were more abundant in the TC group than in the NC group. Additionally, <italic>Eikenella</italic> was more abundant in TT than in NT. Identical results were obtained using the Mann&#x2013;Whitney U test (<xref ref-type="fig" rid="f6">
<bold>Figures&#xa0;6C, D</bold>
</xref>), and <italic>Eikenella</italic> remained significantly more abundant in the TT and TC groups than in the normal groups (<italic>p &lt; 0.05</italic>). The results of TC versus TT and NC versus NT are shown in <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Figures S5~10</bold>
</xref>. Comparison among the four groups revealed that the abundance of some bacteria species was higher in the tumor groups than that in the normal groups, such as <italic>Fusobacterium, Eikenella, Shewanella, and Listeria</italic>, whereas the abundance of <italic>Akkermansia</italic> was lower in the tumor group than that in the normal group. The ratios of some bacteria species in cells to that in tissues, such as <italic>Gemella, Escherichia, Shigella</italic>, and <italic>Blautia</italic>, were remarkable changed (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6E</bold>
</xref>).</p>
<fig id="f6" position="float">
<label>Figure&#xa0;6</label>
<caption>
<p>Analysis of biomarkers. LEfSe analysis is the most common method for differential analysis of microbiota. LDA<italic>&gt;3.0, p&lt; 0.05</italic> were used as a threshold for discriminative features. <bold>(A, B)</bold> Cladogram and LDA score bar graph. <bold>(C, D)</bold> Mann&#x2013;Whitney U test analysis. <bold>(E)</bold> Abundance maps of some of the different species. TC, TT, NC, and NT are the four groups. LEfSe, linear discriminant analysis effect size; NC, normal cell; NT, normal tissue; TC, tumor cell; TT, tumor tissue.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fonc-14-1374769-g006.tif"/>
</fig>
</sec>
<sec id="s3_3_6">
<label>3.3.6</label>
<title>Predictive functional profiling of microbiota</title>
<p>PICRUSt2 was used to explore the potential functional profiles of the CRC microbiota from the four groups. The COG results showed significant differences in the levels of predicted Zn-dependent protease, minimal metalloprotease (MMP)-like domain, dienelactone hydrolase, regulator of RNase E activity, and RraA between the TC and NC groups (<italic>p &lt; 0.05</italic>); Uncharacterized conserved protein YkwD, contains CAP (CSP/antigen 5/PR1) domain, phosphoenolpyruvate carboxykinase, GTP-dependent uncharacterized protein, pyridoxamine 5&#x2019;-phosphate oxidase (PNPOx-like) family, high-affinity nickel permease, HD-GYP domain, and c-di-GMP phosphodiesterase class II (or its inactivated variant) levels were significantly different between the TT and NT groups (<italic>p &lt; 0.05</italic>) (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7A</bold>
</xref>). KO results showed that significant differences in the levels of pilus assembly proteins, CpaE, phoN, acid phosphatase (class A), tight adherence protein C, and putative membrane proteins between the TC and NC groups (<italic>p &lt; 0.05</italic>); formyltetrahydrofolate deformylase, type IV pilus assembly protein PilA, and holin-like protein LrgB levels were significantly different between the TT and NT groups (<italic>p &lt; 0.05</italic>) (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7B</bold>
</xref>). The pathway analysis results showed that only the function of pyruvate fermentation to isobutanol (engineered) was significantly different between the TC and NC groups (<italic>p &lt; 0.05</italic>) (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7C</bold>
</xref>); there was no difference in the functions between the TT and NT groups. The results for other groups are provided in <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Figures S11, S12</bold>
</xref>.</p>
<fig id="f7" position="float">
<label>Figure&#xa0;7</label>
<caption>
<p>Predictive functional profiling of microbiota. PICRUSt2 was performed to explore the potential functional profiling of CRC microbiota from the four groups. <bold>(A)</bold> COG database results. <bold>(B)</bold> KO database results. <bold>(C)</bold> Pathway database results. TC, TT, NC, and NT are the four groups. CRC, colorectal cancer; NC, normal cell; NT, normal tissue; TC, tumor cell; TT, tumor tissue.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fonc-14-1374769-g007.tif"/>
</fig>
</sec>
</sec>
</sec>
<sec id="s4" sec-type="discussion">
<label>4</label>
<title>Discussion</title>
<p>The intestinal microecology plays an important role in the occurrence and development of CRC. In recent years, the development of sequencing technology and bioinformatics has enabled the study of the composition, metabolism, and metagenome of intestinal microecology (<xref ref-type="bibr" rid="B14">14</xref>). The gut microbiota has various effects on the transformation process, progression, and response of tumors to immunotherapy. In addition to multiple functions, such as maintaining the epithelial barrier, immune response, digestive function, and regulating neurotransmitters, the gut microbiota are essential for maintaining homeostasis. Changes in the gut microbiota may interfere with homeostasis, leading to dysregulation of the flora (<xref ref-type="bibr" rid="B15">15</xref>). Microbiota dysregulation manifests as changes in composition, decreases in diversity, and increases in potential pathogenic bacteria (<xref ref-type="bibr" rid="B16">16</xref>, <xref ref-type="bibr" rid="B17">17</xref>). Dysregulated intestinal microecology has &#x201c;pre-carcinogenic properties&#x201d; and contributes to the formation of a tumor-like microenvironment. Intestinal epithelial function impairment, pathogen recognition dysfunction, and abnormal immune response accompanied by the dysregulation of the microbiome lead to a severe reduction in intestinal resistance to colonization by pathogenic microorganisms, thus making it easier for pathogenic bacteria to invade intestinal epithelial tissues and cells (<xref ref-type="bibr" rid="B18">18</xref>). Additionally, the dysregulation of microbiota can promote the proliferation, invasion and metastasis of tumor cells by increasing the release of endotoxin and mediating the M2 polarization of TLR4-dependent macrophages (<xref ref-type="bibr" rid="B19">19</xref>, <xref ref-type="bibr" rid="B20">20</xref>). Furthermore, dysregulated microbiota can have a significant impact on chemotherapy efficacy by affecting drug absorption, decomposition, and toxicity, thus affecting the development and prognosis of CRC (<xref ref-type="bibr" rid="B15">15</xref>).</p>
<p>While many studies have revealed a strong relationship between gut microbiota and CRC, most have focused on fecal microbiota and their metabolites. In 2020, Nejman et&#xa0;al. (<xref ref-type="bibr" rid="B5">5</xref>) studied 1,526 tumors and adjacent normal tissues from seven cancer types: breast, lung, ovarian, pancreatic, melanoma, bone, and brain tissue tumors and found that tumor microbiota are composed of tumor-specific intracellular bacteria. In addition, intracellular bacteria can cause cells to become cancerous by hybridizing their own DNA with that of the host cell (<xref ref-type="bibr" rid="B21">21</xref>). Host cells infected with bacteria undergo changes in metabolism, exhibiting a &#x201c;Warburg-like metabolism&#x201d; similar to that of cancer cells (<xref ref-type="bibr" rid="B22">22</xref>). After bacteria enter cancer tissues or cells, they not only affect chemotherapy efficacy through drug transformation, but also antibiotic resistance through immune evasion (<xref ref-type="bibr" rid="B23">23</xref>). Bacteria colonize tumors in humans, proliferate within them, and modulate immune function, ultimately affecting the survival of patients with cancer and their response to treatment (<xref ref-type="bibr" rid="B24">24</xref>). In 2022, Fu et&#xa0;al. found that intracellular bacteria enhanced the resistance of circulating tumor cells to fluid shear stress by reorganizing the cytoskeleton, improving the survival ability of host cells, and thus promoting tumor metastasis (<xref ref-type="bibr" rid="B25">25</xref>). In summary, intratumoral microbiota play an important and direct role in inducing tumor occurrence, altering cell metabolism, promoting tumor metastasis, and enabling the tumor to evade immune surveillance. The in-depth study and subsequent regulation of intratumoral microbiota may aid the discovery of novel treatment options for patients with cancer (<xref ref-type="bibr" rid="B5">5</xref>).</p>
<p>However, given the low microbial biomass of tumors, characterization of tumor microbiome remains a challenge. Tumors have a high host-to-bacterial DNA ratio, which makes the identification of tumor-associated bacteria extremely complex (<xref ref-type="bibr" rid="B26">26</xref>). Nevertheless, with the improvements in new-generation analytic tools, there has been some progress in the study of intratumoral microbiomes in recent years. Moreover, it has now been proven that each cancer subtype has a unique microbiome, characterized by bacterial communities with specific metabolic functions (<xref ref-type="bibr" rid="B5">5</xref>, <xref ref-type="bibr" rid="B26">26</xref>).</p>
<p>However, in contrast to most tumors, such as breast cancer, which have no direct connection with the external environment, CRC tumor microbes are contaminated by fecal microbes owing to the structure of the gut folds and villi. The low microbial biomass within the tumor, the large and complex intestinal microbiome, and dealing with contamination, are all challenges in detecting tumor-associated microbiota in CRC. These difficulties have limited research on intratumoral microbes in the tumors that are directly related to the external environment. In a study by Nejman et&#xa0;al., they included 22 colorectal tumors to enrich their analysis. Unlike that of other cancers, such as lung and ovarian cancer, the Jaccard index of bacterial species profiles was similar between CRC and normal tissues, and bacteria belonging to the Firmicutes and Bacteroidetes phyla were the most abundant species in colorectal tumors (<xref ref-type="bibr" rid="B5">5</xref>). These results were consistent with those obtained in our tissue-level study and a previous study (<xref ref-type="bibr" rid="B27">27</xref>). In another study, Dohlman et&#xa0;al. distinguished tissue-resident microbiota from contaminants. It was found that the Firmicutes and Bacteroidetes were the most dominant phyla in CRC tissues, and bacteria in CRC tissues were enriched for mucosa-related species (<xref ref-type="bibr" rid="B28">28</xref>). In the four major gastrointestinal cancer microbiome studies, no differences in microbial community profiles were observed between tumor and normal samples (<xref ref-type="bibr" rid="B29">29</xref>).</p>
<p>LPS is a component of the cell wall in Gram-negative bacteria. Most intestinal bacteria are Gram-negative, such as <italic>Escherichia coli</italic>. Bacteria were mainly distributed in the feces found on the surface of intestinal tissues and in epithelial cells and macrophages within the tissues. This is consistent with the findings of a previous report (<xref ref-type="bibr" rid="B4">4</xref>). We attempted to avoid contamination by intestinal fecal bacteria through cleaning the sampling sites and controlling their location. The effects of different sampling sites and number of washes on the detection of CRC intratumoral microbiota were compared. The results showed that different sampling sites and number of washes affected microbiota detection. Although cleaning methods can be controlled, sampling sites varied and led to different microbial distributions. Carcinogenesis is related to the bacterial composition inside the intestinal epithelial cells; however, in the normal intestinal tissues, the contamination of fecal bacteria cannot be avoided due to the tissue structure. Therefore, fecal bacterial contamination cannot be avoided by controlling or cleaning the sampling site in sampling normal tissues.</p>
<p>We established a new method for isolating tumor cell-associated microbiota using cell enzymatic digestion. First, we destroyed the gut folds and villus structures by rough shearing and cell enzymatic hydrolysis to exclude fecal bacterial contamination. In this process, we chose dispase and collagenase, instead of trypsin, in order to destroy the intestinal fold and villus structure, whilst preserving the integrity of the cells. Then, we used a 40-&#x3bc;m cell sieve to filter out the undigested tissue blocks to further avoid contamination of fecal bacteria in the tissue sample. After that, the collected single cell suspension was separated repeatedly with a centrifugal force of 300 x g, to remove the &#x201c;passing bacteria&#x201d; where the binding between the cells are bacteria is weak. This process helps to reduce fecal bacterial contamination. The final microbiome consisted mainly of bacteria from the cells or those that colonized the cells.</p>
<p>Finally, 5R-16s rRNA gene sequencing was used to analyze and compare the differences in CRC-associated microbiota of a total of 44 conventional and enzymatically treated samples from 11 patients. The sequencing results of the tissue samples were consistent with those of previous studies (<xref ref-type="bibr" rid="B5">5</xref>, <xref ref-type="bibr" rid="B27">27</xref>). However, different results were obtained for the cell samples treated with enzymatic hydrolysis. In contrast to the that at the tissue level, Proteobacteria were the most abundant species in the cell samples. The abundance of the Firmicutes phyla was lower in the tumor tissues than that in the normal tissues but was higher in the cell groups. The proportion of Actinobacteria also differed between the two groups. The cell enzymatic digestion method may reduce the contamination of fecal bacteria, enabling low biomass intratumoral microbiota to be obtained. This was also confirmed at the genus level, where more differential bacteria were found by comparing the abundance and tumor/normal ratios of species.</p>
<p>Combining the detection results from tissues can also predict bacterial distributions. When the biomass of a species is high in both tumor tissues and cells, the bacteria are enriched and have invaded the tumor cells. If there is no difference in the species abundance between the tumor and normal tissues, the intracellular biomass increases, indicating that the bacteria have invaded the tumor. Conventional detection methods can only detect the enrichment of bacteria in tissues, but not the movement of bacteria.</p>
<p>Cell enzymatic digestion can also be combined with microbial culture to detect tumor cell-associated microorganisms, which may be more accurate than tissue culture and narrows the scope of future research. However, owing to the limitations of laboratory conditions, this approach has not been investigated and requires further research. There is a scarcity of research and skillset for tumor research in the microbial field. Therefore, it is difficult to conduct further timely research and comparisons between studies. Due to the difficulty of acquiring clinical biopsy tissues, the sample size is relatively small. Further studies will increase the sample size and expand this study, including using animal models to explore further research, which may identify microbial targets and provide help for future CRC diagnosis and intervention.</p>
<p>This study provides a new method to detect CRC cell-associated microbiota. However, our study has some limitations: first, the sample size is small. A total of 12 patient samples were used and only 11 patient samples were used for full analyses. Second, the method of isolating tumor-associated microbiome was the first of its kind that has not been used by other investigators, the reproducibility of which needs to be validated in another cohort of samples. Third, fecal microbiome were not analyzed using the same 5R-16s rRNA-seq because the abundance of fecal microbiome was way more than the tumor-associated microbiome. Future studies shall use more samples and multiple cohorts of patients including at least a testing cohort and a validation cohort to verify the reproducibility of the study. It is desirable to analyze both fecal and tumor-associated microbiomes and compare the microbial populations between them.</p>
</sec>
<sec id="s5" sec-type="conclusions">
<label>5</label>
<title>Conclusion</title>
<p>Different sampling sites and number of washes led to different microbiota detection results. The cell enzymatic digestion method can reduce the coverage of fecal bacteria in intestinal tissues to allow low biomass intratumoral microbiota to be obtained. This is a new method for the isolation of tumor cell-associated microbiota, and in combination with microbial composition analysis of tissues can predict the movement of bacteria. It can also provide additional assistance for generating microbial cultures. Thus, this study provides a new method for studying tumor-associated microbiota.</p>
</sec>
<sec id="s6" sec-type="data-availability">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/<xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Material</bold>
</xref>. Further inquiries can be directed to the corresponding authors.</p>
</sec>
<sec id="s7" sec-type="ethics-statement">
<title>Ethics statement</title>
<p>All participants provided informed consent, and the study was approved by the Medical Ethics Committee of Chengde Medical University (No. 202215). The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants&#x2019; legal guardians/next of kin.</p>
</sec>
<sec id="s8" sec-type="author-contributions">
<title>Author contributions</title>
<p>YZZ: Formal analysis, Writing &#x2013; original draft. YJL: Funding acquisition, Writing &#x2013; original draft. JYP: Methodology, Writing &#x2013; original draft. SKJ: Methodology, Writing &#x2013; review &amp; editing. XYZ: Data curation, Writing &#x2013; original draft. EHZ: Resources, Writing &#x2013; review &amp; editing. YHL: Funding acquisition, Writing &#x2013; review &amp; editing.</p>
</sec>
</body>
<back>
<sec id="s9" sec-type="funding-information">
<title>Funding</title>
<p>The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was funded by the Natural Science Foundation of Hebei Province, China (H2019406073, H2020406008), Central guidance for local scientific and technological development funds of Hebei Province (236Z7701G), the Technology Innovation Guidance Project-Science and Technology Work Conference of Hebei Provincial Department of Science and Technology, Chengde Medical College Discipline Construction Funds, and the University Research Project of the Chengde Medical College (202007, KY202310, KY202402).</p>
</sec>
<ack>
<title>Acknowledgments</title>
<p>We thank the Affiliated Hospital of Chengde Medical College for the support of this project. We would like to thank Editage (<ext-link ext-link-type="uri" xlink:href="http://www.editage.cn">www.editage.cn</ext-link>) for English language editing.</p>
</ack>
<sec id="s10" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec id="s11" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors&#xa0;and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="s12" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fonc.2024.1374769/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fonc.2024.1374769/full#supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="Presentation_1.zip" id="SM1" mimetype="application/zip"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<label>1</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sung</surname> <given-names>H</given-names>
</name>
<name>
<surname>Ferlay</surname> <given-names>J</given-names>
</name>
<name>
<surname>Siegel</surname> <given-names>RL</given-names>
</name>
<name>
<surname>Laversanne</surname> <given-names>M</given-names>
</name>
<name>
<surname>Soerjomataram</surname> <given-names>I</given-names>
</name>
<name>
<surname>Jemal</surname> <given-names>A</given-names>
</name>
<etal/>
</person-group>. <article-title>Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries</article-title>. <source>CA Cancer J Clin</source>. (<year>2021</year>) <volume>71</volume>:<page-range>209&#x2013;49</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.3322/caac.21660</pub-id>
</citation>
</ref>
<ref id="B2">
<label>2</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bastiaanssen</surname> <given-names>T</given-names>
</name>
<name>
<surname>Cowan</surname> <given-names>C</given-names>
</name>
<name>
<surname>Claesson</surname> <given-names>MJ</given-names>
</name>
<name>
<surname>Dinan</surname> <given-names>TG</given-names>
</name>
<name>
<surname>Cryan</surname> <given-names>JF</given-names>
</name>
</person-group>. <article-title>Making sense of &#x2026; the microbiome in psychiatry</article-title>. <source>Int J Neuropsychopharmacol</source>. (<year>2019</year>) <volume>22</volume>:<fpage>37</fpage>&#x2013;<lpage>52</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/ijnp/pyy067</pub-id>
</citation>
</ref>
<ref id="B3">
<label>3</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dougherty</surname> <given-names>MW</given-names>
</name>
<name>
<surname>Jobin</surname> <given-names>C</given-names>
</name>
</person-group>. <article-title>Intestinal bacteria and colorectal cancer: etiology and treatment</article-title>. <source>Gut Microbes</source>. (<year>2023</year>) <volume>15</volume>:<elocation-id>2185028</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/19490976.2023.2185028</pub-id>
</citation>
</ref>
<ref id="B4">
<label>4</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Galeano Ni&#xf1;o</surname> <given-names>JL</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>H</given-names>
</name>
<name>
<surname>LaCourse</surname> <given-names>KD</given-names>
</name>
<name>
<surname>Kempchinsky</surname> <given-names>AG</given-names>
</name>
<name>
<surname>Baryiames</surname> <given-names>A</given-names>
</name>
<name>
<surname>Barber</surname> <given-names>B</given-names>
</name>
<etal/>
</person-group>. <article-title>Effect of the intratumoral microbiota on spatial and cellular heterogeneity in cancer</article-title>. <source>Nature</source>. (<year>2022</year>) <volume>611</volume>:<page-range>810&#x2013;7</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41586-022-05435-0</pub-id>
</citation>
</ref>
<ref id="B5">
<label>5</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nejman</surname> <given-names>D</given-names>
</name>
<name>
<surname>Livyatan</surname> <given-names>I</given-names>
</name>
<name>
<surname>Fuks</surname> <given-names>G</given-names>
</name>
<name>
<surname>Gavert</surname> <given-names>N</given-names>
</name>
<name>
<surname>Zwang</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Geller</surname> <given-names>LT</given-names>
</name>
<etal/>
</person-group>. <article-title>The human tumor microbiome is composed of tumor type-specific intracellular bacteria</article-title>. <source>Science</source>. (<year>2020</year>) <volume>368</volume>:<page-range>973&#x2013;80</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1126/science.aay9189</pub-id>
</citation>
</ref>
<ref id="B6">
<label>6</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sepich-Poore</surname> <given-names>GD</given-names>
</name>
<name>
<surname>Zitvogel</surname> <given-names>L</given-names>
</name>
<name>
<surname>Straussman</surname> <given-names>R</given-names>
</name>
<name>
<surname>Hasty</surname> <given-names>J</given-names>
</name>
<name>
<surname>Wargo</surname> <given-names>JA</given-names>
</name>
<name>
<surname>Knight</surname> <given-names>R</given-names>
</name>
</person-group>. <article-title>The microbiome and human cancer</article-title>. <source>Science</source>. (<year>2021</year>) <volume>371</volume>:<elocation-id>eabc4552</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1126/science.abc4552</pub-id>
</citation>
</ref>
<ref id="B7">
<label>7</label>
<citation citation-type="journal">
<article-title>Tumors appear rife with bacterial lodgers</article-title>. <source>Cancer Discovery</source>. (<year>2020</year>) <volume>10</volume>:<page-range>1085&#x2013;6</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1158/2159-8290.CD-NB2020-062</pub-id>
</citation>
</ref>
<ref id="B8">
<label>8</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fu</surname> <given-names>A</given-names>
</name>
<name>
<surname>Yao</surname> <given-names>B</given-names>
</name>
<name>
<surname>Dong</surname> <given-names>T</given-names>
</name>
<name>
<surname>Cai</surname> <given-names>S</given-names>
</name>
</person-group>. <article-title>Emerging roles of intratumor microbiota in cancer metastasis</article-title>. <source>Trends Cell Biol</source>. (<year>2023</year>) <volume>33</volume>:<page-range>583&#x2013;93</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.tcb.2022.11.007</pub-id>
</citation>
</ref>
<ref id="B9">
<label>9</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname> <given-names>Z</given-names>
</name>
<name>
<surname>Hong</surname> <given-names>LL</given-names>
</name>
<name>
<surname>Ling</surname> <given-names>ZQ</given-names>
</name>
</person-group>. <article-title>Potential role of intratumor bacteria outside the gastrointestinal tract: More than passengers</article-title>. <source>Cancer Med</source>. (<year>2023</year>) <volume>12</volume>:<page-range>16756&#x2013;73</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/cam4.6298</pub-id>
</citation>
</ref>
<ref id="B10">
<label>10</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Poore</surname> <given-names>GD</given-names>
</name>
<name>
<surname>Kopylova</surname> <given-names>E</given-names>
</name>
<name>
<surname>Zhu</surname> <given-names>Q</given-names>
</name>
<name>
<surname>Carpenter</surname> <given-names>C</given-names>
</name>
<name>
<surname>Fraraccio</surname> <given-names>S</given-names>
</name>
<name>
<surname>Wandro</surname> <given-names>S</given-names>
</name>
<etal/>
</person-group>. <article-title>Microbiome analyses of blood and tissues suggest cancer diagnostic approach</article-title>. <source>Nature</source>. (<year>2020</year>) <volume>579</volume>:<page-range>567&#x2013;74</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41586-020-2095-1</pub-id>
</citation>
</ref>
<ref id="B11">
<label>11</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fuks</surname> <given-names>G</given-names>
</name>
<name>
<surname>Elgart</surname> <given-names>M</given-names>
</name>
<name>
<surname>Amir</surname> <given-names>A</given-names>
</name>
<name>
<surname>Zeisel</surname> <given-names>A</given-names>
</name>
<name>
<surname>Turnbaugh</surname> <given-names>PJ</given-names>
</name>
<name>
<surname>Soen</surname> <given-names>Y</given-names>
</name>
<etal/>
</person-group>. <article-title>Combining 16S rRNA gene variable regions enables high-resolution microbial community profiling</article-title>. <source>Microbiome</source>. (<year>2018</year>) <volume>6</volume>:<elocation-id>17</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/s40168-017-0396-x</pub-id>
</citation>
</ref>
<ref id="B12">
<label>12</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Segata</surname> <given-names>N</given-names>
</name>
<name>
<surname>Izard</surname> <given-names>J</given-names>
</name>
<name>
<surname>Waldron</surname> <given-names>L</given-names>
</name>
<name>
<surname>Gevers</surname> <given-names>D</given-names>
</name>
<name>
<surname>Miropolsky</surname> <given-names>L</given-names>
</name>
<name>
<surname>Garrett</surname> <given-names>WS</given-names>
</name>
<etal/>
</person-group>. <article-title>Metagenomic biomarker discovery and explanation</article-title>. <source>Genome Biol</source>. (<year>2011</year>) <volume>12</volume>:<fpage>R60</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/gb-2011-12-6-r60</pub-id>
</citation>
</ref>
<ref id="B13">
<label>13</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Douglas</surname> <given-names>GM</given-names>
</name>
<name>
<surname>Maffei</surname> <given-names>VJ</given-names>
</name>
<name>
<surname>Zaneveld</surname> <given-names>JR</given-names>
</name>
<name>
<surname>Yurgel</surname> <given-names>SN</given-names>
</name>
<name>
<surname>Brown</surname> <given-names>JR</given-names>
</name>
<name>
<surname>Taylor</surname> <given-names>CM</given-names>
</name>
<etal/>
</person-group>. <article-title>PICRUSt2 for prediction of metagenome functions</article-title>. <source>Nat Biotechnol</source>. (<year>2020</year>) <volume>38</volume>:<page-range>685&#x2013;8</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41587-020-0548-6</pub-id>
</citation>
</ref>
<ref id="B14">
<label>14</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Parks</surname> <given-names>DH</given-names>
</name>
<name>
<surname>Tyson</surname> <given-names>GW</given-names>
</name>
<name>
<surname>Hugenholtz</surname> <given-names>P</given-names>
</name>
<name>
<surname>Beiko</surname> <given-names>RG</given-names>
</name>
</person-group>. <article-title>STAMP: statistical analysis of taxonomic and functional profiles</article-title>. <source>Bioinformatics</source>. (<year>2014</year>) <volume>30</volume>:<page-range>3123&#x2013;4</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/bioinformatics/btu494</pub-id>
</citation>
</ref>
<ref id="B15">
<label>15</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Alpert</surname> <given-names>O</given-names>
</name>
<name>
<surname>Begun</surname> <given-names>L</given-names>
</name>
<name>
<surname>Issac</surname> <given-names>T</given-names>
</name>
<name>
<surname>Solhkhah</surname> <given-names>R</given-names>
</name>
</person-group>. <article-title>The brain-gut axis in gastrointestinal cancers</article-title>. <source>J Gastrointest Oncol</source>. (<year>2021</year>) <volume>12</volume>:<page-range>S301&#x2013;10</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.21037/jgo-2019-gi-04</pub-id>
</citation>
</ref>
<ref id="B16">
<label>16</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lloyd-Price</surname> <given-names>J</given-names>
</name>
<name>
<surname>Arze</surname> <given-names>C</given-names>
</name>
<name>
<surname>Ananthakrishnan</surname> <given-names>AN</given-names>
</name>
<name>
<surname>Schirmer</surname> <given-names>M</given-names>
</name>
<name>
<surname>Avila-Pacheco</surname> <given-names>J</given-names>
</name>
<name>
<surname>Poon</surname> <given-names>TW</given-names>
</name>
<etal/>
</person-group>. <article-title>Multi-omics of the gut microbial ecosystem in inflammatory bowel diseases</article-title>. <source>Nature</source>. (<year>2019</year>) <volume>569</volume>:<page-range>655&#x2013;62</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41586-019-1237-9</pub-id>
</citation>
</ref>
<ref id="B17">
<label>17</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kostic</surname> <given-names>AD</given-names>
</name>
<name>
<surname>Xavier</surname> <given-names>RJ</given-names>
</name>
<name>
<surname>Gevers</surname> <given-names>D</given-names>
</name>
</person-group>. <article-title>The microbiome in inflammatory bowel disease: current status and the future ahead</article-title>. <source>Gastroenterology</source>. (<year>2014</year>) <volume>146</volume>:<page-range>1489&#x2013;99</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1053/j.gastro.2014.02.009</pub-id>
</citation>
</ref>
<ref id="B18">
<label>18</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname> <given-names>YZ</given-names>
</name>
<name>
<surname>Li</surname> <given-names>YY</given-names>
</name>
</person-group>. <article-title>Inflammatory bowel disease: pathogenesis</article-title>. <source>World J Gastroenterol</source>. (<year>2014</year>) <volume>20</volume>:<page-range>91&#x2013;9</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.3748/wjg.v20.i1.91</pub-id>
</citation>
</ref>
<ref id="B19">
<label>19</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname> <given-names>R</given-names>
</name>
<name>
<surname>Zhou</surname> <given-names>R</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>H</given-names>
</name>
<name>
<surname>Li</surname> <given-names>W</given-names>
</name>
<name>
<surname>Pan</surname> <given-names>M</given-names>
</name>
<name>
<surname>Yao</surname> <given-names>X</given-names>
</name>
<etal/>
</person-group>. <article-title>Gut microbiota-stimulated cathepsin K secretion mediates TLR4-dependent M2 macrophage polarization and promotes tumor metastasis in colorectal cancer</article-title>. <source>Cell Death Differ</source>. (<year>2019</year>) <volume>26</volume>:<page-range>2447&#x2013;63</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41418-019-0312-y</pub-id>
</citation>
</ref>
<ref id="B20">
<label>20</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rooks</surname> <given-names>MG</given-names>
</name>
<name>
<surname>Garrett</surname> <given-names>WS</given-names>
</name>
</person-group>. <article-title>Gut microbiota, metabolites and host immunity</article-title>. <source>Nat Rev Immunol</source>. (<year>2016</year>) <volume>16</volume>:<page-range>341&#x2013;52</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/nri.2016.42</pub-id>
</citation>
</ref>
<ref id="B21">
<label>21</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dong</surname> <given-names>QL</given-names>
</name>
<name>
<surname>Xing</surname> <given-names>XY</given-names>
</name>
</person-group>. <article-title>Cancer cells arise from bacteria</article-title>. <source>Cancer Cell Int</source>. (<year>2018</year>) <volume>18</volume>:<fpage>205</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/s12935-018-0699-4</pub-id>
</citation>
</ref>
<ref id="B22">
<label>22</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Escoll</surname> <given-names>P</given-names>
</name>
<name>
<surname>Buchrieser</surname> <given-names>C</given-names>
</name>
</person-group>. <article-title>Metabolic reprogramming of host cells upon bacterial infection: Why shift to a Warburg-like metabolism</article-title>. <source>FEBS J</source>. (<year>2018</year>) <volume>285</volume>:<page-range>2146&#x2013;60</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/febs.14446</pub-id>
</citation>
</ref>
<ref id="B23">
<label>23</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Singh</surname> <given-names>R</given-names>
</name>
<name>
<surname>Kumar</surname> <given-names>CS</given-names>
</name>
<name>
<surname>Banerjee</surname> <given-names>M</given-names>
</name>
<name>
<surname>Gupta</surname> <given-names>S</given-names>
</name>
</person-group>. <article-title>A dual drug delivery platform for cancer-bacteria cotargeting</article-title>. <source>ACS Appl Bio Mater</source>. (<year>2019</year>) <volume>2</volume>:<page-range>5032&#x2013;41</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1021/acsabm.9b00724</pub-id>
</citation>
</ref>
<ref id="B24">
<label>24</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kalaora</surname> <given-names>S</given-names>
</name>
<name>
<surname>Nagler</surname> <given-names>A</given-names>
</name>
<name>
<surname>Nejman</surname> <given-names>D</given-names>
</name>
<name>
<surname>Alon</surname> <given-names>M</given-names>
</name>
<name>
<surname>Barbolin</surname> <given-names>C</given-names>
</name>
<name>
<surname>Barnea</surname> <given-names>E</given-names>
</name>
<etal/>
</person-group>. <article-title>Identification of bacteria-derived HLA-bound peptides in melanoma</article-title>. <source>Nature</source>. (<year>2021</year>) <volume>592</volume>:<page-range>138&#x2013;43</page-range>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41586-021-03368-8</pub-id>
</citation>
</ref>
<ref id="B25">
<label>25</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fu</surname> <given-names>A</given-names>
</name>
<name>
<surname>Yao</surname> <given-names>B</given-names>
</name>
<name>
<surname>Dong</surname> <given-names>T</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Yao</surname> <given-names>J</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>Y</given-names>
</name>
<etal/>
</person-group>. <article-title>Tumor-resident intracellular microbiota promotes metastatic colonization in breast cancer</article-title>. <source>Cell</source>. (<year>2022</year>) <volume>185</volume>:<fpage>1356</fpage>&#x2013;<lpage>1372.e26</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.cell.2022.02.027</pub-id>
</citation>
</ref>
<ref id="B26">
<label>26</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Heymann</surname> <given-names>C</given-names>
</name>
<name>
<surname>Bard</surname> <given-names>JM</given-names>
</name>
<name>
<surname>Heymann</surname> <given-names>MF</given-names>
</name>
<name>
<surname>Heymann</surname> <given-names>D</given-names>
</name>
<name>
<surname>Bobin-Dubigeon</surname> <given-names>C</given-names>
</name>
</person-group>. <article-title>The intratumoral microbiome: Characterization methods and functional impact</article-title>. <source>Cancer Lett</source>. (<year>2021</year>) <volume>522</volume>:<fpage>63</fpage>&#x2013;<lpage>79</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.canlet.2021.09.009</pub-id>
</citation>
</ref>
<ref id="B27">
<label>27</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gao</surname> <given-names>Z</given-names>
</name>
<name>
<surname>Guo</surname> <given-names>B</given-names>
</name>
<name>
<surname>Gao</surname> <given-names>R</given-names>
</name>
<name>
<surname>Zhu</surname> <given-names>Q</given-names>
</name>
<name>
<surname>Qin</surname> <given-names>H</given-names>
</name>
</person-group>. <article-title>Microbiota disbiosis is associated with colorectal cancer</article-title>. <source>Front Microbiol</source>. (<year>2015</year>) <volume>6</volume>:<elocation-id>20</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fmicb.2015.00020</pub-id>
</citation>
</ref>
<ref id="B28">
<label>28</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dohlman</surname> <given-names>AB</given-names>
</name>
<name>
<surname>Arguijo Mendoza</surname> <given-names>D</given-names>
</name>
<name>
<surname>Ding</surname> <given-names>S</given-names>
</name>
<name>
<surname>Gao</surname> <given-names>M</given-names>
</name>
<name>
<surname>Dressman</surname> <given-names>H</given-names>
</name>
<name>
<surname>Iliev</surname> <given-names>ID</given-names>
</name>
<etal/>
</person-group>. <article-title>The cancer microbiome atlas: a pan-cancer comparative analysis to distinguish tissue-resident microbiota from contaminants</article-title>. <source>Cell Host Microbe</source>. (<year>2021</year>) <volume>29</volume>:<fpage>281</fpage>&#x2013;<lpage>298.e5</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.chom.2020.12.001</pub-id>
</citation>
</ref>
<ref id="B29">
<label>29</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>J</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>Y</given-names>
</name>
<name>
<surname>Li</surname> <given-names>Z</given-names>
</name>
<name>
<surname>Gao</surname> <given-names>X</given-names>
</name>
<name>
<surname>Huang</surname> <given-names>D</given-names>
</name>
</person-group>. <article-title>Global analysis of microbiota signatures in four major types of gastrointestinal cancer</article-title>. <source>Front Oncol</source>. (<year>2021</year>) <volume>11</volume>:<elocation-id>685641</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fonc.2021.685641</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>