The OCT4 pseudogene POU5F1B is amplified and promotes an aggressive phenotype in gastric cancer

Targeting transposable elements in cancer: developments and opportunities

Transposable elements (TEs), comprising nearly 50% of the human genome, have transitioned from being perceived as "genomic junk" to key players in cancer progression. Contemporary research links TE regulatory disruptions with cancer development, underscoring their therapeutic potential. Advances in long-read sequencing, computational analytics, single-cell sequencing, proteomics, and CRISPR-Cas9 technologies have enriched our understanding of TEs' clinical implications, notably their impact on genome architecture, gene regulation, and evolutionary processes. In cancer, TEs, including long interspersed element-1 (LINE-1), Alus, and long terminal repeat (LTR) elements, demonstrate altered patterns, influencing both tumorigenic and tumor-suppressive mechanisms. TE-derived nucleic acids and tumor antigens play critical roles in tumor immunity, bridging innate and adaptive responses. Given their central role in oncology, TE-targeted therapies, particularly through reverse transcriptase inhibitors and epigenetic modulators, represent a novel avenue in cancer treatment. Combining these TE-focused strategies with existing chemotherapy or immunotherapy regimens could enhance efficacy and offer a new dimension in cancer treatment. This review delves into recent TE detection advancements, explores their multifaceted roles in tumorigenesis and immune regulation, discusses emerging diagnostic and therapeutic approaches centered on TEs, and anticipates future directions in cancer research.

Enhancing breast cancer outcomes with machine learning-driven glutamine metabolic reprogramming signature

Background

This study aims to identify precise biomarkers for breast cancer to improve patient outcomes, addressing the limitations of traditional staging in predicting treatment responses.

Methods

Our analysis encompassed data from over 7,000 breast cancer patients across 14 datasets, which included in-house clinical data and single-cell data from 8 patients (totaling 43,766 cells). We utilized an integrative approach, applying 10 machine learning algorithms in 54 unique combinations to analyze 100 existing breast cancer signatures. Immunohistochemistry assays were performed for empirical validation. The study also investigated potential immunotherapies and chemotherapies.

Results

Our research identified five consistent glutamine metabolic reprogramming (GMR)-related genes from multi-center cohorts, forming the foundation of a novel GMR-model. This model demonstrated superior accuracy in predicting recurrence and mortality risks compared to existing clinical and molecular features. Patients classified as high-risk by the model exhibited poorer outcomes. IHC validation in 30 patients reinforced these findings, suggesting the model's broad applicability. Intriguingly, the model indicates a differential therapeutic response: low-risk patients may benefit more from immunotherapy, whereas high-risk patients showed sensitivity to specific chemotherapies like BI-2536 and ispinesib.

Conclusions

The GMR-model marks a significant leap forward in breast cancer prognosis and the personalization of treatment strategies, offering vital insights for the effective management of diverse breast cancer patient populations.

Breaking paradigms: Long non-coding RNAs forming gene fusions with potential implications in cancer

Long non-coding RNAs (lncRNAs) are non-coding RNAs longer than 200 nucleotides with dynamic regulatory functions. They interact with a wide range of molecules such as DNA, RNA, and proteins to modulate diverse cellular functions through several mechanisms and, if deregulated, they can lead to cancer development and progression. Recently, it has been described that lncRNAs are susceptible to form gene fusions with mRNAs or other lncRNAs, breaking the paradigm of gene fusions consisting mainly of protein-coding genes. However, their biological significance in the tumor phenotype is still uncertain. Therefore, their recent identification opens a new line of research to study their biological role in tumorigenesis, and their potential as biomarkers with clinical relevance or as therapeutic targets. The present study aimed to review the lncRNA fusions identified so far and to know which of them have been associated with a potential function. We address the current challenges to deepen their study as well as the reasons why they represent a future therapeutic window in cancer.

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

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

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSION

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

Chromosomal Instability in Gastric Cancer: Role in Tumor Development, Progression, and Therapy

According to the Cancer Genome Atlas (TCGA), gastric cancers are classified into four molecular subtypes: Epstein-Barr virus-positive (EBV+), tumors with microsatellite instability (MSI), tumors with chromosomal instability (CIN), and genomically stable (GS) tumors. However, the gastric cancer (GC) with chromosomal instability remains insufficiently described and does not have effective markers for molecular and histological verification and diagnosis. The CIN subtype of GC is characterized by chromosomal instability, which is manifested by an increased frequency of aneuploidies and/or structural chromosomal rearrangements in tumor cells. Structural rearrangements in the CIN subtype of GC are not accidental and are commonly detected in chromosomal loci, being abnormal because of specific structural organization. The causes of CIN are still being discussed; however, according to recent data, aberrations in the TP53 gene may cause CIN development or worsen its phenotype. Clinically, patients with the CIN subtype of GC demonstrate poor survival, but receive the maximum benefit from adjuvant chemotherapy. In the review, we consider the molecular mechanisms and possible causes of chromosomal instability in GC, the common rearrangements of chromosomal loci and their impact on the development and clinical course of the disease, as well as the driver genes, their functions, and perspectives on their targeting in the CIN subtype of GC.

Zyxin inhibits the epithelial-mesenchymal transition process in gastric cancer by upregulating SIRT1

Tumor development relies on the stemness of cancer stem cells, which is regulated by environmental cues. Previous studies have shown that zyxin can inhibit the expression of genes for embryonic stem cell status. In the present study, the expression levels of zyxin protein in cancer tissues and adjacent noncancerous tissues from 73 gastric cancer patients with different clinical stages were analyzed by Western blot. We showed that the relative expression levels of zyxin in gastric cancer tissues (cancer tissues/adjacent tissues) were significantly downregulated in advanced clinical stages. Overexpression of zyxin inhibited the stemness and epithelial-mesenchymal transition (EMT) processes in gastric cancer cells. Zyxin also inhibited the proliferation, migration, and invasion but increased the sensitivity of cancer cells to drugs. Overexpression of zyxin in MKN45 cells inhibited tumor growth in nude mice. We show that the interactions between zyxin and SIRT1 led to the upregulation of SIRT1, reduced acetylation levels of histone H3 K9 and K23, decreased transcription levels of SNAI 1/2, and inhibition of the EMT process. This study demonstrated that zyxin negatively regulates the progression of gastric cancer by inhibiting the stemness of cancer stem cells and EMT. Our findings shed new light on the pathogenesis of gastric cancer.

PLAGL2 promotes Snail expression and gastric cancer progression via UCA1/miR-145-5p/YTHDF1 axis

OBJECTIVES

Although great progress has made in gastric cancer (GC) in the past years, the overall 5-year survival rate remains to be low for advanced GC patients. A recent study showed that PLAGL2 was increased in GC and enhanced the proliferation and metastasis of GC. Nevertheless, the underlying mechanism still needs to be investigated.

METHODS

