FOXA1-induced LINC01207 facilitates head and neck squamous cell carcinoma via up-regulation of TNRC6B

Proliferative verrucous and homogeneous Leukoplakias exhibit differential methylation patterns

OBJECTIVE

Proliferative verrucous leukoplakia (PVL) is considered a clinically distinct entity from other oral leucoplakias (OLs) due to its clinical presentation and evolution. However, molecular differences between them remain unclear. We aimed to determine whether there are methylation differences between PVL and other forms of OLs.

MATERIALS AND METHODS

Oral biopsies from 12 patients with PVL, eight patients with homogeneous leucoplakia (HL), and 10 healthy individuals were obtained for a genome-wide DNA methylation analysis via the Infinium EPIC Platform.

RESULTS

A total of 1815 differentially methylated CpGs were found between PVL and HL, with a prominent state of hypermethylation in HL patients. CpGs covered 813 genes with distinct roles, including cell adhesion, extracellular matrix organization, and cell and synaptic signaling. 43% of these genes had been previously described in cancer and associated with prognosis. We developed a multinomial logistic regression model able to differentiate HL, PVL, and control samples. The model had a cross-validated estimate of 73% and included differentially methylated cancer-related genes between the pathological conditions and the healthy donors, including ADNP, BRCA2, CDK13, GNB1, NIN, NUMB, PIK3C2B, PTK2, SHISA4, THSD7B, WWP1, and ZNF292. It also included CpGs covering differentially methylated genes in HL (MEN1 and TNRC6B) and PVL (ACOXL, ADH1B, CAMTA1, CBFA2T3, CPXM2, LRFN2, SORCS2, and SPN).

CONCLUSIONS

PVL and HL present differential methylation patterns that could be linked to their differential clinical behavior. Our findings show the potential of methylation markers and suggest novel diagnostic biomarkers.

LncRNA WAC-AS1 promotes osteosarcoma Metastasis and stemness by sponging miR-5047 to upregulate SOX2

Cancer stemness and osteosarcoma (OS) malignant progression are closely associated. However, the molecular mechanisms underlying this association have not been fully demonstrated. Long noncoding RNAs (lncRNAs) are an intriguing class of widely prevalent endogenous RNAs involved in OS progression, the vast majority of which have not been characterized functionally. Here, we identified tumor promoter lncRNA WAC-AS1 to be highly expressed in OS tumors and associated with worse survival. Further analysis revealed that WAC-AS1 increased tumorsphere formation of OS cells and promoted metastasis, as confirmed by cell proliferation, transwell and wound healing assays. MiR-5047 was identified as a downstream target of WAC-AS1. Subsequently, based on bioinformatics analysis, RIP assay and luciferase reporter assay, SOX2 mRNA was verified as a target of miR-5047. WAC-AS1 enhanced OS cell proliferation and stemness via acting as a ceRNA by binding to miR-5047, thereby increasing SOX2 expression. In addition, SOX2 bound to the promoter region of WAC-AS1 and promoted its transcription, thereby forming a positive feedback loop to regulate OS malignancy. Taken together, our findings show WAC-AS1 is a tumor promoter and a key regulator of OS cell stemness and metastasis via a miR-5047/SOX2 axis.

Competing endogenous RNAs in head and neck squamous cell carcinoma: a review

