POU5F1P1, a putative cancer susceptibility gene, is overexpressed in prostatic carcinoma

Identification of prostate cancer associated genes for diagnosis and prognosis: a modernized in silico approach

Prostate cancer (PCa) ranks as the second leading cause of cancer-related deaths in men. Diagnosing PCa relies on molecular markers known as diagnostic biomarkers, while prognostic biomarkers are used to identify key proteins involved in PCa treatments. This study aims to gather PCa-associated genes and assess their potential as either diagnostic or prognostic biomarkers for PCa. A corpus of 152,064 PCa-related data from PubMed, spanning from May 1936 to December 2020, was compiled. Additionally, 4199 genes associated with PCa terms were collected from the National Center of Biotechnology Information (NCBI) database. The PubMed corpus data was extracted using pubmed.mineR to identify PCa-associated genes. Network and pathway analyses were conducted using various tools, such as STRING, DAVID, KEGG, MCODE 2.0, cytoHubba app, CluePedia, and ClueGO app. Significant marker genes were identified using Random Forest, Support Vector Machines, Neural Network algorithms, and the Cox Proportional Hazard model. This study reports 3062 unique PCa-associated genes along with 2518 corresponding unique PMIDs. Diagnostic markers such as IL6, MAPK3, JUN, FOS, ACTB, MYC, and TGFB1 were identified, while prognostic markers like ACTB and HDAC1 were highlighted in PubMed. This suggests that the potential target genes provided by PubMed data outweigh those in the NCBI database.

Proteomics Can Rise to the Challenge of Pseudogenes' Coding Nature

Throughout the past decade, technological advances in genomics and transcriptomics have revealed pervasive translation throughout mammalian genomes. These putative proteins are usually excluded from proteomics analyses, as they are absent from common protein repositories. A sizable portion of these noncanonical proteins is translated from pseudogenes. Pseudogenes are commonly termed defective copies of coding genes unable to produce proteins. Here, we suggest that proteomics can help in their annotation. First, we define important terms and review specific examples underlining the caveats in pseudogene annotation and their coding potential. Then, we will discuss the challenges inherent to pseudogenes that have thus far rendered complex their confidence in omics data. Finally, we identify recent developments in experimental procedures, instrumentation, and computational methods in proteomics that put the field in a unique position to solve the pseudogene annotation conundrum.

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.

Unravelling similarities and differences in the role of circular and linear PVT1 in cancer and human disease

The plasmacytoma variant translocation 1 (PVT1) is a long non-coding RNA gene involved in human disease, mainly in cancer onset/progression. Although widely analysed, its biological roles need to be further clarified. Notably, functional studies on PVT1 are complicated by the occurrence of multiple transcript variants, linear and circular, which generate technical issues in the experimental procedures used to evaluate its impact on human disease. Among the many PVT1 transcripts, the linear PVT1 (lncPVT1) and the circular hsa_circ_0001821 (circPVT1) are frequently reported to perform similar pathologic and pro-tumorigenic functions when overexpressed. The stimulation of cell proliferation, invasion and drug resistance, cell metabolism regulation, and apoptosis inhibition is controlled through multiple targets, including MYC, p21, STAT3, vimentin, cadherins, the PI3K/AKT, HK2, BCL2, and CASP3. However, some of this evidence may originate from an incorrect evaluation of these transcripts as two separate molecules, as they share the lncPVT1 exon-2 sequence. We here summarise lncPVT1/circPVT1 functions by mainly focusing on shared pathways, pointing out the potential bias that may exist when the biological role of each transcript is analysed. These considerations may improve the knowledge about lncPVT1/circPVT1 and their specific targets, which deserve further studies due to their diagnostic, prognostic, and therapeutic potential.

MYC DNA Methylation in Prostate Tumor Tissue Is Associated with Gleason Score

Increasing evidence suggests a role of epigenetic mechanisms at chromosome 8q24, an important cancer genetic susceptibility region, in prostate cancer. We investigated whether MYC DNA methylation at 8q24 (six CpG sites from exon 3 to the 3′ UTR) in prostate tumor was associated with tumor aggressiveness (based on Gleason score, GS), and we incorporated RNA expression data to investigate the function. We accessed radical prostatectomy tissue for 50 Caucasian and 50 African American prostate cancer patients at the University of Maryland Medical Center, selecting an equal number of GS 6 and GS 7 cases per group. MYC DNA methylation was lower in tumor than paired normal prostate tissue for all six CpG sites (median difference: −14.74 to −0.20 percentage points), and we observed similar results for two nearby sites in The Cancer Genome Atlas (p < 0.0001). We observed significantly lower methylation for more aggressive (GS 7) than less aggressive (GS 6) tumors for three exon 3 sites (for CpG 212 (chr8:128753145), GS 6 median = 89.7%; GS 7 median = 85.8%; p-value = 9.4 × 10⁻⁴). MYC DNA methylation was not associated with MYC expression, but was inversely associated with PRNCR1 expression after multiple comparison adjustment (q-value = 0.04). Findings suggest that prostate tumor MYC exon 3 hypomethylation is associated with increased aggressiveness.

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+.

High Cancer Susceptibility Candidate 8 Expression Is Associated With Poor Prognosis of Pancreatic Adenocarcinoma: Validated Analysis Based on Four Cancer Databases

Objective

The aim of this study was to explore the association between the expression of a long non-coding RNA (lncRNA), cancer susceptibility candidate 8 (CASC8), and pancreatic adenocarcinoma (PAAD).

Materials and Methods