Gene and protein expressions were assessed using RT-qPCR and western blot. The migration, proliferation and invasion of GC cells were examined using scratch assay, CCK-8 assay and Transwell assay, respectively. ChIP-PCR, dual-luciferase assay, RIP-qPCR and CoiP were utilized to confirm the interaction among PLAGL2, UCA1, miR-145-5p and YTHDF1 as well as METTL3, YTHDF1 and eEF-2. A mouse xenograft model was used utilized to further confirm the regulatory network.

RESULTS

PLAGL2 bound to the upstream promoter of UCA1, which regulated YTHDF1 by sponging miR-145-5p. METTL3 can mediate the m6A modification level of Snail. YTHDF1 recognized m6A-modified Snail by interacting with eEF-2 and thus promoted Snail expression, which eventually induced epithelial-mesenchymal transition (EMT) in GC cells and metastasis of GC.

CONCLUSIONS

Overall, our study demonstrates that PLAGL2 enhances Snail expression and GC progression via the UCA1/miR-145-5p/YTHDF1 axis, suggesting that PLAGL2 may become a therapeutic target for GC treatment.

HERV-K (HML-2) insertion polymorphisms in the 8q24.13 region and their potential etiological associations with acute myeloid leukemia

Human endogenous retroviruses (HERVs) are LTR retrotransposons that are present in the human genome. Among them, members of the HERV-K (HML-2) group are suspected to play a role in the development of different types of cancer, including lung, ovarian, and prostate cancer, as well as leukemia. Acute myeloid leukemia (AML) is an important disease that causes 1% of cancer deaths in the United States and has a survival rate of 28.7%. Here, we describe a method for assessing the statistical association between HERV-K (HML-2) transposable element insertion polymorphisms (or TIPs) and AML, using whole-genome sequencing and read mapping using TIP_finder software. Our results suggest that 101 polymorphisms involving HERV-K (HML-2) elements were correlated with AML, with a percentage between 44.4 to 56.6%, most of which (70) were located in the region from 8q24.13 to 8q24.21. Moreover, it was found that the TRIB1, LRATD2, POU5F1B, MYC, PCAT1, PVT1, and CCDC26 genes could be displaced or fragmented by TIPs. Furthermore, a general method was devised to facilitate analysis of the correlation between transposable element insertions and specific diseases. Finally, although the relationship between HERV-K (HML-2) TIPs and AML remains unclear, the data reported in this study indicate a statistical correlation, as supported by the χ² test with p-values < 0.05.

The 8q24 region hosts miRNAs altered in biospecimens of colorectal and bladder cancer patients

BACKGROUND

The 8q24 locus is enriched in cancer-associated polymorphisms and, despite containing relatively few protein-coding genes, it hosts the MYC oncogene and other genetic elements connected to tumorigenesis, including microRNAs (miRNAs). Research on miRNAs may provide insights into the transcriptomic regulation of this multiple cancer-associated region.

MATERIAL AND METHODS

We profiled all miRNAs located in the 8q24 region in 120 colorectal cancer (CRC) patients and 80 controls. miRNA profiling was performed on cancer/non-malignant adjacent mucosa, stool, and plasma extracellular vesicles (EVs), and the results validated with The Cancer Genome Atlas (TCGA) data. To verify if the 8q24-annotated miRNAs altered in CRC were dysregulated in other cancers and biofluids, we evaluated their levels in bladder cancer (BC) cases from the TCGA dataset and in urine and plasma EVs from a set of BC cases and healthy controls.

RESULTS

Among the detected mature miRNAs in the region, 12 were altered between CRC and adjacent mucosa (adj. p < 0.05). Five and four miRNAs were confirmed as dysregulated in the CRC and BC TCGA dataset, respectively. A co-expression analysis of tumor/adjacent tissue data from the CRC group revealed a correlation between the dysregulated miRNAs and CRC-related genes (PVT1 and MYC) annotated in 8q24 region. miR-30d-5p and miR-151a-3p, altered in CRC tissue, were also dysregulated in stool of CRC patients and urine of BC cases, respectively. Functional enrichment of dysregulated miRNA target genes highlighted terms related to TP53-mediated cell cycle control.

CONCLUSIONS

Altered expression of 8q24-annotated miRNAs may be relevant for the initiation and/or progression of cancer.

The alarmin interleukin-33 promotes the expansion and preserves the stemness of Tcf-1+ CD8+ T cells in chronic viral infection

T cell factor 1 (Tcf-1) expressing CD8⁺ T cells exhibit stem-like self-renewing capacity, rendering them key for immune defense against chronic viral infection and cancer. Yet, the signals that promote the formation and maintenance of these stem-like CD8⁺ T cells (CD8⁺SL) remain poorly defined. Studying CD8⁺ T cell differentiation in mice with chronic viral infection, we identified the alarmin interleukin-33 (IL-33) as pivotal for the expansion and stem-like functioning of CD8⁺SL as well as for virus control. IL-33 receptor (ST2)-deficient CD8⁺ T cells exhibited biased end differentiation and premature loss of Tcf-1. ST2-deficient CD8⁺SL responses were restored by blockade of type I interferon signaling, suggesting that IL-33 balances IFN-I effects to control CD8⁺SL formation in chronic infection. IL-33 signals broadly augmented chromatin accessibility in CD8⁺SL and determined these cells' re-expansion potential. Our study identifies the IL-33-ST2 axis as an important CD8⁺SL-promoting pathway in the context of chronic viral infection.

Transposon-activated POU5F1B promotes colorectal cancer growth and metastasis

The treatment of colorectal cancer (CRC) is an unmet medical need in absence of early diagnosis. Here, upon characterizing cancer-specific transposable element-driven transpochimeric gene transcripts (TcGTs) produced by this tumor in the SYSCOL cohort, we find that expression of the hominid-restricted retrogene POU5F1B through aberrant activation of a primate-specific endogenous retroviral promoter is a strong negative prognostic biomarker. Correlating this observation, we demonstrate that POU5F1B fosters the proliferation and metastatic potential of CRC cells. We further determine that POU5F1B, in spite of its phylogenetic relationship with the POU5F1/OCT4 transcription factor, is a membrane-enriched protein that associates with protein kinases and known targets or interactors as well as with cytoskeleton-related molecules, and induces intracellular signaling events and the release of trans-acting factors involved in cell growth and cell adhesion. As POU5F1B is an apparently non-essential gene only lowly expressed in normal tissues, and as POU5F1B-containing TcGTs are detected in other tumors besides CRC, our data provide interesting leads for the development of cancer therapies.

Considering hormone-sensitive cancers as a single disease in the UK biobank reveals shared aetiology

