Pseudogene-expressed RNAs: a new frontier in cancers

Current status and future perspectives in dysregulated miR-492

MicroRNAs (miRNAs) are a class of single-stranded small non-coding RNAs with a length of 21-23 nucleotides. One such miRNA, miR-492, is located in the KRT19 pseudogene 2 (KRT19P2) of chromosome 12q22 and can also be generated from the processing of the KRT19 transcript at chromosome 17q21. Aberrant expression of miR-492 has been observed in cancers of various physiological systems. miR-492 has been shown to target at least 11 protein-coding genes, which are involved in the regulation of cellular behaviors such as growth, cell cycle, proliferation, epithelial- mesenchymal transition (EMT), invasion and migration. The expression of miR-492 can be regulated by both endogenous and exogenous factors. Furthermore, miR-492 is involved in the regulation of several signaling pathways including the PI3K/AKT signaling pathway, WNT/β-catenin signaling pathway, and MAPK signaling pathway. High expression of miR-492 has been closely associated with shorter overall survival in patients with gastric cancer, ovarian cancer, oropharyngeal carcinoma, colorectal cancer, and hepatocellular carcinoma. This study systematically summarizes the related research findings on miR-492, providing potential insights for future investigations.

The pseudogene DUXAP10 contributes to gefitinib resistance in NSCLC by repressing OAS2 expression

Gefitinib, an epidermal growth factor receptor-tyrosine kinase inhibitor (EGFR-TKI),is the currently recommended first-line therapy for advanced EGFR-mutant lung cancer, and understanding the mechanism of resistance is the key to formulating therapeutic strategies for EGFR-TKIs. In this study, we evaluate the expression patterns and potential biological functions of the pseudogene DUXAP10 in gefitinib resistance. We find that pseudogene DUXAP10 expression is significantly upregulated in NSCLC gefitinib-resistant cells and tissues. Gain and loss of function assays reveal that knockdown of DUXAP10 by siRNA reverses gefitinib resistance both in vitro and in vivo. Furthermore, DUXAP10 interacts with the histone methyltransferase enhancer of zeste homolog 2 (EZH2) to repress the expression of 2',5'-oligoadenylate synthetase (OAS2). Overall, our study highlights the pivotal role of DUXAP10 in gefitinib resistance, and the DUXAP10/EZH2/OAS2 axis might be a promising therapeutic target to overcome acquired gefitinib resistance in NSCLC.

Scavenger: A pipeline for recovery of unaligned reads utilising similarity with aligned reads

Read alignment is an important step in RNA-seq analysis as the result of alignment forms the basis for downstream analyses. However, recent studies have shown that published alignment tools have variable mapping sensitivity and do not necessarily align all the reads which should have been aligned, a problem we termed as the false-negative non-alignment problem. Here we present Scavenger, a python-based bioinformatics pipeline for recovering unaligned reads using a novel mechanism in which a putative alignment location is discovered based on sequence similarity between aligned and unaligned reads. We showed that Scavenger could recover unaligned reads in a range of simulated and real RNA-seq datasets, including single-cell RNA-seq data. We found that recovered reads tend to contain more genetic variants with respect to the reference genome compared to previously aligned reads, indicating that divergence between personal and reference genomes plays a role in the false-negative non-alignment problem. Even when the number of recovered reads is relatively small compared to the total number of reads, the addition of these recovered reads can impact downstream analyses, especially in terms of estimating the expression and differential expression of lowly expressed genes, such as pseudogenes.

A Penalization Method for Estimating Heterogeneous Covariate Effects in Cancer Genomic Data

In high-throughput profiling studies, extensive efforts have been devoted to searching for the biomarkers associated with the development and progression of complex diseases. The heterogeneity of covariate effects associated with the outcomes across subjects has been noted in the literature. In this paper, we consider a scenario where the effects of covariates change smoothly across subjects, which are ordered by a known auxiliary variable. To this end, we develop a penalization-based approach, which applies a penalization technique to simultaneously select important covariates and estimate their unique effects on the outcome variables of each subject. We demonstrate that, under the appropriate conditions, our method shows selection and estimation consistency. Additional simulations demonstrate its superiority compared to several competing methods. Furthermore, applying the proposed approach to two The Cancer Genome Atlas datasets leads to better prediction performance and higher selection stability.

Predicting Pseudogene-miRNA Associations Based on Feature Fusion and Graph Auto-Encoder

Pseudogenes were originally regarded as non-functional components scattered in the genome during evolution. Recent studies have shown that pseudogenes can be transcribed into long non-coding RNA and play a key role at multiple functional levels in different physiological and pathological processes. microRNAs (miRNAs) are a type of non-coding RNA, which plays important regulatory roles in cells. Numerous studies have shown that pseudogenes and miRNAs have interactions and form a ceRNA network with mRNA to regulate biological processes and involve diseases. Exploring the associations of pseudogenes and miRNAs will facilitate the clinical diagnosis of some diseases. Here, we propose a prediction model PMGAE (Pseudogene–MiRNA association prediction based on the Graph Auto-Encoder), which incorporates feature fusion, graph auto-encoder (GAE), and eXtreme Gradient Boosting (XGBoost). First, we calculated three types of similarities including Jaccard similarity, cosine similarity, and Pearson similarity between nodes based on the biological characteristics of pseudogenes and miRNAs. Subsequently, we fused the above similarities to construct a similarity profile as the initial representation features for nodes. Then, we aggregated the similarity profiles and associations of nodes to obtain the low-dimensional representation vector of nodes through a GAE. In the last step, we fed these representation vectors into an XGBoost classifier to predict new pseudogene–miRNA associations (PMAs). The results of five-fold cross validation show that PMGAE achieves a mean AUC of 0.8634 and mean AUPR of 0.8966. Case studies further substantiated the reliability of PMGAE for mining PMAs and the study of endogenous RNA networks in relation to diseases.

