Genome-wide multi-omics profiling of the 8p11-p12 amplicon in breast carcinoma

Digital RNA sequencing using unique molecular identifiers enables ultrasensitive RNA mutation analysis

Mutation analysis is typically performed at the DNA level since most technical approaches are developed for DNA analysis. However, some applications, like transcriptional mutagenesis, RNA editing and gene expression analysis, require RNA analysis. Here, we combine reverse transcription and digital DNA sequencing to enable low error digital RNA sequencing. We evaluate yield, reproducibility, dynamic range and error correction rate for seven different reverse transcription conditions using multiplexed assays. The yield, reproducibility and error rate vary substantially between the specific conditions, where the yield differs 9.9-fold between the best and worst performing condition. Next, we show that error rates similar to DNA sequencing can be achieved for RNA using appropriate reverse transcription conditions, enabling detection of mutant allele frequencies <0.1% at RNA level. We also detect mutations at both DNA and RNA levels in tumor tissue using a breast cancer panel. Finally, we demonstrate that digital RNA sequencing can be applied to liquid biopsies, analyzing cell-free gene transcripts. In conclusion, we demonstrate that digital RNA sequencing is suitable for ultrasensitive RNA mutation analysis, enabling several basic research and clinical applications.

NSD3 in Cancer: Unraveling Methyltransferase-Dependent and Isoform-Specific Functions

NSD3 (nuclear receptor-binding SET domain protein 3) is a member of the NSD histone methyltransferase family of proteins. In recent years, it has been identified as a potential oncogene in certain types of cancer. The NSD3 gene encodes three isoforms, the long version (NSD3L), a short version (NSD3S) and the WHISTLE isoforms. Importantly, the NSD3S isoform corresponds to the N-terminal region of the full-length protein, lacking the methyltransferase domain. The chromosomal location of NSD3 is frequently amplified across cancer types, such as breast, lung, and colon, among others. Recently, this amplification has been correlated to a chromothripsis event, that could explain the different NSD3 alterations found in cancer. The fusion proteins containing NSD3 have also been reported in leukemia (NSD3-NUP98), and in NUT (nuclear protein of the testis) midline carcinoma (NSD3-NUT). Its role as an oncogene has been described by modulating different cancer pathways through its methyltransferase activity, or the short isoform of the protein, through protein interactions. Specifically, in this review we will focus on the functions that have been characterized as methyltransferase dependent, and those that have been correlated with the expression of the NSD3S isoform. There is evidence that both the NSD3L and NSD3S isoforms are relevant for cancer progression, establishing NSD3 as a therapeutic target. However, further functional studies are needed to differentiate NSD3 oncogenic activity as dependent or independent of the catalytic domain of the protein, as well as the contribution of each isoform and its clinical significance in cancer progression.

Pan-cancer analysis of genomic and transcriptomic data reveals the prognostic relevance of human proteasome genes in different cancer types

Background

The human proteasome gene family (PSM) consists of 49 genes that play a crucial role in cancer proteostasis. However, little is known about the effect of PSM gene expression and genetic alterations on clinical outcome in different cancer forms.

Methods

Here, we performed a comprehensive pan-cancer analysis of genetic alterations in PSM genes and the subsequent prognostic value of PSM expression using data from The Cancer Genome Atlas (TCGA) containing over 10,000 samples representing up to 33 different cancer types. External validation was performed using a breast cancer cohort and KM plotter with four cancer types.

Results

The PSM genetic alteration frequency was high in certain cancer types (e.g. 67%; esophageal adenocarcinoma), with DNA amplification being most common. Compared with normal tissue, most PSM genes were predominantly overexpressed in cancer. Survival analysis also established a relationship with PSM gene expression and adverse clinical outcome, where PSMA1 and PSMD11 expression were linked to more unfavorable prognosis in ≥ 30% of cancer types for both overall survival (OS) and relapse-free interval (PFI). Interestingly, PSMB5 gene expression was associated with OS (36%) and PFI (27%), and OS for PSMD2 (42%), especially when overexpressed.

Conclusion

These findings indicate that several PSM genes may potentially be prognostic biomarkers and novel therapeutic targets for different cancer forms.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-10079-4.

Pan-cancer analysis identifies BIRC5 as a prognostic biomarker

Background

