Pseudogenes: newly discovered players in human cancer

Unraveling the molecular landscape of osteoarthritis: A comprehensive review focused on the role of non-coding RNAs

Osteoarthritis (OA) poses a significant global health challenge, with its prevalence anticipated to increase in the coming years. This review delves into the emerging molecular biomarkers in OA pathology, focusing on the roles of various molecules such as metabolites, noncoding RNAs (ncRNAs), including microRNAs (miRNAs) and long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs). Advances in omics technologies have transformed biomarker identification, enabling comprehensive analyses of the complex pathways involved in OA pathogenesis. Notably, ncRNAs, especially miRNAs and lncRNAs, exhibit dysregulated expression patterns in OA, presenting promising opportunities for diagnosis and therapy. Additionally, the intricate interplay between epigenetic modifications and OA progression highlights the regulatory role of epigenetics in gene expression dynamics. Genome-wide association studies have pinpointed key genes undergoing epigenetic changes, providing insights into the inflammatory processes and chondrocyte hypertrophy typical of OA. Understanding the molecular landscape of OA, including biomarkers and epigenetic mechanisms, holds significant potential for developing innovative diagnostic tools and therapeutic strategies for OA management.

The physical and evolutionary energy landscapes of devolved protein sequences corresponding to pseudogenes

Protein evolution is guided by structural, functional, and dynamical constraints ensuring organismal viability. Pseudogenes are genomic sequences identified in many eukaryotes that lack translational activity due to sequence degradation and thus over time have undergone "devolution." Previously pseudogenized genes sometimes regain their protein-coding function, suggesting they may still encode robust folding energy landscapes despite multiple mutations. We study both the physical folding landscapes of protein sequences corresponding to human pseudogenes using the Associative Memory, Water Mediated, Structure and Energy Model, and the evolutionary energy landscapes obtained using direct coupling analysis (DCA) on their parent protein families. We found that generally mutations that have occurred in pseudogene sequences have disrupted their native global network of stabilizing residue interactions, making it harder for them to fold if they were translated. In some cases, however, energetic frustration has apparently decreased when the functional constraints were removed. We analyzed this unexpected situation for Cyclophilin A, Profilin-1, and Small Ubiquitin-like Modifier 2 Protein. Our analysis reveals that when such mutations in the pseudogene ultimately stabilize folding, at the same time, they likely alter the pseudogenes' former biological activity, as estimated by DCA. We localize most of these stabilizing mutations generally to normally frustrated regions required for binding to other partners.

Subcellular RNA distribution and its change during human embryonic stem cell differentiation

The spatial localization of RNA within cells is closely related to its function and also involved in cell fate determination. However, the atlas of RNA distribution within cells and dynamic changes during the developmental process are largely unknown. In this study, five subcellular components, including cytoplasmic extract, membrane extract, soluble nuclear extract, chromatin-bound nuclear extract, and cytoskeletal extract, were isolated and the rules of subcellular RNA distribution in human embryonic stem cells (hESCs) and its change during hESC differentiation are summarized for the first time. The overall distribution patterns of coding and non-coding RNAs are revealed. Interestingly, some developmental genes are found to be transcribed but confined to the chromatin in undifferentiated hESC. Unexpectedly, alternative splicing and polyadenylation endow spatial heterogeneity among different isoforms of the same gene. Finally, the dynamic pattern of RNA distribution during hESC differentiation is characterized, which provides new clues for a comprehensive understanding hESC pluripotency and differentiation.

Pathological roles of miRNAs and pseudogene-derived lncRNAs in human cancers, and their comparison as prognosis/diagnosis biomarkers

This review examines and compares the diagnostic and prognostic capabilities of miRNAs and lncRNAs derived from pseudogenes in cancer patients. Additionally, it delves into their roles in cancer pathogenesis. Both miRNAs and pseudogene-derived lncRNAs have undergone thorough investigation as remarkably sensitive and specific cancer biomarkers, offering significant potential for cancer detection and monitoring. . Extensive research is essential to gain a complete understanding of the precise roles these non-coding RNAs play in cancer, allowing the development of novel targeted therapies and biomarkers for improved cancer detection and treatment approaches.

LncRNA FTH1P3: A New Biomarker for Cancer-Related Therapeutic Development

Cancer is a persistent and urgent health problem that affects the entire world. Not long ago, regulatory biomolecules referred to as long noncoding RNAs (lncRNAs) might have value for their innate abundance and stability. These single-stranded RNAs potentially interfere with several physiological and biochemical cellular processes involved in many human pathological situations, particularly cancer diseases. Ferritin heavy chain1 pseudogene 3 (FTH1P3), a lncRNA that is ubiquitously transcribed and belongs to the ferritin heavy chain (FHC) family, represents a novel class of lncRNAs primarily found in oral squamous cell carcinoma. Further research has shown that FTH1P3 is involved in other malignancies such as uveal melanoma, glioma, esophageal squamous cell carcinoma, non-small cell lung cancer, breast cancer, laryngeal squamous cell carcinoma, and cervical cancer. Accordingly, FTH1P3 significantly enhances cancer symptoms, including cell proliferation, invasion, metastasis, chemoresistance, and inhibition of apoptosis through many specific mechanisms. Notably, the clinical data significantly demonstrated the association of FTH1P3 overexpression with poor prognosis and poor overall survival within the examined samples. Here, we summarize all the research published to date (13 articles) on FTH1P3, focusing on the biological function underlying the regulatory mechanism and its possible clinical relevance.

Prediction of Tumor Microenvironment Characteristics and Treatment Response in Lung Squamous Cell Carcinoma by Pseudogene OR7E47P-related Immune Genes

OBJECTIVE

Pseudogenes are initially regarded as nonfunctional genomic sequences, but some pseudogenes regulate tumor initiation and progression by interacting with other genes to modulate their transcriptional activities. Olfactory receptor family 7 subfamily E member 47 pseudogene (OR7E47P) is expressed broadly in lung tissues and has been identified as a positive regulator in the tumor microenvironment (TME) of lung adenocarcinoma (LUAD). This study aimed to elucidate the correlation between OR7E47P and tumor immunity in lung squamous cell carcinoma (LUSC).

METHODS

Clinical and molecular information from The Cancer Genome Atlas (TCGA) LUSC cohort was used to identify OR7E47P-related immune genes (ORIGs) by weighted gene correlation network analysis (WGCNA). Based on the ORIGs, 2 OR7E47P clusters were identified using non-negative matrix factorization (NMF) clustering, and the stability of the clustering was tested by an extreme gradient boosting classifier (XGBoost). LASSO-Cox and stepwise regressions were applied to further select prognostic ORIGs and to construct a predictive model (ORPScore) for immunotherapy. The Botling cohorts and 8 immunotherapy cohorts (the Samstein, Braun, Jung, Gide, IMvigor210, Lauss, Van Allen, and Cho cohorts) were included as independent validation cohorts.

RESULTS

OR7E47P expression was positively correlated with immune cell infiltration and enrichment of immune-related pathways in LUSC. A total of 57 ORIGs were identified to classify the patients into 2 OR7E47P clusters (Cluster 1 and Cluster 2) with distinct immune, mutation, and stromal programs. Compared to Cluster 1, Cluster 2 had more infiltration by immune and stromal cells, lower mutation rates of driver genes, and higher expression of immune-related proteins. The clustering performed well in the internal and 5 external validation cohorts. Based on the 7 ORIGs (HOPX, STX2, WFS, DUSP22, SLFN13, GGCT, and CCSER2), the ORPScore was constructed to predict the prognosis and the treatment response. In addition, the ORPScore was a better prognostic factor and correlated positively with the immunotherapeutic response in cancer patients. The area under the curve values ranged from 0.584 to 0.805 in the 6 independent immunotherapy cohorts.

