Processed pseudogenes acquired somatically during cancer development

MGAT4EP promotes tumor progression and serves as a prognostic marker for breast cancer

Breast cancer remains a global health challenge with varied prognoses despite treatment advancements. Therefore, this study explores the pseudogene MGAT4EP as a potential biomarker and therapeutic target in breast cancer. Using TCGA data and bioinformatics, MGAT4EP was identified as significantly overexpressed in breast cancer tissues and associated with poor prognosis. Multivariate Cox regression confirmed MGAT4EP as important prognostic factor. A clinical prediction model based on MGAT4EP expression showed high accuracy for 1-, 3-, and 5-year survival rates and was translated into a nomogram for clinical application. Functional studies revealed that silencing MGAT4EP via siRNA promoted apoptosis, inhibited migration and invasion in breast cancer cells. RNA-seq, GSEA, and GO analyses linked MGAT4EP to apoptosis and focal adhesion pathways. Notably, knock down of MGAT4EP significantly suppressed tumor growth and metastasis in xenograft and lung metastasis models. Taken together, these findings establish MGAT4EP as an attractive target for metastatic breast cancer and provide a potential a promising therapeutic target for breast cancer treatment.

Structural features of somatic and germline retrotransposition events in humans

BACKGROUND

Transposons are DNA sequences able to move or copy themselves to other genomic locations leading to insertional mutagenesis. Although transposon-derived sequences account for half of the human genome, most elements are no longer transposition competent. Moreover, transposons are normally repressed through epigenetic silencing in healthy adult tissues but become derepressed in several human cancers, with high activity detected in colorectal cancer. Their impact on non-malignant and malignant tissue as well as the differences between somatic and germline retrotransposition remain poorly understood. With new sequencing technologies, including long read sequencing, we can access intricacies of retrotransposition, such as insertion sequence details and nested repeats, that have been previously challenging to characterize.

RESULTS

In this study, we investigate somatic and germline retrotransposition by analyzing long read sequencing from 56 colorectal cancers and 112 uterine leiomyomas. We identified 1495 somatic insertions in colorectal samples, while striking lack of insertions was detected in uterine leiomyomas. Our findings highlight differences between somatic and germline events, such as transposon type distribution, insertion length, and target site preference. Leveraging long-read sequencing, we provide an in-depth analysis of the twin-priming phenomenon, detecting it across transposable element types that remain active in humans, including Alus. Additionally, we detect an abundance of germline transposons in repetitive DNA, along with a relationship between replication timing and insertion target site.

CONCLUSIONS

Our study reveals a stark contrast in somatic transposon activity between colorectal cancers and uterine leiomyomas, and highlights differences between somatic and germline transposition. This suggests potentially different conditions in malignant and non-malignant tissues, as well as in germline and somatic tissues, which could be involved in the transposition process. Long-read sequencing provided important insights into transposon behavior, allowing detailed examination of structural features such as twin priming and nested elements.

Maintenance of long-term transposable element activity through regulation by nonautonomous elements

Transposable elements are DNA sequences that can move and replicate within genomes. Broadly, there are 2 types: autonomous elements, which encode the necessary enzymes for transposition, and nonautonomous elements, which rely on the enzymes produced by autonomous elements for their transposition. Nonautonomous elements have been proposed to regulate the numbers of transposable elements, which is a possible explanation for the persistence of transposition activity over long evolutionary times. However, previous modeling studies indicate that interactions between autonomous and nonautonomous elements usually result in the extinction of one type. Here, we study a stochastic model that allows for the stable coexistence of autonomous and nonautonomous retrotransposons. We determine the conditions for this coexistence and derive an analytical expression for the stationary distribution of their copy numbers, showing that nonautonomous elements regulate stochastic fluctuations and the number of autonomous elements in stationarity. We find that the stationary variances of each element can be expressed as a function of the average copy numbers and their covariance, enabling data comparison and model validation. These results suggest that continued transposition activity of transposable elements, regulated by nonautonomous elements, is a possible evolutionary outcome that could for example explain the long coevolutionary history of autonomous LINE1 and nonautonomous Alu element transposition in the human ancestry.

Competing endogenous RNAs network dysregulation in oral cancer: a multifaceted perspective on crosstalk and competition

Oral cancer progresses from asymptomatic to advanced stages, often involving cervical lymph node metastasis, resistance to chemotherapy, and an unfavorable prognosis. Clarifying its potential mechanisms is vital for developing effective theraputic strategies. Recent research suggests a substantial involvement of non-coding RNA (ncRNA) in the initiation and advancement of oral cancer. However, the underlying roles and functions of various ncRNA types in the growth of this malignant tumor remain unclear. Competing endogenous RNAs (ceRNAs) refer to transcripts that can mutually regulate each other at the post-transcriptional level by vying for shared miRNAs. Networks of ceRNAs establish connections between the functions of protein-coding mRNAs and non-coding RNAs, including microRNA, long non-coding RNA, pseudogenic RNA, and circular RNA, piwi-RNA, snoRNA. A growing body of research has indicated that imbalances in ceRNAs networks play a crucial role in various facets of oral cancer, including development, metastasis, migration, invasion, and inflammatory responses. Hence, delving into the regulatory pathways of ceRNAs in oral cancer holds the potential to advance our understanding of the pathological mechanisms, facilitate early diagnosis, and foster targeted drug development for this malignancy. The present review summarized the fundamental role of ceRNA network, discussed the limitations of current ceRNA applications, which have been improved through chemical modification and carrier delivery as new biomarkers for diagnosis and prognosis is expected to offer a groundbreaking therapeutic approach for individuals with oral cancer.

