Expressed pseudogenes in the transcriptional landscape of human cancers

A prognostic binary classifier comprised of five critical mRNAs stratified pancreatic cancer patients following resection

Background

Pancreatic cancer is characterized by an extremely poor prognosis, even following potentially curative resection. Classical prognostic markers such as histopathological or clinical parameters have limited predictive power. The present study aimed to establish a prognostic model combining mRNA expression data with histopathological and clinical data to better predict survival and stratify pancreatic cancer patients following resection. We pioneered three models in one study and systematically evaluated the clinical benefits of all three models.

Methods

To identify differentially expressed genes in pancreatic cancer, mRNA data from normal (GTEx database) and pancreatic cancer (TCGA database) tissues were used. Survival analysis was carried out to identify prognosis-relevant genes from the identified differentially expressed genes and LASSO regression was used to filter out hub genes. The risk score of several hub genes was calculated according to gene expression and coefficients. Validation was carried out using an independent set of GEO microarray data. Multivariate COX regression was used for identifying independent clinical and pathological risk factors related to patient's survival in the TCGA database and a prognostic model combining mRNA expression data with histopathological and clinical data was established. Another prognostic model using clinicopathological factors from the SEER database was conceived based on multivariate COX regression. NRI (net reclassification improvement) and IDI (integrated discrimination index) were used to compare the predictive capabilities of the different models.

Results

We identified 1589 differentially expressed genes (DEGs) through the comparison of normal and pancreatic cancer tissues, of whom 317 were associated with prognosis(p < 0.05). LASSO regression identified five hub genes, MYEOV, ANXA2P2, MET, CEP55, and KRT7, that were used for the five-mRNA-classifier prognostic model. The classifier could stratify patients into a short and long survival group: 5-year overall survival in the training set (TCGA, 6 % vs 52 %, p < 0.001), test set (TCGA, 18 % vs 55 %,p < 0.01) and external validation set (GEO, 0 % vs 25 %, p < 0.05). Sensitivity analysis showed that the mRNA model (model 1) was better than the clinicopathological no-mRNA model (model 2) in predicting 5-year survival in the TCGA database (AUC: 0.877 vs 0.718, z = 3.165, p < 0.01) and better than the multi-factor prognostic model (model 3) from the SEER database (AUC: 0.754, z = 2.637, p < 0.01). On predictive performance, model 1 improved model 2 (NRI = 0.084, z = 1.288, p = 0.198; IDI = 0.055, z = 1.041,p = 0.298) and model 3 (NRI = 0.167,z = 1.961,p = 0.05; IDI = 0.086, z = 1.427, p = 0.154).

Conclusion

The five-mRNA-classifier is a reliable and feasible instrument to predict the prognosis of pancreatic cancer patients following resection. It might help in patiens counseling and assist clinicians in providing individualized treatment for patients in different risk groups.

Loss to gain: pseudogenes in microorganisms, focusing on eubacteria, and their biological significance

Pseudogenes are defined as "non-functional" copies of corresponding parent genes. The cognition of pseudogenes continues to be refreshed through accumulating and updating research findings. Previous studies have predominantly focused on mammals, but pseudogenes have received relatively less attention in the field of microbiology. Given the increasing recognition on the importance of pseudogenes, in this review, we focus on several aspects of microorganism pseudogenes, including their classification and characteristics, their generation and fate, their identification, their abundance and distribution, their impact on virulence, their ability to recombine with functional genes, the extent to which some pseudogenes are transcribed and translated, and the relationship between pseudogenes and viruses. By summarizing and organizing the latest research progress, this review will provide a comprehensive perspective and improved understanding on pseudogenes in microorganisms. KEY POINTS: • Concept, classification and characteristics, identification and databases, content, and distribution of microbial pseudogenes are presented. • How pseudogenization contribute to pathogen virulence is highlighted. • Pseudogenes with potential functions in microorganisms are discussed.

Structural Features and Physiological Associations of Human 14-3-3ζ Pseudogenes

There are about 14,000 pseudogenes that are mutated or truncated sequences resembling functional parent genes. About two-thirds of pseudogenes are processed, while others are duplicated. Although initially thought dead, emerging studies indicate they have functional and regulatory roles. We study 14-3-3ζ, an adaptor protein that regulates cytokine signaling and inflammatory diseases, including rheumatoid arthritis, cancer, and neurological disorders. To understand how 14-3-3ζ (gene symbol YWHAZ) performs diverse functions, we examined the human genome and identified nine YWHAZ pseudogenes spread across many chromosomes. Unlike the 32 kb exon-to-exon sequence in YWHAZ, all pseudogenes are much shorter and lack introns. Out of six, four YWHAZ exons are highly conserved, but the untranslated region (UTR) shows significant diversity. The putative amino acid sequence of pseudogenes is 78-97% homologous, resulting in striking structural similarities with the parent protein. The OMIM and Decipher database searches revealed chromosomal loci containing pseudogenes are associated with human diseases that overlap with the parent gene. To the best of our knowledge, this is the first report on pseudogenes of the 14-3-3 family protein and their implications for human health. This bioinformatics-based study introduces a new insight into the complexity of 14-3-3ζ's functions in biology.

The network structure of hematopoietic cancers

Hematopoietic cancers (HCs) are a heterogeneous group of malignancies that affect blood, bone marrow and lymphatic system. Here, by analyzing 1960 RNA-Seq samples from three independent datasets, we explored the co-expression landscape in HCs, by inferring gene co-expression networks (GCNs) with four cancer phenotypes (B and T-cell acute leukemia -BALL, TALL-, acute myeloid leukemia -AML-, and multiple myeloma -MM-) as well as non-cancer bone marrow. We characterized their structure (topological features) and function (enrichment analyses). We found that, as in other types of cancer, the highest co-expression interactions are intra-chromosomal, which is not the case for control GCNs. We also detected a highly co-expressed group of overexpressed pseudogenes in HC networks. The four GCNs present only a small fraction of common interactions, related to canonical functions, like immune response or erythrocyte differentiation. With this approach, we were able to reveal cancer-specific features useful for detection of disease manifestations.

Hsa-miR-665 Is a Promising Biomarker in Cancer Prognosis

Biomarkers are biomolecules used to identify or predict the presence of a specific disease or condition. They play an important role in early diagnosis and may be crucial for treatment. MicroRNAs (miRNAs), a group of small non-coding RNAs, are more and more regarded as promising biomarkers for several reasons. Dysregulation of miRNAs has been linked with development of several diseases, including many different types of cancer, and abnormal levels can be present in early stages of tumor development. Because miRNAs are stable molecules secreted and freely circulating in blood and urine, they can be sampled with little or no invasion. Here, we present an overview of the current literature, focusing on the types of cancers for which dysregulation of miR-665 has been associated with disease progression, recurrence, and/or prognosis. It needs to be emphasized that the role of miR-665 sometimes seems ambiguous, in the sense that it can be upregulated in one cancer type and downregulated in another and can even change during the progression of the same cancer. Caution is thus needed before using miR-665 as a biomarker, and extrapolation between different cancer types is not advisable. Moreover, more detailed understanding of the different roles of miR-665 will help in determining its potential as a diagnostic and prognostic biomarker.

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.

Small Open Reading Frame-Encoded Micro-Peptides: An Emerging Protein World

