Abstract
BACKGROUND: The estrogen receptor (ER)-positive subtype of breast cancer (BC) is the most common type of BC. A number of long noncoding RNAs (lncRNAs) play critical roles in cancer biology, including BC. Previous lncRNA profiling studies have focused only on triple-negative BC and HER 2-positive BC, and no studies have specifically focused on lncRNAs in ER-positive BC. In this study, we analyzed the expression profile of the lncRNAs and mRNAs found in this particular subtype of BC for the first time. METHODS: We evaluated lncRNA microarray data from four pairs of primary BC and adjuvant nontumor breast tissues. Then, we screened out the differently expressed genes and measured the correlation of the expression levels of lncRNAs and ERalpha by Pearson's correlation coefficient analysis. We also performed classification and length distribution of the dysregulated lncRNAs. KEGG pathway analysis was used to understand the biological roles of these differently expressed genes. lncRNA-mRNA coexpression networks were constructed. Finally, RT-PCR was employed to validate the microarray analysis findings. RESULTS: We screened out 2,178 differently expressed lncRNAs, and 13 lncRNAs were found to be associated with the ER expression level. Classification analysis showed that most lncRNAs belonged to intergenic lncRNA and were from 400 to 800 nt in length. Chromosome distribution showed that many of the lncRNAs were mapped to chromosome 1. In the pathway analysis, most of the genes were related to cancer-associated behaviors, such as p53 signaling pathway, cell cycle, focal adhesion, and ECM-receptor interaction. lncRNA-mRNA coexpression networks were constructed, and the lncRNAs related to ESR1, BRCA1, and BRCA2 in the two groups were significantly different. The RT-PCR results were consistent with the data obtained from the microarrays. CONCLUSION: These results provide useful information for exploring potential novel biomarkers as diagnosis and therapy targets for the clinical treatment of ER-positive BC.