Our understanding of RNA biology has evolved with recent advances in research from it being a non-functional product to molecules of the genome with specific regulatory functions. Competitive endogenous RNA (ceRNA), which has gained prominence over time as an essential part of post-transcriptional regulatory mechanism, is one such example. The ceRNA biology hypothesis states that coding RNA and non-coding RNA co-regulate each other using microRNA (miRNA) response elements. The ceRNA components include long non-coding RNAs, pseudogene and circular RNAs that exert their effect by interacting with miRNA and regulate the expression level of its target genes. Emerging evidence has revealed that the dysregulation of the ceRNA network is attributed to the pathogenesis of various cancers, including the head and neck squamous cell carcinoma (HNSCC). This is the most prevalent cancer developed from the mucosal epithelium in the lip, oral cavity, larynx and pharynx. Although many efforts have been made to comprehend the cause and subsequent treatment of HNSCC, the morbidity and mortality rate remains high. Hence, there is an urgent need to understand the holistic progression of HNSCC, mediated by ceRNA, that can have immense relevance in identifying novel biomarkers with a defined therapeutic intervention. In this review, we have made an effort to highlight the ceRNA biology hypothesis with a focus on its involvement in the progression of HNSCC. For the identification of such ceRNAs, we have additionally highlighted a number of databases and tools.

Integrated analysis of lncRNA-associated ceRNA network in p16-positive and p16-negative head and neck squamous cell carcinoma

Determination of human papillomavirus (HPV) status has become clinically relevant for head and neck squamous cell carcinoma (HNSCC) patients. p16 immunohistochemistry is one of the recommended methods for classifying HPV status. However, long noncoding RNAs (lncRNAs) and related competing endogenous RNA (ceRNA) networks linked to different p16-status HNSCC are still absent. In the present study, The Cancer Genome Atlas database provided RNA profiles as well as clinical information from 26 p16-positive HNSCC samples, 71 p16-negative HNSCC samples, and 44 adjacent normal control samples. Differentially expressed RNAs (DERNAs) between HNSCC samples and normal samples were identified by limma package in R. Functional enrichment analysis of differentially expressed mRNAs was performed using Clusterprofiler package in R. Survival analysis of DERNAs was carried out by survival package in R. The ceRNA network was constructed using GDCRNATools package in R. A total of 102 lncRNAs, 196 microRNAs (miRNAs), and 2282 mRNAs were identified as p16-positive-specific DERNAs. There were 90 lncRNAs, 153 miRNAs, and 2038 mRNAs were identified as p16-negative-specific DERNAs. Functional enrichment analysis revealed that the differentially expressed mRNAs in the p16-positive and the p16-negative group were mainly enriched in the "DNA replication" and "extracellular matrix -receptor interaction" pathway, respectively. Among the top 25 DERNAs, there were 1 key lncRNA, 1 key miRNA, and 1 key messenger RNA in the p16-positive group and 2 key lncRNAs, 1 key miRNA, and 2 key mRNAs in the p16-negative group were significantly related to the overall survival. Then the ceRNA network in the p16-positive and p16-negative group was constructed. There were 5 lncRNAs, 16 miRNAs, and 66 mRNAs included in the p16-positive group ceRNA network and 1 lncRNA, 4 miRNAs, and 28 mRNAs included in the p16-negative group ceRNA network. Among the RNAs in the ceRNA network, 5 mRNAs were significantly related to the overall survival. Taken together, we revealed the differential RNA expression profiling and the differential ceRNA network in the p16-positive and p16-negative group of HNSCC. Our findings provided a novel insight into this HPV-related cancer and potential biomarkers and therapeutic targets for HNSCC based on p16 status.

A bioinformatic approach of targeting SARS-CoV-2 replication by silencing a conserved alternative reserve of the orf8 gene using host miRNAs

The causative agent of the COVID-19 pandemic, the SARS-CoV-2 virus has yielded multiple relevant mutations, many of which have branched into major variants. The Omicron variant has a huge similarity with the original viral strain (first COVID-19 strain from Wuhan). Among different genes, the highly variable orf8 gene is responsible for crucial host interactions and has undergone multiple mutations and indels. The sequence of the orf8 gene of the Omicron variant is, however, identical with the gene sequence of the wild type. orf8 modulates the host immunity making it easier for the virus to conceal itself and remain undetected. Variants seem to be deleting this gene without affecting the viral replication. While analyzing, we came across the conserved orf7a gene in the viral genome which exhibits a partial sequence homology as well as functional similarity with the SARS-CoV-2 orf8. Hence, we have proposed here in our hypothesis that, orf7a might be an alternative reserve of orf8 present in the virus which was compensating for the lost gene. A computational approach was adopted where we screened various miRNAs targeted against the orf8 gene. These miRNAs were then docked onto the orf8 mRNA sequences. The same set of miRNAs was then used to check for their binding affinity with the orf7a reference mRNA. Results showed that miRNAs targeting the orf8 had favorable shape complementarity and successfully docked with the orf7a gene as well. These findings provide a basis for developing new therapeutic approaches where both orf8 and orf7a can be targeted simultaneously.

