NANOGP1, a tandem duplicate of NANOG, exhibits partial functional conservation in human naïve pluripotent stem cells.

Gene duplication events can drive evolution by providing genetic material for new gene functions, and they create opportunities for diverse developmental strategies to emerge between species. To study the contribution of duplicated genes to human early development, we examined the evolution and function of NANOGP1, a tandem duplicate of the transcription factor NANOG. We found that NANOGP1 and NANOG have overlapping but distinct expression profiles, with high NANOGP1 expression restricted to early epiblast cells and naïve-state pluripotent stem cells. Sequence analysis and epitope-tagging revealed that NANOGP1 is protein coding with an intact homeobox domain. The duplication that created NANOGP1 occurred earlier in primate evolution than previously thought and has been retained only in great apes, whereas Old World monkeys have disabled the gene in different ways, including homeodomain point mutations. NANOGP1 is a strong inducer of naïve pluripotency; however, unlike NANOG, it is not required to maintain the undifferentiated status of human naïve pluripotent cells. By retaining expression, sequence and partial functional conservation with its ancestral copy, NANOGP1 exemplifies how gene duplication and subfunctionalisation can contribute to transcription factor activity in human pluripotency and development.

Terrestrial Birth and Body Size Tune UCP1 Functionality in Seals.

The molecular evolution of the mammalian heater protein UCP1 is a powerful biomarker to understand thermoregulatory strategies during species radiation into extreme climates, such as aquatic life with high thermal conductivity. While fully aquatic mammals lost UCP1, most semiaquatic seals display intact UCP1 genes, apart from large elephant seals. Here, we show that UCP1 thermogenic activity of the small-bodied harbor seal is equally potent compared to terrestrial orthologs, emphasizing its importance for neonatal survival on land. In contrast, elephant seal UCP1 does not display thermogenic activity, not even when translating a repaired or a recently highlighted truncated version. Thus, the thermogenic benefits for neonatal survival during terrestrial birth in semiaquatic pinnipeds maintained evolutionary selection pressure on UCP1 function and were only outweighed by extreme body sizes among elephant seals, fully eliminating UCP1-dependent thermogenesis.

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.

PMS2 or PMS2CL? Characterization of variants detected in the 3' of the PMS2 gene.

PMS2 germline pathogenic variants are one of the major causes for Lynch syndrome and constitutional mismatch repair deficiencies. Variant identification in the 3' region of this gene is complicated by the presence of the pseudogene PMS2CL which shares a high sequence homology with PMS2. Consequently, short-fragment screening strategies (NGS, Sanger) may fail to discriminate variant's gene localization. Using a comprehensive analysis strategy, we assessed 42 NGS-detected variants in 76 patients and found 32 localized on PMS2 while 6 on PMS2CL. Interestingly, four variants were detected in either of them in different patients. Clinical phenotype was well correlated to genotype, making it very helpful in variant assessment. Our findings emphasize the necessity of more specific complementary analyses to confirm the gene origin of each variant detected in different individuals in order to avoid variant misinterpretation. In addition, we characterized two PMS2 genomic alterations involving Alu-mediated tandem duplication and gene conversion. Those mechanisms seemed to be particularly favored in PMS2 which contribute to frequent genomic rearrangements in the 3' region of the gene.

PIN1P1 is activated by CREB1 and promotes gastric cancer progression via interacting with YBX1 and upregulating PIN1.

Long noncoding RNAs (lncRNAs) play critical roles in the carcinogenesis and progression of cancers. However, the role and mechanism of the pseudogene lncRNA PIN1P1 in gastric carcinoma remain unclear. The expression and effects of lncRNA PIN1P1 in gastric cancer were investigated. The transcriptional regulation of CREB1 on PIN1P1 was determined by ChIP and luciferase assays. The mechanistic model of PIN1P1 in gastric cancer was further explored by RNA pull-down, RIP and western blot analysis. PIN1P1 was overexpressed in gastric cancer tissues, and upregulated PIN1P1 predicted poor prognosis in patients. CREB1 was directly combined with the promoter region of PIN1P1 to promote the transcription of PIN1P1. CREB1-mediated enhanced proliferation, migration and invasion could be partially reversed by downregulation of PIN1P1. Overexpressed PIN1P1 promoted the proliferation, migration and invasion of gastric cancer cells, whereas decreased PIN1P1 showed the opposite effects. PIN1P1 directly interacted with YBX1 and promoted YBX1 protein expression, leading to upregulation of PIN1, in which E2F1 may be involved. Silencing of YBX1 during PIN1P1 overexpression could partially rescue PIN1 upregulation. PIN1, the parental gene of PIN1P1, was elevated in gastric cancer tissues, and its upregulation was correlated with poor patient outcomes. PIN1 facilitated gastric cancer cell proliferation, migration and invasion. To sum up, CREB1-activated PIN1P1 could promote gastric cancer progression through YBX1 and upregulating PIN1, suggesting that it is a potential target for gastric cancer.