Small open reading frames (sORFs) are often overlooked features in genomes. In the past, they were labeled as noncoding or "transcriptional noise". However, accumulating evidence from recent years suggests that sORFs may be transcribed and translated to produce sORF-encoded polypeptides (SEPs) with less than 100 amino acids. The vigorous development of computational algorithms, ribosome profiling, and peptidome has facilitated the prediction and identification of many new SEPs. These SEPs were revealed to be involved in a wide range of basic biological processes, such as gene expression regulation, embryonic development, cellular metabolism, inflammation, and even carcinogenesis. To effectively understand the potential biological functions of SEPs, we discuss the history and development of the newly emerging research on sORFs and SEPs. In particular, we review a range of recently discovered bioinformatics tools for identifying, predicting, and validating SEPs as well as a variety of biochemical experiments for characterizing SEP functions. Lastly, this review underlines the challenges and future directions in identifying and validating sORFs and their encoded micropeptides, providing a significant reference for upcoming research on sORF-encoded peptides.

Human VDAC pseudogenes: an emerging role for VDAC1P8 pseudogene in acute myeloid leukemia

BACKGROUND

Voltage-dependent anion selective channels (VDACs) are the most abundant mitochondrial outer membrane proteins, encoded in mammals by three genes, VDAC1, 2 and 3, mostly ubiquitously expressed. As 'mitochondrial gatekeepers', VDACs control organelle and cell metabolism and are involved in many diseases. Despite the presence of numerous VDAC pseudogenes in the human genome, their significance and possible role in VDAC protein expression has not yet been considered.

RESULTS

We investigated the relevance of processed pseudogenes of human VDAC genes, both in physiological and in pathological contexts. Using high-throughput tools and querying many genomic and transcriptomic databases, we show that some VDAC pseudogenes are transcribed in specific tissues and pathological contexts. The obtained experimental data confirm an association of the VDAC1P8 pseudogene with acute myeloid leukemia (AML).

CONCLUSIONS

Our in-silico comparative analysis between the VDAC1 gene and its VDAC1P8 pseudogene, together with experimental data produced in AML cellular models, indicate a specific over-expression of the VDAC1P8 pseudogene in AML, correlated with a downregulation of the parental VDAC1 gene.

PA2G4P4 Promotes Glioma Cell Migration and Tumorigenesis through the PTEN/AKT/mTOR Signaling Pathway

Dysregulation of pseudogene expression is closely related to the progression of various cancers, including glioma. Proliferation-associated 2G4 pseudogene 4 (PA2G4P4) could affect cell viability and apoptosis of glioma cells. However, the specific regulatory mechanism of PA2G4P4 is not clear. In this paper, we found that PA2G4P4 overexpres-sion promoted glioma cell proliferation, migration and cell cycle progression, whereas PA2G4P4 knockdown inhibited cancer progression. Knockdown of PA2G4P4 also suppressed the tumorigenesis of glioma cells in vivo. Furthermore, knockdown of PA2G4 after overexpression of PA2G4P4 decreased the cell viability and migration ability to normal level. The protein level of a tumor suppressor gene phosphatase and tensing homolog (PTEN) was greatly decreased in U87 cells after PA2G4P4 overexpression, while increased after PA2G4 knockdown; on the contrary, the protein levels of P-AKT and P-S6 were obviously induced in U87 cells after PA2G4P4 overexpression, and decreased after PA2G4 knockdown. The cell ability, colony formation ability and cell migration ability were all recovered to normal level by adding an AKT inhibitor MK2206 to the glioma cells, which were induced by PA2G4P4 overexpression. Our results revealed that PA2G4P4 could regulate glioma cell proliferation and migration through PTEN/AKT/mTOR signaling pathway by targeting PA2G4 gene. PA2G4P4 may become a target for glioma treatment.

Improving the genome and proteome annotations of the marine model diatom Thalassiosira pseudonana using a proteogenomics strategy

Diatoms are unicellular eukaryotic phytoplankton that account for approximately 20% of global carbon fixation and 40% of marine primary productivity; thus, they are essential for global carbon biogeochemical cycling and climate. The availability of ten diatom genome sequences has facilitated evolutionary, biological and ecological research over the past decade; however, a complimentary map of the diatom proteome with direct measurements of proteins and peptides is still lacking. Here, we present a proteome map of the model marine diatom Thalassiosira pseudonana using high-resolution mass spectrometry combined with a proteogenomic strategy. In-depth proteomic profiling of three different growth phases and three nutrient-deficient samples identified 9526 proteins, accounting for ~ 81% of the predicted protein-coding genes. Proteogenomic analysis identified 1235 novel genes, 975 revised genes, 104 splice variants and 234 single amino acid variants. Furthermore, our quantitative proteomic analysis experimentally demonstrated that a considerable number of novel genes were differentially translated under different nutrient conditions. These findings substantially improve the genome annotation of T. pseudonana and provide insights into new biological functions of diatoms. This relatively comprehensive diatom proteome catalog will complement available diatom genome and transcriptome data to advance biological and ecological research of marine diatoms.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-022-00161-y.

New insights into circRNA and its mechanisms in angiogenesis regulation in ischemic stroke: a biomarker and therapeutic target

Ischemic stroke accounts for about 71% of strokes worldwide. Due to limited recommended therapeutics for ischemic stroke, more attention is focused on angiogenesis in ischemic stroke. Not long after ischemic stroke, angiogenesis arises and is vital for the prognosis. Various pro-angiogenic, anti-angiogenic factors and their downstream pathways engage in angiogenesis regulation. CircRNAs are differentially expressed after ischemic stroke. Up to now, circRNAs have been found to exert many functions in regulating apoptosis, autophagy, proliferation, and differentiation of neurons and neural stem cells mainly as miRNAs sponges or proteins decoy. Thus, many circRNAs are considered promising biomarkers or therapeutic targets for ischemic stroke. Besides, circRNAs participate in the modulation of endothelial-mesenchymal transition and blood-brain barrier maintenance. Moreover, circRNAs play significant roles in endothelial dysfunction concerning inflammation responses, apoptosis, proliferation, and migration. They correlate with many angiogenesis-related signaling pathways and genes via the circRNA/miRNA/mRNA network. Novel insights into circRNAs significance in angiogenesis regulation in ischemic stroke could be provided for further researches on the clinical application of circRNAs in ischemic stroke.

Evolution and function of developmentally dynamic pseudogenes in mammals

BACKGROUND

Pseudogenes are excellent markers for genome evolution, which are emerging as crucial regulators of development and disease, especially cancer. However, systematic functional characterization and evolution of pseudogenes remain largely unexplored.

RESULTS

To systematically characterize pseudogenes, we date the origin of human and mouse pseudogenes across vertebrates and observe a burst of pseudogene gain in these two lineages. Based on a hybrid sequencing dataset combining full-length PacBio sequencing, sample-matched Illumina sequencing, and public time-course transcriptome data, we observe that abundant mammalian pseudogenes could be transcribed, which contribute to the establishment of organ identity. Our analyses reveal that developmentally dynamic pseudogenes are evolutionarily conserved and show an increasing weight during development. Besides, they are involved in complex transcriptional and post-transcriptional modulation, exhibiting the signatures of functional enrichment. Coding potential evaluation suggests that 19% of human pseudogenes could be translated, thus serving as a new way for protein innovation. Moreover, pseudogenes carry disease-associated SNPs and conduce to cancer transcriptome perturbation.

CONCLUSIONS

Our discovery reveals an unexpectedly high abundance of mammalian pseudogenes that can be transcribed and translated, and these pseudogenes represent a novel regulatory layer. Our study also prioritizes developmentally dynamic pseudogenes with signatures of functional enrichment and provides a hybrid sequencing dataset for further unraveling their biological mechanisms in organ development and carcinogenesis in the future.

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.