m6A Modification-Mediated DUXAP8 Regulation of Malignant Phenotype and Chemotherapy Resistance of Hepatocellular Carcinoma Through miR-584-5p/MAPK1/ERK Pathway Axis

Hepatocellular carcinoma (HCC) has a poor prognosis due to its high malignancy, rapid disease progression, and the presence of chemotherapy resistance. Long-stranded non-coding RNAs (lncRNAs) affect many malignant tumors, including HCC. However, their mechanism of action in HCC remains unclear. This study aimed to clarify the role of DUXAP8 in regulating the malignant phenotype and chemotherapy resistance in HCC. Using an in vivo xenograft tumor model, the regulatory functions and mechanisms of lncRNA DUXAP8 in the progression and response of HCC to chemotherapy were explored. It was found that DUXAP8 was significantly upregulated in a patient-derived xenograft tumor model based on sorafenib treatment, which is usually associated with a relatively poor prognosis in patients. In HCC, DUXAP8 maintained its upregulation in the expression by increasing the stability of m6A methylation-mediated RNA. DUXAP8 levels were positively correlated with the proliferation, migration, invasion, and chemotherapy resistance of HCC in vivo and in vitro. In the mechanistic study, it was found that DUXAP8 competitively binds to miR-584-5p through a competing endogenous RNA (ceRNA) mechanism, thus acting as a molecular sponge for miR-584-5p to regulate MAPK1 expression, which in turn activates the MAPK/ERK pathway. These findings can provide ideas for finding new prognostic indicators and therapeutic targets for patients with HCC.

Competing Endogenous RNA Networks in Glioma

Gliomas are the most common and malignant primary brain tumors. Various hallmarks of glioma, including sustained proliferation, migration, invasion, heterogeneity, radio- and chemo-resistance, contribute to the dismal prognosis of patients with high-grade glioma. Dysregulation of cancer driver genes is a leading cause for these glioma hallmarks. In recent years, a new mechanism of post-transcriptional gene regulation was proposed, i.e., "competing endogenous RNA (ceRNA)." Long non-coding RNAs, circular RNAs, and transcribed pseudogenes act as ceRNAs to regulate the expression of related genes by sponging the shared microRNAs. Moreover, coding RNA can also exert a regulatory role, independent of its protein coding function, through the ceRNA mechanism. In the latest glioma research, various studies have reported that dysregulation of certain ceRNA regulatory networks (ceRNETs) accounts for the abnormal expression of cancer driver genes and the establishment of glioma hallmarks. These achievements open up new avenues to better understand the hidden aspects of gliomas and provide new biomarkers and potential efficient targets for glioma treatment. In this review, we summarize the existing knowledge about the concept and logic of ceRNET and highlight the emerging roles of some recently found ceRNETs in glioma progression.

Long non-coding RNA DUXAP8 promotes the cell proliferation, migration, and invasion of papillary thyroid carcinoma via miR-223-3p mediated regulation of CXCR4

Papillary thyroid carcinoma (PTC) is a differentiated type of thyroid malignancy with a high incidence. Long non-coding RNA (lncRNA) DUXAP8 has been reported to participate in the proliferation, migration, and invasion of several cancer types. However, its association with PTC has not yet been reported. The current study aimed to investigate the role of DUXAP8 in PTC and revealed the underlying mechanisms. The expression of DUXAP8 was knocked down in two PTC cell lines and the effects of DUXAP8 on the PTC biological behavior were examined by cell counting kit-8 (CCK-8), wound healing, and transwell invasion assays. Luciferase reporter assay was used to detect the binding activity between miR-223-3p and DUXAP8. We found that knockdown of DUXAP8 inhibited the proliferation, migration, and invasion of PTC cells. DUXAP8 could sponge miR-223-3p through the specific binding site. CXCR4 was a target of miR-223-3p. The malignant phenotypes of the PTC cells were suppressed by the over-expression of miR-223-3p. Moreover, miR-223-3p inhibition or CXCR4 over-expression partly restored the proliferation, migration, and invasion activities of DUXAP8-downregulated PTC cells. The results evidenced that DUXAP8 acted as an oncogene in PTC, these effects seemed to partly dependent on the miR-223-3p/CXCR4 axis.

The Pseudogene DUXAP8 Promotes Colorectal Cancer Cell Proliferation, Invasion, and Migration by Inducing Epithelial-Mesenchymal Transition Through Interacting with EZH2 and H3K27me3

Background

Colorectal cancer (CRC) is the third leading cause of cancer death worldwide. The long noncoding RNA (lncRNA) DUXAP8 has been reported to play an important role in CRC. This study investigated the mechanism by which this lncRNA regulates CRC progression.

Methods

The levels of lncRNA DUXAP8 in CRC tissues and cell lines were detected by qRT-PCR. We then knocked down or forced overexpression of DUXAP8, and the resultant effect on cell proliferation was determined by the Edu assay and a cell cycle analysis, and the effect on cell apoptosis was determined by flow cytometry. The cell invasion/migration ability and the epithelial-to-mesenchymal transition (EMT) markers were determined by Transwell/wound healing assays and Western blotting. CHIP and RNA pull-down assays were performed to determine the binding of Zeste gene enhancer 2 (EZH2) and trimethylated histone H3 to Lys27 (H3K27me3) in the E-cadherin promoter regions, or to DUXAP8.

Results

The levels of lncRNA DUXAP8 were significantly increased in CRC tissues and CRC cell lines. Knockdown of lncRNA DUXAP8 inhibited cell proliferation and the EMT process, and increased cell apoptosis, and overexpression of lncRNA DUXAP8 had an opposite effect. Both ChIP and RNA pull-down assays showed that the E-cadherin promoter region was bound by H3K27me3 and EZH2, which restrained E-cadherin expression. However, that binding was suppressed and E-cadherin expression was markedly induced by lncRNA DUXAP8 knockdown. Furthermore, lncRNA DUXAP8 could interact with EZH2 and H3K27me3.

Conclusion

