LY6K-AS lncRNA is a lung adenocarcinoma prognostic biomarker and regulator of mitotic progression

LINC00460-FUS-MYC feedback loop drives breast cancer metastasis and doxorubicin resistance

Therapeutic resistance and metastasis largely contribute to mortality from breast cancer and therefore understanding the underlying mechanisms of such remains an urgent challenge. By cross-analysis of TCGA and GEO databases, LINC00460 was identified as an oncogenic long non-coding RNA, highly expressed in Doxorubicin resistant breast cancer. LINC00460 was further demonstrated to promote stem cell-like and epithelial-mesenchymal transition (EMT) characteristics in breast cancer cells. LINC00460 interacts with FUS protein with consequent enhanced stabilization, which further promotes MYC mRNA maturation. LINC00460 expression was transcriptionally enhanced by c-MYC protein, forming a positive feedback loop to promote metastasis and Doxorubicin resistance. LINC00460 depletion in Doxorubicin-resistant breast cancer cells restored sensitivity to Doxorubicin and increased the efficacy of c-MYC inhibitor therapy. Collectively, these findings implicate LINC00460 as a promising prognostic biomarker and potential therapeutic target to overcome Doxorubicin resistance in breast cancer.

Endogenous retroviral solo-LTRs in human genome

Human endogenous retroviruses (HERVs) are derived from the infection and integration of exogenetic retroviruses. HERVs account for 8% of human genome, and the majority of HERVs are solitary LTRs (solo-LTRs) due to homologous recombination. Multiple findings have showed that solo-LTRs could provide an enormous reservoir of transcriptional regulatory sequences involved in diverse biological processes, especially carcinogenesis and cancer development. The link between solo-LTRs and human diseases still remains poorly understood. This review focuses on the regulatory modules of solo-LTRs, which contribute greatly to the diversification and evolution of human genes. More importantly, although inactivating mutations, insertions and deletions have been identified in solo-LTRs, the inherited regulatory elements of solo-LTRs initiate the expression of chimeric lncRNA transcripts, which have been reported to play crucial roles in human health and disease. These findings provide valuable insights into the evolutionary and functional mechanisms underlying the presence of HERVs in human genome. Taken together, in this review, we will present evidences showing the regulatory and encoding capacity of solo-LTRs as well as the significant impact on various aspects of human biology.

The pancancer overexpressed NFYC Antisense 1 controls cell cycle mitotic progression through in cis and in trans modes of action

Antisense RNAs (asRNAs) represent an underappreciated yet crucial layer of gene expression regulation. Generally thought to modulate their sense genes in cis through sequence complementarity or their act of transcription, asRNAs can also regulate different molecular targets in trans, in the nucleus or in the cytoplasm. Here, we performed an in-depth molecular characterization of NFYC Antisense 1 (NFYC-AS1), the asRNA transcribed head-to-head to NFYC subunit of the proliferation-associated NF-Y transcription factor. Our results show that NFYC-AS1 is a prevalently nuclear asRNA peaking early in the cell cycle. Comparative genomics suggests a narrow phylogenetic distribution, with a probable origin in the common ancestor of mammalian lineages. NFYC-AS1 is overexpressed pancancer, preferentially in association with RB1 mutations. Knockdown of NFYC-AS1 by antisense oligonucleotides impairs cell growth in lung squamous cell carcinoma and small cell lung cancer cells, a phenotype recapitulated by CRISPR/Cas9-deletion of its transcription start site. Surprisingly, expression of the sense gene is affected only when endogenous transcription of NFYC-AS1 is manipulated. This suggests that regulation of cell proliferation is at least in part independent of the in cis transcription-mediated effect on NFYC and is possibly exerted by RNA-dependent in trans effects converging on the regulation of G2/M cell cycle phase genes. Accordingly, NFYC-AS1-depleted cells are stuck in mitosis, indicating defects in mitotic progression. Overall, NFYC-AS1 emerged as a cell cycle-regulating asRNA with dual action, holding therapeutic potential in different cancer types, including the very aggressive RB1-mutated tumors.