starBase database was used to perform differential expression, survival, and competing endogenous RNA (ceRNA) network and H19/miR-671 correlation analyses for CASC8 in 178 PAAD samples. Using the cBioPortal database website, we analyzed the alteration in CASC8 expression and its correlation with the overall survival in PAAD. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were also performed using the circlncRNAnet database. Analysis of CASC8 polymorphisms was performed using the UCSC Xena database. Finally, the expression of CASC8 in Chinese PAAD tissues was validated by qPCR.

Results

The expression of CASC8 was observed to be high in 178 PAAD samples [fold change = 8.71, P = 0.0014, false discovery rate (FDR) = 0.04] and was related with poor prognosis, but not in pancreatic neuroendocrine tumor (pNET). CASC8 amplification was noted in 6% of the PAAD patients; however, the gene amplification did not affect the expression of CASC8 but was involved with the overall survival time of PAAD patients. Network analysis indicated that H19 is the ceRNA pair of CASC8 and that CASC8 competitively binds to miR-671 and might participate in the process of epithelial-to-mesenchymal transition (EMT). The correlation analysis showed that CASC8 was significantly negatively correlated with SMAD7. The analysis of CASC8 polymorphism showed that high copy number segment (CNS) of CASC8 is associated with low survival. Validation using PAAD tissues from Chinese patients was consistent with the in silico findings.

Conclusion

CASC8 is specifically expressed at a high level in PAAD and associated with poor prognosis, which might be through its interaction with H19, miR-671, and SMAD7. These results indicate that CASC8 could serve as a novel marker for predicting the prognosis and as a potential target for the therapy of PAAD.

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.

The Octamer-Binding Transcription Factor 4 (OCT4) Pseudogene, POU Domain Class 5 Transcription Factor 1B (POU5F1B), is Upregulated in Cervical Cancer and Down-Regulation Inhibits Cell Proliferation and Migration and Induces Apoptosis in Cervical Cancer Cell Lines

Background

The POU domain class 5 transcription factor 1B (POU5F1B), is a pseudogene that is homologous to octamer-binding transcription factor 4 (OCT4), and is located adjacent to the MYC gene on human chromosome 8q24. POU5F1B has been reported to be transcribed in several types of cancer, but its role in cervical cancer remains unclear. This study aimed to investigate the expression and function of POU5F1B in tissue samples of human cervical cancer and in cervical cancer cell lines in vitro.

Material/Methods

Quantitative reverse transcription polymerase chain reaction (qRT-PCR) was used to quantify POU5F1B expression in cervical cancer tissues and in SiHa, HeLa, CaSki, and C33A human cervical cancer cell lines. Functional in vitro studies included analysis of the effects of POU5F1B expression on cervical cancer cell proliferation, migration, and apoptosis using a Cell Counting Kit-8 (CCK-8) assay, cell migration assays, and flow cytometry. Luciferase activity assays, qRT-PCR, and Western blot were performed to confirm the expression of POU5F1B.

Results

POU5F1B was significantly upregulated in cervical cancer tissues and cell lines. Interference with the expression of POU5F1B significantly inhibited cell proliferation, apoptosis, migration and invasion, and induced apoptosis in vitro. Western blot demonstrated that POU5F1B could modulate the expression of the OCT4 protein.

Conclusions

POU5F1B was upregulated in cervical cancer and down-regulation inhibited cell proliferation and migration and induced apoptosis in cervical cancer cell lines by modulating OCT4. Further studies are required to determine whether POU5F1B might be a diagnostic or prognostic biomarker or therapeutic target in cervical cancer.

Deletion of the murine ortholog of the 8q24 gene desert has anti-cancer effects in transgenic mammary cancer models

Background

The gene desert on human chromosomal band 8q24 harbors multiple genetic variants associated with common cancers, including breast cancer. The locus, including the gene desert and its flanking genes, MYC, PVT1 and FAM84B, is also frequently amplified in human breast cancer. We generated a megadeletion (MD) mouse model lacking 430-Kb of sequence orthologous to the breast cancer-associated region in the gene desert. The goals were to examine the effect of the deletion on mammary cancer development and on transcript level regulation of the candidate genes within the locus.

Methods

The MD allele was engineered using the MICER system in embryonic stem cells and bred onto 3 well-characterized transgenic models for breast cancer, namely MMTV-PyVT, MMTV-neu and C3(1)-TAg. Mammary tumor growth, latency, multiplicity and metastasis were compared between homozygous MD and wild type mice carrying the transgenes. A reciprocal mammary gland transplantation assay was conducted to distinguish mammary cell-autonomous from non-mammary cell-autonomous anti-cancer effects. Gene expression analysis was done using quantitative real-time PCR. Chromatin interactions were evaluated by 3C. Gene-specific patient outcome data were analysed using the METABRIC and TCGA data sets through the cBioPortal website.

Results

Mice homozygous for the MD allele are viable, fertile, lactate sufficiently to nourish their pups, but maintain a 10% lower body weight mainly due to decreased adiposity. The deletion interferes with mammary tumorigenesis in mouse models for luminal and basal breast cancer. In the MMTV-PyVT model the mammary cancer-reducing effects of the allele are mammary cell-autonomous. We found organ-specific effects on transcript level regulation, with Myc and Fam84b being downregulated in mammary gland, prostate and mammary tumor samples. Through analysis using the METABRIC and TCGA datasets, we provide evidence that MYC and FAM84B are frequently co-amplified in breast cancer, but in contrast with MYC, FAM84B is frequently overexpressed in the luminal subtype, whereas MYC activity affect basal breast cancer outcomes.

Conclusion

Deletion of a breast cancer-associated non-protein coding region affects mammary cancer development in 3 transgenic mouse models. We propose Myc as a candidate susceptibility gene, regulated by the gene desert locus, and a potential role for Fam84b in modifying breast cancer development.

Electronic supplementary material