Hormone-related cancers, including cancers of the breast, prostate, ovaries, uterine, and thyroid, globally contribute to the majority of cancer incidence. We hypothesize that hormone-sensitive cancers share common genetic risk factors that have rarely been investigated by previous genomic studies of site-specific cancers. Here, we show that considering hormone-sensitive cancers as a single disease in the UK Biobank reveals shared genetic aetiology. We observe that a significant proportion of variance in disease liability is explained by the genome-wide single nucleotide polymorphisms (SNPs), i.e., SNP-based heritability on the liability scale is estimated as 10.06% (SE 0.70%). Moreover, we find 55 genome-wide significant SNPs for the disease, using a genome-wide association study. Pair-wise analysis also estimates positive genetic correlations between some pairs of hormone-sensitive cancers although they are not statistically significant. Our finding suggests that heritable genetic factors may be a key driver in the mechanism of carcinogenesis shared by hormone-sensitive cancers.

Exploring the effects of genetic variation on gene regulation in cancer in the context of 3D genome structure

Background

Numerous genome-wide association studies (GWAS) conducted to date revealed genetic variants associated with various diseases, including breast and prostate cancers. Despite the availability of these large-scale data, relatively few variants have been functionally characterized, mainly because the majority of single-nucleotide polymorphisms (SNPs) map to the non-coding regions of the human genome. The functional characterization of these non-coding variants and the identification of their target genes remain challenging.

Results

In this communication, we explore the potential functional mechanisms of non-coding SNPs by integrating GWAS with the high-resolution chromosome conformation capture (Hi-C) data for breast and prostate cancers. We show that more genetic variants map to regulatory elements through the 3D genome structure than the 1D linear genome lacking physical chromatin interactions. Importantly, the association of enhancers, transcription factors, and their target genes with breast and prostate cancers tends to be higher when these regulatory elements are mapped to high-risk SNPs through spatial interactions compared to simply using a linear proximity. Finally, we demonstrate that topologically associating domains (TADs) carrying high-risk SNPs also contain gene regulatory elements whose association with cancer is generally higher than those belonging to control TADs containing no high-risk variants.

Conclusions

Our results suggest that many SNPs may contribute to the cancer development by affecting the expression of certain tumor-related genes through long-range chromatin interactions with gene regulatory elements. Integrating large-scale genetic datasets with the 3D genome structure offers an attractive and unique approach to systematically investigate the functional mechanisms of genetic variants in disease risk and progression.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12863-021-01021-x.

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.

Detection and identification of cis-regulatory elements using change-point and classification algorithms

BACKGROUND

Transcriptional regulation is primarily mediated by the binding of factors to non-coding regions in DNA. Identification of these binding regions enhances understanding of tissue formation and potentially facilitates the development of gene therapies. However, successful identification of binding regions is made difficult by the lack of a universal biological code for their characterisation.

RESULTS

We extend an alignment-based method, changept, and identify clusters of biological significance, through ontology and de novo motif analysis. Further, we apply a Bayesian method to estimate and combine binary classifiers on the clusters we identify to produce a better performing composite.

CONCLUSIONS

The analysis we describe provides a computational method for identification of conserved binding sites in the human genome and facilitates an alternative interrogation of combinations of existing data sets with alignment data.

Recurrent Superenhancer of the Oncogene POU5F1B in Colorectal Cancers

Superenhancer usages in single cancer form such as colorectal cancer (CRC) may provide novel efficient targeting candidates. It is unclear whether CRC contains recurrent superenhancers that confer a predisposition to malignancy. We investigated the superenhancer profile of CRC cell line HCT116 and compared it to that of a healthy sigmoid colon. We found that HCT116 had lost most of the normal colon superenhancer activities but gained a new set of tumor-favoring superenhancers that facilitate tumor proliferation, growth signalling, and hypoxia resistance. Inhibiting the superenhancers by JQ-1 treatment had significantly decreased the colony formation capability of HCT116. Then, by comparing the superenhancer genes and robust CRC upregulated genes, we identified a superenhancer associated with a common CRC upregulated oncogene, POU5f1B. POU5f1B overexpression is related to the worse outcome in CRCs. Via performing ChIP-PCR in 35 clinical samples and investigating CRC anti-H3K27ac ChiP-seq public dataset consisting of 36 samples, we further identified that the superenhancer of oncogene POU5F1B is recurrently activated in CRCs, taking 62 and 72 per cent, respectively. Moreover, JQ-1 treatment successfully inhibited the POU5F1B expression in 5 out of 6 POU5F1B superenhancer-positive samples. Therefore, we concluded that the superenhancer activation of POU5F1B contributes partially to its high expression in CRCs, in addition to the well-known gene amplification aetiology.

Construction and Clinical Translation of Causal Pan-Cancer Gene Score Across Cancer Types

Pan-cancer strategy, an integrative analysis of different cancer types, can be used to explain oncogenesis and identify biomarkers using a larger statistical power and robustness. Fine-mapping defines the casual loci, whereas genome-wide association studies (GWASs) typically identify thousands of cancer-related loci and not necessarily have a fine-mapping component. In this study, we develop a novel strategy to identify the causal loci using a pan-cancer and fine-mapping assumption, constructing the CAusal Pan-cancER gene (CAPER) score and validating its performance using internal and external validation on 1,287 individuals and 985 cell lines. Summary statistics of 15 cancer types were used to define 54 causal loci in 15 potential genes. Using the Cancer Genome Atlas (TCGA) training set, we constructed the CAPER score and divided cancer patients into two groups. Using the three validation sets, we found that 19 cancer-related variables were statistically significant between the two CAPER score groups and that 81 drugs had significantly different drug sensitivity between the two CAPER score groups. We hope that our strategies for selecting causal genes and for constructing CAPER score would provide valuable clues for guiding the management of different types of cancers.

The Roles of the Colon Cancer Associated Transcript 2 (CCAT2) Long Non-Coding RNA in Cancer: A Comprehensive Characterization of the Tumorigenic and Molecular Functions

Colon cancer-associated transcript 2 (CCAT2) is an intensively studied lncRNA with important regulatory roles in cancer. As such, cumulative studies indicate that CCAT2 displays a high functional versatility due to its direct interaction with multiple RNA binding proteins, transcription factors, and other species of non-coding RNA, especially microRNA. The definitory mechanisms of CCAT2 are its role as a regulator of the TCF7L2 transcription factor, enhancer of MYC expression, and activator of the WNT/β-catenin pathway, as well as a role in promoting and maintaining chromosome instability through the BOP1–AURKB pathway. Additionally, we highlight how the encompassing rs6983267 SNP has been shown to confer CCAT2 with allele-specific functional and structural particularities, such as the allelic-specific reprogramming of glutamine metabolism. Additionally, we emphasize CCAT2’s role as a competitive endogenous RNA (ceRNA) for multiple tumor suppressor miRNAs, such as miR-4496, miR-493, miR-424, miR-216b, miR-23b, miR-34a, miR-145, miR-200b, and miR-143 and the pro-tumorigenic role of the altered regulatory axis. Additionally, due to its upregulation in tumor tissues, wide distribution across cancer types, and presence in serum samples, we outline CCAT2’s potential as a biomarker and disease indicator and its implications for the development of resistance against current cancer therapy regiments and metastasis.