CONCLUSION

Our study suggests a significant correlation between OR7E47P and TME modulation in LUSC. ORIGs can be applied to molecularly stratify patients, and the ORPScore may serve as a biomarker for clinical decision-making regarding individualized prognostication and immunotherapy.

A novel peptide P1-121aa encoded by STK24P1 regulates vasculogenic mimicry via ELF2 phosphorylation in glioblastoma

Glioblastoma (GBM) is the most common malignant tumor of the central nervous system. Vasculogenic mimicry (VM) is a hematological system composed of tumor cells that exert blood perfusion without relying on vascular endothelial cells. The current poor results of anti-vascular therapy for clinical GBM are associated with the presence of VM; therefore, it is important to investigate VM formation in GBM. Our results demonstrate that STK24P1 encodes P1-121aa with a kinase structural domain, and in vitro kinase assays demonstrated that P1-121aa mediates modification of ELF2 phosphorylation. ChIP and dual luciferase reporter gene assays demonstrated that the transcription factor ELF2 binds to VE-cadherin and the VEGFR2 promoter region, thereby promoting VM formation in glioma cells. P1-121aa, encoded by the pseudogene STK24P1, phosphorylates ELF2 at S107, increasing the stability of the ELF2 protein. ELF2 promotes VEGFR2 and VE-cadherin expression at the transcriptional level, which in turn promotes VM in GBM. This study demonstrates the important roles of STK24P1, P1-121aa, and ELF2 in regulating VM in GBM, which could provide potential targets for GBM treatment.

Pseudogenes in Cancer: State of the Art

Pseudogenes are duplicates of protein-coding genes that have accumulated multiple detrimental alterations, rendering them unable to produce the protein they encode. Initially disregarded as "junk DNA" due to their perceived lack of functionality, research on their biological roles has been hindered by this assumption. Nevertheless, recent focus has shifted towards these molecules due to their abnormal expression in cancer phenotypes. In this review, our objective is to provide a thorough overview of the current understanding of pseudogene formation, the mechanisms governing their expression, and the roles they may play in promoting tumorigenesis.

Multiple ubiquitin E3 ligase genes antagonistically regulate chloroplast-associated protein degradation

The chloroplast is the most prominent member of a diverse group of plant organelles called the plastids, and it is characterized by its vital role in photosynthesis. ^1,2,3 Most of the ∼3,000 different proteins in chloroplasts are synthesized in the cytosol in precursor (preprotein) form, each with a cleavable transit peptide. ^4,5,6,7,8 Preproteins are imported via translocons in the outer and inner envelope membranes of the chloroplast, termed TOC and TIC, respectively. ^{9,10,11,12,13} Discovery of the chloroplast-localized ubiquitin E3 ligase SUPPRESSOR OF PPI1 LOCUS1 (SP1) demonstrated that the nucleocytosolic ubiquitin-proteasome system (UPS) targets the TOC apparatus to dynamically control protein import and chloroplast biogenesis in response to developmental and environmental cues. The relevant UPS pathway is termed chloroplast-associated protein degradation (CHLORAD). ^14,15,16 Two homologs of SP1 exist, SP1-like1 (SPL1) and SPL2, but their roles have remained obscure. Here, we show that SP1 is ubiquitous in the Viridiplantae and that SPL2 and SPL1 appeared early during the evolution of the Viridiplantae and land plants, respectively. Through genetic and biochemical analysis, we reveal that SPL1 functions as a negative regulator of SP1, potentially by interfering with its ability to catalyze ubiquitination. In contrast, SPL2, the more distantly related SP1 homolog, displays partial functional redundancy with SP1. Both SPL1 and SPL2 modify the extent of leaf senescence, like SP1, but do so in diametrically opposite ways. Thus, SPL1 and SPL2 are bona fide CHLORAD system components with negative and positive regulatory functions that allow for nuanced control of this vital proteolytic pathway.

Pseudogenes and the associated ceRNA network as potential prognostic biomarkers for colorectal cancer

Colorectal cancer (CRC) is one of the most common and malignant carcinomas. Many long noncoding RNAs (lncRNAs) have been reported to play important roles in the tumorigenesis of CRC by influencing the expression of some mRNAs via competing endogenous RNA (ceRNA) networks and interacting with miRNAs. Pseudogene is one kind of lncRNA and can act as RNA sponges for miRNAs and regulate gene expression via ceRNA networks. However, there are few studies about pseudogenes in CRC. In this study, 31 differentially expressed (DE) pseudogenes, 17 DE miRNAs and 152 DE mRNAs were identified by analyzing the expression profiles of colon adenocarcinoma obtained from The Cancer Genome Atlas. A ceRNA network was constructed based on these RNAs. Kaplan-Meier analysis showed that 7 pseudogenes, 4 miRNAs and 30 mRNAs were significantly associated with overall survival. Then multivariate Cox regression analysis of the ceRNA-related DE pseudogenes was performed and a 5-pseudogene signature with the greatest prognostic value for CRC was identified. Moreover, the results were validated by the Gene Expression Omnibus database, and quantitative real-time PCR in 113 pairs of CRC tissues and colon cancer cell lines. This study provides a pseudogene-associated ceRNA network, 7 prognostic pseudogene biomarkers, and a 5-pseudogene prognostic risk signature that may be useful for predicting the survival of CRC patients.

Meta-Analytic Comparison of Global RNA Transcriptomes of Acute and Chronic Myeloid Leukemia Cells Reveals Novel Gene Candidates Governing Myeloid Malignancies

Simple Summary

Despite advances in the understanding of genetic risk factors and molecular mechanisms underlying acute myeloid leukemia (AML) and chronic myeloid leukemia (CML), clinical outcomes of current therapies in terms of disease relapse and mortality rate pose a great economic and social burden. To overcome this, the identification of new molecular prognostic biomarkers and pharmacological targets is crucial. Recent studies have suggested that AML and CML may share common pathogenic mechanisms and cellular substrates. To this end, in the present study, global transcriptome profiles of AML and CML at the molecular and cellular level were directly compared using a combination of meta-analysis and modern statistics, and novel candidate genes and specific biological processes associated with the pathogenesis of AML and CML were characterized. Our study significantly improves our current understanding of myeloid leukemia and will help develop new therapeutic targets and biomarkers for disease progression, management and treatment response.

Abstract

Background: Acute myeloid leukemia (AML) and chronic myeloid leukemia (CML) represent a group of hematological malignancies characterized by the pathogenic clonal expansion of leukemic myeloid cells. The diagnosis and clinical outcome of AML and CML are complicated by genetic heterogeneity of disease; therefore, the identification of novel molecular biomarkers and pharmacological targets is of paramount importance. Methods: RNA-seq-based transcriptome data from a total of five studies were extracted from NCBI GEO repository and subjected to an in-depth bioinformatics analysis to identify differentially expressed genes (DEGs) between AML and CML. A systemic literature survey and functional gene ontology (GO) enrichment analysis were performed for the top 100 DEGs to identify novel candidate genes and biological processes associated with AML and CML. Results: LINC01554, PTMAP12, LOC644936, RPS27AP20 and FAM133CP were identified as novel risk genes for AML and CML. GO enrichment analysis showed that DEGs were significantly associated with pre-RNA splicing, reactive oxygen species and glycoprotein metabolism, the cellular endomembrane system, neutrophil migration and antimicrobial immune response. Conclusions: Our study revealed novel biomarkers and specific biological processes associated with AML and CML. Further studies are required to evaluate their value as molecular targets for managing and treating the myeloid malignancies.

Mechanisms Controlling MicroRNA Expression in Tumor