The long non-coding RNA LINC00707 interacts with Smad proteins to regulate TGFβ signaling and cancer cell invasion

BACKGROUND

Long non-coding RNAs (lncRNAs) regulate cellular processes by interacting with RNAs or proteins. Transforming growth factor β (TGFβ) signaling via Smad proteins regulates gene networks that control diverse biological processes, including cancer cell migration. LncRNAs have emerged as TGFβ targets, yet, their mechanism of action and biological role in cancer remain poorly understood.

METHODS

Whole-genome transcriptomics identified lncRNA genes regulated by TGFβ. Protein kinase inhibitors and RNA-silencing, in combination with cDNA cloning, provided loss- and gain-of-function analyses. Cancer cell-based assays coupled to RNA-immunoprecipitation, chromatin isolation by RNA purification and protein screening sought mechanistic evidence. Functional validation of TGFβ-regulated lncRNAs was based on new transcriptomics and by combining RNAscope with immunohistochemical analysis in tumor tissue.

RESULTS

Transcriptomics of TGFβ signaling responses revealed down-regulation of the predominantly cytoplasmic long intergenic non-protein coding RNA 707 (LINC00707). Expression of LINC00707 required Smad and mitogen-activated protein kinase inputs. By limiting the binding of Krüppel-like factor 6 to the LINC00707 promoter, TGFβ led to LINC00707 repression. Functionally, LINC00707 suppressed cancer cell invasion, as well as key fibrogenic and pro-mesenchymal responses to TGFβ, as also attested by RNA-sequencing analysis. LINC00707 also suppressed Smad-dependent signaling. Mechanistically, LINC00707 interacted with and retained Smad proteins in the cytoplasm. Upon TGFβ stimulation, LINC00707 dissociated from the Smad complex, which allowed Smad accumulation in the nucleus. In vivo, LINC00707 expression was negatively correlated with Smad2 activation in tumor tissues.

CONCLUSIONS

LINC00707 interacts with Smad proteins and limits the output of TGFβ signaling, which decreases LINC00707 expression, thus favoring cancer cell invasion. Video Abstract.

Exploring tumor-normal cross-talk with TranNet: Role of the environment in tumor progression

There is a growing awareness that tumor-adjacent normal tissues used as control samples in cancer studies do not represent fully healthy tissues. Instead, they are intermediates between healthy tissues and tumors. The factors that contribute to the deviation of such control samples from healthy state include exposure to the tumor-promoting factors, tumor-related immune response, and other aspects of tumor microenvironment. Characterizing the relation between gene expression of tumor-adjacent control samples and tumors is fundamental for understanding roles of microenvironment in tumor initiation and progression, as well as for identification of diagnostic and prognostic biomarkers for cancers. To address the demand, we developed and validated TranNet, a computational approach that utilizes gene expression in matched control and tumor samples to study the relation between their gene expression profiles. TranNet infers a sparse weighted bipartite graph from gene expression profiles of matched control samples to tumors. The results allow us to identify predictors (potential regulators) of this transition. To our knowledge, TranNet is the first computational method to infer such dependencies. We applied TranNet to the data of several cancer types and their matched control samples from The Cancer Genome Atlas (TCGA). Many predictors identified by TranNet are genes associated with regulation by the tumor microenvironment as they are enriched in G-protein coupled receptor signaling, cell-to-cell communication, immune processes, and cell adhesion. Correspondingly, targets of inferred predictors are enriched in pathways related to tissue remodelling (including the epithelial-mesenchymal Transition (EMT)), immune response, and cell proliferation. This implies that the predictors are markers and potential stromal facilitators of tumor progression. Our results provide new insights into the relationships between tumor adjacent control sample, tumor and the tumor environment. Moreover, the set of predictors identified by TranNet will provide a valuable resource for future investigations.