The online version of this article (10.1186/s12885-018-5109-8) contains supplementary material, which is available to authorized users.

Single nucleotide polymorphisms and sporadic colorectal cancer susceptibility: a field synopsis and meta-analysis

Background

Although mounting non-hereditary colorectal cancer (NHCRC) associated single nucleotide polymorphisms (SNPs) have been observed, no field synopsis and meta-analysis has been conducted through systematically assessing cumulative evidence, during the past 5 years.

Methods

We retrieved the database via the PubMed, Web of Science and Embase gateways to identify publications concerning the associations between SNPs and risk of NHCRC, up to May 1st, 2017. To assess the finding credibility, cumulative evidence was graded based on the Venice criteria. Meta-analysis was also performed for three subgroups including ethnicity (Asian vs Caucasian), primary cancer site (colon vs rectum) and TNM stage (I II vs III IV). Then, we arranged those high quality SNPs into different regions according to their locations on genes to evaluate their functional roles on CRC development.

Results

5114 publications were collected and 1001 of them met our inclusion criteria, which totally included 1788 SNPs in 793 genes or distinct chromosomal loci. Totally, we performed 359 primary and subgroup meta-analyses for 160 SNPs in 96 distinct genes. By utilizing the Venice criteria, we identified 15 high quality SNPs with 25 high credibility significant associations. Furthermore, we artificially divided the high quality SNPs into different groups, based on their SNP loci (exon region, intron region, promoter region, downstream region, non-coding region and intergenic region).

Conclusion

We have identified 15 high quality SNPs which may act as promising genetic biomarkers for clinical NHCRC susceptibility screening and explored their functional roles on the NHCRC development based on their locations on genes.

Association between lncRNA CASC8 polymorphisms and the risk of cancer: a meta-analysis

Objective

To explore the relationship between single-nucleotide polymorphisms (SNPs) in one of the long noncoding RNA (lncRNA), cancer susceptibility candidate 8 (CASC8) gene and the risk of cancer.

Materials and methods

A meta-analysis was conducted to summarize the relationship between common SNPs (rs10505477 and rs7837328) in the lncRNA CASC8 gene and the risk of cancer. The relevant references were retrieved from several authoritative databases. Rigorous inclusion and exclusion criteria were adopted to ensure the credibility of the results. The fixed effects or random effects model was used to calculate the OR and 95% CI. We tested for publication bias.

Results

Fifteen articles containing 20 datasets (24,504 cases and 22,969 controls) were finally included in the meta-analysis. Compared to the individuals carrying the rs10505477 TT genotype, those with the TC or CC genotype had a decreased risk of cancer (TC vs TT: OR 0.876, 95% CI 0.832–0.923, P<0.001; CC vs TT: OR 0.748, 95% CI 0.703–0.795, P<0.001). Allele C of rs10505477 might be a protective factor for decreasing susceptibility to cancer (OR 0.866, 95% CI 0.840–0.893, P<0.001). As for rs7837328, the GA and AA genotypes were associated with increased risks of cancer as compared to the GG genotype (ORs 1.209 and 1.336; 95% CIs 1.127–1.298 and 1.202–1.484, respectively); its A allele could significantly increase the risk of cancer compared with the G allele (OR 1.169, 95% CI 1.114–1.227, P<0.001).

Conclusion

The rs10505477 and rs7837328 polymorphisms might be associated with risk of cancer.

Cumulative Evidence for Relationships Between 8q24 Variants and Prostate Cancer

Multiple independent cancer susceptibility loci at chromosome 8q24 have been identified by GWAS (Genome-wide association studies). Forty six articles including 60,293 cases and 62,971 controls were collected to conduct a meta-analysis to evaluate the associations between 21 variants in 8q24 and prostate cancer risk. Of the 21 variants located in 8q2\5 were significantly associated with the risk of prostate cancer. In particular, both homozygous AA and heterozygous CA genotypes of rs16901979, as well as the AA and CA genotypes of rs1447295, were associated with the risk of prostate cancer. Our study showed that variants in the 8q24 region are associated with prostate cancer risk in this large-scale research synopsis and meta-analysis. Further studies are needed to explore the role of the 8q24 variants involved in the etiology of prostate cancer.

Identification of driver copy number alterations in diverse cancer types and application in drug repositioning

Results from numerous studies suggest an important role for somatic copy number alterations (SCNAs) in cancer progression. Our work aimed to identify the drivers (oncogenes or tumor suppressor genes) that reside in recurrently aberrant genomic regions, including a large number of genes or non-coding genes, which remain a challenge for decoding the SCNAs involved in carcinogenesis. Here, we propose a new approach to comprehensively identify drivers, using 8740 cancer samples involving 18 cancer types from The Cancer Genome Atlas (TCGA). On average, 84 drivers were revealed for each cancer type, including protein-coding genes, long non-coding RNAs (lncRNA) and microRNAs (miRNAs). We demonstrated that the drivers showed significant attributes of cancer genes, and significantly overlapped with known cancer genes, including MYC, CCND1 and ERBB2 in breast cancer, and the lncRNA PVT1 in multiple cancer types. Pan-cancer analyses of drivers revealed specificity and commonality across cancer types, and the non-coding drivers showed a higher cancer-type specificity than that of coding drivers. Some cancer types from different tissue origins were found to converge to a high similarity because of the significant overlap of drivers, such as head and neck squamous cell carcinoma (HNSC) and lung squamous cell carcinoma (LUSC). The lncRNA SOX2-OT, a common driver of HNSC and LUSC, showed significant expression correlation with the oncogene SOX2. In addition, because some drivers are common in multiple cancer types and have been targeted by known drugs, we found that some drugs could be successfully repositioned, as validated by the datasets of drug response assays in cell lines. Our work reported a new method to comprehensively identify drivers in SCNAs across diverse cancer types, providing a feasible strategy for cancer drug repositioning as well as novel findings regarding cancer-associated non-coding RNA discovery.

