Organotypic three-dimensional (3D) cell cultures more accurately mimic the characteristics of solid tumors in vivo in comparison with traditional two-dimensional (2D) monolayer cell models. Currently, studies on the regulation of long non-coding RNAs (lncRNAs) have not been explored in breast cancer cells cultured in 3D microenvironments. In the present research, we studied the expression and potential roles of lncRNAs in estrogen receptor-positive luminal B subtype BT-474 breast cancer cells grown over extracellular matrix proteins-enriched 3D cultures. Global expression profiling using DNA microarrays identifies 290 upregulated and 183 downregulated lncRNAs in 3D cultures relative to 2D condition. Using a co-expression analysis approach of lncRNAs and mRNAs pairs expressed in the same experimental conditions, we identify hundreds of regulatory axes modulating genes involved in cancer hallmarks, such as responses to estrogens, cell proliferation, hypoxia, apical junctions, and resistance to endocrine therapy. In addition, we identified 102 lncRNAs/mRNA correlations in 3D cultures, which were similar to those reported in TCGA datasets obtained from luminal B breast cancer patients. Interestingly, we also found a set of mRNAs transcripts co-expressed with LINC00847 and CTD-2566J3.1 lncRNAs, which were predictors of pathologic complete response and overall survival. Finally, both LINC00847 and CTD -2566J3.1 were co-expressed with essential genes for cancer genetic dependencies, such as FOXA1 y GINS2. Our experimental and predictive findings show that co-expressed lncRNAs/mRNAs pairs exhibit a high degree of similarity with those found in luminal B breast cancer patients, suggesting that they could be adequate pre-clinical tools to identify not only biomarkers related to endocrine therapy response and PCR, but to understand the biological behavior of cancer cells in 3D microenvironments.
Keywords: breast cancer; co-regulation networks; lncRNAs; organotypic 3D cultures; therapy response.