Hsa_circ_0031608: A Potential Modulator of VSMC Phenotype in the Rupture of Intracranial Aneurysms

Background and Purpose

Phenotypic modulation of vascular smooth muscle cells (VSMCs) plays an important role in the development of intracranial aneurysms (IAs). Growing evidence has demonstrated that circular RNAs (circRNAs) may serve as a potential modulator of VSMC phenotype in various vascular diseases. This study aimed to assess the potential function of circRNAs in the rupture of IAs and VSMC phenotypic modulation.

Methods

Using surgically dissected human ruptured (n = 8) and unruptured (n = 8) IA lesions, differentially expressed circRNAs were screened by transcriptomic sequencing and verified using qRT-PCR. Based on the screened circRNA, we predicted and screened the combined miRNA and downstream mRNAs to construct circRNA-miRNA-mRNA networks. Further in vitro experiments were performed to investigate the relationship between the validated circRNA and the phenotypic switching of VSMCs.

Results

We found 1,373 differentially expressed genes in ruptured versus unruptured aneurysms. The top five dysregulated circRNAs were selected for qRT-PCR validation. We found hsa_circ_0031608 was both highly expressed in ruptured IAs and pro-inflammatory transformation of VSMCs. Then, a regulatory circRNA-miRNA-mRNA with one circRNA node, six miRNA nodes, and 84 mRNA nodes was constructed. GO analysis and KEGG pathway enrichment analysis were performed on mRNAs in the network. Then, a PPI network was built based on these mRNAs and five hub genes were identified (FOXO3, DICER1, CCND2, IGF1R, and TNRC6B) by the cytoHubba plugin in Cytoscape software. In vitro, overexpression of hsa_circ_0031608 influenced the expression of VSMC phenotypic markers validated by qPCR and Western blotting. Furthermore, hsa_circ_0031608 promoted the migration and proliferation capacity of VSMCs.

Conclusion

hsa_circ_0031608 regulated the phenotypic modulation of VSMCs and played an important role in the rupture of IAs. The specific mechanism should be further studied and confirmed.

Eight Differential miRNAs in DN Identified by Microarray Analysis as Novel Biomarkers

BACKGROUND

Diabetic nephropathy (DN) is the common cause of renal diseases such as end-stage renal disease (ESRD) and chronic kidney disease (CKD). Various diagnostic applications and treatment methods are used for clinical but remain some prognosis issues. To avoid morbidity and mortality related to DN, early detection of disease complications as well as targeted therapeutic strategies is essential. Considerable evidence indicates that non-coding RNA plays a vital role in the biological processes of various diseases, used as biomarkers and therapeutic targets. And the most known ncRNAs are the microRNAs (miRNAs), long noncoding RNAs (lncRNAs), and circular RNAs (circRNAs).

MATERIALS AND METHODS

Our study aimed to identify potential prognostic ncRNAs involved in DN by bioinformatics analysis and validated expression levels through quantitative polymerase chain reaction (qPCR) and GEO database. Our research focuses on differential expression miRNAs (DEmiRNAs) in DN and their interactions with critical genes.

RESULTS

We identified 8 up-regulated DEmiRNAs, including miR-103a-2-5p, miR-297, miR-548x-3p, miR-604, miR-644a, miR-1256, miR-3911 and miR-5047 finally. We further validated these miRNAs in a murine model.