Our data indicated that lncRNA DUXAP8 could induce the progression of CRC by negatively regulating E-cadherin via interaction with EZH2 and H3K27me3. These findings suggest lncRNA DUXAP8 as target for treating CRC.

Interferon-Mediated Long Non-Coding RNA Response in Macrophages in the Context of HIV

Interferons play a critical role in the innate immune response against a variety of pathogens, such as HIV-1. Recent studies have shown that long non-coding genes are part of a reciprocal feedforward/feedback relationship with interferon expression. They presumably contribute to the cell type specificity of the interferon response, such as the phenotypic and functional transition of macrophages throughout the immune response. However, no comprehensive understanding exists today about the IFN–lncRNA interplay in macrophages, also a sanctuary for latent HIV-1. Therefore, we completed a poly-A+ RNAseq analysis on monocyte-derived macrophages (MDMs) treated with members of all three types of IFNs (IFN-α, IFN-ε, IFN-γ or IFN-λ) and on macrophages infected with HIV-1, revealing an extensive non-coding IFN and/or HIV-1 response. Moreover, co-expression correlation with mRNAs was used to identify important (long) non-coding hub genes within IFN- or HIV-1-associated gene clusters. This study identified and prioritized IFN related hub lncRNAs for further functional validation.

Identification of male infertility-related long non-coding RNAs and their functions based on a competing endogenous RNA network

OBJECTIVE

To identify male infertility-related long non-coding (lnc)RNAs and an lncRNA-related competing endogenous (ce)RNA network.

METHODS

Expression data including 13 normospermic and eight teratozoospermic samples from postmortem donors were downloaded from the GEO database (GSE6872). The limma R package was used to discriminate dysregulated lncRNA and micro (m)RNA profiles. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses of differentially expressed (DE) mRNAs were performed using the clusterProfiler R package. The ceRNA network of dysregulated genes was visualized by Cytoscape.

RESULTS

A total of 101 DE lncRNAs and 1722 mRNAs were identified as male infertility-specific RNAs with thresholds of |log₂FoldChange| >2.0 and adjusted P-value <0.05. GO and KEGG pathways were analyzed for DE mRNAs. Gene set enrichment analysis revealed that DE genes were enriched in embryonic skeletal system development and cytokine-cytokine receptor interactions. A ceRNA network was constructed with 26 key lncRNAs, 33 microRNAs, and 133 mRNAs. DE lncRNAs in male sterility were mainly associated with transferring phosphorus-containing groups and complexes of histone methyltransferases, methyltransferases, PcG proteins, and serine/threonine protein kinases.

CONCLUSION

This provides a novel perspective to study lncRNA-related ceRNA networks in male infertility and assist in identifying new potential biomarkers for diagnostic purposes.

Functional characterization of long noncoding RNAs

PURPOSE OF REVIEW

Mounting evidence suggests that long noncoding RNAs (lncRNAs) are essential regulators of gene expression. Although few lncRNAs have been the subject of detailed molecular and functional characterization, it is believed that lncRNAs play an important role in tissue homeostasis and development. In fact, gene expression profiling studies reveal lncRNAs are developmentally regulated in a tissue-type and cell-type specific manner. Such findings have brought significant attention to their potential contribution to disease cause. The current review summarizes recent studies of lncRNAs in the heart.

RECENT FINDINGS

lncRNA discovery has largely been driven by the implementation of next generation sequencing technologies. To date, such technologies have contributed to the identification of tens of thousands of distinct lncRNAs in humans -- accounting for a large majority of all RNA sequences transcribed across the human genome. Although the functions of these lncRNAs remain largely unknown, gain-of-function and loss-of-function studies (in vivo and in vitro) have uncovered a number of mechanisms by which lncRNAs regulate gene expression and protein function. Such mechanisms have been stratified according to three major functional categories: RNA sponges (RNA-mediated sequestration of free miRNAs; e.g. H19, MEG3, and MALAT1); transcription-modulating lncRNAs (RNA influences regulatory factor recruitment by binding to histone modifiers or transcription factors; e.g. CAIF, MANTIS, and NEAT1); and translation-modulating lncRNAs (RNA modifies protein function via directly interacting with a protein itself or binding partners; e.g. Airn, CCRR, and ZFAS1).

SUMMARY

Recent studies strongly suggest that lncRNAs function via binding to macromolecules (e.g. genomic DNA, miRNAs, or proteins). Thus, lncRNAs constitute an additional mode by which cells regulate gene expression.

The LGMN pseudogene promotes tumor progression by acting as a miR-495-3p sponge in glioblastoma

Pseudogenes, which are long noncoding RNAs that originate from protein-coding genes, have been suggested to play important roles in disease. Although studies have revealed high expression of legumain (LGMN) in many types of tumors, the regulation of LGMN remains largely unknown. Here, we found that a novel LGMN pseudogene (LGMNP1) was upregulated in glioblastoma (GBM) tissues and high LGMNP1 expression in GBM cells enhanced proliferation and invasion. Biochemical analysis showed that cytoplasmic LGMNP1 functionally targeted miR-495-3p in a manner involving an RNA-induced silencing complex. Dual-luciferase reporter assays demonstrated that LGMN was a target of miR-495-3p, and LGMN was upregulated and positively correlated with LGMNP1 in GBM. Moreover, miR-495-3p was downregulated and negatively correlated with LGMNP1 in GBM tissues. Notably, the tumor-promoting effects of LGMNP1 upregulation could be alleviated by miR-495-3p mimics. Furthermore, GBM cells overexpressing LGMNP1 exhibited more aggressive tumor progression and elevated LGMN expression in vivo. Thus, our data illustrate that LGMNP1 exerts its oncogenic activity, at least in part, as a competitive endogenous RNA (ceRNA) that elevates LGMN expression by sponging miR-495-3p. CeRNA-mediated miRNA sequestration might be a novel therapeutic strategy in GBM.

