Copy number loss or silencing of apoptosis-effector genes in cancer

Neuroblastoma and Glioblastoma Cases With Amplified Oncogenes Have Reduced Numbers of Tumor-Resident Adaptive Immune Receptor Recombinations

PURPOSE

In certain cancers, oncogene amplification is correlated with an immunologically cold or noninflamed, tumor immune microenvironment (TIME) and a worse prognosis, for example, in the case of MYCN-amplified neuroblastoma (NBL). However, for other cancer types, the relationship between oncogene amplification and immune response is more complicated or unresolved. One such cancer is glioblastoma multiforme (GBM), in which the epidermal growth factor receptor (EGFR) oncogene is commonly amplified. Unlike MYCN-amplified NBL, EGFR-amplified GBM has not been shown to correlate with a distinct survival probability.

METHODS

Given this contrasting state for NBL and GBM, we sought to apply a genomics approach to evaluating the immune response for cases with gene amplification.

RESULTS

Our results confirmed and added further specificity to the cold TIME of MYCN-amplified NBL. Moreover, we demonstrated a novel state of immunologically cold EGFR-amplified GBM tumors.

CONCLUSION

This approach to using copy number variation and immune receptor recombination read recovery levels to assess gene amplification and TIME, respectively, may be particularly efficient for the rapid evaluation of many other cancer types.

CNV assessments associated with outcome distinctions for adult and pediatric cancers: Loss of BRCA1 in neuroblastoma associates with a lower survival probability

DNA copy number variations (CNV) are common in cancer development, however, CNV detection approaches that include assessments of small CNVs, for example, due to locally misaligned sister chromatid exchanges, have not been substantially applied. Using such approaches, CNVs have been detected, in the cancer setting, for regulatory elements common to both proliferation and apoptosis effector genes, but no linkage has yet been made to cancer patient clinical data. Thus, we hypothesized that copy number losses, including local copy number losses, of specific apoptosis effector genes would be associated with reduced survival. Both whole genome and whole exome files were processed for validations and consistency. Results indicated lower late-stage survival for multiple myeloma cases representing reduced BAD and CASP3 copies, as well as for lung adenocarcinoma cases representing reduced BAX and CASP3 copies. Results also indicated that neuroblastoma (NBL) cases representing reduced copies of CASP9 and BRCA1 had reduced overall survival probabilities, with the BRCA1 results being particularly notable due to previous reports of BRCA1 inactivating mutations in NBL. Overall, novel approaches to assessing CNVs offers the promise of establishing patient risk stratifications and of identifying single genes or other small spaces in the genome where a CNV may be linked to specific outcomes.

An age-based, RNA expression paradigm for survival biomarker identification for pediatric neuroblastoma and acute lymphoblastic leukemia

Background

Pediatric cancer survival rates overall have been improving, but neuroblastoma (NBL) and acute lymphoblastic leukemia (ALL), two of the more prevalent pediatric cancers, remain particularly challenging. One issue not yet fully addressed is distinctions attributable to age of diagnosis.

Methods

In this report, we verified a survival difference based on diagnostic age for both pediatric NBL and pediatric ALL datasets, with younger patients surviving longer for both diseases. We identified several gene expression markers that correlated with age, along a continuum, and then used a series of age-independent survival metrics to filter these initial correlations.

Results

For pediatric NBL, we identified 2 genes that are expressed at a higher level in lower surviving patients with an older diagnostic age; and 4 genes that are expressed at a higher level in longer surviving patients with a younger diagnostic age. For pediatric ALL, we identified 3 genes expressed at a higher level in lower surviving patients with an older diagnostic age; and 17 genes expressed at a higher level in longer surviving patients with a younger diagnostic age.

Conclusions

This process implicated pan-chromosome effects for chromosomes 11 and 17 in NBL; and for the X chromosome in ALL.

Electronic supplementary material

The online version of this article (10.1186/s12935-019-0790-5) contains supplementary material, which is available to authorized users.

Elucidating feed-forward apoptosis signatures in breast cancer datasets: Higher FOS expression associated with a better outcome

Overstimulation of pro-proliferative pathways and high level expression of pro-proliferative transcription factors (TFs) can lead to apoptosis. This is likely due to TF binding sites for pro-proliferative TFs common to pro-proliferative and pro-apoptosis-effector genes. Certain clinical datasets have indicated that molecular markers associated with higher proliferation rates lead to improved outcomes for patients with cancer. These observations have been extensively assessed on a general basis, however there has been little work dissecting feed-forward apoptosis signaling pathways that may represent specific distinctions between a pro-proliferative mechanism and a pro-apoptotic mechanism in samples from patients with cancer. Using The Cancer Genome Atlas datasets and bioinformatic approaches, the present study reports that higher FOS expression levels, along with higher FOS target apoptosis-effector gene expression, is associated with an increased survival, while higher POU2F1 expression is associated with a reduced survival (average difference of 25.9 months survival). In summary, in the datasets examined FOS represents an apoptosis-driver and high POU2F1 represents a driver mechanism for cancer development.

