The GENCODE pseudogene resource

Construction of a prognostic model via WGCNA combined with the LASSO algorithm for stomach adenocarcinoma patients

Objective

This study aimed to identify prognostic signatures to predict the prognosis of patients with stomach adenocarcinoma (STAD), which is necessary to improve poor prognosis and offer possible treatment strategies for STAD patients.

Methods

The overlapping genes between the key model genes that were screened by the weighted gene co-expression network analysis (WGCNA) and differentially expressed genes (DEGs) whose expression was different with significance between normal and tumor tissues were extracted to serve as co-expression genes. Then, enrichment analysis was performed on these genes. Furthermore, the least absolute shrinkage and selection operator (LASSO) regression was performed to screen the hub genes among overlapping genes. Finally, we constructed a model to explore the influence of polygenic risk scores on the survival probability of patients with STAD, and interaction effect and mediating analyses were also performed.

Results

DEGs included 2,899 upregulated genes and 2,896 downregulated genes. After crossing the DEGs and light-yellow module genes that were obtained by WGCNA, a total of 39 overlapping genes were extracted. The gene enrichment analysis revealed that these genes were enriched in the prion diseases, biosynthesis of unsaturated fatty acids, RNA metabolic process, hydrolase activity, etc. PIP5K1P1, PTTG3P, and SNORD15B were determined by LASSO-Cox. The prognostic prediction of the three-gene model was established. The Cox regression analysis showed that the comprehensive risk score for three genes was an independent prognosis factor.

Conclusion

PIP5K1P1, PTTG3P, and SNORD15B are related to the prognosis and overall survival of patients. The three-gene risk model constructed has independent prognosis predictive ability for STAD.

The pseudogene GBP1P1 suppresses influenza A virus replication by acting as a protein decoy for DHX9

Recently, substantial evidence has demonstrated that pseudogene-derived long noncoding RNAs (lncRNAs) as regulatory RNAs have been implicated in basic physiological processes and disease development through multiple modes of functional interaction with DNA, RNA, and proteins. Here, we report an important role for GBP1P1, the pseudogene of guanylate-binding protein 1, in regulating influenza A virus (IAV) replication in A549 cells. GBP1P1 was dramatically upregulated after IAV infection, which is controlled by JAK/STAT signaling. Functionally, ectopic expression of GBP1P1 in A549 cells resulted in significant suppression of IAV replication. Conversely, silencing GBP1P1 facilitated IAV replication and virus production, suggesting that GBP1P1 is one of the interferon-inducible antiviral effectors. Mechanistically, GBP1P1 is localized in the cytoplasm and functions as a sponge to trap DHX9 (DExH-box helicase 9), which subsequently restricts IAV replication. Together, these studies demonstrate that GBP1P1 plays an important role in antagonizing IAV replication.IMPORTANCELong noncoding RNAs (lncRNAs) are extensively expressed in mammalian cells and play a crucial role as regulators in various biological processes. A growing body of evidence suggests that host-encoded lncRNAs are important regulators involved in host-virus interactions. Here, we define a novel function of GBP1P1 as a decoy to compete with viral mRNAs for DHX9 binding. We demonstrate that GBP1P1 induction by IAV is mediated by JAK/STAT activation. In addition, GBP1P1 has the ability to inhibit IAV replication. Importantly, we reveal that GBP1P1 acts as a decoy to bind and titrate DHX9 away from viral mRNAs, thereby attenuating virus production. This study provides new insight into the role of a previously uncharacterized GBP1P1, a pseudogene-derived lncRNA, in the host antiviral process and a further understanding of the complex GBP network.

Combining single-cell profiling and functional analysis explores the role of pseudogenes in human early embryonic development

More and more studies have demonstrated that pseudogenes possess coding ability, and the functions of their transcripts in the development of diseases have been partially revealed. However, the role of pseudogenes in maintenance of normal physiological states and life activities has long been neglected. Here, we identify pseudogenes that are dynamically expressed during human early embryogenesis, showing different expression pattern from that of adult tissues. We explore the expression correlation between pseudogenes and the parent genes, part due to their shared gene regulatory elements or the potential regulation network between them. The essential role of three pseudogenes, PI4KAP1, TMED10P1, and FBXW4P1, in maintaining self-renewal of human embryonic stem cells is demonstrated. We further find that the three pseudogenes might perform their regulatory functions by binding to proteins or microRNAs. The pseudogene-related single-nucleotide polymorphisms are significantly associated with human congenital disease, further illustrating their importance during early embryonic development. Overall, this study is an excavation and exploration of functional pseudogenes during early human embryonic development, suggesting that pseudogenes are not only capable of being specifically activated in pathological states, but also play crucial roles in the maintenance of normal physiological states.

The annotation of GBA1 has been concealed by its protein-coding pseudogene GBAP1

Mutations in GBA1 cause Gaucher disease and are the most important genetic risk factor for Parkinson's disease. However, analysis of transcription at this locus is complicated by its highly homologous pseudogene, GBAP1. We show that >50% of short RNA-sequencing reads mapping to GBA1 also map to GBAP1. Thus, we used long-read RNA sequencing in the human brain, which allowed us to accurately quantify expression from both GBA1 and GBAP1. We discovered significant differences in expression compared to short-read data and identify currently unannotated transcripts of both GBA1 and GBAP1. These included protein-coding transcripts from both genes that were translated in human brain, but without the known lysosomal function-yet accounting for almost a third of transcription. Analyzing brain-specific cell types using long-read and single-nucleus RNA sequencing revealed region-specific variations in transcript expression. Overall, these findings suggest nonlysosomal roles for GBA1 and GBAP1 with implications for our understanding of the role of GBA1 in health and disease.

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.

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.

