scCircle-seq unveils the diversity and complexity of extrachromosomal circular DNAs in single cells

Extrachromosomal circular DNAs (eccDNAs) have emerged as important intra-cellular mobile genetic elements that affect gene copy number and exert in trans regulatory roles within the cell nucleus. Here, we describe scCircle-seq, a method for profiling eccDNAs and unraveling their diversity and complexity in single cells. We implement and validate scCircle-seq in normal and cancer cell lines, demonstrating that most eccDNAs vary largely between cells and are stochastically inherited during cell division, although their genomic landscape is cell type-specific and can be used to accurately cluster cells of the same origin. eccDNAs are preferentially produced from chromatin regions enriched in H3K9me3 and H3K27me3 histone marks and are induced during replication stress conditions. Concomitant sequencing of eccDNAs and RNA from the same cell uncovers the absence of correlation between eccDNA copy number and gene expression levels, except for a few oncogenes, including MYC, contained within a large eccDNA in colorectal cancer cells. Lastly, we apply scCircle-seq to one prostate cancer and two breast cancer specimens, revealing cancer-specific eccDNA landscapes and a higher propensity of eccDNAs to form in amplified genomic regions. scCircle-seq is a scalable tool that can be used to dissect the complexity of eccDNAs across different cell and tissue types, and further expands the potential of eccDNAs for cancer diagnostics.

Genomic dissection and mutation-specific target discovery for breast cancer PIK3CA hotspot mutations

Background

Recent advancements in high-throughput genomics and targeted therapies have provided tremendous potential to identify and therapeutically target distinct mutations associated with cancers. However, to date the majority of targeted therapies are used to treat all functional mutations within the same gene, regardless of affected codon or phenotype.

Results

In this study, we developed a functional genomic analysis workflow with a unique isogenic cell line panel bearing two distinct hotspot PIK3CA mutations, E545K and H1047R, to accurately identify targetable differences between mutations within the same gene. We performed RNA-seq and ATAC-seq and identified distinct transcriptomic and epigenomic differences associated with each PIK3CA hotspot mutation. We used this data to curate a select CRISPR knock out screen to identify mutation-specific gene pathway vulnerabilities. These data revealed AREG as a E545K-preferential target that was further validated through in vitro analysis and publicly available patient databases.

Conclusions

Using our multi-modal genomics framework, we discover distinct differences in genomic regulation between PIK3CA hotspot mutations, suggesting the PIK3CA mutations have different regulatory effects on the function and downstream signaling of the PI3K complex. Our results demonstrate the potential to rapidly uncover mutation specific molecular targets, specifically AREG and a proximal gene regulatory region, that may provide clinically relevant therapeutic targets. The methods outlined provide investigators with an integrative strategy to identify mutation-specific targets for the treatment of other oncogenic mutations in an isogenic system.

Alternative Tissue Fixation Protocols Dramatically Reduce the Impact of DNA Artifacts, Unraveling the Interpretation of Clinical Comprehensive Genomic Profiling

Formalin-fixed paraffin-embedded (FFPE) samples represent the cornerstone of tissue-based analysis in precision medicine. Targeted next-generation sequencing panels are routinely used to analyze a limited number of genes to guide treatment decision-making for advanced-stage patients. The number and complexity of genetic alterations to be investigated are rapidly growing; in several instances, a comprehensive genomic profiling analysis is needed. The poor quality of genetic material extracted from FFPE samples may impact the feasibility/reliability of sequencing data. We sampled 9 colorectal cancers to allow 4 parallel fixations: (1) neutral buffered formalin (NBF), (2) acid-deprived formalin fixation (ADF), (3) precooled ADF (coldADF), and (4) glyoxal acid free (GAF). DNA extraction, fragmentation analysis, and sequencing by 2 large next-generation sequencing panels (OCAv3 and TSO500) followed. We comprehensively analyzed library and sequencing quality controls and the quality of sequencing results. Libraries from coldADF samples showed significantly longer reads than the others with both panels. ADF-derived and coldADF-derived libraries showed the lowest level of noise and the highest levels of uniformity with the OCAv3 panel, followed by GAF and NBF samples. The data uniformity was confirmed by the TSO500 results, which also highlighted the best performance in terms of the total region sequenced for the ADF and coldADF samples. NBF samples had a significantly smaller region sequenced and displayed a significantly lower number of evaluable microsatellite loci and a significant increase in single-nucleotide variations compared with other protocols. Mutational signature 1 (aging and FFPE artifact related) showed the highest (37%) and lowest (17%) values in the NBF and coldADF samples, respectively. Most of the identified genetic alterations were shared by all samples in each lesion. Five genes showed a different mutational status across samples and/or panels: 4 discordant results involved NBF samples. In conclusion, acid-deprived fixatives (GAF and ADF) guarantee the highest DNA preservation/sequencing performance, thus allowing more complex molecular profiling of tissue samples.