Pseudogene Annexin A2 Pseudogene 1 Contributes to Hepatocellular Carcinoma Progression by Modulating Its Parental Gene ANXA2 via miRNA-376a-3p

BACKGROUND

Pseudogenes are defined as key regulators in cancer initiation and progression. But their biological function and clinical significance in hepatocellular carcinoma (HCC) remain to be elucidated. In the current study, we identified a novel pseudogene, Annexin A2 pseudogene 1 (ANXA2P1), in HCC and explored its underlining molecular mechanism.

METHODS AND RESULTS

We analyzed the expression pattern of ANXA2P1 in a TCGA dataset and an HCC sample cohort and evaluated its clinical significance. The biological effects on HCC cells proliferation, migration, invasion and epithelial-mesenchymal transition (EMT) process were assessed by Cell Counting Kit-8 assay, Transwell assay and Western blot, respectively. The ANXA2P1/miR-376a-3p/ANXA2 axis was determined by bioinformatics analysis and dual-luciferase reporter assays. ANXA2P1 exerted as an oncogene that was significantly overexpressed in HCC tissues and was associated with disease progression and unfavorable prognosis of HCC patients. ANXA2P1 knockdown suppressed cell growth, cell migration and invasion and reversed EMT phenotype in HCC. Mechanistically, ANXA2P1 acts as a competing endogenous RNA for miR-376a-3p, thereby leading to the upregulation of its cognate gene ANXA2.

CONCLUSIONS

ANXA2P1/miR-376a-3p/ANXA2 axis plays an important role in the progression of HCC. Our findings may provide valuable therapeutic target for treating HCC.

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.

The ELF3-regulated lncRNA UBE2CP3 is over-stabilized by RNA-RNA interactions and drives gastric cancer metastasis via miR-138-5p/ITGA2 axis

LncRNAs play essential roles in tumorigenesis and tumor progression. Pseudogene UBE2CP3 is an antisense intronic lncRNA. However, the biological function of UBE2CP3 in gastric cancer (GC) remains unknown. In this study, we revealed that lncRNA UBE2CP3 was aberrantly upregulated in multiple independent gastric cancer cohorts, and its overexpression was clinically associated with poor prognosis in GC. UBE2CP3 was mainly located in cytoplasm and promoted migratory and invasive capacities of GC cells in vitro and in vivo. Mechanismly, a novel dysregulated ceRNA network UB2CP3/miR-138-5p/ITGA2 was identified in GC by transcriptome sequencing. Furthermore, rescue assay further confirmed that UBE2CP3 mainly promoted GC progression through miR-138-5p/ITGA2 axis. More importantly, our data proved that UBE2CP3/IGFBP7 could form an RNA duplex, thereby directly interacting with the ILF3 protein. In turn, this RNA-RNA interaction between IGFBP7 mRNA and UBE2CP3 mediated by ILF3 protein plays an essential role in protecting the mRNA stability of UBE2CP3. In addition, transcription factor ELF3 was identified to be a direct repressor of lncRNA UBE2CP3 in GC. Taken together, overexpression of UBE2CP3 promotes tumor progression via cascade amplification of ITGA2 upregulation in GC. Our finding has revealed that the dysregulation of UBE2CP3 is probably due to the downregulation of ELF3 and/or the overexpression of IGFBP7 mRNA in GC. Our findings reveal, for the first time, that UBE2CP3 plays crucial a role in GC progression by modulating miR-138-5p/ITGA2 axis, suggesting that UBE2CP3 may serve as a potential therapeutic target in GC.

LncRNA CASC21 induces HGH1 to mediate colorectal cancer cell proliferation, migration, EMT and stemness

Colorectal cancer (CRC) is one of the leading causes of cancer-related deaths worldwide. Long non-coding RNAs (lncRNAs) have been increasingly reported to serve vital parts in malignancies including CRC. Although cancer susceptibility 21 (CASC21) has been uncovered to play a part in CRC, its mechanism still needs further explanation. Thus, our study aimed to further explore the influence and mechanism of CASC21 in CRC progression. Quantitative real-time RT-PCR and western blot were performed to detect gene expression; a series of functional assays were performed to investigate the effect of CASC21 on CRC cells; in vivo tumour growth was evaluated via the nude mice xenograft model. The results revealed that CASC21 facilitated CRC cell proliferation, migration, epithelial-mesenchymal transition (EMT) and stemness. In addition, CASC21 was co-expressed with and bound to transcription factor POU5F1B (POU class 5 homeobox 1B). CASC21 recruited POU5F1B to HGH1 promoter to activate the transcription of HGH1 homolog. Also, CASC21 served as a competitive endogenous RNA (ceRNA) to up-regulate HGH1 via endogenously sponging miR-485-5p. Moreover, HGH1 overexpression counteracted the suppression of CASC21 deficiency on CRC tumour growth. In summary, our study indicated that CASC21 enhanced the expression of HGH1 to promote the malignancy of CRC by recruiting POU5F1B and sponging miR-485-5p, suggesting a key role of CASC21 in CRC progression.

Genetic association studies of alterations in protein function expose recessive effects on cancer predisposition

The characterization of germline genetic variation affecting cancer risk, known as cancer predisposition, is fundamental to preventive and personalized medicine. Studies of genetic cancer predisposition typically identify significant genomic regions based on family-based cohorts or genome-wide association studies (GWAS). However, the results of such studies rarely provide biological insight or functional interpretation. In this study, we conducted a comprehensive analysis of cancer predisposition in the UK Biobank cohort using a new gene-based method for detecting protein-coding genes that are functionally interpretable. Specifically, we conducted proteome-wide association studies (PWAS) to identify genetic associations mediated by alterations to protein function. With PWAS, we identified 110 significant gene-cancer associations in 70 unique genomic regions across nine cancer types and pan-cancer. In 48 of the 110 PWAS associations (44%), estimated gene damage is associated with reduced rather than elevated cancer risk, suggesting a protective effect. Together with standard GWAS, we implicated 145 unique genomic loci with cancer risk. While most of these genomic regions are supported by external evidence, our results also highlight many novel loci. Based on the capacity of PWAS to detect non-additive genetic effects, we found that 46% of the PWAS-significant cancer regions exhibited exclusive recessive inheritance. These results highlight the importance of recessive genetic effects, without relying on familial studies. Finally, we show that many of the detected genes exert substantial cancer risk in the studied cohort determined by a quantitative functional description, suggesting their relevance for diagnosis and genetic consulting.

Prognostic Nomogram of Prognosis-Related Genes and Clinicopathological Characteristics to Predict the 5-Year Survival Rate of Colon Cancer Patients