MicroRNAs (miRNAs) are widely present in many organisms and regulate the expression of genes in various biological processes such as cell differentiation, metabolism, and development. Numerous studies have shown that miRNAs are abnormally expressed in tumor tissues and are closely associated with tumorigenesis. MiRNA-based cancer gene therapy has consistently shown promising anti-tumor effects and is recognized as a new field in cancer treatment. So far, some clinical trials involving the treatment of malignancies have been carried out; however, studies of miRNA-based cancer gene therapy are still proceeding slowly. Therefore, furthering our understanding of the regulatory mechanisms of miRNA can bring substantial benefits to the development of miRNA-based gene therapy or other combination therapies and the clinical outcome of patients with cancer. Recent studies have revealed that the aberrant expression of miRNA in tumors is associated with promoter sequence mutation, epigenetic alteration, aberrant RNA modification, etc., showing the complexity of aberrant expression mechanisms of miRNA in tumors. In this paper, we systematically summarized the regulation mechanisms of miRNA expression in tumors, with the aim of providing assistance in the subsequent elucidation of the role of miRNA in tumorigenesis and the development of new strategies for tumor prevention and treatment.

The pseudogenes of eukaryotic translation elongation factors (EEFs): Role in cancer and other human diseases

The eukaryotic translation elongation factors (EEFs), i.e. EEF1A1, EEF1A2, EEF1B2, EEF1D, EEF1G, EEF1E1 and EEF2, are coding-genes that play a central role in the elongation step of translation but are often altered in cancer. Less investigated are their pseudogenes. Recently, it was demonstrated that pseudogenes have a key regulatory role in the cell, especially via non-coding RNAs, and that the aberrant expression of ncRNAs has an important role in cancer development and progression. The present review paper, for the first time, collects all that published about the EEFs pseudogenes to create a base for future investigations. For most of them, the studies are in their infancy, while for others the studies suggest their involvement in normal cell physiology but also in various human diseases. However, more investigations are needed to understand their functions in both normal and cancer cells and to define which can be useful biomarkers or therapeutic targets.

A Novel Pseudogene Methylation Signature to Predict Temozolomide Outcome in Non-G-CIMP Glioblastomas

Objective

Alterations in the methylation state of pseudogenes may serve as clinically useful biomarkers of glioblastomas (GBMs) that do not have glioma-CpG island methylator phenotype (G-CIMP).

Methods

Non-G-CIMP GBM datasets were included for evaluation, and a RISK-score signature was determined from the methylation state of pseudogene loci. Both bioinformatic and experimental analyses were performed for biological validation.

Results

By integrating clinical information with DNA methylation microarray data, we screened a panel of eight CpGs from discovery cohorts of non-G-CIMP GBMs. Each CpG could accurately and independently predict the prognosis of patients under a treatment regime that combined radiotherapy (RT) and temozolomide (TMZ). The 8-CpG signature appeared to show opposite prognostic correlations between patients treated with RT/TMZ and those treated with RT monotherapy. The analyses further indicated that this signature had predictive value for TMZ efficacy because different survival benefits between RT/TMZ and RT therapies were observed in each risk subgroup. The incorporation of other risk factors, such as age and O-6-methylguanine-DNA methyltransferase (MGMT) promoter methylation status, with our pseudogene methylation signature could provide precise risk classification. In vitro experimental data revealed that two locus-specific pseudogenes (ZNF767P and CLEC4GP1) may modulate TMZ resistance via distinct mechanisms in GBM cells.

Conclusion

The biologically and clinically relevant RISK-score signature, based on pseudogene methylation loci, may offer information for predicting TMZ responses of non-G-CIMP GBMs, that is independent from, but complementary to, MGMT-based approaches.

Pseudogene AK4P1 promotes pancreatic ductal adenocarcinoma progression through relieving miR-375-mediated YAP1 degradation

Pseudogenes have been reported to play oncogenic or tumor-suppressive roles in cancer progression. However, the molecular mechanism of most pseudogenes in pancreatic ductal adenocarcinoma (PDAC) remains unknown. Herein, we characterized a novel pseudogene-miRNA-mRNA network associated with PDAC progression using bioinformatics analysis. After screening by dreamBase and GEPIA, 12 up-regulated and 7 down-regulated differentially expressed pseudogenes (DEPs) were identified. According to survival analysis, only elevated AK4P1 indicated a poor prognosis for PDAC patients. Moreover, we found that AK4 acts as a cognate gene of AK4P1 and also predicts worse survival for PDAC patients. Furthermore, 32 miRNAs were predicted to bind to AK4P1 by starBase, among which miR-375 was identified as the most potential binding miRNA of AK4P1. A total of 477 potential target genes of miR-375 were obtained by miRNet, in which 49 hub genes with node degree ≥ 20 were identified by STRING. Subsequent analysis for hub genes demonstrated that YAP1 may be a functional downstream target of AK4P1. To confirmed the above findings, microarray, and qRT-PCR assay revealed that YAP1 was dramatically upregulated in both PDAC cells and tissues. Functional experiments showed that knockdown of YAP1 significantly suppressed PDAC cells growth, increased apoptosis, and decreased the ability of invasion. In conclusion, amplification of AK4P1 may fuel the onset and development of PDAC by targeting YAP1 through competitively binding to miR-375, and serve as a promising biomarker and therapeutic target for PDAC.

Virtual Gene Concept and a Corresponding Pragmatic Research Program in Genetical Data Science

Mendel proposed an experimentally verifiable paradigm of particle-based heredity that has been influential for over 150 years. The historical arguments have been reflected in the near past as Mendel's concept has been diversified by new types of omics data. As an effect of the accumulation of omics data, a virtual gene concept forms, giving rise to genetical data science. The concept integrates genetical, functional, and molecular features of the Mendelian paradigm. I argue that the virtual gene concept should be deployed pragmatically. Indeed, the concept has already inspired a practical research program related to systems genetics. The program includes questions about functionality of structural and categorical gene variants, about regulation of gene expression, and about roles of epigenetic modifications. The methodology of the program includes bioinformatics, machine learning, and deep learning. Education, funding, careers, standards, benchmarks, and tools to monitor research progress should be provided to support the research program.

ceRNAs in Cancer: Mechanism and Functions in a Comprehensive Regulatory Network

Noncoding RNAs have been shown with powerful ability in post-transcriptional regulation, enabling intertwined RNA crosstalk and global molecular interaction in a large amount of dysfunctional conditions including cancer. Competing endogenous RNAs (ceRNAs) are those competitively binding with shared microRNAs (miRNAs), freeing their counterparts from miRNA-induced degradation, thus actively influencing and connecting with each other. Constantly updated analytical approaches boost outstanding advancement achieved in this burgeoning hotspot in multilayered intracellular communication, providing new insights into pathogenesis and clinical treatment. Here, we summarize the mechanisms and correlated factors under this RNA interplay and deregulated transcription profile in neoplasm and tumor progression, underscoring the great significance of ceRNAs for diagnostic values, monitoring biomarkers, and prognosis evaluation in cancer.

Pseudogene Transcripts in Head and Neck Cancer: Literature Review and In Silico Analysis