Sequence Alignment/Map format: a comprehensive review of approaches and applications

The Sequence Alignment/Map (SAM) format file is the text file used to record alignment information. Alignment is the core of sequencing analysis, and downstream tasks accept mapping results for further processing. Given the rapid development of the sequencing industry today, a comprehensive understanding of the SAM format and related tools is necessary to meet the challenges of data processing and analysis. This paper is devoted to retrieving knowledge in the broad field of SAM. First, the format of SAM is introduced to understand the overall process of the sequencing analysis. Then, existing work is systematically classified in accordance with generation, compression and application, and the involved SAM tools are specifically mined. Lastly, a summary and some thoughts on future directions are provided.

LncRNA PTTG3P promotes tumorigenesis and metastasis of NSCLC by binding with ILF3 to maintain mRNA stability and form a positive feedback loop with E2F1

Non-small cell lung cancer (NSCLC) is a highly lethal disease worldwide. We found the pseudogene-derived lncRNA PTTG3P is upregulated in NSCLC and associated with larger tumor size, advanced staging, and poor prognosis. This study investigated the oncogenic roles and mechanisms of PTTG3P in NSCLC. We demonstrate that PTTG3P promoted NSCLC cell proliferation, migration, tumorigenesis, and metastasis while inhibiting apoptosis in vitro and in vivo. Mechanistically, PTTG3P formed an RNA-protein complex with ILF3 to maintain MAP2K6 and E2F1 mRNA stability, two oncogenic factors involved in NSCLC progression. RNA-seq revealed MAP2K6 and E2F1 were downregulated upon PTTG3P knockdown. RIP and RNA stability assays showed PTTG3P/ILF3 interaction stabilized MAP2K6 and E2F1 transcripts. Interestingly, E2F1 transcriptionally upregulated PTTG3P by binding its promoter, forming a positive feedback loop. Knockdown of E2F1 or PTTG3P attenuated their mutual regulatory effects on cell growth and migration. Thus, a PTTG3P/ILF3/E2F1 axis enhances oncogene expression to promote NSCLC pathogenesis. Our study reveals PTTG3P exerts oncogenic functions in NSCLC via mRNA stabilization and a feedback loop, highlighting its potential as a prognostic biomarker and therapeutic target.

Genome annotation: From human genetics to biodiversity genomics

Within the next decade, the genomes of 1.8 million eukaryotic species will be sequenced. Identifying genes in these sequences is essential to understand the biology of the species. This is challenging due to the transcriptional complexity of eukaryotic genomes, which encode hundreds of thousands of transcripts of multiple types. Among these, a small set of protein-coding mRNAs play a disproportionately large role in defining phenotypes. Due to their sequence conservation, orthology can be established, making it possible to define the universal catalog of eukaryotic protein-coding genes. This catalog should substantially contribute to uncovering the genomic events underlying the emergence of eukaryotic phenotypes. This piece briefly reviews the basics of protein-coding gene prediction, discusses challenges in finalizing annotation of the human genome, and proposes strategies for producing annotations across the eukaryotic Tree of Life. This lays the groundwork for obtaining the catalog of all genes-the Earth's code of life.

The X-factor in ART: does the use of assisted reproductive technologies influence DNA methylation on the X chromosome?

BACKGROUND

Assisted reproductive technologies (ART) may perturb DNA methylation (DNAm) in early embryonic development. Although a handful of epigenome-wide association studies of ART have been published, none have investigated CpGs on the X chromosome. To bridge this knowledge gap, we leveraged one of the largest collections of mother-father-newborn trios of ART and non-ART (natural) conceptions to date to investigate sex-specific DNAm differences on the X chromosome. The discovery cohort consisted of 982 ART and 963 non-ART trios from the Norwegian Mother, Father, and Child Cohort Study (MoBa). To verify our results from the MoBa cohort, we used an external cohort of 149 ART and 58 non-ART neonates from the Australian 'Clinical review of the Health of adults conceived following Assisted Reproductive Technologies' (CHART) study. The Illumina EPIC array was used to measure DNAm in both datasets. In the MoBa cohort, we performed a set of X-chromosome-wide association studies ('XWASs' hereafter) to search for sex-specific DNAm differences between ART and non-ART newborns. We tested several models to investigate the influence of various confounders, including parental DNAm. We also searched for differentially methylated regions (DMRs) and regions of co-methylation flanking the most significant CpGs. Additionally, we ran an analogous model to our main model on the external CHART dataset.

RESULTS

In the MoBa cohort, we found more differentially methylated CpGs and DMRs in girls than boys. Most of the associations persisted after controlling for parental DNAm and other confounders. Many of the significant CpGs and DMRs were in gene-promoter regions, and several of the genes linked to these CpGs are expressed in tissues relevant for both ART and sex (testis, placenta, and fallopian tube). We found no support for parental DNAm-dependent features as an explanation for the observed associations in the newborns. The most significant CpG in the boys-only analysis was in UBE2DNL, which is expressed in testes but with unknown function. The most significant CpGs in the girls-only analysis were in EIF2S3 and AMOT. These three loci also displayed differential DNAm in the CHART cohort.

