http://orcid.org/0000-0001-8375-7313
Correction
4 Dec 2018: Shih AJ, Menzin A, Whyte J, Lovecchio J, Liew A, et al. (2018) Correction: Identification of grade and origin specific cell populations in serous epithelial ovarian cancer by single cell RNA-seq. PLOS ONE 13(12): e0208778. https://doi.org/10.1371/journal.pone.0208778 View correction
Figures
Abstract
Here we investigated different cell populations within ovarian cancer using single-cell RNA seq: fourteen samples from nine patients with differing grades (high grade, low grade and benign) as well as different origin sites (primary and metastatic tumor site, ovarian in origin and fallopian in origin). We were able to identify sixteen distinct cell populations with specific cells correlated to high grade tumors, low grade tumors, benign and one population unique to a patient with a breast cancer relapse. Furthermore the proportion of these populations changes from primary to metastatic in a shift from mainly epithelial cells to leukocytes with few cancer epithelial cells in the metastases. Differential gene expression shows myeloid lineage cells are the primary cell group expressing soluble factors in primary samples while fibroblasts do so in metastatic samples. The leukocytes that were captured did not seem to be suppressed through known pro-tumor cytokines from any of the cell populations. Single cell RNA-seq is necessary to de-tangle cellular heterogeneity for better understanding of ovarian cancer progression.
Citation: Shih AJ, Menzin A, Whyte J, Lovecchio J, Liew A, Khalili H, et al. (2018) Identification of grade and origin specific cell populations in serous epithelial ovarian cancer by single cell RNA-seq. PLoS ONE 13(11): e0206785. https://doi.org/10.1371/journal.pone.0206785
Editor: Sandra Orsulic, Cedars-Sinai Medical Center, UNITED STATES
Received: August 21, 2018; Accepted: October 18, 2018; Published: November 1, 2018
Copyright: © 2018 Shih et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: The data discussed in this publication have been deposited in NCBI’s Gene Expression Omnibus and are accessible through GEO Series accession number GSE118828 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE118828).
Funding: This study was funded by the following private funders to The Feinstein Institute: Linda and Dr. Myron Teitelbaum Foundation, Monter Cancer Center at Northwell Health Cancer Institute, Moms Who Kick, Manhasset Women’s Coalition Against Breast Cancer, Henry Gabbay, Stand Up 2 Cancer, and the John L. Lovecchio, MD Gynecologic Oncology Research Fund. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Ovarian cancer is the 5th leading cause of cancer deaths in women living in the United States [1]. In 2017 alone, over 22,000 women were diagnosed with ovarian cancer and approximately 14,000 died from their disease. Epithelial ovarian cancers (EOCs) are broken down into four histological subgroups: serous, mucinous, endometroid and clear cell [2]. Serous ovarian cancers are the most common, comprising ~50% of cases [3] and can be further subdivided into high grade serous ovarian cancer (HGSOC) and low grade serous ovarian cancer (LGSOC) at ~90% and ~10% respectively [2].
The treatment options for HGSOC and LGSOC are similar with two standards of care, 1) Primary cytoreductive surgery (PCS) followed by adjuvant chemotherapy and 2) Neoadjuvant chemotherapy (NACT) followed by interval cytoreductive surgery (ICS) and adjuvant chemotherapy [4]. For both types of treatments, a combination chemotherapy consisting of platinum-based and taxane-based drugs (e.g. carboplatin and paclitaxel) is administered [5].
Each EOC type has a different prognosis with HGSOC having the worst outcome with the highest mortality rate of all the gynecological cancers. HGSOC is highly curable (>90% 5 year survival) if diagnosed early as local disease when the cancer is confined to the ovaries; however, most women are diagnosed with advanced stage metastatic disease. For these women, cure rates are exceptionally low: less than 25% of patients with late stage HGSOC will live more than 5 years [6]. Although HGSOC is initially sensitive to chemotherapy, almost invariably, relapse occurs followed by drug resistant progressive disease [7,8]. Overall 5-year survival of LGSOC is higher at 75% although this is dependent on having no residual disease following surgery as LGSOC has greater chemoresistance [9].
There has been no significant decrease in mortality rates in almost 30 years [10]. There is a critical need to improve our understanding of the underlying mechanisms leading to drug resistant ovarian cancer and identify potentially actionable therapeutic targets.
Ovarian cancer is a complex disease with significant tumor heterogeneity and as such there has been little success in identifying actionable targets. A meta-analysis of gene expression data from 1251 HGSOC tumors did not identify a collective prognostic gene expression signature [11]. Recently, Patch et al [12] performed a comprehensive whole genome analysis (DNA, RNA, miRNA, CNV, methylation) on 92 HGSOC tumors associated with different levels of drug response (i.e. refractory, resistant, or sensitive). Even with this in-depth analysis, no actionable targets or prognostic molecular profiles were identified. Given the rarity of LGSOC and its chemoresistance, most studies have been focused on recurrent disease using chemotherapy or hormone replacement therapy with mixed results [13]. Overall, analyses of bulk tumor tissue at the DNA and RNA levels have failed to provide results of significant clinical value.
Single-cell RNA-seq (scRNA-seq) allows for the quantitative and qualitative analysis of cell composition in complex tissues without a priori knowledge of the cell populations present. Several thousand genes can be quantitated simultaneously at the individual cell level. Using this approach we sought to identify commonalities and differences in cell composition of tumor samples from women with differing grades of serous epithelial ovarian cancer. Previous studies have examined ovarian cancer at the single cell level [14,15]; here we expand by assaying many more cells as well as identifying cell type specific differential expression.
Materials and methods
Subject recruitment, sample acquisition and sample processing
Women scheduled for surgery to evaluate a suspicious pelvic mass were recruited through the Tissue Donation Program at The Feinstein Institute. Pathological discard tissue, primary tumor and metastatic lesions (when available), were obtained at time of surgery and frozen for later analysis. Tissue was minced and frozen in 40% FBS, 40% RPMI and 20% DMSO. This freezing protocol routinely yields greater than 85% viable cells. This study was approved by the Institutional Review Board (IRB) of Northwell Health. See Table 1 for study subject characteristics.
Single-cell suspensions
Thawed tissue samples were digested and dissociated into single-cell suspensions using Miltenyi’s gentleMACS Octo dissociator and tumor dissociation kit following manufacturer’s instructions. Single cells were collected by straining digested tissue through a MACS Smart strainer (70um), washed, then layered over a Ficoll gradient to remove red blood cells and debris. To enrich the tumor cell fraction for live cells, dead cells were excluded using Dead Cell Removal kit and MS columns (Miltenyi Biotec) and remaining viable cells were prepared for scRNA-seq.
Single-cell RNA-seq (scRNA-seq)
Using the BioRad droplet digital SEQ Single-cell Isolator and the Illumina SureCell Whole Transcriptome Analysis 3’ library prep kit, scRNA-seq was performed on isolated cells from ovarian tumor tissue samples following the manufacturer’s instructions. Briefly single-cells were encapsulated, lysed then barcoded within each droplet. Following first and second strand cDNA synthesis, Illumina’s Nextera technology was used to generate a library for NGS. Final libraries were assessed and quantified using a High Sensitivity DNA chip and a 2100 BioAnalyzer (Agilent Technologies) prior to sequencing on a NexSeq 500 high output flow cell.
Data analysis
Sample data was demultiplexed and raw sequence files generated using bcl2fastq2 app in Illumina’s cloud application, Basespace, followed by the SureCell RNA Single-Cell application to extract the cell barcodes and the Unique Molecular Identifiers (UMIs). The raw sequence files were aligned to the hg19 human reference genome and gene expression was quantified. A knee plot of the genic UMIs per cell barcode was used to set thresholds for cells for downstream analysis. The data discussed in this publication have been deposited in NCBI’s Gene Expression Omnibus [16] and are accessible through GEO Series accession number GSE118828 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE118828).
The R [17] software package Seurat was used for further analysis [18,19]. Genes were initially filtered on expression in at least three cells and each cell needed to have at least 200 genes expressed. The entire dataset was then log-normalized with several factors regressed out: total number of UMIs and cell cycle scores. A principal component analysis of the most variable genes was performed and an elbow plot was used to select the principal components (PCs) capturing the most variance in the dataset. These PCs were used as edge weights in an unsupervised graph-based clustering to identify cell clusters. T-distributed stochastic neighbor embedding (tSNE) was used for visualization of the cell clusters. Expression levels of cell-type specific markers were used to determine the putative identities of each cell cluster. All colors used in plots were from Paul Tol’s color schemes and templates [20].
Once the cell clusters were identified, differential expression was done using a likelihood-ratio test tuned for single-cell expression, which assumes an underlying negative-binomial distribution as well as accounts for drop out using a zero-inflation regression. We defined differentially expressed genes as those with an adjusted p-value less than 0.05, an average log2 fold change greater than 1. For potential biomarkers we put further thresholds of being expressed in at least 50% of cells of the defining cluster with at least 50% difference of detectable expression of cells to the comparison group.
Results
Patient description
A total of 2911 cells were captured by scRNA-seq on nine patients with fourteen samples, see Table 1. All nine patients had a sample from the primary site, four had samples from a metastatic site and one had a sample from a normal ovary; all of the metastatic samples were from omentum proximal to the primary site. Five of the patients were diagnosed with HGSOC (with one being fallopian in origin), two were diagnosed with LGSOC, one had a cancer that was peritoneal in origin and one had a benign cystadenoma. Average age of diagnosis for the patients was 62.4 ± 6.5 years. All the tumors were diagnosed as stage III or greater with the exception of the benign cystadenoma and one of the LGSOC, which was at stage IA. Two of the patients (HG3 and NA1) had a previous history of breast cancer with one of them receiving neoadjuvant chemotherapy. The self-reported ethnicity had a breakdown of four White, one Asian, one Black/African-American, one with multiple races and two that did not self-report. Hematoxylin and eosin staining showed no visible lymphocyte infiltration in either tumor or metastatic lesions.
Cell Identification
Pooled scRNA-seq data from all samples (tumor, normal, benign) identified sixteen distinct major clusters of cells that correlated with tumor origin, tumor status or known cellular gene expression patterns (Fig 1). The major clusters broadly divided into epithelial/mesothelial cells (EPCAM/CD326, KRT), lymphocytes (PTPRC/CD45, CD3E, CD19, MS4A1/CD20), endothelial cells (PECAM-1/CD31, CD34), fibroblasts (ACTA2, DCN, ACTB) and stromal cells (THY1/CD90, ENG/CD105, VIM, CD44), among others (S1 Fig). Pooled primary tumor data alone resulted in a large epithelial cell cluster that could be further subdivided into four groups, one group unique to HG3, one shared by HGSOC primary tumor samples, one shared by both LGSOC primary tumor samples and one cluster enriched in the ovarian tumor of fallopian origin HG2F. As expected normal ovary tissue displayed a significant presence of fibroblasts and stromal cells with no evidence of epithelial cells. ScRNA-seq data from a benign ovarian malignancy exhibited a unique cell cluster with epithelial markers, but is mesothelial in origin with other cell types divided among other clusters.
CD45 positive lymphocytes separated into five clusters with CD4+ and CD8+ T-cells (CD3D, CD4 and CD8, respectively), B-cells (CD19 and MS4A1/CD20 positive) and two clusters of myeloid lineage (CD14 positive). The myeloid clusters appeared to originate either from primary tumors or the metastatic lesion. Only the PECAM1/CD31 positive cluster of endothelial progenitor cells overlapped with all tissue samples analyzed.
Cell population differences
Each subject sample had differing relative proportions of identified cells types, see Fig 2 for both percentages and numbers of cell types identified. An unsupervised hierarchical clustering of populations of cells from all subjects resulted in five groups. There were three major groups, which were predominantly from HGSOC primary tumor samples, HGSOC metastatic tumor samples and LGSOC samples. There were single outlier groups from benign, a primary tumor sample from HG3 and the normal and peritoneal samples clustering together. Surprisingly, the neoadjuvant treated primary tumor sample clustered with metastatic lesions. The outlier groups of benign and normal ovary samples were predictably comprised mostly of benign mesothelial cells as well as normal fibroblasts and stromal cells, respectively.
Top is a dendrogram grouping similar samples by cell populations. Middle is the total number of cells captured. Bottom are bar plots of each cell type. Number following colon in legend is number of cells captured for each identified cell type.
A definitive shift in cell populations between HGSOC primary tumors and corresponding HGSOC metastatic samples was observed (Fig 2). Primary tumor samples were predominantly epithelial cells (68.3% epithelial cells versus 11.1% lymphocytes), whereas corresponding metastatic samples were predominantly immune cells (66.2% lymphocytes versus 10.5% epithelial cells) with lesser amounts of epithelial cells, fibroblasts and stromal cells. Shifting from mainly epithelial to mainly lymphocytes seems to drive two of the big clusters (HGSOC primary and HGSOC metastatic) as well as the placement of the neoadjuvant treated sample in HGSOC metastatic: treatment killed off tumor epithelial cells. Clustering of LGSOC samples are based on the presence of LGSOC epithelial cells, which interestingly were present in a large quantity in the metastatic LGSOC sample. Two metastatic samples, HG2F-M and HG3-M, had significantly higher amounts of B-cells compared to all other samples (p-value < 2.2E-16 for both via Fisher’s test).
Tumor epithelial cell expression
Using Seurat, potential biomarker discovery for all cancer epithelial cells versus every other cell type was performed (Table 2). Twelve differentially expressed genes were identified which included WFDC2, two members of the claudin family (CLDN3 and CLDN4) as well as three keratin genes (KRT8, KRT18, KRT19).
Examining each tumor epithelial cluster versus the other three showed two clusters with potential biomarkers. HGSOC-F epithelial cells had unique identifiers with the top three expressed genes being TPPP3, C20orf85 and FOXJ1 (Table 3). The epithelial cell cluster unique to patient HG3 from this tumor had eleven differentially expressed genes with highest expression found in SST, TFF3 and PIGR. There were no distinct gene profiles for either of the HGSOC or LGSOC ovarian tumor epithelial cells, all highly expressed genes were also expressed within other epithelial clusters.
Non-epithelial cell expression
Comparison of normal fibroblasts and stromal cells to their tumor counterparts showed a striking upregulation of collagen genes including COL1A1, COL1A2, COL3A1, COL4A1, COL4A2 COL5A1, COL5A2, COL6A1 and COL6A3 as well as collagen support genes, SPARC, SERPINH1 and SERPINE1 (Table 4). COL4A1, COL4A2, SPARC, SERPINH1, and SERPINE1 were also expressed within the endothelial progenitor cluster. Cancer fibroblast and stromal cells had increased expression of MMP and MMP-related genes: MMP2, MMP11 and TIMP1. While many secreted factors were expressed by myeloid lineage cells (primary and metastatic), metastatic fibroblasts and cancer stromal cells, as expected relatively few were observed in tumor epithelial cells (S1 Table).
Doing a differential expression between the two clusters of cancer fibroblasts, we found metastatic fibroblasts were found to express higher levels of soluble factors compared to primary fibroblasts including, CXCL12, S100A6, S100A10, SFRP2,SFRP4,IGF1, CXCL14, ANGPTL4 and IL6 (Table 5). Metastatic fibroblasts also had increased expression of two complement cascade genes C3, CFB and an inhibitor of C1 gene, SERPING1. Tumor fibroblasts had increased expression of C7. Comparison between the primary tumor myeloid lineage cells and the metastatic myeloid lineage cells showed elevated levels of CC2, CC3, CC4, CXCL8 and TNF in the primary tumor.
Comparison of cell populations to TCGA HGSOC sub-types
The Cancer Genome Atlas (TCGA) collected 489 HGSOC samples and performed mRNA expression, microRNA expression, promoter methylation and DNA copy number assays [21]. Clustering of the samples based on both mRNA and microRNA expression showed four stable subtypes within HGSOC that were labeled differentiated, immunoreactive, mesenchymal and proliferative.
Using genes sets defining the four HGSOC subtypes, we used the average expression of groups of cells to classify samples, first by patient sample ID then by identified cell clusters. The threshold of classification was defined as two standard deviations above the mean of the combined group comparisons. By patient sample ID, only two patients were identified with HGSOC sub-types (Table 6), LG2-P as differentiated and LG2-M as mesenchymal, which are LGSOC samples.
However using identified cell clusters, there was significantly more correlation with the prescribed function of each sub-type (Table 6). Benign epithelial cells and LGSOC epithelial cells were identified as the differentiated subtype. Primary and metastatic myeloid lineage cells showed increased expression of the immunoresponsive subtype. Primary fibroblasts, metastatic fibroblasts and cancer stromal cells were identified as the mesenchymal subtype while normal fibroblasts and normal stromal cells had no similarities to any of the HGSOC subtypes. None of the clusters showed enrichment of proliferative subtype. Here we see that myeloid lineage defines immunoresponsive subtype, while tumor fibroblasts and stromal cells define mesenchymal, and epithelial cells define differentiated. The proliferative subtype may either be a cancer type or a rare cell population that was not captured in this study.
Discussion
With scRNA-seq, we can investigate intra and inter tumor heterogeneity gene expression differences at cellular resolution. Given that epithelial cells are the predominant dysregulated cell type in epithelial ovarian cancer, we first wanted to compare these cells to all others and their potential relationship to cancer progression. The genes we identified at the single cell level in all cancer epithelial cells are concordant with previous bulk RNA seq studies. WFDC2 has been shown to be a biomarker for ovarian cancer [22] and overexpression promotes ovarian tumor growth [23]. Similarly, SCGB2A1 is expressed in all ovarian cancers [24]. CD24 has prognostic value in survival as a marker for cancer stem cells [25–28] whereas CLDN3 and CLDN4 have been shown to regulate the epithelial to mesenchymal transition [29]. These expression profiles imply that in ovarian cancer, epithelial cells are transitioning into a mesenchymal state, which is a hallmark of cancer.
Although all the malignant tumors had been classified as serous ovarian cancer, there are specific subsets of cancerous epithelial cells that correlate with pathological findings. Specific expression patterns of subsets of tumor epithelial cells expressed functional genes related to their origin. HGSOC-F epithelial cells biomarker TPPP3 is a tubulin polymerization protein that has been previously implicated in colorectal cancer [30], while the other three genes are necessary for proper cilia function. FOXJ1 is a differentiating factor for ciliated epithelial cells in the neonatal oviduct [31], RSPH1 is necessary for proper cilia function [32] and ZMYND10 affects dynenin motor function in cilia [33]. Given that these cells are fallopian in origin, their expression of cilia related genes is not unexpected as cilia are necessary in normal fallopian function. However these genes may also provide extra motility and transitioning to a mesenchymal state for the subtype of HGSOC that is fallopian in origin.
Two genes identified in HG3 tumor epithelial cells are broadly implicated in many cancers with worse prognoses, PIGR [34,35] and AGR2 [36,37]. They are the only cells in this dataset that express two soluble factors, a cytokine CXCL17 which has been found to promote tumor progression in other cancers [38–40] and a hormone SST, which suppresses other growth factors and has been considered as a potential pharmaceutical target in cancers [41,42]. HG3 is the only patient that had a previous cancer diagnosis without neoadjuvant chemotherapy: a breast malignancy more than ten years prior to her ovarian cancer diagnosis. Whether this cancer is a late recurrence of breast cancer or different in context of increased cancer risk, the biomarkers here imply a more aggressive phenotype with different activating factors.
The idea of pre-metastatic niche was postulated by Stephen Paget in 1889 as a tumor “seed” finding its appropriate “soil” [43]. Formation of a pre-metastatic niche involves tumor associated cells to create a hospitable microenvironment for tumor cells through specific secreted factors [44]. Secreted factors were mainly being produced by cells of myeloid lineage and metastatic fibroblasts especially compared to epithelial cells in primary tumor samples. Given there were few epithelial cells in metastatic samples it is unclear if they are producing many secreted factors. Therefore myeloid lineage cells are crucial in both primary and metastatic ovarian tumors, with fibroblasts playing a more important role in metastatic growth than in primary.
Collagen has a complex role in cancer progression. Initially thought to provide a barrier to tumor invasion, collagen has been shown to remodel the ECM to help promote angiogenesis, tumor invasion and migration [45]. Furthermore, matrix metalloproteases work to degrade the ECM to increase tumor invasion, particularly in metastasis [46]. Increased expression of both types of genes in ECM remodeling is thought to play a role in fibrosis/desmoplasia that is seen in many types of cancer [47]. We saw significant desmoplasia in pathology slides that was increased in metastatic samples compared to primary samples. High expression of several collagens and MMP-related genes in our dataset was primarily in cancer fibroblasts and stromal cells, not directly from the cancer epithelial cells themselves. This supports the theory that tumor stroma plays an important role in the tumor development [48]. In addition, there was high expression of COL4A1, COL4A2, SERPINH1 and SPARC in endothelial cells. Type IV-collagen has previously been shown to be necessary in angiogenic functions of endothelial cells [49,50], with SERPINE1 implicated in maintaining ECM of arterial walls [51]. SPARC and SERPINH1 play important roles in type-IV collagen synthesis and function [52]. Concurrent expression of these genes in four different clusters of cells imply collagen IV is crucial in the recruitment of endothelial cells for angiogenesis. Furthermore, as endothelial cells are the only cluster of cells that are represented in normal, primary tumor and metastatic tumor this implies there is no cell or class state switch in endothelial cells and these cells are functioning normally given their microenvironment and stimuli.
Creating a pro-tumor environment also requires avoiding immune surveillance from B-cells and T-cells, which can have both negative regulation and positive regulation of cancer, depending on specific tumor microenvironment (B-cells reviewed in [53], T-cells reviewed in [54]). Regulatory B-cells are expected to promote tumor growth by suppressing immune response through IL-10 and STAT3. Our cluster of B-cells do not seem to be regulatory B-cells, as they are CD38- and CD27- and with most similar marker expression as mature naive. More importantly, they do not express any IL-10 or STAT3. We also only see minor sporadic expression of markers for T-reg cells (PDCD1/PD-1, FOXP3), with little expression of pro-tumor cytokines, IL-10 and CCL22. Presence of immune cells in metastatic tumors is paradoxical, especially in light of no expression of pro-tumor factors. In fact we actually see expression of interferon-Ɣ in CD8 cells which is thought to by anti-tumor [55,56]. Recent research has shown a pro-tumor role of interferon-Ɣ through expression of PD-L1 [57], but expression of PD-L1 is not evident in our dataset.
One critical component of managing the immune system is the complement system, which regulates the innate immune system and plays a role in tumor progression [58,59]. Metastatic fibroblast cells overexpress C3, which is one of the key components that must be activated in the complement system. In murine models C3 has roles in tumor growth [60] and angiogenesis [61]. Furthermore, ascitic fluid taken from ovarian cancer patients showed surface deposition of C3 with increased CFH (which is only mildly overexpressed in these cells in our dataset) [62,63]. The other genes overexpressed in these cells are inhibitors, CFB and SERPING1. Adding to this, metastatic fibroblasts are also the only cells expressing CXCL12, which along with CXCR4 play a pivotal role across many cancers [64–66], including ovarian cancer [67,68]. We see CXCR4 expressed only in B-cells and T-cells. Metastatic fibroblasts also express two S100 family proteins, S100A10 and S100A6, which are important in creating inflammation for tumor growth and metastases [69,70]. Taken together it appears that metastatic fibroblasts are creating an inflammatory environment through CXCL12 and S100 family proteins while suppressing immune response through complement system inactivation.
Conclusions
Our results show four major malignant epithelial cell types seen in nine patients diagnosed with either high grade or low grade serous epithelial ovarian cancer. Gene profiles were found for all four in aggregate or for two clusters individually, particularly HGSOC that is fallopian in origin and one in a patient with breast cancer recurrence. These may serve either as surveillance option, e.g. CA125 in ovarian cancer [71] or treatment, e.g. HER2 inhibitors in breast cancer [72,73]. Secondly, we found that primary tumor epithelial cells secrete much less factors in comparison to cells of myeloid lineage or fibroblasts. Between primary and metastatic tumors, fibroblasts secrete increasing levels of necessary factors to fuel metastatic growth, with the caveat that we captured only a few metastatic epithelial cells and may see more secreted factors with a larger sample size. Metastatic fibroblasts could potentially be tumor epithelial cells that have fully transformed (as many of the tumor epithelial expressed genes imply), but metastatic fibroblasts we observed do not express PROM1/CD133, a well described cancer stem cell marker [74]. Third, once ovarian cancer metastasizes, relatively few tumor cells are required to create and maintain a favorable microenvironment as evidenced by the shift of cell populations from epithelial cells to lymphocytes. Fourth, B-cells and T-cells do not appear to be suppressed through known pro-tumor cytokines, but could be suppressed indirectly through complement pathways. Finally, using previous data from ovarian cancers analyzed by TCGA, we find cancer subtype to be correlated with specific cell types.
Only through high-resolution studies at single cell resolution are we able to identify and quantify heterogeneity within tumors. Our ability to study patient derived primary and corresponding metastatic lesions using high-throughput single cell analysis represents a unique opportunity to study ovarian cancer without a priori knowledge of tumor and stromal cell inter-relationships. Thus single cell assessment of patient samples can provide critical information needed to understand ovarian cancer progression.
Supporting information
S1 Fig. Coloring of tSNE plot by known markers defining each cluster of cells.
https://doi.org/10.1371/journal.pone.0206785.s001
(TIFF)
S1 Table. Gene expression of secreted factors increased in each cluster.
https://doi.org/10.1371/journal.pone.0206785.s002
(XLS)
S2 Table. All genes upregulated in each cluster.
https://doi.org/10.1371/journal.pone.0206785.s003
(XLS)
Acknowledgments
We are exceptionally grateful to patients and their families.
We also gratefully acknowledge the research coordinators and staff members of the Northwell Health Tissue Donation Program for facilitating the collection of tissue samples.
References
- 1. Siegel RL, Miller KD, Jemal A. Cancer Statistics, 2017. CA Cancer J Clin. 2017;67: 7–30. pmid:28055103
- 2. Seidman JD, Horkayne-Szakaly I, Haiba M, Boice CR, Kurman RJ, Ronnett BM. The histologic type and stage distribution of ovarian carcinomas of surface epithelial origin. Int J Gynecol Pathol. 2004;23: 41–44. pmid:14668549
- 3. Stuart GCE, Kitchener H, Bacon M, duBois A, Friedlander M, Ledermann J, et al. 2010 Gynecologic Cancer InterGroup (GCIG) consensus statement on clinical trials in ovarian cancer: report from the Fourth Ovarian Cancer Consensus Conference. Int J Gynecol Cancer. 2011;21: 750–755. pmid:21543936
- 4. Wright AA, Bohlke K, Armstrong DK, Bookman MA, Cliby WA, Coleman RL, et al. Neoadjuvant chemotherapy for newly diagnosed, advanced ovarian cancer: Society of Gynecologic Oncology and American Society of Clinical Oncology Clinical Practice Guideline. Gynecol Oncol. 2016;143: 3–15. pmid:27650684
- 5. Morgan RJ Jr, Armstrong DK, Alvarez RD, Bakkum-Gamez JN, Behbakht K, Chen L-M, et al. Ovarian Cancer, Version 1.2016, NCCN Clinical Practice Guidelines in Oncology. J Natl Compr Canc Netw. 2016;14: 1134–1163. pmid:27587625
- 6.
Howlader N, Noone AM, Krapcho M, Miller D, Bishop K, Kosary CL, Yu M, Ruhl J, Tatalovich Z, Mariotto A, Lewis DR, Chen HS, Feuer EJ. SEER Cancer Statistics Review, 1975–2014. Cronin KA, editor. Bethesda, MD: National Cancer Institute.; based on November 2016 SEER data submission, posted to the SEER web site, April 2017.
- 7. Cooke SL, Brenton JD. Evolution of platinum resistance in high-grade serous ovarian cancer. Lancet Oncol. 2011;12: 1169–1174. pmid:21742554
- 8. McGuire WP, Hoskins WJ, Brady MF, Kucera PR, Partridge EE, Look KY, et al. Cyclophosphamide and cisplatin compared with paclitaxel and cisplatin in patients with stage III and stage IV ovarian cancer. N Engl J Med. 1996;334: 1–6. pmid:7494563
- 9. Kaldawy A, Segev Y, Lavie O, Auslender R, Sopik V, Narod SA. Low-grade serous ovarian cancer: A review. Gynecol Oncol. 2016;143: 433–438. pmid:27581327
- 10. Parmar MKB, Ledermann JA, Colombo N, du Bois A, Delaloye J-F, Kristensen GB, et al. Paclitaxel plus platinum-based chemotherapy versus conventional platinum-based chemotherapy in women with relapsed ovarian cancer: the ICON4/AGO-OVAR-2.2 trial. Lancet. 2003;361: 2099–2106. pmid:12826431
- 11. Waldron L, Haibe-Kains B, Culhane AC, Riester M, Ding J, Wang XV, et al. Comparative meta-analysis of prognostic gene signatures for late-stage ovarian cancer. J Natl Cancer Inst. 2014;106. pmid:24700801
- 12. Patch A-M, Christie EL, Etemadmoghadam D, Garsed DW, George J, Fereday S, et al. Whole-genome characterization of chemoresistant ovarian cancer. Nature. 2015;521: 489–494. pmid:26017449
- 13. Della Pepa C, Tonini G, Santini D, Losito S, Pisano C, Di Napoli M, et al. Low Grade Serous Ovarian Carcinoma: From the molecular characterization to the best therapeutic strategy. Cancer Treat Rev. 2015;41: 136–143. pmid:25573350
- 14. Winterhoff BJ, Maile M, Mitra AK, Sebe A, Bazzaro M, Geller MA, et al. Single cell sequencing reveals heterogeneity within ovarian cancer epithelium and cancer associated stromal cells. Gynecol Oncol. 2017;144: 598–606. pmid:28111004
- 15. Gonzalez VD, Samusik N, Chen TJ, Savig ES, Aghaeepour N, Quigley DA, et al. Commonly Occurring Cell Subsets in High-Grade Serous Ovarian Tumors Identified by Single-Cell Mass Cytometry. Cell Rep. 2018;22: 1875–1888. pmid:29444438
- 16. Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30: 207–210. pmid:11752295
- 17.
R Core Team. R: A language and environment for statistical computing [Internet]. Vienna, Austria; 2017. https://www.R-project.org/
- 18. Satija R, Farrell JA, Gennert D, Schier AF, Regev A. Spatial reconstruction of single-cell gene expression data. Nat Biotechnol. 2015;33: 495–502. pmid:25867923
- 19. Macosko EZ, Basu A, Satija R, Nemesh J, Shekhar K, Goldman M, et al. Highly Parallel Genome-wide Expression Profiling of Individual Cells Using Nanoliter Droplets. Cell. 2015;161: 1202–1214. pmid:26000488
- 20.
Tol P. Paul Tol’s Notes, Colour schemes and templates [Internet]. May 2018. https://personal.sron.nl/~pault/
- 21. Cancer Genome Atlas Research Network. Integrated genomic analyses of ovarian carcinoma. Nature. 2011;474: 609–615. pmid:21720365
- 22. Hellström I, Raycraft J, Hayden-Ledbetter M, Ledbetter JA, Schummer M, McIntosh M, et al. The HE4 (WFDC2) protein is a biomarker for ovarian carcinoma. Cancer Res. 2003;63: 3695–3700. pmid:12839961
- 23. Moore RG, Hill EK, Horan T, Yano N, Kim K, MacLaughlan S, et al. HE4 (WFDC2) gene overexpression promotes ovarian tumor growth. Sci Rep. 2014;4: 3574. pmid:24389815
- 24. Bellone S, Tassi R, Betti M, English D, Cocco E, Gasparrini S, et al. Mammaglobin B (SCGB2A1) is a novel tumour antigen highly differentially expressed in all major histological types of ovarian cancer: implications for ovarian cancer immunotherapy. Br J Cancer. 2013;109: 462–471. pmid:23807163
- 25. Lim J, Lee K-M, Shim J, Shin I. CD24 regulates stemness and the epithelial to mesenchymal transition through modulation of Notch1 mRNA stability by p38MAPK. Arch Biochem Biophys. 2014;558: 120–126. pmid:24977325
- 26. Wang G, Zhang Z, Ren Y. TROP-1/Ep-CAM and CD24 are potential candidates for ovarian cancer therapy. Int J Clin Exp Pathol. 2015;8: 4705–4714. pmid:26191160
- 27. Kristiansen G, Denkert C, Schlüns K, Dahl E, Pilarsky C, Hauptmann S. CD24 is expressed in ovarian cancer and is a new independent prognostic marker of patient survival. Am J Pathol. 2002;161: 1215–1221. pmid:12368195
- 28. Gao M-Q, Choi Y-P, Kang S, Youn JH, Cho N-H. CD24+ cells from hierarchically organized ovarian cancer are enriched in cancer stem cells. Oncogene. 2010;29: 2672–2680. pmid:20190812
- 29. Lin X, Shang X, Manorek G, Howell SB. Regulation of the Epithelial-Mesenchymal Transition by Claudin-3 and Claudin-4. PLoS One. 2013;8: e67496. pmid:23805314
- 30. Ye K, Li Y, Zhao W, Wu N, Liu N, Li R, et al. Knockdown of Tubulin Polymerization Promoting Protein Family Member 3 inhibits cell proliferation and invasion in human colorectal cancer. J Cancer. 2017;8: 1750–1758. pmid:28819371
- 31. Okada A, Ohta Y, Brody SL, Watanabe H, Krust A, Chambon P, et al. Role of foxj1 and estrogen receptor alpha in ciliated epithelial cell differentiation of the neonatal oviduct. J Mol Endocrinol. 2004;32: 615–625. pmid:15171704
- 32. Onoufriadis A, Shoemark A, Schmidts M, Patel M, Jimenez G, Liu H, et al. Targeted NGS gene panel identifies mutations in RSPH1 causing primary ciliary dyskinesia and a common mechanism for ciliary central pair agenesis due to radial spoke defects. Hum Mol Genet. 2014;23: 3362–3374. pmid:24518672
- 33. Kobayashi D, Asano-Hoshino A, Nakakura T, Nishimaki T, Ansai S, Kinoshita M, et al. Loss of zinc finger MYND-type containing 10 (zmynd10) affects cilia integrity and axonemal localization of dynein arms, resulting in ciliary dysmotility, polycystic kidney and scoliosis in medaka (Oryzias latipes). Dev Biol. 2017;430: 69–79. pmid:28823919
- 34. Qi X, Li X, Sun X. Reduced expression of polymeric immunoglobulin receptor (pIgR) in nasopharyngeal carcinoma and its correlation with prognosis. Tumour Biol. 2016;37: 11099–11104. pmid:26910773
- 35. Yue X, Ai J, Xu Y, Chen Y, Huang M, Yang X, et al. Polymeric immunoglobulin receptor promotes tumor growth in hepatocellular carcinoma. Hepatology. 2017;65: 1948–1962. pmid:28073159
- 36. Pohler E, Craig AL, Cotton J, Lawrie L, Dillon JF, Ross P, et al. The Barrett’s antigen anterior gradient-2 silences the p53 transcriptional response to DNA damage. Mol Cell Proteomics. 2004;3: 534–547. pmid:14967811
- 37. Tian S-B, Tao K-X, Hu J, Liu Z-B, Ding X-L, Chu Y-N, et al. The prognostic value of AGR2 expression in solid tumours: a systematic review and meta-analysis. Sci Rep. 2017;7: 15500. pmid:29138453
- 38. Guo YJ, Zhou YJ, Yang XL, Shao ZM, Ou ZL. The role and clinical significance of the CXCL17-CXCR8 (GPR35) axis in breast cancer. Biochem Biophys Res Commun. 2017;493: 1159–1167. pmid:28943434
- 39. Ohlsson L, Hammarström M-L, Lindmark G, Hammarström S, Sitohy B. Ectopic expression of the chemokine CXCL17 in colon cancer cells. Br J Cancer. 2016;114: 697–703. pmid:26889977
- 40. Li L, Yan J, Xu J, Liu C-Q, Zhen Z-J, Chen H-W, et al. CXCL17 expression predicts poor prognosis and correlates with adverse immune infiltration in hepatocellular carcinoma. PLoS One. 2014;9: e110064. pmid:25303284
- 41. Sun L-C, Coy DH. Somatostatin receptor-targeted anti-cancer therapy. Curr Drug Deliv. 2011;8: 2–10. pmid:21034425
- 42. Rai U, Thrimawithana TR, Valery C, Young SA. Therapeutic uses of somatostatin and its analogues: Current view and potential applications. Pharmacol Ther. 2015;152: 98–110. pmid:25956467
- 43. Paget S. THE DISTRIBUTION OF SECONDARY GROWTHS IN CANCER OF THE BREAST. Lancet. 1889;133: 571–573.
- 44. Peinado H, Zhang H, Matei IR, Costa-Silva B, Hoshino A, Rodrigues G, et al. Pre-metastatic niches: organ-specific homes for metastases. Nat Rev Cancer. 2017;17: 302–317. pmid:28303905
- 45. Fang M, Yuan J, Peng C, Li Y. Collagen as a double-edged sword in tumor progression. Tumour Biol. 2014;35: 2871–2882. pmid:24338768
- 46. Deryugina EI, Quigley JP. Matrix metalloproteinases and tumor metastasis. Cancer Metastasis Rev. 2006;25: 9–34. pmid:16680569
- 47. Rybinski B, Franco-Barraza J, Cukierman E. The wound healing, chronic fibrosis, and cancer progression triad. Physiol Genomics. 2014;46: 223–244. pmid:24520152
- 48. Davidson B, Trope CG, Reich R. The Role of the Tumor Stroma in Ovarian Cancer. Front Oncol. 2014;4. pmid:24860785
- 49. Bahramsoltani M, Slosarek I, De Spiegelaere W, Plendl J. Angiogenesis and collagen type IV expression in different endothelial cell culture systems. Anat Histol Embryol. 2014;43: 103–115. pmid:23551189
- 50. Dunér P, Gonçalves I, Grufman H, Edsfeldt A, To F, Nitulescu M, et al. Increased aldehyde-modification of collagen type IV in symptomatic plaques—a possible cause of endothelial dysfunction. Atherosclerosis. 2015;240: 26–32. pmid:25746374
- 51. Ruigrok YM, Rinkel GJE, van’t Slot R, Wolfs M, Tang S, Wijmenga C. Evidence in favor of the contribution of genes involved in the maintenance of the extracellular matrix of the arterial wall to the development of intracranial aneurysms. Hum Mol Genet. 2006;15: 3361–3368. pmid:17038484
- 52. Kuo DS, Labelle-Dumais C, Gould DB. COL4A1 and COL4A2 mutations and disease: insights into pathogenic mechanisms and potential therapeutic targets. Hum Mol Genet. 2012;21: R97–110. pmid:22914737
- 53. Sarvaria A, Madrigal JA, Saudemont A. B cell regulation in cancer and anti-tumor immunity. Cell Mol Immunol. 2017;14: 662–674. pmid:28626234
- 54. Santoiemma PP, Powell DJ Jr. Tumor infiltrating lymphocytes in ovarian cancer. Cancer Biol Ther. 2015;16: 807–820. pmid:25894333
- 55. Burke F, Smith PD, Crompton MR, Upton C, Balkwill FR. Cytotoxic response of ovarian cancer cell lines to IFN-gamma is associated with sustained induction of IRF-1 and p21 mRNA. Br J Cancer. 1999;80: 1236–1244. pmid:10376977
- 56. Marth C, Helmberg M, Mayer I, Fuith LC, Daxenbichler G, Dapunt O. Effects of biological response modifiers on ovarian carcinoma cell lines. Anticancer Res. 1989;9: 461–467. pmid:2502063
- 57. Abiko K, Matsumura N, Hamanishi J, Horikawa N, Murakami R, Yamaguchi K, et al. IFN-γ from lymphocytes induces PD-L1 expression and promotes progression of ovarian cancer. Br J Cancer. 2015;112: 1501–1509. pmid:25867264
- 58. Pio R, Corrales L, Lambris JD. The role of complement in tumor growth. Adv Exp Med Biol. 2014;772: 229–262. pmid:24272362
- 59. Afshar-Kharghan V. The role of the complement system in cancer. J Clin Invest. 2017;127: 780–789. pmid:28248200
- 60. Cho MS, Vasquez HG, Rupaimoole R, Pradeep S, Wu S, Zand B, et al. Autocrine effects of tumor-derived complement. Cell Rep. 2014;6: 1085–1095. pmid:24613353
- 61. Nunez-Cruz S, Gimotty PA, Guerra MW, Connolly DC, Wu Y-Q, DeAngelis RA, et al. Genetic and pharmacologic inhibition of complement impairs endothelial cell function and ablates ovarian cancer neovascularization. Neoplasia. 2012;14: 994–1004. pmid:23226093
- 62. Junnikkala S, Hakulinen J, Jarva H, Manuelian T, Bjørge L, Bützow R, et al. Secretion of soluble complement inhibitors factor H and factor H-like protein (FHL-1) by ovarian tumour cells. Br J Cancer. 2002;87: 1119–1127. pmid:12402151
- 63. Bjørge L, Hakulinen J, Vintermyr OK, Jarva H, Jensen TS, Iversen OE, et al. Ascitic complement system in ovarian cancer. Br J Cancer. 2005;92: 895–905. pmid:15726105
- 64. Shen B, Zheng M-Q, Lu J-W, Jiang Q, Wang T-H, Huang X-E. CXCL12-CXCR4 promotes proliferation and invasion of pancreatic cancer cells. Asian Pac J Cancer Prev. 2013;14: 5403–5408. pmid:24175834
- 65. Sun Y, Mao X, Fan C, Liu C, Guo A, Guan S, et al. CXCL12-CXCR4 axis promotes the natural selection of breast cancer cell metastasis. Tumour Biol. 2014;35: 7765–7773. pmid:24810923
- 66. Sun X, Cheng G, Hao M, Zheng J, Zhou X, Zhang J, et al. CXCL12 / CXCR4 / CXCR7 chemokine axis and cancer progression. Cancer Metastasis Rev. 2010;29: 709–722. pmid:20839032
- 67. Guo Q, Gao B-L, Zhang X-J, Liu G-C, Xu F, Fan Q-Y, et al. CXCL12-CXCR4 Axis Promotes Proliferation, Migration, Invasion, and Metastasis of Ovarian Cancer. Oncol Res. 2014;22: 247–258. pmid:26629936
- 68. Popple A, Durrant LG, Spendlove I, Rolland P, Scott IV, Deen S, et al. The chemokine, CXCL12, is an independent predictor of poor survival in ovarian cancer. Br J Cancer. 2012;106: 1306–1313. pmid:22415233
- 69. Bresnick AR, Weber DJ, Zimmer DB. S100 proteins in cancer. Nat Rev Cancer. 2015;15: 96–109. pmid:25614008
- 70. Lukanidin E, Sleeman JP. Building the niche: the role of the S100 proteins in metastatic growth. Semin Cancer Biol. 2012;22: 216–225. pmid:22381352
- 71. Bast RC Jr, Feeney M, Lazarus H, Nadler LM, Colvin RB, Knapp RC. Reactivity of a monoclonal antibody with human ovarian carcinoma. J Clin Invest. 1981;68: 1331–1337. pmid:7028788
- 72. Slamon DJ, Clark GM, Wong SG, Levin WJ, Ullrich A, McGuire WL. Human breast cancer: correlation of relapse and survival with amplification of the HER-2/neu oncogene. Science. 1987;235: 177–182. pmid:3798106
- 73. Slamon DJ, Leyland-Jones B, Shak S, Fuchs H, Paton V, Bajamonde A, et al. Use of chemotherapy plus a monoclonal antibody against HER2 for metastatic breast cancer that overexpresses HER2. N Engl J Med. 2001;344: 783–792. pmid:11248153
- 74. Burgos-Ojeda D, Rueda BR, Buckanovich RJ. Ovarian cancer stem cell markers: prognostic and therapeutic implications. Cancer Lett. 2012;322: 1–7. pmid:22334034