IntroductionThe high incidence of breast cancer (BC) prompted us to explore more factors that might affect its occurrence, development, treatment, and also recurrence. Dysregulation of cholesterol metabolism has been widely observed in BC; however, the detailed role of how cholesterol metabolism affects chemo-sensitivity, and immune response, as well as the clinical outcome of BC is unknown. MethodsWith Mendelian randomization (MR) analysis, the potential causal relationship between genetic variants of cholesterol and BC risk was assessed first. Then we analyzed 73 cholesterol homeostasis-related genes (CHGs) in BC samples and their expression patterns in the TCGA cohort with consensus clustering analysis, aiming to figure out the relationship between cholesterol homeostasis and BC prognosis. Based on the CHG analysis, we established a CAG_score used for predicting therapeutic response and overall survival (OS) of BC patients. Furthermore, a machine learning method was adopted to accurately predict the prognosis of BC patients by comparing multi-omics differences of different risk groups. ResultsWe observed that the alterations in plasma cholesterol appear to be correlative with the venture of BC (MR Egger, OR: 0.54, 95% CI: 0.35-0.84, p