Loss of KLK4::KLKP1 pseudogene expression by RNA chromogenic in-situ hybridization is associated with PTEN loss and increased risk of biochemical recurrence in a cohort of middle eastern men with prostate cancer

BACKGROUND

KLK4::KLKP1 fusion is a recently described pseudogene that is enriched in prostate cancer (PCa). This new biomarker has not been characterized in the Middle Eastern population.

OBJECTIVE

To establish the incidence and prognostic value of KLK4::KLKP1 fusion in a cohort of Middle Eastern men with PCa and explore the relationship of this marker to other relevant biomarkers (PTEN, ERG, SPINK1).

DESIGN, SETTING, AND PARTICIPANTS

We interrogated a cohort of 340 Middle Eastern men with localized PCa treated by radical prostatectomy between 2005 and 2015. KLK4::KLKP1 fusion status was assessed by RNA Chromogenic in situ hybridization (CISH) and correlated to pathological and clinical parameters.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

RNA-CISH expression of KLK4::KLKP1 was correlated with prognostic factors, ERG, PTEN, and SPINK1 expression, and biochemical recurrence (BCR) following prostatectomy.

RESULTS AND LIMITATIONS

51.7% of patient samples showed positive KLK4::KLKP1 expression; more commonly in cores of PCa (38%) versus non-cancer (20.6%) (p < 0.0001) and in lower Gleason Grade Group tumors (1-3) vs (4-5). KLK4::KLKP1 expression positively correlated with ERG positivity and inversely associated with PTEN loss. No significant association was found with SPINK1 expression, seminal vesicle invasion, positive surgical margin, pathological stage, or patient age (< 50 or ≥ 50). The association between PTEN loss and BCR increased when combined with KLK4::KLKP1 negativity (HR 2.31, CI 1.03-5.20, p = 0.042).

CONCLUSIONS

KLK4::KLKP1 expression is more common in this cohort of Middle Eastern men than has been reported in North American men. It is associated with ERG positivity and inversely correlated with PTEN loss. In isolation, KLK4::KLKP1 expression was not significantly associated with clinical outcome or pathological parameters. However, its expression is associated with certain molecular subtypes (ERG-positive, PTEN-intact) and as we demonstrate may help further stratify the risk of recurrence within these groups.

TOMM40 genetic variants associated with healthy aging and longevity: a systematic review

Introduction

Healthy aging relies on mitochondrial functioning because this organelle provides energy and diminishes oxidative stress. Single nucleotide polymorphisms (SNPs) in TOMM40, a critical gene that produces the outer membrane protein TOM40 of mitochondria, have been associated with mitochondrial dysfunction and neurodegenerative processes. Yet it is not clear whether or how the mitochondria may impact human longevity. We conducted this review to ascertain which SNPs have been associated with markers of healthy aging.

Methods

Using the PRISMA methodology, we conducted a systematic review on PubMed and Embase databases to identify associations between TOMM40 SNPs and measures of longevity and healthy aging.

Results

Twenty-four articles were selected. The TOMM40 SNPs rs2075650 and rs10524523 were the two most commonly identified and studied SNPs associated with longevity. The outcomes associated with the TOMM40 SNPs were changes in BMI, brain integrity, cognitive functions, altered inflammatory network, vulnerability to vascular risk factors, and longevity.

Discussions

Our systematic review identified multiple TOMM40 SNPs potentially associated with healthy aging. Additional research can help to understand mechanisms in aging, including resilience, prevention of disease, and adaptation to the environment.

Pseudogenes and Liquid Phase Separation in Epigenetic Expression

Pseudogenes have been considered as non-functional genes. However, peptides and long non-coding RNAs produced by pseudogenes are expressed in different tumors. Moreover, the dysregulation of pseudogenes is associated with cancer, and their expressions are higher in tumors compared to normal tissues. Recent studies show that pseudogenes can influence the liquid phase condensates formation. Liquid phase separation involves regulating different epigenetic stages, including transcription, chromatin organization, 3D DNA structure, splicing, and post-transcription modifications like m⁶A. Several membrane-less organelles, formed through the liquid phase separate, are also involved in the epigenetic regulation, and their defects are associated with cancer development. However, the association between pseudogenes and liquid phase separation remains unrevealed. The current study sought to investigate the relationship between pseudogenes and liquid phase separation in cancer development, as well as their therapeutic implications.

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.

The role of lncRNA-mediated ceRNA regulatory networks in pancreatic cancer

Non-coding RNAs (ncRNAs), which occupy the vast majority of human transcripts are known for their inability to encode proteins. NcRNAs consist of a diverse range of RNA species, including long non-coding RNAs (lncRNAs), which have significant meaning for epigenetic modification, post-transcriptional regulation of target genes, molecular interference, etc. The dysregulation of ncRNAs will mediate the pathogenesis of diverse human diseases, like cancer. Pancreatic cancer, as one of the most lethal malignancies in the digestive system that is hard to make a definite diagnosis at an early clinicopathological stage with a miserable prognosis. Therefore, the identification of potential and clinically applicable biomarker is momentous to improve the overall survival rate and positively ameliorate the prognosis of patients with pancreatic carcinoma. LncRNAs as one kind of ncRNAs exert multitudinous biological functions, and act as molecular sponges, relying on microRNA response elements (MREs) to competitively target microRNAs (miRNAs), thereby attenuating the degradation or inhibition of miRNAs to their own downstream protein-coding target genes, also thus regulating the initiation and progression of neoplasms. LncRNAs, which emerge aforementioned function are called competing endogenous RNAs (ceRNAs). Consequently, abundant research of lncRNAs as potential biomarkers is of critical significance for the molecular diagnosis, targeted therapy, as well as prognosis monitoring of pancreatic cancer.

Identification of Type 2 Diabetes Biomarkers From Mixed Single-Cell Sequencing Data With Feature Selection Methods

Diabetes is the most common disease and a major threat to human health. Type 2 diabetes (T2D) makes up about 90% of all cases. With the development of high-throughput sequencing technologies, more and more fundamental pathogenesis of T2D at genetic and transcriptomic levels has been revealed. The recent single-cell sequencing can further reveal the cellular heterogenicity of complex diseases in an unprecedented way. With the expectation on the molecular essence of T2D across multiple cell types, we investigated the expression profiling of more than 1,600 single cells (949 cells from T2D patients and 651 cells from normal controls) and identified the differential expression profiling and characteristics at the transcriptomics level that can distinguish such two groups of cells at the single-cell level. The expression profile was analyzed by several machine learning algorithms, including Monte Carlo feature selection, support vector machine, and repeated incremental pruning to produce error reduction (RIPPER). On one hand, some T2D-associated genes (MTND4P24, MTND2P28, and LOC100128906) were discovered. On the other hand, we revealed novel potential pathogenic mechanisms in a rule manner. They are induced by newly recognized genes and neglected by traditional bulk sequencing techniques. Particularly, the newly identified T2D genes were shown to follow specific quantitative rules with diabetes prediction potentials, and such rules further indicated several potential functional crosstalks involved in T2D.

Role of hepatitis B virus in development of hepatocellular carcinoma: Focus on covalently closed circular DNA