The BIRC5 gene encodes for the Survivin protein, which is a member of the inhibitor of apoptosis family. Survivin is found in humans during fetal development, but generally not in adult cells thereafter. Previous studies have shown that Survivin is abundant in most cancer cells, thereby making it a promising target for anti-cancer drugs and a potential prognostic tool.

Methods

To assess genetic alterations and mutations in the BIRC5 gene as well as BIRC5 co-expression with other genes, genomic and transcriptomic data were downloaded via cBioPortal for approximately 9000 samples from The Cancer Genome Atlas (TCGA) representing 33 different cancer types and 11 pan-cancer organ systems, and validated using the ICGC Data Portal and COSMIC. TCGA BIRC5 RNA sequencing data from 33 different cancer types and matching normal tissue samples for 16 cancer types were downloaded from Broad GDAC Firehose and validated using breast cancer microarray data from our previous work and data sets from the GENT2 web-based tool. Survival data were analyzed with multivariable Cox proportional hazards regression analysis and validated using KM plotter for breast-, ovarian-, lung- and gastric cancer.

Results

Although genetic alterations in BIRC5 were not common in cancer, BIRC5 expression was significantly higher in cancer tissue compared to normal tissue in the 16 different cancer types. For 14/33 cancer types, higher BIRC5 expression was linked to worse overall survival (OS, 4/14 after adjusting for both age and tumor grade and 10/14 after adjusting only for age). Interestingly, higher BIRC5 expression was associated with better OS in lung squamous cell carcinoma and ovarian serous cystadenocarcinoma. Higher BIRC5 expression was also linked to shorter progressive-free interval (PFI) for 14/33 cancer types (4/14 after adjusting for both age and tumor grade and 10/14 after adjusting only for age). External validation showed that high BIRC5 expression was significantly associated with worse OS for breast-, lung-, and gastric cancer.

Conclusions

Our findings suggest that BIRC5 overexpression is associated with the initiation and progression of several cancer types, and thereby a promising prognostic biomarker.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09371-0.

Regulation of mRNA Translation by Hormone Receptors in Breast and Prostate Cancer

Breast and prostate cancer are the second and third leading causes of death amongst all cancer types, respectively. Pathogenesis of these malignancies is characterised by dysregulation of sex hormone signalling pathways, mediated by the estrogen receptor-α (ER) in breast cancer and androgen receptor (AR) in prostate cancer. ER and AR are transcription factors whose aberrant function drives oncogenic transcriptional programs to promote cancer growth and progression. While ER/AR are known to stimulate cell growth and survival by modulating gene transcription, emerging findings indicate that their effects in neoplasia are also mediated by dysregulation of protein synthesis (i.e., mRNA translation). This suggests that ER/AR can coordinately perturb both transcriptional and translational programs, resulting in the establishment of proteomes that promote malignancy. In this review, we will discuss relatively understudied aspects of ER and AR activity in regulating protein synthesis as well as the potential of targeting mRNA translation in breast and prostate cancer.

Breast cancers, mammary stem cells, and cancer stem cells, characteristics, and hypotheses

The cellular heterogeneity of breast cancers still represents a major therapeutic challenge. The latest genomic studies have classified breast cancers in distinct clusters to inform the therapeutic approaches and predict clinical outcomes. The mammary epithelium is composed of luminal and basal cells, and this seemingly hierarchical organization is dependent on various stem cells and progenitors populating the mammary gland. Some cancer cells are conceptually similar to the stem cells as they can self-renew and generate bulk populations of nontumorigenic cells. Two models have been proposed to explain the cell of origin of breast cancer and involve either the reprogramming of differentiated mammary cells or the dysregulation of mammary stem cells or progenitors. Both hypotheses are not exclusive and imply the accumulation of independent mutational events. Cancer stem cells have been isolated from breast tumors and implicated in the development, metastasis, and recurrence of breast cancers. Recent advances in single-cell sequencing help deciphering the clonal evolution within each breast tumor. Still, few clinical trials have been focused on these specific cancer cell populations.

Amplification of 8p11.23 in cancers and the role of amplicon genes