Prospective study of DNA methylation at chromosome 8q24 in peripheral blood and prostate cancer risk

BACKGROUND

Chromosome 8q24 has emerged as an important genetic susceptibility region for several cancers, including prostate cancer; however, little is known about the contribution of DNA methylation in this region to risk.

METHODS

We prospectively evaluated DNA methylation at 8q24 in relation to prostate cancer using pre-diagnostic blood samples from 694 prostate cancer cases (including 172 aggressive cases) and 703 controls in the Prostate, Lung, Colorectal and Ovarian Cancer Screening Trial. We used logistic regression to estimate odds ratios and 95% confidence intervals.

RESULTS

Although none remained significant after adjustment for multiple testing (q>0.05), of the 50 CpG sites meeting quality control, we identified 8 sites that were nominally associated with prostate cancer (Ptrend<0.05), including 6 correlated (Spearman ρ: 0.20-0.52) sites in POU5F1B and 2 intergenic sites (most significant site: Chr8:128428897 in POU5F1B, Ptrend=0.01). We also identified two correlated (ρ=0.39) sites in MYC (Chr8:128753187 and Chr8:128753154) that were associated with aggressive (Ptrend=0.02 and 0.03), but not non-aggressive disease (Ptrend=0.70 and 0.20; Pheterogeneity=0.01 and 4.6 × 10-3). These findings persisted after adjustment for the top 8q24 prostate cancer variants in our study.

CONCLUSIONS

Although requiring replication, our findings provide some evidence that 8q24 DNA methylation levels may be associated with prostate cancer risk.

Upregulation of FAM84B during prostate cancer progression

Although the FAM84B gene lies within chromosome 8q24, a locus frequently altered in prostate cancer (PC), its alteration during prostate tumorigenesis has not been well studied. We report here FAM84B upregulation in DU145 cell-derived prostate cancer stem-like cells (PCSLCs) and DU145 cell-produced lung metastases compared to subcutaneous xenograft tumors. FAM84B protein was detected in bone metastases and primary PCs. Nanostring examination of 7 pairs of tumor adjacent normal and PC tissues revealed elevations in FAM84B mRNA levels in all carcinomas. Furthermore, through analysis of FAM84B expression using large datasets within the Gene Expression Omnibus and OncomineTM database, we demonstrate significant increases in FAM84B mRNA in 343 primary PCs versus 181 normal tissues, and elevations in the FAM84B gene copy number (GCN) in 171 primary PCs versus 61 normal tissues. While FAM84B was not detected at higher levels via immunohistochemistry in high grade (Gleason score/GS 8-10) tumors compared to GS6-7 PCs, analyses of FAM84B mRNA and GCN using datasets within the cBioPortal database demonstrated FAM84B upregulation in 12% (67/549) of primary PCs and 18% (73/412) of metastatic castration resistant PCs (mCRPCs), and GCN increases in 4.8% (26/546) of primary PCs and 26% (121/467) of mCRPCs, revealing an association of the aforementioned changes with CRPC development. Of note, an increase in FAM84B expression was observed in xenograft CRPCs produced by LNCaP cells. Furthermore, FAM84B upregulation and GCN increases correlate with decreases in disease free survival and overall survival. Collectively, we demonstrate a novel association of FAM84B with PC tumorigenesis and CRPC progression.

Fine-scale mapping of 8q24 locus identifies multiple independent risk variants for breast cancer

Previous genome-wide association studies among women of European ancestry identified two independent breast cancer susceptibility loci represented by single nucleotide polymorphisms (SNPs) rs13281615 and rs11780156 at 8q24. A fine-mapping study across 2.06 Mb (chr8:127,561,724-129,624,067, hg19) in 55,540 breast cancer cases and 51,168 controls within the Breast Cancer Association Consortium was conducted. Three additional independent association signals in women of European ancestry, represented by rs35961416 (OR = 0.95, 95% CI = 0.93-0.97, conditional p = 5.8 × 10(-6) ), rs7815245 (OR = 0.94, 95% CI = 0.91-0.96, conditional p = 1.1 × 10(-6) ) and rs2033101 (OR = 1.05, 95% CI = 1.02-1.07, conditional p = 1.1 × 10(-4) ) were found. Integrative analysis using functional genomic data from the Roadmap Epigenomics, the Encyclopedia of DNA Elements project, the Cancer Genome Atlas and other public resources implied that SNPs rs7815245 in Signal 3, and rs1121948 in Signal 5 (in linkage disequilibrium with rs11780156, r(2)  = 0.77), were putatively functional variants for two of the five independent association signals. The results highlighted multiple 8q24 variants associated with breast cancer susceptibility in women of European ancestry.

Expression of evolutionarily novel genes in tumors

The evolutionarily novel genes originated through different molecular mechanisms are expressed in tumors. Sometimes the expression of evolutionarily novel genes in tumors is highly specific. Moreover positive selection of many human tumor-related genes in primate lineage suggests their involvement in the origin of new functions beneficial to organisms. It is suggested to consider the expression of evolutionarily young or novel genes in tumors as a new biological phenomenon, a phenomenon of TSEEN (tumor specifically expressed, evolutionarily novel) genes.

Clinical Significance of Long Non-Coding RNA CASC8 rs10505477 Polymorphism in Lung Cancer Susceptibility, Platinum-Based Chemotherapy Response, and Toxicity