Chronic infection with hepatitis B virus (HBV) remains a major global health problem, especially in developing countries. It may lead to prolonged liver damage, fibrosis, cirrhosis, and hepatocellular carcinoma. Persistent chronic HBV infection is related to host immune response and the stability of the covalently closed circular DNA (cccDNA) in human hepatocytes. In addition to being essential for viral transcription and replication, cccDNA is also suspected to play a role in persistent HBV infections or hepatitis relapses since cccDNA is very stable in non-dividing human hepatocytes. Understanding the pathogenicity and oncogenicity of HBV components would be essential in the development of new diagnostic tools and treatment strategies. This review summarizes the role and molecular mechanisms of HBV cccDNA in hepatocyte transformation and hepatocarcinogenesis and current efforts to its detection and targeting.

Non-Coding RNAs Are Brokers in Breast Cancer Interactome Networks and Add Discrimination Power between Subtypes

Despite the power of high-throughput genomics, most non-coding RNA (ncRNA) biotypes remain hard to identify, characterize, and validate. This is a clear indication that intensive next-generation sequencing research has led to great efficiency and accuracy in detecting ncRNAs, but not in their functionalization. Computational scientists continue to support the discovery process by spotting significant data features (expression or mutational profiles), elucidating phenotype uncertainty, and delineating complex regulation landscapes for biological pathways and pathophysiological processes. With reference to transcriptome regulation dynamics in cancer, this work introduces a novel network-driven inference approach designed to reveal the potential role of computationally identified ncRNAs in discriminating between breast cancer (BC) subtypes beyond the traditional gene expression signatures. As heterogeneity cast in the subtypes is a characteristic of most cancers, the proposed approach is generalizable beyond BC. Expression profiles of a wide transcriptome spectrum were obtained for a number of BC patients (and controls) listed in TCGA and processed with RNA-Seq. The well-known PAM50 subtype signature was available for the samples and used to move from differentially expressed transcript profiles to subtype-specific biclusters associating gene patterns with patients. Co-expressed gene networks were then generated and annotations were provided, focusing on the biclusters with basal and luminal signatures. These were used to build template maps, i.e., networks in which to embed the ncRNAs and contextually functionalize them based on their interactors. This inference approach is able to assess the influence of ncRNAs at the level of BC subtype. Network topology was considered through the brokerage measure to account for disruptiveness effects induced by the removal of nodes corresponding to ncRNAs. Equivalently, it is shown that ncRNAs can act as brokers of network interactome dynamics, and removing them allows the refinement of subtype-related characteristics previously obtained by gene signatures only. The results of the study elucidate the role of pseudogenes in two major BC subtypes, considering the contextual annotations. Put into a wider perspective, ncRNA brokers may help predictive functionalization studies targeted to new disease phenotypes, for instance those linked to the tumor microenvironment or metabolism, or those specifically involving metastasis. Overall, the approach may represent an in silico prioritization strategy toward the systems identification of new diagnostic and prognostic biomarkers.

Proteogenomic Analysis of Breast Cancer Transcriptomic and Proteomic Data, Using De Novo Transcript Assembly: Genome-Wide Identification of Novel Peptides and Clinical Implications

We have carried out proteogenomic analysis of the breast cancer transcriptomic and proteomic data, available at The Clinical Proteomic Tumor Analysis Consortium resource, to identify novel peptides arising from alternatively spliced events as well as other noncanonical expressions. We used a pipeline that consisted of de novo transcript assembly, six frame-translated custom database, and a combination of search engines to identify novel peptides. A portfolio of 4,387 novel peptide sequences initially identified was further screened through PepQuery validation tool (Clinical Proteomic Tumor Analysis Consortium), which yielded 1,558 novel peptides. We considered the dataset of 1,558 validated through PepQuery to understand their functional and clinical significance, leaving the rest to be further verified using other validation tools and approaches. The novel peptides mapped to the known gene sequences as well as to genomic regions yet undefined for translation, 580 novel peptides mapped to known protein-coding genes, 147 to non–protein-coding genes, and 831 belonged to novel translational sequences. The novel peptides belonging to protein-coding genes represented alternatively spliced events or 5′ or 3′ extensions, whereas others represented translation from pseudogenes, long noncoding RNAs, or novel peptides originating from uncharacterized protein-coding sequences—mostly from the intronic regions of known genes. Seventy-six of the 580 protein-coding genes were associated with cancer hallmark genes, which included key oncogenes, transcription factors, kinases, and cell surface receptors. Survival association analysis of the 76 novel peptide sequences revealed 10 of them to be significant, and we present a panel of six novel peptides, whose high expression was found to be strongly associated with poor survival of patients with human epidermal growth factor receptor 2–enriched subtype. Our analysis represents a landscape of novel peptides of different types that may be expressed in breast cancer tissues, whereas their presence in full-length functional proteins needs further investigations.

•

Novel protein variants and peptides from noncoding sequences are rapidly emerging.

•

Mining of mass spectrometry data using proteogenomic analysis reveals such entities.

•

Novel peptides from coding and noncoding sequences identified in breast cancer.

•

Novel peptides mapped to cancer hallmark genes in breast cancer.

•

Panel of novel peptides with prognostic potential found for HER2-enriched subtype.

Short open reading frames (sORFs) and microproteins: an update on their identification and validation measures

A short open reading frame (sORFs) constitutes ≤ 300 bases, encoding a microprotein or sORF-encoded protein (SEP) which comprises ≤ 100 amino acids. Traditionally dismissed by genome annotation pipelines as meaningless noise, sORFs were found to possess coding potential with ribosome profiling (RIBO-Seq), which unveiled sORF-based transcripts at various genome locations. Nonetheless, the existence of corresponding microproteins that are stable and functional was little substantiated by experimental evidence initially. With recent advancements in multi-omics, the identification, validation, and functional characterisation of sORFs and microproteins have become feasible. In this review, we discuss the history and development of an emerging research field of sORFs and microproteins. In particular, we focus on an array of bioinformatics and OMICS approaches used for predicting, sequencing, validating, and characterizing these recently discovered entities. These strategies include RIBO-Seq which detects sORF transcripts via ribosome footprints, and mass spectrometry (MS)-based proteomics for sequencing the resultant microproteins. Subsequently, our discussion extends to the functional characterisation of microproteins by incorporating CRISPR/Cas9 screen and protein–protein interaction (PPI) studies. Our review discusses not only detection methodologies, but we also highlight on the challenges and potential solutions in identifying and validating sORFs and their microproteins. The novelty of this review lies within its validation for the functional role of microproteins, which could contribute towards the future landscape of microproteomics.

Beyond Genes: Inclusion of Alternative Splicing and Alternative Polyadenylation to Assess the Genetic Architecture of Predisposition to Voluntary Alcohol Consumption in Brain of the HXB/BXH Recombinant Inbred Rat Panel

Post transcriptional modifications of RNA are powerful mechanisms by which eukaryotes expand their genetic diversity. For instance, researchers estimate that most transcripts in humans undergo alternative splicing and alternative polyadenylation. These splicing events produce distinct RNA molecules, which in turn yield distinct protein isoforms and/or influence RNA stability, translation, nuclear export, and RNA/protein cellular localization. Due to their pervasiveness and impact, we hypothesized that alternative splicing and alternative polyadenylation in brain can contribute to a predisposition for voluntary alcohol consumption. Using the HXB/BXH recombinant inbred rat panel (a subset of the Hybrid Rat Diversity Panel), we generated over one terabyte of brain RNA sequencing data (total RNA) and identified novel splice variants (via StringTie) and alternative polyadenylation sites (via aptardi) to determine the transcriptional landscape in the brains of these animals. After establishing an analysis pipeline to ascertain high quality transcripts, we quantitated transcripts and integrated genotype data to identify candidate transcript coexpression networks and individual candidate transcripts associated with predisposition to voluntary alcohol consumption in the two-bottle choice paradigm. For genes that were previously associated with this trait (e.g., Lrap, Ift81, and P2rx4) (Saba et al., Febs. J., 282, 3556–3578, Saba et al., Genes. Brain. Behav., 20, e12698), we were able to distinguish between transcript variants to provide further information about the specific isoforms related to the trait. We also identified additional candidate transcripts associated with the trait of voluntary alcohol consumption (i.e., isoforms of Mapkapk5, Aldh1a7, and Map3k7). Consistent with our previous work, our results indicate that transcripts and networks related to inflammation and the immune system in brain can be linked to voluntary alcohol consumption. Overall, we have established a pipeline for including the quantitation of alternative splicing and alternative polyadenylation variants in the transcriptome in the analysis of the relationship between the transcriptome and complex traits.

