Abstract
Human endogenous retroviruses (HERVs) are the remnants of ancient retroviral infections integrated into the human genome. Although most HERVs are silenced or rendered inactive by various regulatory mechanisms, they retain the potential to influence the nearby genes. We analyzed the regulatory map of 91 HERV-Ks on neighboring genes in human breast cancer and investigated the impact of HERV-Ks on the tumor microenvironment (TME) and prognosis of breast cancer. Nine RNA-seq datasets were obtained from GEO and NCBI SRA. Differentially expressed genes and HERV-Ks were analyzed using DESeq2. Validation of high-risk prognostic candidate genes using TCGA data. These included Overall survival (multivariate Cox regression model), immune infiltration analysis (TIMER), tumor mutation burden (maftools), and drug sensitivity analysis (GSCA). A total of 88 candidate genes related to breast cancer prognosis were screened, of which CD48, SLAMF7, SLAMF1, IGLL1, IGHA1, and LRRC8A were key genes. Functionally, these six key genes were significantly enriched in some immune function-related pathways, which may be associated with poor prognosis for breast cancer (p = 0.00016), and the expression levels of these genes were significantly correlated with the sensitivity of breast cancer treatment-related drugs. Mechanistically, they may influence breast cancer development by modulating the infiltration of various immune cells into the TME. We further experimentally validated these genes to confirm the results obtained from bioinformatics analysis. This study represents the first report on the regulatory potential of HERV-K in the neighboring breast cancer genome. We identified three key HERV-Ks and five neighboring genes that hold promise as novel targets for future interventions and treatments for breast cancer.