Expression profiling of luminal B breast tumor in Indian women

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.