Identification of New Prognostic Genes and Construction of a Prognostic Model for Lung Adenocarcinoma

Lung adenocarcinoma (LUAD) is a rapidly progressive malignancy, and its mortality rate is very high. In this study, we aimed at finding novel prognosis-related genes and constructing a credible prognostic model to improve the prediction for LUAD patients. Differential gene expression, mutant subtype, and univariate Cox regression analyses were conducted with the dataset from the Cancer Genome Atlas (TCGA) database to screen for prognostic features. These features were employed in the following multivariate Cox regression analysis and the produced prognostic model included the stage and expression of SMCO2, SATB2, HAVCR1, GRIA1, and GALNT4, as well as mutation subtypes of TP53. The exactness of the model was confirmed by an overall survival (OS) analysis and disease-free survival (DFS) analysis, which indicated that patients in the high-risk group had a poorer prognosis compared to those in the low-risk group. The area under the receiver operating characteristic curve (AUC) was 0.793 in the training group and 0.779 in the testing group. The AUC of tumor recurrence was 0.778 in the training group and 0.815 in the testing group. In addition, the number of deceased patients increased as the risk scores raised. Furthermore, the knockdown of prognostic gene HAVCR1 suppressed the proliferation of A549 cells, which supports our prognostic model that the high expression of HAVCR1 predicts poor prognosis. Our work created a reliable prognostic risk score model for LUAD and provided potential prognostic biomarkers.

H3K4me3-related lncRNAs signature and comprehensive analysis of H3K4me3 regulating tumor immunity in lung adenocarcinoma

BACKGROUD

The role of epigenetic modifications in tumorigenesis has been widely reported. However, the role and mechanism of H3K4me3 modification in lung adenocarcinoma (LUAD) are rarely reported systematically. We, therefore, sought to analyze the characteristics of LUAD associated with H3K4me3 modification, build an H3K4me3-lncRNAs score model to predict the prognosis of patients with LUAD and clarify the potential value of H3K4me3 in immunotherapy of LUAD.

METHODS

We evaluated H3K4me3-lncRNA patterns and H3K4me3-lncRNA scores of 477 LUAD samples based on 53 lncRNAs closely correlated to H3K4me3 regulators and comprehensive analyzed the role of these patterns in tumorigenesis and tumor immunity. Using Gene set variation analysis (GSVA), we systematically evaluated the H3K4me3 level of every sample and deeply analyzed the effect of H3K4me3 on the prognosis of LUAD. In addition, we included two independent immunotherapy cohorts to study the impact of high H3K4me3 score on the prognosis of patients. We also used an independent cohort with 52 matched paraffin specimens of LUAD to verify the impact of high H3K3me3 expression on the prognosis of patients.

RESULTS

We identified three H3K4me3-lncRNA patterns that exhibited specific immune characteristics. Characterized by immunosuppressive and increased TGFβ-mediated epithelial-mesenchymal transition (EMT), patients with high H3K4me3-lncRNA score had a poor overall survival and decreased H3K4me3 score. H3K4me3 score was significantly positively correlated with CD4+T-cell and CD8+T-cell activation, programmed cell death and immune checkpoints (ICs) expression, and was negatively correlated with MYC pathway, TP53 pathway, and cell proliferation. Patients with high H3K4me3 score showed elevated expression of ICs, potentiated CD4 T-cell and CD8 T-cell activation, increased programmed cell death, and suppressed cell proliferation and TGFβ-mediated EMT. Patients with high H3K4me3 score and high expression of CTLA4, ICOS, TIGIT, PDCD1LG2, IDO1, CD274, PDCD1, LAG3, or HAVCR2 had the best survival advantage. Two independent immunotherapy cohorts verified that patients with high H3K4me3 score showed an increased inflamed tumor microenvironment (TME) phenotype and enhanced anti-PD-1/L1 immunotherapy response. Immunohistochemistry (IHC) data from 52 matched paraffin specimens of LUAD confirmed that the protein level of H3K4me3 in tumor was significantly lower than that of paracancerous tissues and H3K4me3 brought significant survival benefits to patients with LUAD.