ANXA2P2: A Potential Immunological and Prognostic Signature in Ovarian Serous Cystadenocarcinoma via Pan-Carcinoma Synthesis

Background

Although the effect of pseudogene ANXA2P2 on some tumors has been reported in a few literatures, the therapeutic potential and prognostic value of ANXA2P2 in ovarian serous cystadenocarcinoma (OV) have not been elucidated.

Methods

The correlation for ANXA2P2 expression patterns to prognostic characteristics, tumor immune microenvironment, immune cell infiltration level, tumor mutation burden (TMB), tumor microsatellite instability (MSI), drug sensitivity, and pathway function enrichment were investigated in pan-carcinoma via TCGA and GTEx databases. Subsequently, the role of ANXA2P2 expression levels in the pathway enrichments and prognosis prediction in OV were further explored using weighted correlation network analysis (WGCNA) analysis, gene mutation analysis, and risk-independent prognostic analysis.

Results

ANXA2P2 was frequently overexpressed in a variety of tumors compared with normal tissues. The correlation analysis for prognostic characteristics, tumor immune microenvironment, immune cell infiltration level, TMB, MSI, drug sensitivity, and pathway function enrichment revealed that ANXA2P2 expression patterns might deal a significant impact on the pathogenesis, development, and prognosis of various tumors. Then, GSVA, GSEA, WGCNA, gene mutation, and independent prognostic analysis for OV have indicated that high expression in ANXA2P2 could be mostly enriched in TNF-α signaling-via-NF-κB, epithelial-mesenchymal transition, apical junction, IL-6-JAK STAT3 signaling, etc., which were also proved to act as crucial factors on tumorigenesis, development, invasion, and metastasis. The mutation of TP53 (94%), TTN (24%), and CSMD3 (9%) in the biological process of tumor had been confirmed by relevant studies. Finally, the independent prognostic analysis demonstrated that ANXA2P2 expression in OV contributes greatly to the dependability of 3- and 5-year survival prediction.

Conclusion

In summary, our findings might provide a helpful foundation for prospective explorative researches, afford new strategies for the clinical treatment, deal prognosis prediction, and give new hope for OV patients.

PΨFinder: a practical tool for the identification and visualization of novel pseudogenes in DNA sequencing data

BACKGROUND

Processed pseudogenes (PΨgs) are disabled gene copies that are transcribed and may affect expression of paralogous genes. Moreover, their insertion in the genome can disrupt the structure or the regulatory region of a gene, affecting its expression level. These events have been identified as occurring mutations during cancer development, thus being able to identify PΨgs and their location will improve their impact on diagnostic testing, not only in cancer but also in inherited disorders.

RESULTS

We have implemented PΨFinder (P-psy-finder), a tool that identifies PΨgs, annotates known ones and predicts their insertion site(s) in the genome. The tool screens alignment files and provides user-friendly summary reports and visualizations. To demonstrate its applicability, we scanned 218 DNA samples from patients screened for hereditary colorectal cancer. We detected 423 PΨgs distributed in 96% of the samples, comprising 7 different parent genes. Among these, we confirmed the well-known insertion site of the SMAD4-PΨg within the last intron of the SCAI gene in one sample. While for the ubiquitous CBX3-PΨg, present in 82.6% of the samples, we found it reversed inserted in the second intron of the C15ORF57 gene.

CONCLUSIONS

PΨFinder is a tool that can automatically identify novel PΨgs from DNA sequencing data and determine their location in the genome with high sensitivity (95.92%). It generates high quality figures and tables that facilitate the interpretation of the results and can guide the experimental validation. PΨFinder is a complementary analysis to any mutational screening in the identification of disease-causing mutations within cancer and other diseases.

The World of Pseudogenes: New Diagnostic and Therapeutic Targets in Cancers or Still Mystery Molecules?

Pseudogenes were once considered as “junk DNA”, due to loss of their functions as a result of the accumulation of mutations, such as frameshift and presence of premature stop-codons and relocation of genes to inactive heterochromatin regions of the genome. Pseudogenes are divided into two large groups, processed and unprocessed, according to their primary structure and origin. Only 10% of all pseudogenes are transcribed into RNAs and participate in the regulation of parental gene expression at both transcriptional and translational levels through senseRNA (sRNA) and antisense RNA (asRNA). In this review, about 150 pseudogenes in the different types of cancers were analyzed. Part of these pseudogenes seem to be useful in molecular diagnostics and can be detected in various types of biological material including tissue as well as biological fluids (liquid biopsy) using different detection methods. The number of pseudogenes, as well as their function in the human genome, is still unknown. However, thanks to the development of various technologies and bioinformatic tools, it was revealed so far that pseudogenes are involved in the development and progression of certain diseases, especially in cancer.

N6 -methyladenosine-Mediated Upregulation of WTAPP1 Promotes WTAP Translation and Wnt Signaling to Facilitate Pancreatic Cancer Progression

Pseudogenes may play important roles in cancer. Here, we explore the mechanism and function of a pseudogene WTAPP1 in the progress of pancreatic ductal adenocarcinoma (PDAC). WTAPP1 RNA was significantly elevated in PDAC and was associated with poor prognosis in patients. Overexpression of WTAPP1 RNA promoted PDAC proliferation and invasiveness in vitro and in vivo. Mechanistically, N ⁶-methyladenosine (m⁶A) modification stabilized WTAPP1 RNA via CCHC-type zinc finger nucleic-acid binding protein (CNBP), resulting in increased levels of WTAPP1 RNA in PDAC cells. Excessive WTAPP1 RNA bound its protein-coding counterpart WT1-associated protein (WTAP) mRNA and recruited more EIF3 translation initiation complex to promote WTAP translation. Increased WTAP protein enhanced the activation of Wnt signaling and provoked the malignant phenotypes of PDAC. Decreasing WTAPP1 RNA significantly suppressed the in vivo growth and metastasis of PDAC cell lines and patient-derived xenografts. These results indicate that m⁶A-mediated increases in WTAPP1 expression promote PDAC progression and thus may serve as a therapeutic target. SIGNIFICANCE: This study reveals how aberrant m⁶A modification of the WTAPP1 pseudogene results in increased translation of its protein-coding counterpart to promote Wnt signaling, which contributes to pancreatic cancer progression.

Processed pseudogenes: A substrate for evolutionary innovation: Retrotransposition contributes to genome evolution by propagating pseudogene sequences with rich regulatory potential throughout the genome