Fatty acid synthase as a new therapeutic target for HER2-positive gastric cancer

PURPOSE

Trastuzumab is an HER2-specific agent approved as the gold-standard therapy for advanced HER2-positive (HER2+) gastric cancer (GC), but the high rate and rapid appearance of resistance limit its clinical efficacy, resulting in the need to identify new vulnerabilities. Defining the drivers influencing HER2+ cancer stem cell (CSC) maintenance/survival could represent a clinically useful strategy to counteract tumor growth and therapy resistance. Accumulating evidence show that targeting crucial metabolic hubs, as the fatty acid synthase (FASN), may be clinically relevant.

METHODS

FASN protein and transcript expression were examined by WB and FACS and by qRT-PCR and GEP analyses, respectively, in trastuzumab-sensitive and trastuzumab-resistant HER2+ GC cell lines cultured in adherent (2D) or gastrosphere promoting (3D) conditions. Molecular data were analyzed in silico in public HER2+ GC datasets. The effectiveness of the FASN inhibitor TVB3166 to overcome anti-HER2 therapy resistance was tested in vitro in gastrospheres forming efficiency bioassays and in vivo in mice bearing trastuzumab-resistant GC cells.

RESULTS

We compared the transcriptome profiles of HER2+ GC cells cultured in 2D versus 3D conditions finding a significant enrichment of FASN in 3D cultures. FASN upregulation significantly correlated with high stemness score and poor prognosis in HER2+ GC cases. TVB3166 treatment significantly decreased GCSCs in all cell targets. HER2 and FASN cotargeting significantly decreased the capability to form gastrospheres versus monotherapy and reduced the in vivo growth of trastuzumab-resistant GC cells.

CONCLUSION

Our findings indicate that cotargeting HER2 and FASN increase the benefit of anti-HER2 therapy representing a new opportunity for metabolically combating trastuzumab-resistant HER2+ GC.

HER2 Copy Number and Resistance Mechanisms in Patients with HER2-positive Advanced Gastric Cancer Receiving Initial Trastuzumab-based Therapy in JACOB Trial

PURPOSE

In JACOB trial, pertuzumab added to trastuzumab-chemotherapy did not significantly improve survival of patients with HER2-positive metastatic gastric cancer, despite 3.3 months increase versus placebo. HER2 copy-number variation (CNV) and AMNESIA panel encompassing primary resistance alterations (KRAS/PIK3CA/MET mutations, KRAS/EGFR/MET amplifications) may improve patients' selection for HER2 inhibition.

EXPERIMENTAL DESIGN

In a post hoc analysis of JACOB on 327 samples successfully sequenced by next-generation sequencing (NGS; Oncomine Focus DNA), HER2 CNV, HER2 expression by IHC, and AMNESIA were correlated with overall response rate (ORR), progression-free survival (PFS), and overall survival (OS) by univariable/multivariable models.

RESULTS

Median HER2 CNV was 4.7 (interquartile range, 2.2-16.9). HER2 CNV-high versus low using the median as cutoff was associated with longer median PFS (10.5 vs. 6.4 months; HR = 0.48; 95% confidence interval: 0.38-0.62; P < 0.001) and OS (20.3 vs. 13.0 months; HR = 0.54; 0.42-0.72; P < 0.001). Combining HER2 CNV and IHC improved discriminative ability, with better outcomes restricted to HER2-high/HER2 3+ subgroup. AMNESIA positivity was found in 51 (16%), with unadjusted HR = 1.35 (0.98-1.86) for PFS; 1.43 (1.00-2.03) for OS.In multivariable models, only HER2 CNV status remained significant for PFS (P < 0.001) and OS (P = 0.004). Higher ORR was significantly associated with IHC 3+ [61% vs. 34% in 2+; OR = 3.11 (1.89-5.17)] and HER2-high [59% vs. 43% in HER2-low; OR = 1.84 (1.16-2.94)], with highest OR in the top CNV quartile. These biomarkers were not associated with treatment effect of pertuzumab.

CONCLUSIONS

HER2 CNV-high assessed by NGS may be associated with better ORR, PFS, and OS in a JACOB subgroup, especially if combined with HER2 3+. The negative prognostic role of AMNESIA requires further clinical validation.

