Responses to neoadjuvant chemoradiotherapy (nCRT) and therapy-related toxicities in rectal cancer vary among patients. To provide the individualized therapeutic option for each patient, predictive markers of therapeutic responses and toxicities are in critical need. We aimed to identify the association of gut microbiome with and its potential predictive value for therapeutic responses and toxicities. In the present study, we collected fecal microbiome samples from patients with rectal cancer at treatment initiation and just after nCRT. Taxonomic profiling 16S ribosomal RNA gene sequencing was performed on all samples. Patients were classified as responders versus non-responders. Patients were grouped into no or mild diarrhea and severe diarrhea. STAMP and high-dimensional class comparisons linear discriminant analysis of effect size (LEfSe) were used to compare the compositional differences between groups. Phylogenetic Investigation of Communities by Reconstruction of Unobserved States (PICRUSt) was utilized to predict differences in metabolic function between groups. Ten patients were classified as responders and 12 patients were classified as non-responders. Fourteen patients experienced no or mild diarrhea and 8 patients experienced severe diarrhea. Several bacteria taxa with significantly different relative abundances before and after nCRT were identified. Similarly, several baseline bacteria taxa and predicted pathways with significantly different relative abundances between responders and non-responders or between patients no or mild diarrhea and severe diarrhea were identified. Specifically, was identified as enriched in responders and several bacteria taxa in the order etc. were identified as enriched in non-responders. Pathways including fatty acid metabolism were predicted to be enriched in responders. In addition, , , and etc. were identified as enriched in patients with no or mild diarrhea. Pathways including primary bile acid biosynthesis were predicted to be enriched in patients with no or mild diarrhea. Together, the microbiota and pathway markers identified in this study may be utilized to predict the therapeutic responses and therapy-related toxicities of nCRT in patients with rectal cancer. More patient data is needed to verify the current findings and the results of metagenomic, metatranscriptomic, and metabolomic analyses will further mine key biomarkers at the compositional and functional level.