Long noncoding RNA FOXD3-AS1 promotes colon adenocarcinoma progression and functions as a competing endogenous RNA to regulate SIRT1 by sponging miR-135a-5p

More and more documents have proved that the abnormal expression of long noncoding RNAs (lncRNAs) are correlated with the initiation and progression of colorectal cancer (CRC). lncRNA FOXD3-AS1 has been reported in glioma for its oncogenic property. According to the survival analysis of The Cancer Genome Atlas database, FOXD3-AS1 upregulation implied lower survival rate of patients with CRC. Quantitative real-time polymerase chain reaction showed the overexpression of FOXD3-AS1 in both CRC tissues and cells. The Kaplan-Meier method demonstrated the prognostic value of FOXD3-AS1 for patients with CRC. To explore the effect of FOXD3-AS1 on CRC progression, loss-of-function experiments were carried out, whose results indicated that knockdown of FOXD3-AS1 suppressed cell proliferation, migration, and invasion, inhibited cell cycle and promoted cell apoptosis in vitro. In vivo experiments affirmed that FOXD3-AS1 affected tumor growth. FOXD3-AS1 expression was enriched in the cytoplasm of CRC cells. Mechanism experiments revealed that FOXD3-AS1 served as a competing endogenous RNA to upregulate SIRT1 by sponging miR-135a-5p. In addition, SIRT1 silencing also restrained cell proliferation and motility. Rescue assays revealed the biological function of FOXD3-AS1/miR-135a-5p/SIRT1 axis in CRC progression. In conclusion, FOXD3-AS1 promotes CRC progression by regulating miR-135a-5p/SIRT1 axis, shedding lights on the way to CRC treatments.

Identification of Long Non-Coding RNA and Circular RNA Expression Profiles in Atrial Fibrillation

BACKGROUND

Long non-coding RNA (lncRNA) and circular RNA (circRNA) have both been found to play important roles in cardiovascular diseases, including myocardial infarction, heart failure, and atherosclerosis. However, the role of lncRNA and circRNA in atrial fibrillation (AF) has rarely been investigated. This study aimed to identify lncRNA and circRNA expression profiles in AF patients.

METHODS

Atrial tissues from seven patients with AF and seven matched controls were collected. The lncRNA and circRNA expression profiles of atrial tissues were identified using Hiseq/Proton RNA sequencing. Validation was performed by reverse transcription quantitative real-time PCR (qRT-PCR) on 35 pairs of AF patients and controls. Gene Ontology (GO) categories and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed. A competing endogenous RNA (ceRNA) network was constructed.

RESULTS

A total of 557 lncRNAs and 280 circRNAs were significantly differentially expressed with fold change >1.5 (p<0.05). An lncRNA Voltage Dependent Anion Channel 2 Pseudogene 2 (VDAC2P2) and two circRNAs chr13_41887361_41865736_-21625 and chr13_100368574_100301460_-67114 were validated, using qRT-PCR, to have significantly different expression levels. GO and KEGG pathway analysis showed that some pathways such as ribosome and chromatin modification, Rap1 signalling and cardiac muscle contraction were involved in the pathogenesis of AF. Competing endogenous RNAs were predicted based on constructional network analysis. The LncRNA-miRNA-mRNA and circRNA-miRNA-mRNA networks were constructed by co-expressing lncRNA/circRNA and mRNAs, which were competitively combined with miRNAs.

CONCLUSION

This study characterised lncRNA and circRNA expression and their interaction with mRNA and miRNA in AF.

Classification of early and late stage liver hepatocellular carcinoma patients from their genomics and epigenomics profiles

Background

Liver Hepatocellular Carcinoma (LIHC) is one of the major cancers worldwide, responsible for millions of premature deaths every year. Prediction of clinical staging is vital to implement optimal therapeutic strategy and prognostic prediction in cancer patients. However, to date, no method has been developed for predicting the stage of LIHC from the genomic profile of samples.

Methods

The Cancer Genome Atlas (TCGA) dataset of 173 early stage (stage-I), 177 late stage (stage-II, Stage-III and stage-IV) and 50 adjacent normal tissue samples for 60,483 RNA transcripts and 485,577 methylation CpG sites, was extensively analyzed to identify the key transcriptomic expression and methylation-based features using different feature selection techniques. Further, different classification models were developed based on selected key features to categorize different classes of samples implementing different machine learning algorithms.

Results

In the current study, in silico models have been developed for classifying LIHC patients in the early vs. late stage and cancerous vs. normal samples using RNA expression and DNA methylation data. TCGA datasets were extensively analyzed to identify differentially expressed RNA transcripts and methylated CpG sites that can discriminate early vs. late stages and cancer vs. normal samples of LIHC with high precision. Naive Bayes model developed using 51 features that combine 21 CpG methylation sites and 30 RNA transcripts achieved maximum MCC (Matthew’s correlation coefficient) 0.58 with an accuracy of 78.87% on the validation dataset in discrimination of early and late stage. Additionally, the prediction models developed based on 5 RNA transcripts and 5 CpG sites classify LIHC and normal samples with an accuracy of 96–98% and AUC (Area Under the Receiver Operating Characteristic curve) 0.99. Besides, multiclass models also developed for classifying samples in the normal, early and late stage of cancer and achieved an accuracy of 76.54% and AUC of 0.86.

Conclusion

