Circulating tumor cells (CTCs) have potential utility as a surrogate biomarker of tumor biology via a liquid biopsy. The aim of this study was to evaluate if the nCounter NanoString assay could be used for accurate gene expression profiling of CTCs using the PAM50 research-use-only CodeSet. Analysis was performed on CTCs isolated by the ANGLE Parsortix system from healthy blood spiked with the breast cancer cell lines Hs578T, SkBr3, MDA-MB-231 or MCF7. Using cell lines as gold standard positive controls and Parsortix processed blood without spiking (unspiked) as negative controls, we found an average of 12 significantly differentially expressed genes among spiked samples versus unspiked controls. We validated our findings with the NanoStringDiff differential expression statistical method. The NanoString recommended targeted pre-amplification introduced false positive results due to pre-amplification bias, and the amplification of non-cancer genes from normal leukocytes confounded gene expression profiling of CTCs. Pre-amplification bias is a concern for other similar assays that may be used as discovery tools or target validation of transcripts of interest in gene expression profiling of CTCs. We recommend the use of an unspiked negative control when evaluating CTC technologies regarding gene expression profiling. Given that the molecular profiling of CTCs as a liquid biopsy may have clinical ramifications for potential treatment selection in future clinical trials, our study emphasizes cautious consideration of pre-analytical variables such as amplification bias in the context of liquid biopsy studies.

Whole exome sequencing (WES), targeted gene panel sequencing and single nucleotide polymorphism (SNP) arrays are increasingly used for the identification of actionable alterations that are critical to cancer care. Here, we compared The Cancer Genome Atlas (TCGA) and the Genomics Evidence Neoplasia Information Exchange (GENIE) breast cancer genomic datasets (array and next generation sequencing (NGS) data) in detecting genomic alterations in clinically relevant genes. We performed an in silico analysis to determine the concordance in the frequencies of actionable mutations and copy number alterations/aberrations (CNAs) in the two most common breast cancer histologies, invasive lobular and invasive ductal carcinoma. We found that targeted sequencing identified a larger number of mutational hotspots and clinically significant amplifications that would have been missed by WES and SNP arrays in many actionable genes such as PIK3CA, EGFR, AKT3, FGFR1, ERBB2, ERBB3 and ESR1. The striking differences between the number of mutational hotspots and CNAs generated from these platforms highlight a number of factors that should be considered in the interpretation of array and NGS-based genomic data for precision medicine. Targeted panel sequencing was preferable to WES to define the full spectrum of somatic mutations present in a tumor.

Background

We characterized the whole transcriptome of circulating tumor cells (CTCs) in stage II-III breast cancer to evaluate correlations with primary tumor biology.

Methods

CTCs were isolated from peripheral blood (PB) via immunomagnetic enrichment followed by fluorescence-activated cell sorting (IE/FACS). CTCs, PB, and fresh tumors were profiled using RNA-seq. Formalin-fixed, paraffin-embedded (FFPE) tumors were subjected to RNA-seq and NanoString PAM50 assays with risk of recurrence (ROR) scores.

Results

CTCs were detected in 29/33 (88%) patients. We selected 21 cases to attempt RNA-seq (median number of CTCs = 9). Sixteen CTC samples yielded results that passed quality-control metrics, and these samples had a median of 4,311,255 uniquely mapped reads (less than PB or tumors). Intrinsic subtype predicted by comparing estrogen receptor (ER), progesterone receptor (PR), and human epidermal growth factor receptor 2 (HER2) versus PAM50 for FFPE tumors was 85% concordant. However, CTC RNA-seq subtype assessed by the PAM50 classification genes was highly discordant, both with the subtype predicted by ER/PR/HER2 and by PAM50 tumors. Two patients died of metastatic disease, both of whom had high ROR scores and high CTC counts. We identified significant genes, canonical pathways, upstream regulators, and molecular interaction networks comparing CTCs by various clinical factors. We also identified a 75-gene signature with highest expression in CTCs and tumors taken together that was prognostic in The Cancer Genome Atlas and Molecular Taxonomy of Breast Cancer International Consortium datasets.

Conclusion

It is feasible to use RNA-seq of CTCs in non-metastatic patients to discover novel tumor biology characteristics.

Purpose

The potential utility of circulating tumor cells (CTCs) as liquid biopsies is of great interest. We hypothesized that CTC capture using EpCAM based gating is feasible for most breast cancer subtypes.

Results

Cancer cells could be recovered from all intrinsic subtypes of breast cancer with IE/FACS, however, claudin-low cell lines showed very low capture rates compared to the four other groups (p = 0.03). IE/FACS detection of CTC mimic cells was time sensitive, emphasizing controlling for pre-analytic variables in CTC studies. Median fluorescent intensity for flow cytometry and RNA flow cell type characterization were highly correlated, predicting for CTC isolation across molecular subtypes. RNA-Seq of IE/FACS sorted single cell equivalents showed high correlation compared to bulk cell lines, and distinct gene expression signatures compared to PB.

Materials and methods

Ten cell lines representing all major subtypes of breast cancer were spiked (as CTC mimics) into and recovered from peripheral blood (PB) using immunomagnetic enrichment followed by fluorescence-activated cell sorting (IE/FACS). Flow cytometry and RNA flow were used to quantify the expression of multiple breast cancer related markers of interest. Two different RNA-Seq technologies were used to analyze global gene expression of recovered sorted cells compared to bulk cell lines and PB.

Conclusions

EpCAM based IE/FACS detected and captured a portion of spiked cells from each of the 10 cell lines representing all breast cancer subtypes, including basal-like but not claudin-low cancers. The assay allows for the isolation of high quality RNA suitable for accurate RNA-Seq of heterogeneous rare cell populations.

Neuroblastoma is the most common extracranial solid tumor of childhood. While MYCN and mutant anaplastic lymphoma kinase (ALK^F1174L) cooperate in tumorigenesis, how ALK contributes to tumor formation remains unclear. Here, we used a human stem cell-based model of neuroblastoma. Mis-expression of ALK^F1174L and MYCN resulted in shorter latency compared to MYCN alone. MYCN tumors resembled adrenergic, while ALK/MYCN tumors resembled mesenchymal, neuroblastoma. Transcriptomic analysis revealed enrichment in focal adhesion signaling, particularly the extracellular matrix genes POSTN and FN1 in ALK/MYCN tumors. Patients with ALK-mutant tumors similarly demonstrated elevated levels of POSTN and FN1. Knockdown of POSTN, but not FN1, delayed adhesion and suppressed proliferation of ALK/MYCN tumors. Furthermore, loss of POSTN reduced ALK-dependent activation of WNT signaling. Reciprocally, inhibition of the WNT pathway reduced expression of POSTN and growth of ALK/MYCN tumor cells. Thus, ALK drives neuroblastoma in part through a feedforward loop between POSTN and WNT signaling.