Long non-coding RNA (lncRNA) CASC8 rs10505477 polymorphism has been identified to be related to risk of many kinds of cancers, such as colorectal cancer, gastric cancer, and invasive ovarian cancer, and it may be involved in the prognosis of gastric cancer patients who have received platinum-based chemotherapy after surgical treatment. So far, there is no study investigating the clinical significance of lncRNA CASC8 rs10505477 in lung cancer susceptibility and treatment. In this study, we genotyped 498 lung cancer patients and 213 healthy control subjects to explore the correlation between the rs10505477 polymorphism and lung cancer risk in a Chinese population. Among the 498 patients, 467 were selected for the chemotherapy response and toxicity study. We found that the single nucleotide polymorphisms (SNP) rs10505477 was greatly related to lung cancer risk in male and adenocarcinoma subgroups in recessive model (adjusted OR = 0.51, 95%CI = 0.29–0.90, p = 0.02; adjusted OR = 0.52, 95%CI = 0.30–0.89, p = 0.02, respectively). It was also closely correlated with platinum-based chemotherapy response in dominant model (adjusted OR = 1.58, 95%CI = 1.05–2.39, p = 0.03). Additionally, we observed that CASC8 rs10505477 polymorphism was significantly relevant to severe hematologic toxicity in non-small-cell lung cancer (NSCLC) subgroup in dominant model (adjusted OR = 0.59, 95%CI = 0.35–0.98, p = 0.04) and in additive model (adjusted OR = 0.62, 95%CI = 0.43–0.90, p = 0.01). Furthermore, it was found that rs10505477 polymorphism was greatly associated with gastrointestinal toxicity in SCLC and cisplatin subgroups in dominant model (adjusted OR = 7.82, 95%CI = 1.36–45.07, p = 0.02; adjusted OR = 1.94, 95%CI = 1.07–3.53, p = 0.03, respectively). Thus, lncRNA CASC8 rs10505477 could serve as a possible risk marker for diagnosing lung cancer, and could be used to forecast the response and toxicity of platinum-based treatment in lung cancer patients.

An examination of clinical differences between carriers and non-carriers of chromosome 8q24 risk alleles in a New Zealand Caucasian population with prostate cancer

Background. Prostate cancer makes up approximately 15% of all cancers diagnosed in men in developed nations and approximately 4% of cases in developing nations. Although it is clear that prostate cancer has a genetic component and single nucleotide polymorphisms (SNPs) can contribute to prostate cancer risk, detecting associations is difficult in multi-factorial diseases, as environmental and lifestyle factors also play a role. In this study, specific clinical characteristics, environmental factors and genetic risk factors were assessed for interaction with prostate cancer.

Methods. A total of 489 prostate cancer cases and 427 healthy controls were genotyped for SNPs found on chromosome 8q24 and a genetic risk score was calculated. In addition the SNPs were tested for an association with a number of clinical and environmental factors.

Results. Age and tobacco use were positively associated, whilst alcohol consumption was negatively associated with prostate cancer risk. The following SNPs found on chromosome 8q24 were statistically significantly associated with prostate cancer: rs10086908, rs16901979; rs1447295and rs4242382. No association between Gleason score and smoking status, or between Gleason score and genotype were detected.

Conclusion. A genetic risk score was calculated based on the 15 SNPs tested and found to be significantly associated with prostate cancer risk. Smoking significantly contributed to the risk of developing prostate cancer, and this risk was further increased by the presence of four SNPs in the 8q24 chromosomal region.

Pseudogene-expressed RNAs: a new frontier in cancers

Over the past decade, the importance of non-protein-coding functional elements in the human genome has emerged from the water and been identified as a key revelation in post-genomic biology. Since the completion of the ENCODE (Encyclopedia of DNA Elements) and FANTOM (Functional Annotation of Mammals) project, tens of thousands of pseudogenes as well as numerous long non-coding RNA (lncRNA) genes were identified. However, while pseudogenes were initially regarded as non-functional relics littering the human genome during evolution, recent studies have revealed that they play critical roles at multiple levels in diverse physiological and pathological processes, especially in cancer through parental-gene-dependent or parental-gene-independent regulation. Herein, we review the current knowledge of pseudogenes and synthesize the nascent evidence for functional properties and regulatory modalities exerted by pseudogene-transcribed RNAs in human cancers and prospect the potential as molecular signatures in cancer reclassification and tailored therapy.

Analysis of Prostate Cancer Susceptibility Variants in South African Men: Replicating Associations on Chromosomes 8q24 and 10q11

Genome-wide association studies (GWAS) have implicated single nucleotide polymorphisms (SNPs) on chromosomes 2p15, 6q25, 7p15.2, 7q21, 8q24, 10q11, 10q26, 11q13, 17q12, 17q24, 19q13, and Xp11, with prostate cancer (PCa) susceptibility and/or tumour aggressiveness, in populations of African, European, and Asian ancestry. The objective of this study was to confirm these associations in South African Mixed Ancestry and White men. We evaluated 17 prioritised GWAS SNPs in South African cases (331 Mixed Ancestry and 155 White) and controls (178 Mixed Ancestry and 145 White). The replicated SNP associations for the different South African ethnic groups were rs7008482 (8q24) (p = 2.45 × 10⁻⁵), rs6983267 (8q24) (p = 4.48 × 10⁻⁷), and rs10993994 (10q11) (p = 1.40 × 10⁻³) in Mixed Ancestry men and rs10993994 (p = 1.56 × 10⁻⁹) in White men. No significant associations were observed for the analyses stratified by disease aggressiveness in the individual and the combined population group analysis. The present study demonstrates that a number of known PCa susceptibility variants may contribute to disease susceptibility in South African men. Larger genetic investigations extended to other South African population groups are warranted to confirm the role of these and other SNPs in disease susceptibility.