CONCLUSIONS

Genes that co-localized with the significant CpGs and DMRs associated with ART are implicated in several key biological processes (e.g., neurodevelopment) and disorders (e.g., intellectual disability and autism). These connections are particularly compelling in light of previous findings indicating that neurodevelopmental outcomes differ in ART-conceived children compared to those naturally conceived.

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.

Global DNA and protein interactomes of FLT1P1 (Fms-related tyrosine kinase 1 pseudogene 1) revealed its molecular regulatory functions associated with preeclampsia

BACKGROUND

Preeclampsia (PE) is one of the most serious pregnancy complications with unknown pathogenesis. Emerging evidence has demonstrated that Fms-related tyrosine kinase 1 (FLT1) is highly involved in PE development. As a pseudogene of FLT1, FLT1P1 increased in PE samples. However, its functions remain largely unknown.

METHODS AND RESULTS

In this study, co-expression analysis was performed to identify the potential target genes of FTL1P1. Then chromatin isolation using RNA purification (ChIRP) method was employed to explore the interactomes of FLT1P1, including interacting with DNA fragments and proteins. We found that in PE samples, both FLT1P1 and FLT1 were highly expressed and closely correlated. ChIRP-protein data revealed that FLT1P1 interacts with translation- and transcription-related proteins, including 4 transcription factors (TFs). ChIRP-DNA analysis revealed that FLT1P1 preferentially interacted with DNA fragments downstream of transcription start sites (TSSs). Functional analysis of its interacting genes revealed that they were enriched in transcriptional regulation and apoptosis-related pathways. Twenty-six TFs, including CREB1 and SRF, were extracted from the potential FLT1P1-interacting gene sets and were potential targets of FLT1P1. CREB1 could bind to FLT1 promoter, and was negatively correlated with FLT1 at the expression level, making it a potential regulator of FLT1.

CONCLUSIONS

Our study extensively investigated the interactome profiles of FLT1P1, especially the prompter region of TF gene CREB1, and revealed the potential molecular regulatory mechanisms of FLT1 expression in PE samples. Our results provide a novel view of PE pathogenesis, and suggest that FLT1P1 could serve as a potential therapeutic target in PE diagnosis and treatment.

Competing endogenous RNA (ceRNA) networks in Parkinson's disease: A systematic review

Parkinson's disease (PD) is a distinctive clinical syndrome with several causes and clinical manifestations. Aside from an infectious cause, PD is a rapidly developing neurological disorder with a global rise in frequency. Notably, improved knowledge of molecular pathways and the developing novel diagnostic methods may result in better therapy for PD patients. In this regard, the amount of research on ceRNA axes is rising, highlighting the importance of these axes in PD. CeRNAs are transcripts that cross-regulate one another via competition for shared microRNAs (miRNAs). These transcripts may be either coding RNAs (mRNAs) or non-coding RNAs (ncRNAs). This research used a systematic review to assess validated loops of ceRNA in PD. The Prisma guideline was used to conduct this systematic review, which entailed systematically examining the articles of seven databases. Out of 309 entries, forty articles met all criteria for inclusion and were summarized in the appropriate table. CeRNA axes have been described through one of the shared vital components of the axes, including lncRNAs such as NEAT1, SNHG family, HOTAIR, MALAT1, XIST, circRNAs, and lincRNAs. Understanding the multiple aspects of this regulatory structure may aid in elucidating the unknown causal causes of PD and providing innovative molecular therapeutic targets and medical fields.

Genetics in prenatal diagnosis

The options for prenatal genetic testing have evolved rapidly in the past decade, and advances in sequencing technology now allow genetic diagnoses to be made down to the single-base-pair level, even before the birth of the child. This offers women the opportunity to obtain information regarding the foetus, thereby empowering them to make informed decisions about their pregnancy. As genetic testing becomes increasingly available to women, clinician knowledge and awareness of the options available to women is of great importance. Additionally, comprehensive pretest and posttest genetic counselling about the advantages, pitfalls and limitations of genetic testing should be provided to all women. This review article aims to cover the range of genetic tests currently available in prenatal screening and diagnosis, their current applications and limitations in clinical practice as well as what the future holds for prenatal genetics.

Non-coding RNAs in stroke pathology, diagnostics, and therapeutics

Ischemic stroke is a leading cause of death and disability worldwide. Methods to alleviate functional deficits after ischemic stroke focus on restoration of cerebral blood flow to the affected area. However, pharmacological or surgical methods such as thrombolysis and thrombectomy have a narrow effective window. Harnessing and manipulating neurochemical processes of recovery may provide an alternative to these methods. Recently, non-coding RNA (ncRNA) have been increasingly investigated for their contributions to the pathology of diseases and potential for diagnostic and therapeutic applications. Here we will review several ncRNA - H19, MALAT1, ANRIL, NEAT1, pseudogenes, small nucleolar RNA, piwi-interacting RNA and circular RNA - and their involvement in stroke pathology. We also examine these ncRNA as potential diagnostic biomarkers, particularly in circulating blood, and as targets for therapeutic interventions. An important aspect of this is a discussion of potential methods of treatment delivery to allow for targeting of interventions past the blood-brain barrier, including lipid nanoparticles, polymer nanoparticles, and viral and non-viral vectors. Overall, several long non-coding RNA (lncRNA) discussed here have strong implications for the development of pathology and functional recovery after ischemic stroke. LncRNAs H19 and ANRIL show potential as diagnostic biomarkers, while H19 and MALAT1 may prove to be effective therapeutics for both minimising damage as well as promoting recovery. Other ncRNA have also been implicated in ischemic stroke but are currently too poorly understood to make inferences for diagnosis or treatment. Whilst the field of ncRNAs is relatively new, significant work has already highlighted that ncRNAs represent a promising novel investigative tool for understanding stroke pathology, could be used as diagnostic biomarkers, and as targets for therapeutic interventions.

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.