Chromosomal rearrangements and instability caused by the LINE-1 retrotransposon

LINE-1 (L1) retrotransposition is widespread in many cancers, especially those with a high burden of chromosomal rearrangements. However, whether and to what degree L1 activity directly impacts genome integrity is unclear. Here, we apply whole-genome sequencing to experimental models of L1 expression to comprehensively define the spectrum of genomic changes caused by L1. We provide definitive evidence that L1 expression frequently and directly causes both local and long-range chromosomal rearrangements, small and large segmental copy-number alterations, and subclonal copy-number heterogeneity due to ongoing chromosomal instability. Mechanistically, all these alterations arise from DNA double-strand breaks (DSBs) generated by L1-encoded ORF2p. The processing of ORF2p-generated DSB ends prior to their ligation can produce diverse rearrangements of the target sequences. Ligation between DSB ends generated at distal loci can generate either stable chromosomes or unstable dicentric, acentric, or ring chromosomes that undergo subsequent evolution through breakage-fusion bridge cycles or DNA fragmentation. Together, these findings suggest L1 is a potent mutagenic force capable of driving genome evolution beyond simple insertions.

The reconstruction of evolutionary dynamics of processed pseudogenes indicates deep silencing of "retrobiome" in naked mole rat

Approximately half of mammalian genomes are occupied by retrotransposons, highly repetitive interspersed genetic elements expanded through the mechanism of reverse transcription. The evolution of this "retrobiome" involved a series of explosive amplifications, presumably associated with high mutation rates, interspersed with periods of silencing. A by-product of retrotransposon activity is the formation of processed pseudogenes (PPGs)-intron-less, promoter-less DNA copies of messenger RNA (mRNA). We examined the proportion of PPGs with varying degrees of deviation from their ancestor mRNAs as an indicator of the intensity of retrotranspositions at different times in the past. Our analysis revealed a high proportion of "young'' (recently acquired) PPGs in the DNA of mice and rats, indicating significant retrobiome activity during the recent evolution of these species. The ongoing process of new PPG entries in mouse germ line DNA was confirmed by identifying diversity in PPG content within the single strain of mice, C57BL/6. In contrast, the highly abundant PPGs of the naked mole rat (NMR) exhibited substantial deviation from their mRNAs, with a near-complete lack of PPGs without mutations, indicative of the silencing of the retrobiome in the most recent evolutionary past, preceded by a period of high activity. This distinctive feature of the NMR genome was confirmed through the analysis of a broad range of mammalian species. The peculiar evolutionary dynamics of PPGs in the NMR, an organism with exceptional longevity and resistance to cancer, may reflect the role played by the retrobiome in aging and cancer.

Coding, or non-coding, that is the question

The advent of high-throughput sequencing uncovered that our genome is pervasively transcribed into RNAs that are seemingly not translated into proteins. It was also found that non-coding RNA transcripts outnumber canonical protein-coding genes. This mindboggling discovery prompted a surge in non-coding RNA research that started unraveling the functional relevance of these new genetic units, shaking the classic definition of "gene". While the non-coding RNA revolution was still taking place, polysome/ribosome profiling and mass spectrometry analyses revealed that peptides can be translated from non-canonical open reading frames. Therefore, it is becoming evident that the coding vs non-coding dichotomy is way blurrier than anticipated. In this review, we focus on several examples in which the binary classification of coding vs non-coding genes is outdated, since the same bifunctional gene expresses both coding and non-coding products. We discuss the implications of this intricate usage of transcripts in terms of molecular mechanisms of gene expression and biological outputs, which are often concordant, but can also surprisingly be discordant. Finally, we discuss the methodological caveats that are associated with the study of bifunctional genes, and we highlight the opportunities and challenges of therapeutic exploitation of this intricacy towards the development of anticancer therapies.

Transposon DNA sequences facilitate the tissue-specific gene transfer of circulating tumor DNA between human cells

The exchange of genes between cells is known to play an important physiological and pathological role in many organisms. We show that circulating tumor DNA (ctDNA) facilitates cell-specific gene transfer between human cancer cells and explain part of the mechanisms behind this phenomenon. As ctDNA migrates into the nucleus, genetic information is transferred. Cell targeting and ctDNA integration require ERVL, SINE or LINE DNA sequences. Chemically manufactured AluSp and MER11C sequences replicated multiple myeloma (MM) ctDNA cell targeting and integration. Additionally, we found that ctDNA may alter the treatment response of MM and pancreatic cancer models. This study shows that retrotransposon DNA sequences promote cancer gene transfer. However, because cell-free DNA has been detected in physiological and other pathological conditions, our findings have a broader impact than just cancer. Furthermore, the discovery that transposon DNA sequences mediate tissue-specific targeting will open up a new avenue for the delivery of genes and therapies.

LINE-1 retrotransposition and its deregulation in cancers: implications for therapeutic opportunities

