Gut mucus layer degradation is associated with aggressive cervical cancer phenotype

Gut mucus layer degradation is associated with aggressive cervical cancer phenotype | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Help Center Sign In Submit a Preprint Cite Share Download PDF Research article Gut mucus layer degradation is associated with aggressive cervical cancer phenotype Lauren E. Colbert This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-796525/v1 This work is licensed under a CC BY 4.0 License Status: Under Review Version 1 posted You are reading this latest preprint version Abstract Background Gut microbiome community composition differs between cervical cancer (CC) patients and healthy controls, and increased gut diversity is associated with improved outcomes after treatment. We proposed that functions of specific microbial species adjoining the mucus layer may directly impact the biology of CC. Results In this study, we examined metagenomes of rectal swabs in 41 CC patients using whole-genome shotgun sequencing and found a significant association between molecular functions encoded by the metagenomes with markers of aggressive cancer including initial tumor size and stage. Profiling of the molecular function abundances and their distributions identified 2 microbial communities co-existing in each metagenome but with distinct metabolism and taxonomic structures. Community A ( Clostridia and Proteobacteria predominant) was characterized by high activity of pathways involved in stress response, mucus glycan degradation and utilization of degradation byproducts. This community was prevalent in larger, advanced stage tumors. Conversely, community B ( Bacteroidia predominant) was characterized by fast growth, active oxidative phosphorylation, and production of vitamins. This community was prevalent in small, early-stage tumors. Conclusions Based on these results, we propose that increased mucus layer degradation is associated with a more aggressive cervical cancer phenotype. Metagenomics mucus layer gut cervical cancer tumor size bacteria Bacteroidales Clostridiales Figures Figure 1 Figure 2 Figure 3 Figure 4 Background Cervical Cancer (CC) is the fourth most common cancer in women worldwide [ 1 ]. The vast majority of these cases are caused by the human papillomavirus (HPV). Most individuals infected by the virus do not develop the cancer, as the viral infection is cleared before the virus integrates in the genome, expressing the cell cycle disrupting proteins E6 and E7[ 2 ]. Prior reports [ 3 ] have identified connections between the local (cervical/vaginal) microbiome and HPV infection and between cervical intraepithelial neoplasia and CC development. We previously reported differences between CC patients and healthy controls [ 4 ] and associations of gut diversity and composition with survival of CC patients based on 16Sv4 rRNA sequencing [ 5 ]. The effect of gut microbial functions on development and progression of CC in human is unknown. In a murine model of HPV cancer, we recently identified a significant association between gut and vaginal microbiomes and impact on the development and progression of CC [ 6 ]. Causal relationships are well established between the gut microbiome composition and the development of other cancers, such as colon cancer [ 7 ]. Additionally, a growing body of evidence links microorganisms to the efficacy of cancer therapies [ 8 – 11 ]. In most studies, bacterial DNA extracted from fecal samples is used as proxy for evaluation of the gut microbiome composition. However, the microbial communities adjoining and populating the outer mucus layer of the intestine are different from stool as a whole [ 12 ] and may be more relevant in modulating immune function. Indeed, composition of bacteria within the stool primarily reflect the dietary habits of patients, while microbial communities adjoining the mucus layer [ 13 , 14 ] also protect the gut epithelial surface that holds many important immune and metabolic functions. The difference may be more pronounced when the mucus layer is somehow disturbed as a result of a disease [ 15 ]. The mucus layer consists of an outer mucus layer, which is a habitat for commensal bacteria, and a smaller inner layer, which is attached to the epithelial cells and is lacking commensal bacteria [ 16 ]. This protective coatings are the crucial interface between the host and microorganisms [ 17 ] and are mainly (~ 89%) comprised of glycans (polysaccharides) attached by O-glycosylation to Muc2 (mucin 2) protein [ 18 ]. Galactose and N-acetylgalactosamine are major components of mucus glycan in normal human descending colon, although the structure varies in other parts of the intestine [ 19 ]. The mucus glycans provide attachment sites for bacteria [ 20 ] and select species that can influence the mucus layer structure and the host [ 21 ] by degrading mucin glycans, metabolizing products of the degradation, and producing metabolites that affect the host, positively or negatively. According to recent studies [ 22 , 23 ], mucins have also potent beneficial properties and can modulate microbial phenotypes suppressing quorum sensing, biofilm formation, and secretion of toxins. Thus, the microbial community adjoining and populating the outer mucus layer occupies a rather distinct environment within the intestine and can be affected not only by diet, but also by the mucus barrier and by the host. Previous bacterial 16S rRNA gene sequencing studies found that stool and rectal swab microbiotas from the same subject were similar and that rectal swabs can be used as a proxy for fecal samples [ 13 , 14 , 24 ]; However, in light of the findings referenced above, the mucus layer of the intestine, including rectum, can significantly modify phenotype of the adjoining microbial community. This means that even if the structure of the community is similar between stool and swabs, metabolic functions implemented by the same organisms in these different environments may be different. In addition, the use of Whole Genome Shotgun sequencing ( WGS ) instead of 16S RNA can provide finer resolution of the microbial community structure and functions and can reveal features not found by 16S. In this study we propose that molecular functions encoded by metagenomes, further referred as MA (Mucus Associated) because they populate the mucus layer or adjoin to it, are different among CC patients and likely associate with their clinical characteristics. We therefore explore metabolic characteristics of microbial communities populating the mucus layer in rectums of 41 CC patients. Rectal swabs were used to sample the communities and then to sequence them by WGS . Results Study cohort Patients were enrolled on an IRB approved (MDACC 2014 − 0543) study for longitudinal sampling of the gut microbiome (Table S1). For all 41 patients enrolled on study, median age was 49 (range 29 to 72 years), and median BMI was 28.6 (range 17.5 to 46.7). Most patients were Non-Hispanic White (44%) or Hispanic (44%) and had squamous carcinoma (74%) and stage II disease (54%). Median tumor size was 5.4cm (range 1.8cm to 11.5cm). More information on clinicopathological characteristics of the patients is provided in Table S2. Functional richness of MA metagenome associates with larger tumor size and with advanced CC stage We first explored associations of total number of predicted genes and functions with patient and tumor characteristics (Table S3, Fig. 1 a). A higher numbers of unique functions was positively correlated with increased tumor size (R = 0.41, P = 0.008) (Fig. 1 b). Importantly, a large number of unique KEGG orthologous groups (KOs) was significantly associated not only with large tumor size (P = 0.02), but also with advanced stages (III/IV) of cervical cancer (P = 0.005; Fig. 1 a). Unsupervised hierarchical clustering reveals microbial communities with 2 distinct profiles of molecular functions Next, we explored whether specific genes or functions were associated with patient and tumor characteristics using unsupervised hierarchical clustering of abundances. We selected 2,396 KOs that were found in greater than 30 samples. These KOs clustered into 2 large clusters and 4 sub-clusters (Fig. 1 c). In general, Cluster 2 was comprised of smaller tumors (“ST cluster”), when compared with Cluster 1 (“LT cluster”; t-test P = 0.06). Within the LT cluster, Sub-cluster 1.1 was dominated by younger patients (t-test P = 0.05) with large tumors (t-test P = 0.03). The KEGG pathway analysis of the ST cluster and the LT cluster using KEGG mapper [ 38 ] also identified a set of pathways enriched in each of the clusters (Fig. 1 c). The ST Cluster was enriched with pathways associated with rapid microbial cell proliferation, including ribosome biogenesis (Additional File 1: Figure S3a), DNA repair (Additional File 1: Figure S3b), oxidative phosphorylation (Additional File 1: Figure S3c), and synthesis of peptidoglycan and lipopolysaccharides ( Fig. 1 d). Almost all enzymes involved in synthesis of biotin (vitamin B7) were also found in the ST cluster ( Fig. 1 d). The LT cluster (Cluster 1 in Fig. 1 c) was enriched with pathways associated with bacterial stress response, as indicated by activation of quorum sensing (Additional File 1: Figure S4a), sporulation, and degradation of glycan. The glycan degradation pathways included high activity of the KEGG pathway, “Other glycan degradation” (Fig. 1 e). In addition, pathways involved in utilization of glycan degradation byproducts, such as fructose and mannose (Additional File 1: Figure S5a), and especially galactose (Additional File 1: Figure S5b) were enriched. Galactose is the major constituent (~ 85%) of glycans in normal human gastric mucus [ 44 ] and was also enriched in the LT cluster. Both branches of the pentose phosphate pathway, the oxidative branch maintaining redox balance in stress conditions and non-oxidative branch that supply glycolysis with intermediates derived from pentoses [ 45 ], were also more active in the LT cluster (Additional File 1: Figure S6). Most enzymes involved in utilization of ethanolamine (Fig. 1 e), a breakdown product of human and bacterial cell membrane [ 46 ], and enzymes involved in production of ornithine (Fig. 1 e), a precursor for synthesis of polyamines, such as putrescine [ 47 , 48 ], were found enriched in the LT cluster. The biological processes of flagellar assembly and chemotaxis (Additional File 1: Figure S7) were also up-regulated in the LT cluster, revealing the importance of mobility for bacteria residing in the metagenomes [ 49 ]. Supervised comparison of metagenomes associated with largest and smallest tumors To confirm the association of large tumor size with glycan degrading microbial communities, we divided patient samples into “Large Tumor” (LT-group) and “Small Tumor” (ST-group) groups. Two groups of metagenomes were selected to quantify abundances of the communities. The LT-group included 14 metagenomes of patients with the largest tumors, and the ST-group had 14 metagenomes of patients with the smallest tumors (Fig. 2 a). The LT-group were also significantly more enriched in patients with stage III/IV of CC then ST-group (chi-squared contingency table tests p-values 0.02), since tumor size is roughly correlated with stage in the 2009 FIGO staging system used in the study. The diversity and richness of molecular functions was significantly higher in LT versus ST-group (Fig. 2 b and Fig. 2 c respectively). There were 1,496 differentially abundant KOs between the groups with most of the KOs (75%) were more abundant in the LT-group. The density plots of the differentially abundant KOs were also different between the groups. The KO abundances were significantly more scattered around the mean (standard deviation is 2.6) in LT-group when compared with ST-group (standard deviation is 1.3). The KO abundances mean value was almost twice as low in LT-group than in ST-group (Fig. 2 d and Fig. 2 e respectively). The observed density plots indicated that 2 distinct populations (sets) of molecular functions, referred to as LT-abundant and ST-abundant, might coincide in each metagenome but be most representative in either LT- or ST-group of samples. This conclusion was confirmed by density plots of all KOs found in LT-group (Fig. 2 f) and in ST-group (Fig. 2 g). By fitting parameters of the normal distributions from Fig. 2 d and Fig. 2 e to density plots in Fig. 2 f and 2 g using multiple linear regression, we found that although both, LT- and ST- abundant functions were present in LT- and ST-group of metagenomes, the fitting parameter of LT-abundant population was increased in LT-group by 1.2 and the fitting parameter of ST-abundant population was decreased by 1.5 in the group. CAZyme Enzyme families associated with small tumors To associate the LT-abundant KO population with the glycan degradation we compared abundances of carbohydrate-active enzymes (CAZyme) families in this population and in the ST- abundant KO population using CAZy database [39]. The database classifies CAZymes according to their functions, such as synthesis of complex carbohydrates or their hydrolysis. Glycoside hydrolases GH, which are enzymes involved in hydrolysis of glycosidic bonds between carbohydrates or between a carbohydrate and a non-carbohydrate moiety, were highly enriched in LT abundant KO population (Fig. 2h). There were 16 GH among KOs in the population, and none in the ST abundant population, which included only 3 Glycosyl Transferases (GT); the enzymes that are mainly involved in biosynthesis of disaccharides, oligosaccharides, and polysaccharides. We further compared ST- and LT-abundant KO populations in terms of abundances of glycosidases, enzymes annotated by EC 3.2.1- involved in hydrolyzes of O- and S-glycosyl compounds. There were 32 glycosidases in 28 samples, and many of them were significantly more abundant in metagenomes of patients with large tumors (Fig 2i). Only one enzyme, lysozyme, was found to be significantly more abundant in the small tumors group of patients. Lysozyme is a known component of two-component cell lysis cassette in bacteriophages [50]. Biological processes associated with small and large tumors The KEGG pathway analysis of differentially abundant KOs was further used to identify specific biological processes represented by LT-abundant and ST-abundant KOs (Fig. 2j). The representative pathways revealed by the analysis were consistent with results of the supervised analysis of 41 samples. Very active proliferation in ST-group of metagenomes were indicated by a high pathway enrichment score (ST- versus LT-abundant population) for DNA replication, ribosome, t-RNA biosynthesis, homologous recombination and mismatch repair. Stress response and degradation of the mucus layer in LT-group were indicated by enrichment of biofilm formation pathway, galactose metabolism (Fig. 2j, Additional File 1: Figure S8), and by a significant increase in the abundance of trehalose-specific PTS components (Fig. 2k) among LT- versus ST-abundant KOs. Ornithine biosynthesis (Additional File 1: Figure S9), putrescine and spermidine metabolisms (Additional File 1: Figure S8; Fig. 2k) were also enriched in LT- abundant KOs, while enzymes involved in vitamin B12 production were enriched among ST-abundant KOs (Fig2 k). Association of KEGG pathway with clinicopathological characteristics and recurrence free survival (RFS) None of the studied clinical characteristics showed significant (p<0.05) association with RFS. This analysis was limited by the low number of events in the study cohort. Among KEGG pathways enriched in either LT- or ST-group of MA metagenomes, only a trend of negative association with RFS was found for high activity of glycan degradation pathway and of Phosphotransferase system (Fig. 3e and Fig. 3f). Both pathways were enriched in LT-group of MA metagenomes. Vice versa, significant positive association with RFS was found for high activity of Ribosome biogenesis and DNA repair pathways enriched in ST- group of MA metagenomes (Fig. 3g and 3h). Associations with RFS were also observed for stage and positive nodal status (Additional File 1: Figure S10A and S10B). These positive associations are well documented in the literature as very significant in studies of large number of CC patients (>10,000) [51]. In this study, patients with positive nodes also had significantly increased activity of glycan degradation pathway (p=0.005) and ornithine biosynthesis (p=0.05) in MA metagenomes (Fig.3a and Fig. 3b). Activity of the pathways was also significantly higher in patients with stage III/IV tumors (Fig. 3c and Fig. 3d), while activity of ribosome biogenesis was higher in metagenomes of stage I/II patients (Additional File 1: Figure S11). Taxonomic associations with tumor size Difference in taxonomic structure of MA metagenomes in LT-group versus ST-group can be seen at different taxonomic levels (Fig. 4a, Additional File 1: Figure S12-S14) with dominance of Phylum Firmicutes ( Class Clostridia and Order Clostridiales) and Phylum Proteobacteria in LT- group and Phylum Bacteroidetes (Class Bacteroidia, Families Prevotellaceae and Ruminococcaceae in ST-group. Like the differentially abundant molecular functions, the density plots of differentially abundant species were also different between the groups (Fig. 4b) with a greater median abundance of species in ST-group. No difference, however, was seen between LT- and ST-groups of metagenomes in terms of species richness, diversity, and evenness (Fig. 4c). Taxonomic annotation of most abundant contigs (OTUs) in each group and comparison of their abundances between groups revealed that many of the contigs belonged to Order Bacteroidales, and they were significantly more abundant in ST-group. There were 5 putative species of Genus Porphyromonas (P. endodontalis, P. sp COT_239OH1446, P. bennonis, P. sp HMSC065F10 and P. uenonis ) and 2 species of Genus Prevotella ( P. Timonensis and P. Buccalis ) among Bacteroidales (Fig. 4d, Additional File 2: Table S3). None of the genera was significantly abundant in LT-group. Contigs annotated by Class Clostiridia and other levels of taxonomy in the Class were abundant in L-group including 12 putative species of Family Ruminococcaceae. Only 1 species of Ruminococcaceae ( AF41_9 ) was found more abundant in ST group. The observations were consistent with results obtained by comparison of LT- and ST-groups using LEfSe (Fig. 4e and Additional File 1: Figure S15). The latter analysis has also identified less abundant taxa enriched in one of the groups, such as Class Tissierellia, Genera Ezakiella, Murdochiella and Hungatella, in ST-group, as well as Families Ruminococcaceae, Lachnospieraceae, Enterobacteriaceae , and Order Enterobacterales in LT group (Additional File 1: Figure S16). Taxonomic associations with pathways To identify taxa utilizing the glycan degradation pathway in each metagenome we selected contigs that encode enzymes of the pathway and quantified abundances of taxa representing these contigs. Overlay of the enzyme abundances across metagenomes with the relative abundances of taxa involved in glycan degradation (Fig. 4f) showed that major taxa degrading glycan are different between MA metagenomes of patients with small and large tumors. Clostridiales is more likely to encode enzymes of the glycan degradation in the LT-group, while Bacteroidales is the major taxon that encode the enzymes in ST-group. The result is consistent with the difference in taxonomic structure of MA metagenomes described in the previous paragraph. Further correlation analysis revealed a significant positive association of the glycan degradation pathway score with Bacteroidales, Clostridiales, and with non-classified taxa. A significant negative correlation was found with Actinomycetales, Chlamydiales, and Tissierellales (Additional File 1: Figure S17). In general, there was significant variation across taxa involved in glycan degradation. The variation was especially dramatic in MA metagenomes of patients with large tumors. Discussion Herein, we have described significant differences in the gut metagenomes adjoining the mucus layer (MA metagenome) between cervical cancer patients with small, early stage tumors and large, advanced stage tumors at the time of treatment. Specifically, MA metagenomes of patients with small tumors were enriched in molecular functions associated with biosynthesis and rapid proliferation. By contrast, MA metagenomes of patients with large tumors were enriched with functions associated with bacterial stress response and degradation of glycan. In addition, the utilization of ethanolamine and functions related to bacterial cells mobility were more abundant among young patients with large tumors and advanced stage (Fig. 1 c), suggesting particularly aggressive tumor biology. Activity of glycan degradation and of ornithine biosynthesis were also significantly associated with the advance tumor stage and the presence of cancer cells in lymph nodes, while activity of the ribosome biogenesis significantly associated with low stage tumors and better recurrence-free survival. This dominance may lead to significant changes in the metabolic environment of the intestine that directly impact tumor growth and progression. It is possible, this aggressive tumor biology is fueled by differing molecular functions in MA metagenome. We have described two distinct microbial consortia, which we refer to here as proliferating and mucus degrading (Fig. 4 g). Both consortia co-exist; although patients with larger tumor size and advanced stage show an obvious dominance of the mucus degrading consortia over proliferating. The tumor promoting effects of glycan degradation and utilization of major degradation products, such as galactose and other pentoses [ 19 ], may help drive this aggressive tumor growth. Glycan degradation leads to production of metabolites with known tumor growth promoting effects and with resistance to radio- and chemotherapy, such as ornithine [ 52 ] and ceramides [ 53 , 54 ]. The metabolites are produced as by-products of glycan degradation (Fig. 1 e) and can be further metabolized either by bacteria or by host and give oncogenic effect after their metabolization. Ceramide itself, for example, is a powerful tumor suppressor, but products of its metabolism are potent tumor survival factors associated with resistance to cancer therapies [ 53 ]. Ornithine can be used by the consortium or by the host to synthesize the polyamines putrescine and spermidine (Additional File 1: Figure S8), which cause tumorigenic transformation and tumor progression [ 55 ]. If the gut metabolites cross the gut mucosa into blood through the mucus layer, they may promote systemic proliferation of tumor cells. In addition, the intensive ethanolamine utilization (Fig. 1 e) by the LT-group metagenomes, especially pronounced in younger patients (Fig. 1 c), also suggests a more pathogenic environment in the intestine. The set of enzymes involved in the pathway are similar to the eut operon in Salmonella enterica serovar Typhimurium [ 56 ], a known gastrointestinal pathogen. Many other species that contain the eut genes are also pathogens [ 46 ], because ethanolamine is derived from the membrane phospholipid phosphatidylethanolamine, an important component of all bacterial and eukaryotic cells. The intensive utilization of ethanolamine may indicate a degradation of the colonic epithelium [ 46 ] that secretes peptides inhibiting bacterial penetration into the inner colonic mucus layer and blocking bacterial mobility [ 17 , 57 ]. Indeed, the activity of pathways associated with mobility, such as chemotaxis and flagella assembly, are also enhanced in MA metagenomes of the LT-group, particularly in sub-cluster of younger patients (Fig. 1 c). Conversely, the increased synthesis of peptidoglycan and lipopolysaccharides by the proliferating consortia may improve the intestinal microenvironment and the immune response in patients with smaller tumors [ 58 ]. The increased production of vitamins B7 (biotin) and B12 (cobalamin) may also be important for immune cell function and reduction of cellular oxidative stress [ 59 – 61 ]. Overall, these metabolites may have a tumor suppressive effect. The findings suggest a potential therapeutic intervention in CC by shifting the balance in MA metagenome from mucus degrading consortium to proliferating one for suppression of the tumor growth and enhancement responses to treatment. This study does not make clear what comes first, the biological mechanism of mucus degradation or tumor progression. It is likely that both processes are tightly related and influence each other. Cervical cancer typically develops over years. Differences in the microbial flora of a large tumor may reflect changes in the tumor microenvironment directly affected by local tumor progression, or gut metagenome functions may allow increased tumor growth via systemic factors released into the bloodstream. Alternatively, tumor size may be a surrogate for other aggressive tumor biology, such as hypoxia or necrosis, which increase as tumors grow. Although this study associates large tumor size with degradation of mucus layer and reveal by-products of the degradation with known tumor promoting effects, it doesn’t provide direct evidence for an association of tumor-promoting metabolites with tumor size. It is possible that tumor progression and associated changes of organismal processes may be primary and drive changes in the mucus layer microenvironment and in molecular functions of the MA community through production of cytokines and metabolites [ 17 ] or by changing the mucin glycans that are also important host signals selecting microorganisms and making them less pathogenic [ 23 ]. Indeed, the study reveals stress-associated changes in the mucus layer microenvironment that are indicated by reduced proliferation, and by increased activity of quorum sensing and sporulation in mucus degrading consortium. Importantly, taxonomic structure differed between groups only when functions and pathways were considered, suggesting taxonomic comparison alone is not adequate. Further studies using mouse or cell culture models may be necessary to decipher biological mechanisms underlying the discovered association between CC and mucus layer degradation. In addition, because potential mechanisms underlying the association of the tumor size and mucus layer degradation are not specific to CC, it is very likely that similar associations between the tumor size and structure and function of MA metagenomes may be seen in other cancer types. Shifting the balance between proliferating and mucus degrading consortia in MA metagenomes may directly affect cancer therapy and drug toxicity. The mucus layer plays a very important protective role in the intestine [ 62 ], therefore its unbalanced degradation can impair sensitivity to or ability to resist toxic effects of drugs. Recently, a growing body of evidence linked microorganisms to cancer therapy efficacy and toxicity [ 8 – 10 ]. Our previous study of metagenomes in fecal samples of melanoma patients with different response to anti-PD-1 therapy [ 10 ] reveals that degradation associated pathways are enriched in non-responders, while biosynthetic pathways are enriched in responders. We were not able to associate catabolic pathways with the mucus layer degradation in this study; however, we can speculate that the degradation can be responsible for the enrichment of catabolic processes in non-responding patients. The small number of samples and the use of fecal samples instead of rectal swabs might complicate identification of the mucus layer degradation in the study. It is likely that fecal samples are not optimal for the evaluation because they represent a different environment in the intestine. Further studies using WGS are necessary to compare fecal and swab samples from the same patient when mucus layer degradation is evaluated on association of patient’s response to cancer therapies, cancer stage, or patient survival. Conclusions Overall, these findings suggest potential interventions related to mucus enhancement that could improve outcomes and decrease cervical cancer progression. Further mechanistic studies are important to explore if the disturbed balance between proliferating and mucus degrading consortia can be corrected by diet, probiotics, antibiotics, or other interventions. Methods Patient samples collection and processing Forty-one patients treated at the University of Texas, M.D. Anderson Cancer Center or the Lyndon B. Johnson Clinic (LBJ) at Harris Health with a diagnosis of CC participated in the study. The patients were enrolled on an Institutional Review Board (IRB) approved prospective protocol. Informed consent was obtained to collect rectal swabs. All patient samples were acquired prior to receiving any treatment. Samples were collected from each out of 41 patients by a clinician performing rectal exams using a matrix-designed quick-release Isohelix swab. The swabs were initially stored in 20ul of proteinase K and 400 μl of lysis buffer (Isohelix) within 1 h of sample collection and then were frozen and kept at -80°C. WGS sequencing and metagenome assembly Whole Genome Shotgun sequencing was performed on genomic bacterial DNA (gDNA), which was extracted to maximize bacterial DNA yield from specimens while keeping background amplification to a minimum [25, 26]. Libraries were constructed from each sample using the KAPA Hyper Prep Kit (Kapa Biosystems, Wilmington, MA, USA) and sequenced using the Illumina HiSeqX platform with the 2 x 150 bp paired-end read protocol. Sequencing reads were derived from raw BCL files which were retrieved from the sequencer and called into fastqs by Casava v1.8.3 (Illumina). The appropriate read preparation steps, such as quality control, trimming and filtering, and host DNA removal prior to further analysis, were performed using an in-house pipeline (Additional file 1: Figure S1b; Additional file 2: Table S1). Briefly, paired-end raw sequence reads were filtered and trimmed using BBMap [27]. The trimmed reads were mapped to a hg38 reference database (GCA_000001405.28) using bowtie2 [28] to remove host contamination. The cleaned reads were then assembled to longer sequences (contigs) using both MEGAHIT [29] and metaSPAdes [30]. The assembled contigs were filtered, to remove those that were smaller than 1000bp, and binned using MetaBAT2 [31]. The assembled and binned contigs were used for gene predictions by Prodigal [32]. Annotation of the genes by KEGG ortholog groups (KOs) was implemented by KofamKOALA [33], and taxonomic classification of the contigs was done by CAT and BAT [34]. The read coverage of each assembled contig was calculated by aligning the cleaned reads directly to the contig using BBMap and by counting the mapped read by featureCounts [35]. The read coverage, GC content, taxonomic and functional annotations for each gene/contig were summarized using a Perl script. All the software tools were running with default parameters if not specified. Versions and sources of the software tools or packages used in the pipeline are listed in Additional file 2: Table S1. Output of the assembly pipeline was a set of assembled contigs for each sample, their taxonomic annotation and read coverage, gene predictions for each contig, and functional annotation of each gene by KEGG Orthologous Group (KO) if found. Computational analysis of assembled genomes Annotations of assembled contigs for all samples were aggregated into 1 table, referred to as the Metagenome Function Abundance (MFA) table (Additional file 1: Figure S1b). Each column in the table represents a metagenome, and each row represents a predicted known molecular function annotated by KO. Thus, each cell in the MFA table has a quantity of the specified (KO id) molecular function in the specified sample (Sample id). Quantification of the molecular functions by MFA table is explained using a toy example provided in Additional file 1: Figure S2. The MFA table was normalized using total number of reads in each sample and then multiplied by 1,000,000. Further analysis of the normalized MFA table included unsupervised and supervised methods, annotation by biological processes and pathways, and integration with clinicopathological characteristics of the patients (Additional file 1: Figure S1c). Unsupervised analysis of MFA table The MFA table was filtered to select KOs found in at least 30 patients, log2-transformed and centered using the median, and then hierarchically clustered by the open source clustering software [36] with default parameters. The inferred clusters of samples were tested for association with clinical information including age, CC stage and tumor size, using fisher.test() and wilcox.test() functions in R. The inferred clusters of KOs were searched for overlap with known KEGG referenced pathways and modules using the “Search&Color Pathway” tool in KEGG mapper [37, 38] with further manual curation of the results. Supervised analysis of MFA table The analysis was used to find differentially abundant KOs in MA metagenomes of patients with large versus small tumors. Samples for the analysis were selected by sorting all 41 samples by tumor size of the CC patients and assigning top 14 samples with largest tumors to LT-group and bottom 14 samples with smallest tumors to ST-group for the comparison. All KOs found in at least 1 sample were included in the analysis. Fisher’s test and Mann Whitney tests were used to find differentially abundant KOs between the group with p-value < 0.05 (either test) without adjustment. The KO was considered enriched in LT-group if Fisher’s test p-values 0.05 but the Wilcoxon test p-value <0.05 than the enrichment was inferred by difference in mean abundances between ST- and LT-groups. The differentially abundant KOs were searched for overlap with known KEGG referenced pathways and modules as described above. The overlapping KOs were used to infer the pathways enrichment score calculated as ratio of difference between percentage of KOs overlapped with the pathway in ST- and LT-groups to sum of the percentages. Only top scored pathways that include 9 and more KOs for LT-group and 4 and more KOs for ST-group are considered. Annotation of the differentially abundant KOs by carbohydrate-active enzymes (CAZy) families was implemented using the mapping table between the KO ID and CAZy family ID [39]. The table was downloaded from KEGG on 26 October 2019. The differentially abundant KOs were further used to quantify relative abundances of microbial subpopulations expressing them in ST and in LT group of samples. The decomposition of cell populations, referred to as ST-dominant and LT-dominant, was implemented based on characteristics of density plots of KO The differentially abundant KOs in ST and LT group were further used to construct 2 normal curves for ST-abundant and LT-abundant KOs. Parameters of the curves, mean and standard deviation, were assumed to characterize molecular functions expressed by the microbial community dominating in either ST-group or LT-group, although each of the groups is comprised of both communities. To quantify abundances of the communities in each group, we fitted the 2 normal curves to density plots of all KOs found in the two groups using multiple linear regression. The obtained fitting coefficients were used to quantify the dominance of each set in metagenomes of each group. Survival analysis The analysis was used to associate the Recurrence Free Survival (RFS) probability with clinicopathological characteristics and KEGG pathway [40] . Activity of each pathway found to be associated with the tumor size was quantified for each patient using the mean value of log2-transformed normalized abundances of KOs involved in the pathway. The mean value is referred to as the pathway activity score. To generate Kaplan-Meier plots for the pathway activity score, we categorized the score based on a cut-off set at the first quartile. Observations for each variable falling within the first quarter were labeled as “Low” and those greater than the first quartile cut-off were labelled as “High”. The analysis was implemented for each pathway identified as differentially enriched in either small or large tumor group metagenomes. We used WHO standards to set cut-offs to categorize BMI into “Underweight/Normal weight” vs. “Overweight/Obese”. We used a cut-off set at the median for other continuous clinicopathological characteristics, namely, age and tumor size. Cox proportional hazards regression analysis [41] was used to evaluate predictive value of KEGG pathways and clinical characteristics on RFS time. The R packages ‘survival’ and ‘survminer’ were used to compute the survival curves, and to visualize them as Kaplan-Meier plots. Taxonomic characterization of the cohort All predicted contigs annotated by Length (L), Depths (D), and by taxonomy were used to quantify abundances of the unique set of species in 14 samples of LT group and in 14 samples of ST group. For each sample, depth of the contig was multiplied by its length and the values were summarized for all contigs that belong to the same taxa to quantify the taxa abundance. The values were normalized using the sum of the values for all contigs identified in the sample. The OTU (operating taxonomic units) table was created by merging the abundances of all identified taxa for all 28 samples. Two approaches were used to find differentially abundant putative taxa between LT and ST group. In the first approach, the OTUs were selected by Fisher’s exact test (for rare OTUs) and Mann Whitney tests (for common OTUs) with p-value cutoff < 0.05 (either test) without adjustment. In the other approach, the Linear discriminant analysis (LDA) Effect Size (LEfSe) [42] available as a Galaxy [43] module was used to determine the taxa that are differentially abundant between LT and ST groups. The analysis was run with default parameters, except the threshold for the logarithmic LDA score for discriminative taxa. The threshold was set to 2.5 instead of 2. To reveal taxa involved in glycan degradation pathway in each metagenome we have selected all contigs that encode enzymes involved in the pathway. The set of unique taxa from taxonomic annotation of the contigs was considered as the representative taxa of the pathway. To create an OTU table of the representative taxa, each OTU were quantified by multiplying the length and the depth of each selected contig and by summarizing the obtained values for all contigs encoding the OTU. The 100% stacked area plot was used to visualize taxonomic structure of the community involved in the pathway across all studied samples. The Pearson correlation coefficients were used to evaluate association of the log2-transformed score of the KEGG pathway with the abundance of each identified taxon. Declarations Ethics approval and consent to participate Informed written consent was obtained from eligible patients willing to participate in our study. This study was under a protocol approved by the UT M.D. Anderson Cancer Center Institutional Review Board (MDACC 2014–0543) Consent for publication Not applicable. Availability of data and materials The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request. Raw sequencing reads generated by WGS are available in the BioProject database under accession PRJNA702617(https://www.ncbi.nlm.nih.gov/bioproject/PRJNA702617). Competing interests The authors declare that they have no competing interests. Funding This research was supported in part by the Radiological Society of North America Resident/Fellow Award (to L.E.C.); by the National Institutes of Health (NIH) through MD Anderson’s Cancer Center Support Grant P30CA016672 and through T32 grant CA101642-14 (T.T.S.); and by The University of Texas MD Anderson Cancer Center HPV-related Cancers Moonshot (L.E.C., A. K.). The funding sources were not involved in the development of the research hypothesis, study design, data analysis, or manuscript writing. Authors’ Contributions LEC, AHK, and TVK designed the study and wrote the paper with contribution from XW, JZ, AF, PJE, AJ, LLL, KMS, NA, RRJ, MBEA, and LR. TS, TTS, KYC, EL, MAK, GB, JY, JFP, MPM, and PO performed the experiment and sequencing. XW, JZ, TVK, XS and HS developed the computation pipeline for metagenome assembly and annotation. TVK, XW, and MBEA processed the data. All authors read and approved the final manuscript. Acknowledgements We sincerely thank the patients who volunteered to participate in this study. We acknowledge the Moonshot program at the University of Texas MD Anderson Cancer Center for financial support of the study and High-Performance Computing at the University of Texas MD Anderson Cancer Center for providing computational resources that have contributed to the research results reported in this paper. References Vu M, Yu J, Awolude OA, Chuang L. Cervical cancer worldwide. Curr Probl Cancer. 2018;42:457–65. Szymonowicz KA, Chen JJ. Biological and clinical aspects of HPV-related cancers. Cancer Biology Medicine. 2020;17:864–78. Curty G, de Carvalho PS, Soares MA. The Role of the Cervicovaginal Microbiome on the Genesis and as a Biomarker of Premalignant Cervical Intraepithelial Neoplasia and Invasive Cervical Cancer. International Journal of Molecular Sciences 2020, 21. Sims TT, Colbert LE, Zheng J, Delgado Medrano AY, Hoffman KL, Ramondetta L, Jazaeri A, Jhingran A, Schmeler KM, Daniel CR, Klopp A. Gut microbial diversity and genus-level differences identified in cervical cancer patients versus healthy controls. Gynecol Oncol. 2019;155:237–44. Sims TT, El Alam MB, Karpinets TV, Dorta-Estremera S, Hegde VL, Nookala S, Yoshida-Court K, Wu XG, Biegert GWG, Medrano AYD, et al: Gut microbiome diversity is an independent predictor of survival in cervical cancer patients receiving chemoradiation. Communications Biology 2021, 4. Karpinets TV, Solley TN, Mikkelson MD, Dorta-Estremera S, Nookala SS, Medrano AYD, Petrosino JF, Mezzari MP, Zhang J, Futreal PA, et al. Effect of Antibiotics on Gut and Vaginal Microbiomes Associated with Cervical Cancer Development in Mice. Cancer Prev Res (Phila). 2020;13:997–1006. Alhinai EA, Walton GE, Commane DM. The Role of the Gut Microbiota in Colorectal Cancer Causation. Int J Mol Sci 2019, 20. Zitvogel L, Galluzzi L, Viaud S, Vetizou M, Daillere R, Merad M, Kroemer G. Cancer and the gut microbiota: an unexpected link. Sci Transl Med. 2015;7:271ps271. Alexander JL, Wilson ID, Teare J, Marchesi JR, Nicholson JK, Kinross JM. Gut microbiota modulation of chemotherapy efficacy and toxicity. Nat Rev Gastroenterol Hepatol. 2017;14:356–65. Gopalakrishnan V, Spencer CN, Nezi L, Reuben A, Andrews MC, Karpinets TV, Prieto PA, Vicente D, Hoffman K, Wei SC, et al. Gut microbiome modulates response to anti-PD-1 immunotherapy in melanoma patients. Science. 2018;359:97–103. Proenca JT, Barral DC, Gordo I. Commensal-to-pathogen transition: One-single transposon insertion results in two pathoadaptive traits in Escherichia coli -macrophage interaction. Sci Rep. 2017;7:4504. Lu HP, Lai YC, Huang SW, Chen HC, Hsieh CH, Yu HT. Spatial heterogeneity of gut microbiota reveals multiple bacterial communities with distinct characteristics. Scientific Reports 2014, 4. Bassis CM, Moore NM, Lolans K, Seekatz AM, Weinstein RA, Young VB, Hayden MK, Program CDCPE. Comparison of stool versus rectal swab samples and storage conditions on bacterial community profiles. BMC Microbiol. 2017;17:78. Reyman M, van Houten MA, Arp K, Sanders EAM, Bogaert D. Rectal swabs are a reliable proxy for faecal samples in infant gut microbiota research based on 16S-rRNA sequencing. Sci Rep. 2019;9:16072. Mentella MC, Scaldaferri F, Pizzoferrato M, Gasbarrini A, Miggiano GAD. Nutrition, IBD and Gut Microbiota: A Review. Nutrients 2020, 12. Johansson ME, Phillipson M, Petersson J, Velcich A, Holm L, Hansson GC. The inner of the two Muc2 mucin-dependent mucus layers in colon is devoid of bacteria. Proc Natl Acad Sci U S A. 2008;105:15064–9. Schroeder BO. Fight them or feed them: how the intestinal mucus layer manages the gut microbiota. Gastroenterology Report. 2019;7:3–12. Holmen Larsson JM, Thomsson KA, Rodriguez-Pineiro AM, Karlsson H, Hansson GC. Studies of mucus in mouse stomach, small intestine, and colon. III. Gastrointestinal Muc5ac and Muc2 mucin O-glycan patterns reveal a regiospecific distribution. Am J Physiol Gastrointest Liver Physiol. 2013;305:G357–63. Capon C, Maes E, Michalski JC, Leffler H, Kim YS. Sd(a)-antigen-like structures carried on core 3 are prominent features of glycans from the mucin of normal human descending colon. Biochem J. 2001;358:657–64. Juge N. Microbial adhesins to gastrointestinal mucus. Trends Microbiol. 2012;20:30–9. Cai R, Cheng C, Chen J, Xu X, Ding C, Gu B. Interactions of commensal and pathogenic microorganisms with the mucus layer in the colon. Gut Microbes. 2020;11:680–90. Kavanaugh NL, Zhang AQ, Nobile CJ, Johnson AD, Ribbeck K. Mucins suppress virulence traits of Candida albicans. mBio. 2014;5:e01911. Wheeler KM, Carcamo-Oyarce G, Turner BS, Dellos-Nolan S, Co JY, Lehoux S, Cummings RD, Wozniak DJ, Ribbeck K. Mucin glycans attenuate the virulence of Pseudomonas aeruginosa in infection. Nat Microbiol. 2019;4:2146–54. Short MI, Hudson R, Besasie BD, Reveles KR, Shah DP, Nicholson S, Johnson-Pais TL, Weldon K, Lai Z, Leach RJ, et al. Comparison of rectal swab, glove tip, and participant-collected stool techniques for gut microbiome sampling. BMC Microbiol. 2021;21:26. Human Microbiome Project C. A framework for human microbiome research. Nature. 2012;486:215–21. Human Microbiome Project C. Structure, function and diversity of the healthy human microbiome. Nature. 2012;486:207–14. 2020 BhsnpbfAM. Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9. Li D, Liu CM, Luo R, Sadakane K, Lam TW. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics. 2015;31:1674–6. Nurk S, Meleshko D, Korobeynikov A, Pevzner PA. metaSPAdes: a new versatile metagenomic assembler. Genome Res. 2017;27:824–34. Kang DD, Li F, Kirton E, Thomas A, Egan R, An H, Wang Z. MetaBAT 2: an adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies. PeerJ. 2019;7:e7359. Hyatt D, Chen GL, Locascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics. 2010;11:119. Aramaki T, Blanc-Mathieu R, Endo H, Ohkubo K, Kanehisa M, Goto S, Ogata H. KofamKOALA: KEGG Ortholog assignment based on profile HMM and adaptive score threshold. Bioinformatics. 2020;36:2251–2. von Meijenfeldt FAB, Arkhipova K, Cambuy DD, Coutinho FH, Dutilh BE. Robust taxonomic classification of uncharted microbial sequences and bins with CAT and BAT. Genome Biol. 2019;20:217. Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30:923–30. de Hoon MJ, Imoto S, Nolan J, Miyano S. Open source clustering software. Bioinformatics. 2004;20:1453–4. Kanehisa M, Sato Y, Kawashima M, Furumichi M, Tanabe M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 2016;44:D457–62. Kanehisa M, Sato Y. KEGG Mapper for inferring cellular functions from protein sequences. Protein Sci. 2020;29:28–35. Lombard V, Ramulu HG, Drula E, Coutinho PM, Henrissat B. The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 2014;42:D490–5. Li JCA. Modeling survival data: Extending the Cox model. Sociological Methods Research. 2003;32:117–20. Cox DR. Regression Models and Life-Tables. Journal of the Royal Statistical Society Series B-Statistical Methodology. 1972;34:187-+. Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, Huttenhower C. Metagenomic biomarker discovery and explanation. Genome Biology 2011, 12. Afgan E, Baker D, Batut B, van den Beek M, Bouvier D, Cech M, Chilton J, Clements D, Coraor N, Gruning BA, et al. The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2018 update. Nucleic Acids Res. 2018;46:W537–44. Clamp JR, Cooper B, Creeth JM, Ene D, Barrett J, Gough M. The Presence of Polysaccharide in Normal Human Gastric Mucus. Biochem J. 1983;215:421–3. Stincone A, Prigione A, Cramer T, Wamelink MMC, Campbell K, Cheung E, Olin-Sandoval V, Gruning NM, Kruger A, Alam MT, et al. The return of metabolism: biochemistry and physiology of the pentose phosphate pathway. Biol Rev. 2015;90:927–63. Kaval KG, Garsin DA. Ethanolamine Utilization in Bacteria. mBio 2018, 9. Pegg AE. Mammalian Polyamine Metabolism and Function. Iubmb Life. 2009;61:880–94. Vazquez A, Kamphorst JJ, Markert E, Schug ZT, Tardito S, Gottlieb E. Cancer metabolism at a glance. J Cell Sci. 2016;129:3367–73. Matilla MA, Krell T. The effect of bacterial chemotaxis on host infection and pathogenicity. Fems Microbiology Reviews. 2018;42:40–67. Srividhya KV, Krishnaswamy S. Subclassification and targeted characterization of prophage-encoded two-component cell lysis cassette. J Biosci. 2007;32:979–90. McComas KN, Torgeson AM, Ager BJ, Hellekson C, Burt LM, Maurer KA, Werner TL, Gaffney DK. The variable impact of positive lymph nodes in cervical cancer: Implications of the new FIGO staging system. Gynecol Oncol. 2020;156:85–92. Geck RC, Foley JR, Murray Stewart T, Asara JM, Casero RA Jr, Toker A. Inhibition of the polyamine synthesis enzyme ornithine decarboxylase sensitizes triple-negative breast cancer cells to cytotoxic chemotherapy. J Biol Chem. 2020;295:6263–77. Yura Y, Masui A, Hamada M. Inhibitors of Ceramide- and Sphingosine-Metabolizing Enzymes as Sensitizers in Radiotherapy and Chemotherapy for Head and Neck Squamous Cell Carcinoma. Cancers 2020, 12. Morad SAF, Cabot MC. Ceramide-orchestrated signalling in cancer cells. Nat Rev Cancer. 2013;13:51–65. Murray-Stewart TR, Woster PM, Casero RA. Targeting polyamine metabolism for cancer therapy and prevention. Biochem J. 2016;473:2937–53. Roof DM, Roth JR. Ethanolamine utilization in Salmonella typhimurium. J Bacteriol. 1988;170:3855–63. Bergstrom JH, Birchenough GMH, Katona G, Schroeder BO, Schutte A, Ermund A, Johansson MEV, Hansson GC. Gram-positive bacteria are held at a distance in the colon mucus by the lectin-like protein ZG16. Proc Natl Acad Sci USA. 2016;113:13833–8. Tang C, Ding R, Sun J, Liu J, Kan J, Jin C. The impacts of natural polysaccharides on intestinal microbiota and immune responses - a review. Food Funct. 2019;10:2290–312. Kung JT, Mackenzie CG, Talmage DW. The requirement for biotin and fatty acids in the cytotoxic T-cell response. Cell Immunol. 1979;48:100–10. Peterson CT, Rodionov DA, Osterman AL, Peterson SN. B Vitamins and Their Role in Immune Regulation and Cancer. Nutrients 2020, 12. Agrawal S, Agrawal A, Said HM. Biotin deficiency enhances the inflammatory response of human dendritic cells. American Journal of Physiology-Cell Physiology. 2016;311:C386–91. Camilleri M. Leaky gut: mechanisms, measurement and clinical implications in humans. Gut. 2019;68:1516–26. Supplementary Files AdditionalFile1SupplementaryFigures.pdf Supplemental Figures [See Additional File1_Supplementary Figures] AdditionalFile2SupplementaryTables.xlsx Supplemental Tables [See Additional File2_Supplementary Tables] Cite Share Download PDF Status: Under Review Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-796525","acceptedTermsAndConditions":true,"allowDirectSubmit":false,"archivedVersions":[],"articleType":"Research article","associatedPublications":[],"authors":[{"id":44698184,"identity":"d17608f2-61fb-4de5-b398-6fa09bea2368","order_by":0,"name":"Lauren E. Colbert","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA9klEQVRIiWNgGAWjYBACA/YDjIdBDH4gm4GBjRgtPAkMIC0Skg1Ea5FQgGgxOEC8FuYDhwsq7tUZXzu88XNBmU2+fHvvAYaPe2rxaGFLODzjTLGE2e20YukZ59IsN5w5l8A449lxPFp4DA7ztiUAteQYSPO2HTYwkMgxYOY5cIyAln8JEsazc4x/87b9N5CfQVCLDlBLQ4KEgXSOGdCWAwYMN8BaanBr4UkD+uVYguSM22ll1jznkg0MzpwxODjjwAGcWuzbDx98XFCTwM8/O3nzbZ4yOwP59h7DBx8O1OHUgh0ArThMohYgINWWUTAKRsEoGMYAADY9VM4weVi7AAAAAElFTkSuQmCC","orcid":"","institution":"","correspondingAuthor":true,"submittingAuthor":false,"prefix":"","firstName":"Lauren","middleName":"E.","lastName":"Colbert","suffix":""}],"badges":[],"createdAt":"2021-08-09 10:20:41","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-796525/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-796525/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":16876853,"identity":"268c4675-e442-4c60-8d23-f84c2759834a","added_by":"auto","created_at":"2021-12-30 16:52:44","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":705081,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eAssociation of molecular functions of MA metagenomes with KEGG pathways and clinical characteristics of CC patients. a \u003c/strong\u003eLarge number of unique functions is associated with larger tumor size and advance CC stage. \u003cstrong\u003eb \u003c/strong\u003eLinear relationship between the number of\u003cstrong\u003e \u003c/strong\u003eunique functions and the tumor size. \u003cstrong\u003ec \u003c/strong\u003eUnsupervised hierarchical clustering of metagenomes in terms of abundances of molecular functions. The heat map of the clustering shows abundances of 2396 KOs as rows and 41 MA metagenomes as columns. Boxes with enriched KEGG pathways are colored by red if the clustered molecular functions are more abundant in metagenomes of patients with Large tumors (Cluster 1) and in green if the molecular functions are more abundant in metagenomes of patients with Small tumors (Cluster 2). \u003cstrong\u003ed \u003c/strong\u003eSynthesis of peptidoglycan, lipopolysaccharides and biotin are enriched among KOs that are abundant in Cluster 2 (Small tumor) metagenomes. \u003cstrong\u003ee\u003c/strong\u003e Glycan degradation, ethanolamine utilization and production of ornithine are enriched among KOs that are abundant in Cluster 1 (Large tumor) metagenomes. Enzymes involved in the pathways (c and b) have red borders if they are found among KOs of the cluster.\u003c/p\u003e","description":"","filename":"ScreenShot20211230at11.22.09AM.png","url":"https://assets-eu.researchsquare.com/files/rs-796525/v1/e99e66e16308a0090a28f82b.png"},{"id":16876983,"identity":"68f3859d-7675-4f23-9a68-8447710ad811","added_by":"auto","created_at":"2021-12-30 16:55:44","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":773751,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eSupervised comparison of Largest tumors (LT-group) and Smallest tumors (ST-group) of MA metagenomes.\u003c/strong\u003e \u003cstrong\u003ea\u003c/strong\u003e Grouping of samples into ST- and LT-groups for supervised comparison. \u003cstrong\u003eb,c\u003c/strong\u003e Greater richness and diversity of molecular functions in LT-group of metagenomes. \u003cstrong\u003ed,e \u003c/strong\u003eDensity plots of\u003cstrong\u003e \u003c/strong\u003edifferentially abundant KOs in LT- (c) and ST- (d) groups of metagenomes. Parameters of the normal distributions were used to find amplitudes of LT- and ST-abundant molecular functions among all KOs in each tumor group. \u003cstrong\u003ef,g \u003c/strong\u003eDensity plots of all KOs identified in LT- (e) and ST- (f) groups of metagenomes. \u003cstrong\u003eh \u003c/strong\u003eCarbohydrate Active enZymes (CAZymes) identified in LT- and ST- groups. \u003cstrong\u003ei\u003c/strong\u003e Glycoside hydrolases involved in degradation of glycans (EC 3.2.1.-) are differentially enriched in LT- versus ST- group of metagenomes. \u003cstrong\u003ej\u003c/strong\u003e KEGG pathways enriched among KOs differentially abundant between ST- and LT- groups of metagenomes. k Some KEGG pathways differentially abundant between the groups. Phosphotransferase and putresine utilization pathway are more abundant in LT-group, while vitamin B12 production is more abundant in ST-group.\u0026nbsp;\u003c/p\u003e","description":"","filename":"ScreenShot20211230at11.22.30AM.png","url":"https://assets-eu.researchsquare.com/files/rs-796525/v1/332ac43b01d183ca8c1950cb.png"},{"id":16876850,"identity":"3929b237-04de-4502-aa82-bfc0a6a96548","added_by":"auto","created_at":"2021-12-30 16:52:44","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":327411,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eAssociation of KEGG pathway and clinicopathological characteristics. a \u003c/strong\u003eIncreased activity of Glycan degradation pathway and \u003cstrong\u003eb\u003c/strong\u003e Ornithine biosynthesis in metagenomes of patients with positive nodes on imaging. Categorization of patients into the negative (red color) and positive (green color) groups according to absence or presence of cancer cells in lymph nodes respectively revealed that 6 out of 12 patients with negative nodes had very low activity of Glycan degradation pathway. \u003cstrong\u003ec \u003c/strong\u003eIncreased activity of Glycan degradation pathway and \u003cstrong\u003ed\u003c/strong\u003e Ornithine biosynthesis in metagenomes of patients with stage III/IV cervical cancer. \u003cstrong\u003ee, f \u003c/strong\u003eNegative association of Phosphotransferase system and of \u003cstrong\u003e\u0026nbsp;\u003c/strong\u003eGlycan degradation pathway with recurrence free survival (RFS). Both pathways were enriched in metagenomes of patients with large tumors. The pathway activity was quantified in each metagenome by the mean value of log2-transformed normalized abundances of KOs involved in the pathway; the value is referred to as the pathway activity score, which was categorized as high in the first quartile and low in the rest. \u003cstrong\u003eg, h\u003c/strong\u003e Significant positive association of Ribosome biogenesis and of DNA repair pathways with RRS. Both pathways were enriched in metagenomes of patients with small tumors.\u003c/p\u003e","description":"","filename":"ScreenShot20211230at11.22.43AM.png","url":"https://assets-eu.researchsquare.com/files/rs-796525/v1/78b99c391f5b49c0dedf540f.png"},{"id":16876982,"identity":"1faa8b3c-de76-4785-85a6-d8f5008713dc","added_by":"auto","created_at":"2021-12-30 16:55:44","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":709580,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eTaxonomic profiling and biological models of MA metagenomes in LT- and ST-groups. a\u003c/strong\u003e Taxonomic structure of metagenomes at the \u003cem\u003ePhyla\u003c/em\u003e level. \u003cstrong\u003eb\u003c/strong\u003e Density plots of differentially abundant species in LT- and ST-groups.\u003cstrong\u003e \u003c/strong\u003eFisher’s test and Mann Whitney tests were used to find differentially abundant OTUs with p-value \u0026lt; 0.05 (either test) without adjustment. \u003cstrong\u003ec \u003c/strong\u003eRichness, diversity, and evenness of OTUs in LT and ST-groups.\u003cstrong\u003e d \u003c/strong\u003eTaxonomic annotation of most abundant contigs in LT- and ST-groups. \u003cstrong\u003ee\u003c/strong\u003e Differentially abundant taxa between ST- and LT-\u003cstrong\u003e \u003c/strong\u003egroups identified by LEfSe. \u003cstrong\u003ef \u003c/strong\u003eCladogram of differentially abundant taxa. \u003cstrong\u003e\u0026nbsp;g \u003c/strong\u003eBiological\u003cstrong\u003e \u003c/strong\u003emodels of Mucus Adjoining (MA) bacterial consortium in LT- and ST-groups.\u003c/p\u003e","description":"","filename":"ScreenShot20211230at11.23.00AM.png","url":"https://assets-eu.researchsquare.com/files/rs-796525/v1/8bed58f3e8397c939bf5c659.png"},{"id":16877013,"identity":"8e0412b3-275d-4112-8889-b7254f469451","added_by":"auto","created_at":"2021-12-30 16:55:47","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":525528,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-796525/v1/385a14d1-e584-45ec-ad67-355dc5d9e346.pdf"},{"id":16876855,"identity":"49bb467e-0ee6-4f7a-a455-d19f91aeed02","added_by":"auto","created_at":"2021-12-30 16:52:44","extension":"pdf","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":4431550,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eSupplemental Figures [See Additional File1_Supplementary Figures]\u003c/strong\u003e\u003c/p\u003e","description":"","filename":"AdditionalFile1SupplementaryFigures.pdf","url":"https://assets-eu.researchsquare.com/files/rs-796525/v1/db92839388ab3c75efc39164.pdf"},{"id":16876851,"identity":"62ffdf79-fc3b-4808-83e4-4404ea4b8a08","added_by":"auto","created_at":"2021-12-30 16:52:44","extension":"xlsx","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":27988,"visible":true,"origin":"","legend":"\u003cp\u003e\u003cstrong\u003eSupplemental Tables [See Additional File2_Supplementary Tables]\u003c/strong\u003e\u003c/p\u003e","description":"","filename":"AdditionalFile2SupplementaryTables.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-796525/v1/8f97ff2b48eb3e233b49e9ea.xlsx"}],"financialInterests":"","formattedTitle":"Gut mucus layer degradation is associated with aggressive cervical cancer phenotype","fulltext":[{"header":"Background","content":"\u003cp\u003eCervical Cancer (CC) is the fourth most common cancer in women worldwide [\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e]. The vast majority of these cases are caused by the human papillomavirus (HPV). Most individuals infected by the virus do not develop the cancer, as the viral infection is cleared before the virus integrates in the genome, expressing the cell cycle disrupting proteins E6 and E7[\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e]. Prior reports [\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e] have identified connections between the local (cervical/vaginal) microbiome and HPV infection and between cervical intraepithelial neoplasia and CC development. We previously reported differences between CC patients and healthy controls [\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e] and associations of gut diversity and composition with survival of CC patients based on 16Sv4 rRNA sequencing [\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e]. The effect of gut microbial functions on development and progression of CC in human is unknown. In a murine model of HPV cancer, we recently identified a significant association between gut and vaginal microbiomes and impact on the development and progression of CC [\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e]. Causal relationships are well established between the gut microbiome composition and the development of other cancers, such as colon cancer [\u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e]. Additionally, a growing body of evidence links microorganisms to the efficacy of cancer therapies [\u003cspan additionalcitationids=\"CR9 CR10\" citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eIn most studies, bacterial DNA extracted from fecal samples is used as proxy for evaluation of the gut microbiome composition. However, the microbial communities adjoining and populating the outer mucus layer of the intestine are different from stool as a whole [\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e] and may be more relevant in modulating immune function. Indeed, composition of bacteria within the stool primarily reflect the dietary habits of patients, while microbial communities adjoining the mucus layer [\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e, \u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e] also protect the gut epithelial surface that holds many important immune and metabolic functions. The difference may be more pronounced when the mucus layer is somehow disturbed as a result of a disease [\u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eThe mucus layer consists of an outer mucus layer, which is a habitat for commensal bacteria, and a smaller inner layer, which is attached to the epithelial cells and is lacking commensal bacteria [\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e]. This protective coatings are the crucial interface between the host and microorganisms [\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e] and are mainly (~\u0026thinsp;89%) comprised of glycans (polysaccharides) attached by O-glycosylation to Muc2 (mucin 2) protein [\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e]. Galactose and N-acetylgalactosamine are major components of mucus glycan in normal human descending colon, although the structure varies in other parts of the intestine [\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e]. The mucus glycans provide attachment sites for bacteria [\u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e20\u003c/span\u003e] and select species that can influence the mucus layer structure and the host [\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e] by degrading mucin glycans, metabolizing products of the degradation, and producing metabolites that affect the host, positively or negatively. According to recent studies [\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e, \u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e], mucins have also potent beneficial properties and can modulate microbial phenotypes suppressing quorum sensing, biofilm formation, and secretion of toxins. Thus, the microbial community adjoining and populating the outer mucus layer occupies a rather distinct environment within the intestine and can be affected not only by diet, but also by the mucus barrier and by the host.\u003c/p\u003e \u003cp\u003ePrevious bacterial 16S rRNA gene sequencing studies found that stool and rectal swab microbiotas from the same subject were similar and that rectal swabs can be used as a proxy for fecal samples [\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e, \u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e, \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e]; However, in light of the findings referenced above, the mucus layer of the intestine, including rectum, can significantly modify phenotype of the adjoining microbial community. This means that even if the structure of the community is similar between stool and swabs, metabolic functions implemented by the same organisms in these different environments may be different. In addition, the use of Whole Genome Shotgun sequencing (\u003cem\u003eWGS\u003c/em\u003e) instead of 16S RNA can provide finer resolution of the microbial community structure and functions and can reveal features not found by 16S.\u003c/p\u003e \u003cp\u003eIn this study we propose that molecular functions encoded by metagenomes, further referred as MA (Mucus Associated) because they populate the mucus layer or adjoin to it, are different among CC patients and likely associate with their clinical characteristics. We therefore explore metabolic characteristics of microbial communities populating the mucus layer in rectums of 41 CC patients. Rectal swabs were used to sample the communities and then to sequence them by \u003cem\u003eWGS\u003c/em\u003e.\u003c/p\u003e"},{"header":"Results","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003eStudy cohort\u003c/h2\u003e \u003cp\u003ePatients were enrolled on an IRB approved (MDACC 2014\u0026thinsp;\u0026minus;\u0026thinsp;0543) study for longitudinal sampling of the gut microbiome (Table S1). For all 41 patients enrolled on study, median age was 49 (range 29 to 72 years), and median BMI was 28.6 (range 17.5 to 46.7). Most patients were Non-Hispanic White (44%) or Hispanic (44%) and had squamous carcinoma (74%) and stage II disease (54%). Median tumor size was 5.4cm (range 1.8cm to 11.5cm). More information on clinicopathological characteristics of the patients is provided in Table S2.\u003c/p\u003e \u003cp\u003e \u003cb\u003eFunctional richness of MA metagenome associates with larger tumor size and with advanced CC stage\u003c/b\u003e \u003c/p\u003e \u003cp\u003eWe first explored associations of total number of predicted genes and functions with patient and tumor characteristics (Table S3, Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ea). A higher numbers of unique functions was positively correlated with increased tumor size (R\u0026thinsp;=\u0026thinsp;0.41, P\u0026thinsp;=\u0026thinsp;0.008) (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eb). Importantly, a large number of unique KEGG orthologous groups (KOs) was significantly associated not only with large tumor size (P\u0026thinsp;=\u0026thinsp;0.02), but also with advanced stages (III/IV) of cervical cancer (P\u0026thinsp;=\u0026thinsp;0.005; Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ea).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003e \u003cb\u003eUnsupervised hierarchical clustering reveals microbial communities with 2 distinct profiles of molecular functions\u003c/b\u003e \u003c/p\u003e \u003cp\u003eNext, we explored whether specific genes or functions were associated with patient and tumor characteristics using unsupervised hierarchical clustering of abundances. We selected 2,396 KOs that were found in greater than 30 samples. These KOs clustered into 2 large clusters and 4 sub-clusters (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ec). In general, Cluster 2 was comprised of smaller tumors (\u0026ldquo;ST cluster\u0026rdquo;), when compared with Cluster 1 (\u0026ldquo;LT cluster\u0026rdquo;; t-test P\u0026thinsp;=\u0026thinsp;0.06). Within the LT cluster, Sub-cluster 1.1 was dominated by younger patients (t-test P\u0026thinsp;=\u0026thinsp;0.05) with large tumors (t-test P\u0026thinsp;=\u0026thinsp;0.03).\u003c/p\u003e \u003cp\u003eThe KEGG pathway analysis of the ST cluster and the LT cluster using KEGG mapper [\u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e38\u003c/span\u003e] also identified a set of pathways enriched in each of the clusters (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ec). The ST Cluster was enriched with pathways associated with rapid microbial cell proliferation, including ribosome biogenesis (Additional File 1: Figure S3a), DNA repair (Additional File 1: Figure S3b), oxidative phosphorylation (Additional File 1: Figure S3c), and synthesis of peptidoglycan and lipopolysaccharides \u003cb\u003e(\u003c/b\u003eFig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ed). Almost all enzymes involved in synthesis of biotin (vitamin B7) were also found in the ST cluster \u003cb\u003e(\u003c/b\u003eFig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ed).\u003c/p\u003e \u003cp\u003eThe LT cluster (Cluster 1 in Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ec) was enriched with pathways associated with bacterial stress response, as indicated by activation of quorum sensing (Additional File 1: Figure S4a), sporulation, and degradation of glycan. The glycan degradation pathways included high activity of the KEGG pathway, \u0026ldquo;Other glycan degradation\u0026rdquo; (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ee). In addition, pathways involved in utilization of glycan degradation byproducts, such as fructose and mannose (Additional File 1: Figure S5a), and especially galactose (Additional File 1: Figure S5b) were enriched. Galactose is the major constituent (~\u0026thinsp;85%) of glycans in normal human gastric mucus [\u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e] and was also enriched in the LT cluster. Both branches of the pentose phosphate pathway, the oxidative branch maintaining redox balance in stress conditions and non-oxidative branch that supply glycolysis with intermediates derived from pentoses [\u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e45\u003c/span\u003e], were also more active in the LT cluster (Additional File 1: Figure S6). Most enzymes involved in utilization of ethanolamine (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ee), a breakdown product of human and bacterial cell membrane [\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e], and enzymes involved in production of ornithine (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ee), a precursor for synthesis of polyamines, such as putrescine [\u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e47\u003c/span\u003e, \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e48\u003c/span\u003e], were found enriched in the LT cluster. The biological processes of flagellar assembly and chemotaxis (Additional File 1: Figure S7) were also up-regulated in the LT cluster, revealing the importance of mobility for bacteria residing in the metagenomes [\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e].\u003c/p\u003e \u003cp\u003e \u003cb\u003eSupervised comparison of metagenomes associated with largest and smallest tumors\u003c/b\u003e \u003c/p\u003e \u003cp\u003eTo confirm the association of large tumor size with glycan degrading microbial communities, we divided patient samples into \u0026ldquo;Large Tumor\u0026rdquo; (LT-group) and \u0026ldquo;Small Tumor\u0026rdquo; (ST-group) groups. Two groups of metagenomes were selected to quantify abundances of the communities. The LT-group included 14 metagenomes of patients with the largest tumors, and the ST-group had 14 metagenomes of patients with the smallest tumors (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ea). The LT-group were also significantly more enriched in patients with stage III/IV of CC then ST-group (chi-squared contingency table tests p-values 0.02), since tumor size is roughly correlated with stage in the 2009 FIGO staging system used in the study. The diversity and richness of molecular functions was significantly higher in LT versus ST-group (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eb and Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ec respectively). There were 1,496 differentially abundant KOs between the groups with most of the KOs (75%) were more abundant in the LT-group. The density plots of the differentially abundant KOs were also different between the groups. The KO abundances were significantly more scattered around the mean (standard deviation is 2.6) in LT-group when compared with ST-group (standard deviation is 1.3). The KO abundances mean value was almost twice as low in LT-group than in ST-group (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ed and Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ee respectively). The observed density plots indicated that 2 distinct populations (sets) of molecular functions, referred to as LT-abundant and ST-abundant, might coincide in each metagenome but be most representative in either LT- or ST-group of samples. This conclusion was confirmed by density plots of all KOs found in LT-group (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ef) and in ST-group (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eg). By fitting parameters of the normal distributions from Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ed and Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ee to density plots in Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ef and \u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eg using multiple linear regression, we found that although both, LT- and ST- abundant functions were present in LT- and ST-group of metagenomes, the fitting parameter of LT-abundant population was increased in LT-group by 1.2 and the fitting parameter of ST-abundant population was decreased by 1.5 in the group.\u003c/p\u003e\u003cp\u003e\u003cstrong\u003eCAZyme Enzyme families associated with small tumors\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; To associate the LT-abundant KO population with the glycan degradation we compared abundances of carbohydrate-active enzymes (CAZyme) families in this population and in the ST- abundant KO population using CAZy database [39]. \u0026nbsp;The database classifies CAZymes according to their functions, such as synthesis of complex carbohydrates or their hydrolysis.\u003c/p\u003e\n\u003cp\u003e\u0026nbsp;Glycoside hydrolases GH, which are enzymes involved in hydrolysis of glycosidic bonds between carbohydrates or between a carbohydrate and a non-carbohydrate moiety, were highly enriched in LT abundant KO population (Fig. 2h). There were 16 GH among KOs in the population, and none in the ST abundant population, which included only 3 Glycosyl Transferases (GT); the enzymes that are mainly involved in biosynthesis of disaccharides, oligosaccharides, and polysaccharides. We further compared ST- and LT-abundant KO populations in terms of abundances of glycosidases, enzymes annotated by EC 3.2.1- involved in hydrolyzes of O- and S-glycosyl compounds. \u0026nbsp;There were 32 glycosidases in 28 samples, and many of them were significantly more abundant in metagenomes of patients with large tumors (Fig 2i). Only one enzyme, lysozyme, was found to be significantly more abundant in the small tumors group of patients. Lysozyme is a known component of two-component cell lysis cassette in bacteriophages [50].\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eBiological processes associated with small and large tumors\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe KEGG pathway analysis of differentially abundant KOs was further used to identify specific biological processes represented by LT-abundant and ST-abundant KOs (Fig. 2j). The representative pathways revealed by the analysis were consistent with results of the supervised analysis of 41 samples. Very active proliferation in ST-group of metagenomes were indicated by a high pathway enrichment score (ST- versus LT-abundant population) for DNA replication, ribosome, t-RNA biosynthesis, homologous recombination and mismatch repair. Stress response and degradation of the mucus layer in LT-group were indicated by enrichment of biofilm formation pathway, galactose metabolism (Fig. 2j, Additional File 1: Figure S8), and by a significant increase in the abundance of trehalose-specific PTS components (Fig. 2k) among LT- versus ST-abundant KOs. Ornithine biosynthesis (Additional File 1: Figure S9), putrescine and spermidine metabolisms (Additional File 1: Figure S8; Fig. 2k) were also enriched in LT- abundant KOs, while enzymes involved in vitamin B12 production were enriched among ST-abundant KOs (Fig2 k).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAssociation of KEGG pathway with clinicopathological characteristics and recurrence free survival (RFS)\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNone of the studied clinical characteristics showed significant (p\u0026lt;0.05) association with RFS. This analysis was limited by the low number of events in the study cohort. \u0026nbsp;Among KEGG pathways enriched in either LT- or ST-group of MA metagenomes, only a trend of negative association with RFS was found for high activity of glycan degradation pathway and of Phosphotransferase system (Fig. 3e and Fig. 3f). Both pathways were enriched in LT-group of MA metagenomes. Vice versa, significant positive association with RFS was found for high activity of Ribosome biogenesis and DNA repair pathways enriched in ST- group of MA metagenomes (Fig. 3g and 3h).\u003cstrong\u003e\u0026nbsp;\u003c/strong\u003eAssociations with RFS were also observed for stage and positive nodal status (Additional File 1: Figure S10A and S10B). These positive associations are well documented in the literature as very significant in studies of large number of CC patients (\u0026gt;10,000) [51]. In this study, patients with positive nodes also had significantly increased activity of glycan degradation pathway (p=0.005) and ornithine biosynthesis (p=0.05) in MA metagenomes (Fig.3a and Fig. 3b). Activity of the pathways was also significantly higher in patients with stage III/IV tumors (Fig. 3c and Fig. 3d), while activity of ribosome biogenesis was higher in metagenomes of stage I/II patients (Additional File 1: Figure S11).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eTaxonomic associations with tumor size\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eDifference in taxonomic structure of MA metagenomes in LT-group versus ST-group can be seen at different taxonomic levels (Fig. 4a, Additional File 1: Figure S12-S14) with dominance of Phylum \u003cem\u003eFirmicutes (\u003c/em\u003eClass\u003cem\u003e\u0026nbsp;Clostridia\u0026nbsp;\u003c/em\u003eand\u003cem\u003e\u0026nbsp;\u003c/em\u003eOrder\u003cem\u003e\u0026nbsp;Clostridiales)\u003c/em\u003e and Phylum \u003cem\u003eProteobacteria\u003c/em\u003e in LT- group and Phylum \u003cem\u003eBacteroidetes\u003c/em\u003e (Class\u003cem\u003e\u0026nbsp;Bacteroidia,\u0026nbsp;\u003c/em\u003eFamilies\u003cem\u003e\u0026nbsp;Prevotellaceae\u0026nbsp;\u003c/em\u003eand \u003cem\u003eRuminococcaceae\u003c/em\u003e in ST-group. Like the differentially abundant molecular functions, the density plots of differentially abundant species were also different between the groups (Fig. 4b) with a greater median abundance of species in ST-group. No difference, however, was seen between LT- and ST-groups of metagenomes in terms of species richness, diversity, and evenness (Fig. 4c).\u003cstrong\u003e\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eTaxonomic annotation of most abundant contigs (OTUs) in each group and comparison of their abundances between groups revealed that many of the contigs belonged to Order\u003cem\u003e\u0026nbsp;Bacteroidales,\u003c/em\u003e and\u003cem\u003e\u0026nbsp;\u003c/em\u003ethey\u003cem\u003e\u0026nbsp;\u003c/em\u003ewere significantly more abundant in ST-group. There were 5 putative species of Genus\u003cem\u003e\u0026nbsp;Porphyromonas (P.\u003c/em\u003e \u003cem\u003eendodontalis, P. sp COT_239OH1446, P. bennonis, P. sp HMSC065F10 and P. uenonis\u003c/em\u003e) and 2 species of Genus\u003cem\u003e\u0026nbsp;Prevotella\u003c/em\u003e (\u003cem\u003eP. Timonensis and P. Buccalis\u003c/em\u003e) among \u003cem\u003eBacteroidales\u003c/em\u003e (Fig. 4d, Additional File 2: Table S3). None of the genera was significantly abundant in LT-group. Contigs annotated by Class\u003cem\u003e\u0026nbsp;Clostiridia\u003c/em\u003e and other levels of taxonomy in the \u003cem\u003eClass\u003c/em\u003e were abundant in L-group including 12 putative species of Family \u003cem\u003eRuminococcaceae.\u0026nbsp;\u003c/em\u003eOnly 1 species of \u003cem\u003eRuminococcaceae\u003c/em\u003e (\u003cem\u003eAF41_9\u003c/em\u003e) was found more abundant in ST group. The observations\u0026nbsp;were consistent with results obtained by comparison of LT- and ST-groups using LEfSe (Fig. 4e and\u0026nbsp;Additional File 1: Figure S15). The latter analysis has also identified less abundant taxa enriched in one of the groups, such as Class\u003cem\u003e\u0026nbsp;Tissierellia,\u0026nbsp;\u003c/em\u003eGenera\u003cem\u003e\u0026nbsp;Ezakiella, Murdochiella and Hungatella,\u003c/em\u003e in ST-group, as well as Families\u003cem\u003e\u0026nbsp;Ruminococcaceae, Lachnospieraceae, Enterobacteriaceae\u003c/em\u003e, and Order\u003cem\u003e\u0026nbsp;Enterobacterales\u003c/em\u003e in LT group (Additional File 1: Figure S16).\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eTaxonomic associations with pathways\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eTo identify taxa utilizing the glycan degradation pathway in each metagenome we selected contigs that encode enzymes of the pathway and quantified abundances of taxa representing these contigs. Overlay of the enzyme abundances across metagenomes with the relative abundances of taxa involved in glycan degradation (Fig. 4f) showed that major taxa degrading glycan are different between MA metagenomes of patients with small and large tumors. \u003cem\u003eClostridiales\u003c/em\u003e is more likely to encode enzymes of the glycan degradation in the LT-group, while \u003cem\u003eBacteroidales\u0026nbsp;\u003c/em\u003eis the major\u003cem\u003e\u0026nbsp;\u003c/em\u003etaxon that encode the enzymes in ST-group. The result is consistent with the\u0026nbsp;difference in taxonomic structure of MA metagenomes described in the previous paragraph. Further correlation analysis revealed a significant positive association of the glycan degradation pathway score with \u003cem\u003eBacteroidales, Clostridiales, and\u0026nbsp;\u003c/em\u003ewith\u003cem\u003e\u0026nbsp;\u003c/em\u003enon-classified taxa. \u0026nbsp;A significant negative correlation was found with \u003cem\u003eActinomycetales, Chlamydiales,\u003c/em\u003e and \u003cem\u003eTissierellales\u003c/em\u003e (Additional File 1: Figure S17). In general, there was significant variation across taxa involved in glycan degradation. The variation was especially dramatic in MA metagenomes of patients with large tumors. \u0026nbsp;\u003c/p\u003e"},{"header":"Discussion","content":"\u003cp\u003eHerein, we have described significant differences in the gut metagenomes adjoining the mucus layer (MA metagenome) between cervical cancer patients with small, early stage tumors and large, advanced stage tumors at the time of treatment. Specifically, MA metagenomes of patients with small tumors were enriched in molecular functions associated with biosynthesis and rapid proliferation. By contrast, MA metagenomes of patients with large tumors were enriched with functions associated with bacterial stress response and degradation of glycan. In addition, the utilization of ethanolamine and functions related to bacterial cells mobility were more abundant among young patients with large tumors and advanced stage (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ec), suggesting particularly aggressive tumor biology. Activity of glycan degradation and of ornithine biosynthesis were also significantly associated with the advance tumor stage and the presence of cancer cells in lymph nodes, while activity of the ribosome biogenesis significantly associated with low stage tumors and better recurrence-free survival. This dominance may lead to significant changes in the metabolic environment of the intestine that directly impact tumor growth and progression.\u003c/p\u003e \u003cp\u003eIt is possible, this aggressive tumor biology is fueled by differing molecular functions in MA metagenome. We have described two distinct microbial consortia, which we refer to here as proliferating and mucus degrading (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eg). Both consortia co-exist; although patients with larger tumor size and advanced stage show an obvious dominance of the mucus degrading consortia over proliferating. The tumor promoting effects of glycan degradation and utilization of major degradation products, such as galactose and other pentoses [\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e], may help drive this aggressive tumor growth. Glycan degradation leads to production of metabolites with known tumor growth promoting effects and with resistance to radio- and chemotherapy, such as ornithine [\u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e52\u003c/span\u003e] and ceramides [\u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e53\u003c/span\u003e, \u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e54\u003c/span\u003e]. The metabolites are produced as by-products of glycan degradation (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ee) and can be further metabolized either by bacteria or by host and give oncogenic effect after their metabolization. Ceramide itself, for example, is a powerful tumor suppressor, but products of its metabolism are potent tumor survival factors associated with resistance to cancer therapies [\u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e53\u003c/span\u003e]. Ornithine can be used by the consortium or by the host to synthesize the polyamines putrescine and spermidine (Additional File 1: Figure S8), which cause tumorigenic transformation and tumor progression [\u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e55\u003c/span\u003e]. If the gut metabolites cross the gut mucosa into blood through the mucus layer, they may promote systemic proliferation of tumor cells. In addition, the intensive ethanolamine utilization (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ee) by the LT-group metagenomes, especially pronounced in younger patients (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ec), also suggests a more pathogenic environment in the intestine. The set of enzymes involved in the pathway are similar to the \u003cem\u003eeut\u003c/em\u003e operon in \u003cem\u003eSalmonella enterica serovar Typhimurium\u003c/em\u003e [\u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e56\u003c/span\u003e], a known gastrointestinal pathogen. Many other species that contain the \u003cem\u003eeut\u003c/em\u003e genes are also pathogens [\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e], because ethanolamine is derived from the membrane phospholipid phosphatidylethanolamine, an important component of all bacterial and eukaryotic cells. The intensive utilization of ethanolamine may indicate a degradation of the colonic epithelium [\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e] that secretes peptides inhibiting bacterial penetration into the inner colonic mucus layer and blocking bacterial mobility [\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e, \u003cspan citationid=\"CR57\" class=\"CitationRef\"\u003e57\u003c/span\u003e]. Indeed, the activity of pathways associated with mobility, such as chemotaxis and flagella assembly, are also enhanced in MA metagenomes of the LT-group, particularly in sub-cluster of younger patients (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003ec).\u003c/p\u003e \u003cp\u003eConversely, the increased synthesis of peptidoglycan and lipopolysaccharides by the proliferating consortia may improve the intestinal microenvironment and the immune response in patients with smaller tumors [\u003cspan citationid=\"CR58\" class=\"CitationRef\"\u003e58\u003c/span\u003e]. The increased production of vitamins B7 (biotin) and B12 (cobalamin) may also be important for immune cell function and reduction of cellular oxidative stress [\u003cspan additionalcitationids=\"CR60\" citationid=\"CR59\" class=\"CitationRef\"\u003e59\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e61\u003c/span\u003e]. Overall, these metabolites may have a tumor suppressive effect.\u003c/p\u003e \u003cp\u003eThe findings suggest a potential therapeutic intervention in CC by shifting the balance in MA metagenome from mucus degrading consortium to proliferating one for suppression of the tumor growth and enhancement responses to treatment.\u003c/p\u003e \u003cp\u003eThis study does not make clear what comes first, the biological mechanism of mucus degradation or tumor progression. It is likely that both processes are tightly related and influence each other. Cervical cancer typically develops over years. Differences in the microbial flora of a large tumor may reflect changes in the tumor microenvironment directly affected by local tumor progression, or gut metagenome functions may allow increased tumor growth via systemic factors released into the bloodstream. Alternatively, tumor size may be a surrogate for other aggressive tumor biology, such as hypoxia or necrosis, which increase as tumors grow. Although this study associates large tumor size with degradation of mucus layer and reveal by-products of the degradation with known tumor promoting effects, it doesn\u0026rsquo;t provide direct evidence for an association of tumor-promoting metabolites with tumor size. It is possible that tumor progression and associated changes of organismal processes may be primary and drive changes in the mucus layer microenvironment and in molecular functions of the MA community through production of cytokines and metabolites [\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e] or by changing the mucin glycans that are also important host signals selecting microorganisms and making them less pathogenic [\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e]. Indeed, the study reveals stress-associated changes in the mucus layer microenvironment that are indicated by reduced proliferation, and by increased activity of quorum sensing and sporulation in mucus degrading consortium.\u003c/p\u003e \u003cp\u003eImportantly, taxonomic structure differed between groups only when functions and pathways were considered, suggesting taxonomic comparison alone is not adequate. Further studies using mouse or cell culture models may be necessary to decipher biological mechanisms underlying the discovered association between CC and mucus layer degradation. In addition, because potential mechanisms underlying the association of the tumor size and mucus layer degradation are not specific to CC, it is very likely that similar associations between the tumor size and structure and function of MA metagenomes may be seen in other cancer types.\u003c/p\u003e \u003cp\u003eShifting the balance between proliferating and mucus degrading consortia in MA metagenomes may directly affect cancer therapy and drug toxicity. The mucus layer plays a very important protective role in the intestine [\u003cspan citationid=\"CR62\" class=\"CitationRef\"\u003e62\u003c/span\u003e], therefore its unbalanced degradation can impair sensitivity to or ability to resist toxic effects of drugs. Recently, a growing body of evidence linked microorganisms to cancer therapy efficacy and toxicity [\u003cspan additionalcitationids=\"CR9\" citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e]. Our previous study of metagenomes in fecal samples of melanoma patients with different response to anti-PD-1 therapy [\u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e] reveals that degradation associated pathways are enriched in non-responders, while biosynthetic pathways are enriched in responders. We were not able to associate catabolic pathways with the mucus layer degradation in this study; however, we can speculate that the degradation can be responsible for the enrichment of catabolic processes in non-responding patients. The small number of samples and the use of fecal samples instead of rectal swabs might complicate identification of the mucus layer degradation in the study. It is likely that fecal samples are not optimal for the evaluation because they represent a different environment in the intestine. Further studies using WGS are necessary to compare fecal and swab samples from the same patient when mucus layer degradation is evaluated on association of patient\u0026rsquo;s response to cancer therapies, cancer stage, or patient survival.\u003c/p\u003e"},{"header":"Conclusions","content":"\u003cp\u003eOverall, these findings suggest potential interventions related to mucus enhancement that could improve outcomes and decrease cervical cancer progression. Further mechanistic studies are important to explore if the disturbed balance between proliferating and mucus degrading consortia can be corrected by diet, probiotics, antibiotics, or other interventions.\u003c/p\u003e "},{"header":"Methods","content":"\u003cp\u003e\u003cstrong\u003ePatient samples collection and processing\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eForty-one patients treated at the University of Texas, M.D. Anderson Cancer Center or the Lyndon B. Johnson Clinic (LBJ) at Harris Health with a diagnosis of CC participated in the study. The patients were enrolled on an Institutional Review Board (IRB) approved prospective protocol. \u0026nbsp;Informed consent was obtained to collect rectal swabs. \u0026nbsp;All patient samples were acquired prior to receiving any treatment. \u0026nbsp;\u003c/p\u003e\n\u003cp\u003eSamples were collected from each out of 41 patients by a clinician performing rectal exams using a matrix-designed quick-release Isohelix swab. The swabs were initially stored in 20ul of proteinase K and 400 \u0026mu;l of lysis buffer (Isohelix) within 1 h of sample collection and then were frozen and kept at -80\u0026deg;C.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eWGS sequencing and metagenome assembly\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWhole Genome Shotgun sequencing was performed on genomic bacterial DNA (gDNA), which was extracted to maximize bacterial DNA yield from specimens while keeping background amplification to a minimum [25, 26]. Libraries were constructed from each sample using the KAPA Hyper Prep Kit (Kapa Biosystems, Wilmington, MA, USA) and sequenced using the Illumina HiSeqX platform with the 2 x 150 bp paired-end read protocol. Sequencing reads were derived from raw BCL files which were retrieved from the sequencer and called into fastqs by Casava v1.8.3 (Illumina). The appropriate read preparation steps, such as quality control, trimming and filtering, and host DNA removal prior to further analysis, were performed using an in-house pipeline (Additional file 1: Figure S1b; Additional file 2: Table S1). Briefly, paired-end raw sequence reads were filtered and trimmed using BBMap [27]. The trimmed reads were mapped to a hg38 reference database (GCA_000001405.28) using bowtie2 [28] to remove host contamination. The cleaned reads were then assembled to longer sequences (contigs) using both MEGAHIT [29] and metaSPAdes [30]. The assembled contigs were filtered, to remove those that were smaller than 1000bp, and binned using MetaBAT2 [31]. The assembled and binned contigs were used for gene predictions by Prodigal [32]. Annotation of the genes by KEGG ortholog groups (KOs) was implemented \u0026nbsp;by KofamKOALA [33], and taxonomic classification of the contigs was done by CAT and BAT [34]. The read coverage of each assembled contig was calculated by aligning the cleaned reads directly to the contig using BBMap and by counting the mapped read by featureCounts [35]. The read coverage, GC content, taxonomic and functional annotations for each gene/contig were summarized using a Perl script. \u0026nbsp;All the software tools were running with default parameters if not specified. Versions and sources of the software tools or packages used in the pipeline are listed in Additional file 2: Table S1. Output of the assembly pipeline was a set of assembled contigs for each sample, their taxonomic annotation and read coverage, gene predictions for each contig, and functional annotation of each gene by KEGG Orthologous Group (KO) if found.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eComputational analysis of assembled genomes\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAnnotations of assembled contigs for all samples were aggregated into 1 table, referred to as the Metagenome Function Abundance (MFA) table (Additional file 1: Figure S1b). Each column in the table represents a metagenome, and each row represents a predicted known molecular function annotated by KO. Thus, each cell in the MFA table has a quantity of the specified (KO id) molecular function in the specified sample (Sample id). Quantification of the molecular functions by MFA table is explained using a toy example provided in Additional file 1: Figure S2. The MFA table was normalized using total number of reads in each sample and then multiplied by 1,000,000. Further analysis of the normalized MFA table included unsupervised and supervised methods, annotation by biological processes and pathways, and integration with clinicopathological characteristics of the patients (Additional file 1: Figure S1c).\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eUnsupervised analysis of MFA table\u003c/em\u003e\u003c/p\u003e\n\u003cp\u003eThe MFA table was filtered to select KOs found in at least 30 patients, log2-transformed and centered using the median, and then hierarchically clustered by the open source clustering software [36] with default parameters. The inferred clusters of samples were tested for association with clinical information including age, CC stage and tumor size, using fisher.test() and wilcox.test() functions in R. The inferred clusters of KOs were searched for overlap with known KEGG referenced pathways and modules using the \u0026ldquo;Search\u0026amp;Color Pathway\u0026rdquo; tool in KEGG mapper [37, 38] with further manual curation of the results.\u003cem\u003e\u0026nbsp;\u003c/em\u003e\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eSupervised analysis of MFA table\u003c/em\u003e\u003c/p\u003e\n\u003cp\u003eThe analysis was used to find differentially abundant KOs in MA metagenomes of patients with large versus small tumors. Samples for the analysis were selected by sorting all 41 samples by tumor size of the CC patients and assigning top 14 samples with largest tumors to LT-group and bottom 14 samples with smallest tumors to ST-group for the comparison. All KOs found in at least 1 sample were included in the analysis. Fisher\u0026rsquo;s test and Mann Whitney tests were used to find differentially abundant KOs between the group with p-value \u0026lt; 0.05 (either test) without adjustment. The KO was considered enriched in LT-group if Fisher\u0026rsquo;s test p-values \u0026lt; 0.05 and the KO is more common in the group. If the Fisher\u0026rsquo;s test p-value \u0026gt;0.05 but the Wilcoxon test p-value \u0026lt;0.05 than the enrichment was inferred by difference in mean abundances between ST- and LT-groups. The differentially abundant KOs were searched for overlap with known KEGG referenced pathways and modules as described above. The overlapping KOs were used to infer the pathways enrichment score calculated as ratio of difference between percentage of KOs overlapped with the pathway in ST- and LT-groups to sum of the percentages. Only top scored pathways that include 9 and more KOs for LT-group and 4 and more KOs for ST-group are considered.\u003c/p\u003e\n\u003cp\u003eAnnotation of \u0026nbsp;the differentially abundant KOs by \u0026nbsp;carbohydrate-active enzymes (CAZy) families \u0026nbsp;was implemented using the mapping table between the KO ID and \u0026nbsp;CAZy family ID\u0026nbsp;[39]. The table was downloaded from KEGG on 26 October 2019.\u003c/p\u003e\n\u003cp\u003eThe differentially abundant KOs were further used to quantify relative abundances of microbial subpopulations expressing them in ST and in LT group of samples. The decomposition of cell populations, referred to as ST-dominant and LT-dominant, was implemented based on characteristics of density plots of KO\u003c/p\u003e\n\u003cp\u003eThe differentially abundant KOs in ST and LT group were further used to construct 2 normal curves for ST-abundant and LT-abundant KOs. \u0026nbsp;Parameters of the curves, mean and standard deviation, were assumed to characterize molecular functions expressed by the microbial community dominating in either ST-group or LT-group, although each of the groups is comprised of both communities. To quantify abundances of the communities in each group, we fitted the 2 normal curves to density plots of all KOs found in the two groups using multiple linear regression. The obtained fitting coefficients were used to quantify the dominance of each set in metagenomes of each group.\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eSurvival analysis\u003c/em\u003e\u003c/p\u003e\n\u003cp\u003eThe analysis was used to associate the Recurrence Free Survival (RFS) probability with clinicopathological characteristics and KEGG pathway \u003cem\u003e[40]\u003c/em\u003e\u003cem\u003e.\u0026nbsp;\u003c/em\u003eActivity \u003cem\u003eof each\u0026nbsp;\u003c/em\u003epathway found to be associated with the tumor size was quantified for each patient using the mean value of log2-transformed normalized abundances of KOs involved in the pathway. The mean value is referred to as the pathway activity score.\u0026nbsp;To generate Kaplan-Meier plots for the pathway activity score, we categorized the score based on a cut-off set at the first quartile. Observations for each variable falling within the first quarter were labeled as \u0026ldquo;Low\u0026rdquo; and those greater than the first quartile cut-off were labelled as \u0026ldquo;High\u0026rdquo;. The analysis was implemented for each pathway identified as differentially enriched in either small or large tumor group metagenomes.\u0026nbsp;We used WHO standards to set cut-offs to categorize BMI into \u0026ldquo;Underweight/Normal weight\u0026rdquo; vs. \u0026ldquo;Overweight/Obese\u0026rdquo;. We used a cut-off set at the median for other continuous clinicopathological characteristics, namely, age and tumor size. Cox proportional hazards\u003cstrong\u003e\u0026nbsp;\u003c/strong\u003eregression analysis [41] was used to evaluate predictive value of KEGG pathways and clinical characteristics on RFS time. The R packages \u0026lsquo;survival\u0026rsquo; and \u0026lsquo;survminer\u0026rsquo; were used to compute the survival curves, and to visualize them as Kaplan-Meier plots.\u003c/p\u003e\n\u003cp\u003e\u003cem\u003eTaxonomic characterization of the cohort\u0026nbsp;\u003c/em\u003e\u003c/p\u003e\n\u003cp\u003eAll predicted contigs annotated by Length (L), Depths (D), and by taxonomy were used to quantify abundances of the unique set of species in 14 samples of LT group and in 14 samples of ST group. For each sample, depth of the contig was multiplied by its length and the values were summarized for all contigs that belong to the same taxa to quantify the taxa abundance.\u003c/p\u003e\n\u003cp\u003eThe values were normalized using the sum of the values for all contigs identified in the sample.\u003c/p\u003e\n\u003cp\u003eThe OTU (operating taxonomic units) table was created by merging the abundances of all identified taxa for all 28 samples. Two approaches were used to find\u0026nbsp;differentially abundant putative taxa between LT and ST group. In the first approach, the OTUs were selected by \u0026nbsp;Fisher\u0026rsquo;s exact test (for rare OTUs) and Mann Whitney tests (for common OTUs) with p-value cutoff \u0026lt; 0.05 (either test) without adjustment. In the other approach, the\u0026nbsp;Linear discriminant analysis (LDA) Effect Size (LEfSe) \u0026nbsp;[42] available as a Galaxy [43] module was used to determine the taxa that are differentially abundant between LT and ST groups. The analysis was run with default parameters, except the threshold\u0026nbsp;for the logarithmic LDA score for discriminative taxa. The threshold was set to 2.5 instead of 2.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eTo reveal taxa involved in glycan degradation pathway in each metagenome we have selected all contigs that encode enzymes involved in the pathway. The set of unique taxa from taxonomic annotation of the contigs was considered as the representative taxa of the pathway. To create an OTU table of the representative taxa, each OTU were quantified by multiplying the length and the depth of each selected contig and by summarizing the obtained values for all contigs encoding the OTU. The 100% stacked area plot was used to visualize taxonomic structure of the community involved in the pathway across all studied samples. The Pearson correlation coefficients were used to evaluate association of the log2-transformed score of the KEGG pathway with the abundance of each identified taxon.\u0026nbsp;\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eEthics approval and consent to participate\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eInformed written consent was obtained from eligible patients willing to participate in our study.\u003c/p\u003e\n\u003cp\u003eThis study was under a protocol approved by the UT M.D. Anderson Cancer Center Institutional Review Board (MDACC 2014\u0026ndash;0543)\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eConsent for publication\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eNot applicable.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAvailability of data and materials\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request. Raw sequencing reads generated by WGS are available in the BioProject database under accession PRJNA702617(https://www.ncbi.nlm.nih.gov/bioproject/PRJNA702617).\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCompeting interests\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe authors declare that they have no competing interests.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eFunding\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThis research was supported in part by the Radiological Society of North America Resident/Fellow Award (to L.E.C.); by the National Institutes of Health (NIH) through MD Anderson\u0026rsquo;s Cancer Center Support Grant P30CA016672 and through T32 grant CA101642-14 (T.T.S.); and by The University of Texas MD Anderson Cancer Center HPV-related Cancers Moonshot (L.E.C., A. K.). The funding sources were not involved in the development of the research hypothesis, study design, data analysis, or manuscript writing.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthors\u0026rsquo; Contributions\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eLEC, AHK, and TVK designed the study and wrote the paper with contribution from XW, JZ, AF, PJE, AJ, LLL, KMS, NA, RRJ, MBEA, and LR.\u003c/p\u003e\n\u003cp\u003eTS, TTS, \u0026nbsp;KYC, EL, MAK, GB, JY, JFP, MPM, \u0026nbsp;and PO performed the experiment and sequencing.\u003c/p\u003e\n\u003cp\u003eXW, JZ, TVK, XS and HS developed the computation pipeline for metagenome assembly and annotation.\u003c/p\u003e\n\u003cp\u003eTVK, XW, and MBEA processed the data.\u003c/p\u003e\n\u003cp\u003eAll authors read and approved the final manuscript.\u003cstrong\u003e\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgements\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eWe sincerely thank the patients who volunteered to participate in this study.\u0026nbsp;\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eWe acknowledge the Moonshot program at the University of Texas MD Anderson Cancer Center for financial support of the study and High-Performance Computing at the University of Texas MD Anderson Cancer Center for providing computational resources that have contributed to the research results reported in this paper. \u0026nbsp;\u0026nbsp;\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n \u003cli\u003e\u003cspan\u003eVu M, Yu J, Awolude OA, Chuang L. Cervical cancer worldwide. Curr Probl Cancer. 2018;42:457\u0026ndash;65.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eSzymonowicz KA, Chen JJ. Biological and clinical aspects of HPV-related cancers. Cancer Biology Medicine. 2020;17:864\u0026ndash;78.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eCurty G, de Carvalho PS, Soares MA. The Role of the Cervicovaginal Microbiome on the Genesis and as a Biomarker of Premalignant Cervical Intraepithelial Neoplasia and Invasive Cervical Cancer. \u003cem\u003eInternational Journal of Molecular Sciences\u003c/em\u003e 2020, 21.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eSims TT, Colbert LE, Zheng J, Delgado Medrano AY, Hoffman KL, Ramondetta L, Jazaeri A, Jhingran A, Schmeler KM, Daniel CR, Klopp A. Gut microbial diversity and genus-level differences identified in cervical cancer patients versus healthy controls. Gynecol Oncol. 2019;155:237\u0026ndash;44.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eSims TT, El Alam MB, Karpinets TV, Dorta-Estremera S, Hegde VL, Nookala S, Yoshida-Court K, Wu XG, Biegert GWG, Medrano AYD, et al: Gut microbiome diversity is an independent predictor of survival in cervical cancer patients receiving chemoradiation. \u003cem\u003eCommunications Biology\u003c/em\u003e 2021, 4.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eKarpinets TV, Solley TN, Mikkelson MD, Dorta-Estremera S, Nookala SS, Medrano AYD, Petrosino JF, Mezzari MP, Zhang J, Futreal PA, et al. Effect of Antibiotics on Gut and Vaginal Microbiomes Associated with Cervical Cancer Development in Mice. Cancer Prev Res (Phila). 2020;13:997\u0026ndash;1006.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eAlhinai EA, Walton GE, Commane DM. The Role of the Gut Microbiota in Colorectal Cancer Causation. \u003cem\u003eInt J Mol Sci\u003c/em\u003e 2019, 20.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eZitvogel L, Galluzzi L, Viaud S, Vetizou M, Daillere R, Merad M, Kroemer G. Cancer and the gut microbiota: an unexpected link. Sci Transl Med. 2015;7:271ps271.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eAlexander JL, Wilson ID, Teare J, Marchesi JR, Nicholson JK, Kinross JM. Gut microbiota modulation of chemotherapy efficacy and toxicity. Nat Rev Gastroenterol Hepatol. 2017;14:356\u0026ndash;65.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eGopalakrishnan V, Spencer CN, Nezi L, Reuben A, Andrews MC, Karpinets TV, Prieto PA, Vicente D, Hoffman K, Wei SC, et al. Gut microbiome modulates response to anti-PD-1 immunotherapy in melanoma patients. Science. 2018;359:97\u0026ndash;103.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eProenca JT, Barral DC, Gordo I. Commensal-to-pathogen transition: One-single transposon insertion results in two pathoadaptive traits in Escherichia coli -macrophage interaction. Sci Rep. 2017;7:4504.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eLu HP, Lai YC, Huang SW, Chen HC, Hsieh CH, Yu HT. Spatial heterogeneity of gut microbiota reveals multiple bacterial communities with distinct characteristics. \u003cem\u003eScientific Reports\u003c/em\u003e 2014, 4.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eBassis CM, Moore NM, Lolans K, Seekatz AM, Weinstein RA, Young VB, Hayden MK, Program CDCPE. Comparison of stool versus rectal swab samples and storage conditions on bacterial community profiles. BMC Microbiol. 2017;17:78.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eReyman M, van Houten MA, Arp K, Sanders EAM, Bogaert D. Rectal swabs are a reliable proxy for faecal samples in infant gut microbiota research based on 16S-rRNA sequencing. Sci Rep. 2019;9:16072.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eMentella MC, Scaldaferri F, Pizzoferrato M, Gasbarrini A, Miggiano GAD. Nutrition, IBD and Gut Microbiota: A Review. \u003cem\u003eNutrients\u003c/em\u003e 2020, 12.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eJohansson ME, Phillipson M, Petersson J, Velcich A, Holm L, Hansson GC. The inner of the two Muc2 mucin-dependent mucus layers in colon is devoid of bacteria. Proc Natl Acad Sci U S A. 2008;105:15064\u0026ndash;9.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eSchroeder BO. Fight them or feed them: how the intestinal mucus layer manages the gut microbiota. Gastroenterology Report. 2019;7:3\u0026ndash;12.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eHolmen Larsson JM, Thomsson KA, Rodriguez-Pineiro AM, Karlsson H, Hansson GC. Studies of mucus in mouse stomach, small intestine, and colon. III. Gastrointestinal Muc5ac and Muc2 mucin O-glycan patterns reveal a regiospecific distribution. Am J Physiol Gastrointest Liver Physiol. 2013;305:G357\u0026ndash;63.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eCapon C, Maes E, Michalski JC, Leffler H, Kim YS. Sd(a)-antigen-like structures carried on core 3 are prominent features of glycans from the mucin of normal human descending colon. Biochem J. 2001;358:657\u0026ndash;64.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eJuge N. Microbial adhesins to gastrointestinal mucus. Trends Microbiol. 2012;20:30\u0026ndash;9.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eCai R, Cheng C, Chen J, Xu X, Ding C, Gu B. Interactions of commensal and pathogenic microorganisms with the mucus layer in the colon. Gut Microbes. 2020;11:680\u0026ndash;90.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eKavanaugh NL, Zhang AQ, Nobile CJ, Johnson AD, Ribbeck K. Mucins suppress virulence traits of Candida albicans. mBio. 2014;5:e01911.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eWheeler KM, Carcamo-Oyarce G, Turner BS, Dellos-Nolan S, Co JY, Lehoux S, Cummings RD, Wozniak DJ, Ribbeck K. Mucin glycans attenuate the virulence of Pseudomonas aeruginosa in infection. Nat Microbiol. 2019;4:2146\u0026ndash;54.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eShort MI, Hudson R, Besasie BD, Reveles KR, Shah DP, Nicholson S, Johnson-Pais TL, Weldon K, Lai Z, Leach RJ, et al. Comparison of rectal swab, glove tip, and participant-collected stool techniques for gut microbiome sampling. BMC Microbiol. 2021;21:26.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eHuman Microbiome Project C. A framework for human microbiome research. Nature. 2012;486:215\u0026ndash;21.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eHuman Microbiome Project C. Structure, function and diversity of the healthy human microbiome. Nature. 2012;486:207\u0026ndash;14.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003e2020 BhsnpbfAM.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eLangmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357\u0026ndash;9.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eLi D, Liu CM, Luo R, Sadakane K, Lam TW. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics. 2015;31:1674\u0026ndash;6.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eNurk S, Meleshko D, Korobeynikov A, Pevzner PA. metaSPAdes: a new versatile metagenomic assembler. Genome Res. 2017;27:824\u0026ndash;34.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eKang DD, Li F, Kirton E, Thomas A, Egan R, An H, Wang Z. MetaBAT 2: an adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies. PeerJ. 2019;7:e7359.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eHyatt D, Chen GL, Locascio PF, Land ML, Larimer FW, Hauser LJ. Prodigal: prokaryotic gene recognition and translation initiation site identification. BMC Bioinformatics. 2010;11:119.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eAramaki T, Blanc-Mathieu R, Endo H, Ohkubo K, Kanehisa M, Goto S, Ogata H. KofamKOALA: KEGG Ortholog assignment based on profile HMM and adaptive score threshold. Bioinformatics. 2020;36:2251\u0026ndash;2.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003evon Meijenfeldt FAB, Arkhipova K, Cambuy DD, Coutinho FH, Dutilh BE. Robust taxonomic classification of uncharted microbial sequences and bins with CAT and BAT. Genome Biol. 2019;20:217.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eLiao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30:923\u0026ndash;30.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003ede Hoon MJ, Imoto S, Nolan J, Miyano S. Open source clustering software. Bioinformatics. 2004;20:1453\u0026ndash;4.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eKanehisa M, Sato Y, Kawashima M, Furumichi M, Tanabe M. KEGG as a reference resource for gene and protein annotation. Nucleic Acids Res. 2016;44:D457\u0026ndash;62.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eKanehisa M, Sato Y. KEGG Mapper for inferring cellular functions from protein sequences. Protein Sci. 2020;29:28\u0026ndash;35.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eLombard V, Ramulu HG, Drula E, Coutinho PM, Henrissat B. The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 2014;42:D490\u0026ndash;5.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eLi JCA. Modeling survival data: Extending the Cox model. Sociological Methods Research. 2003;32:117\u0026ndash;20.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eCox DR. Regression Models and Life-Tables. Journal of the Royal Statistical Society Series B-Statistical Methodology. 1972;34:187-+.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eSegata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, Huttenhower C. Metagenomic biomarker discovery and explanation. \u003cem\u003eGenome Biology\u003c/em\u003e 2011, 12.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eAfgan E, Baker D, Batut B, van den Beek M, Bouvier D, Cech M, Chilton J, Clements D, Coraor N, Gruning BA, et al. The Galaxy platform for accessible, reproducible and collaborative biomedical analyses: 2018 update. Nucleic Acids Res. 2018;46:W537\u0026ndash;44.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eClamp JR, Cooper B, Creeth JM, Ene D, Barrett J, Gough M. The Presence of Polysaccharide in Normal Human Gastric Mucus. Biochem J. 1983;215:421\u0026ndash;3.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eStincone A, Prigione A, Cramer T, Wamelink MMC, Campbell K, Cheung E, Olin-Sandoval V, Gruning NM, Kruger A, Alam MT, et al. The return of metabolism: biochemistry and physiology of the pentose phosphate pathway. Biol Rev. 2015;90:927\u0026ndash;63.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eKaval KG, Garsin DA. Ethanolamine Utilization in Bacteria. \u003cem\u003emBio\u003c/em\u003e 2018, 9.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003ePegg AE. Mammalian Polyamine Metabolism and Function. Iubmb Life. 2009;61:880\u0026ndash;94.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eVazquez A, Kamphorst JJ, Markert E, Schug ZT, Tardito S, Gottlieb E. Cancer metabolism at a glance. J Cell Sci. 2016;129:3367\u0026ndash;73.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eMatilla MA, Krell T. The effect of bacterial chemotaxis on host infection and pathogenicity. Fems Microbiology Reviews. 2018;42:40\u0026ndash;67.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eSrividhya KV, Krishnaswamy S. Subclassification and targeted characterization of prophage-encoded two-component cell lysis cassette. J Biosci. 2007;32:979\u0026ndash;90.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eMcComas KN, Torgeson AM, Ager BJ, Hellekson C, Burt LM, Maurer KA, Werner TL, Gaffney DK. The variable impact of positive lymph nodes in cervical cancer: Implications of the new FIGO staging system. Gynecol Oncol. 2020;156:85\u0026ndash;92.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eGeck RC, Foley JR, Murray Stewart T, Asara JM, Casero RA Jr, Toker A. Inhibition of the polyamine synthesis enzyme ornithine decarboxylase sensitizes triple-negative breast cancer cells to cytotoxic chemotherapy. J Biol Chem. 2020;295:6263\u0026ndash;77.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eYura Y, Masui A, Hamada M. Inhibitors of Ceramide- and Sphingosine-Metabolizing Enzymes as Sensitizers in Radiotherapy and Chemotherapy for Head and Neck Squamous Cell Carcinoma. \u003cem\u003eCancers\u003c/em\u003e 2020, 12.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eMorad SAF, Cabot MC. Ceramide-orchestrated signalling in cancer cells. Nat Rev Cancer. 2013;13:51\u0026ndash;65.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eMurray-Stewart TR, Woster PM, Casero RA. Targeting polyamine metabolism for cancer therapy and prevention. Biochem J. 2016;473:2937\u0026ndash;53.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eRoof DM, Roth JR. Ethanolamine utilization in Salmonella typhimurium. J Bacteriol. 1988;170:3855\u0026ndash;63.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eBergstrom JH, Birchenough GMH, Katona G, Schroeder BO, Schutte A, Ermund A, Johansson MEV, Hansson GC. Gram-positive bacteria are held at a distance in the colon mucus by the lectin-like protein ZG16. Proc Natl Acad Sci USA. 2016;113:13833\u0026ndash;8.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eTang C, Ding R, Sun J, Liu J, Kan J, Jin C. The impacts of natural polysaccharides on intestinal microbiota and immune responses - a review. Food Funct. 2019;10:2290\u0026ndash;312.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eKung JT, Mackenzie CG, Talmage DW. The requirement for biotin and fatty acids in the cytotoxic T-cell response. Cell Immunol. 1979;48:100\u0026ndash;10.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003ePeterson CT, Rodionov DA, Osterman AL, Peterson SN. B Vitamins and Their Role in Immune Regulation and Cancer. \u003cem\u003eNutrients\u003c/em\u003e 2020, 12.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eAgrawal S, Agrawal A, Said HM. Biotin deficiency enhances the inflammatory response of human dendritic cells. American Journal of Physiology-Cell Physiology. 2016;311:C386\u0026ndash;91.\u003c/span\u003e\u003c/li\u003e\n \u003cli\u003e\u003cspan\u003eCamilleri M. Leaky gut: mechanisms, measurement and clinical implications in humans. Gut. 2019;68:1516\u0026ndash;26.\u003c/span\u003e\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":true,"hideJournal":false,"highlight":"","institution":"","isAcceptedByJournal":true,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":true,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"bmc-cancer","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"bcan","sideBox":"Learn more about [BMC Cancer](http://bmccancer.biomedcentral.com/)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/bcan/default.aspx","title":"BMC Cancer","twitterHandle":"BMC_series","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"em","reportingPortfolio":"BMC Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true},"keywords":"Metagenomics, mucus layer, gut, cervical cancer, tumor size, bacteria, Bacteroidales, Clostridiales","lastPublishedDoi":"10.21203/rs.3.rs-796525/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-796525/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003ch2\u003eBackground\u003c/h2\u003e \u003cp\u003eGut microbiome community composition differs between cervical cancer (CC) patients and healthy controls, and increased gut diversity is associated with improved outcomes after treatment. We proposed that functions of specific microbial species adjoining the mucus layer may directly impact the biology of CC.\u003c/p\u003e\u003ch2\u003eResults\u003c/h2\u003e \u003cp\u003eIn this study, we examined metagenomes of rectal swabs in 41 CC patients using whole-genome shotgun sequencing and found a significant association between molecular functions encoded by the metagenomes with markers of aggressive cancer including initial tumor size and stage. Profiling of the molecular function abundances and their distributions identified 2 microbial communities co-existing in each metagenome but with distinct metabolism and taxonomic structures. Community A (\u003cem\u003eClostridia\u003c/em\u003e and \u003cem\u003eProteobacteria\u003c/em\u003e predominant) was characterized by high activity of pathways involved in stress response, mucus glycan degradation and utilization of degradation byproducts. This community was prevalent in larger, advanced stage tumors. Conversely, community B (\u003cem\u003eBacteroidia\u003c/em\u003e predominant) was characterized by fast growth, active oxidative phosphorylation, and production of vitamins. This community was prevalent in small, early-stage tumors.\u003c/p\u003e\u003ch2\u003eConclusions\u003c/h2\u003e \u003cp\u003eBased on these results, we propose that increased mucus layer degradation is associated with a more aggressive cervical cancer phenotype.\u003c/p\u003e","manuscriptTitle":"Gut mucus layer degradation is associated with aggressive cervical cancer phenotype","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2021-12-30 16:52:42","doi":"10.21203/rs.3.rs-796525/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"bmc-cancer","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":false,"externalIdentity":"bcan","sideBox":"Learn more about [BMC Cancer](http://bmccancer.biomedcentral.com/)","snPcode":"","submissionUrl":"https://www.editorialmanager.com/bcan/default.aspx","title":"BMC Cancer","twitterHandle":"BMC_series","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"em","reportingPortfolio":"BMC Series","inReviewEnabled":true,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"39323024-7176-4f39-b426-e7a5595c00ce","owner":[],"postedDate":"December 30th, 2021","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"under-review","subjectAreas":[],"tags":[],"updatedAt":"2022-08-12T07:23:08+00:00","versionOfRecord":[],"versionCreatedAt":"2021-12-30 16:52:42","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-796525","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-796525","identity":"rs-796525","version":["v1"]},"buildId":"GqpaHPwrfC8PjnIFayRh5","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}