Non Coding RNAs as Regulators of Wnt/β-Catenin and Hippo Pathways in Arrhythmogenic Cardiomyopathy

Arrhythmogenic cardiomyopathy (ACM) is an inherited cardiomyopathy histologically characterized by the replacement of myocardium by fibrofatty infiltration, cardiomyocyte loss, and inflammation. ACM has been defined as a desmosomal disease because most of the mutations causing the disease are located in genes encoding desmosomal proteins. Interestingly, the instable structures of these intercellular junctions in this disease are closely related to a perturbed Wnt/β-catenin pathway. Imbalance in the Wnt/β-catenin signaling and also in the crosslinked Hippo pathway leads to the transcription of proadipogenic and profibrotic genes. Aiming to shed light on the mechanisms by which Wnt/β-catenin and Hippo pathways modulate the progression of the pathological ACM phenotype, the study of non-coding RNAs (ncRNAs) has emerged as a potential source of actionable targets. ncRNAs comprise a wide range of RNA species (short, large, linear, circular) which are able to finely tune gene expression and determine the final phenotype. Some share recognition sites, thus referred to as competing endogenous RNAs (ceRNAs), and ensure a coordinating action. Recent cancer research studies regarding the key role of ceRNAs in Wnt/β-catenin and Hippo pathways modulation pave the way to better understanding the molecular mechanisms underlying ACM.

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.

Estrogen distinctly regulates transcription and translation of lncRNAs and pseudogenes in breast cancer cells

Estrogen drives key transcriptional changes in breast cancer and stimulates breast cancer cells' growth with multiple mechanisms to coordinate transcription and translation. In addition to protein-coding transcripts, estrogen can regulate long non-coding RNA (lncRNA) transcripts, plus diverse non-coding RNAs including antisense, enhancer, and intergenic. LncRNA genes comprise the majority of human genes. The accidental, or regulated, translation of their short open reading frames by ribosomes remains a controversial topic. Here we report for the first time an integrated analysis of RNA abundance and ribosome occupancy level, using Ribo-seq combined with RNA-Seq, in the estrogen-responsive, estrogen receptor α positive, human breast cancer cell model MCF7, before and after hormone treatment. Translational profiling can determine, in an unbiased manner, which fraction of the genome is actually translated into proteins, as well as resolving whether transcription and translation respond concurrently, or differentially, to estrogen treatment. Our data showed specific transcripts more robustly detected in RNA-Seq than in the ribosome-profiling data, and vice versa, suggesting distinct gene-specific estrogen responses at the transcriptional and the translational level, respectively. Here, we showed that estrogen stimulation affects the expression levels of numerous lncRNAs, but not their association with ribosomes, and that most lncRNAs are not ribosome-bound. For the first time, we also demonstrated the transcriptional and translational response of expressed pseudogenes to estrogen, pointing to new perspectives for drug-target development in breast cancer in the future.

The Role of Transposable Elements of the Human Genome in Neuronal Function and Pathology

Transposable elements (TEs) have been extensively studied for decades. In recent years, the introduction of whole-genome and whole-transcriptome approaches, as well as single-cell resolution techniques, provided a breakthrough that uncovered TE involvement in host gene expression regulation underlying multiple normal and pathological processes. Of particular interest is increased TE activity in neuronal tissue, and specifically in the hippocampus, that was repeatedly demonstrated in multiple experiments. On the other hand, numerous neuropathologies are associated with TE dysregulation. Here, we provide a comprehensive review of literature about the role of TEs in neurons published over the last three decades. The first chapter of the present review describes known mechanisms of TE interaction with host genomes in general, with the focus on mammalian and human TEs; the second chapter provides examples of TE exaptation in normal neuronal tissue, including TE involvement in neuronal differentiation and plasticity; and the last chapter lists TE-related neuropathologies. We sought to provide specific molecular mechanisms of TE involvement in neuron-specific processes whenever possible; however, in many cases, only phenomenological reports were available. This underscores the importance of further studies in this area.

mRNA Vaccines: Why Is the Biology of Retroposition Ignored?