Background: The Cancer Genome Atlas (TCGA) has established a genome-wide gene expression profile, increasing our understanding of the impact of tumor heredity on clinical outcomes. The aim of this study was to construct a nomogram using data from the TCGA regarding prognosis-related genes and clinicopathological characteristics to predict the 5-years survival rate of colon cancer (CC) patients. Methods: Kaplan-Meier and Cox regression analyses were used to identify genes associated with the 5-years survival rate of CC patients. Cox regression was used to analyze the relationship between the clinicopathological features and prognostic genes and overall survival rates in patients with CC and to identify independent risk factors for the prognosis of CC patients. A nomogram for predicting the 5-years survival rate of CC patients was constructed by R software. Results: A total of eight genes (KCNJ14, CILP2, ATP6V1G2, GABRD, RIMKLB, SIX2, PLEKHA8P1, and MPP2) related to the 5-years survival of rate CC patients were identified. Age, stage, and PLEKHA8P1 were independent risk factors for the 5-years survival rate in patients with CC. The accuracy, sensitivity and specificity of the nomogram model constructed by age, TNM staging, and PLEKHA8P1 for predicting the 5-years survival of rate CC patients were 83.3, 83.97, and 85.79%, respectively. Conclusion: The nomogram can correctly predict the 5-year survival rate of patients with CC, thus aiding the individualized decision-making process for patients with CC.

Markers and Reporters to Reveal the Hierarchy in Heterogeneous Cancer Stem Cells

A subpopulation within cancer, known as cancer stem cells (CSCs), regulates tumor initiation, chemoresistance, and metastasis. At a closer look, CSCs show functional heterogeneity and hierarchical organization. The present review is an attempt to assign marker profiles to define the functional heterogeneity and hierarchical organization of CSCs, based on a series of single-cell analyses. The evidences show that analogous to stem cell hierarchy, self-renewing Quiescent CSCs give rise to the Progenitor CSCs with limited proliferative capacity, and later to a Progenitor-like CSCs, which differentiates to Proliferating non-CSCs. Functionally, the CSCs can be tumor-initiating cells (TICs), drug-resistant CSCs, or metastasis initiating cells (MICs). Although there are certain marker profiles used to identify CSCs of different cancers, molecules like CD44, CD133, ALDH1A1, ABCG2, and pluripotency markers [Octamer binding transcriptional factor 4 (OCT4), SOX2, and NANOG] are used to mark CSCs of a wide range of cancers, ranging from hematological malignancies to solid tumors. Our analysis of the recent reports showed that a combination of these markers can demarcate the heterogeneous CSCs in solid tumors. Reporter constructs are widely used for easy identification and quantification of marker molecules. In this review, we discuss the suitability of reporters for the widely used CSC markers that can define the heterogeneous CSCs. Since the CSC-specific functions of CD44 and CD133 are regulated at the post-translational level, we do not recommend the reporters for these molecules for the detection of CSCs. A promoter-based reporter for ABCG2 may also be not relevant in CSCs, as the expression of the molecule in cancer is mainly regulated by promoter demethylation. In this context, a dual reporter consisting of one of the pluripotency markers and ALDH1A1 will be useful in marking the heterogeneous CSCs. This system can be easily adapted to high-throughput platforms to screen drugs for eliminating CSCs.

Co-amplification of genes in chromosome 8q24: a robust prognostic marker in hepatocellular carcinoma

Background

Hepatocellular carcinoma (HCC) is a leading cause of tumor-associated death worldwide, owing to its high 5-year postoperative recurrence rate and inter-individual heterogeneity. Thus, a prognostic model is urgently needed for patients with HCC. Several researches have reported that copy number amplification of the 8q24 chromosomal region is associated with low survival in many cancers. In the present work, we set out to construct a multi-gene model for prognostic prediction in HCC.

Methods

RNA sequencing and copy number variant data of tumor tissue samples of HCC from The Cancer Genome Atlas (n=328) were used to identify differentially expressed messenger RNAs of genes located on the chromosomal 8q24 region by the Wilcox test. Univariate Cox and Lasso-Cox regression analyses were carried out for the screening and construction of a prognostic multi-gene signature in The Cancer Genome Atlas cohort (n=119). The multi-gene signature was validated in a cohort from the International Cancer Genome Consortium (n=240). A nomogram for prognostic prediction was built, and the underpinning molecular mechanisms were studied by Gene Set Enrichment Analysis.

Results

We successfully established a 7-gene prognostic signature model to predict the prognosis of patients with HCC. Using the model, we divided individuals into high-risk and low-risk sets, which showed a significant difference in overall survival in the training dataset (HR =0.17, 95% CI: 0.1-0.28; P<0.001) and in the testing dataset (HR = 0.42, 95% CI: 0.23-0.74; P=0.002). Multivariate Cox regression analysis showed the signature to be an independent prognostic factor of HCC survival. A nomogram including the prognostic signature was constructed and showed a better predictive performance in short-term (1 and 3 years) than in long-term (5 years) survival. Furthermore, Gene Set Enrichment Analysis identified several pathways of significance, which may aid in explaining the underlying molecular mechanism.

Conclusions

Our 7-gene signature is a reliable prognostic marker for HCC, which may provide meaningful information for therapeutic customization and treatment-related decision making.

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

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

The Landscape of the Heritable Cancer Genome

Genome-wide association studies (GWAS) have found hundreds of single-nucleotide polymorphisms (SNP) associated with increased risk of cancer. However, the amount of heritable risk explained by SNPs is limited, leaving most of the cancer heritability unexplained. Tumor sequencing projects have shown that causal mutations are enriched in genic regions. We hypothesized that SNPs located in protein coding genes and nearby regulatory regions could explain a significant proportion of the heritable risk of cancer. To perform gene-level heritability analysis, we developed a new method, called Bayesian Gene Heritability Analysis (BAGHERA), to estimate the heritability explained by all genotyped SNPs and by those located in genic regions using GWAS summary statistics. BAGHERA was specifically designed for low heritability traits such as cancer and provides robust heritability estimates under different genetic architectures. BAGHERA-based analysis of 38 cancers reported in the UK Biobank showed that SNPs explain at least 10% of the heritable risk for 14 of them, including late onset malignancies. We then identified 1,146 genes, called cancer heritability genes (CHG), explaining a significant proportion of cancer heritability. CHGs were involved in hallmark processes controlling the transformation from normal to cancerous cells. Importantly, 60 of them also harbored somatic driver mutations, and 27 are tumor suppressors. Our results suggest that germline and somatic mutation information could be exploited to identify subgroups of individuals at higher risk of cancer in the broader population and could prove useful to establish strategies for early detection and cancer surveillance. SIGNIFICANCE: This study describes a new statistical method to identify genes associated with cancer heritability in the broader population, creating a map of the heritable cancer genome with gene-level resolution.See related commentary by Bader, p. 2586.

Microarray-based Analysis of Genes, Transcription Factors, and Epigenetic Modifications in Lung Cancer Exposed to Nitric Oxide

BACKGROUND/AIM

Nitric oxide (NO) is recognized as an important biological mediator that exerts several human physiological functions. As its nature is an aqueous soluble gas that can diffuse through cells and tissues, NO can affect cell signaling, the phenotype of cancer and modify surrounding cells. The variety of effects of NO on cancer cell biology has convinced researchers to determine the defined mechanisms of these effects and how to control this mediator for a better understanding as well as for therapeutic gain.

