Abstract
PURPOSE: In this study, we aimed at profiling of luminal B breast cancer specific gene expression pattern in Indian women using mRNA-seq and validation based on TCGA expression data. METHODS: RNA isolated from luminal B tumor and adjacent normal tissues was used for library construction and sequencing. Reference-based assemblies of these reads were used for differential gene expression analysis using DeSeq2. The DEGs were evaluated using TCGA expression data. Kaplan-Meier survival method was used to evaluate association between genes showing luminal B specific differential expression pattern and breast cancer prognosis and statistical significance was assessed using log-rank test. Alternate splicing analysis was done using rmats. RESULTS: Differential expression analysis identified 2371 differentially expressed genes (DEGs) in luminal B breast tumors in comparison with adjacent normal tissues of Indian Women. Of them, 1692 DEGs were validated using TCGA luminal B paired samples. Integration of this data with the DEGs obtained by comparative analysis of unpaired luminal B with luminal A unpaired samples from TCGA resulted in 291 DEGs showing luminal B specific expression pattern. Further, 26 genes of prognostic value were identified. Differential splicing analysis between luminal B tumors and adjacent normal tissues in our cohort led to the identification of 687 genes showing significant differential alternate splicing events. CONCLUSION: This study profiled gene expression pattern of luminal B tumors of Indian women and identified 26 key genes of prognostic value for luminal B breast cancer. This study also profiled differential alternate splicing and identified important alternate splicing events in luminal B breast cancer.