The major advantage of mRNA vaccines over more conventional approaches is their potential for rapid development and large-scale deployment in pandemic situations. In the current COVID-19 crisis, two mRNA COVID-19 vaccines have been conditionally approved and broadly applied, while others are still in clinical trials. However, there is no previous experience with the use of mRNA vaccines on a large scale in the general population. This warrants a careful evaluation of mRNA vaccine safety properties by considering all available knowledge about mRNA molecular biology and evolution. Here, I discuss the pervasive claim that mRNA-based vaccines cannot alter genomes. Surprisingly, this notion is widely stated in the mRNA vaccine literature but never supported by referencing any primary scientific papers that would specifically address this question. This discrepancy becomes even more puzzling if one considers previous work on the molecular and evolutionary aspects of retroposition in murine and human populations that clearly documents the frequent integration of mRNA molecules into genomes, including clinical contexts. By performing basic comparisons, I show that the sequence features of mRNA vaccines meet all known requirements for retroposition using L1 elements-the most abundant autonomously active retrotransposons in the human genome. In fact, many factors associated with mRNA vaccines increase the possibility of their L1-mediated retroposition. I conclude that is unfounded to a priori assume that mRNA-based therapeutics do not impact genomes and that the route to genome integration of vaccine mRNAs via endogenous L1 retroelements is easily conceivable. This implies that we urgently need experimental studies that would rigorously test for the potential retroposition of vaccine mRNAs. At present, the insertional mutagenesis safety of mRNA-based vaccines should be considered unresolved.

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.

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.

Functions of lncRNA DUXAP8 in non-small cell lung cancer

Non-small cell lung cancer (NSCLC) poses a serious threat to public health due to its significant morbidity and mortality rates. The processes of NSCLC formation and development are quite complex and involve numerous regulatory biomolecules. Long non-coding RNAs (lncRNAs) have attracted attention since they have been found to play critical roles in the tumorigenesis of various human malignancies. Recently, double homeobox A pseudogene 8 (DUXAP8) was identified as an oncogenic lncRNA that is overexpressed in different tumor types. In NSCLC, high expression of DUXAP8 is associated with poor prognosis in patients. The regulatory mechanism underlying the oncogenic effects of DUXAP8 can be divided into transcriptional level and post-transcriptional level. DUXAP8 promotes proliferation, epithelial-mesenchymal transition, and aerobic glycolysis in NSCLC cells. Moreover, DUXAP8 shows potential for the diagnosis and treatment of NSCLC. Herein, we review the molecular mechanisms underlying the DUXAP8-mediated phenotypes of NSCLC as well as its potential clinical applications.

CRISPR/Cas9-induced gene conversion between ATAD3 paralogs

Paralogs and pseudogenes are abundant within the human genome, and can mediate non-allelic homologous recombination (NAHR) or gene conversion events. The ATAD3 locus contains three paralogs situated in tandem, and is therefore prone to NAHR-mediated deletions and duplications associated with severe neurological phenotypes. To study this locus further, we aimed to generate biallelic loss-of-function variants in ATAD3A by CRISPR/Cas9 genome editing. Unexpectedly, two of the generated clones underwent gene conversion, as evidenced by replacement of the targeted sequence of ATAD3A by a donor sequence from its paralog ATAD3B. We highlight the complexity of CRISPR/Cas9 design, end-product formation, and recombination repair mechanisms for CRISPR/Cas9 delivery as a nucleic acid molecular therapy when targeting genes that have paralogs or pseudogenes, and advocate meticulous evaluation of resultant clones in model organisms. In addition, we suggest that endogenous gene conversion may be used to repair missense variants in genes with paralogs or pseudogenes.

Between the Devil and the Deep Blue Sea: Non-Coding RNAs Associated with Transmissible Cancers in Tasmanian Devil, Domestic Dog and Bivalves

Currently there are nine known examples of transmissible cancers in nature. They have been observed in domestic dog, Tasmanian devil, and six bivalve species. These tumours can overcome host immune defences and spread to other members of the same species. Non-coding RNAs (ncRNAs) are known to play roles in tumorigenesis and immune system evasion. Despite their potential importance in transmissible cancers, there have been no studies on ncRNA function in this context to date. Here, we present possible applications of the CRISPR/Cas system to study the RNA biology of transmissible cancers. Specifically, we explore how ncRNAs may play a role in the immortality and immune evasion ability of these tumours.

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.

Pseudogene-Mediated Gene Conversion After CRISPR-Cas9 Editing Demonstrated by Partial CD33 Conversion with SIGLEC22P

Although gene editing workflows typically consider the possibility of off-target editing, pseudogene-directed homology repair has not, to our knowledge, been reported previously. Here, we employed a CRISPR-Cas9 strategy for targeted excision of exon 2 in CD33 in U937 human monocyte cell line. Candidate clonal cell lines were screened by using a clinically relevant antibody known to label the IgV domain encoded by exon 2 (P67.6, gemtuzumab). In addition to the anticipated deletion of exon 2, we also found unexpected P67.6-negative cell lines, which had apparently retained CD33 exon 2. Sequencing revealed that these lines underwent gene conversion from the nearby SIGLEC22P pseudogene during homology repair that resulted in three missense mutations relative to CD33. Ectopic expression studies confirmed that the P67.6 epitope is dependent upon these amino acids. In summation, we report that pseudogene-directed homology repair can lead to aberrant CRISPR gene editing.

Multiomics Analysis of Neuroblastoma Cells Reveals a Diversity of Malignant Transformations