Once considered nonfunctional, pseudogene transcripts are now known to provide valuable information for cancer susceptibility, including head and neck cancer (HNC), a serious health problem worldwide, with about 50% unimproved overall survival over the last decades. The present review focuses on the role of pseudogene transcripts involved in HNC risk and prognosis. We combined current literature and in silico analyses from The Cancer Genome Atlas (TCGA) database to identify the most deregulated pseudogene transcripts in HNC and their genetic variations. We then built a co-expression network and performed gene ontology enrichment analysis to better understand the pseudogenes’ interactions and pathways in HNC. In the literature, few pseudogenes have been studied in HNC. Our in silico analysis identified 370 pseudogene transcripts associated with HNC, where SPATA31D5P, HERC2P3, SPATA31C2, MAGEB6P1, SLC25A51P1, BAGE2, DNM1P47, SPATA31C1, ZNF733P and OR2W5 were found to be the most deregulated and presented several genetic alterations. NBPF25P, HSP90AB2P, ZNF658B and DPY19L2P3 pseudogenes were predicted to interact with 12 genes known to participate in HNC, DNM1P47 was predicted to interact with the TP53 gene, and HLA-H pseudogene was predicted to interact with HLA-A and HLA-B genes. The identified pseudogenes were associated with cancer biology pathways involving cell communication, response to stress, cell death, regulation of the immune system, regulation of gene expression, and Wnt signaling. Finally, we assessed the prognostic values of the pseudogenes with the Kaplan–Meier Plotter database, and found that expression of SPATA31D5P, SPATA31C2, BAGE2, SPATA31C1, ZNF733P and OR2W5 pseudogenes were associated with patients’ survival. Due to pseudogene transcripts’ potential for cancer diagnosis, progression, and as therapeutic targets, our study can guide new research to HNC understanding and development of new target therapies.

Studying the Oncosuppressive Functions of PTENP1 as a ceRNA

PTENP1 is a processed pseudogene of the tumour suppressor phosphatase and tensin homolog deleted on chromosome 10 (PTEN). It functions posttranscriptionally to regulate PTEN by acting as a sponge for microRNAs that target PTEN. PTENP1 therefore functions as a competitive endogenous RNA (ceRNA), competing with PTEN for binding of microRNAs (miRNA) and thereby modulating PTEN cellular abundance. Studies of the overexpression of PTENP1 all confirm its oncosuppressive function to be mediated through the suppression of cell proliferation, induction of apoptosis, and inhibition of cell migration and invasion of cancer cells of differing types. These oncosuppressive functions are a direct consequence of miRNA binding by PTENP1 and the subsequent liberation of PTEN from miRNA induced suppression. In this chapter, we will focus initially on the description of a high efficiency transient transfection method to introduce and overexpress PTENP1 in the cell type of interest, followed by accurate methodologies to measure transfection efficiency by flow cytometry. We will then continue to describe two methods to analyze cell proliferation, namely the CCK-8 assay and Click-iT^® EdU assay. Due to commonalities in the manifestation of the oncosuppressive effects of PTENP1, mediated through its role as a ceRNA, the methods presented in this chapter will have wide applicability to a variety of different cell types.

Pseudogene legumain promotes thyroid carcinoma progression via the microRNA-495/autophagy pathway

The pseudogene legumain (LGMN) has been reported to regulate cancer cell biology. However, the role of LGMN in thyroid carcinoma remains unknown. Herein, Cell Counting Kit 8 and Transwell assays were performed to evaluate cellular proliferation and invasion capacity, respectively. In addition, a tube formation assay was performed to assess HUVEC angiogenesis. The results showed that LGMN depletion attenuated cellular proliferation, invasion and tube formation ability, and that LGMN expression was dysregulated in thyroid carcinoma tumors. Furthermore, patients with high LGMN expression levels exhibited a lower overall survival rate than those with low expression levels. LGMN and microRNA (miR)-495 modulated the expression levels of autophagy-related gene 3 (ATG3) and p62. Finally, ATG3 overexpression rescued the LGMN-regulated thyroid carcinoma phenotype. In conclusion, LGMN was found to promote thyroid carcinoma progression via the miR-495/autophagy axis, thus providing novel insights for understanding the pathogenesis of thyroid carcinoma.

Pseudogene ACTBP2 increases blood-brain barrier permeability by promoting KHDRBS2 transcription through recruitment of KMT2D/WDR5 in Aβ1-42 microenvironment

The blood–brain barrier (BBB) has a vital role in maintaining the homeostasis of the central nervous system (CNS). Changes in the structure and function of BBB can accelerate Alzheimer’s disease (AD) development. β-Amyloid (Aβ) deposition is the major pathological event of AD. We elucidated the function and possible molecular mechanisms of the effect of pseudogene ACTBP2 on the permeability of BBB in Aβ_1–42 microenvironment. BBB model treated with Aβ_1–42 for 48 h were used to simulate Aβ-mediated BBB dysfunction in AD. We proved that pseudogene ACTBP2, RNA-binding protein KHDRBS2, and transcription factor HEY2 are highly expressed in ECs that were obtained in a BBB model in vitro in Aβ_1–42 microenvironment. In Aβ_1–42-incubated ECs, ACTBP2 recruits methyltransferases KMT2D and WDR5, binds to KHDRBS2 promoter, and promotes KHDRBS2 transcription. The interaction of KHDRBS2 with the 3′UTR of HEY2 mRNA increases the stability of HEY2 and promotes its expression. HEY2 increases BBB permeability in Aβ_1–42 microenvironment by transcriptionally inhibiting the expression of ZO-1, occludin, and claudin-5. We confirmed that knocking down of Khdrbs2 or Hey2 increased the expression levels of ZO-1, occludin, and claudin-5 in APP/PS1 mice brain microvessels. ACTBP2/KHDRBS2/HEY2 axis has a crucial role in the regulation of BBB permeability in Aβ_1–42 microenvironment, which may provide a novel target for the therapy of AD.

Regulatory Potential of Competing Endogenous RNAs in Myotonic Dystrophies

Non-coding RNAs (ncRNAs) have been reported to be implicated in cell fate determination and various human diseases. All ncRNA molecules are emerging as key regulators of diverse cellular processes; however, little is known about the regulatory interaction among these various classes of RNAs. It has been proposed that the large-scale regulatory network across the whole transcriptome is mediated by competing endogenous RNA (ceRNA) activity attributed to both protein-coding and ncRNAs. ceRNAs are considered to be natural sponges of miRNAs that can influence the expression and availability of multiple miRNAs and, consequently, the global mRNA and protein levels. In this review, we summarize the current understanding of the role of ncRNAs in two neuromuscular diseases, myotonic dystrophy type 1 and 2 (DM1 and DM2), and the involvement of expanded CUG and CCUG repeat-containing transcripts in miRNA-mediated RNA crosstalk. More specifically, we discuss the possibility that long repeat tracts present in mutant transcripts can be potent miRNA sponges and may affect ceRNA crosstalk in these diseases. Moreover, we highlight practical information related to innovative disease modelling and studying RNA regulatory networks in cells. Extending knowledge of gene regulation by ncRNAs, and of complex regulatory ceRNA networks in DM1 and DM2, will help to address many questions pertinent to pathogenesis and treatment of these disorders; it may also help to better understand general rules of gene expression and to discover new rules of gene control.

The lnc-CTSLP8 upregulates CTSL1 as a competitive endogenous RNA and promotes ovarian cancer metastasis

Background

Ovarian cancer is highly lethal and has a poor prognosis due to metastasis. Long non-coding RNAs (lncRNAs) are key regulators of tumor development, but their role in ovarian cancer metastasis remains unclear.

Methods

The expression of lnc-CTSLP8 in ovarian cancer was analyzed in public databases (TCGA and GEO) and validated via qRT-PCR. Lnc-CTSLP8 overexpression and knockout cell lines were constructed using a lentiviral vector and the CRISP/Cas9 system. Cell proliferation, colony formation, migration, and invasion were analyzed. An ovarian orthotopic tumor mouse model was used for the in vivo study. Changes in autophagosomes, autolysosomes, and mitochondria in ovarian cancer cells were observed via transmission electron microscopy. EMT markers were detected by immunoblotting and immunofluorescence assays. RNA immunoprecipitation, RNA pull-down, and dual luciferase reporter assays were performed to confirm the interaction between lnc-CTSLP8 and miR-199a-5p.

Results