RICTOR amplification identifies a subgroup in small cell lung cancer and predicts response to drugs targeting mTOR

Small cell lung cancer (SCLC) is an aggressive cancer that represents ~15% of all lung cancers. Currently there are no targeted therapies to treat SCLC. Our genomic analysis of a metastatic SCLC cohort identified recurrent RICTOR amplification. Here, we examine the translational potential of this observation. RICTOR was the most frequently amplified gene observed (~14% patients), and co-amplified with FGF10 and IL7R on chromosome 5p13. RICTOR copy number variation correlated with RICTOR protein expression in SCLC cells. In parallel, cells with RICTOR copy number (CN) gain showed increased sensitivity to three mTOR inhibitors, AZD8055, AZD2014 and INK128 in cell growth assays, with AZD2014 demonstrating the best inhibition of downstream signaling. SCLC cells with RICTOR CN gain also migrated more rapidly in chemotaxis and scratch wound assays and were again more sensitive to mTOR inhibitors. The overall survival in SCLC patients with RICTOR amplification was significantly decreased (p = 0.021). Taken together, our results suggest that SCLC patients with RICTOR amplification may constitute a clinically important subgroup because of their potential response to mTORC1/2 inhibitors.

Identification of specific feed-forward apoptosis mechanisms and associated higher survival rates for low grade glioma and lung squamous cell carcinoma

The mechanisms of cell proliferation due to the overexpression of certain transcription factors (TFs) have been well documented in the cancer setting. However, many of these same TFs have pro-apoptotic effects, particularly when expressed or activated at high levels, a process referred to as feed-forward apoptosis (FFA). To determine whether cancers could be stratified on the basis of specific FFA signatures, RNASeq data representing samples from the cancer genome atlas were analyzed, revealing that high expression of the pro-proliferative TFs, MYC and YY1, is associated with a favorable outcome in low-grade glioma (LGG) and lung squamous cell carcinoma (LUSC), respectively. Analysis of the RNASeq data also led to the identification of specific apoptosis-effector genes whose expression levels correlate with increased survival rates, for both LGG and LUSC. Although FFA has been demonstrated as a general effect in cancer, in this report, for the first time, results identify specific TFs and their responsive effector genes that distinguish subsets of cancer samples undergoing more or less of a FFA process in a way that is associated with distinct patient survival rates.

Impact of SNPs on CpG Islands in the MYC and HRAS oncogenes and in a wide variety of tumor suppressor genes: A multi-cancer approach

Single nucleotide polymorphisms (SNPs) that occur within CpG Islands may lead to increased hypermethylation if a SNP allele has the potential to form a CpG dinucleotide, as well as potentially lead to hypomethylation if a SNP allele eliminates a CpG dinucleotide. We analyzed CpG-related SNP allele frequencies in whole genome sequences (WGS) across 5 TCGA cancer datasets, thereby exploiting a more recent appreciation for signaling pathway degeneracy in cancer. The cancer data sets were analyzed for SNPs in CpG islands associated with the oncogenes, HRAS and MYC, and in the CpG islands associated with the tumor suppressor genes, APC, DCC, and RB1. We determined that one SNP allele (rs3824120) in a CpG island associated with MYC which eliminated a CpG was more common in the cancer datasets than in the 100Genomes databases (p < 0.01). For HRAS, 2 SNP alleles (rs112690925, rs7939028) that created CpG's occurred significantly less frequently in the cancer data sets than in the general SNP databases (e.g., rs7939028, p < 0.0002, in comparison with AllSNPs(142)). Also, one SNP allele (rs4940177) that created a CpG in a CpG island associated with the DCC tumor suppressor gene, was more common in the cancer datasets (p < 0.0007). To understand a broader picture of the potential of SNP alleles to create CpG's in CpG islands of tumor suppressor genes, we developed a scripted algorithm to assess the SNP alleles associated with the CpG islands of 43 tumor suppressor genes. The following tumor suppressor genes have the possibility of significant, percent increases in their CpG counts, depending on which SNP allele(s) is present: VHL, BRCA1, BRCA2, CHEK2, PTEN and RB1.

Stratifying melanoma and breast cancer TCGA datasets on the basis of the CNV of transcription factor binding sites common to proliferation- and apoptosis-effector genes

