BackgroundGut microbiome (GM) was observed to be associated with the incidence of Hirschsprung disease (HD). However, the effect and mechanism of GM in HD is still unclear. To investigate the relationship between GM and HD and the effect of metabolites as mediators, a bidirectional two-step Mendelian randomization (MR) study was conducted.MethodsThe study selected instrument variables (IVs) from summary-level genome-wide association studies (GWAS). The MiBioGen consortium provided the GWAS data for GM, while the GWAS data for metabolites and HD were obtained from the GWAS Catalog consortium. Two-sample MR analyses were performed to estimate bidirectional correlations between IVs associated with GM and HD. Then, genetic variants related to 1,400 metabolite traits were selected for further mediation analyses using the Product method.ResultsThis study found that seven genus bacteria had a significant causal relationship with the incidence of HD but not vice versa. 27 metabolite traits were significantly correlated with HD. After combining the significant results, three significant GM-metabolites-HD lines have been identified. In the Peptococcus-Stearoyl sphingomyelin (d18:1/18:0)-HD line, the Stearoyl sphingomyelin (d18:1/18:0) levels showed a mediation proportion of 14.5%, while in the Peptococcus-lysine-HD line, the lysine levels had a mediation proportion of 12.9%. Additionally, in the Roseburia-X-21733-HD line, the X-21733 levels played a mediation proportion of 23.5%.ConclusionOur MR study indicates a protective effect of Peptococcus on HD risk that is partially mediated through serum levels of stearoyl sphingomyelin (d18:1/18:0) and lysine, and a risk effect of Roseburia on HD that is partially mediated by X-21733 levels. These findings could serve as novel biomarkers and therapeutic targets for HD.