A novel pseudogene, lnc-CTSLP8, was identified in ovarian cancer, with significantly elevated expression in metastatic tumor tissues compared to primary ovarian tumors. When overexpressed, lnc-CTSLP8 promoted ovarian cancer in vitro and in vivo by acting as a sponge for miR-199a-5p. Autophagy and EMT in ovarian cancer were also enhanced by lnc-CTSLP8. Mechanistically, lnc-CTSLP8 upregulated CTSL1 as a competitive endogenous RNA and exhibited oncogenic effects. Moreover, CTSL1 inhibitor treatment and miR-199a-5p overexpression abrogated the effects of lnc-CTSLP8 overexpression.

Conclusions

lnc-CTSLP8 acts as a ceRNA in ovarian cancer and represents a potential therapeutic target for metastatic ovarian cancer.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13046-021-01957-z.

Identification of key pseudogenes in nasopharyngeal carcinoma based on RNA-Seq analysis

BACKGROUND

Nasopharyngeal carcinoma (NPC) is a malignant head and neck tumor, and more than 70% of new cases are in East and Southeast Asia. However, association between NPC and pseudogenes playing important roles in genesis of multiple tumor types is still not clear and needs to be investigated.

METHODS

Using RNA-Sequencing (RNA-seq) technology, we analyzed pseudogene expression in 13 primary NPC and 6 recurrent NPC samples as well as their paracancerous counterparts. Quantitative PCR was used to validate the differentially expressed pseudogenes.

RESULTS

We found 251 differentially expressed pseudogenes including 73 up-regulated and 178 down-regulated ones between primary NPC and paracancerous tissues. Enrichment analysis of gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway were conducted to filter out the key pseudogenes. We reported that pseudogenes from cytochrome P450 (CYP) family, such as CYP2F2P, CYP2G1P, CYP4F24P, CYP2B7P and CYP2G2P were significantly down-regulated in NPC compared to paracancerous tissues, while IGHV1OR15-2, IGHV3-11, FCGR1CP and IGHV3-69-1 belonging to Fc gamma receptors were significantly up-regulated. CYP2B7P, CYP2F2P and CYP4F26P were enriched in arachidonic acid metabolism pathway. The qRT-PCR analysis validated the lower expression of pseudogenes CYP2F2P and CYP2B7P in NPC tissues and cell lines compared to paracancerous tissues and normal human nasopharyngeal epithelial cell line. CYP2B7P overexpression weakened migratory and invasive capacity of NPC cell line. Moreover, the expression pattern of those pseudogenes in recurrent NPC tissues was different from the primary NPC.

CONCLUSION

This study suggested the role of pseudogenes in tumorigenesis and progression, potentially functioning as therapeutic targets to NPC.

Identification of dysregulated competing endogenous RNA networks in glioblastoma: A way toward improved therapeutic opportunities

Glioblastoma is recognized as one of the leading causes of death worldwide. Although there have been considerable advancements in understanding the causative molecular mechanisms of this malignancy, effective therapeutic strategies are still in limited use. It has been revealed that non-coding RNAs (ncRNAs) play critical roles in glioblastoma development, while interactions between the regulatory molecules such as long ncRNAs (lncRNAs), microRNAs (miRNAs), transcribed pseudogenes, and circular RNAs (circRNAs) remain to be fully deciphered. Over the recent years, researchers have discovered a new category of RNA molecules called competing endogenous RNA (ceRNA). This kind of RNA can contribute to molecular interactions in the form of ceRNA networks (ceRNETs). Multiple lines of evidence have demonstrated that dysregulation of various ceRNA networks is involved in glioblastoma development. Therefore, gaining insights into these dysregulations might offer potential for the early diagnosis of glioblastoma patients and identification of efficient therapeutic targets. In this review, we provide an overview of recent discoveries on ceRNA networks and the involvement of dysregulated networks in posing limitations to temozolomide therapy. We also describe signaling pathways relevant to the progression of glioblastoma.

Non-Coding RNAs in Cancer Diagnosis and Therapy: Focus on Lung Cancer

Over the last several decades, clinical evaluation and treatment of lung cancers have largely improved with the classification of genetic drivers of the disease, such as EGFR, ALK, and ROS1. There are numerous regulatory factors that exert cellular control over key oncogenic pathways involved in lung cancers. In particular, non-coding RNAs (ncRNAs) have a diversity of regulatory roles in lung cancers such that they have been shown to be involved in inducing proliferation, suppressing apoptotic pathways, increasing metastatic potential of cancer cells, and acquiring drug resistance. The dysregulation of various ncRNAs in human cancers has prompted preclinical studies examining the therapeutic potential of restoring and/or inhibiting these ncRNAs. Furthermore, ncRNAs demonstrate tissue-specific expression in addition to high stability within biological fluids. This makes them excellent candidates as cancer biomarkers. This review aims to discuss the relevance of ncRNAs in cancer pathology, diagnosis, and therapy, with a focus on lung cancer.

Structural characterization and duplication modes of pseudogenes in plants

We identified and characterized the pseudogene complements of five plant species: four dicots (Arabidopsis thaliana, Vitis vinifera, Populus trichocarpa and Phaseolus vulgaris) and one monocot (Oryza sativa). Retroposition was considered of modest importance for pseudogene formation in all investigated species except V. vinifera, which showed an unusually high number of retro-pseudogenes in non coding genic regions. By using a pipeline for the classification of sequence duplicates in plant genomes, we compared the relative importance of whole genome, tandem, proximal, transposed and dispersed duplication modes in the pseudo and functional gene complements. Pseudogenes showed higher tendencies than functional genes to genomic dispersion. Dispersed pseudogenes were prevalently fragmented and showed high sequence divergence at flanking regions. On the contrary, those deriving from whole genome duplication were proportionally less than expected based on observations on functional loci and showed higher levels of flanking sequence conservation than dispersed pseudogenes. Pseudogenes deriving from tandem and proximal duplications were in excess compared to functional loci, probably reflecting the high evolutionary rate associated with these duplication modes in plant genomes. These data are compatible with high rates of sequence turnover at neutral sites and double strand break repairs mediated duplication mechanisms.

Genome-wide analysis of pseudogenes reveals HBBP1's human-specific essentiality in erythropoiesis and implication in β-thalassemia

The human genome harbors 14,000 duplicated or retroposed pseudogenes. Given their functionality as regulatory RNAs and low conservation, we hypothesized that pseudogenes could shape human-specific phenotypes. To test this, we performed co-expression analyses and found that pseudogene exhibited tissue-specific expression, especially in the bone marrow. By incorporating genetic data, we identified a bone-marrow-specific duplicated pseudogene, HBBP1 (η-globin), which has been implicated in β-thalassemia. Extensive functional assays demonstrated that HBBP1 is essential for erythropoiesis by binding the RNA-binding protein (RBP), HNRNPA1, to upregulate TAL1, a key regulator of erythropoiesis. The HBBP1/TAL1 interaction contributes to a milder symptom in β-thalassemia patients. Comparative studies further indicated that the HBBP1/TAL1 interaction is human-specific. Genome-wide analyses showed that duplicated pseudogenes are often bound by RBPs and less commonly bound by microRNAs compared with retropseudogenes. Taken together, we not only demonstrate that pseudogenes can drive human evolution but also provide insights on their functional landscapes.

Strategies to Study the Functions of Pseudogenes in Mouse Models of Cancer

Aberrant expression of pseudogenes has been observed in many cancer types. Deregulated pseudogenes engage in a multitude of biological processes at the DNA, RNA, and protein levels and eventually facilitate disease progression. To investigate pseudogene functions in cancer, cell lines and cell line transplantation models have been widely used. However, cancer biology is best studied in the context of an intact organism. Here, we present various strategies to investigate pseudogenes in genetically engineered mouse models and discuss advantages and disadvantages of the different approaches.

Global gene methylation profiling of common warts caused by human papillomaviruses infection