CONCLUSIONS

We build an H3K4me3-lncRNAs score model to predict the prognosis of patients with LUAD. More importantly, this study revealed characteristics of H3K4me3 modification in LUAD and clarified the important potential role of H3K4me3 on tumor immunotherapy and patients' survival.

LINC00668 promoted non-small lung cancer progression by miR-518c-3p/TRIP4 axis

BACKGROUND

Non-small lung cancer ranks first in the cancer-related death of all malignant tumors. Exploring novel biological targets is of great significance for diagnosis and therapy of NSCLC.

OBJECTIVE

In this study, we aimed to explore the effect of LINC00668 on the biological functions of NSCLC cells and the underlying mechanism.

METHODS

RT-qPCR assays and western blot assays were utilized to estimate the relative gene expression at mRNA and protein levels, respectively. CCK8, colony formation, wound healing, transwell, and cell apoptosis assays were employed to assess cell function. IHC and FISH assays were used to determine the gene expression in NSCLC tissues. RIP and dual-luciferase assays were conducted to validate the combination between LINC00668 and miR-518c-3p. The correlation of expression between miR-518c-3p and LINC00668 or TRIP4 was determined by Pearson correlation analysis.

RESULTS

LINC00668 was aberrantly upregulated in NSCLC tumor tissues and cell lines. Inhibition of LINC00668 significantly suppressed tumor proliferation, migration, invasion and promoted cell apoptosis. Mechanistically, LINC00668 could bind to miR-518c-3p, thus targeting the 3'UTR of TRIP4. TRIP4 overexpression rescued the weakened cell function mediated by LINC00668 silencing.

CONCLUSIONS

LINC00668 acted as an oncogene in NSCLC progression through miR-518c-3p/TRIP4 axis. Our study disclosed a new mechanism of LINC00668 functioned in NSCLC and may give a deeper insight of the targeted therapy of NSCLC in the future.

Gene amplification-driven lncRNA SNHG6 promotes tumorigenesis via epigenetically suppressing p27 expression and regulating cell cycle in non-small cell lung cancer

Long non-coding RNAs (lncRNAs) have been validated to play essential roles in non-small cell lung carcinoma (NSCLC) progression. In this study, through systematically screening GSE33532 and GSE29249 from Gene Expression Omnibus (GEO) database and bioinformatics analysis, we found the significant upregulation of SNHG6 in NSCLC. The activation of SNHG6 was driven by copy number amplification and high expression of SNHG6 indicated a poor prognosis. Functionally, the knockdown of SNHG6 inhibited NSCLC cell proliferation, migration, and suppressed the G1/S transition of the cell cycle. SNHG6 overexpression had the opposite effects. Mechanically, SNHG6 recruited EZH2 to the promoter region of p27 and increased H3K27me3 enrichment, thus epigenetically repressing the expression of p27, regulating the cell cycle, and promoting tumorigenesis of NSCLC. SNHG6 silencing restrained tumor growth in vivo and suppressed the expressions of cell cycle-related proteins in the G1/S transition. In conclusion, our study uncovered a novel mechanism of SNHG6 activation and its function. SNHG6 can be considered a potential target for the diagnosis and treatment of NSCLC in the future.

A prognostic and therapeutic hallmark developed by the integrated profile of basement membrane and immune infiltrative landscape in lung adenocarcinoma