Neuroblastoma (NB) is a pediatric cancer of the developing sympathetic nervous system that exhibits significant variation in the stage of differentiation and cell composition of tumors. Global loss of DNA methylation and genomic 5-hydroxymethylcytosine (5hmC) is a hallmark of human cancers. Here, we used our recently developed single-base resolution approaches, hmTOP-seq and uTOP-seq, for construction of 5hmC maps and identification of large partially methylated domains (PMDs) in different NB cell subpopulations. The 5hmC profiles revealed distinct signatures characteristic to different cell lineages and stages of malignant transformation of NB cells in a conventional and oxygen-depleted environment, which often occurs in tumors. The analysis of the cell-type-specific PMD distribution highlighted differences in global genome organization among NB cells that were ascribed to the same lineage identity by transcriptomic networks. Collectively, we demonstrated a high informativeness of the integrative epigenomic and transcriptomic research and large-scale genome structure in investigating the mechanisms that regulate cell identities and developmental stages of NB cells. Such multiomics analysis, as compared with mutational studies, open new ways for identification of novel disease-associated features which bring prognostic and therapeutic value in treating this aggressive pediatric disease.

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.

Gucy2d selectively marks inhibitory dynorphin neurons in the spinal dorsal horn but is dispensable for pain and itch sensitivity

Introduction

Inhibitory neurons in the spinal dorsal horn can be classified based on expression of neurochemical marker genes. However, these marker genes are often expressed throughout the central nervous system, which poses challenges for manipulating genetically identified spinal neurons without undesired off-target effects.

Objectives

We investigated whether Gucy2d, previously identified as a highly selective marker of dynorphin-lineage neurons in the dorsal horn, is expressed in other locations within the adult mouse spinal cord, dorsal root ganglia (DRG), or brain. In addition, we sought to molecularly characterize Gucy2d-expressing dorsal horn neurons and investigate whether the disruption of Gucy2d gene expression affects sensitivity to itch or pain.

Methods

In situ hybridization experiments assessed Gucy2d mRNA expression in the adult mouse spinal cord, DRG, and brain, and its colocalization with Pax2, Bhlhb5, and Pde2a in dorsal horn neurons. Knockout mice lacking Gucy2d expression were compared with littermate controls to assess sensitivity to chloroquine-induced itch and dry skin-mediated chronic itch, as well as heat, cold, or mechanical stimuli.

Results

Gucy2d is selectively expressed in dynorphin-lineage neurons in lamina I-III of the adult mouse spinal cord but not in the brain or DRG. Spinal Gucy2d-expressing neurons are inhibitory neurons that also express the transcription factor Bhlhb5 and the cGMP-dependent phosphodiesterase Pde2a. Gucy2d knockout mice did not exhibit altered responses to itch or pain.

Conclusions

The selective expression of Gucy2d within a subpopulation of inhibitory dorsal horn neurons may yield a means to selectively manipulate inhibitory signaling at the level of the spinal cord without effects on the brain.

Lack of CpG islands in human unitary pseudogenes and its implication

CpG islands (CGIs) are aggregation of CpG dinucleotides in the promoters of mammalian genes. These CGIs are present in almost all the housekeeping genes and some tissue-specific genes in the mammalian genome. Extensive research has been done on the prevalence and role of CGIs in protein-coding genes. However, little is known about CGIs in pseudogenes. In the current research project, we focused on CGIs in three main classes of pseudogenes e.g., duplicated pseudogenes (DPGs), processed pseudogenes (PPGs), and unitary pseudogenes (UPGs). We discovered a predominant absence of CGIs in the promoters of all three pseudogenes. We also compared the CGI profile of these pseudogenes with their parent genes and found that unitary pseudogenes (UPGs) differ from the DPGs and PPGs in the sense that in the latter, lack of CGIs is a consequential event while in UPGs, this lack of CGIs in their promoters is not a result of pseudogenization process. We also discussed the implication of the results obtained from this comparison. To our knowledge, this is the first-ever study highlighting this aspect of UPGs throwing new insights into the evolution of genome in general and especially in the context of pseudogenes.

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.

PTENP1-AS contributes to BRAF inhibitor resistance and is associated with adverse clinical outcome in stage III melanoma

BRAF inhibitors (BRAFi) selectively target oncogenic BRAFV600E/K and are effective in 80% of advanced cutaneous malignant melanoma cases carrying the V600 mutation. However, the development of drug resistance limits their clinical efficacy. Better characterization of the underlying molecular processes is needed to further improve treatments. We previously demonstrated that transcription of PTEN is negatively regulated by the PTEN pseudogene antisense RNA, PTENP1-AS, and here we investigated the impact of this transcript on clinical outcome and BRAFi resistance in melanoma. We observed that increased expression levels of PTENP1-AS in BRAFi resistant cells associated with enrichment of EZH2 and H3K27me3 at the PTEN promoter, consequently reducing the expression levels of PTEN. Further, we showed that targeting of the PTENP1-AS transcript sensitized resistant cells to BRAFi treatment and that high expression of PTENP1-AS in stage III melanoma correlated with poor survival. Collectively, the data presented here show that PTENP1-AS is a promising target for re-sensitizing cells to BRAFi and also a possible prognostic marker for clinical outcome in stage III melanoma.

Long Non-coding RNA RUNDC3A-AS1 Promotes Lung Metastasis of Thyroid Cancer via Targeting the miR-182-5p/ADAM9