Infection with the human papillomaviruses (HPV) often involves the epigenetic modification of the host genome. Despite its prevalence among the population, host genome methylation in HPV-induced warts is not clearly understood. In this study, genome-wide methylation profiling was carried out on paired healthy skin and wart samples in order to investigate the effects that benign HPV infection has on gene methylation status. To overcome this gap in knowledge, paired wart (n = 12) and normal skin (n = 12) samples were obtained from Arab males in order to perform DNA extraction and subsequent genome-wide methylation profiling on the Infinium Methylation EPIC Bead Chip microarray. Analysis of differential methylation revealed a clear pattern of discrimination between the wart and normal skin samples. In warts, the most differentially methylated (DM) genes included long non-coding RNAs (AC005884, AL049646.2, AC126121.2, AP001790.1, and AC107959.3), microRNAs (MIR374B, MIR596, MIR1255B1, MIR26B, and MIR196A2),snoRNAs (SNORD114-22, SNORD70, and SNORD114-31), pseudogenes (AC069366.1, RNU4ATAC11P, AC120057.1, NANOGP3, AC106038.2, TPT1P2, SDC4P, PKMP3, and VN2R3P), and protein-coding genes (AREG, GJB2, C12orf71, AC020909.2, S100A8, ZBED2, FABP7, and CYSLTR1). In addition, pathway analysis revealed that, among the most differentially methylated genes, STAT5A, RARA, MEF2D, MAP3K8, and THRA were the common regulators. It can be observed that HPV-induced warts involve a clear and unique epigenetic alteration to the host genome.

Pseudogenes as Biomarkers and Therapeutic Targets in Human Cancers

Pseudogenes are commonly labeled as "junk DNA" given their perceived nonfunctional status. However, the advent of large-scale genomics projects prompted a revisit of pseudogene biology, highlighting their key functional and regulatory roles in numerous diseases, including cancers. Integrative analyses of cancer data have shown that pseudogenes can be transcribed and even translated, and that pseudogenic DNA, RNA, and proteins can interfere with the activity and function of key protein coding genes, acting as regulators of oncogenes and tumor suppressors. Capitalizing on the available clinical research, we are able to get an insight into the spread and variety of pseudogene biomarker and therapeutic potential. In this chapter, we describe pseudogenes that fulfill their role as diagnostic or prognostic biomarkers, both as unique elements and in collaboration with other genes or pseudogenes. We also report that the majority of prognostic pseudogenes are overexpressed and exert an oncogenic role in colorectal, liver, lung, and gastric cancers. Finally, we highlight a number of pseudogenes that can establish future therapeutic avenues.

The HMGA1-pseudogene7 shows oncogenic activity in vivo

We have recently reported that transgenic mice overexpressing the HMGA1-pseudogene7 develop hematological neoplasia marked by monoclonal B-cell populations, and diagnosed as Diffuse Large B-cell Lymphoma. These findings prove the HMGA1-pseudogene7 oncogenic role in vivo.

The development and controversy of competitive endogenous RNA hypothesis in non-coding genes

As a momentous post-transcriptional regulator, microRNAs (miRNAs) are attracting more and more attention. The classical miRNAs regulated mechanism shows it binds to the targets' 3'UTR thus play the role in post-transcription. Meanwhile, single miRNA can target multiple genes, so those should compete to bind that miRNA. Vice versa, single gene can sponge mass of miRNAs as well. Thus the competitive endogenous RNAs (ceRNAs) hypothesis was put forward in 2011. The ceRNA hypothesis has made huge achievements, in particular in non-coding genes, which including long non-coding RNAs (lncRNAs), circle RNAs (circRNAs) and pseudogenes, even viral transcripts. It also contributed greatly to epigenetics development. However, an increasing number of controversies have occurred with applause. Based on this situation, this review introduces something in detail about the ceRNAs hypothesis achieved in lncRNAs, circRNAs, pseudogenes and viral transcripts, respectively. Meanwhile, it also covers controversy of the ceRNAs hypothesis.

The Interplay between Long Noncoding RNAs and Proteins of the Epigenetic Machinery in Ovarian Cancer

Comprehensive large-scale sequencing and bioinformatics analyses have uncovered a myriad of cancer-associated long noncoding RNAs (lncRNAs). Aberrant expression of lncRNAs is associated with epigenetic reprogramming during tumor development and progression, mainly due to their ability to interact with DNA, RNA, or proteins to regulate gene expression. LncRNAs participate in the control of gene expression patterns during development and cell differentiation and can be cell and cancer type specific. In this review, we described the potential of lncRNAs for clinical applications in ovarian cancer (OC). OC is a complex and heterogeneous disease characterized by relapse, chemoresistance, and high mortality rates. Despite advances in diagnosis and treatment, no significant improvements in long-term survival were observed in OC patients. A set of lncRNAs was associated with survival and response to therapy in this malignancy. We manually curated databases and used bioinformatics tools to identify lncRNAs implicated in the epigenetic regulation, along with examples of direct interactions between the lncRNAs and proteins of the epigenetic machinery in OC. The resources and mechanisms presented herein can improve the understanding of OC biology and provide the basis for further investigations regarding the selection of novel biomarkers and therapeutic targets.

Comprehensive analysis of LDHAP5 pseudogene expression and potential pathogenesis in ovarian serous cystadenocarcinoma

Background

We aimed to identify differentially expressed pseudogenes and explore their potential functions in four types of common gynecological malignancies (e.g., cervical squamous cell carcinoma, ovarian serous cystadenocarcinoma, uterine corpus endometrial carcinoma, and uterine carcinosarcoma) using bioinformatics technology.

Materials and methods

We identified up-regulated and down-regulated pseudogenes and built a pseudogene-miRNA-mRNA regulatory network through public datasets to explore their potential functions in carcinogenesis and cancer prognosis.

Results

Among the 63 up-regulated pseudogenes identified, LDHAP5 demonstrated the greatest potential as a candidate pseudogene due to its significant association with poor overall survival in ovarian serous cystadenocarcinoma. KEGG pathway analysis revealed that LDHAP5 showed significant enrichment in MicroRNAs in cancer, Pathway in cancer and PI3K-AKT signaling pathway. Further analysis revealed that EGFR was the potential target mRNA of LDHAP5, which may play an important role in ovarian serous cystadenocarcinoma.

Conclusions

LDHAP5 was associated with the occurrence and prognosis of ovarian serous cystadenocarcinoma, and thus shows potential as a novel therapeutic target against such cancer.

Genomic Indexing by Somatic Gene Recombination of mRNA/ncRNA - Does It Play a Role in Genomic Mosaicism, Memory Formation, and Alzheimer's Disease?

Recent evidence indicates that genomic individuality of neurons, characterized by DNA-content variation, is a common if not universal phenomenon in the human brain that occurs naturally but can also show aberrancies that have been linked to the pathomechanism of Alzheimer’s disease and related neurodegenerative disorders. Etiologically, this genomic mosaic has been suggested to arise from defects of cell cycle regulation that may occur either during brain development or in the mature brain after terminal differentiation of neurons. Here, we aim to draw attention towards another mechanism that can give rise to genomic individuality of neurons, with far-reaching consequences. This mechanism has its origin in the transcriptome rather than in replication defects of the genome, i.e., somatic gene recombination of RNA. We continue to develop the concept that somatic gene recombination of RNA provides a physiological process that, through integration of intronless mRNA/ncRNA into the genome, allows a particular functional state at the level of the individual neuron to be indexed. By insertion of defined RNAs in a somatic recombination process, the presence of specific mRNA transcripts within a definite temporal context can be “frozen” and can serve as an index that can be recalled at any later point in time. This allows information related to a specific neuronal state of differentiation and/or activity relevant to a memory trace to be fixed. We suggest that this process is used throughout the lifetime of each neuron and might have both advantageous and deleterious consequences.

