Combinatorial epigenetic regulation of non-coding RNAs has profound effects on oncogenic pathways in breast cancer subtypes

N6-methyladenosine (m6A)-circHECA from secondary hair follicle of cashmere goats: identification, regulatory network and expression regulated potentially by methylation of its host gene promoter

OBJECTIVE

The objective of this study was to identify the N6-methyladenosine (m6A)- circHECA molecule in secondary hair follicles (SHFs) of cashmere goats, and generate its potential regulatory network, as well as explore the potential relationship between transcriptional pattern of m6A-circHECA and promoter methylation of its host gene (HECA).

METHODS

The validation of circHECA m6A sites was performed using methylation immunoprecipitation (Me-RIP) along with reverse transcription-quantitative polymerase chain reaction (RT-qPCR) technique. The nucleus and cytoplasm localizations of m6AcircHECA were performed using SHF stem cells of cashmere goats with RT-qPCR analysis. Based on in-silico analysis, the regulatory networks of m6A-circHECA were generated with related signal pathway enrichment. The methylation level of promoter region of m6A-circHECA host gene (HECA) was assessed by the bisulfite sequencing PCR (BSPPCR) technique.

RESULTS

The m6A-circHECA was confirmed to contain four m6A modification sites including m6A-213, m6A-297, m6A-780, and m6A-927, and it was detected mainly in cytoplasm of the SHF stem cells of cashmere goats. The integrated regulatory network analysis showed directly or indirectly complex regulatory relationships between m6A-circHECA of cashmere goats and its potential target molecules: miRNAs, mRNAs, and proteins. The regulatory network and pathway enrichment indicated that m6A-circHECA might play multiple roles in the SHF physiology process of cashmere goats through directly or indirectly interacting or regulating its potential target molecules. A higher methylation level of promoter region of HECA gene in SHFs of cashmere goats might cause the lower expression of m6A-circHECA.

CONCLUSION

The m6A-circHECA might play multiple roles in SHF physiology process of cashmere goats through miRNA mediated pathways along with directly or indirectly interaction with its target proteins. The promoter methylation of m6A-circHECA host gene (HECA) most likely was implicated in its expression inhibition in SHFs of cashmere goats.

Novel histone acetylation-related lncRNA signature for predicting prognosis and tumor microenvironment in esophageal carcinoma

Histone acetylation is one of the most common epigenetic modifications and plays a crucial role in tumorigenesis. However, the prognostic significance of histone acetylation-related lncRNAs (HARlncRNAs) in esophageal carcinoma (ESCA) is not well understood. A total of 653 differentially expressed lncRNAs (DElncRNAs) were identified between 162 ESCA tissues and 11 normal tissues in the TCGA database, and 7 of them were correlated with acetylation regulators. We employed univariate Cox regression analysis, combining it with clinical prognosis information, to select 3 prognostic-related HARlncRNAs for further analysis. Subsequently, we used LASSO regression analysis to construct a risk signature for ESCA and identified C21orf62-AS1 and SSTR5.AS1 as potential biomarkers for the prognosis of ESCA patients. Based on the risk score calculated using the risk signature, we categorized patients into high- and low-risk groups. We identified the risk score as an independent risk factor and validated it in the training, test, and GSE53624 datasets. Additionally, patients categorized by their risk scores exhibited distinct immune statuses, tumor mutation burdens, responses to immunotherapy, and drug sensitivities.

Insights into the role of long non-coding RNAs in DNA methylation mediated transcriptional regulation

DNA methylation is one of the most important epigenetic mechanisms that governing regulation of gene expression, aberrant DNA methylation patterns are strongly associated with human malignancies. Long non-coding RNAs (lncRNAs) have being discovered as a significant regulator on gene expression at the epigenetic level. Emerging evidences have indicated the intricate regulatory effects between lncRNAs and DNA methylation. On one hand, transcription of lncRNAs are controlled by the promoter methylation, which is similar to protein coding genes, on the other hand, lncRNA could interact with enzymes involved in DNA methylation to affect the methylation pattern of downstream genes, thus regulating their expression. In addition, circular RNAs (circRNAs) being an important class of noncoding RNA are also found to participate in this complex regulatory network. In this review, we summarize recent research progress on this crosstalk between lncRNA, circRNA, and DNA methylation as well as their potential functions in complex diseases including cancer. This work reveals a hidden layer for gene transcriptional regulation and enhances our understanding for epigenetics regarding detailed mechanisms on lncRNA regulatory function in human cancers.