Processed pseudogenes may serve as a genetic reservoir for evolutionary innovation. Here, we argue that through the activity of long interspersed element-1 retrotransposons, processed pseudogenes disperse coding and noncoding sequences rich with regulatory potential throughout the human genome. While these sequences may appear to be non-functional, a lack of contemporary function does not prohibit future development of biological activity. Here, we discuss the dynamic evolution of certain processed pseudogenes into coding and noncoding genes and regulatory elements, and their implication in wide-ranging biological and pathological processes. Also see the video abstract here: https://youtu.be/iUY_mteVoPI.

CircRNAs in gastric cancer: current research and potential clinical implications

Gastric cancer (GC) has a dismal prognosis and is also one of the most commonly diagnosed malignancies worldwide. circRNAs are covalently closed circular RNA molecules without 5'-cap and a 3'-tail, currently listed among the broad noncoding RNA family. circRNAs participate in a variety of pathophysiological processes relevant to human diseases, especially malignancies, including GC. Compelling evidence has shown that circRNAs can function by sponging miRNAs, interacting with RNA binding proteins, and encoding proteins or peptides. Yet, our current understanding of these RNA circles remains very limited. Here, we overview the biogenesis, characteristics, functions, and degradation of circRNAs. Moreover, we give an account of the circRNAs that have been linked with GC, describing their functions and mechanisms of action in the context of GC. Next, we discuss the potential value of circRNAs as diagnostic or prognostic GC biomarkers and summarize future prospects, important questions, and challenges of circRNA-based therapeutics.

Pseudogenes and their potential functions in hematopoiesis

Pseudogenes are DNA regions comprising defective copies of functional genes, the majority of which were generated by RNA- or DNA-level duplications. They exist across almost all forms of life and account for about one-quarter of the annotated genes in the human genome. Although these have been considered nonfunctional for decades, a growing number of pseudogenes have been found to be transcribed and to play crucial regulatory roles. Accumulating evidence indicates that they regulate gene expression through molecular interactions at the protein, RNA, and DNA levels. However, pseudogenes are often excluded in multiple genomewide analyses and functional screening, and their biological activities remain to be systematically disclosed. Here, we summarize the features of and progress of research on pseudogenes, in addition to discussing what is unknown about these genetic elements. Our previous findings, together with evidence of their poor conservation, prompted us to propose that pseudogenes may contribute to primate- or human-specific regulation, especially in hematopoiesis.

Systematic functional interrogation of human pseudogenes using CRISPRi

BACKGROUND

The human genome encodes over 14,000 pseudogenes that are evolutionary relics of protein-coding genes and commonly considered as nonfunctional. Emerging evidence suggests that some pseudogenes may exert important functions. However, to what extent human pseudogenes are functionally relevant remains unclear. There has been no large-scale characterization of pseudogene function because of technical challenges, including high sequence similarity between pseudogene and parent genes, and poor annotation of transcription start sites.

RESULTS

To overcome these technical obstacles, we develop an integrated computational pipeline to design the first genome-wide library of CRISPR interference (CRISPRi) single-guide RNAs (sgRNAs) that target human pseudogene promoter-proximal regions. We perform the first pseudogene-focused CRISPRi screen in luminal A breast cancer cells and reveal approximately 70 pseudogenes that affect breast cancer cell fitness. Among the top hits, we identify a cancer-testis unitary pseudogene, MGAT4EP, that is predominantly localized in the nucleus and interacts with FOXA1, a key regulator in luminal A breast cancer. By enhancing the promoter binding of FOXA1, MGAT4EP upregulates the expression of oncogenic transcription factor FOXM1. Integrative analyses of multi-omic data from the Cancer Genome Atlas (TCGA) reveal many unitary pseudogenes whose expressions are significantly dysregulated and/or associated with overall/relapse-free survival of patients in diverse cancer types.

CONCLUSIONS

Our study represents the first large-scale study characterizing pseudogene function. Our findings suggest the importance of nuclear function of unitary pseudogenes and underscore their underappreciated roles in human diseases. The functional genomic resources developed here will greatly facilitate the study of human pseudogene function.

Pseudogenes: Four Decades of Discovery

A pseudogene is defined as a genomic DNA sequence that looks like a mutated or truncated version of a known functional gene. Nearly four decades since their first discovery it has been estimated that between ~12,000 and ~20,000 pseudogenes exist in the human genome. Early efforts to characterize functions for pseudogenes were unsuccessful, thus they were considered functionless relics of evolutionary selection, junk DNA or genetic fossils. Remarkably, an increasing number of pseudogenes have been reported to be expressed as RNA transcripts above and beyond levels considered accidental or spurious transcription. There is emerging evidence that some expressed pseudogene transcripts have biological functions and should be defined as a subclass of functional long noncoding RNAs (lncRNA). In this introductory chapter, I briefly summarize the history and the current knowledge of pseudogenes, and highlight the emerging functions of some pseudogenes in human biology and disease. This second iteration of Pseudogenes in Methods in Molecular Biology highlights new methodological approaches to investigate this intriguing family of lncRNAs and the extent of their biological function.

Enhancing Gonadotrope Gene Expression Through Regulatory lncRNAs

The world of long non-coding RNAs (lncRNAs) has opened up massive new prospects in understanding the regulation of gene expression. Not only are there seemingly almost infinite numbers of lncRNAs in the mammalian cell, but they have highly diverse mechanisms of action. In the nucleus, some are chromatin-associated, transcribed from transcriptional enhancers (eRNAs) and/or direct changes in the epigenetic landscape with profound effects on gene expression. The pituitary gonadotrope is responsible for activation of reproduction through production and secretion of appropriate levels of the gonadotropic hormones. As such, it exemplifies a cell whose function is defined through changes in developmental and temporal patterns of gene expression, including those that are hormonally induced. Roles for diverse distal regulatory elements and eRNAs in gonadotrope biology have only just begun to emerge. Here, we will present an overview of the different kinds of lncRNAs that alter gene expression, and what is known about their roles in regulating some of the key gonadotrope genes. We will also review various screens that have detected differentially expressed pituitary lncRNAs associated with changes in reproductive state and those whose expression is found to play a role in gonadotrope-derived nonfunctioning pituitary adenomas. We hope to shed light on this exciting new field, emphasize the open questions, and encourage research to illuminate the roles of lncRNAs in various endocrine systems.

Diverse heterochromatin-associated proteins repress distinct classes of genes and repetitive elements

Heterochromatin, typically marked by histone H3 trimethylation at lysine 9 (H3K9me3) or lysine 27 (H3K27me3), represses different protein-coding genes in different cells, as well as repetitive elements. The basis for locus specificity is unclear. Previously, we identified 172 proteins that are embedded in sonication-resistant heterochromatin (srHC) harbouring H3K9me3. Here, we investigate in humans how 97 of the H3K9me3-srHC proteins repress heterochromatic genes. We reveal four groups of srHC proteins that each repress many common genes and repeat elements. Two groups repress H3K9me3-embedded genes with different extents of flanking srHC, one group is specific for srHC genes with H3K9me3 and H3K27me3, and one group is specific for genes with srHC as the primary feature. We find that the enhancer of rudimentary homologue (ERH) is conserved from Schizosaccharomyces pombe in repressing meiotic genes and, in humans, now represses other lineage-specific genes and repeat elements. The study greatly expands our understanding of H3K9me3-based gene repression in vertebrates.

Dealing with Pseudogenes in Molecular Diagnostics in the Next Generation Sequencing Era