CONCLUSION

Identifying these up-regulated genes and elucidating these miRNAs regulatory network will contribute to a better understanding of the molecular mechanism of DN and how they can be used as new biomarkers and potential therapeutic targets for DN.

Long Noncoding Competing Endogenous RNA Networks in Pancreatic Cancer

Pancreatic cancer (PC) is a highly malignant disease characterized by insidious onset, rapid progress, and poor therapeutic effects. The molecular mechanisms associated with PC initiation and progression are largely insufficient, hampering the exploitation of novel diagnostic biomarkers and development of efficient therapeutic strategies. Emerging evidence recently reveals that noncoding RNAs (ncRNAs), including long ncRNAs (lncRNAs) and microRNAs (miRNAs), extensively participate in PC pathogenesis. Specifically, lncRNAs can function as competing endogenous RNAs (ceRNAs), competitively sequestering miRNAs, therefore modulating the expression levels of their downstream target genes. Such complex lncRNA/miRNA/mRNA networks, namely, ceRNA networks, play crucial roles in the biological processes of PC by regulating cell growth and survival, epithelial-mesenchymal transition and metastasis, cancer stem cell maintenance, metabolism, autophagy, chemoresistance, and angiogenesis. In this review, the emerging knowledge on the lncRNA-associated ceRNA networks involved in PC initiation and progression will be summarized, and the potentials of the competitive crosstalk as diagnostic, prognostic, and therapeutic targets will be comprehensively discussed.

Long non-coding RNA LINC01207 promotes cell proliferation and migration but suppresses apoptosis and autophagy in oral squamous cell carcinoma by the microRNA-1301-3p/lactate dehydrogenase isoform A axis

Long noncoding RNAs (lncRNAs) have been reported to participate in the progression of various cancers, including oral squamous cell carcinoma (OSCC). This study aims to find out whether lncRNA LINC01207 regulates the progression of OSCC. Reverse transcription quantitative polymerase chain reaction (RT-qPCR) was conducted to evaluate gene expression in OSCC cells and tissues. Cell viability, proliferation, migration, apoptosis, and autophagy were detected using Cell Counting Kit-8 (CCK-8), colony formation, Transwell assays, flow cytometry, and western blot analysis. Luciferase reporter and RNA immunoprecipitation (RIP) assays were conducted to assess the interactions among genes. We found that LINC01207 was overexpressed in OSCC cells and tissues. LINC01207 silencing inhibited OSCC cell proliferation and migration but promoted apoptosis and autophagy, and LINC01207 overexpression had an opposite result. LINC01207 interacted with microRNA-1301-3p (miR-1301-3p) while lactate dehydrogenase isoform A (LHDA) was targeted by miR1301-3p. Effects caused by LINC01207 downregulation on OSCC cells were reversed by overexpression of LDHA. Overall, LINC01207 promotes OSCC progression via the miR-1301-3p/LDHA axis

LINC01207 is up-regulated in gastric cancer tissues and promotes disease progression by regulating miR-671-5p/DDX5 axis

LINC01207 is involved in the progression of some cancers. This study was designed to delve into the biological function and mechanism of LINC01207 in gastric cancer. qPCR was adopted to examine the expression levels of LINC01207, miR-671-5p, DDX5 mRNA in gastric cancer tissues and cells. After LINC01207 was overexpressed or depleted, MTT and BrdU assays were conducted to detect cell proliferation. Transwell assay was employed to detect cell migration and invasion. Western blot was used to detect the expression of DDX5 protein in cells. Bioinformatics analysis, luciferase reporter assay and RNA pull-down assay were performed to predict and validate the binding site between miR-671-5p and LINC01207 or DDX5. LINC01207, DDX5 mRNA were up-regulated in gastric cancer, while miR-671-5p was down-regulated; high expression of LINC01207 and transfection of miR-671-5p inhibitors facilitated the proliferation of gastric cancer cells; however, knocking down LINC01207 and the overexpression of miR-671-5p mimics had opposite biological effects. LINC01207 and miR-671-5p were interacted and miR-671-5p was negatively regulated by LINC01207. MiR-671-5p could reverse the function of LINC01207. DDX5 was a downstream target of miR-671-5p and was positively modulated by LINC01207. LINC01207 promotes the proliferation and metastasis of gastric cancer cells by regulating miR-671-5p/DDX5 axis.