Epigenetic silencing of Oct4 by a complex containing SUV39H1 and Oct4 pseudogene lncRNA

Pseudogene-derived, long non-coding RNAs (lncRNAs) act as epigenetic regulators of gene expression. Here we present a panel of new mouse Oct4 pseudogenes and demonstrate that the X-linked Oct4 pseudogene Oct4P4 critically impacts mouse embryonic stem cells (mESCs) self-renewal. Sense Oct4P4 transcription produces a spliced, nuclear-restricted lncRNA that is efficiently upregulated during mESC differentiation. Oct4P4 lncRNA forms a complex with the SUV39H1 HMTase to direct the imposition of H3K9me3 and HP1α to the promoter of the ancestral Oct4 gene, located on chromosome 17, leading to gene silencing and reduced mESC self-renewal. Targeting Oct4P4 expression in primary mouse embryonic fibroblasts causes the re-acquisition of self-renewing features of mESC. We demonstrate that Oct4P4 lncRNA plays an important role in inducing and maintaining silencing of the ancestral Oct4 gene in differentiating mESCs. Our data introduces a sense pseudogene–lncRNA-based mechanism of epigenetic gene regulation that controls the cross-talk between pseudogenes and their ancestral genes.

Pseudogene derived long non-coding RNAs (lncRNAs) can regulate the expression of their ancestral genes. Here, the authors show that the Oct4 pseudogene OctP4 lncRNA plays an important role in inducing and maintaining silencing of the ancestral Oct4 gene in differentiating mouse embryonic stem cells.

Involvement of long non-coding RNA in colorectal cancer: From benchtop to bedside (Review)

Colorectal cancer (CRC) is one of the greatest threats to public health. Recent advances in whole-genome transcriptome analyses have enabled the identification of numerous members of a novel class of non-coding (nc)RNA, long ncRNA (lncRNA), which is broadly defined as RNA molecules that are >200 nt in length and lacking an open reading frame. In the present review, all lncRNAs associated with CRC are briefly summarized, with a particular focus on their potential roles as clinical biomarkers. CRC-associated lncRNAs involved in the underlying mechanisms of CRC progression are also initially included. This should benefit the development of novel markers and effective therapeutic targets for patients with CRC.

DNA methylation levels at chromosome 8q24 in peripheral blood are associated with 8q24 cancer susceptibility loci

Chromosome 8q24 has emerged as an important region for genetic susceptibility to various cancers, but little is known about the contribution of DNA methylation at 8q24. To evaluate variability in DNA methylation levels at 8q24 and the relationship with cancer susceptibility single nucleotide polymorphisms (SNPs) in this region, we quantified DNA methylation levels in peripheral blood at 145 CpG sites nearby 8q24 cancer susceptibility SNPs or MYC using pyrosequencing among 80 Caucasian men in the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial. For the 60 CpG sites meeting quality control, which also demonstrated temporal stability over a 5-year period, we calculated pairwise Spearman correlations for DNA methylation levels at each CpG site with 42 8q24 cancer susceptibility SNPs. To account for multiple testing, we adjusted P values into q values reflecting the false discovery rate (FDR). In contrast to the MYC CpG sites, most sites nearby the SNPs demonstrated good reproducibility, high methylation levels, and moderate-high between-individual variation. We observed 10 statistically significant (FDR < 0.05) CpG site-SNP correlations. These included correlations between an intergenic CpG site at Chr8:128393157 and the prostate cancer SNP rs16902094 (ρ = -0.54; P = 9.7 × 10(-7); q = 0.002), a PRNCR1 CpG site at Chr8:128167809 and the prostate cancer SNP rs1456315 (ρ = 0.52; P = 1.4 × 10(-6); q = 0.002), and two POU5F1B CpG sites and several prostate/colorectal cancer SNPs (for Chr8:128498051 and rs6983267, ρ = 0.46; P = 2.0 × 10(-5); q = 0.01). This is the first report of correlations between blood DNA methylation levels and cancer susceptibility SNPs at 8q24, suggesting that DNA methylation at this important susceptibility locus may contribute to cancer risk.

Clinical significance of POU5F1P1 rs10505477 polymorphism in Chinese gastric cancer patients receving cisplatin-based chemotherapy after surgical resection

This study aimed to investigate the association between POU class5 homeobox 1 pseudogene 1 gene (POU5F1P1) rs10505477 polymorphism and the prognosis of Chinese gastric cancer patients, who received cisplatin-based chemotherapy after surgical resection. POU5F1P1 rs10505477 was genotyped using the SNaPshot method in 944 gastric cancer patients who received gastrectomy. The association of rs10505477 G > A polymorphism with the progression and prognosis in gastric cancer patients was statistically analyzed using the SPSS version 18.0 for Windows. The results reveal that rs10505477 polymorphism has a negatively effect on the overall survival of gastric cancer patients in cisplatin-based chemotherapy subgroup (HR = 1.764, 95% CI = 1.069–2.911, p = 0.023). Our preliminary study indicates for the first time that POU5F1P1 rs10505477 is correlated with survival of gastric cancer patients who receving cisplatin-based chemotherapy after gastrectomy. Further studies are warranted to investigate the mechanism and to verify our results in different populations.

A genetic variant in pseudogene E2F3P1 contributes to prognosis of hepatocellular carcinoma

