Patients and sample collection
32 female participants were assigned into 2 groups: group NH, meaning chronic radiative proctitis (CRP) patients with no or mild hematochezia (hematochezia frequency < twice per week and/or no fresh blood observed under proctoscope) (n = 16); and group H, meaning CRP patients with hematochezia (hematochezia frequency > twice per week and/or fresh blood was found under proctoscope) (n = 13) (Table 1). Two patients presented with rectal ulcers (diameter >1mm) and/or bacterial infection were excluded. Age of the 30 patients ranged from 28y to 71y. Patients have received total dose of 48 to 54 Gy irradiation as for radiotherapy of uterus or cervical carcinoma but have not received antibiotic treatment for at least 3 months. 3 samples from healthy adults were collected and analyzed in parallel.
Diversity of gut microbiota in CRP patients.
To investigate the ecological complexity of the gut microbiota of CRP patients, alpha diversity analysis within each group was generated using Community richness calculators Chao1 estimator and ace estimator as well as Community diversity calculators the Simpson index and the Shannon index in MOTHUR. As shown in Figure 1, the Community richness is lower in CRP patients with hematochezia (Mean of Chao1 = 1065; Mean of ace = 1738) as compare with CRP patients without hematochezia (Mean of Chao1 = 943.0; Mean of ace = 1488) but not statistically different (p > 0.1). Nevertheless, intragroup diversities are similar between the two groups. This result suggests that the overall complexity of gut microbiota may not be a sensitive indicator for hematochezia of CRP.
Differential gut microbiota of CRP patients with hematochezia
Venn diagram was generated to illustrate the assigned OTUs in each group. About 70% of the OTUs were shared in CRP patients while 27% and 29% of the taxonomic units were unique in H group and NH group, respectively (Figure 2a). Ordination analysis including Principal Component Analysis (PCA), Principal Co-ordinates Analysis (PCoA) and Non-Metric Multi-Dimensional Scaling (NMDS) were conducted to measure the difference in bacterial community composition in two-dimension. Inter-individual variation accounted for a part of variance, which was consistent with previous studies8,15. The ordination plots demonstrated a difference among healthy controls, CRP patients and CRP patients with hematochezia, suggesting that a proportion of bacteria population is related to the hematochezia symptom of CRP (Figure 2b, Figure S1-2).
To further study the microbiota composition, the OUT matrixes were visualized by Krona (Figure S3). As one of the most common bacterial components of the gut microbiome, Bacteroides was the most abundant bacteria in both of the CRP groups and that in healthy controls (Figure 2c, Figure S3). Composition of Bacteroides was higher in the NH group than that in H group. More Escherichia/Shigella, Fusobacterium and Prevotella were found in gut microbiota of CRP patients with hematochezia as compare to the NH group.
The linear discriminant analysis on effective size (LefSe) generated a global view of gut microbiota and revealed 8 bacteria taxa with differential relative abundance in both NH and H groups (LDA score > 2). The Verrucomicrobia phylum, Mogibacterium which is a species of the Clostridiales Incertae Sedis XIII genus, and Allisonella were identified from the H group. Akkermansia is the most abundant genus within the Verrucomicrobias (account for >99% of the Verrucomicrobias; 2% of the microbiota in the NH group; 4% in the H group) while the other two genera of Verrucomicrobias accounted for less than 0.02% of the microbiota. Our result also showed that the compositions of the Enterobacteriales family and the Porphyromonadaceae genus were lower in the H group (Figure 3, Figure S4). The linear discriminant analysis result was consistent with the Krona plots, as the Escherichia/Shigella is the major species of Enterobacteriales and the Parabacteroides belongs to the Porphyromonadaceaes.
Microbial functional dysbiosis in Hematochezia group
As the composition of microbiota community altered, functional and metabolic changes could be predicted based on previous molecular studies16,17. The Fe-S protein (COG1600), the Glutathione peroxidase (GPx, COG0386) and the Glutaredoxin-related protein (COG4545) were significantly increased in Hematochezia group (Figure 4, Table S1). The GPx is an important enzyme for detoxin and clearage of hydrogen peroxide. Notably, the metabolism of Arachidonic Acid (ARA) was also higher in Hematochezia group according to the annotation of the KEGG database (Figure S5). ARA is not only an important composition of cell membrane but also crucial for biosynthesis of prostaglandin (PGE2 and PGI2) which are commonly known as inflammatory mediators.