Reduction of Paraoxonase Expression Followed by Inactivation across Independent Semiaquatic Mammals Suggests Stepwise Path to Pseudogenization.

Convergent adaptation to the same environment by multiple lineages frequently involves rapid evolutionary change at the same genes, implicating these genes as important for environmental adaptation. Such adaptive molecular changes may yield either change or loss of protein function; loss of function can eliminate newly deleterious proteins or reduce energy necessary for protein production. We previously found a striking case of recurrent pseudogenization of the Paraoxonase 1 (Pon1) gene among aquatic mammal lineages-Pon1 became a pseudogene with genetic lesions, such as stop codons and frameshifts, at least four times independently in aquatic and semiaquatic mammals. Here, we assess the landscape and pace of pseudogenization by studying Pon1 sequences, expression levels, and enzymatic activity across four aquatic and semiaquatic mammal lineages: pinnipeds, cetaceans, otters, and beavers. We observe in beavers and pinnipeds an unexpected reduction in expression of Pon3, a paralog with similar expression patterns but different substrate preferences. Ultimately, in all lineages with aquatic/semiaquatic members, we find that preceding any coding-level pseudogenization events in Pon1, there is a drastic decrease in expression, followed by relaxed selection, thus allowing accumulation of disrupting mutations. The recurrent loss of Pon1 function in aquatic/semiaquatic lineages is consistent with a benefit to Pon1 functional loss in aquatic environments. Accordingly, we examine diving and dietary traits across pinniped species as potential driving forces of Pon1 functional loss. We find that loss is best associated with diving activity and likely results from changes in selective pressures associated with hypoxia and hypoxia-induced inflammation.

CSPG4P12 polymorphism served as a susceptibility marker for esophageal cancer in Chinese population.

BACKGROUND: Chondroitin sulfate proteoglycan 4 pseudogene 12 (CSPG4P12) has been implicated in the pathogenesis of various cancers. This study aimed to evaluate the association of the CSPG4P12 polymorphism with esophageal squamous cell carcinoma (ESCA) risk and to explore the biological impact of CSPG4P12 expression on ESCA cell behavior. METHODS: A case-control study was conducted involving 480 ESCA patients and 480 healthy controls to assess the association between the rs8040855 polymorphism and ESCA risk. The CSPG4P12 rs8040855 genotype was identified using the TaqMan-MGB probe method. Logistic regression model was used to evaluate the association of CSPG4P12 SNP with the risk of ESCA by calculating the odds ratios (OR) and 95% confidence intervals (95%CI ). The effects of CSPG4P12 overexpression on cell proliferation, migration, and invasion were examined in ESCA cell lines. Co-expressed genes were identified via the CBioportal database, with pathway enrichment analyzed using SangerBox. The binding score of CSPG4P12 to P53 was calculated using RNA protein interaction prediction (RPISeq). Additionally, Western Blot analysis was performed to investigate the impact of CSPG4P12 overexpression on the P53/PI3K/AKT signaling pathway. RESULTS: The presence of at least one rs8040855 G allele was associated with a reduced susceptibility to ESCA compared to the CC genotype (OR = 0.51, 95%CI = 0.28-0.93, P = 0.03). Stratification analysis revealed that the CSPG4P12 rs8040855 C allele significantly decreased the risk of ESCA among younger individuals (≤ 57 years) and non-drinkers (OR = 0.31, 95%CI = 0.12-0.77, P = 0.01; OR = 0.42, 95%CI=0.20-0.87, P = 0.02, respectively). CSPG4P12 expression was found to be downregulated in ESCA tissues compared to adjacent normal tissues. Overexpression of CSPG4P12 in ESCA cells inhibited their proliferation, migration, and invasion capabilities. Furthermore, Western Blot analysis indicated that CSPG4P12 overexpression led to a reduction in PI3K and p-AKT protein expression levels. P53 silencing rescues the inhibitory effect of CSPG4P12 on p-AKT. CONCLUSION: The CSPG4P12 rs8040855 variant is associated with reduced ESCA risk and the overexpression of CSPG4P12 inhibited the migration and invasion of ESCA cells by P53/PI3K/AKT pathway. These findings suggest that CSPG4P12 may serve as a novel biomarker for ESCA susceptibility and a potential target for therapeutic intervention.