Long non-coding RNAs (lncRNAs) have been identified as influential indicators in variety of malignancies. Among which, LncRNA RUNDC3A-AS1 is reported to upregulate in thyroid cancer. However, the expression pattern and the pathological function of lncRNA RUNDC3A-AS1 in thyroid cancer is unclear. In this study, we examined the expression levels of lncRNA RUNDC3A-AS1 in the thyroid cancer tissues and cell lines via RT-qPCR analysis. The effects of RUNDC3A-AS1 on thyroid cancer cell metastasis were detected by transwell chamber assay, scratch assay in vitro and lung metastasis model in vivo. The results indicated that RUNDC3A-AS1 was highly expressed in the thyroid cancer tissues and cell lines. Functionally, knockdown of RUNDC3A-AS1 could repress the migration and invasion of thyroid cancer cells in vitro, and inhibit thyroid cancer metastasis to lung in vivo. Mechanistically, RUNDC3A-AS1 served as an inhibitor of miR-182-5p in tumor tissues and cell lines. RUNDC3A-AS1 inhibited the expression of miR-182-5p to increase the expression level of ADAM9, thus further aggravating the malignancy of thyroid cancer. Therefore, the RUNDC3A-AS1/miR-182-5p/ADAM9 axis may be a potential therapeutic target for the treatment of thyroid cancer metastasis.

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.

Kalirin-RAC controls nucleokinetic migration in ADRN-type neuroblastoma

The migrational propensity of neuroblastoma is affected by cell identity, but the mechanisms behind the divergence remain unknown. Using RNAi and time-lapse imaging, we show that ADRN-type NB cells exhibit RAC1- and kalirin-dependent nucleokinetic (NUC) migration that relies on several integral components of neuronal migration. Inhibition of NUC migration by RAC1 and kalirin-GEF1 inhibitors occurs without hampering cell proliferation and ADRN identity. Using three clinically relevant expression dichotomies, we reveal that most of up-regulated mRNAs in RAC1- and kalirin-GEF1-suppressed ADRN-type NB cells are associated with low-risk characteristics. The computational analysis shows that, in a context of overall gene set poverty, the upregulomes in RAC1- and kalirin-GEF1-suppressed ADRN-type cells are a batch of AU-rich element-containing mRNAs, which suggests a link between NUC migration and mRNA stability. Gene set enrichment analysis-based search for vulnerabilities reveals prospective weak points in RAC1- and kalirin-GEF1-suppressed ADRN-type NB cells, including activities of H3K27- and DNA methyltransferases. Altogether, these data support the introduction of NUC inhibitors into cancer treatment research.

Competing Endogenous RNA Networks in Glioma

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

Next generation sequencing for HLA loci in full heritage Pima Indians of Arizona, Part II: HLA-A, -B, and -C with selected non-classical loci at 4-field resolution from whole genome sequences

While the samples and data from the Pima Indians of the Gila River Indian Community have been included in many international HLA workshops and conferences and have been the focus of numerous population reports and the source of novel alleles at the classical HLA loci, they have not been studied for the non-classical loci. In order to expand our HLA-disease association studies, we typed over 300 whole genome sequences from full Pima heritage members, controlled for first degree relationship, and employed recently developed computer algorithms to resolve HLA alleles. Both classical-HLA-A, -B, and -C- and non-classical- HLA-E, -F, -G, -J, -L, -W, -Y, -DPA2, -DPB2, -DMA, -DMB, -DOA, -DRB2, -DRB9, TAP1- loci were typed at the 4-field level of resolution. We present allele and selected haplotype frequencies, test the genotype distributions for population structure, discuss the issues that are created for tests of Hardy-Weinberg equilibrium over the four sample spaces of high resolution HLA typing, and address the implications for the evolution of non-classical pseudogenes that are no longer expressed in a phenotype subject to natural selection.

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.

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.

Methods to Study Protein-Binding to Pseudogene Transcripts

One of the most commonly described biological feature of processed pseudogenes is the ability to influence the expression of their parental coding genes. As evidenced in several studies, the high sequence similarity between these RNA pairs sets up a certain level of competition for posttranscriptional regulators, including, among others, RNA-binding proteins (RBPs). RBPs may affect, positively or negatively, the stability of bound mRNAs, so that, if an overexpressed pseudogene competes with its homologous coding gene, the downstream protein synthesis would change, with potential pathological consequences. Given these premises, a rigorous and comprehensive understanding of interactions between pseudogene-parental gene RNA pairs and RBPs could provide further insights into the biological bases of complex diseases, such as cancer, cardiovascular disease, and type 2 diabetes, identifying novel predictive and/or prognostic biomarkers.Herein, we detail easily adaptable protocols of plasmid-based molecular cloning and RNA-electrophoretic mobility shift assay (EMSA) used in our laboratory for determining the interaction between a cytoplasmatic stabilizing protein (αCP1) and the pseudogene-parental gene RNA pair HMGA1-p /HMGA1. We also offer a general overview of RNA immunoprecipitation procedures and present novel bioinformatic tools for predicting RBPs binding sites on pseudogene transcripts.

Pseudogenes: A Novel Source of Trans-Acting Antisense RNAs

Several recent studies support a functional role for pseudogenes, a copy of a parent gene that has lost protein-coding potential, which was for a long time thought to represent only "junk" DNA. Several hundreds of pseudogenes have now been reported as transcribed RNAs in a large variety of tissues and tumor types. Most studies have focused on pseudogenes expressed in sense direction, relative to their protein-coding gene counterpart, but some reports suggest that pseudogenes can be also transcribed as antisense RNAs (asRNAs). Key regulatory genes, such as PTEN and OCT4, have in fact been reported to be under the regulation of pseudogene-expressed asRNAs. Here, we review what is known about pseudogene-expressed asRNAs, we discuss the functional role that these transcripts may have in gene regulation and we summarize the techniques that are available to study them.