Copy number alterations are widespread in cancer genomes and are part of the genomic instability underlying the pathogenesis of neoplastic diseases. Recurrent copy number alterations of specific chromosomal loci may result in gains of oncogenes or losses of tumor suppressor genes and become entrenched in the genomic framework of certain types of cancers. The locus at chromosome 8p11.23 presents recurrent amplifications most commonly in squamous lung carcinomas, breast cancers, squamous esophageal carcinomas, and urothelial carcinomas. Amplification is rare in other cancers. The amplified segment involves several described oncogenes that may promote cancer cell survival and proliferation, as well as less well characterized genes that could also contribute to neoplastic processes. Genes proposed to be "drivers" in 8p11.23 amplifications include ZNF703, FGFR1 and PLPP5. Additional genes in the locus that could be functionally important in neoplastic networks include co-chaperone BAG4, lysine methyltransferase NSD3, ASH2L, a member of another methyltransferase complex, MLL and the mRNA processing and translation regulators LSM1 and EIF4EBP1. In this paper, genes located in the amplified segment of 8p11.23 will be examined for their role in cancer and data arguing for their importance for cancers with the amplification will be presented.

Inflammatory Breast Cancer: Clinical Implications of Genomic Alterations and Mutational Profiling

Simple Summary

Inflammatory breast cancer (IBC) is an aggressive disease with high mortality rates. Nowadays, there is no targeted treatment for this tumor type. Based on this context, we investigated the molecular profile of this disease by using well-established methodologies (high-resolution microarray platform, targeted next-generation sequencing, and immunohistochemistry) that have proven potential to unveil cancer biomarkers. We found alterations related to IBC aggressiveness and metastasis (gains of MDM4, losses of CHL1, and high homologous recombination deficiency scores), and worse overall survival (variants in HR and mismatch repair genes). We also compared the mutational profiling of our cases with literature data, which includes both non-IBC and IBC cases, validating our findings. Overall, we describe genetic alterations with the potential to be used as prognostic or predictive biomarkers and ultimately improve IBC patients’ care.

Abstract

Inflammatory breast cancer (IBC) is a rare and aggressive type of breast cancer whose molecular basis is poorly understood. We performed a comprehensive molecular analysis of 24 IBC biopsies naïve of treatment, using a high-resolution microarray platform and targeted next-generation sequencing (105 cancer-related genes). The genes more frequently affected by gains were MYC (75%) and MDM4 (71%), while frequent losses encompassed TP53 (71%) and RB1 (58%). Increased MYC and MDM4 protein expression levels were detected in 18 cases. These genes have been related to IBC aggressiveness, and MDM4 is a potential therapeutic target in IBC. Functional enrichment analysis revealed genes associated with inflammatory regulation and immune response. High homologous recombination (HR) deficiency scores were detected in triple-negative and metastatic IBC cases. A high telomeric allelic imbalance score was found in patients having worse overall survival (OS). The mutational profiling was compared with non-IBC (TCGA, n = 250) and IBC (n = 118) from four datasets, validating our findings. Higher frequency of TP53 and BRCA2 variants were detected compared to non-IBC, while PIKC3A showed similar frequency. Variants in mismatch repair and HR genes were associated with worse OS. Our study provided a framework for improved diagnosis and therapeutic alternatives for this aggressive tumor type.

Tumorigenic effects of TLX overexpression in HEK 293T cells

BACKGROUND

The human orphan receptor TLX (NR2E1) is a key regulator of neurogenesis, adult stem cell maintenance, and tumorigenesis. However, little is known about the genetic and transcriptomic events that occur following TLX overexpression in human cell lines.

AIMS

Here, we used cytogenetics and RNA sequencing to investigate the effect of TLX overexpression with an inducible vector system in the HEK 293T cell line.

METHODS AND RESULTS

Conventional spectral karyotyping was used to identify chromosomal abnormalities, followed by fluorescence in situ hybridization (FISH) analysis on chromosome spreads to assess TLX DNA copy number. Illumina paired-end whole transcriptome sequencing was then performed to characterize recurrent genetic variants (single nucleotide polymorphisms (SNPs) and indels), expressed gene fusions, and gene expression profiles. Lastly, flow cytometry was used to analyze cell cycle distribution. Intriguingly, we show that upon transfection with a vector containing the human TLX gene (eGFP-hTLX), an isochromosome forms on the long arm of chromosome 6, thereby resulting in DNA gain of the TLX locus (6q21) and upregulation of TLX. Induction of the eGFP-hTLX vector further increased TLX expression levels, leading to G0-G1 cell cycle arrest, genetic aberrations, modulation of gene expression patterns, and crosstalk with other nuclear receptors (AR, ESR1, ESR2, NR1H4, and NR3C2). We identified a 49-gene signature associated with central nervous system (CNS) development and carcinogenesis, in addition to potentially cancer-driving gene fusions (LARP1-CNOT8 and NSL1-ZDBF2) and deleterious genetic variants (frameshift insertions in the CTSH, DBF4, POSTN, and WDR78 genes).