Basement membranes (BMs) are specialised extracellular matrices that maintain cellular integrity and resist the breaching of carcinoma cells for metastases while regulating tumour immunity. The tumour immune microenvironment (TME) is essential for tumour growth and the response to and benefits from immunotherapy. In this study, the BM score and TME score were constructed based on the expression signatures of BM-related genes and the presence of immune cells in lung adenocarcinoma (LUAD), respectively. Subsequently, the BM-TME classifier was developed with the combination of BM score and TME score for accurate prognostic prediction. Further, Kaplan-Meier survival estimation, univariate Cox regression analysis and receiver operating characteristic curves were used to cross-validate and elucidate the prognostic prediction value of the BM-TME classifier in several cohorts. Findings from functional annotation analysis suggested that the potential molecular regulatory mechanisms of the BM-TME classifier were closely related to the cell cycle, mitosis and DNA replication pathways. Additionally, the guiding value of the treatment strategy of the BM-TME classifier for LUAD was determined. Future clinical disease management may benefit from the findings of our research.

Z-DNA-Containing Long Terminal Repeats of Human Endogenous Retrovirus Families Provide Alternative Promoters for Human Functional Genes

Transposable elements (TEs) account for approximately 45% of the human genome. TEs have proliferated randomly and integrated into functional genes during hominoid radiation. They appear as right-handed B-DNA double helices and slightly elongated left-handed Z-DNAs. Human endogenous retrovirus (HERV) families are widely distributed in human chromosomes at a ratio of 8%. They contain a 5'-long terminal repeat (LTR)-gag-pol-env-3'-LTR structure. LTRs contain the U3 enhancer and promoter region, transcribed R region, and U5 region. LTRs can influence host gene expression by acting as regulatory elements. In this review, we describe the alternative promoters derived from LTR elements that overlap Z-DNA by comparing Z-hunt and DeepZ data for human functional genes. We also present evidence showing the regulatory activity of LTR elements containing Z-DNA in GSDML. Taken together, the regulatory activity of LTR elements with Z-DNA allows us to understand gene function in relation to various human diseases.

Transmembrane Protein-Based Risk Model and H3K4me3 Modification Characteristics in Lung Adenocarcinoma

The role and mechanism of transmembrane proteins (TMEMs) in tumorigenesis remain unclear. Based on 4 independent cohorts containing 1,208 cases, we identified 3 TMEMs (TMEM273, TMEM164, and TMEM125), which were used to construct a risk model to predict the prognosis of LUAD. The two patterns based on the risk score exhibited a high degree of consistency with the characteristics of immune cell infiltration and epigenetic distribution. Patients with a low-risk score, characterized by an increased activation of immunity, H3K4me3 modification, tumor cell apoptosis, chemokine secretion, and TMB, had better disease-free survival (DFS) and overall survival (OS). Obvious immunosuppression, increased epithelial–mesenchymal transition, a low H3K4me3 level, shortened cell cycle, and accelerated cell division manifested in high-risk patients, with poorer DFS and OS. The model showed a better prognostic value than the tumor immune dysfunction and exclusion score. Correlation analysis told us that patients with high scores were suitable for treatment with CD276 inhibitors for their higher levels of CD276 expression. The risk score had a strong negative correlation with HAVCR2 and ICOS among patients with EGFR-WT, KRAS-WT, STK11-WT, or TP53-MUT, and patients with these mutation types with low scores were suitable for treatment with HAVCR2 or ICOS inhibitors. This work comprehensively analyzed the role and mechanism of TMEMs in LUAD and revealed the characteristics of histone methylation modification. The TMEM-based signature gave us deep insight into immune cell infiltration profiles and provided an individualized immunotherapy strategy.

Overexpression of lncRNA A2M-AS1 inhibits cell growth and aggressiveness via regulating the miR-587/bone morphogenetic protein 3 axis in lung adenocarcinoma