Presence of pseudogenes is a dreadful issue in next generation sequencing (NGS), because their contamination can interfere with the detection of variants in the genuine gene and generate false positive and false negative variants.In this chapter we focus on issues related to the application of NGS strategies for analysis of genes with pseudogenes in a clinical setting. The degree to which a pseudogene impacts the ability to accurately detect and map variants in its parent gene depends on the degree of similarity (homology) with the parent gene itself. Hereby, target enrichment and mapping strategies are crucial factors to avoid "contaminating" pseudogene sequences. For target enrichment, we describe advantages and disadvantages of PCR- and capture-based strategies. For mapping strategies, we discuss crucial parameters that need to be considered to accurately distinguish sequences of functional genes from pseudogenic sequences. Finally, we discuss some examples of genes associated with Mendelian disorders, for which interesting NGS approaches are described to avoid interference with pseudogene sequences.

Limbic Expression of mRNA Coding for Chemoreceptors in Human Brain-Lessons from Brain Atlases

Animals strongly rely on chemical senses to uncover the outside world and adjust their behaviour. Chemical signals are perceived by facial sensitive chemosensors that can be clustered into three families, namely the gustatory (TASR), olfactory (OR, TAAR) and pheromonal (VNR, FPR) receptors. Over recent decades, chemoreceptors were identified in non-facial parts of the body, including the brain. In order to map chemoreceptors within the encephalon, we performed a study based on four brain atlases. The transcript expression of selected members of the three chemoreceptor families and their canonical partners was analysed in major areas of healthy and demented human brains. Genes encoding all studied chemoreceptors are transcribed in the central nervous system, particularly in the limbic system. RNA of their canonical transduction partners (G proteins, ion channels) are also observed in all studied brain areas, reinforcing the suggestion that cerebral chemoreceptors are functional. In addition, we noticed that: (i) bitterness-associated receptors display an enriched expression, (ii) the brain is equipped to sense trace amines and pheromonal cues and (iii) chemoreceptor RNA expression varies with age, but not dementia or brain trauma. Extensive studies are now required to further understand how the brain makes sense of endogenous chemicals.

Positive natural selection of N6-methyladenosine on the RNAs of processed pseudogenes

BACKGROUND

Canonical nonsense-mediated decay (NMD) is an important splicing-dependent process for mRNA surveillance in mammals. However, processed pseudogenes are not able to trigger NMD due to their lack of introns. It is largely unknown whether they have evolved other surveillance mechanisms.

RESULTS

Here, we find that the RNAs of pseudogenes, especially processed pseudogenes, have dramatically higher m6A levels than their cognate protein-coding genes, associated with de novo m6A peaks and motifs in human cells. Furthermore, pseudogenes have rapidly accumulated m6A motifs during evolution. The m6A sites of pseudogenes are evolutionarily younger than neutral sites and their m6A levels are increasing, supporting the idea that m6A on the RNAs of pseudogenes is under positive selection. We then find that the m6A RNA modification of processed, rather than unprocessed, pseudogenes promotes cytosolic RNA degradation and attenuates interference with the RNAs of their cognate protein-coding genes. We experimentally validate the m6A RNA modification of two processed pseudogenes, DSTNP2 and NAP1L4P1, which promotes the RNA degradation of both pseudogenes and their cognate protein-coding genes DSTN and NAP1L4. In addition, the m6A of DSTNP2 regulation of DSTN is partially dependent on the miRNA miR-362-5p.

CONCLUSIONS

Our discovery reveals a novel evolutionary role of m6A RNA modification in cleaning up the unnecessary processed pseudogene transcripts to attenuate their interference with the regulatory network of protein-coding genes.

Long-read cDNA sequencing identifies functional pseudogenes in the human transcriptome

Pseudogenes are gene copies presumed to mainly be functionless relics of evolution due to acquired deleterious mutations or transcriptional silencing. Using deep full-length PacBio cDNA sequencing of normal human tissues and cancer cell lines, we identify here hundreds of novel transcribed pseudogenes expressed in tissue-specific patterns. Some pseudogene transcripts have intact open reading frames and are translated in cultured cells, representing unannotated protein-coding genes. To assess the biological impact of noncoding pseudogenes, we CRISPR-Cas9 delete the nucleus-enriched pseudogene PDCL3P4 and observe hundreds of perturbed genes. This study highlights pseudogenes as a complex and dynamic component of the human transcriptional landscape.

Increased expression of peptides from non-coding genes in cancer proteomics datasets suggests potential tumor neoantigens

Neoantigen-based immunotherapy has yielded promising results in clinical trials. However, it is limited to tumor-specific mutations, and is often tailored to individual patients. Identifying suitable tumor-specific antigens is still a major challenge. Previous proteogenomics studies have identified peptides encoded by predicted non-coding sequences in human genome. To investigate whether tumors express specific peptides encoded by non-coding genes, we analyzed published proteomics data from five cancer types including 933 tumor samples and 275 matched normal samples and compared these to data from 31 different healthy human tissues. Our results reveal that many predicted non-coding genes such as DGCR9 and RHOXF1P3 encode peptides that are overexpressed in tumors compared to normal controls. Furthermore, from the non-coding genes-encoded peptides specifically detected in cancers, we predict a large number of “dark antigens” (neoantigens from non-coding genomic regions), which may provide an alternative source of neoantigens beyond standard tumor specific mutations.

Rong Xiang et al. analyze the expression of non-coding genes encoded peptides in publicly-available proteomics data from five cancer types and matched controls. They identify peptides from non-coding genes including DGCR9 and RHOXF1P3 that are upregulated in tumors compared to controls, suggesting that non-coding gene-encoded peptides may be a source of neoantigens in some cancers.

Circulating Biomarkers for Early Stage Non-Small Cell Lung Carcinoma Detection: Supplementation to Low-Dose Computed Tomography

Lung cancer is currently the leading cause of cancer death in both developing and developed countries. Given that lung cancer has poor prognosis in later stages, it is essential to achieve an early diagnosis to maximize patients’ overall survival. Non-small cell lung cancer (NSCLC) is the most common form of primary lung cancer in both smokers and non-smokers. The current standard screening method, low‐dose computed tomography (LDCT), is the only radiological method that demonstrates to have mortality benefits across multiple large randomized clinical trials (RCT). However, these RCTs also found LDCT to have a significant false positive rate that results in unnecessary invasive biopsies being performed. Due to the lack of both sensitive and specific screening methods for the early detection of lung cancer, there is an urgent need for alternative minimally or non-invasive biomarkers that may provide diagnostic, and/or prognostic information. This has led to the identification of circulating biomarkers that can be readily detectable in blood and have been extensively studied as prognosis markers. Circulating microRNA (miRNA) in particular has been investigated for these purposes as an augmentation to LDCT, or as direct diagnosis of lung cancer. There is, however, a lack of consensus across the studies on which miRNAs are the most clinically useful. Besides miRNA, other potential circulating biomarkers include circulating tumor cells (CTCs), circulating tumor DNA (ctDNAs) and non-coding RNAs (ncRNAs). In this review, we provide the current outlook of several of these biomarkers for the early diagnosis of NSCLC.

Circular RNAs and Cardiovascular Regeneration

circular RNAs (circRNAs) are a type of non-coding RNAs that are widely present in eukaryotic cells. They have the characteristics of stable structure, high abundance, and cell or tissue specific expression. circRNAs are single-stranded RNAs that are covalently back spliced to form closed circular loops. They may participate in gene expression and regulation through a variety of action modes. circRNAs can encode proteins or function by acting as miRNA sponges for protein translation. Since 2016, a growing number of research studies have shown that circRNAs play important role in the pathogenesis of cardiovascular disease. With the construction of circRNA database, the differential expression of circRNAs in the heart tissue samples from different species and the gradual elucidation of its mode of action in disease may become an ideal diagnosis biomarker and an effective therapeutic target. What can be expected surely has a broader application prospect. In this review, we summarize recent publications on circRNA biogenesis, expression profiles, functions, and the most recent studies of circRNAs in the field of cardiovascular diseases with special emphasis on cardiac regeneration.