CONCLUSION

Taken together, these findings illustrate that TLX may play a pivotal role in tumorigenesis via genomic instability and perturbation of cancer-related processes.

Transcriptome analysis reveals overlap in fusion genes in a phase I clinical cohort of TNBC and HGSOC patients treated with buparlisib and olaparib

PURPOSE

Fusion genes can be therapeutically relevant if they result in constitutive activation of oncogenes or repression of tumor suppressors. However, the prevalence and role of fusion genes in female cancers remain largely unexplored. Here, we investigate the fusion gene landscape in triple-negative breast cancer (TNBC) and high-grade serous ovarian cancer (HGSOC), two subtypes of female cancers with high molecular similarity but limited treatment options at present.

METHODS

RNA-seq was utilized to identify fusion genes in a cohort of 18 TNBC and HGSOC patients treated with the PI3K inhibitor buparlisib and the PARP inhibitor olaparib in a phase I clinical trial (NCT01623349). Differential gene expression analysis was performed to assess the function of fusion genes in silico. Finally, these findings were correlated with the reported clinical outcomes.

RESULTS

A total of 156 fusion genes was detected, whereof 44/156 (28%) events occurred in more than one patient. Low recurrence across samples indicated that the majority of fusion genes were private passenger events. The long non-coding RNA MALAT1 was involved in 97/156 (62%) fusion genes, followed in prevalence by MUC16, FOXP1, WWOX and XIST. Gene expression of FOXP1 was significantly elevated in patients with vs. without FOXP1 fusion (P= 0.02). From a clinical perspective, FOXP1 fusions were associated with a favorable overall survival.

CONCLUSIONS

In summary, this study provides the first characterization of fusion genes in a cohort of TNBC and HGSOC patients. An improved mechanistic understanding of fusion genes will support the future identification of innovative therapeutic approaches for these challenging diseases.

Genome-wide DNA copy number profiling and bioinformatics analysis of ovarian cancer reveals key genes and pathways associated with distinct invasive/migratory capabilities

Ovarian cancer (OC) metastasis presents major hurdles that must be overcome to improve patient outcomes. Recent studies have demonstrated copy number variations (CNVs) frequently contribute to alterations in oncogenic drivers. The present study used a CytoScan HD Array to analyse CNVs and loss of heterozygosity (LOH) in the entire genomes of 6 OC patients and human OC cell lines to determine the genetic target events leading to the distinct invasive/migratory capacities of OC. The results showed that LOH at Xq11.1 and Xp21.1 and gains at 8q21.13 were novel, specific CNVs. Ovarian cancer-related CNVs were then screened by bioinformatics analysis. In addition, transcription factors-target gene interactions were predicted with information from PASTAA analysis. As a result, six genes (i.e., GAB2, AKT1, EGFR, COL6A3, UGT1A1 and UGT1A8) were identified as strong candidates by integrating the above data with gene expression and clinical outcome data. In the transcriptional regulatory network, 4 known cancer-related transcription factors (TFs) interacted with 6 CNV-driven genes. The protein/DNA arrays revealed 3 of these 4 TFs as potential candidate gene-related transcription factors in OC. We then demonstrated that these six genes can serve as potential biomarkers for OC. Further studies are required to elucidate the pathogenesis of OC.

The Identification and Analysis of mRNA-lncRNA-miRNA Cliques From the Integrative Network of Ovarian Cancer

Ovarian cancer is one of the leading causes of cancer mortality in women. Since little clinical symptoms were shown in the early period of ovarian cancer, most patients were found in phases III-IV or with abdominal metastasis when diagnosed. The lack of effective early diagnosis biomarkers makes ovarian cancer difficult to screen. However, in essence, the fundamental problem is we know very little about the regulatory mechanisms during tumorigenesis of ovarian cancer. There are emerging regulatory factors, such as long noncoding RNAs (lncRNAs) and microRNAs (miRNAs), which have played important roles in cancers. Therefore, we analyzed the RNA-seq profiles of 407 ovarian cancer patients. An integrative network of 20,424 coding RNAs (mRNAs), 10,412 lncRNAs, and 742 miRNAs were construed with variance inflation factor (VIF) regression method. The mRNA-lncRNA-miRNA cliques were identified from the network and analyzed. Such promising cliques showed significant correlations with survival and stage of ovarian cancer and characterized the complex sponge regulatory mechanism, suggesting their contributions to tumorigenicity. Our results provided novel insights of the regulatory mechanisms among mRNAs, lncRNAs, and miRNAs and highlighted several promising regulators for ovarian cancer detection and treatment.