MATERIALS AND METHODS

We used bioinformatics and pharmacological experiments to elucidate the potential regulation and underlying mechanisms of NO in non-small a lung cancer cell model.

RESULTS

Using microarrays, we identified a total of 151 NO-regulated genes (80 up-regulated genes, 71 down-regulated genes) with a strong statistically significant difference compared to untreated controls. Among these, the genes activated by a factor of more than five times were: DCBLD2, MGC24975, RAB40AL, PER3, RCN1, MRPL51, PTTG1, KLF5, NFIX. On the other hand, the expression of RBMS2, PDP2, RBAK, ORMDL2, GRPEL2, ZNF514, MTHFD2, POLR2D, RCBTB1, JOSD1, RPS27, GPR4 genes were significantly decreased by a factor of more than five times. Bioinformatics further revealed that NO exposure of lung cancer cells resulted in a change in transcription factors (TFs) and epigenetic modifications (histone modification and miRNA). Interestingly, NO treatment was shown to potentiate cancer stem cell-related genes and transcription factors Oct4, Klf4, and Myc.

CONCLUSION

Through this comprehensive approach, the present study illustrated the scheme of how NO affects molecular events in lung cancer cells.

8q24.21 Locus: A Paradigm to Link Non-Coding RNAs, Genome Polymorphisms and Cancer

The majority of the human genome is comprised of non-protein-coding genes, but the relevance of non-coding RNAs in complex diseases has yet to be fully elucidated. One class of non-coding RNAs is long non-coding RNAs or lncRNAs, many of which have been identified to play a range of roles in transcription and translation. While the clinical importance of the majority of lncRNAs have yet to be identified, it is puzzling that a large number of disease-associated genetic variations are seen in lncRNA genes. The 8q24.21 locus is rich in lncRNAs and very few protein-coding genes are located in this region. Interestingly, the 8q24.21 region is also a hot spot for genetic variants associated with an increased risk of cancer. Research focusing on the lncRNAs in this area of the genome has indicated clinical relevance of lncRNAs in different cancers. In this review, we summarise the lncRNAs in the 8q24.21 region with respect to their role in cancer and discuss the potential impact of cancer-associated genetic polymorphisms on the function of lncRNAs in initiation and progression of cancer.

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.

A novel Oct4/Pou5f1-like non-coding RNA controls neural maturation and mediates developmental effects of ethanol

Prenatal ethanol exposure can result in loss of neural stem cells (NSCs) and decreased brain growth. Here, we assessed whether a noncoding RNA (ncRNA) related to the NSC self-renewal factor Oct4/Pou5f1, and transcribed from a processed pseudogene locus on mouse chromosome 9 (mOct4pg9), contributed to the loss of NSCs due to ethanol. Mouse fetal cortical-derived NSCs, cultured ex vivo to mimic the early neurogenic environment of the fetal telencephalon, expressed mOct4pg9 ncRNA at significantly higher levels than the parent Oct4/Pou5f1 mRNA. Ethanol exposure increased expression of mOct4pg9 ncRNA, but decreased expression of Oct4/Pou5f1. Gain- and loss-of-function analyses indicated that mOct4pg9 overexpression generally mimicked effects of ethanol exposure, resulting in increased proliferation and expression of transcripts associated with neural maturation. Moreover, mOct4pg9 associated with Ago2 and with miRNAs, including the anti-proliferative miR-328-3p, whose levels were reduced following mOct4pg9 overexpression. Finally, mOct4pg9 inhibited Oct4/Pou5f1-3'UTR-dependent protein translation. Consistent with these observations, data from single-cell transcriptome analysis showed that mOct4pg9-expressing progenitors lack Oct4/Pou5f1, but instead overexpress transcripts for increased mitosis, suggesting initiation of transit amplification. Collectively, these data suggest that the inhibitory effects of ethanol on brain development are explained, in part, by a novel ncRNA which promotes loss of NSC identity and maturation.

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.

Nintedanib promotes antitumour immunity and shows antitumour activity in combination with PD-1 blockade in mice: potential role of cancer-associated fibroblasts

BACKGROUND

Cancer-associated fibroblasts (CAFs) in the tumour microenvironment (TME) suppress antitumour immunity, and the tyrosine kinase inhibitor nintedanib has antifibrotic effects.

METHODS

We performed a preclinical study to evaluate whether nintedanib might enhance antitumour immunity by targeting CAFs and thereby improve the response to immune checkpoint blockade (ICB).

RESULTS

Whereas nintedanib did not suppress the growth of B16-F10 melanoma cells in vitro, it prolonged survival in a syngeneic mouse model of tumour formation by these cells, suggestive of an effect on the TME without direct cytotoxicity. Gene expression profiling indeed showed that nintedanib influenced antitumour immunity and fibrosis. Tumoural infiltration of CD8⁺ T cells and granzyme B production were increased by nintedanib, and its antitumour activity was attenuated by antibody-mediated depletion of these cells, indicating that nintedanib suppressed tumour growth in a CD8⁺ T cell-dependent manner. Moreover, nintedanib inhibited the proliferation and activation of fibroblasts. Finally, the combination of nintedanib with ICB showed enhanced antitumour efficacy in B16-F10 tumour-bearing mice.

CONCLUSIONS

Our results suggest that nintedanib targeted CAFs and thereby attenuated the immunosuppressive nature of the TME and promoted the intratumoural accumulation and activation of CD8⁺ T cells, with these effects contributing to enhanced antitumour activity in combination with ICB.

HPV-16 E2/E6 and POU5F1B as Biomarkers to Determine Cervical High-Grade Squamous Lesions and More

Purpose

Human papillomavirus-16 (HPV-16) is the most carcinogenic HPV genotype. This study aimed to evaluate the clinical value of POU5F1B and HPV-16-E2/E6 by cervical cytology specimens to predict the cervical intraepithelial neoplasia two grade and more (CIN2+).

Methods

Finally, 248 patients with HPV-16 single infection were enrolled. Using cytology specimen by real-time quantitative PCR (qPCR), POU5F1B mRNA and HPV-16-E2/E6 were detected. The relationship of POU5F1B, HPV-16-E2/E6 and CIN2+ were analyzed, and the optimal cut-off values of POU5F1B and HPV-16-E2/E6 to predict CIN2+ were calculated.

Results

The mean HPV-16-E2/E6 decreased significantly with cervical lesions development, especially compared with CIN2+ (p<0.05). And the POU5F1B demonstrated higher expression in CIN2+ than that of normal cervical tissue and CIN1 (p<0.05). What is more, POU5F1B was negatively correlated with HPV-16-E2/E6. It demonstrated that the area under the receiver operating characteristic curve (AUC) for POU5F1B (0.9058) was higher than that for HPV-16-E2/E6 (0.8983), and the sensitivity and specificity of POU5F1B in the diagnosis of CIN2+ were higher than HPV-E2/E6. Furthermore, it demonstrated that the POU5F1B had the highest odds ratio (OR= 16.84; 95% CI (8.00–35.46)) for the detection of CIN 2+.