Long interspersed element 1 (LINE-1) is the only protein-coding transposon that is active in humans. LINE-1 propagates in the genome using RNA intermediates via retrotransposition. This activity has resulted in LINE-1 sequences occupying approximately one-fifth of our genome. Although most copies of LINE-1 are immobile, ∼100 copies are retrotransposition-competent. Retrotransposition is normally limited via epigenetic silencing, DNA repair, and other host defense mechanisms. In contrast, LINE-1 overexpression and retrotransposition are hallmarks of cancers. Here, we review mechanisms of LINE-1 regulation and how LINE-1 may promote genetic heterogeneity in tumors. Finally, we discuss therapeutic strategies to exploit LINE-1 biology in cancers.

Precise characterization of somatic complex structural variations from tumor/control paired long-read sequencing data with nanomonsv

We present our novel software, nanomonsv, for detecting somatic structural variations (SVs) using tumor and matched control long-read sequencing data with a single-base resolution. The current version of nanomonsv includes two detection modules, Canonical SV module, and Single breakend SV module. Using tumor/control paired long-read sequencing data from three cancer and their matched lymphoblastoid lines, we demonstrate that Canonical SV module can identify somatic SVs that can be captured by short-read technologies with higher precision and recall than existing methods. In addition, we have developed a workflow to classify mobile element insertions while elucidating their in-depth properties, such as 5' truncations, internal inversions, as well as source sites for 3' transductions. Furthermore, Single breakend SV module enables the detection of complex SVs that can only be identified by long-reads, such as SVs involving highly-repetitive centromeric sequences, and LINE1- and virus-mediated rearrangements. In summary, our approaches applied to cancer long-read sequencing data can reveal various features of somatic SVs and will lead to a better understanding of mutational processes and functional consequences of somatic SVs.

Rilotumumab Resistance Acquired by Intracrine Hepatocyte Growth Factor Signaling

Drug resistance is a long-standing impediment to effective systemic cancer therapy and acquired drug resistance is a growing problem for molecularly-targeted therapeutics that otherwise have shown unprecedented successes in disease control. The hepatocyte growth factor (HGF)/Met receptor pathway signaling is frequently involved in cancer and has been a subject of targeted drug development for nearly 30 years. To anticipate and study specific resistance mechanisms associated with targeting this pathway, we engineered resistance to the HGF-neutralizing antibody rilotumumab in glioblastoma cells harboring autocrine HGF/Met signaling, a frequent abnormality of this brain cancer in humans. We found that rilotumumab resistance was acquired through an unusual mechanism comprising dramatic HGF overproduction and misfolding, endoplasmic reticulum (ER) stress-response signaling and redirected vesicular trafficking that effectively sequestered rilotumumab and misfolded HGF from native HGF and activated Met. Amplification of MET and HGF genes, with evidence of rapidly acquired intron-less, reverse-transcribed copies in DNA, was also observed. These changes enabled persistent Met pathway activation and improved cell survival under stress conditions. Point mutations in the HGF pathway or other complementary or downstream growth regulatory cascades that are frequently associated with targeted drug resistance in other prevalent cancer types were not observed. Although resistant cells were significantly more malignant, they retained sensitivity to Met kinase inhibition and acquired sensitivity to inhibition of ER stress signaling and cholesterol biosynthesis. Defining this mechanism reveals details of a rapidly acquired yet highly-orchestrated multisystem route of resistance to a selective molecularly-targeted agent and suggests strategies for early detection and effective intervention.

svaRetro and svaNUMT: modular packages for annotating retrotransposed transcripts and nuclear integration of mitochondrial DNA in genome sequencing data

Nuclear integration of mitochondrial genomes and retrocopied transcript insertion are biologically important but often-overlooked aspects of structural variant (SV) annotation. While tools for their detection exist, these typically rely on reanalysis of primary data using specialised detectors rather than leveraging calls from general purpose structural variant callers. Such reanalysis potentially leads to additional computational expense and does not take advantage of advances in general purpose structural variant calling. Here, we present svaRetro and svaNUMT; R packages that provide functions for annotating novel genomic events, such as nonreference retrocopied transcripts and nuclear integration of mitochondrial DNA. The packages were developed to work within the Bioconductor framework. We evaluate the performance of these packages to detect events using simulations and public benchmarking datasets, and annotate processed transcripts in a public structural variant database. svaRetro and svaNUMT provide modular, SV-caller agnostic tools for downstream annotation of structural variant calls.

Phenotypic mutations contribute to protein diversity and shape protein evolution

Errors in DNA replication generate genetic mutations, while errors in transcription and translation lead to phenotypic mutations. Phenotypic mutations are orders of magnitude more frequent than genetic ones, yet they are less understood. Here, we review the types of phenotypic mutations, their quantifications, and their role in protein evolution and disease. The diversity generated by phenotypic mutation can facilitate adaptive evolution. Indeed, phenotypic mutations, such as ribosomal frameshift and stop codon readthrough, sometimes serve to regulate protein expression and function. Phenotypic mutations have often been linked to fitness decrease and diseases. Thus, understanding the protein heterogeneity and phenotypic diversity caused by phenotypic mutations will advance our understanding of protein evolution and have implications on human health and diseases.

Pseudogenes and Liquid Phase Separation in Epigenetic Expression

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

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

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

Structural Variation in Cancer: Role, Prevalence, and Mechanisms