Eukaryotic initiation factor 4E-binding protein as an oncogene in breast cancer

BACKGROUND

Eukaryotic Initiation Factor 4E-Binding Protein (EIF4EBP1, 4EBP1) is overexpressed in many human cancers including breast cancer, yet the role of 4EBP1 in breast cancer remains understudied. Despite the known role of 4EBP1 as a negative regulator of cap-dependent protein translation, 4EBP1 is predicted to be an essential driving oncogene in many cancer cell lines in vitro, and can act as a driver of cancer cell proliferation. EIF4EBP1 is located within the 8p11-p12 genomic locus, which is frequently amplified in breast cancer and is known to predict poor prognosis and resistance to endocrine therapy.

METHODS

Here we evaluated the effect of 4EBP1 targeting using shRNA knock-down of expression of 4EBP1, as well as response to the mTORC targeted drug everolimus in cell lines representing different breast cancer subtypes, including breast cancer cells with the 8p11-p12 amplicon, to better define a context and mechanism for oncogenic 4EBP1.

RESULTS

Using a genome-scale shRNA screen on the SUM panel of breast cancer cell lines, we found 4EBP1 to be a strong hit in the 8p11 amplified SUM-44 cells, which have amplification and overexpression of 4EBP1. We then found that knock-down of 4EBP1 resulted in dramatic reductions in cell proliferation in 8p11 amplified breast cancer cells as well as in other luminal breast cancer cell lines, but had little or no effect on the proliferation of immortalized but non-tumorigenic human mammary epithelial cells. Kaplan-Meier analysis of EIF4EBP1 expression in breast cancer patients demonstrated that overexpression of this gene was associated with reduced relapse free patient survival across all breast tumor subtypes.

CONCLUSIONS

These results are consistent with an oncogenic role of 4EBP1 in luminal breast cancer and suggests a role for this protein in cell proliferation distinct from its more well-known role as a regulator of cap-dependent translation.

Transcriptomic and genomic profiling of early-stage ovarian carcinomas associated with histotype and overall survival

Ovarian cancer is the most lethal gynecological malignancy in the western world. Despite recent efforts to characterize ovarian cancer using molecular profiling, few targeted treatment options are currently available. Here, we examined genetic variants, fusion transcripts, SNP genotyping, and gene expression patterns for early-stage (I and II) ovarian carcinomas (n=96) in relation to clinicopathological characteristics and clinical outcome, thereby identifying novel genetic features of ovarian carcinomas. Furthermore, mutation frequencies of specific genetic variants and/or their gene expression patterns were associated with histotype and overall survival, e.g. SLC28A2 (mucinous ovarian carcinoma histotype), ARCN1 (low expression in 0-2 year survival group), and tumor suppressor MTUS1 (mutation status and overall survival). The long non-coding RNA MALAT1 was identified as a highly promiscuous fusion transcript in ovarian carcinoma. Moreover, gene expression deregulation for 23 genes was associated with tumor aggressiveness. Taken together, the novel biomarkers identified here may improve ovarian carcinoma subclassification and patient stratification according to histotype and overall survival.

Mutational game changer: Chromothripsis and its emerging relevance to cancer

In recent years, the paradigm that genomic abnormalities in cancer cells arise through progressive accumulation of mutational events has been challenged by the discovery of single catastrophic events. One such phenomenon termed chromothripsis, involving massive chromosomal rearrangements arising all at once, has emerged as a major mutational game changer. The strong interest in this process stems from its widespread association with a range of cancer types and its potential as a mutational driver. In this review, we first describe chromothripsis detection and incidence in cancers. We then explore recently proposed underlying mechanistic origins, which explain the curious observations of the highly localised nature of the rearrangements on chromothriptic chromosomes. Detection of chromothriptic patterns following incorporation of single chromosomes into micronuclei or following telomere attrition have greatly contributed to our understanding of the reasons behind this chromosomal restriction. These underlying cellular events have been found to be participants in the tumourigenic process, strongly suggesting a potential role for chromothripsis in cancer development. Thus, we discuss potential implications of chromothripsis for cancer progression and therapy.