Emerging roles of circular RNAs in non‑small cell lung cancer (Review)

Circular RNAs (circRNAs) are a class of novel endogenous transcripts with limited protein‑coding abilities. CircRNAs have been demonstrated to function as critical regulators of tumor development and distant metastasis through binding to microRNAs (miRNAs) and interacting with RNA‑binding proteins, thereby regulating transcription and translation. Emerging evidence has illustrated that certain circRNAs can serve as biomarkers for diagnosis and prognosis of cancer, and/or serve as potential therapeutic targets. Expression of functional circRNAs is commonly dysregulated in cancer and this is correlated with advanced Tumor‑Node‑Metastasis stage, lymph node status, distant metastasis, poor differentiation and shorter overall survival of cancer patients. Recently, an increasing number of studies have shown that circRNAs are closely associated with NSCLC. Functional experiments have revealed that circRNAs are intricately associated with the pathological progression of NSCLC. The present review provides an overview of the regulatory effect of circRNAs in the development and progression of NSCLC, taking into consideration various physiological and pathological processes, such as proliferation, apoptosis, invasion and migration, and their potential value as biomarkers and therapeutic targets.

Identification of Transcriptomic Differences between Lower Extremities Arterial Disease, Abdominal Aortic Aneurysm and Chronic Venous Disease in Peripheral Blood Mononuclear Cells Specimens

Several human tissues are investigated in studies of molecular biomarkers associated with diseases development. Special attention is focused on the blood and its components due to combining abundant information about systemic responses to pathological processes as well as high accessibility. In the current study, transcriptome profiles of peripheral blood mononuclear cells (PBMCs) were used to compare differentially expressed genes between patients with lower extremities arterial disease (LEAD), abdominal aortic aneurysm (AAA) and chronic venous disease (CVD). Gene expression patterns were generated using the Ion S5XL next-generation sequencing platform and were analyzed using DESeq2 and UVE-PLS methods implemented in R programming software. In direct pairwise analysis, 21, 58 and 10 differentially expressed genes were selected from the comparison of LEAD vs. AAA, LEAD vs. CVD and AAA vs. CVD patient groups, respectively. Relationships between expression of dysregulated genes and age, body mass index, creatinine levels, hypertension and medication were identified using Spearman rank correlation test and two-sided Mann-Whitney U test. The functional analysis, performed using DAVID website tool, provides potential implications of selected genes in pathological processes underlying diseases studied. Presented research provides new insight into differences of pathogenesis in LEAD, AAA and CVD, and selected genes could be considered as potential candidates for biomarkers useful in diagnosis and differentiation of studied diseases.

Circular RNAs: Their Role in the Pathogenesis and Orchestration of Breast Cancer

As one of the most frequently occurring malignancies in women, breast cancer (BC) is still an enormous threat to women all over the world. The high mortality rates in BC patients are associated with BC recurrence, metastatic progression to distant organs, and therapeutic resistance. Circular RNAs (circRNAs), belonging to the non-coding RNAs (ncRNAs), are connected end to end to form covalently closed single-chain circular molecules. CircRNAs are widely found in different species and a variety of human cells, with the features of diversity, evolutionary conservation, stability, and specificity. CircRNAs are emerging important participators in multiple diseases, including cardiovascular disease, inflammation, and cancer. Recent studies have shown that circRNAs are involved in BC progress by regulating gene expression at the transcriptional or post-transcriptional level via binding to miRNAs then inhibiting their function, suggesting that circRNAs may be potential targets for early diagnosis, treatment, and prognosis of BC. Herein, in this article, we have reviewed and summarized the current studies about the biogenesis, features, and functions of circRNAs. More importantly, we emphatically elucidate the pivotal functions and mechanisms of circRNAs in BC growth, metastasis, diagnosis, and drug resistance. Deciphering the complex networks, especially the circRNA-miRNA target gene axis, will endow huge potentials in developing therapeutic strategies for combating BC.

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.

Pseudogenes in Cardiovascular Disease

Cardiovascular disease is the main disease that affects human life span. In recent years, the disease has been increasingly addressed at the molecular levels, for example, pseudogenes are now known to be involved in the pathogenesis and development of cardiovascular diseases. Pseudogenes are non-coding homologs of protein-coding genes and were once called “junk gene.” Since they are highly homologous to their functional parental genes, it is somewhat difficult to distinguish them. With the development of sequencing technology and bioinformatics, pseudogenes have become readily identifiable. Recent studies indicate that pseudogenes are closely related to cardiovascular diseases. This review provides an overview of pseudogenes and their roles in the pathogenesis of cardiovascular diseases. This new knowledge adds to our understanding of cardiovascular disease at the molecular level and will help develop new biomarkers and therapeutic approaches designed to prevent and treat the disease.

Cancer, Retrogenes, and Evolution

This review summarizes the knowledge about retrogenes in the context of cancer and evolution. The retroposition, in which the processed mRNA from parental genes undergoes reverse transcription and the resulting cDNA is integrated back into the genome, results in additional copies of existing genes. Despite the initial misconception, retroposition-derived copies can become functional, and due to their role in the molecular evolution of genomes, they have been named the “seeds of evolution”. It is convincing that retrogenes, as important elements involved in the evolution of species, also take part in the evolution of neoplastic tumors at the cell and species levels. The occurrence of specific “resistance mechanisms” to neoplastic transformation in some species has been noted. This phenomenon has been related to additional gene copies, including retrogenes. In addition, the role of retrogenes in the evolution of tumors has been described. Retrogene expression correlates with the occurrence of specific cancer subtypes, their stages, and their response to therapy. Phylogenetic insights into retrogenes show that most cancer-related retrocopies arose in the lineage of primates, and the number of identified cancer-related retrogenes demonstrates that these duplicates are quite important players in human carcinogenesis.

Long non-coding antisense RNA HYOU1-AS is essential to human breast cancer development through competitive binding hnRNPA1 to promote HYOU1 expression

Triple negative breast cancer (TNBC) has poor prognosis due to lack of biomarker and therapeutic target. Emerging research has revealed long noncoding RNAs (lncRNAs) are involved in breast cancer progression, but their functions and regulatory mechanisms remain poorly understood, especially in TNBC. In this study, we performed lncRNA microarray analysis of five TNBC samples and their matched normal tissues, and discovered a number of differentially expressed lncRNAs. We identified an antisense lncRNA, HYOU1-AS, which is transcribed from the opposite strand of the hypoxia up-regulated 1 (HYOU1) gene, enriched in the nucleus and highly expressed in TNBC. HYOU1-AS knockdown could inhibit the proliferation and migration of the TNBC MDA-MB-231 cells, and reduce their xenograft tumor formation in nude mice. In mechanistic studies, we found that HYOU1-AS could promote the expression of HYOU1, a proliferative gene, through competitively binding to hnRNPA1, an RNA-binding protein, to relieve its post-transcriptional inhibition of the HYOU1 mRNA. Consistently, increased HYOU1 levels correlated with poor clinical outcomes of breast cancer patients based on our study of the TCGA database. Overall, our data indicated that the lncRNA HYOU1-AS promoted TNBC progression through upregulating HYOU1.

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.