To the Editor: Breast cancer is one of the malignant diseases that cause death in women and is a severe threat to women's health. With the progress of medical treatment, there are many methods to treat breast cancer, such as drug therapy and hormone therapy. Among them, molecular targeted therapy has dramatically improved the treatment effect of breast cancer. Therefore, it is vital to find important molecular markers.[1] Breast cancer can be divided into four subtypes: triple-negative (TN), lumA, lumB, and HER2+.[2] HER2+ breast cancer accounts for 15–20% of breast cancers, with a higher grade, a more aggressive phenotype, and a worse prognosis.[3]
The main objective of this study was to screen for key genes that appear only in the HER2+ subtype of breast cancer relative to the other three subtypes, including basal, lumA, and lumB. Screening for these key genes can help determine how HER2+ subtype differs from other subtypes in terms of pathogenesis, then treat them in a more targeted manner.
The dataset was obtained from the The Cancer Genome Atlas (TCGA) database (https://portal.gdc.cancer.gov/) and included 47 HER2+ samples, 87 basal samples, 291 lumA samples, 118 lumB samples, and 97 samples from normal tissues. The clinical features of HER2+ samples are shown in [Supplementary Table 1, https://links.lww.com/CM9/B268]. Differential gene analysis was done for each of the four subtypes vs. normal samples using the limma package in the R language. According to the screening condition |logFC| >2.0, P <0.05, a total of 2063 differentially expressed genes were obtained from HER2+ samples, of which 897 genes were upregulated, and 1166 genes were downregulated [Figure 1A]. The numbers of differential genes obtained for basal, lumA, and lumB under the same screening conditions were 2093, 1365, and 2084, respectively [Figure 1B].
Figure 1: (A) Volcano map of HER2+ breast cancer samples. (B) Volcano map of basal, lumA and lumB breast cancer samples. And the Venn diagram of genes among DEGs of three subgroups (basal, lumA, lumB) and HER2+. (C) Hub genes of HER2+ breast cancer samples. (D) Protein-protein interaction networks of key genes of HER2+ breast cancer samples. (E) GO and KEGG analysis for the hub genes. (F) Validation of expression levels of the three hub genes (UTS2, DRD4, GLP1R) among basal, HER2+, lumA, lumB, and normal tissues from the TCGA database. (G,H) OS analysis of UTS2, DRD4, and GLP1R. (I) Immunohistochemistry. The left panel shows the expression of UTS2 in normal tissues and the right panel shows the expression of UTS2 in breast cancer. (J) Nomogram. Prediction of patient survival at 2 years and 4 years by three factors: UST2 expression (median UST2 expression as a cut-off value, with 0 representing low expression and 1 representing high expression), age, and staging. (K) Calibration maps. The closer the prediction curve is to the actual curve, the more accurate the prediction is. GO: Gene Ontology; OS: Overall survival; HR: Hazard ratio; KEGG: Kyoto Encyclopedia of Genes and Genomes; TCGA: The Cancer Genome Atlas.
Weighted correlation network analysis (WGCNA) of HER2+ differential genes was done, and four modules were obtained. Among them, the turquoise module was the most correlated with HER2+, with a correlation of 0.94. The turquoise module was a key co-expression module with a gene count of 760 [Figure 1C].
Protein-protein interaction networks (PPIs) were constructed for key genes of key modules and visualized using Cytoscape software. The top 100 key genes were screened using the MCC algorithm of the CytoHubba plugin [Figure 1D].
Gene Ontology (GO) analysis revealed that key genes were mainly enriched in ion transmembrane transport, collagen-containing extracellular matrix, transporter complex, and extracellular matrix structural constituent. The main Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were neuroactive ligand-receptor interaction and cytokine-cytokine receptor interaction [Figure 1E].
The HER2+ differential genes were compared with basal, lumA, and lumB, respectively, to screen key genes that distinguish HER2+ subtypes from other subtypes. As a result, DRD4, UTS2, and GLP1R were obtained as key genes present only in the HER2+ subtype [Figure 1B]. DRD4 and UTS2 were upregulated, while GLP1R was downregulated.
The expression values of DRD4, UTS2, and GLP1R were compared in the four subtypes and the normal samples, respectively. The results showed that the expression of UTS2 was significantly higher in HER2+ subtype than in other subtypes and normal samples. In contrast, the expression values of GLP1R in HER2+ were lower than those in other subtypes and normal samples [Figure 1F].
The PROGgeneV2 (http://genomics.jefferson.edu/proggene/) database was used to see the impact of UTS2, DRD4, and GLP1R expression on the prognosis of HER2+ patients. The data used for survival analysis was GSE6130 from the Gene Expression Omnibus (GEO) database. The clinical features of GSE6130 are shown in Supplementary Table 2, https://links.lww.com/CM9/B268. Results showed that patients with low DRD4 and GLP1R expression had more down survival times than those with high expression, while patients with high UTS2 expression had even lower survival rates [Figure 1G, H]. Similarly, looking at the effect of UTS2, DRD4, and GLP1R on survival in breast cancer patients through the Gene Expression Profiling Interactive Analysis (GEPIA) (http://gepia.cancer-pku.cn/) database yielded the same results as above [Figure 1G, H].
The above results suggest that UTS2 is a key gene associated with HER2+ prognosis. The accuracy of the results was further verified by using the Human Protein Atlas (HPA) database (http://www.proteinatlas.org/) to view its expression in normal breast tissue and breast cancer tissue, which showed that the protein expression level of UTS2 was significantly higher in breast cancer than in normal tissue [Figure 1I].
Finally, clinical information of HER2+ patients, including stage and age, combined with UTS2 expression values were used to construct Cox regression models to predict patient survival at two years and four years. UTS2 was divided into two groups of high and low expression using median values, with zero representing low expression and one representing high expression. The results revealed that patients with low UTS2 expression had a lower score and a better prognosis [Figure 1J]. Also, the accuracy of the outcome prediction showed that the results were reliable [Figure 1K].
All statistical analyses were performed with R software (version 3.6.2, Institute for Statistics and Mathematics, Vienna, Austria;https://www.r-project.org/). Continuous variables are represented using means, and categorical variables are represented using counts. Kaplan–Meier survival rates were compared between groups using the log-rank statistic. The centerpiece of the HPA is its unique antibody collection for mapping the entire human proteome by immunohistochemistry and immunocytochemistry. Multivariate analysis was used to construct a nomogram predicting patient survival. Calibration plots were also used to determine the nomogram's prognostic value. P <0.05 indicated a statistically significant difference.
The above analysis shows that UTS2 is the key gene that distinguishes HER2+ from other subtypes and that high expression of UTS2 leads to poorer patient prognosis. UTS2 (urotensin 2) has three (splice variants), 251 orthologues, and 1 Ensembl protein family member. It is a protein coding gene. Diseases associated with UTS2 include portal hypertension and congestive heart failure. Among its related pathways are RET signaling and signaling by GPCR.[4] Gene Ontology (GO) annotations related to this gene include signaling receptor binding and hormone activity.[5] Chen et al[6] constructed an immunogenetic model of colon cancer using genes such as UTS2 better to predict patient survival at 5 years. Wang et al[7] used multivariate Cox proportional hazards analysis to construct a prognostic model of immunogenetic correlates, including eight genes such as UTS2, which better predicted patient survival. UTS2, an immune-related gene, is highly associated with the prognosis of HER2+ breast cancer and is expected to play an important role in immunotherapy of HER2+.
Conflicts of interest
None.
References
1. Solanki M, Visscher D. Pathology of breast cancer in the last half century. Hum Pathol 2020;95: 137–148. doi: 10.1016/j.humpath. 2019.09.007.
2. Lee A, Moon BI, Kim TH. BRCA1/BRCA2 pathogenic variant breast cancer: Treatment and prevention strategies. Ann Lab Med 2020;40: 114–121. doi: 10.3343/alm.2020.40.2.114.
3. Ross JS, Slodkowska EA, Symmans WF, Pusztai L, Ravdin PM, Hortobagyi GN. The HER-2 receptor and breast cancer: Ten years of targeted anti-HER-2 therapy and personalized medicine. Oncologist 2009;14: 320–368. doi: 10.1634/theoncologist.2008-0230.
4. Ozan B, Demiryürek S, Safdar M, Inanc Y, Demiryürek AT. Lack of association between urotensin-II (UTS2) gene polymorphisms (Thr21Met and Ser89Asn) and migraine. Bosn J Basic Med Sci 2017;17: 268–273. doi: 10.17305/bjbms.2017.2138.
5. Langham RG, Kelly DY. Urotensin II and the kidney. Curr Opin Nephrol Hypertens 2013;22: 107–112. doi: 10.1097/MNH.0b013e32835b6d57.
6. Chen H, Luo J, Guo J. Development and validation of a five-immune gene prognostic risk model in colon cancer. BMC Cancer 2020;20: 395. doi: 10.1186/s12885-020-06799-0.
7. Li F, Guo P, Dong K, Guo P, Wang H, Lv X. Identification of key biomarkers and potential molecular mechanisms in renal cell carcinoma by bioinformatics analysis. J Comput Biol 2019;26: 1278–1295. doi: 10.1089/cmb.2019.0145.