Our study reveals stage prediction of LIHC samples with high accuracy based on the genomics and epigenomics profiling is a challenging task in comparison to the classification of cancerous and normal samples. Comprehensive analysis, differentially expressed RNA transcripts, methylated CpG sites in LIHC samples and prediction models are available from CancerLSP (http://webs.iiitd.edu.in/raghava/cancerlsp/).

Low Expression of Pseudogene POU5F1B Affects Diagnosis and Prognosis in Acute Myeloid Leukemia (AML)

Background

The transcription factor Oct-4 is necessary for maintaining pluripotency and self-renewal of embryonic stem cells, and POU5F1B is a processed pseudogene of Oct-4 with coding capacity. The purpose of this study is to evaluate the expression and clinical implication of POU5F1B in AML.

Material/Methods

The expression of the POU5F1B transcript was evaluated in 175 newly diagnosed AML patients and 39 healthy controls by use of real-time quantitative PCR (RQ-PCR).

Results

POU5F1B was underexpressed in AML compared with controls (P<0.001). The receiver operating characteristic (ROC) curve revealed that the POU5F1B transcript level was able to differentiate AML patients from healthy individuals (AUC=0.682). In non-APL AML patients, the POU5F1B^low group had significantly higher WBC than the POU5F1B^high group (20.2×10⁹vs. 4.6×10⁹ L⁻¹, P=0.021). Among whole-cohort AML, non-APL AML, and intermediate-risk AML, POU5F1B^high patients had obviously higher complete remission (CR) rates than POU5F1B^low patients (P=0.012, P=0.012 and P=0.027). In addition, Kaplan-Meier analysis demonstrated better overall survival (OS, P=0.019, P=0.007 and P=0.046, respectively) in POU5F1B^high patients compared with POU5F1B^low patients. Furthermore, in multivariate survival analysis, POU5F1B was independently associated with OS in non-APL AML patients and intermediate-risk AML as a favorable prognostic factor.

Conclusions

POU5F1B was frequently underexpressed in AML, and might contribute to the diagnosis and prognosis of AML.

Identification of a Five-Pseudogene Signature for Predicting Survival and Its ceRNA Network in Glioma

Background: Glioma is the most common primary brain tumor with a dismal prognosis. It is urgent to develop novel molecular biomarkers and conform to individualized schemes. Methods: Differentially expressed pseudogenes between low grade glioma (LGG) and glioblastoma multiforme (GBM) were identified in the training cohort. Least absolute shrinkage and selection operator (LASSO) regression and multivariate Cox proportional hazards regression analyses were used to select pseudogenes associated with prognosis of glioma. A risk signature was constructed based on the selected pseudogenes for predicting the survival of glioma patients. A pseudogene-miRNA-mRNA regulatory network was established and visualized using Cytoscape 3.5.1. Gene Oncology (GO) and signaling pathway analyses were performed on the targeted genes to investigate functional roles of the risk signature. Results: Five pseudogenes (ANXA2P2, EEF1A1P9, FER1L4, HILS1, and RAET1K) correlating with glioma survival were selected and used to establish a risk signature. Time-dependent receiver operating characteristic (ROC) curves revealed that the risk signature could accurately predict the 1, 3, and 5-year survival of glioma patients. GO and signaling pathway analyses showed that the risk signature was involved in regulation of proliferation, migration, angiogenesis, and apoptosis in glioma. Conclusions: In this study, a risk signature with five pseudogenes was constructed and shown to accurately predict 1-, 3-, and 5-year survival for glioma patient. The risk signature may serve as a potential target against glioma.

Pseudogenes as Functionally Significant Elements of the Genome

Pseudogene is a gene copy that has lost its original function. For a long time, pseudogenes have been considered as "junk DNA" that inevitably arises as a result of ongoing evolutionary process. However, experimental data obtained during recent years indicate this understanding of the nature of pseudogenes is not entirely correct, and many pseudogenes perform important genetic functions. In the review, we have addressed classification of pseudogenes, methods of their detection in the genome, and the problem of their evolutionary conservatism and prevalence among species belonging to different taxonomic groups in the light of modern data. The mechanisms of gene expression regulation by pseudogenes and the role of pseudogenes in pathogenesis of various human diseases are discussed.

DUXAP8, a pseudogene derived lncRNA, promotes growth of pancreatic carcinoma cells by epigenetically silencing CDKN1A and KLF2

BACKGROUND

Recent studies highlight pseudogene derived long non-coding RNAs (lncRNAs) as key regulators of cancer biology. However, few of them have been well characterized in pancreatic cancer. Here, we aimed to identify the association between pseudogene derived lncRNA DUXAP8 and growth of pancreatic cancer cells.

METHODS

We screened for pseudogene derived lncRNAs associated with human pancreatic cancer by comparative analysis of three independent datasets from GEO. Quantitative real-time reverse transcription polymerase chain reaction was used to assess the relative expression of DUXAP8 in pancreatic cancer tissues and cells. Loss-of-function approaches were used to investigate the potential functional roles of DUXAP8 in pancreatic cancer cell proliferation and apoptosis in vitro and in vivo. RNA immunoprecipitation, chromosome immunoprecipitation assay and rescue experiments were performed to analyze the association of DUXAP8 with target proteins and genes in pancreatic cancer cells.

RESULTS

Pancreatic cancer tissues had significantly higher DUXAP8 levels than paired adjacent normal tissues. High DUXAP8 expression was associated with a larger tumor size, advanced pathological stage and shorter overall survival of pancreatic cancer patients. Moreover, silencing DUXAP8 expression by siRNA or shRNA inhibited pancreatic cancer cell proliferation and promoted apoptosis in vitro and in vivo. Mechanistic analyses indicated that DUXAP8 regulates PC cell proliferation partly through downregulation of tumor suppressor CDKN1A and KLF2 expression.

CONCLUSION

Our results suggest that tumor expression of pseudogene derived lncRNA DUXAP8 plays an important role in pancreatic cancer progression. DUXAP8 may serve as a candidate biomarker and represent a novel therapeutic target of pancreatic cancer.

TUSC2P suppresses the tumor function of esophageal squamous cell carcinoma by regulating TUSC2 expression and correlates with disease prognosis

BACKGROUND

Pseudogenes are RNA transcripts with high homology with its parent protein-coding genes. Although pseudogenes lost the ability to produce protein, it still exert import biological function, and play important role in the pathogenesis of a wide varity of tumors; However, the role of pseudogenes in esophageal squamous cell carcinoma (ESCC) is poorly understood.

METHODS

TUSC2P function in ESCC were explored using both in vitro and in vivo experiments cell proliferation, invasion and apoptosis assay was performed to evaluated the effect of TUSC2P on the tumor biology of ESCC. Expression of relative genes was assessed by quantitative real-time PCR (qRT-PCR) and western blotting in EC109 and TE-1 cell, as well as ESCC patients. 3'UTR luciferase assay was used to confirm the direct binding of miRNAs with TUSC2 and TUSC2P 3'UTR. Relation betweenTUSC2P, TUSC2 and ESCC prognosis was predicted by survival analysis (n = 56).

RESULTS

Pseudogene TUSC2P was down regulated in ESCC tissues compared with paired normal adjacent tissues, and the expression of TUSC2P was significantly correlated with survivalof ESCC patients. Over expression of TUSC2P in EC109 and TE-1 cells resulted in altered expression of TUSC2, thus inhibited proliferation, invasion and promoted apoptosis. Dual luciferase assay demonstrated that TUSC2P 3'UTR decoyed miR-17-5p, miR-520a-3p, miR-608, miR-661 from binding to TUSC2.

CONCLUSIONS

TUSC2P can suppresses the tumor function of esophageal squamous cell carcinoma by regulating TUSC2 expression and may also serve as a prognostic factor for ESCC patients.

Pan-cancer analysis of somatic mutations and transcriptomes reveals common functional gene clusters shared by multiple cancer types

To discover functional gene clusters across cancers, we performed a systematic pan-cancer analysis of 33 cancer types. We identified genes that were associated with somatic mutations and were the cores of a co-expression network. We found that multiple cancer types have relatively exclusive hub genes individually; however, the hub genes cooperate with each other based on their functional relationship. When we built a protein-protein interaction network of hub genes and found nine functional gene clusters across cancer types, the gene clusters divided not only the region of the network map, but also the function of the network by their distinct roles related to the development and progression of cancer. This functional relationship between the clusters and cancers was underpinned by the high expression of module genes and enrichment of programmed cell death, and known candidate cancer genes. In addition to protein-coding hub genes, non-coding hub genes had a possible relationship with cancer. Overall, our approach of investigating cancer genes enabled finding pan-cancer hub genes and common functional gene clusters shared by multiple cancer types based on the expression status of the primary tumour and the functional relationship of genes in the biological network.

PTEN/PTENP1: 'Regulating the regulator of RTK-dependent PI3K/Akt signalling', new targets for cancer therapy

Regulation of the PI-3 kinase (PI3K)/Akt signalling pathway is essential for maintaining the integrity of fundamental cellular processes, cell growth, survival, death and metabolism, and dysregulation of this pathway is implicated in the development and progression of cancers. Receptor tyrosine kinases (RTKs) are major upstream regulators of PI3K/Akt signalling. The phosphatase and tensin homologue (PTEN), a well characterised tumour suppressor, is a prime antagonist of PI3K and therefore a negative regulator of this pathway. Loss or inactivation of PTEN, which occurs in many tumour types, leads to overactivation of RTK/PI3K/Akt signalling driving tumourigenesis. Cellular PTEN levels are tightly regulated by a number of transcriptional, post-transcriptional and post-translational regulatory mechanisms. Of particular interest, transcription of the PTEN pseudogene, PTENP1, produces sense and antisense transcripts that exhibit post-transcriptional and transcriptional modulation of PTEN expression respectively. These additional levels of regulatory complexity governing PTEN expression add to the overall intricacies of the regulation of RTK/PI-3 K/Akt signalling. This review will discuss the regulation of oncogenic PI3K signalling by PTEN (the regulator) with a focus on the modulatory effects of the sense and antisense transcripts of PTENP1 on PTEN expression, and will further explore the potential for new therapeutic opportunities in cancer treatment.

Long non-coding RNA MT1DP shunts the cellular defense to cytotoxicity through crosstalk with MT1H and RhoC in cadmium stress

Metallothioneins (MTs) are known to protect cells against oxidative stress, especially providing protection against cadmium (Cd) toxicity in hepatocytes. There are various gene variants and pseudogenes for MTs; however, there is little understanding on the functions of those non-coding MT members that are known to be expressed as long non-coding RNAs (lncRNAs) nowadays. Different from most protein-coding MT members, MT1DP was here found that remarkably induced to provoke cytotoxicity in hepatocytes in response to Cd treatment. MT1DP exerted such a pro-apoptotic function in Cd-treated hepatocytes through interacting with two partners: RhoC and MT1H. On one hand, MT1DP interacted with RhoC protein to increase the latter’s stability by preventing lysosome-dependent protein degradation. Therefore, upon Cd stress, MT1DP/RhoC complex was quickly reinforced to activate RhoC-CCN1/2-AKT signaling and potentiate Ca²⁺ influx, leading to enhanced Cd uptake and elevated Cd toxicity. On the other hand, MT1H, a protein-coding member of the MT family with little known function, was found to quickly respond to Cd exposure along with MT1DP. Mechanistically, MT1H and MT1DP were uncovered to mutually protect each other through a reciprocal ceRNA mechanism, building up a positive feedback loop to enforce MT1DP-conducted signaling upon Cd exposure. Moreover, MT1DP was found to contribute much more to the activation of RhoC-CCN1/2-AKT signaling than MT1H. Considered together, we here unveiled a mystery whether a pseudogene within the MT family, MT1DP, has actual biological functions in regulating Cd-induced cellular defense. Our findings unearthed an important role of pseudogene MT1DP in calibrating the cellular machinery to switch the cellular defense to cytotoxicity through crosslinking an interplay between its two partners, namely MT1H and RhoC, under cadmium stress.

RetrogeneDB-a database of plant and animal retrocopies

For a long time, retrocopies were considered ‘junk DNA’, but numerous studies have shown that retrocopies may gain functionality and become so-called retrogenes. Retrogenes may code fully functional proteins that coexist with parental gene products or may even replace them. Retrocopies may also function as regulatory RNAs and, for example, become a source of small interfering RNAs, act as trans natural antisense transcripts or as alternative targets for miRNAs. Numerous researchers have emphasized that retrogenes play a crucial role in various organisms’ developmental stages and diseases. Despite the ever-growing evidence of the importance of retrocopies, resources dedicated to retroposition are very limited. Here, we report an update of the RetrogeneDB, which, to the best of our knowledge, is the largest database dedicated to retrocopies. It provides annotations of 86 458 retrocopies in 62 animal and 37 plant species. The database contains information about the retrocopies’ localization, open reading frame conservation, expression, RNA Polymerase II activity and the alternative transcription start site studies. Orthologous relationships between retrogenes were also determined, which made retrocopy conservation studies much more valuable. Additionally, based on the RNA-Seq data from the Geuvadis project, the expression levels of retrocopies were estimated in a total of 50 individuals from 5 human populations. The information is now presented in a new, more user-friendly web interface, with easy access to the source data, which may be used for the downstream analysis. RetrogeneDB is freely available at http://yeti.amu.edu.pl/retrogenedb.

Database URL:http://yeti.amu.edu.pl/retrogenedb

Secondary database URL:http://rhesus.amu.edu.pl/retrogenedb

Epithelial-to-mesenchymal transition in gallbladder cancer: from clinical evidence to cellular regulatory networks

Gallbladder cancer (GBC), with late diagnosis, rapid disease progression and early metastasis, is a highly aggressive malignant tumor found worldwide. Patients with GBC have poor survival, low curative resection rates and early recurrence. For such a lethal tumor, uncovering the mechanisms and exploring new strategies to prevent tumor progression and metastasis are critically important. Epithelial-to-mesenchymal transition (EMT) has a prominent role in the early steps of tumor progression and metastasis by initiating polarized epithelial cell transition into motile mesenchymal cells. Accumulating evidence suggests that EMT can be modulated by the cooperation of multiple mechanisms affecting common targets. Signaling pathways, transcriptional and post-transcriptional regulation and epigenetic alterations are involved in the stepwise EMT regulatory network in GBC. Loss of epithelial markers, acquisition of mesenchymal markers and dysregulation of EMT-inducing transcription factors (EMT-TFs) have been observed and are associated with the clinicopathology and prognosis of GBC patients. Therefore, EMT may be a detectable and predictable event for predicting GBC progression and metastasis in the clinic. In this review, we will provide an overview of EMT from the clinical evidence to cellular regulatory networks that have been studied thus far in clinical and basic GBC studies.

Pseudogenes of annexin A2, novel prognosis biomarkers for diffuse gliomas

Diffuse gliomas is a kind of common malignant primary brain tumor. Pseudogenes have multilayered biological function in the progression of human cancers. In this study, Differentially Expressed Pseudogenes (DEPs) between glioblastomas and non-tumor controls were found by bioinformatics analysis, of which the annexin A2 pseudogenes (ANXA2P1, ANXA2P2 and ANXA2P3) were significantly up-regulated, along with the parent gene annexin A2 (ANXA2). Among four glioblastoma subtypes, ANXA2P1 and ANXA2P2 were preferentially expressed in mesenchymal subtype and less expressed in proneural subtype. Meanwhile, Pearson’s correlation analysis revealed that the expression level of ANXA2 was positively correlated with ANXA2 pseudogenes expression. Then, the expression patterns of ANXA2 and its pseudogenes were validated in diffuse glioma specimens (n=99) and non-tumor tissues (n=12) by quantitative real-time PCR (qRT-PCR). Additionally, Kaplan–Meier analysis revealed that highly expressed ANXA2 and annexin A2 pseudogenes were associated with the poor survival outcome of glioma patients. Cox regression analyses suggested that ANXA2, ANXA2P1 and ANXA2P2 were the independent prognosis factors for gliomas. Furthermore, down-regulation of ANXA2 and ANXA2 pseudogenes might contribute to the improvement of patients’ survival who received chemotherapy and radiotherapy. These results demonstrated that ANXA2 pseudogenes and ANXA2 could be used as the novel biomarkers for diagnosis, prognosis and target therapy of gliomas.

The pseudogene derived from long non-coding RNA DUXAP10 promotes colorectal cancer cell growth through epigenetically silencing of p21 and PTEN

Recently, substantial evidence has demonstrated that pseudogene derived lncRNAs are crucial regulators of cancer development and progression. DUXAP10,a pseudogene derived long non-coding RNA(lncRNA), is overexpression in colorectal cancer (CRC), but its expression pattern, biological function and underlying mechanism in CRC is still undetermined. In this study, we observed that DUXAP10 was up-regulated in CRC tissues which was positively correlated with advanced pathological stages, larger tumor sizes and lymph node metastasis. Additionally, knockdown of DUXAP10 inhibited cell proliferation, induced cell apoptosis and increase the number of G0/G1 cells significantly in the HCT116 and SW480 cell lines. Moreover, DUXAP10 silencing inhibited tumor growth in vivo. Further mechanism study showed that, by binding to histone demethylase lysine-specific demethylase 1 (LSD1), DUXAP10 promote CRC cell growth and reduced cell apoptosis through silencing the expression of p21 and phosphatase and tensin homolog (PTEN) tumor suppressor. Our findings suggested that the pseudogene-derived from lncRNA DUXAP10 promotes the biological progression of CRC and is likely to be a potential therapeutic target for CRC intervention.

Identifying large indels in targeted next generation sequencing assays for myeloid neoplasms: a cautionary tale of the ZRSR1 pseudogene

Targeted next generation sequencing platforms have been increasingly utilised for identification of novel mutations in myeloid neoplasms, such as acute myeloid leukaemia (AML), and hold great promise for use in routine clinical diagnostics. In this study, we evaluated the utility of an open source variant caller in detecting large indels in a targeted sequencing of AML samples. While we found that this bioinformatics pipeline has the potential to accurately capture large indels (>20 bp) in patient samples, we highlighted the pitfall of a confounding ZRSR1 pseudogene that led to an erroneous ZRSR2 variant call. We further discuss possible clinical implications of the ZRSR1 pseudogene in myeloid neoplasms based on its molecular features. Knowledge of the confounding ZRSR1 pseudogene in ZRSR2 sequencing assays could be particularly important in AML diagnostics because the detection of ZRSR2 in AML patients is highly specific for an s-AML diagnosis.

Protein-Coding Genes' Retrocopies and Their Functions

Transposable elements, often considered to be not important for survival, significantly contribute to the evolution of transcriptomes, promoters, and proteomes. Reverse transcriptase, encoded by some transposable elements, can be used in trans to produce a DNA copy of any RNA molecule in the cell. The retrotransposition of protein-coding genes requires the presence of reverse transcriptase, which could be delivered by either non-long terminal repeat (non-LTR) or LTR transposons. The majority of these copies are in a state of “relaxed” selection and remain “dormant” because they are lacking regulatory regions; however, many become functional. In the course of evolution, they may undergo subfunctionalization, neofunctionalization, or replace their progenitors. Functional retrocopies (retrogenes) can encode proteins, novel or similar to those encoded by their progenitors, can be used as alternative exons or create chimeric transcripts, and can also be involved in transcriptional interference and participate in the epigenetic regulation of parental gene expression. They can also act in trans as natural antisense transcripts, microRNA (miRNA) sponges, or a source of various small RNAs. Moreover, many retrocopies of protein-coding genes are linked to human diseases, especially various types of cancer.

The pseudogene-derived long noncoding RNA SFTA1P is down-regulated and suppresses cell migration and invasion in lung adenocarcinoma

Pseudogenes were once considered to be genomic fossils without biological function. Interestingly, recent evidence showed that a lot of pseudogenes are transcribed in human cancers, and their alterations contribute to multiple cancer development and progression. It is apparent that many pseudogenes transcribe noncoding RNAs and contribute to the role noncoding genome plays in human cancers. On this basis, some pseudogene transcripts are currently ranked among the classes of long noncoding RNAs. In this study, we identified a new pseudogene-derived long noncoding RNA termed SFTA1P by analyzing the microarray data of non-small cell lung cancer from Gene Expression Omnibus datasets. We found that SFTA1P expression was significantly decreased in non-small cell lung cancer tissues compared with normal tissues in non-small cell lung cancer microarray data. Moreover, decreased SFTA1P expression is only correlated with lung adenocarcinoma patients' poor survival time but not with lung squamous cell carcinoma patients' survival. In addition, gain-of-function studies including growth curves, migration, invasion assays, and in vivo studies were performed to verify the tumor suppressor role of SFTA1P in non-small cell lung cancer. Finally, the potential underlying pathways involved in SFTA1P were investigated by analyzing the SFTA1P-correlated genes in The Cancer Genome Atlas lung adenocarcinoma and normal tissues RNA sequencing data. Taken together, these findings demonstrate that pseudogene-derived long noncoding RNA SFTA1P exerts the tumor suppressor functions in human lung adenocarcinoma. Our investigation reveals the novel roles of pseudogene in lung adenocarcinoma, which may serve as a new target for lung adenocarcinoma diagnosis and therapy.

The Pseudogene DUXAP8 Promotes Non-small-cell Lung Cancer Cell Proliferation and Invasion by Epigenetically Silencing EGR1 and RHOB

Recently, the non-protein-coding functional elements in the human genome have been identified as key regulators in postgenomic biology, and a large number of pseudogenes as well as long non-coding RNAs (lncRNAs) have been found to be transcribed in multiple human cancers. However, only a small proportion of these pseudogenes has been functionally characterized. In this study, we screened for pseudogenes associated with human non-small-cell lung cancer (NSCLC) by comparative analysis of several independent datasets from the GEO. We identified a transcribed pseudogene named DUXAP8 that is upregulated in tumor tissues. Patients with higher DUXAP8 expression exhibited shorter survival, suggesting DUXAP8 as a new candidate prognostic marker for NSCLC patients. Knockdown of DUXAP8 impairs cell growth, migration, and invasion, and induces apoptosis both in vitro and in vivo. Mechanistically, DUXAP8 represses the tumor suppressors EGR1 and RHOB by recruiting histone demethylase LSD1 and histone methyltransferase EZH2, thereby promoting cell proliferation, migration, and invasion. These findings indicate that the pseudogene DUXAP8 may act as an oncogene in NSCLC by silencing EGR1 and RHOB transcription by binding with EZH2 and LSD1, which may offer a novel therapeutic target for this disease.

Pseudogenes regulate parental gene expression via ceRNA network

The concept of competitive endogenous RNA (ceRNA) was first proposed by Salmena and colleagues. Evidence suggests that pseudogene RNAs can act as a 'sponge' through competitive binding of common miRNA, releasing or attenuating repression through sequestering miRNAs away from parental mRNA. In theory, ceRNAs refer to all transcripts such as mRNA, tRNA, rRNA, long non-coding RNA, pseudogene RNA and circular RNA, because all of them may become the targets of miRNA depending on spatiotemporal situation. As binding of miRNA to the target RNA is not 100% complementary, it is possible that one miRNA can bind to multiple target RNAs and vice versa. All RNAs crosstalk through competitively binding to miRNAvia miRNA response elements (MREs) contained within the RNA sequences, thus forming a complex regulatory network. The ratio of a subset of miRNAs to the corresponding number of MREs determines repression strength on a given mRNA translation or stability. An increase in pseudogene RNA level can sequester miRNA and release repression on the parental gene, leading to an increase in parental gene expression. A massive number of transcripts constitute a complicated network that regulates each other through this proposed mechanism, though some regulatory significance may be mild or even undetectable. It is possible that the regulation of gene and pseudogene expression occurring in this manor involves all RNAs bearing common MREs. In this review, we will primarily discuss how pseudogene transcripts regulate expression of parental genes via ceRNA network and biological significance of regulation.