Somatic rearrangements resulting in genomic structural variation drive malignant phenotypes by altering the expression or function of cancer genes. Pan-cancer studies have revealed that structural variants (SVs) are the predominant class of driver mutation in most cancer types, but because they are difficult to discover, they remain understudied when compared with point mutations. This review provides an overview of the current knowledge of somatic SVs, discussing their primary roles, prevalence in different contexts, and mutational mechanisms. SVs arise throughout the life history of cancer, and 55% of driver mutations uncovered by the Pan-Cancer Analysis of Whole Genomes project represent SVs. Leveraging the convergence of cell biology and genomics, we propose a mechanistic classification of somatic SVs, from simple to highly complex DNA rearrangement classes. The actions of DNA repair and DNA replication processes together with mitotic errors result in a rich spectrum of SV formation processes, with cascading effects mediating extensive structural diversity after an initiating DNA lesion has formed. Thanks to new sequencing technologies, including the sequencing of single-cell genomes, open questions about the molecular triggers and the biomolecules involved in SV formation as well as their mutational rates can now be addressed.

Expected final online publication date for the Annual Review of Genomics and Human Genetics, Volume 23 is October 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

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.

Normal and Neoplastic Growth Suppression by the Extended Myc Network

Among the first discovered and most prominent cellular oncogenes is MYC, which encodes a bHLH-ZIP transcription factor (Myc) that both activates and suppresses numerous genes involved in proliferation, energy production, metabolism and translation. Myc belongs to a small group of bHLH-ZIP transcriptional regulators (the Myc Network) that includes its obligate heterodimerization partner Max and six "Mxd proteins" (Mxd1-4, Mnt and Mga), each of which heterodimerizes with Max and largely opposes Myc's functions. More recently, a second group of bHLH-ZIP proteins (the Mlx Network) has emerged that bears many parallels with the Myc Network. It is comprised of the Myc-like factors ChREBP and MondoA, which, in association with the Max-like member Mlx, regulate smaller and more functionally restricted repertoires of target genes, some of which are shared with Myc. Opposing ChREBP and MondoA are heterodimers comprised of Mlx and Mxd1, Mxd4 and Mnt, which also structurally and operationally link the two Networks. We discuss here the functions of these "Extended Myc Network" members, with particular emphasis on their roles in suppressing normal and neoplastic growth. These roles are complex due to the temporal- and tissue-restricted expression of Extended Myc Network proteins in normal cells, their regulation of both common and unique target genes and, in some cases, their functional redundancy.

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.

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

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

mRNA vaccines: Why and how they should be modified

The COVID-19 pandemic has stimulated the production of different therapeutic approaches for the resolution of coronavirus infections. On one hand, nanobiomolecules have been proposed as bait material for viruses,1,2 on the other hand unconventional messenger RNA vaccines have been produced like SARS-CoV-2 mRNA vaccines (BioNTech/Pfizer BNT162b2 and Moderna mRNA-1273). [...]

Noncoding RNAs Serve as the Deadliest Universal Regulators of all Cancers

Numerous cancer drivers have been identified, but they are specific to a given cancer type and condition; universal cancer drivers and universal cancer mechanisms still remain largely unclear. Here, we identified the deadliest universal drivers for all cancers via developing algorithms to analyze massive RNAseqs and clinical data from The Cancer Genome Atlas (TCGA). In general, noncoding RNAs primarily serve as the most important inducers and suppressors for all types of cancers. In particular, pseudogenes are primary inducers, and specifically the antisense RNA RP11-335K5.2 serves as the most universal cancerous driver, independently of the cancer type and condition. Therefore, noncoding RNAs, instead of proteins as conventionally thought, primarily drive cancer, which establishes a novel field for future cancer research and therapy.

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.

Genomic Mosaicism Formed by Somatic Variation in the Aging and Diseased Brain

Over the past 20 years, analyses of single brain cell genomes have revealed that the brain is composed of cells with myriad distinct genomes: the brain is a genomic mosaic, generated by a host of DNA sequence-altering processes that occur somatically and do not affect the germline. As such, these sequence changes are not heritable. Some processes appear to occur during neurogenesis, when cells are mitotic, whereas others may also function in post-mitotic cells. Here, we review multiple forms of DNA sequence alterations that have now been documented: aneuploidies and aneusomies, smaller copy number variations (CNVs), somatic repeat expansions, retrotransposons, genomic cDNAs (gencDNAs) associated with somatic gene recombination (SGR), and single nucleotide variations (SNVs). A catch-all term of DNA content variation (DCV) has also been used to describe the overall phenomenon, which can include multiple forms within a single cell's genome. A requisite step in the analyses of genomic mosaicism is ongoing technology development, which is also discussed. Genomic mosaicism alters one of the most stable biological molecules, DNA, which may have many repercussions, ranging from normal functions including effects of aging, to creating dysfunction that occurs in neurodegenerative and other brain diseases, most of which show sporadic presentation, unlinked to causal, heritable genes.

Loss of MGA repression mediated by an atypical polycomb complex promotes tumor progression and invasiveness

MGA, a transcription factor and member of the MYC network, is mutated or deleted in a broad spectrum of malignancies. As a critical test of a tumor suppressive role, we inactivated Mga in two mouse models of non-small cell lung cancer using a CRISPR-based approach. MGA loss significantly accelerated tumor growth in both models and led to de-repression of non-canonical Polycomb ncPRC1.6 targets, including genes involved in metastasis and meiosis. Moreover, MGA deletion in human lung adenocarcinoma lines augmented invasive capabilities. We further show that MGA-MAX, E2F6, and L3MBTL2 co-occupy thousands of promoters and that MGA stabilizes these ncPRC1.6 subunits. Lastly, we report that MGA loss also induces a pro-growth effect in human colon organoids. Our studies establish MGA as a bona fide tumor suppressor in vivo and suggest a tumor suppressive mechanism in adenocarcinomas resulting from widespread transcriptional attenuation of MYC and E2F target genes mediated by MGA-MAX associated with a non-canonical Polycomb complex.