Conclusion

HPV-16-E2/E6≤0.6471 or POU5F1B≥1.0310 in cervical exfoliated cells can be used as a reliable predictor of CIN2+. POU5F1B can be used as a new auxiliary biomarker to determine the HPV infection status and a reliable predictor of CIN2+. The expression of POU5F1B≥1.0310 had the highest OR for the detection of CIN2+.

Exploring the effects of genetic variation on gene regulation in cancer in the context of 3D genome structure.. [DOI: 10.1101/2020.10.06.328567] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Numerous genome-wide association studies (GWAS) conducted to date revealed genetic variants associated with various diseases, including breast and prostate cancers. Despite the availability of these large-scale data, relatively few variants have been functionally characterized, mainly because the majority of single-nucleotide polymorphisms (SNPs) map to the non-coding regions of the human genome. The functional characterization of these non-coding variants and the identification of their target genes remain challenging. In this communication, we explore the potential functional mechanisms of non-coding SNPs by integrating GWAS with the high-resolution chromosome conformation capture (Hi-C) data for breast and prostate cancers. We show that more genetic variants map to regulatory elements through the 3D genome structure than the 1D linear genome lacking physical chromatin interactions. Importantly, the association of enhancers, transcription factors, and their target genes with breast and prostate cancers tends to be higher when these regulatory elements are mapped to high-risk SNPs through spatial interactions compared to simply using a linear proximity. Finally, we demonstrate that topologically associating domains (TADs) carrying high-risk SNPs also contain gene regulatory elements whose association with cancer is generally higher than those belonging to control TADs containing no high-risk variants. Our results suggest that many SNPs may contribute to the cancer development by affecting the expression of certain tumor-related genes through long-range chromatin interactions with gene regulatory elements. Integrating large-scale genetic datasets with the 3D genome structure offers an attractive and unique approach to systematically investigate the functional mechanisms of genetic variants in disease risk and progression.

Targeting of Lung Cancer Stem Cell Self-Renewal Pathway by a Small Molecule Verrucarin J

Despite considerable advances made in understanding of lung cancer biology, there has been meek improvement in lung cancer treatment outcome with 4% to 5% increase in 5-year survival rates in the last four decades. Underlying problem of lung cancer recurrence and poor prognosis is attributed to the presence of cancer stem cells (CSCs) which possess the potential to differentiate, proliferate and trigger chemo-resistance, tumor progression and metastasis, despite initial elimination of the tumor. To address specific targeting of CSCs, we investigated the effects of a small molecule Verrucarin J (VJ) on lung cancer cell lines A549 and H1793. VJ significantly inhibited cell proliferation of both cell lines, with IC₅₀ values of approximately 10 nM for A549 and 20 nM for H1793 respectively after 48 h of treatment. A549 cell line when treated with VJ, induced cell apoptosis with concomitant down regulation of key CSC specific genes- ALDH1, LGR5, OCT4 and CD133 in a dose-dependent manner. To delineate the molecular mechanism by which VJ targets lung cancer cells and CSCs, we determined the effects of VJ on CSC self-renewal pathways Wnt1/β-catenin and Notch1. Treatment of A549 cell line with VJ inhibited significantly both the signalling pathways, suggesting inhibition of expression of CSC genes by VJ through the inhibition of CSC self-renewal signalling pathways. Taken together, our results suggest that VJ may serve as a potent anticancer drug to target cancer cells and CSCs.

PWAS: proteome-wide association study-linking genes and phenotypes by functional variation in proteins

We introduce Proteome-Wide Association Study (PWAS), a new method for detecting gene-phenotype associations mediated by protein function alterations. PWAS aggregates the signal of all variants jointly affecting a protein-coding gene and assesses their overall impact on the protein's function using machine learning and probabilistic models. Subsequently, it tests whether the gene exhibits functional variability between individuals that correlates with the phenotype of interest. PWAS can capture complex modes of heritability, including recessive inheritance. A comparison with GWAS and other existing methods proves its capacity to recover causal protein-coding genes and highlight new associations. PWAS is available as a command-line tool.

A three-gene signature might predict prognosis in patients with acute myeloid leukemia