Pan-Cancer Methylated Dysregulation of Long Non-coding RNAs Reveals Epigenetic Biomarkers

Different cancer types not only have common characteristics but also have their own characteristics respectively. The mechanism of these specific and common characteristics is still unclear. Pan-cancer analysis can help understand the similarities and differences among cancer types by systematically describing different patterns in cancers and identifying cancer-specific and cancer-common molecular biomarkers. While long non-coding RNAs (lncRNAs) are key cancer modulators, there is still a lack of pan-cancer analysis for lncRNA methylation dysregulation. In this study, we integrated lncRNA methylation, lncRNA expression and mRNA expression data to illuminate specific and common lncRNA methylation patterns in 23 cancer types. Then, we screened aberrantly methylated lncRNAs that negatively regulated lncRNA expression and mapped them to the ceRNA relationship for further validation. 29 lncRNAs were identified as diagnostic biomarkers for their corresponding cancer types, with lncRNA AC027601 was identified as a new KIRC-associated biomarker, and lncRNA ACTA2-AS1 was regarded as a carcinogenic factor of KIRP. Two lncRNAs HOXA-AS2 and AC007228 were identified as pan-cancer biomarkers. In general, the cancer-specific and cancer-common lncRNA biomarkers identified in this study may aid in cancer diagnosis and treatment.

HFS-SLPEE: A Novel Hierarchical Feature Selection and Second Learning Probability Error Ensemble Model for Precision Cancer Diagnosis

The emergence of high-throughput RNA-seq data has offered unprecedented opportunities for cancer diagnosis. However, capturing biological data with highly nonlinear and complex associations by most existing approaches for cancer diagnosis has been challenging. In this study, we propose a novel hierarchical feature selection and second learning probability error ensemble model (named HFS-SLPEE) for precision cancer diagnosis. Specifically, we first integrated protein-coding gene expression profiles, non-coding RNA expression profiles, and DNA methylation data to provide rich information; afterward, we designed a novel hierarchical feature selection method, which takes the CpG-gene biological associations into account and can select a compact set of superior features; next, we used four individual classifiers with significant differences and apparent complementary to build the heterogeneous classifiers; lastly, we developed a second learning probability error ensemble model called SLPEE to thoroughly learn the new data consisting of classifiers-predicted class probability values and the actual label, further realizing the self-correction of the diagnosis errors. Benchmarking comparisons on TCGA showed that HFS-SLPEE performs better than the state-of-the-art approaches. Moreover, we analyzed in-depth 10 groups of selected features and found several novel HFS-SLPEE-predicted epigenomics and epigenetics biomarkers for breast invasive carcinoma (BRCA) (e.g., TSLP and ADAMTS9-AS2), lung adenocarcinoma (LUAD) (e.g., HBA1 and CTB-43E15.1), and kidney renal clear cell carcinoma (KIRC) (e.g., IRX2 and BMPR1B-AS1).

Comprehensive landscape of epigenetic-dysregulated lncRNAs reveals a profound role of enhancers in carcinogenesis in BC subtypes