Transcription factors that activate both proliferation- and apoptosis-effector genes, along with a number of related observations, have led to a proposal for a feed forward mechanism of activating the two gene classes, whereby a certain concentration of a transcription factor activates the proliferation-effector genes and a higher concentration of the transcription factor activates the apoptosis-effector genes. We reasoned that this paradigm of regulation could lead to, in the cancer setting, a selection for relatively reduced copy numbers of apoptosis-effector gene, transcription factor binding sites (TFBS). Thus, the aim of this investigation was to examine the DNA sequencing read depths of TFBS for a set of proliferation- and apoptosis-effector genes, normalized to the read depths found in matching blood samples, as provided by the cancer genome atlas (TCGA); and thereby document copy number differences among these TFBS. We determined that the melanoma and breast cancer, TCGA datasets could be divided into three categories: (i) no detectable copy number variation for the proliferation- and apoptosis-effector, shared TFBS; (ii) a relative increase in the copy number of proliferation-effector gene TFBS, compared with the copy number of the apoptosis-effector gene TFBS; and (iii) a relative decrease in the number of proliferation-effector gene TFBS. Thus, we conclude that changes in the relative copies of the shared TFBS, for proliferation- and apoptosis-effector genes, have the potential of impacting tumor cell proliferative and apoptotic capacities.

TCGA: Increased oncoprotein coding region mutations correlate with a greater expression of apoptosis-effector genes and a positive outcome for stomach adenocarcinoma

Oncogene mutations are primarily thought to facilitate uncontrolled cell growth. However, overexpression of oncoproteins likely leads to apoptosis in a feed forward mechanism, whereby a certain level of oncoprotein leads to the activation of pro-proliferation effector genes and higher levels lead to activation of pro-apoptotic effector genes. TCGA STAD barcodes having no oncoprotein coding region mutations represented reduced expression of the apoptosis-effector genes compared with barcodes with multiple oncoprotein coding region mutations. Furthermore, STAD barcodes in a "no-subsequent tumor" group, representing 224 samples, and in a "positive outcome" group, had more oncoprotein coding regions mutated, on average, than barcodes of the new tumor and negative outcome groups, respectively. BRAF, CTNNB1, KRAS and MTOR coding region mutations (as a group) had the strongest association with the no-subsequent tumor group. Tumor suppressor coding region mutations were also correlated with no-subsequent tumor. These results are consistent with an oncoprotein-mediated, feed-forward mechanism of apoptosis in patients. Importantly, the no-subsequent tumor group also had more overall mutations. This result leads to considerations of unhealthy cells or cells with more neo-antigens for immune rejection. However, a probabilistic aspect of mutagenesis is also consistent with more oncoprotein and tumor suppressor protein mutations, in cases of more overall mutations, and thus a higher likelihood of activation of feed forward apoptosis pathways.

Gastric cancer and gene copy number variation: emerging cancer drivers for targeted therapy

Gastric cancer (GC) is among the most common malignancy in the world with poor prognosis and limited treatment options. It has been established that gastric carcinogenesis is caused by a complex interaction between host and environmental factors. Copy number variation (CNV) refers to a form of genomic structural variation that results in abnormal gene copy numbers, including gene amplification, gain, loss and deletion. DNA CNV is an important influential factor for the expression of both protein-coding and non-coding genes, affecting the activity of various signaling pathways. CNV arises as a result of preferential selection that favors cancer development, and thus, targeting the amplified 'driver genes' in GC may provide novel opportunities for personalized therapy. The detection of CNVs in chromosomal or mitochondrial DNA from tissue or blood samples may assist the diagnosis, prognosis and targeted therapy of GC. In this review, we discuss the recent CNV discoveries that shed light on the molecular pathogenesis of GC, with a specific emphasis on CNVs that display diagnostic, prognostic or therapeutic significances in GC.

Tumor suppressor genes are larger than apoptosis-effector genes and have more regions of active chromatin: Connection to a stochastic paradigm for sequential gene expression programs

Apoptosis- and proliferation-effector genes are substantially regulated by the same transactivators, with E2F-1 and Oct-1 being notable examples. The larger proliferation-effector genes have more binding sites for the transactivators that regulate both sets of genes, and proliferation-effector genes have more regions of active chromatin, i.e, DNase I hypersensitive and histone 3, lysine-4 trimethylation sites. Thus, the size differences between the 2 classes of genes suggest a transcriptional regulation paradigm whereby the accumulation of transcription factors that regulate both sets of genes, merely as an aspect of stochastic behavior, accumulate first on the larger proliferation-effector gene "traps," and then accumulate on the apoptosis effector genes, thereby effecting sequential activation of the 2 different gene sets. As IRF-1 and p53 levels increase, tumor suppressor proteins are first activated, followed by the activation of apoptosis-effector genes, for example during S-phase pausing for DNA repair. Tumor suppressor genes are larger than apoptosis-effector genes and have more IRF-1 and p53 binding sites, thereby likewise suggesting a paradigm for transcription sequencing based on stochastic interactions of transcription factors with different gene classes. In this report, using the ENCODE database, we determined that tumor suppressor genes have a greater number of open chromatin regions and histone 3 lysine-4 trimethylation sites, consistent with the idea that a larger gene size can facilitate earlier transcriptional activation via the inclusion of more transactivator binding sites.