A unique case of hyperammonemia due to CA5A deficiency: Impact of coexisting gene mutations, pseudogene, and microdeletion.

Carbonic anhydrase 5A (CA5A) belongs to a family of carbonic anhydrases which are zinc metalloenzymes involved in the reversible hydration of CO2 to bicarbonate. Mutations in CA5A are very rare and known to cause Carbonic anhydrase 5A deficiency (CA5AD), an autosomal recessive inborn error of metabolism characterized clinically by acute onset of encephalopathy in infancy or early childhood. CA5A also has two very identical pseudogenes whose interference may result in compromised accuracy in targeted sequencing. We report a unique case of CA5AD caused by compound heterozygous variant (NM_001739.2: c.721G>A: p.Glu241Lys & NM_001739.2: c.619-3420_c.774 + 502del4078bp) in an infant in order to expand the phenotypic spectrum and underscore the impact of pseudogenes, which can introduce complexities in molecular genetic analysis.

Molecular basis of FAAH-OUT-associated human pain insensitivity.

Chronic pain affects millions of people worldwide and new treatments are needed urgently. One way to identify novel analgesic strategies is to understand the biological dysfunctions that lead to human inherited pain insensitivity disorders. Here we report how the recently discovered brain and dorsal root ganglia-expressed FAAH-OUT long non-coding RNA (lncRNA) gene, which was found from studying a pain-insensitive patient with reduced anxiety and fast wound healing, regulates the adjacent key endocannabinoid system gene FAAH, which encodes the anandamide-degrading fatty acid amide hydrolase enzyme. We demonstrate that the disruption in FAAH-OUT lncRNA transcription leads to DNMT1-dependent DNA methylation within the FAAH promoter. In addition, FAAH-OUT contains a conserved regulatory element, FAAH-AMP, that acts as an enhancer for FAAH expression. Furthermore, using transcriptomic analyses in patient-derived cells we have uncovered a network of genes that are dysregulated from disruption of the FAAH-FAAH-OUT axis, thus providing a coherent mechanistic basis to understand the human phenotype observed. Given that FAAH is a potential target for the treatment of pain, anxiety, depression and other neurological disorders, this new understanding of the regulatory role of the FAAH-OUT gene provides a platform for the development of future gene and small molecule therapies.

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.

Evidence That Aquaporin 11 (AQP11) in the Spiny Dogfish (Squalus acanthias) May Represent a Pseudogene.

Various attempts to amplify an AQP11 cDNA from tissues of the spiny dogfish (Squalus acanthias) were made. Two pairs of deoxy-inosine-containing degenerate primers were designed based on conserved amino acid sequences from an AQP11 alignment. These primers yielded some faint bands from gill cDNA that were sequenced. Blast searches with the sequences showed they were not AQP11. An elasmobranch AQP11 nucleotide sequence alignment was produced to identify conserved regions to make further degenerate primers. One primer pair produced a short 148 bp fragment showing particularly strong amplification in gill and intestine. It was sequenced and represented a piece of the AQP11 gene. However, as the fragment may have resulted from contaminating genomic DNA (in total RNA used to make cDNA), 5' and 3' RACE were performed to amplify the two ends of the putative cDNA. Furthermore, 5' and 3' RACE amplifications depend on the presence of a 5' cap nucleotide and a poly A tail, respectively on the putative AQP11 mRNA. Hence, successful amplification was only possible from cDNA and not genomic DNA. Nested RACE amplifications were performed using gill and intestinal RACE cDNA, but none of the DNA fragments sequenced were AQP11. Consequently, the spiny dogfish AQP11 gene may represent a pseudogene.

Paleozoic Protein Fossils Illuminate the Evolution of Vertebrate Genomes and Transposable Elements.