Certain pseudogenes may regulate their protein-coding cousins by competing for miRNAs and play an active biological role in cancer. However, few studies have focused on the association of genetic variations in pseudogenes with cancer prognosis. We selected six potentially functional single nucleotide polymorphisms (SNPs) in cancer-related pseudogenes, and performed a case-only study to assess the association between those SNPs and the prognosis of hepatocellular carcinoma (HCC) in 331 HBV-positive HCC patients without surgical treatment. Log-rank test and Cox proportional hazard models were used for survival analysis. We found that the A allele of rs9909601 in E2F3P1 was significantly associated with a better prognosis compared with the G allele [adjusted hazard ratio (HR)  =  0.69, 95% confidence interval (CI)  =  0.56–0.86, P  =  0.001]. Additionally, this protective effect was more predominant for patients without chemotherapy and transcatheter hepatic arterial chemoembolization (TACE) treatment. Interestingly, we also detected a statistically significant multiplicative interaction between genotypes of rs9909601 and chemotherapy or TACE status on HCC survival (P for multiplicative interaction < 0.001). These findings indicate that rs9909601 in the pseudogene E2F3P1 may be a genetic marker for HCC prognosis in Chinese.

An expressed retrogene of the master embryonic stem cell gene POU5F1 is associated with prostate cancer susceptibility

Genetic association studies of prostate and other cancers have identified a major risk locus at chromosome 8q24. Several independent risk variants at this locus alter transcriptional regulatory elements, but an affected gene and mechanism for cancer predisposition have remained elusive. The retrogene POU5F1B within the locus has a preserved open reading frame encoding a homolog of the master embryonic stem cell transcription factor Oct4. We find that 8q24 risk alleles are expression quantitative trait loci correlated with reduced expression of POU5F1B in prostate tissue and that predicted deleterious POU5F1B missense variants are also associated with risk of transformation. POU5F1 is known to be self-regulated by the encoded Oct4 transcription factor. We further observe that POU5F1 expression is directly correlated with POU5F1B expression. Our results suggest that a pathway critical to self-renewal of embryonic stem cells may also have a role in the origin of cancer.

OCT4 as a target of miR-34a stimulates p63 but inhibits p53 to promote human cell transformation

Human cell transformation is a key step for oncogenic development, which involves multiple pathways; however, the mechanism remains unclear. To test our hypothesis whether cell oncogenic transformation shares some mechanisms with the process of reprogramming non-stem cells to induced pluripotent stem cells (iPSC), we studied the relationship among the key factors for promoting or inhibiting iPSC in radiation-transformed human epithelial cell lines derived from different tissues (lung, breast and colon). We unexpectedly found that p63 and OCT4 were highly expressed (accompanied by low expressed p53 and miR-34a) in all transformed cell lines examined when compared with their non-transformed counterparts. We further elucidated the relationship of these factors: the 3p strand of miR-34a directly targeted OCT4 by binding to the 3′ untranslated region (3′-UTR) of OCT4 and, OCT4, in turn, stimulated p63 but inhibited p53 expression by binding to a specific region of the p63 or p53 promoter. Moreover, we revealed that the effects of OCT4 on promoting cell oncogenic transformation were by affecting p63 and p53. These results support that a positive loop exists in human cells: OCT4 upregulation as a consequence of inhibition of miR-34a, promotes p63 but suppresses p53 expression, which further stimulates OCT4 upregulation by downregulating miR-34a. This functional loop contributes significantly to cell transformation and, most likely, also to the iPSC process.

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

POU5F1B (POU domain class 5 transcription factor 1B), a processed pseudogene that is highly homologous to OCT4, was recently shown to be transcribed in cancer cells, but its clinical relevance and biological function have remained unclear. We now show that POU5F1B, which is located adjacent to MYC on human chromosome 8q24, is frequently amplified in gastric cancer (GC) cell lines. POU5F1B, but not OCT4, was also found to be expressed at a high level in GC cell lines and clinical specimens. In addition, the DNA copy number and mRNA abundance for POU5F1B showed a positive correlation in both cancer cell lines and GC specimens. Overexpression of POU5F1B in GC cells promoted colony formation in vitro as well as both tumorigenicity and tumor growth in vivo, and these effects were enhanced in the additional presence of MYC overexpression. Furthermore, knockdown of POU5F1B expression with a short hairpin RNA confirmed a role for the endogenous pseudogene in the promotion of cancer cell growth in vitro and tumor growth in vivo. POU5F1B overexpression induced upregulation of various growth factors in GC cells as well as exhibited mitogenic, angiogenic and antiapoptotic effects in GC xenografts. Finally, amplification of POU5F1B was detected in 17 (12%) of 145 cases of GC and was a significant predictor of poor prognosis in patients with stage IV disease. In conclusion, we found that the POU5F1B pseudogene is amplified and expressed at a high level in, as well as confers an aggressive phenotype on, GC, and that POU5F1B amplification is associated with a poor prognosis in GC patients.

Prognostication of OCT4 isoform expression in prostate cancer

Cancer stem cells (CSCs) refer to a subset of tumor cells that self-renew and affect tumor heterogeneity. This model has attracted considerable interest in recent years due to its implications in the prognosis and clinical management of cancer because CSCs mediate the occurrence, growth, and recurrence of tumors. OCT4 is central to embryonic stem cell self-renewal and differentiation into specific lineages and encodes two chief isoforms that are generated by alternative splicing--OCT4A and OCT4B. Their function in prostate cancer (PCa) is unknown. The prognostic function of OCT4 isoforms in PCa samples was examined by immunohistochemistry (IHC) and sensitivity and specificity of the antibodies used were evaluated by molecular biology techniques. Biochemical and pathological data and specimens from 193 patients with PCa were evaluated retrospectively. IHC, western blot, immunofluorescence, and automated image analysis were also performed. IHC was performed on a tissue microarray, and western blot and immunofluorescence were performed using the PCa cell line DU-145. IHC expression of OCT4 isoforms correlated with biochemical and pathological parameters, particularly biochemical recurrence-free survival (BCRFS). Patients with higher levels of OCT4B had lower Gleason scores and decreased likelihood of experiencing biochemical recurrence (BR). OCT4A(+) OCT4B(-) patients had the shortest BCRFS, and positivity for OCT4B expression was an independent prognostic factor for BCRFS in the multivariate analysis. We conclude that the expression of OCT4B is a strong marker of good prognosis, and its presence is associated with a decreased likelihood of BR. Thus, OCT4B might represent a powerful clinical prognostic biomarker for PCa patients.