CTCF-induced upregulation of LINC01207 promotes gastric cancer progression via miR-1301-3p/PODXL axis

BACKGROUND

Long non coding RNAs (lncRNAs) have been validated to be involved in the complicated biological processes during tumor progression. LINC01207 has been identified as an oncogene in several cancer types. However, the function of LINC01207 and its underlying molecular mechanism in gastric cancer (GC) are poorly understood.

METHODS

The expression level of LINC01207, miR-1301-3p and PODXL mRNA was detected in GC tissues and cells by RT-qPCR. The level of PODXL protein was examined by western blot. Colony formation assay, EdU assay, TUNEL assay, caspase-3 activity test and transwell assays were carried out to analyze the effect of LINC01207 on GC cell proliferation, apoptosis, migration and invasion. The interaction between RNAs was confirmed by luciferase reporter assay, RNA pull-down assay and RIP assay.

RESULTS

LINC01207 was expressed at high level in GC tissues and cells. Silencing of LINC01207 impaired GC cell proliferation, migration and invasion but promoted cell apoptosis. Mechanistically, LINC01207 acted as a ceRNA by sponging miR-1301-3p to upregulate PODXL. Besides, miR-1301-3p silencing or PODXL overexpression could abolish the inhibitory effect of LINC01207 knockdown on GC cell growth and migration. CCCTC-binding factor (CTCF) could transcriptionally activate LINC01207 in GC cells.

CONCLUSIONS

CTCF-induced activation of LINC01207 contributes to GC progression through regulating miR-1301-3p/PODXL axis.

Identification of MYLIP gene and miRNA-802 involved in the growth and metastasis of cervical cancer cells

BACKGROUND

The dysregulation of microRNA-802 (miR-802) has crucial roles in cancer progression. Nevertheless, the bio-function of miR-802 in cervical cancer remains unclear.

OBJECTIVE

Hence, we illuminated the potential roles of miR-802 in cervical cancer cell growth, migration, and invasion.

METHODS

The levels of miR-802 and myosin regulatory light chain interacting protein (MYLIP) were measured using qRT-PCR assay. The potential effects of miRNA-802 on cervical cancer cell proliferation and metastatic phenotypes were determined using CCK-8, colony formation, wound healing and Transwell invasion assays. MYLIP was validated as a downstream target gene of miRNA-802 using bioinformatics analysis tool and luciferase report gene assay. The impact of miR-802 on the growth of cervical cancer cell in vivo was analyzed using xenograft model. The expression of MYLIP was measured by western blotting and immunohistochemistry (IHC).

RESULTS

MiRNA-802 was distinctly down-regulated in cervical cancer cells as well as clinical cervical cancer samples. Upregulation of miRNA-802 significantly inhibited the growth and aggressiveness of cervical cancer cell. Additional, MYLIP was a functional target of miR-802. MYLIP was ovrerexpressed in cervical cancer and MYLIP level was negatively associated with the level of miR-802. Overexpression of MYLIP eliminated the inhibitory effects of miR-802 on growth and metastatic-related traits of cervical cancer cell. In vivo, miR-802 also markedly reduced the tumor growth of cervical cancer cell and decreased the expression of MYLIP.

CONCLUSIONS

MiR-802 inhibits the growth and metastatic-related phenotypes of cervical cancer cell through targeting MYLIP.