Genomes hold a treasure trove of protein fossils: Fragments of formerly protein-coding DNA, which mainly come from transposable elements (TEs) or host genes. These fossils reveal ancient evolution of TEs and genomes, and many fossils have been exapted to perform diverse functions important for the host's fitness. However, old and highly degraded fossils are hard to identify, standard methods (e.g. BLAST) are not optimized for this task, and few Paleozoic protein fossils have been found. Here, a recently optimized method is used to find protein fossils in vertebrate genomes. It finds Paleozoic fossils predating the amphibian/amniote divergence from most major TE categories, including virus-related Polinton and Gypsy elements. It finds 10 fossils in the human genome (eight from TEs and two from host genes) that predate the last common ancestor of all jawed vertebrates, probably from the Ordovician period. It also finds types of transposon and retrotransposon not found in human before. These fossils have extreme sequence conservation, indicating exaptation: some have evidence of gene-regulatory function, and they tend to lie nearest to developmental genes. Some ancient fossils suggest "genome tectonics," where two fragments of one TE have drifted apart by up to megabases, possibly explaining gene deserts and large introns. This paints a picture of great TE diversity in our aquatic ancestors, with patchy TE inheritance by later vertebrates, producing new genes and regulatory elements on the way. Host-gene fossils too have contributed anciently conserved DNA segments. This paves the way to further studies of ancient protein fossils.

Pseudofinder: Detection of Pseudogenes in Prokaryotic Genomes.

Prokaryotic genomes are usually densely packed with intact and functional genes. However, in certain contexts, such as after recent ecological shifts or extreme population bottlenecks, broken and nonfunctional gene fragments can quickly accumulate and form a substantial fraction of the genome. Identification of these broken genes, called pseudogenes, is a critical step for understanding the evolutionary forces acting upon, and the functional potential encoded within, prokaryotic genomes. Here, we present Pseudofinder, an open-source software dedicated to pseudogene identification and analysis in bacterial and archaeal genomes. We demonstrate that Pseudofinder's multi-pronged, reference-based approach can detect a wide variety of pseudogenes, including those that are highly degraded and typically missed by gene-calling pipelines, as well newly formed pseudogenes containing only one or a few inactivating mutations. Additionally, Pseudofinder can detect genes that lack inactivating substitutions but experiencing relaxed selection. Implementation of Pseudofinder in annotation pipelines will allow more precise estimations of the functional potential of sequenced microbes, while also generating new hypotheses related to the evolutionary dynamics of bacterial and archaeal genomes.

Pseudogene TDGF1P3 regulates the proliferation and metastasis of colorectal cancer cells via the miR-338-3p-PKM2 axis.

Research in the past decade has revealed significant roles of pseudogenes in colorectal cancer (CRC). Here, the role of teratocarcinoma-derived growth factor 1 pseudogene 3 (TDGF1P3) in regulating the proliferation and invasion of CRC cells was investigated; in addition, its downstream targets were analyzed, and the underlying mechanisms were elucidated. TDGF1P3 was determined to be upregulated in CRC cells and tissues. Silencing TDGF1P3 substantially repressed cell proliferation, migration, and invasion in vitro. Similarly, in vivo assays showed that TDGF1P3 knockdown attenuated tumor growth in nude mice. Mechanistic investigations revealed that TDGF1P3 directly bound to miR-338-3p, thereby preventing miR-338-3p from binding to its target mRNA pyruvate kinase M2 (PKM2). Functional rescue tests indicated that TDGF1P3 regulates CRC cell proliferation and invasion by restraining the miR-338-3p-PKM2 axis. Thus, these data illustrated that TDGF1P3 exerts its oncogenic activity by upregulating PKM2 via competitively binding miR-338-3p, which may be a therapeutic target for CRC.

The diagnostic and prognostic value of pseudogene SIGLEC17P in lung adenocarcinoma and a preliminary functional study.

Among malignant tumors, lung adenocarcinoma (LUAD) is the leading cause of death worldwide. This study explored the diagnostic, prognostic value, and preliminary functional verification of sialic acid binding Ig like lectin 17, pseudogene (SIGLEC17P) in LUAD. Prognostic lncRNAs for LUAD were identified by The Cancer Genome Atlas and quantitative real-time PCR (qRT-PCR) was used to detect the expression of SIGLEC17P in LUAD and paracarcinoma tissues. Subsequently, lentiviral vectors were used to overexpress SIGLEC17P in A549 and H1299 cells. The effects of SIGLEC17P overexpression on the proliferation, migration, and invasiveness of LUAD cells (A549 and H1299) were evaluated by Cell Counting Kit-8, wound healing, and transwell migration assays, respectively. Bioinformatics analyses were performed to reveal the potential pathways in which SIGLEC17P is involved in LUAD. qRT-PCR results revealed low SIGLEC17P expression in LUAD tissues and a significant association with the N stage, T stage, and tumor node metastasis stage. Furthermore, the receiver operating characteristic curve demonstrated a reliable diagnostic value. The proliferation, migration, and invasion of LUAD cells were inhibited by overexpression of SIGLEC17P. Bioinformatics analyses suggested that SIGLEC17P might exert antioncogenic effects in LUAD through the mir-20-3p/ADH1B or mir-4476-5p/DPYSL axis. In summary, our results revealed that SIGLEC17P acts as a prognostic biomarker, independent prognostic factor, and potential therapeutic target for patients with LUAD.

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.

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.