Lung adenocarcinoma (LUAD) is a malignant tumor that occurs in the lungs. Numerous reports have substantiated the participation of long non-coding RNAs (lncRNAs) in the tumorigenesis of LUAD. Previously, lncRNA alpha-2-macroglobulin antisense RNA 1 (A2M-AS1) was confirmed to be an important regulator in the biological processes of LUAD and dysregulation of A2M-AS1 was associated with non-small cell lung cancer (NSCLC) progression. However, the precise mechanism of A2M-AS1 in LUAD has not been elucidated. Therefore, our study was designed to investigate the detailed molecular mechanism of A2M-AS1 in LUAD. Herein, the expression of lncRNA A2M-AS1, microRNA (miRNA) miR-587, and bone morphogenetic protein 3 (BMP3) in LUAD cell lines and tissues were detected by real-time quantitative polymerase chain reaction (RT-qPCR) and western blotting. The viability, proliferation, migration and invasion of LUAD cells were tested by cell counting kit-8 (CCK-8), colony formation and Transwell assays. In vivo tumor growth was investigated by xenograft animal experiment. Interactions among A2M-AS1, miR-587 and BMP3 were measured by RNA pulldown and luciferase reporter assays. In this study, A2M-AS1 was downregulated in LUAD tissues and cells and related to poor prognosis in LUAD patients. A2M-AS1 overexpression suppressed LUAD cell proliferation, migration and invasion in vitro and inhibited tumor growth in vivo. Mechanistically, A2M-AS1 directly bound with miR-587 to promote BMP3 expression in LUAD cells. Low expression of BMP3 was found in LUAD tissues and cells and was closely correlated with poor prognosis in LUAD patients. BMP3 deficiency reserved the inhibitory influence of A2M-AS1 overexpression on LUAD cell behaviors. Overall, A2M-AS1 inhibits cell growth and aggressiveness via regulating the miR-587/BMP3 axis in LUAD.

Development and Validation of Lactate Metabolism-Related lncRNA Signature as a Prognostic Model for Lung Adenocarcinoma

BACKGROUND

Lung cancer has been a prominent research focus in recent years due to its role in cancer-related fatalities globally, with lung adenocarcinoma (LUAD) being the most prevalent histological form. Nonetheless, no signature of lactate metabolism-related long non-coding RNAs (LMR-lncRNAs) has been developed for patients with LUAD. Accordingly, we aimed to develop a unique LMR-lncRNA signature to determine the prognosis of patients with LUAD.

METHOD

The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases were utilized to derive the lncRNA expression patterns. Identification of LMR-lncRNAs was accomplished by analyzing the co-expression patterns between lncRNAs and LMR genes. Subsequently, the association between lncRNA levels and survival outcomes was determined to develop an effective signature. In the TCGA cohort, Cox regression was enlisted to build an innovative signature consisting of three LMR-lncRNAs, which was validated in the GEO validation cohort. GSEA and immune infiltration analysis were conducted to investigate the functional annotation of the signature and the function of each type of immune cell.

RESULTS

Fourteen differentially expressed LMR-lncRNAs were strongly correlated with the prognosis of patients with LUAD and collectively formed a new LMR-lncRNA signature. The patients could be categorized into two cohorts based on their LMR-lncRNA signatures: a low-risk and high-risk group. The overall survival of patients with LUAD in the high-risk group was considerably lower than those in the low-risk group. Using Cox regression, this signature was shown to have substantial potential as an independent prognostic factor, which was further confirmed in the GEO cohort. Moreover, the signature could anticipate survival across different groups based on stage, age, and gender, among other variables. This signature also correlated with immune cell infiltration (including B cells, neutrophils, CD4⁺ T cells, CD8⁺ T cells, etc.) as well as the immune checkpoint blockade target CTLA-4.

CONCLUSION

We developed and verified a new LMR-lncRNA signature useful for anticipating the survival of patients with LUAD. This signature could give potentially critical insight for immunotherapy interventions in patients with LUAD.

LRMP Associates With Immune Infiltrates and Acts as a Prognostic Biomarker in Lung Adenocarcinoma