The Tumor-Specific Expression of L1 Retrotransposons Independently Correlates with Time to Relapse in Hormone-Negative Breast Cancer Patients

BACKGROUND

Long-Interspersed Nuclear Element (L1) retrotransposons are silenced in healthy tissues but unrepressed in cancer. Even if L1 reactivation has been associated with reduced overall survival in breast cancer (BC) patients, a comprehensive correlation with clinicopathological features is still missing.

METHODS

Using quantitative, reverse-transcription PCR, we assessed L1 mRNA expression in 12 BC cells, 210 BC patients and in 47 normal mammary tissues. L1 expression was then correlated with molecular and clinicopathological data.

RESULTS

We identified a tumor-exclusive expression of L1s, absent in normal mammary cells and tissues. A positive correlation between L1 expression and tumor dedifferentiation, lymph-node involvement and increased immune infiltration was detected. Molecular subtyping highlighted an enrichment of L1s in basal-like cells and cancers. By exploring disease-free survival, we identified L1 overexpression as an independent biomarker for patients with a high risk of recurrence in hormone-receptor-negative BCs.

CONCLUSIONS

Overall, L1 reactivation identified BCs with aggressive features and patients with a worse clinical fate.

MiR-100 is a predictor of endocrine responsiveness and prognosis in patients with operable luminal breast cancer

PURPOSE

Overexpression of miR-100 in stem cells derived from basal-like breast cancers causes loss of stemness, induction of luminal breast cancer markers and response to endocrine therapy. We, therefore, explored miR-100 as a novel biomarker in patients with luminal breast cancer.

METHODS

miR-100 expression was studied in 90 patients with oestrogen-receptor-positive/human-epidermal growth factor receptor 2-negative breast cancer enrolled in a prospective study of endocrine therapy given either preoperatively, or for the treatment of de novo metastatic disease. Response was defined as a Ki67 ≤2.7% after 21±3 days of treatment. The prognostic role of miR-100 expression was evaluated in the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) and The Cancer Genome Atlas (TCGA) breast cancer datasets. Additionally, we explored the correlation between miR-100 and the expression its targets reported as being associated with endocrine resistance. Finally, we evaluated whether a signature based on miR-100 and its target genes could predict the luminal A molecular subtype.

RESULTS

Baseline miR-100 was significantly anticorrelated with baseline and post-treatment Ki67 (p<0.001 and 0.004, respectively), and independently associated with response to treatment (OR 3.329, p=0.047). In the METABRIC dataset, high expression of miR-100 identified women with luminal A tumours treated with adjuvant endocrine therapy with improved overall survival (HR 0.55, p<0.001). miR-100 was negatively correlated with PLK1, FOXA1, mTOR and IGF1R expression, potentially explaining its prognostic effect. Finally, a miR-100-based signature developed in patients enrolled in the prospective study outperformed Ki67 alone in predicting the luminal A phenotype.

CONCLUSIONS

Our findings suggest that miR-100 should be further explored as a biomarker in patients with luminal breast cancer.

Thyroid hormone inhibits hepatocellular carcinoma progression via induction of differentiation and metabolic reprogramming

BACKGROUND & AIMS

Only limited therapeutic options are currently available for hepatocellular carcinoma (HCC), making the development of effective alternatives essential. Based on the recent finding that systemic or local hypothyroidism is associated with HCC development in humans and rodents, we investigated whether the thyroid hormone triiodothyronine (T3) could inhibit the progression of HCCs.

METHODS

Different rat and mouse models of hepatocarcinogenesis were investigated. The effect of T3 on tumorigenesis and metabolism/differentiation was evaluated by transcriptomic analysis, quantitative reverse transcription PCR, immunohistochemistry, and enzymatic assay.

RESULTS

A short treatment with T3 caused a shift in the global expression profile of the most aggressive preneoplastic nodules towards that of normal liver. This genomic reprogramming preceded the disappearance of nodules and involved reprogramming of metabolic genes, as well as pro-differentiating transcription factors, including Kruppel-like factor 9, a target of the thyroid hormone receptor β (TRβ). Treatment of HCC-bearing rats with T3 strongly reduced the number and burden of HCCs. Reactivation of a local T3/TRβ axis, a switch from Warburg to oxidative metabolism and loss of markers of poorly differentiated hepatocytes accompanied the reduced burden of HCC. This effect persisted 1 month after T3 withdrawal, suggesting a long-lasting effect of the hormone. The antitumorigenic effect of T3 was further supported by its inhibitory activity on cell growth and the tumorigenic ability of human HCC cell lines.

