AbstractMajor depressive disorder (MDD) is a debilitating mental disease, but its underlying molecular mechanisms remain obscure. Our previously established model of naturally occurring depression-like (DL) behaviors in Macaca fascicularis, which is characterized by microbiota-gut–brain (MGB) axis disturbances, can be used to interrogate how a disturbed gut ecosystem may impact the molecular pathology of MDD. Here, gut metagenomics were used to characterize how gut virus and bacterial species, and associated metabolites, change in depression-like monkey model. We identified a panel of 33 gut virus and 14 bacterial species that could discriminate the depression-like from control macaques. In addition, using lipidomic analyses of central and peripheral samples obtained from these animals, we found that the DL macaque were characterized by alterations in the relative abundance, carbon-chain length, and unsaturation degree of 1,2-diacylglyceride (DG) in the prefrontal cortex (PFC), in a brain region-specific manner. In addition, lipid-reaction analysis identified more active and inactive lipid pathways in PFC than in amygdala or hippocampus, with DG being a key nodal player in these lipid pathways. Significantly, co-occurrence network analysis showed that the DG levels may be relevant to the onset of negative emotions behaviors in PFC. Together our findings suggest that altered DG levels and structure in the PFC are hallmarks of the DL macaque, thus providing a new framework for understanding the gut microbiome’s role in depression.