Background: Lymphoid-restricted membrane protein (LRMP) is an endoplasmic reticulum-associated protein that is expressed in a developmentally regulated manner in both B and T cell lineages. However, the role of LRMP in the growth, prognosis and immune infiltration in lung adenocarcinoma (LUAD) remains unclear.

Method: The expression levels of LRMP mRNA in tumor and normal tissues were analyzed using Tumor Immune Estimation Resource 2.0 (TIMER 2.0) and Gene Expression Profiling Interactive Analysis 2 (GEPIA 2). LRMP protein expression was examined using the Human Protein Atlas. In vitro experiments, including qRT-PCR Western blot and immunohistochemistry staining were also performed to investigate LRMP expression. GEPIA2 and Kaplan-Meier plotter databases were used to analyze the clinical prognostic significance of LRMP. To further confirm the underlying function of LRMP, the data were analyzed using gene set enrichment analysis. Moreover, we also constructed plasmids to overexpress LRMP and explored the effect of LRMP in A549 cell line. Additionally, Tumor Immune single-cell Hub was used to investigate the distribution of LRMP in the LUAD immune microenvironment; TIMER and CIBERSORT were used to investigate the relationships among LRMP, LRMP co-expressed genes, and tumor-infiltrating immune cells; Finally, the correlations between LRMP and immune checkpoints were analyzed using TIMER 2.0.

Results: The expression of LRMP was significantly lower in LUAD tissues and cell lines. High LRMP expression is associated with a better prognosis in patients with LUAD. In vitro experimental studies demonstrated that overexpression of LRMP could decrease the proliferation, migration and invasion in A549 cells, and downregulated multiple oncogenic signaling pathways, including p-STAT3, p-PI3K-p-AKT, p-MEK and EMT pathways. GSEA results showed that immuno-related and cell adhesion pathways were enriched in samples with high LRMP expression. LRMP and its co-expressed genes were positively correlated with various tumor-infiltrating immune cells and their markers. Additionally, LRMP positively correlated with immune checkpoints.

Conclusions: Our data suggest that LRMP may act as a tumor suppressor gene and indicates a better prognosis. Moreover, LRMP is associated with immune infiltrates which may be involved in immunotherapy response in LUAD. Further studies are needed to validate these findings.

Aberrant expression of SNHG12 contributes to N, N-dimethylformamide-induced hepatic apoptosis both in short-term and long-term DMF exposure

N, N-Dimethylformamide (DMF) can cause liver damage in occupationally exposed workers, but the molecular mechanism of DMF-induced liver damage has not been fully elucidated. Researches have proved that lncRNA plays a major function in chemical-induced liver toxicity and can be used as a biomarker and therapeutic target for liver injury. In order to verify that lncRNA also participates in DMF-induced liver damage, we treated HL-7702 cells with 75 or 150 mM DMF, and obtained lncRNA expression profiles through high-throughput sequencing. Among the differentially expressed lncRNAs, lncRNA SNHG12 was proved to be significantly downregulated in DMF-treated HL-7702 cells and participate in DMF-mediated apoptosis, even under long-term low-dose DMF exposure (5-10 mM, 8 weeks). In addition, according to bioinformatics analysis, miR-218-5p is expected to be a potential target of SNHG12, which was verified by the dual luciferase reporter assay in HEK293FT cells. MiR-218-5p mimic can induce apoptosis in HL-7702 cells. Among the predicted targets of miR-218-5p, protein kinase C epsilon (PRKCE) was reported to be involved in apoptosis, and was indeed downregulated by miR-218-5p mimic in our study. Further experiments showed that changes of the expression of SNHG12 can affect the expression of PRKCE. In the epidemiological study of occupational population, we also found that SNHG12 was downregulated in the serum exosomes of workers exposed to DMF. These results indicated that SNHG12 can mediate DMF-induced apoptosis of HL-7702 cells through miR-218-5p/PRKCE pathway.