Aberrant expression of long non-coding RNAs (lncRNA) is associated with altered DNA methylation and histone modifications during carcinogenesis. However, identifying epigenetically dysregulated lncRNAs and characterizing their functional mechanisms in different cancer subtypes are still major challenges for cancer studies. In this study, we systematically analyzed the epigenetic alterations of lncRNAs at important regulatory elements in three breast cancer subtypes. We identified 87, 691, and 1,197 epigenetically dysregulated lncRNAs in luminal, basal, and claudin-low subtypes of breast cancer, respectively. The landscape of epigenetically dysregulated lncRNAs at enhancer elements revealed that epigenetic changes of the majority of lncRNAs occurred in a subtype-specific manner and contributed to subtype-specific biological functions. We identified six acetylation of lysine 27 on histone H3 (H3K27ac)-dysregulated lncRNAs and three DNA methylation-dysregulated lncRNAs (CTC-303L1.2, RP11-738B7.1, and SLC26A4-AS1) as prognostic biomarkers of basal subtype. These lncRNAs were involved in immune response-related biological functions. Treatment of the basal breast cancer cell line MDA-MB-468 with CREBBP/EP300 bromodomain inhibitors downregulated H3K27 acetylation levels and caused a decrease in the expression of five H3K27ac-dysregulated lncRNAs (LINC00393, KB-1836B5.1, RP1-140K8.5, AC005162.1, and AC020916.2) and inhibition of the growth of breast cancer cells. One epigenetically dysregulated lncRNA (LINC01983) and four lncRNA regulators (UCA1, RP11-221J22.2, RP11-221J22.1, and RP1-212P9.3) were identified as prognostic biomarkers of the luminal molecular subtype of breast cancer by controlling the tumor necrosis factor (TNF) signaling pathway, T helper (Th)17 cell differentiation, and T cell migration. Finally, our results highlighted a profound role of enhancer-related H3K27ac-dysregulated lncRNAs, DNA methylation-dysregulated lncRNAs, and lncRNA regulators in breast cancer subtype carcinogenesis and their potential prognostic value.

Revealing Epigenetic Factors of circRNA Expression by Machine Learning in Various Cellular Contexts

Circular RNAs (circRNAs) have been identified as naturally occurring RNAs that are highly represented in the eukaryotic transcriptome. Although a large number of circRNAs have been reported, the underlying regulatory mechanism of circRNAs biogenesis remains largely unknown. Here, we integrated in-depth multi-omics data including epigenome, transcriptome, and non-coding RNA and identified candidate circRNAs in six cellular contexts. Next, circRNAs were divided into two classes (high versus low) with different expression levels. Machine learning models were constructed that predicted circRNA expression levels based on 11 different histone modifications and host gene expression. We found that the models achieve great accuracy in predicting high versus low expressed circRNAs. Furthermore, the expression levels of host genes of circRNAs, H3k36me3, H3k79me2, and H4k20me1 contributed greatly to the classification models in six cellular contexts. In summary, all these results suggest that epigenetic modifications, particularly histone modifications, can effectively predict expression levels of circRNAs.

•

CircRNAs exhibit specific expression in various cellular contexts

•

High and low expressed circRNAs exhibit different biological functions

•

Histone modifications are significantly correlated with circRNAs expression

•

Machine learning models were constructed for predicting circRNAs expression

Blocking circ_0000520 Suppressed Breast Cancer Cell Growth, Migration and Invasion Partially via miR-1296/SP1 Axis Both in vitro and in vivo

Background

Breast cancer (BCa) is an overwhelming malignant tumor mainly in women globally. Circular RNAs (circRNAs) are a special type of noncoding RNAs involved in competing endogenous RNA (ceRNA) network, a classic molecular mechanism of the tumorigenesis of human cancers, including BCa. Here, we intended to explore the role and mechanism of hsa_circ_0000520 (circ_0000520) in BCa cells.

Methods

Expression of circ_0000520, miRNA-1296-5p (miR-1296) and specificity protein 1 (SP1) was measured by real time-quantitative PCR and Western blotting. Cell growth was measured by cell counting kit-8, colony formation assay and flow cytometry method. Cell migration and invasion were assessed by transwell assays and Western blotting. Tumor growth was determined by xenograft models. The direct interaction among circ_0000520, miR-1296 and SP1 was confirmed by dual-luciferase reporter assay and RNA pull-down assay.

Results

circ_0000520 was upregulated in BCa tumors and cell lines (T47D, MCF7, MDA-MB-231, BT549, and SKBR3), and circ_0000520 high expression was associated with poor overall survival. Blocking circ_0000520 suppressed cell viability, colony formation, migration and invasion, but promoted cell cycle arrest and apoptosis rate in MDA-MB-231 and MCF7 cells. circ_0000520 could directly regulate miR-1296 expression, and SP1 was a novel target for miR-1296. Moreover, the anti-tumor role of circ_0000520 silencing was abrogated by miR-1296 downregulation or SP1 restoration. Notably, tumor growth of MDA-MB-231 cells in mice was restrained by circ_0000520 deletion.

