Gut microbiota and metabolites have been identified to contribute to the pathogenesis of functional constipation (FC); however, the underlying mechanism(s) have not been elucidated, and the relationship between the gut microbiota and metabolites in FC has received limited attention in the literature. 16S rDNA sequencing and non-targeted metabolomic detection based on liquid chromatography-mass spectrometry (LC-MS/MS) technologies were combined to analyze the altered gut microbiome and metabolic profile of fecal samples from FC patients and healthy individuals (healthy control; HC). The richness and diversity of gut microbiota significantly (p < 0.01) increased in FC patients. Compared to the HC group, 18 genera, including Intestinibacter, Klebsiella, and Akkermansia, exhibited statistically significant changes (p < 0.05). Metabolic analysis showed that metabolic profiles were also markedly altered with 79 metabolites, such as (-)-caryophyllene oxide, chenodeoxycholic acid, and biliverdin, indicating significant inter-group differences (p < 0.05). Besides, the primary bile acid biosynthesis, as well as the metabolic profile of porphyrin and chlorophyll, were the most dominant enriched pathways (FDR < 0.01), in which chenodeoxycholic acid and biliverdin were significantly enriched, respectively. Correlation analysis demonstrated a strong relationship between 10 genera and 19 metabolites (r > 0.6, FDR < 0.05), and notably, Intestinibacter showed a negative correlation with biliverdin (FDR < 0.001), which highlighted the interplay of the gut microbiota and metabolites in the pathogenesis of FC. Our research describes the characteristics of the gut microbiota and metabolic profiles and the correlation between the gut microbiota and metabolites in FC patients. This may contribute to the understanding of the underlying mechanisms involved in FC pathogenesis and may provide novel insights into therapeutic interventions.