CONCLUSIONS

Collectively, these findings suggest that reactivation of the T3/TRβ axis induces differentiation of neoplastic cells towards a more benign phenotype and that T3 or its analogs, particularly agonists of TRβ, could be useful tools in HCC therapy.

LAY SUMMARY

Hepatocellular carcinoma (HCC) represents an important challenge for global health. Recent findings showed that systemic or local hypothyroidism is associated with HCC development. In rat models, we showed that administration of the thyroid hormone T3 impaired HCC progression, even when given at late stages. This is relevant from a translational point of view as HCC is often diagnosed at an advanced stage when it is no longer amenable to curative treatments. Thyroid hormones and/or thyromimetics could be useful for the treatment of patients with HCC.

Clustered protocadherins methylation alterations in cancer

Background

Clustered protocadherins (PCDHs) map in tandem at human chromosome 5q31 and comprise three multi-genes clusters: α-, β- and γ-PCDH. The expression of this cluster consists of a complex mechanism involving DNA hub formation through DNA-CCTC binding factor (CTCF) interaction. Methylation alterations can affect this interaction, leading to transcriptional dysregulation. In cancer, clustered PCDHs undergo a mechanism of long-range epigenetic silencing by hypermethylation.

Results

In this study, we detected frequent methylation alterations at CpG islands associated to these clustered PCDHs in all the solid tumours analysed (colorectal, gastric and biliary tract cancers, pilocytic astrocytoma), but not hematologic neoplasms such as chronic lymphocytic leukemia. Importantly, several altered CpG islands were associated with CTCF binding sites. Interestingly, our analysis revealed a hypomethylation event in pilocytic astrocytoma, suggesting that in neuronal tissue, where PCDHs are highly expressed, these genes become hypomethylated in this type of cancer. On the other hand, in tissues where PCDHs are lowly expressed, these CpG islands are targeted by DNA methylation. In fact, PCDH-associated CpG islands resulted hypermethylated in gastrointestinal tumours.

Conclusions

Our study highlighted a strong alteration of the clustered PCDHs methylation pattern in the analysed solid cancers and suggested these methylation aberrations in the CpG islands associated with PCDH genes as powerful diagnostic biomarkers.

Electronic supplementary material

The online version of this article (10.1186/s13148-019-0695-0) contains supplementary material, which is available to authorized users.

Stromal contribution to the colorectal cancer transcriptome

Recent studies identified a poor-prognosis stem/serrated/mesenchymal (SSM) transcriptional subtype of colorectal cancer (CRC). We noted that genes upregulated in this subtype are also prominently expressed by stromal cells, suggesting that SSM transcripts could derive from stromal rather than epithelial cancer cells. To test this hypothesis, we analyzed CRC expression data from patient-derived xenografts, where mouse stroma supports human cancer cells. Species-specific expression analysis showed that the mRNA levels of SSM genes were mostly due to stromal expression. Transcriptional signatures built to specifically report the abundance of cancer-associated fibroblasts (CAFs), leukocytes or endothelial cells all had significantly higher expression in human CRC samples of the SSM subtype. High expression of the CAF signature was associated with poor prognosis in untreated CRC, and joint high expression of the stromal signatures predicted resistance to radiotherapy in rectal cancer. These data show that the distinctive transcriptional and clinical features of the SSM subtype can be ascribed to its particularly abundant stromal component.

Storage temperature effects on blood orange fruit quality

Orange fruits of two blood varieties (Tarocco and Moro) were stored at 8 degrees C and 22 degrees C for 85 and 106 days, respectively, and analyzed periodically for standard quality parameters (total soluble solids, total acidity, ascorbic acid, juice yield, and rind color) and sensory influencing parameters (anthocyanins, and total and free hydroxycinnamic acids). A decrease in total acidity (TA) and juice yield during storage was observed for both cultivars; total soluble solids (TSS) increased only in the Tarocco oranges stored at 8 degrees C. The increase in TSS observed for Tarocco and the simultaneous decrease in TA in both varieties resulted in a higher maturity index (TSS/TA) for the two cultivars. No loss of vitamin C was noted in Tarocco orange at either temperature, whereas a sharp reduction in vitamin C occurred in the first 50 days of storage for Moro. A significant increase in anthocyanin content was observed in Tarocco and Moro stored at 8 degrees C. Overlong storage induces extensive hydrolysis of hydroxycinnamic derivatives to free acids in Moro orange and these, in turn, could develop the malodorous vinylphenols.