Conclusion

circ_0000520 knockdown could suppress cell growth, migration and invasion both in vitro and in vivo through regulating miR-1296/SP1 pathway.

Long Non-Coding RNA NEAT1 Promotes the Proliferation, Migration, and Metastasis of Human Breast-Cancer Cells by Inhibiting miR-146b-5p Expression

Background

Breast cancer (BC) is the most commonly diagnosed cancer in women. Tumor recurrence and metastasis are the key causes of death in BC patients. Long non-coding RNA (lncRNA) is closely associated with BC progression. lncRNA nuclear-enriched abundant transcript (NEAT)1 has been reported to regulate the proliferation and mobility of several types of cancer cells. However, how lncRNA NEAT1 affects the proliferation and invasion of BC cells is not known.

Methods

Quantitative real time-polymerase chain reaction (qRT-PCR) was used to measure expression of lncRNA NEAT1 and microRNA (miR)-146b-5p in BC tissues and cell lines. Cell Counting Kit (CCK)-8, cell colony-formation, wound-healing, and Transwell™ assays were undertaken to determine the effects of lncRNA NEAT1 and miR-146b-5p on progression of BC cells. The interaction between lncRNA NEAT1 and miR-146b-5p was examined by luciferase reporter, RNA-binding protein immunoprecipitation (RIP), and RNA-pulldown assays.

Results

Expression of lncRNA NEAT1 was upregulated in BC tissues and cell lines. High expression of lncRNA NEAT1 predicted poor overall survival in BC patients. Silencing of expression of lncRNA NEAT1 inhibited epithelial–mesenchymal transition (EMT) and suppressed the proliferation, migration and invasion of BC cells. Ectopic expression of lncRNA NEAT1 induced EMT and promoted BC progression. Mechanistic investigations revealed that miR-146b-5p was a direct target of lncRNA NEAT1, and its expression was correlated negatively with expression of lncRNA NEAT1 in BC tissues.

Conclusion

lncRNA NEAT1 could (i) serve as a novel prognostic marker for BC and (ii) be a potential therapeutic target for BC.

Systematic review regulatory principles of non-coding RNAs in cardiovascular diseases

Cardiovascular diseases (CVDs) continue to be a major cause of morbidity and mortality, and non-coding RNAs (ncRNAs) play critical roles in CVDs. With the recent emergence of high-throughput technologies, including small RNA sequencing, investigations of CVDs have been transformed from candidate-based studies into genome-wide undertakings, and a number of ncRNAs in CVDs were discovered in various studies. A comprehensive review of these ncRNAs would be highly valuable for researchers to get a complete picture of the ncRNAs in CVD. To address these knowledge gaps and clinical needs, in this review, we first discussed dysregulated ncRNAs and their critical roles in cardiovascular development and related diseases. Moreover, we reviewed >28 561 published papers and documented the ncRNA-CVD association benchmarking data sets to summarize the principles of ncRNA regulation in CVDs. This data set included 13 249 curated relationships between 9503 ncRNAs and 139 CVDs in 12 species. Based on this comprehensive resource, we summarized the regulatory principles of dysregulated ncRNAs in CVDs, including the complex associations between ncRNA and CVDs, tissue specificity and ncRNA synergistic regulation. The highlighted principles are that CVD microRNAs (miRNAs) are highly expressed in heart tissue and that they play central roles in miRNA-miRNA functional synergistic network. In addition, CVD-related miRNAs are close to one another in the functional network, indicating the modular characteristic features of CVD miRNAs. We believe that the regulatory principles summarized here will further contribute to our understanding of ncRNA function and dysregulation mechanisms in CVDs.

LncMAP: Pan-cancer atlas of long noncoding RNA-mediated transcriptional network perturbations