HMGA1-pseudogene7 transgenic mice develop B cell lymphomas

We have recently identified and characterized two pseudogenes (HMGA1P6 and HMGA1P7) of the HMGA1 gene, which has a critical role in malignant cell transformation and cancer progression. HMGA1P6 and HMGAP17 act as microRNA decoy for HMGA1 and other cancer-related genes upregulating their protein levels. We have previously shown that they are upregulated in several human carcinomas, and their expression positively correlates with a poor prognosis and an advanced cancer stage. To evaluate in vivo oncogenic activity of HMGA1 pseudogenes, we have generated a HMGA1P7 transgenic mouse line overexpressing this pseudogene. By a mean age of 12 months, about 50% of the transgenic mice developed splenomegaly and accumulation of lymphoid cells in several body compartments. For these mice FACS and immunohistochemical analyses suggested the diagnosis of B-cell lymphoma that was further supported by clonality analyses and RNA expression profile of the pathological tissues of the HMGA1P7 transgenic tissues. Therefore, these results clearly demonstrate the oncogenic activity of HMGA1 pseudogenes in vivo.

MYC-regulated pseudogene HMGA1P6 promotes ovarian cancer malignancy via augmenting the oncogenic HMGA1/2

Pseudogenes have long been considered as nonfunctional genomic sequences. Recent studies have shown that they can potentially regulate the expression of protein-coding genes and are dysregulated in diseases including cancer. However, the potential roles of pseudogenes in ovarian cancer have not been well studied. Here we characterized the pseudogene expression profile in HGSOC (high-grade serous ovarian carcinoma) by microarray. We identified 577 dysregulated pseudogenes and most of them were up-regulated (538 of 577). HMGA1P6 (High mobility group AT-hook 1 pseudogene 6) was one of the overexpressed pseudogenes and its expression was inversely correlated with patient survival. Mechanistically, HMGA1P6 promoted ovarian cancer cell malignancy by acting as a ceRNA (competitive endogenous RNA) that led to enhanced HMGA1 and HMGA2 expression. Importantly, HMGA1P6 was transcriptionally activated by oncogene MYC in ovarian cancer. Our findings reveal that MYC may contribute to oncogenesis through transcriptional regulation of pseudogene HMGA1P6 in ovarian cancer.

Pseudogene-Derived lncRNAs and Their miRNA Sponging Mechanism in Human Cancer

Pseudogenes, abundant in the human genome, are traditionally considered as non-functional “junk genes.” However, recent studies have revealed that pseudogenes act as key regulators at DNA, RNA or protein level in diverse human disorders (including cancer), among which pseudogene-derived long non-coding RNA (lncRNA) transcripts are extensively investigated and has been reported to be frequently dysregulated in various types of human cancer. Growing evidence demonstrates that pseudogene-derived lncRNAs play important roles in cancer initiation and progression by serving as competing endogenous RNAs (ceRNAs) through competitively binding to shared microRNAs (miRNAs), thus affecting both their cognate genes and unrelated genes. Herein, we retrospect those current findings about expression, functions and potential ceRNA mechanisms of pseudogene-derived lncRNAs in human cancer, which may provide us with some crucial clues in developing potential targets for cancer therapy in the future.

Finding relationships among biological entities

Confusion over the concepts of “relationships” and “similarities” lies at the heart of many battles over the direction and intent of research projects. Here is a short story that demonstrates the difference between the two concepts: You look up at the clouds, and you begin to see the shape of a lion. The cloud has a tail, like a lion’s tale, and a fluffy head, like a lion’s mane. With a little imagination the mouth of the lion seems to roar down from the sky. You have succeeded in finding similarities between the cloud and a lion. If you look at a cloud and you imagine a tea kettle producing a head of steam and you recognize that the physical forces that create a cloud and the physical forces that produced steam from a heated kettle are the same, then you have found a relationship. Most popular classification algorithms operate by grouping together data objects that have similar properties or values. In so doing, they may miss finding the true relationships among objects. Traditionally, relationships among data objects are discovered by an intellectual process. In this chapter, we will discuss the scientific gains that come when we classify biological entities by relationships, not by their similarities.

Re-recognition of pseudogenes: From molecular to clinical applications

Pseudogenes were initially regarded as "nonfunctional" genomic elements that did not have protein-coding abilities due to several endogenous inactivating mutations. Although pseudogenes are widely expressed in prokaryotes and eukaryotes, for decades, they have been largely ignored and classified as gene "junk" or "relics". With the widespread availability of high-throughput sequencing analysis, especially omics technologies, knowledge concerning pseudogenes has substantially increased. Pseudogenes are evolutionarily conserved and derive primarily from a mutation or retrotransposon, conferring the pseudogene with a "gene repository" role to store and expand genetic information. In contrast to previous notions, pseudogenes have a variety of functions at the DNA, RNA and protein levels for broadly participating in gene regulation to influence the development and progression of certain diseases, especially cancer. Indeed, some pseudogenes have been proven to encode proteins, strongly contradicting their "trash" identification, and have been confirmed to have tissue-specific and disease subtype-specific expression, indicating their own value in disease diagnosis. Moreover, pseudogenes have been correlated with the life expectancy of patients and exhibit great potential for future use in disease treatment, suggesting that they are promising biomarkers and therapeutic targets for clinical applications. In this review, we summarize the natural properties, functions, disease involvement and clinical value of pseudogenes. Although our knowledge of pseudogenes remains nascent, this field deserves more attention and deeper exploration.

Long noncoding RNA ARHGAP27P1 inhibits gastric cancer cell proliferation and cell cycle progression through epigenetically regulating p15 and p16

Long noncoding RNAs (lncRNAs) have emerged as important regulators in the development and progression of gastric cancer (GC). ARHGAP27P1 is a pseudogene-derived lncRNA, and it has been found to be associated with GC in our preliminary study, but this association has not been studied further. Herein, we confirmed that ARHGAP27P1 was significantly downregulated in GC tissues, plasma and cells. Low expression of ARHGAP27P1 was closely associated with advanced TNM stage, increased invasion depth and lymphatic metastasis. Low ARHGAP27P1 expression also predicted a poor prognosis in GC patients. Functionally, overexpression of ARHGAP27P1 inhibited proliferation, invasion, and migration in GC cells, while silencing of ARHGAP27P1 showed the opposite effects. Mechanistic investigations showed that ARHGAP27P1 had a key role in G0/G1 arrest. We further demonstrated that ARHGAP27P1 was associated with Jumonji-domain containing 3 (JMJD3) and that this association was required for the demethylation of H3K27me3, thereby epigenetically activating expression of p15, p16 and p57. Moreover, knockdown of JMJD3, p15, or p16 consistently reversed the inhibitory effects of ARHGAP27P1 in cell proliferation and cell cycle progression. Taken together, these results suggest that lncRNA ARHGAP27P1, as a novel cell cycle regulator, may serve as a potential target for GC prevention and treatment in human GC.

A disparate role of RP11-424C20.2/UHRF1 axis through control of tumor immune escape in liver hepatocellular carcinoma and thymoma