Deletion of a pseudogene within a fragile site triggers the oncogenic expression of the mitotic CCSER1 gene

The oncogenic role of common fragile sites (CFS), focal and pervasive gaps in the cancer genome arising from replicative stress, remains controversial. Exploiting the TCGA dataset, we found that in most CFS the genes residing within the associated focal deletions are down-regulated, including proteins involved in tumour immune recognition. In a subset of CFS, however, the residing genes are surprisingly overexpressed. Within the most frequent CFS in this group, FRA4F, which is deleted in up to 18% of cancer cases and harbours the CCSER1 gene, we identified a region which includes an intronic, antisense pseudogene, TMSB4XP8. TMSB4XP8 focal ablation or transcriptional silencing elicits the overexpression of CCSER1, through a cis-acting mechanism. CCSER1 overexpression increases proliferation and triggers centrosome amplifications, multinuclearity, and aberrant mitoses. Accordingly, FRA4F is associated in patient samples to mitotic genes deregulation and genomic instability. As a result, cells overexpressing CCSER1 become sensitive to the treatment with aurora kinase inhibitors. Our findings point to a novel tumourigenic mechanism where focal deletions increase the expression of a new class of "dormant" oncogenes.

Giant cells: Linking McClintock's heredity to early embryogenesis and tumor origin throughout millennia of evolution on Earth

The "life code" theory postulates that egg cells, which are giant, are the first cells in reproduction and that damaged or aged giant somatic cells are the first cells in tumorigenesis. However, the hereditary basis for giant cells remains undefined. Here I propose that stress-induced genomic reorganization proposed by Nobel Laureate Barbara McClintock may represent the underlying heredity for giant cells, referred to as McClintock's heredity. Increase in cell size may serve as a response to environmental stress via switching proliferative mitosis to intranuclear replication for reproduction. Intranuclear replication activates McClintock's heredity to reset the genome following fertilization for reproduction or restructures the somatic genome for neoplastic transformation via formation of polyploid giant cancer cells (PGCCs). The genome-based McClintock heredity functions together with gene-based Mendel's heredity to regulate the genomic stability at two different stages of life cycle or tumorigenesis. Thus, giant cells link McClintock's heredity to both early embryogenesis and tumor origin. Cycling change in cell size together with ploidy number switch may represent the most fundamental mechanism on how both germ and soma for coping with environmental stresses for the survival across the tree of life which evolved over millions of years on Earth.

How Chaotic Is Genome Chaos?

Simple Summary

Cancer genomes can undergo major restructurings involving many chromosomal locations at key stages in tumor development. This restructuring process has been designated “genome chaos” by some authors. In order to examine how chaotic cancer genome restructuring may be, the cell and molecular processes for DNA restructuring are reviewed. Examination of the action of these processes in various cancers reveals a degree of specificity that indicates genome restructuring may be sufficiently reproducible to enable possible therapies that interrupt tumor progression to more lethal forms.

Abstract

Cancer genomes evolve in a punctuated manner during tumor evolution. Abrupt genome restructuring at key steps in this evolution has been called “genome chaos.” To answer whether widespread genome change is truly chaotic, this review (i) summarizes the limited number of cell and molecular systems that execute genome restructuring, (ii) describes the characteristic signatures of DNA changes that result from activity of those systems, and (iii) examines two cases where genome restructuring is determined to a significant degree by cell type or viral infection. The conclusion is that many restructured cancer genomes display sufficiently unchaotic signatures to identify the cellular systems responsible for major oncogenic transitions, thereby identifying possible targets for therapies to inhibit tumor progression to greater aggressiveness.

Higher Rates of Processed Pseudogene Acquisition in Humans and Three Great Apes Revealed by Long-Read Assemblies

LINE-1-mediated retrotransposition of protein-coding mRNAs is an active process in modern humans for both germline and somatic genomes. Prior works that surveyed human data mostly relied on detecting discordant mappings of paired-end short reads, or exon junctions contained in short reads. Moreover, there have been few genome-wide comparisons between gene retrocopies in great apes and humans. In this study, we introduced a more sensitive and accurate method to identify processed pseudogenes. Our method utilizes long-read assemblies, and more importantly, is able to provide full-length retrocopy sequences as well as flanking regions which are missed by short-read based methods. From 22 human individuals, we pinpointed 40 processed pseudogenes that are not present in the human reference genome GRCh38 and identified 17 pseudogenes that are in GRCh38 but absent from some input individuals. This represents a significantly higher discovery rate than previous reports (39 pseudogenes not in the reference genome out of 939 individuals). We also provided an overview of lineage-specific retrocopies in chimpanzee, gorilla, and orangutan genomes.

Competing endogenous RNAs in lung cancer