Gene regulatory network perturbations contribute to the development and progression of cancer, however, molecular determinants that mediate transcriptional perturbations remain a fundamental challenge for cancer biology. We show that transcriptional perturbations are widely mediated by long noncoding RNAs (lncRNAs) via integration of genome-wide transcriptional regulation with paired lncRNA and gene expression profiles. Systematic construction of an LncRNA Modulator Atlas in Pan-cancer (LncMAP) reveals distinct types of lncRNA regulatory molecules, which are expressed in multiple tissues, exhibit higher conservation. Strikingly, cancers with similar tissue origin share lncRNA modulators which perturb the regulation of cell cycle and immune response-related functions. Furthermore, we identified a large number of pan-cancer lncRNA modulators with potential clinical significance, which are differentially expressed in cancer or are strongly correlated with drug sensitivity across cell lines. Further stratification of cancer patients based on lncRNA-mediated transcriptional perturbations identifies subtypes with distinct survival rates. Finally, we made a user-friendly web interface available for exploring lncRNA-mediated transcriptional perturbations across cancer types. Our study provides a systems-level dissection of lncRNA-mediated regulatory perturbations in cancer, and also presents a valuable tool and resource for investigating the function of lncRNAs in cancer.

Long non-coding RNA LINC00968 reduces cell proliferation and migration and angiogenesis in breast cancer through up-regulation of PROX1 by reducing hsa-miR-423-5p

Background: Breast cancer (BC) is a common invasive malignancy in women with unclear etiology. A recent study suggested that long non-coding RNA (lncRNA), LINC00968 had a tumor-promoting effect in cancer. However, the role of LINC00968 in BC remains unclear. Therefore, we conducted the present study to determine the effect of LINC00968 in BC and its underlying mechanism. Methods: The expression of LINC00968 and hsa-miR-423-5p in BC tissues and cells was determined using reverse transcription quantitative polymerase chain reaction and western blot analysis. Dual luciferase reporter, RNA pull-down and RNA immunoprecipitation assays were used to determine the relationship among LINC00968, PROX1 and hsa-miR-423-5p. Gain- and loss-function approaches were utilized to examine the effects of LINC00968, PROX1 and hsa-miR-423-5p on cell proliferation, migration, tube formation in vitro; and tumor growth and angiogenesis in vivo. Results: LINC00968 expression reduced while hsa-miR-423-5p increased in BC tissues relative to adjacent normal tissues. Overexpression of LINC00968 was observed to inhibit BC cell proliferation, migration and tube formation abilities in vitro as well as tumor growth in vivo through inhibition of hsa-miR-423-5p. And hsa-miR-423-5p mediated BC cellular functions and tumor growth through down-regulating PROX1. LINC00968 was identified as a competing endogenous RNA to upregulate PROX1 by downregulating hsa-miR-423-5p. More importantly, it was found that LINC00968 increased PROX1 expression in vivo in a concentration-dependent manner. Conclusion: Taken together, this study suggests that LINC00968 inhibits the progression of BC through impeding hsa-miR-423-5p-mediated PROX1 inhibition. LINC00968 may be a potential therapeutic target for BC therapy that warrants further studies.

Targeting microRNA/UHRF1 pathways as a novel strategy for cancer therapy

Ubiquitin-like containing plant homeodomain and RING finger domains 1 (UHRF1) is an anti-apoptotic protein involved in the silencing of several tumor suppressor genes (TSGs) through epigenetic modifications including DNA methylation and histone post-translational alterations, and also epigenetic-independent mechanisms. UHRF1 overexpression is observed in a number of solid tumors and hematological malignancies, and is considered a primary mechanism in inhibiting apoptosis. UHRF1 exerts its inhibitory activity on TSGs by binding to functional domains and therefore influences several epigenetic actors including DNA methyltransferase, histone deacetylase 1, histone acetyltransferase Tat-interacting protein 60 and histone methyltransferases G9a and Suv39H1. UHRF1 is considered to control a large macromolecular protein complex termed epigenetic code replication machinery, in order to maintain epigenetic silencing of TSGs during cell division, thus enabling cancer cells to escape apoptosis. MicroRNAs (miRNAs) are able to regulate the expression of its target gene by functioning as either an oncogene or a tumor suppressor. In the present review, the role of tumor suppressive miRNAs in the regulation of UHRF1, and the importance of targeting the microRNA/UHRF1 pathways in order to induce the reactivation of silenced TSGs and subsequent apoptosis are discussed.