A Pathogenic Variant Reclassified to the Pseudogene PMS2P1 in a Patient with Suspected Hereditary Cancer.

The PMS2 gene is involved in DNA repair by the mismatch repair pathway. Deficiencies in this mechanism have been associated with Lynch Syndrome (LS), which is characterized by a high risk for colorectal, endometrial, ovarian, breast, and other cancers. Germinal pathogenic variants of PMS2 are associated with up to 5% of all cases of LS. The prevalence is overestimated for the existence of multiple homologous pseudogenes. We report the case of a 44-year-old woman diagnosed with breast cancer at 34 years without a relevant cancer family history. The presence of pathogenic variant NM_000535.7:c.1A > T, (p.Met1Leu) in PMS2 was determined by next-generation sequencing analysis with a panel of 322 cancer-associated genes and confirmed by capillary sequencing in the patient. The variant was determined in six family members (brothers, sisters, and a son) and seven non-cancerous unrelated individuals. Analysis of the amplified region showed high homology of PMS2 with five of its pseudogenes. We determined that the variant is associated with the PMS2P1 pseudogene following sequence alignment analysis. We propose considering the variant c.1A > T, (p.Met1Leu) in PMS2 for reclassification as not hereditary cancer-related, given the impact on the diagnosis and treatment of cancer patients and families carrying this variant.

Long noncoding RNAs: biogenesis, regulation, function, and their emerging significance in toxicology.

The repertoire of regulatory noncoding RNAs (ncRNAs) has been enriched by the inclusion of long noncoding RNA (lncRNA) that are longer than 200 nt. Some of the currently known lncRNAs, were reported in the 1990s before the term lncRNA was introduced. These lncRNAs have diverse regulatory functions including regulation of transcription via interactions with proteins and RNAs, chromatin remodeling, translation, posttranslational protein modification, protein trafficking and cell signaling. Predictably, the dysregulation of lncRNA expression due to exposure to toxicants may precipitate adverse health consequences. Dysregulation of lncRNAs has also been implicated in various adverse human health outcomes. There is an increasing agreement that lncRNA expression profiling data needs to be closely examined to determine whether their altered expression can be used as biomarkers of toxicity as well as adverse human health outcomes. This review summarizes the biogenesis, regulation, function of lncRNA and their emerging significance in toxicology and disease conditions. Because our understanding of the lncRNA-toxicity relationship is still evolving, this review discusses this developing field using some examples.

Biological role and diagnostic utility of ribosomal protein L23a pseudogene 53 in cutaneous melanoma.

Background: Skin melanoma is one of the deadliest types of skin cancer and develops from melanocytes. The genetic aberrations in protein-coding genes are well characterized, but little is known about changes in non-coding RNAs (ncRNAs) such as pseudogenes. Ribosomal protein pseudogenes (RPPs) have been described as the largest group of pseudogenes which are dispersed in the human genome. Materials and methids: We looked deeply at the role of one of them, ribosomal protein L23a pseudogene 53 (RPL23AP53), and its potential diagnostic use. The expression level of RPL23AP53 was profiled in melanoma cell lines using real time quantitative reverse transcription polymerase chain reaction (qRT-PCR) and analyzed based on the Cancer Genome Atlas (TCGA) data depending on BRAF status and clinicopathological parameters. Cellular phenotype, which was associated with RPL23AP53 levels, was described based on the REACTOME pathway browser, Gene Set Enrichment Analysis (GSEA) analysis as well as Immune and ESTIMATE Scores. Results: We indicted in vitro changes in RPL23AP53 level depending on a cell line, and based on in silico analysis of TCGA samples demonstrated significant differences in RPL23AP53 expression between primary and metastatic melanoma, as well as correlation between RPL23AP53 and overall survival. No differences depending on BRAF status were observed. RPL23AP53 is associated with several signaling pathways and cellular processes. Conclusions: This study showed that patients with higher expression of RPL23AP53 displayed changed infiltration of lymphocytes, macrophages, and neutrophils compared to groups with lower expression of RPL23AP53. RPL23AP53 pseudogene is differently expressed in melanoma compared with normal tissue and its expression is associated with cellular proliferation. Thus, it may be considered as an indicator of patients' survival and a marker for the immune profile assessment.