Competing endogenous RNAs (ceRNAs) containing microRNA response elements can competitively interact with microRNA via miRNA response elements, which can combine non-coding RNAs with protein-coding RNAs through complex ceRNA networks. CeRNAs include non-coding RNAs (long non-coding RNAs, circular RNAs, and transcribed pseudogenes) and protein-coding RNAs (mRNAs). Molecular interactions in ceRNA networks can coordinate many biological processes; however, they may also lead to ceRNA network imbalance and thus contribute to cancer occurrence when disturbed. Recent studies indicate that many dysregulated RNAs derived from lung cancer may function as ceRNAs to regulate multitudinous biological functions for lung cancer, including tumor cell proliferation, apoptosis, growth, invasion, migration, and metastasis. This study therefore reviewed the research progress in the field of non-coding and protein-coding RNAs as ceRNAs in lung cancer, and highlighted validated ceRNAs involved in biological lung cancer functions. Furthermore, the roles of ceRNAs as novel prognostic and diagnostic biomarkers were also discussed. Interpreting the involvement of ceRNAs networks in lung cancer will provide new insight into cancer pathogenesis and treatment strategies.

Pseudogenes in Cardiovascular Disease

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

Cancer, Retrogenes, and Evolution

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

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

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

Pseudogenes as Biomarkers and Therapeutic Targets in Human Cancers

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

Brain cell somatic gene recombination and its phylogenetic foundations

A new form of somatic gene recombination (SGR) has been identified in the human brain that affects the Alzheimer's disease gene, amyloid precursor protein (APP). SGR occurs when a gene sequence is cut and recombined within a single cell's genomic DNA, generally independent of DNA replication and the cell cycle. The newly identified brain SGR produces genomic complementary DNAs (gencDNAs) lacking introns, which integrate into locations distinct from germline loci. This brief review will present an overview of likely related recombination mechanisms and genomic cDNA-like sequences that implicate evolutionary origins for brain SGR. Similarities and differences exist between brain SGR and VDJ recombination in the immune system, the first identified SGR form that now has a well-defined enzymatic machinery. Both require gene transcription, but brain SGR uses an RNA intermediate and reverse transcriptase (RT) activity, which are characteristics shared with endogenous retrotransposons. The identified gencDNAs have similarities to other cDNA-like sequences existing throughout phylogeny, including intron-less genes and inactive germline processed pseudogenes, with likely overlapping biosynthetic processes. gencDNAs arise somatically in an individual to produce multiple copies; can be functional; appear most frequently within postmitotic cells; have diverse sequences; change with age; and can change with disease state. Normally occurring brain SGR may represent a mechanism for gene optimization and long-term cellular memory, whereas its dysregulation could underlie multiple brain disorders and, potentially, other diseases like cancer. The involvement of RT activity implicates already Food and Drug Administration-approved RT inhibitors as possible near-term interventions for managing SGR-associated diseases and suggest next-generation therapeutics targeting SGR elements.

PseudoChecker: an integrated online platform for gene inactivation inference

The rapid expansion of high-quality genome assemblies, exemplified by ongoing initiatives such as the Genome-10K and i5k, demands novel automated methods to approach comparative genomics. Of these, the study of inactivating mutations in the coding region of genes, or pseudogenization, as a source of evolutionary novelty is mostly overlooked. Thus, to address such evolutionary/genomic events, a systematic, accurate and computationally automated approach is required. Here, we present PseudoChecker, the first integrated online platform for gene inactivation inference. Unlike the few existing methods, our comparative genomics-based approach displays full automation, a built-in graphical user interface and a novel index, PseudoIndex, for an empirical evaluation of the gene coding status. As a multi-platform online service, PseudoChecker simplifies access and usability, allowing a fast identification of disruptive mutations. An analysis of 30 genes previously reported to be eroded in mammals, and 30 viable genes from the same lineages, demonstrated that PseudoChecker was able to correctly infer 97% of loss events and 95% of functional genes, confirming its reliability. PseudoChecker is freely available, without login required, at http://pseudochecker.ciimar.up.pt.

Our Conflict with Transposable Elements and Its Implications for Human Disease

Our genome is a historic record of successive invasions of mobile genetic elements. Like other eukaryotes, we have evolved mechanisms to limit their propagation and minimize the functional impact of new insertions. Although these mechanisms are vitally important, they are imperfect, and a handful of retroelement families remain active in modern humans. This review introduces the intrinsic functions of transposons, the tactics employed in their restraint, and the relevance of this conflict to human pathology. The most straightforward examples of disease-causing transposable elements are germline insertions that disrupt a gene and result in a monogenic disease allele. More enigmatic are the abnormal patterns of transposable element expression in disease states. Changes in transposon regulation and cellular responses to their expression have implicated these sequences in diseases as diverse as cancer, autoimmunity, and neurodegeneration. Distinguishing their epiphenomenal from their pathogenic effects may provide wholly new perspectives on our understanding of disease.

Pan-cancer analysis of whole genomes identifies driver rearrangements promoted by LINE-1 retrotransposition

About half of all cancers have somatic integrations of retrotransposons. Here, to characterize their role in oncogenesis, we analyzed the patterns and mechanisms of somatic retrotransposition in 2,954 cancer genomes from 38 histological cancer subtypes within the framework of the Pan-Cancer Analysis of Whole Genomes (PCAWG) project. We identified 19,166 somatically acquired retrotransposition events, which affected 35% of samples and spanned a range of event types. Long interspersed nuclear element (LINE-1; L1 hereafter) insertions emerged as the first most frequent type of somatic structural variation in esophageal adenocarcinoma, and the second most frequent in head-and-neck and colorectal cancers. Aberrant L1 integrations can delete megabase-scale regions of a chromosome, which sometimes leads to the removal of tumor-suppressor genes, and can induce complex translocations and large-scale duplications. Somatic retrotranspositions can also initiate breakage-fusion-bridge cycles, leading to high-level amplification of oncogenes. These observations illuminate a relevant role of L1 retrotransposition in remodeling the cancer genome, with potential implications for the development of human tumors.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