The identification of effective signatures is crucial to predict the prognosis of acute myeloid leukemia (AML). The investigation aimed to identify a new signature for AML prognostic prediction by using the three-gene expression (octamer-binding transcription factor 4 (OCT4), POU domain type 5 transcription factor 1B (POU5F1B) and B-cell-specific Moloney murine leukemia virus integration site-1 pseudogene 1 (BMI1P1). The expressions of genes were obtained from our previous study. Only the specimens in which three genes were all expressed were included in this research. A three-gene signature was constructed by the multivariate Cox regression analyses to divide patients into high-risk and low-risk groups. Receiver operating characteristic (ROC) analysis of the three-gene signature (area under ROC curve (AUC) = 0.901, 95% CI: 0.821–0.981, P<0.001) indicated that it was a more valuable signature for distinguishing between patients and controls than any of the three genes. Moreover, white blood cells (WBCs, P=0.004), platelets (PLTs, P=0.017), percentage of blasts in bone marrow (BM) (P=0.011) and complete remission (CR, P=0.027) had significant differences between two groups. Furthermore, high-risk group had shorter leukemia-free survival (LFS) and overall survival (OS) than low-risk group (P=0.026; P=0.006), and the three-gene signature was a prognostic factor. Our three-gene signature for prognosis prediction in AML may serve as a prognostic biomarker.

The Oncogenic Potential of the Centromeric Border Protein FAM84B of the 8q24.21 Gene Desert

FAM84B is a risk gene in breast and prostate cancers. Its upregulation is associated with poor prognosis of prostate cancer, breast cancer, and esophageal squamous cell carcinoma. FAM84B facilitates cancer cell proliferation and invasion in vitro, and xenograft growth in vivo. The FAM84B and Myc genes border a 1.2 Mb gene desert at 8q24.21. Co-amplification of both occurs in 20 cancer types. Mice deficient of a 430 Kb fragment within the 1.2 Mb gene desert have downregulated FAM84B and Myc expressions concurrent with reduced breast cancer growth. Intriguingly, Myc works in partnership with other oncogenes, including Ras. FAM84B shares similarities with the H-Ras-like suppressor (HRASLS) family over their typical LRAT (lecithin:retinal acyltransferase) domain. This domain contains a catalytic triad, H23, H35, and C113, which constitutes the phospholipase A_1/2 and O-acyltransferase activities of HRASLS1-5. These enzymatic activities underlie their suppression of Ras. FAM84B conserves H23 and H35 but not C113 with both histidine residues residing within a highly conserved motif that FAM84B shares with HRASLS1-5. Deletion of this motif abolishes FAM84B oncogenic activities. These properties suggest a collaboration of FAM84B with Myc, consistent with the role of the gene desert in strengthening Myc functions. Here, we will discuss recent research on FAM84B-derived oncogenic potential.

Role of Oct3/4 in Cervical Cancer Tumorigenesis

Cervical cancer (CC) is the fourth most common type of cancer that affects women. Compared to other types of cancer, CC has a high mortality rate in women worldwide. Several factors contribute to the development of CC, but persistent high-risk human papillomavirus infection is the main etiologic agent associated with the development of CC. Moreover, several studies reported that alterations in the expression of transcription factors present in a small subpopulation of cells within tumors called cancer stem cells (CSCs), which contribute to the development of CC by promoting tumorigenicity and metastasis. These transcription factors affect self-renewal and maintenance of pluripotency and differentiation in stem cells. OCT3/4 belongs to the family of transcription factors with the POU domain. It consists of five exons and can be edited by alternative splicing into three main transcripts: OCT3/4A, OCT3/4B, and OCT3/4B1. The OCT3/4 expression in CSCs promotes carcinogenesis and the development of malignant tumors, and the loss of expression leads to the loss of self-renewal and proliferation and favors apoptosis. This review describes the main roles of OCT3/4 in CC and its importance in several biological processes that contribute to the development of CC and may serve as molecular targets to improve prognosis of CC.

The Complex Network between MYC Oncogene and microRNAs in Gastric Cancer: An Overview

Despite the advancements in cancer treatments, gastric cancer is still one of the leading causes of death worldwide. In this context, it is of great interest to discover new and more effective ways of treating this disease. Accumulated evidences have demonstrated the amplification of 8q24.21 region in gastric tumors. Furthermore, this is the region where the widely known MYC oncogene and different microRNAs are located. MYC deregulation is key in tumorigenesis in various types of tissues, once it is associated with cell proliferation, survival, and drug resistance. microRNAs are a class of noncoding RNAs that negatively regulate the protein translation, and which deregulation is related with gastric cancer development. However, little is understood about the interactions between microRNAs and MYC. Here, we overview the MYC role and its relationship with the microRNAs network in gastric cancer aiming to identify potential targets useful to be used in clinic, not only as biomarkers, but also as molecules for development of promising therapies.

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

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

Long non-coding RNA CCAT2 promotes oncogenesis in triple-negative breast cancer by regulating stemness of cancer cells

Triple-negative breast cancers (TNBC) are more aggressive due to lacking receptors for hormone therapy and maintaining stemness features in cancer cells. Herein we found long non-coding RNA CCAT2 overexpressed specially in TNBC, and in breast cancer stem cells (BCSC) as well. Enforced overexpression and targeted knockdown demonstrated the oncogenic function of CCAT2 both in vitro and in vivo. CCAT2 promoted the expression of stemness markers including OCT4, Nanog and KLF4, increased mammosphere formation and induced ALDH+ cancer stem cell population in TNBC. A chromosomally adjacent gene OCT4-PG1, as a pseudogene of OCT4, was upregulated by CCAT2, and positively regulated the stemness features of TNBC cells. miR-205 was identified as a target gene of CCAT2 in TNBC. Point-mutation in CCAT2 impaired the sponge inhibition of miR-205. Overexpression of miR-205 rescued the oncogenic phenotypes induced by CCAT2. In addition, Notch2, as a target gene of miR-205, was downregulated by miR-205 and upregulated by CCAT2 in TNBC. Collectively, the current study revealed a novel function of CCAT2 in promoting tumor initiation and progression in TNBC through upregulating OCT4-PG1 expression and activating Notch signaling. These findings not only demonstrated a lncRNA-based therapeutic strategy in treatment of TNBC, but also added a node to the regulatory network of CCAT2 that controls aggressiveness of breast cancer stem cells.

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.

The Contribution of Genetic Variants to the Risk of Papillary Thyroid Carcinoma in the Kazakh Population: Study of Common Single Nucleotide Polymorphisms and Their Clinicopathological Correlations

OBJECTIVE

Risk for developing papillary thyroid carcinoma (PTC), the most common endocrine malignancy, is thought to be mediated by lifestyle, environmental exposures and genetic factors. Recent progress in the genome-wide association studies of thyroid cancer leads to the identification of several genetic variants conferring risk to this malignancy across different ethnicities. We set out to elucidate the impact of selected single nucleotide polymorphisms (SNPs) on PTC risk and to evaluate clinicopathological correlations of these genetic variants in the Kazakh population for the first time.

METHODS

Eight SNPs were genotyped in 485 patients with PTC and 1,008 healthy control Kazakh subjects. The association analysis and multivariable modeling of PTC risk by the genetic factors, supplemented with rigorous statistical validation, were performed.

RESULT

Five of the eight SNPs: rs965513 (FOXE1/PTCSC2, P = 1.3E-16), rs1867277 (FOXE1 5'UTR, P = 7.5E-06), rs2439302 (NRG1 intron 1, P = 4.0E-05), rs944289 (PTCSC3/NKX2-1, P = 4.5E-06) and rs10136427 (BATF upstream, P = 9.8E-03) were significantly associated with PTC. rs966423 (DIRC3, P = 0.07) showed a suggestive association. rs7267944 (DHX35) was associated with PTC risk in males (P = 0.02), rs1867277 (FOXE1) conferred the higher risk in subjects older than 55 years (P = 7.0E-05), and rs6983267 (POU5F1B/CCAT2) was associated with pT3-T4 tumors (P = 0.01). The contribution of genetic component (unidirectional independent effects of rs965513, rs944289, rs2439302 and rs10136427 adjusted for age and sex) to PTC risk in the analyzed series was estimated to be 30-40%.

CONCLUSION

Genetic factors analyzed in the present work display significant association signals with PTC either on the whole group analysis or in particular clinicopathological groups and account for about one-third of the risk for PTC in the Kazakh population.

DUX4 Suppresses MHC Class I to Promote Cancer Immune Evasion and Resistance to Checkpoint Blockade

Advances in cancer immunotherapies make it critical to identify genes that modulate antigen presentation and tumor-immune interactions. We report that DUX4, an early embryonic transcription factor that is normally silenced in somatic tissues, is re-expressed in diverse solid cancers. Both cis-acting inherited genetic variation and somatically acquired mutations in trans-acting repressors contribute to DUX4 re-expression in cancer. Although many DUX4 target genes encode self-antigens, DUX4-expressing cancers were paradoxically characterized by reduced markers of anti-tumor cytolytic activity and lower major histocompatibility complex (MHC) class I gene expression. We demonstrate that DUX4 expression blocks interferon-γ-mediated induction of MHC class I, implicating suppressed antigen presentation in DUX4-mediated immune evasion. Clinical data in metastatic melanoma confirmed that DUX4 expression was associated with significantly reduced progression-free and overall survival in response to anti-CTLA-4. Our results demonstrate that cancers can escape immune surveillance by reactivating a normal developmental pathway and identify a therapeutically relevant mechanism of cell-intrinsic immune evasion.