The gut microbiota has been demonstrated to play a significant role in the pathogenesis of Parkinson's disease (PD). However, conflicting findings regarding specific microbial species have been reported, possibly due to confounding factors within human populations. Herein, our current study investigated the interaction between the gut microbiota and host in a non-human primate (NHP) PD model induced by 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) using a multi-omic approach and a self-controlled design. Our transcriptomic sequencing of peripheral blood leukocytes (PBL) identified key genes involved in pro-inflammatory cytokine dysregulation, mitochondrial function regulation, neuroprotection activation, and neurogenesis associated with PD, such as IL1B, ATP1A3, and SLC5A3. The metabolomic profiles in serum and feces consistently exhibited significant alterations, particularly those closely associated with inflammation, mitochondrial dysfunctions and neurodegeneration in PD, such as TUDCA, ethylmalonic acid, and L-homophenylalanine. Furthermore, fecal metagenome analysis revealed gut dysbiosis associated with PD, characterized by a significant decrease in alpha diversity and altered commensals, particularly species such as Streptococcus, Butyrivibrio, and Clostridium. Additionally, significant correlations were observed between PD-associated microbes and metabolites, such as sphingomyelin and phospholipids. Importantly, PDPC significantly reduced in both PD monkey feces and serum, exhibiting strong correlation with PD-associated genes and microbes, such as SLC5A3 and Butyrivibrio species. Moreover, such multi-omic differential biomarkers were linked to the clinical rating scales of PD monkeys. Our findings provided novel insights into understanding the potential role of key metabolites in the host-microbiota interaction involved in PD pathogenesis.