Re-recognition of pseudogenes: From molecular to clinical applications

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

PseudoFuN: Deriving functional potentials of pseudogenes from integrative relationships with genes and microRNAs across 32 cancers

Background

Long thought “relics” of evolution, not until recently have pseudogenes been of medical interest regarding regulation in cancer. Often, these regulatory roles are a direct by-product of their close sequence homology to protein-coding genes. Novel pseudogene-gene (PGG) functional associations can be identified through the integration of biomedical data, such as sequence homology, functional pathways, gene expression, pseudogene expression, and microRNA expression. However, not all of the information has been integrated, and almost all previous pseudogene studies relied on 1:1 pseudogene–parent gene relationships without leveraging other homologous genes/pseudogenes.

Results

We produce PGG families that expand beyond the current 1:1 paradigm. First, we construct expansive PGG databases by (i) CUDAlign graphics processing unit (GPU) accelerated local alignment of all pseudogenes to gene families (totaling 1.6 billion individual local alignments and >40,000 GPU hours) and (ii) BLAST-based assignment of pseudogenes to gene families. Second, we create an open-source web application (PseudoFuN [Pseudogene Functional Networks]) to search for integrative functional relationships of sequence homology, microRNA expression, gene expression, pseudogene expression, and gene ontology. We produce four “flavors” of CUDAlign-based databases (>462,000,000 PGG pairwise alignments and 133,770 PGG families) that can be queried and downloaded using PseudoFuN. These databases are consistent with previous 1:1 PGG annotation and also are much more powerful including millions of de novo PGG associations. For example, we find multiple known (e.g., miR-20a-PTEN-PTENP1) and novel (e.g., miR-375-SOX15-PPP4R1L) microRNA-gene-pseudogene associations in prostate cancer. PseudoFuN provides a “one stop shop” for identifying and visualizing thousands of potential regulatory relationships related to pseudogenes in The Cancer Genome Atlas cancers.

Conclusions

Thousands of new PGG associations can be explored in the context of microRNA-gene-pseudogene co-expression and differential expression with a simple-to-use online tool by bioinformaticians and oncologists alike.

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

Background

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

Methods

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

Results

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

Conclusion

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

Retroposed copies of RET gene: a somatically acquired event in medullary thyroid carcinoma

BACKGROUND

Different pathogenic germline mutations in the RET oncogene are identified in MEN 2, a hereditary syndrome characterized by medullary thyroid carcinoma (MTC) and other endocrine tumors. Although genetic predisposition is recognized, not all RET mutation carriers will develop the disease during their lifetime or, likewise, RET mutation carriers belonging to the same family may present clinical heterogeneity. It has been suggested that a single germline mutation might not be sufficient for development of MEN 2-associated tumors and a somatic bi-allelic alteration might be required. Here we investigated the presence of somatic second hit mutation in the RET gene in MTC.

METHODS

We integrated Multiplex Ligation-dependent Probe Amplification (MLPA) and whole exome sequencing (WES) to search for copy number alteration (CNA) in the RET gene in MTC samples and medullary thyroid cell lines (TT and MZ-CR-1). We next found reads spanning exon-exon boundaries on RET, an indicative of retrocopy. We subsequently searched for RET retrocopies in the human reference genome (GRCh37) and in the 1000 Genomes Project data, by looking for reads reporting joined exons in the RET locus or distinct genomic regions. To determine RET retrocopy specificity and recurrence, DNA isolated from sporadic and MEN 2-associated MTC (n = 37), peripheral blood (n = 3) and papillary thyroid carcinomas with RET fusion (n = 10) samples were tested using PCR-sequencing methodology.

RESULTS

Through MLPA we have found evidence of CNA in the RET gene in MTC samples and MTC cell lines. WES analysis reinforced the presence of the CNA and hinted for a retroposed copy of RET not found in the human reference genome and 1.000 Genomes Project. Extended analysis confirmed the presence of a somatic MTC-related retrocopy of RET in both sporadic and hereditary tumors. We further unveiled a recurrent (28%) novel point mutation (p.G548 V) found exclusively in the retrocopy of RET. The mutation was also found in cDNA of mutated samples, suggesting it might be functional.

CONCLUSION

We here report a somatic specific RET retroposed copy in MTC samples and cell lines. Our results support the idea that generation of retrocopies in somatic cells is likely to contribute to MTC genesis and progression.

Pseudogene RACGAP1P activates RACGAP1/Rho/ERK signalling axis as a competing endogenous RNA to promote hepatocellular carcinoma early recurrence

Accumulating evidence has indicated crucial roles for pseudogenes in human cancers. However, the roles played by pseudogenes in the pathogenesis of HCC, particularly HCC early recurrence, still incompletely elucidated. Herein, we identify a novel early recurrence related pseudogene RACGAP1P which was significantly upregulated in HCC and was associated with larger tumour size, advanced clinical stage, abnormal AFP level and shorter survival time. In vitro and in vivo experiments have shown that RACGAP1P is a prerequisite for the development of malignant characteristics of HCC cells, including cell growth and migration. Mechanistic investigations indicated that RACGAP1P elicits its oncogenic activity as a ceRNA to sequestrate miR-15-5p from its endogenous target RACGAP1, thereby leading to the upregulation of RACGAP1 and the activation of RhoA/ERK signalling. These results may provide new insights into the functional crosstalk of the pseudogene/miRNA/parent-gene genetic network during HCC early relapse and may contribute to improving the clinical intervention for this subset of HCC patients.

