Comprehensive transcriptomic analysis of heat shock proteins in the molecular subtypes of human breast cancer

Heat Shock Proteins (HSPs), a family of genes with key roles in proteostasis, have been extensively associated with cancer behaviour. However, the HSP family is quite large and many of its members have not been investigated in breast cancer (BRCA), particularly in relation with the current molecular BRCA classification. In this work, we performed a comprehensive transcriptomic study of the HSP gene family in BRCA patients from both The Cancer Genome Atlas (TCGA) and the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) cohorts discriminating the BRCA intrinsic molecular subtypes.

We examined gene expression levels of 1097 BRCA tissue samples retrieved from TCGA and 1981 samples of METABRIC, focusing mainly on the HSP family (95 genes). Data were stratified according to the PAM50 gene expression (Luminal A, Luminal B, HER2, Basal, and Normal-like). Transcriptomic analyses include several statistical approaches: differential gene expression, hierarchical clustering and survival analysis.

Of the 20,531 analysed genes we found that in BRCA almost 30% presented deregulated expression (19% upregulated and 10% downregulated), while of the HSP family 25% appeared deregulated (14% upregulated and 11% downregulated) (|fold change| > 2 comparing BRCA with normal breast tissues). The study revealed the existence of shared HSP genes deregulated in all subtypes of BRCA while other HSPs were deregulated in specific subtypes. Many members of the Chaperonin subfamily were found upregulated while three members (BBS10, BBS12 and CCTB6) were found downregulated. HSPC subfamily had moderate increments of transcripts levels. Various genes of the HSP70 subfamily were upregulated; meanwhile, HSPA12A and HSPA12B appeared strongly downregulated. The strongest downregulation was observed in several HSPB members except for HSPB1. DNAJ members showed heterogeneous expression pattern. We found that 23 HSP genes correlated with overall survival and three HSP-based transcriptional profiles with impact on disease outcome were recognized.

We identified shared and specific HSP genes deregulated in BRCA subtypes. This study allowed the recognition of HSP genes not previously associated with BRCA and/or any cancer type, and the identification of three clinically relevant clusters based on HSPs expression patterns with influence on overall survival.