Expressed pseudogenes in the transcriptional landscape of human cancers

Pseudogene transcripts can provide a novel tier of gene regulation through generation of endogenous siRNAs or miRNA-binding sites. Characterization of pseudogene expression, however, has remained confined to anecdotal observations due to analytical challenges posed by the extremely close sequence similarity with their counterpart coding genes. Here, we describe a systematic analysis of pseudogene "transcription" from an RNA-Seq resource of 293 samples, representing 13 cancer and normal tissue types, and observe a surprisingly prevalent, genome-wide expression of pseudogenes that could be categorized as ubiquitously expressed or lineage and/or cancer specific. Further, we explore disease subtype specificity and functions of selected expressed pseudogenes. Taken together, we provide evidence that transcribed pseudogenes are a significant contributor to the transcriptional landscape of cells and are positioned to play significant roles in cellular differentiation and cancer progression, especially in light of the recently described ceRNA networks. Our work provides a transcriptome resource that enables high-throughput analyses of pseudogene expression.

The 8q24 gene desert: an oasis of non-coding transcriptional activity

Understanding the functional effects of the wide-range of aberrant genetic characteristics associated with the human chromosome 8q24 region in cancer remains daunting due to the complexity of the locus. The most logical target for study remains the MYC proto-oncogene, a prominent resident of 8q24 that was first identified more than a quarter of a century ago. However, many of the amplifications, translocation breakpoints, and viral integration sites associated with 8q24 are often found throughout regions surrounding large expanses of the MYC locus that include other transcripts. In addition, chr.8q24 is host to a number of single nucleotide polymorphisms associated with cancer risk. Yet, the lack of a direct correlation between cancer risk alleles and MYC expression has also raised the possibility that MYC is not always the target of these genetic associations. The 8q24 region has been described as a "gene desert" because of the paucity of functionally annotated genes located within this region. Here we review the evidence for the role of other loci within the 8q24 region, most of which are non-coding transcripts, either in concert with MYC or independent of MYC, as possible candidate gene targets in malignancy.

RNA-Seq of human neurons derived from iPS cells reveals candidate long non-coding RNAs involved in neurogenesis and neuropsychiatric disorders

Genome-wide expression analysis using next generation sequencing (RNA-Seq) provides an opportunity for in-depth molecular profiling of fundamental biological processes, such as cellular differentiation and malignant transformation. Differentiating human neurons derived from induced pluripotent stem cells (iPSCs) provide an ideal system for RNA-Seq since defective neurogenesis caused by abnormalities in transcription factors, DNA methylation, and chromatin modifiers lie at the heart of some neuropsychiatric disorders. As a preliminary step towards applying next generation sequencing using neurons derived from patient-specific iPSCs, we have carried out an RNA-Seq analysis on control human neurons. Dramatic changes in the expression of coding genes, long non-coding RNAs (lncRNAs), pseudogenes, and splice isoforms were seen during the transition from pluripotent stem cells to early differentiating neurons. A number of genes that undergo radical changes in expression during this transition include candidates for schizophrenia (SZ), bipolar disorder (BD) and autism spectrum disorders (ASD) that function as transcription factors and chromatin modifiers, such as POU3F2 and ZNF804A, and genes coding for cell adhesion proteins implicated in these conditions including NRXN1 and NLGN1. In addition, a number of novel lncRNAs were found to undergo dramatic changes in expression, one of which is HOTAIRM1, a regulator of several HOXA genes during myelopoiesis. The increase we observed in differentiating neurons suggests a role in neurogenesis as well. Finally, several lncRNAs that map near SNPs associated with SZ in genome wide association studies also increase during neuronal differentiation, suggesting that these novel transcripts may be abnormally regulated in a subgroup of patients.

Large-scale exploration of gene-gene interactions in prostate cancer using a multistage genome-wide association study

Recent genome-wide association studies have identified independent susceptibility loci for prostate cancer that could influence risk through interaction with other, possibly undetected, susceptibility loci. We explored evidence of interaction between pairs of 13 known susceptibility loci and single nucleotide polymorphisms (SNP) across the genome to generate hypotheses about the functionality of prostate cancer susceptibility regions. We used data from Cancer Genetic Markers of Susceptibility: Stage I included 523,841 SNPs in 1,175 cases and 1,100 controls; Stage II included 27,383 SNPs in an additional 3,941 cases and 3,964 controls. Power calculations assessed the magnitude of interactions our study is likely to detect. Logistic regression was used with alternative methods that exploit constraints of gene-gene independence between unlinked loci to increase power. Our empirical evaluation demonstrated that an empirical Bayes (EB) technique is powerful and robust to possible violation of the independence assumption. Our EB analysis identified several noteworthy interacting SNP pairs, although none reached genome-wide significance. We highlight a Stage II interaction between the major prostate cancer susceptibility locus in the subregion of 8q24 that contains POU5F1B and an intronic SNP in the transcription factor EPAS1, which has potentially important functional implications for 8q24. Another noteworthy result involves interaction of a known prostate cancer susceptibility marker near the prostate protease genes KLK2 and KLK3 with an intronic SNP in PRXX2. Overall, the interactions we have identified merit follow-up study, particularly the EPAS1 interaction, which has implications not only in prostate cancer but also in other epithelial cancers that are associated with the 8q24 locus.