The Function of Non-Coding RNAs in Lung Cancer Tumorigenesis

Lung cancer is the most prevalent and deadliest cancer worldwide. A significant part of lung cancer studies is dedicated to the expression alterations of non-coding RNAs. The non-coding RNAs are transcripts that cannot be translated into proteins. While the study of microRNAs and siRNAs in lung cancer received a lot of attention over the last decade, highly efficient therapeutic option or the diagnostic methods based on non-coding RNAs are still lacking. Because of this, it is of utmost importance to direct future research on lung cancer towards analyzing other RNA types for which the currently available data indicates that are essential at modulating lung tumorigenesis. Through our review of studies on this subject, we identify the following non-coding RNAs as tumor suppressors: ts-46, ts-47, ts-101, ts-53, ts-3676, ts-4521 (tRNA fragments), SNORD116-26, HBII-420, SNORD15A, SNORA42 (snoRNAs), piRNA-like-163, piR-35127, the piR-46545 (piRNAs), CHIAP2, LOC100420907, RPL13AP17 (pseudogenes), and uc.454 (T-UCR). We also found non-coding RNAs with tumor-promoting function: tRF-Leu-CAG, tRNA-Leu, tRNA-Val (tRNA fragments), circ-RAD23B, circRNA 100146, circPVT1, circFGFR3, circ_0004015, circPUM1, circFLI1, circABCB10, circHIPK3 (circRNAs), SNORA42, SNORA3, SNORD46, SNORA21, SNORD28, SNORA47, SNORD66, SNORA68, SNORA78 (snoRNAs), piR-65, piR-34871, piR-52200, piR651 (piRNAs), hY4 5’ fragments (YRNAs), FAM83A-AS1, WRAP53, NKX2-1-AS1 (NATs), DUXAP8, SFTA1P (pseudogene transcripts), uc.338, uc.339 (T-UCRs), and hTERC.

Pseudogenes as Functionally Significant Elements of the Genome

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

Making a virtue of necessity: the pleiotropic role of human endogenous retroviruses in cancer

Like all other mammals, humans harbour an astonishing number of endogenous retroviruses (ERVs), as well as other retroelements, embedded in their genome. These remnants of ancestral germline infection with distinct exogenous retroviruses display various degrees of open reading frame integrity and replication capability. Modern day exogenous retroviruses, as well as the infectious predecessors of ERVs, are demonstrably oncogenic. Further, replication-competent ERVs continue to cause cancers in many other species of mammal. Moreover, human cancers are characterized by transcriptional activation of human endogenous retroviruses (HERVs). These observations conspire to incriminate HERVs as causative agents of human cancer. However, exhaustive investigation of cancer genomes suggests that HERVs have entirely lost the ability for re-infection and thus the potential for insertional mutagenic activity. Although there may be non-insertional mechanisms by which HERVs contribute to cancer development, recent evidence also uncovers potent anti-tumour activities exerted by HERV replication intermediates or protein products. On balance, it appears that HERVs, despite their oncogenic past, now represent potential targets for immune-mediated anti-tumour mechanisms.

This article is part of the themed issue ‘Human oncogenic viruses’.

Network analysis of pseudogene-gene relationships: from pseudogene evolution to their functional potentials

Pseudogenes are fossil relatives of genes. Pseudogenes have long been thought of as "junk DNAs", since they do not code proteins in normal tissues. Although most of the human pseudogenes do not have noticeable functions, ∼20% of them exhibit transcriptional activity. There has been evidence showing that some pseudogenes adopted functions as lncRNAs and work as regulators of gene expression. Furthermore, pseudogenes can even be "reactivated" in some conditions, such as cancer initiation. Some pseudogenes are transcribed in specific cancer types, and some are even translated into proteins as observed in several cancer cell lines. All the above have shown that pseudogenes could have functional roles or potentials in the genome. Evaluating the relationships between pseudogenes and their gene counterparts could help us reveal the evolutionary path of pseudogenes and associate pseudogenes with functional potentials. It also provides an insight into the regulatory networks involving pseudogenes with transcriptional and even translational activities.In this study, we develop a novel approach integrating graph analysis, sequence alignment and functional analysis to evaluate pseudogene-gene relationships, and apply it to human gene homologs and pseudogenes. We generated a comprehensive set of 445 pseudogene-gene (PGG) families from the original 3,281 gene families (13.56%). Of these 438 (98.4% PGG, 13.3% total) were non-trivial (containing more than one pseudogene). Each PGG family contains multiple genes and pseudogenes with high sequence similarity. For each family, we generate a sequence alignment network and phylogenetic trees recapitulating the evolutionary paths. We find evidence supporting the evolution history of olfactory family (both genes and pseudogenes) in human, which also supports the validity of our analysis method. Next, we evaluate these networks in respect to the gene ontology from which we identify functions enriched in these pseudogene-gene families and infer functional impact of pseudogenes involved in the networks. This demonstrates the application of our PGG network database in the study of pseudogene function in disease context.