The immune system is critical in modulating cancer progression. Pseudogenes are a special type of long non-coding RNAs that regulate different tumorigenic processes. However, the potential roles of pseudogenes in tumor-immune interaction remain largely unclear. Here, we reported that pseudogene RP11-424C20.2 and its parental gene UHRF1 were frequently up-regulated and positively correlated in liver hepatocellular carcinoma (LIHC) and thymoma (THYM), but associated with distinct clinical outcomes. We further found that RP11-424C20.2 may act as a competing endogenous RNA (ceRNA) to increase UHRF1 expression through sponging miR-378a-3p. Functional enrichment analysis showed a strong association of UHRF1 with immune-related biological processes. We also observed that UHRF1 expression significantly correlated with immune infiltration, and different types of tumor-infiltrating immune cells displayed different impacts on clinical outcomes. Furthermore, UHRF1 expression in LIHC and THYM showed an opposite correlation with biomarkers from monocyte, dendritic cell, Th1 and T cell exhaustion. Mechanism investigations revealed that RP11-424C20.2/UHRF1 axis regulated immune escape of LIHC and THYM at least partly through IFN-γ-mediated CLTA-4 and PD-L1 pathway. These findings demonstrate a disparate role of RP11-424C20.2/UHRF1 axis in LIHC and THYM via regulating immune infiltrates, and also indicate a therapeutic value for UHRF1 inhibitors in combination with anti-PD-L1/CLTA-4 blockade.

Pseudogene Associated Recurrent Gene Fusion in Prostate Cancer

We present the functional characterization of a pseudogene associated recurrent gene fusion in prostate cancer. The fusion gene KLK4-KLKP1 is formed by the fusion of the protein coding gene KLK4 with the noncoding pseudogene KLKP1. Screening of a cohort of 659 patients (380 Caucasian American; 250 African American, and 29 patients from other races) revealed that the KLK4-KLKP1 is expressed in about 32% of prostate cancer patients. Correlative analysis with other ETS gene fusions and SPINK1 revealed a concomitant expression pattern of KLK4-KLKP1 with ERG and a mutually exclusive expression pattern with SPINK1, ETV1, ETV4, and ETV5. Development of an antibody specific to KLK4-KLKP1 fusion protein confirmed the expression of the full-length KLK4-KLKP1 protein in prostate tissues. The in vitro and in vivo functional assays to study the oncogenic properties of KLK4-KLKP1 confirmed its role in cell proliferation, cell invasion, intravasation, and tumor formation. Presence of strong ERG and AR binding sites located at the fusion junction in KLK4-KLKP1 suggests that the fusion gene is regulated by ERG and AR. Correlative analysis of clinical data showed an association of KLK4-KLKP1 with lower preoperative PSA values and in young men (<50 years) with prostate cancer. Screening of patient urine samples showed that KLK4-KLKP1 can be detected noninvasively in urine. Taken together, we present KLK4-KLKP1 as a class of pseudogene associated fusion transcript in cancer with potential applications as a biomarker for routine screening of prostate cancer.

Deregulation of H19 is associated with cervical carcinoma

CACX is one of the most common cancer affecting women world-wide. Here, expression microarray analysis revealed 8 over-expressed transcribed pseudogenes (GBP1P1, HLA-DRB6, HLA-H, SLC6A10P, NAPSB, KRT16P2, PTTG3P and RNF126P1), down-regulated 7 lincRNAs (H19, MIR100HG, MEG3, DIO3OS, HOXA11-AS, CD27-AS1 and EPB41L4A-AS) and 6 snoRNAs (SNORD97, SNORD3A, SNORD3C, SNORD3D, SNORA12 and SCARNA9) as DEncGs (log2 fold-change ≥ ±1.0) in CACX. Consequently, down-regulation of lincRNA MEG3 and over-expression of pseudogenes, GBP1P1 and PTTG3P in the microarray analysis were found concordant with the real-time quantitative PCR results upon validation. Then, Ingenuity® Pathway analysis (IPA®) analysis with deregulated DEncGs identified functionally important gene, H19. Further, validation (n = 52) of expression confirmed frequent downregulation of H19 with significant association with its deletion (LOH) and promoter methylation (n = 128) in CACX. Moreover, clinicopathological analysis found Indian CACX patients (n = 26) with alterations of H19 by deletion or, promoter methylation with concomitant low expression have poor prognosis.

Genetic variants of PARP4 gene and PARP4P2 pseudogene in patients with multiple primary tumors including thyroid cancer

Recently, the PARP4 gene has been identified as a possible susceptibility gene of primary thyroid and breast cancers. We analyzed PARP4 in 53 patients with multiple primary cancers including a thyroid cancer (TC), in 74 patients with TC alone, and in 88 healthy donors. Two PARP4 intronic variants within the IVS29 (c.3543 + 44T > C) and the IVS22 (c.2758 + 9G > A) were found only in the two patient groups. Moreover, we found a rare variant (r.522C > A) within a PARP4 pseudogene (PARP4P2) in one patient with four primary tumors, and with a familial cancer history. PARP4 mRNA was absent in all primary tumors and matched normal tissues, whereas the pseudogene variant transcript was always expressed. Consistently, immunostaining for PARP4 protein was negative at nuclear level in all tissues, thus suggesting that PARP4P2 pseudogene variant could alter its regulatory role on PARP4, inducing the down-regulation of PARP4 expression at both tumor and normal tissues level. In conclusion, germline intronic PARP4 variants could be a risk factor for the development of TC, and PARP4P2 pseudogene variations associated with PARP4 down-regulation may confer susceptibility to develop multiple metachronous cancers.

The long noncoding RNA TPTE2P1 promotes the viability of colorectal cancer cells

Long noncoding RNAs (lncRNAs) have important functions in tumor development and progression, including colorectal cancer (CRC), but their roles are not completely understood. In this study, the roles of the lncRNA transmembrane phosphoinositide 3-phosphatase and tensin homolog 2 pseudogene 1 (TPTE2P1), previously implicated in gallbladder cancer cell migration and invasion, were evaluated in CRC. In particular, quantitative polymerase chain reaction was used to quantify TPTE2P1 levels in tumor tissues and cell lines. The association between TPTE2P1 and survival was analyzed using the online tool OncoLnc. The 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide assays, colony formation assays, and flow cytometry were used to evaluate the effects of TPTE2P1 on viability, cell cycle progression, and apoptosis. Signaling pathway proteins were quantitated by Western blot analysis. Finally, the role of TPTE2P1 was analyzed in vivo using mouse models. TPTE2P1 levels were significantly higher in CRC tissues than in adjacent normal tissues. Higher expression was associated with a poor survival rate. The silencing of TPTE2P1 led to cell cycle arrest at the S phase and thereby inhibited cell viability. TPTE2P1 knockdown also caused cancer cell apoptosis via the activation of the apoptosis regulator (BCL2)/caspase 3 signaling cascade. In addition, the inhibition of TPTE2P1 had suppressive effects on tumors in vivo. TPTE2P1 is upregulated in CRC and plays essential roles in the regulation of cell viability in vitro and tumor formation in vivo.

High Expression of Pseudogene PTTG3P Indicates a Poor Prognosis in Human Breast Cancer

Pseudogenes play pivotal roles in tumorigenesis. Previous studies have suggested that pituitary tumor-transforming 3, pseudogene (PTTG3P), serves as an oncogene in human cancers. However, its expression pattern, biological function, and underlying mechanism in breast cancer remain unknown. In this study, we demonstrated an elevated expression of PTTG3P in breast cancer and discovered that PTTG3P expression correlated negatively with estrogen receptor (ER) and progesterone receptor (PR) status, but linked positively to basal-like status, triple-negative breast cancer status, Nottingham prognostic index (NPI), and Scarff-Bloom-Richardson grade. High expression of PTTG3P was also found to be associated with a poor prognosis of breast cancer. To explore the potential mechanisms of PTTG3P, a PTTG3P-microRNA (miRNA)-mRNA regulatory network was established. Co-expressed genes of PTTG3P were also obtained. Enrichment analysis for these co-expressed genes revealed that they were significantly enriched in mitotic nuclear division and cell cycle. Subsequent research on mechanism of PTTG3P indicated that its expression correlated positively with PTTG1 expression. However, no significant expression correlation between PTTG3P and PTTG2 was observed. Taken together, our findings suggest that increased expression of pseudogene PTTG3P may be used as a promising prognostic biomarker and novel therapeutic target for breast cancer.