Competing Endogenous RNA Networks as Biomarkers in Neurodegenerative Diseases

Protein aggregation is classically considered the main cause of neuronal death in neurodegenerative diseases (NDDs). However, increasing evidence suggests that alteration of RNA metabolism is a key factor in the etiopathogenesis of these complex disorders. Non-coding RNAs are the major contributor to the human transcriptome and are particularly abundant in the central nervous system, where they have been proposed to be involved in the onset and development of NDDs. Interestingly, some ncRNAs (such as lncRNAs, circRNAs and pseudogenes) share a common functionality in their ability to regulate gene expression by modulating miRNAs in a phenomenon known as the competing endogenous RNA mechanism. Moreover, ncRNAs are found in body fluids where their presence and concentration could serve as potential non-invasive biomarkers of NDDs. In this review, we summarize the ceRNA networks described in Alzheimer's disease, Parkinson's disease, multiple sclerosis, amyotrophic lateral sclerosis and spinocerebellar ataxia type 7, and discuss their potential as biomarkers of these NDDs. Although numerous studies have been carried out, further research is needed to validate these complex interactions between RNAs and the alterations in RNA editing that could provide specific ceRNET profiles for neurodegenerative disorders, paving the way to a better understanding of these diseases.

Identification of a 15-pseudogene based prognostic signature for predicting survival and antitumor immune response in breast cancer

Pseudogenes are noncoding RNAs that have been revealed to play critical roles in oncogenesis and tumor progression. However, their functional roles have not been comprehensively clarified in breast cancer. Here, we systematically analyzed the RNA sequencing data of 13931 pseudogenes in 775 breast cancer patients from The Cancer Genome Atlas dataset, and ultimately identified 15 prognostic pseudogenes by univariate Cox proportional hazard regression. A risk score model was constructed based on the prognostic pseudogenes via LASSO analysis and dichotomized patients into low- and high-risk subgroups. Patients in the high-risk group had a significantly shorter overall survival than those in the low-risk group. The prognostic value of these 15 pseudogenes and the risk score model were further validated in the European Genome-Phenome Archive dataset. Furthermore, we performed consensus clustering of the 15 prognostic pseudogenes and found that their expression pattern was significantly associated with tumor malignancy and host antitumor immune response, in terms of infiltrating immune cell compositions, antigen presenting genes expression, cytolytic activity and T-cell exhausted markers. This study indicated that these 15 prognostic pseudogenes were significantly correlated with tumor malignancy and host antitumor immune response in breast cancer, and might serve as potential targets for immunotherapy.

FTH1 Pseudogenes in Cancer and Cell Metabolism

Ferritin, the principal intracellular iron-storage protein localized in the cytoplasm, nucleus, and mitochondria, plays a major role in iron metabolism. The encoding ferritin genes are members of a multigene family that includes some pseudogenes. Even though pseudogenes have been initially considered as relics of ancient genes or junk DNA devoid of function, their role in controlling gene expression in normal and transformed cells has recently been re-evaluated. Numerous studies have revealed that some pseudogenes compete with their parental gene for binding to the microRNAs (miRNAs), while others generate small interference RNAs (siRNAs) to decrease functional gene expression, and still others encode functional mutated proteins. Consequently, pseudogenes can be considered as actual master regulators of numerous biological processes. Here, we provide a detailed classification and description of the structural features of the ferritin pseudogenes known to date and review the recent evidence on their mutual interrelation within the complex regulatory network of the ferritin gene family.

Mutations in G Protein-Coupled Receptors: Mechanisms, Pathophysiology and Potential Therapeutic Approaches

There are approximately 800 annotated G protein-coupled receptor (GPCR) genes, making these membrane receptors members of the most abundant gene family in the human genome. Besides being involved in manifold physiologic functions and serving as important pharmacotherapeutic targets, mutations in 55 GPCR genes cause about 66 inherited monogenic diseases in humans. Alterations of nine GPCR genes are causatively involved in inherited digenic diseases. In addition to classic gain- and loss-of-function variants, other aspects, such as biased signaling, trans-signaling, ectopic expression, allele variants of GPCRs, pseudogenes, gene fusion, and gene dosage, contribute to the repertoire of GPCR dysfunctions. However, the spectrum of alterations and GPCR involvement is probably much larger because an additional 91 GPCR genes contain homozygous or hemizygous loss-of-function mutations in human individuals with currently unidentified phenotypes. This review highlights the complexity of genomic alteration of GPCR genes as well as their functional consequences and discusses derived therapeutic approaches. SIGNIFICANCE STATEMENT: With the advent of new transgenic and sequencing technologies, the number of monogenic diseases related to G protein-coupled receptor (GPCR) mutants has significantly increased, and our understanding of the functional impact of certain kinds of mutations has substantially improved. Besides the classical gain- and loss-of-function alterations, additional aspects, such as biased signaling, trans-signaling, ectopic expression, allele variants of GPCRs, uniparental disomy, pseudogenes, gene fusion, and gene dosage, need to be elaborated in light of GPCR dysfunctions and possible therapeutic strategies.