Genetics and Expression Profile of the Tubulin Gene Superfamily in Breast Cancer Subtypes and Its Relation to Taxane Resistance

Taxanes are a class of chemotherapeutic agents that inhibit cell division by disrupting the mitotic spindle through the stabilization of microtubules. Most breast cancer (BC) tumors show resistance against taxanes partially due to alterations in tubulin genes. In this project we investigated tubulin isoforms in BC to explore any correlation between tubulin alterations and taxane resistance. Genetic alteration and expression profiling of 28 tubulin isoforms in 6714 BC tumor samples from 4205 BC cases were analyzed. Protein-protein, drug-protein and alterations neighbor genes in tubulin pathways were examined in the tumor samples. To study correlation between promoter activity and expression of the tubulin isoforms in BC, we analyzed the ChIP-seq enrichment of active promoter histone mark H3K4me3 and mRNA expression profile of MCF-7, ZR-75-30, SKBR-3 and MDA-MB-231 cell lines. Potential correlation between tubulin alterations and taxane resistance, were investigated by studying the expression profile of taxane-sensitive and resistant BC tumors also the MDA-MB-231 cells acquired resistance to paclitaxel. All genomic data were obtained from public databases. Results showed that TUBD1 and TUBB3 were the most frequently amplified and deleted tubulin genes in the BC tumors respectively. The interaction analysis showed physical interactions of alpha-, beta- and gamma-tubulin isoforms with each other. The most of FDA-approved tubulin inhibitor drugs including taxanes target only beta-tubulins. The analysis also revealed sex tubulin-interacting neighbor proteins including ENCCT3, NEK2, PFDN2, PTP4A3, SDCCAG8 and TBCE which were altered in at least 20% of the tumors. Three of them are tubulin-specific chaperons responsible for tubulin protein folding. Expression of tubulin genes in BC cell lines were correlated with H3K4me3 enrichment on their promoter chromatin. Analyzing expression profile of BC tumors and tumor-adjacent normal breast tissues showed upregulation of TUBA1A, TUBA1C, TUBB and TUBB3 and downregulation of TUBB2A, TUBB2B, TUBB6, TUBB7P pseudogene, and TUBGCP2 in the tumor tissues compared to the normal breast tissues. Analyzing taxane-sensitive versus taxane-resistant tumors revealed that expression of TUBB3 and TUBB6 was significantly downregulated in the taxane-resistant tumors. Our results suggest that downregulation of tumor beta III- and beta V-tubulins is correlated with taxane resistance in BC. Based on our results, we conclude that aberrant protein folding of tubulins due to mutation and/or dysfunction of tubulin-specific chaperons may be potential mechanisms of taxane resistance. Thus, we propose studying the molecular pathology of tubulin mutations and folding in BC and their impacts on taxane resistance.