Microfluidic System-based Time-course Tracking of Physical Proximity Between Cells and Its Effect on Gene Expression for Elucidating Live Single Cancer-immune Cell Interactions

doi:10.21203/rs.3.rs-1138231/v1

Download PDF

Article

Microfluidic System-based Time-course Tracking of Physical Proximity Between Cells and Its Effect on Gene Expression for Elucidating Live Single Cancer-immune Cell Interactions

https://doi.org/10.21203/rs.3.rs-1138231/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted

You are reading this latest preprint version

Cell-cell communication and physical interactions play a vital role in cancer initiation, homeostasis, progression, and immune response. Here, we report a system that combines live capture of different cell types, co-incubation, time-lapse imaging, and gene expression profiling of doublets using a microfluidic integrated fluidic circuit (IFC) that enables measurement of physical distances between cells and the associated transcriptional profiles due to cell-cell interactions. The temporal variations in natural killer (NK) - triple-negative breast cancer (TNBC) cell distances were tracked and compared with terminally profiled cellular transcriptomes. The results showed the time-bound activities of regulatory modules and alluded to the existence of transcriptional memory. Our experimental and bioinformatic approaches serve as a proof of concept for interrogating live cell interactions at doublet resolution, which can be applied across different cancers and cell types.

Computational Biology

Cancer Biology

Cell-cell interactions

NK cells

single cell analysis

single cell genomics

transcriptome profiling

triple-negative breast cancer

Cell-cell communication sustains the multicellular organism as an integral unit via direct physical interactions, surface receptor-ligand interaction, cell signaling from adjacent cells, nearby cells, or even distant organs^1,2. The investigation of cell-cell interaction in the tumor microenvironment (TME) is one of the barriers to understanding cancer progression and identifying new therapeutic targets³. Despite the advances in high-throughput microscopy and single-cell based molecular analysis, tools to precisely quantify live cell-cell interactions are lacking or, more recently, characterized using combinations of spatial -omics techniques on fresh-frozen or formalin-fixed complex tissues and live cell analyses^4–8. We developed a microfluidic workflow involving capture and co-incubation of live single stromal/cancer cells or doublets using the single-cell dosing mRNA-seq integrated fluidic circuit (IFC) system (Fluidigm^®), which provides both spatial and transcriptional cell-cell interactions. To demonstrate the performance for the quantification of the cell-cell interaction, we applied our platform for natural killer (NK) - triple-negative breast cancer (TNBC) cancer-immune doublets (CIDs).

TNBC was chosen as a model due to its aggressive nature and lack of estrogen receptor, progesterone receptor, and human epidermal growth factor receptor 2 that limit the use of targeted therapies^9,10. However, promising responses have been seen with immunotherapy¹¹, which is emerging as an important component of cancer treatment^12,13. Clinically promising discoveries exploit specific features of tumor-immune cell crosstalk involving immunosuppression and anti-tumoral response signaling^14–17. Among immune cells, NK cells can efficiently kill multiple neighboring cells with oncogenic transformation of surface markers^18,19. NK cell activation, combined with their capacity to enhance antibody responses, supports NK cells’ role as anticancer agents²⁰. It is speculated that genetically engineered endogenous NK cells can exert tumor immunosurveillance and influence tumor growth^21,22. However, heterogeneity is ubiquitous in human cancer making the selection of personalized treatment/therapy a challenge. It is expected that NK cell heterogeneity further contributes to NK-tumor crosstalk dynamics with differential modulation of their cytotoxic response, triggering tumor death when the balance between activation and inhibitory protein levels are considered^23–25. It is essential, therefore, to characterize the molecular level cues emanating from single NK cells when they encounter tumor cells. Thus, a single-cell platform for NK-cancer cell interaction measurement is much needed to study NK-cancer immunotherapy.

To better understand NK-tumor cell interactions, we present a microfluidic workflow involving capture and co-incubation of single NK and cancer cells (CIDs) using the Polaris^TM Single-Cell Dosing mRNA Seq IFC²⁶ (Fluidigm). The doublets captured in the microfluidic chambers were tracked for cell-cell distances, using time-lapse imaging. After 13 hours of incubation with the exchange of growth medium at a defined interval of time (5 hours), the cells were subjected to single-cell RNA sequencing (scRNA-seq). This offered a total of 290 transcriptomes, including single NK and cancer cells as control and NK-cancer cell doublets (CIDs). Unsupervised clustering analysis of the scRNA-seq expression profiles revealed heterogeneity in the TNBC cell line. We also observed a small number of NK killing events among the co-incubated CIDs, which allowed us to characterize the gene expression signature associated with NK-mediated lysis event of the companion cancer cells. We correlated the hourly computed cell-cell distances with terminally profiled gene expression vectors. We noted the existence of transcriptional memory, governed by precise regulatory modules active in a time-bound manner. In addition, we investigated the ligand-protein pairs interactions, which provided cues for the inflated activity of CD24/SIGLEC10 and ANXA1/EGFR in the tumor-NK doublets and supporting a previously described interaction between CD24/SIGLEC10 as a potent immunotherapy target for ovarian and triple negative breast cancer²⁷.

To study the interactions between NK and TNBC cells, we captured single NK-92MI cells (interleukin-2 independent Natural Killer cell line), single MDA-MB-231 cells, and cancer-immune doublets (CIDs, one NK-92MI cell and one MDA-MB-231) and incubated the cells for 13 hours using the Fluidigm Polaris system^{26,28− 31} (Fig. 1A, Fig S1). Time-lapse images of CIDs were captured every hour to measure the distance between the CIDs. Following incubation, single cells and doublets were processed for RNA-sequencing. The integrated fluidic circuit can perform on-chip multi step chemistry that includes cell capture, co-incubation, lysis, reverse transcription, and cDNA amplification²⁹. Single-cell and CIDs were analyzed for expression of marker genes that distinguished the two cell types. We identified previously known markers for both NK-92MI cells and MDA-MB-231 cells. In the case of NK cells, we observed strong expression of NK cell marker genes KLRD1, LAIR1, CCR6 and TNFRSF9. Single cancer cells showed TNBC marker genes HMGA1, ANKRD11, and TACSTD2 (Fig. 1B) when analyzed using SCANPY³² toolkit for analyzing single-cell gene expression data. Further to validate our cell type annotations, we took advantage of the imaging capability of the Polaris system. Prior to cell selection on the microfluidic IFC, the NK and cancer cells were stained with CellTracker™ Deep Red Dye and CellTracker™ Orange CMRA Dye respectively. The z-score normalized intensities of NK and cancer cell channels were subjected to 2D visualization (plotted using R package scatterD3) revealing the grouping of cells as per their annotations based on cell staining (Fig S2). Unsupervised clustering with Seurat v3 and associated Uniform Manifold Approximation and Projection (UMAP) based 2D visualization³³ of the transcriptomes revealed two separate clusters, which were primarily dominated by clonal heterogeneity of the cancer cell line (Fig. 1C). The single NK cells shared one of the clusters, cluster 1, with a cancer cell sub-group. Upon performing principal component analysis (PCA) exclusively on the transcriptomes from this cluster, we noted spatial segregation of the NK cells (Fig. 1C). To further delineate the cancer cell heterogeneity, we performed unsupervised clustering of single cancer cells separately, which also resulted in two distinct clusters (Fig S3A), featuring differentially expressed genes (Fig S3B).

Distance Tracking of Cancer-Immune Cell Interactions Over Time Shows Transcriptional Memory

We tracked CIDs for dynamic changes in the distance between cancer and the companion NK cells under incubation in the same microfluidic chamber over 13 hours (Fig. 1A). At the end of the 13th hour, transcriptome sequencing was performed on the CIDs (n = 102). This allowed us to determine the association between terminally measured gene expression with cancer/immune cell distances measured across different time points. Transcriptomes profiled at the end of the 13th hour featured transcripts that correlate significantly with CID distance measurements across all the time points (Fig. 2A, 2B). Surprisingly, we identified time-bound activities of at least three distinct gene modules (Fig. 2B). At time point 6 (after five-hour incubation and a change in culture medium), we noted a new set of genes (module 2; M2) expressed. However, we did not observe a similar shift in gene expression at later culture medium change time points. We performed module wise transcription factor activity analysis using ShinyGO³⁴ and RcisTarget³⁵, which inferred the regulatory role of three potential transcription factors, namely BRCA1, YY1, and THAP1. Among these, BRCA1 is predicted by ShinyGO, whereas YY1 and THAP1 along with their candidate target genes by RcisTarget (Fig. 2C).

Cancer-Immune Cell Killing Events

The microfluidic system allowed us to design experiments that tracked NK cell killing events and the associated transcriptomic signatures. We observed ten NK cell killing events across a total of 132 CID interactions. Differential expression analysis between CIDs subgroups showed unique gene expression signatures between NK killing and non-killing events. Among the 187 up-regulated genes in the NK killing group were CASP8, SH3BP2, IGF-1, CNPY1 and LMO1 (Fig. 1D). We applied PROGgene V2, a tool for prognostic implications measurement of genes³⁶, to the Cancer Genome Atlas (Breast Invasive Carcinoma; TCGA-BRCA patients), and analyzed the impact of these 187 up-regulated genes on overall survival. 164 out of 187 up-regulated genes overlapped with TCGA-BRCA dataset from PROGgene V2. We observed a subtle survival advantage in the patients having a higher mean expression of this combined gene signature (Hazard ratio (HR)= 0.11, p<0.05) (Fig. 3)

TNBC is a highly aggressive form of breast cancer with limited targeted treatment options. However, immunotherapy has recently turned out to be a promising strategy in the clinical management of the disease. In addition to current widely used T cell-mediated immunotherapy, research is also focused on NK cells that perform key roles in innate immune response in cancer. As such, comprehensive characterization of interacting NK and cancer cells at single-cell resolution might unravel actionable biomarkers and pathways involved in tumor growth and progression.

Here, we performed gene expression analysis of NK-cancer cell doublets resulting from the interacting distance between doublets measured over time. Statistical and bioinformatic investigation of the data enabled identification of gene expression signatures associated with successful and unsuccessful cancer cell lysis (killing events). Using TCGA gene expression profile data, we further confirmed that the identified signatures are linked to patient survival. Using single cell analysis, we also noted that the heterogeneity of the two cancer subclones confounded the process of dimension reduction by overshadowing the distinction of the NK cells’ identity.

Cell interaction distances and transcriptional memory

We noted activation of exclusive, exposure time-bound, gene regulatory modules governing cell-cell physical proximity during the life cycle of NK-cancer cell interaction under co-incubation. We also noted an association between CID cell distance and their terminally profiled transcriptome in live cells. Similarly, Gide and colleagues reported a potential association between cancer/immune cell proximity with anti-PD-1 therapy response in melanoma patients³⁷. This underscores the importance of cell-cell distance as an informative parameter to understand immunosurveillance and response in cancer. We observed changes in gene expression modules over time that correlated to the distance between live cells, suggesting the existence of transcriptional memory. Recent work has highlighted the controlled synthesis and degradation of mRNA transcripts as a major regulatory strategy influencing cell fate decisions³⁸. Transcriptional memory is a phenomenon that allows cells to retain reversible memory to respond to similar stimuli encountered in the future³⁹. At time point 6 (after five-hour incubation and a change in culture medium), we noted a new set of genes expressed. During the initial hours of co-incubation, the cells might undergo the formation of lytic synapses as reported previously^1,40, post which their interactions might take a drift.

In our study, we observed a regulatory role of three potential transcription factors, namely BRCA1, YY1, and THAP1, in transcriptional memory. BRCA1 is a well-studied tumor suppressor gene and has known implications in breast and ovarian cancers. Among the remaining two, YY1 promotes oncogenic activities in breast cancer⁴¹, whereas THAP1 plays a key role in DNA repair and is also found to be overexpressed in breast cancers⁴². One of the genes from module 1 is TRAP1, whose overexpression is involved in promoting breast tumor growth. On the other hand, it also suppresses metastasis by regulating mitochondrial dynamics⁴³. Another gene MELK from this module promotes TNBC proliferation. Targeting MELK can result in cell cycle arrest by reducing cyclin B1 and increasing p27 and p-JNK⁴⁴. EYA2 is also involved in promoting breast cancer proliferation. Its overexpression results in an increase of proliferative markers cyclin E, PCNA, and EGFR⁴⁵.

Cell-cell signaling is a major component of cancer-immune cell interactions. We used gene expressions as a surrogate for ligand-protein activities. We focused on cancer-specific ligand-protein interactions featured in the iTALK ligand-protein pairs database⁴⁶. To estimate the extent of ligand-protein coordination, we computed Pearson’s correlation coefficient for molecule pairs across CIDs and the other single cells (Fig. 4). We observed an elevated correlation between ANXA1 and EGFR in CIDs. ANXA1 is involved in mediating endocytosis of the EGFR receptor ANXA1/S100A11 complex^47–49, and mediates cell-cell communication via exosomal EGFR⁵⁰. We also observed a higher correlation between CD24/SIGLEC10 in the CIDs. Similarly, higher coordination between EGFR and HSP90AA1 was observed in CIDs.

We observed preferentially elevated co-expression between CD24 (receptor) - SIGLEC10 (ligand) transcripts among NK-cancer doublets as compared to single NK and cancer cells. The coordination among this ligand-protein pair in CIDs suggests the importance of the CD24/SIGLEC10 axis in NK/TNBC cells. A past study reported the association of SIGLEC10 in impeding NK cell function and poor patient survival in hepatocellular carcinoma (HCC); CD24/SIGLEC10 axis may thus be involved in regulating NK cell function⁵¹. In another study, the possibility of CD24/SIGLEC10 interaction in TNBC under the exposure of tumor-associated macrophages (TAMs) was reported. It has been observed that CD24 and SIGLEC10 are overexpressed in different tumor types and the TAMs, respectively²⁷. Targeting this interaction may be therapeutically important. Another ligand-protein pair observed was the coordination between EGFR and HSP90AA1. HSP90AA1 is critical for maintaining the stability and function of its client protein EGFR. This stabilization promotes pathogenesis in breast, head and neck cancer^52,53. This occurs via EMT and tumor migration activating signaling pathways in MDA-MB-231 cells⁵⁴.

Cancer-Immune Cell Killing Events

When the CIDs were analyzed for cell killing events, we noted a gene expression signature that included CASP8, SH3BP2, IGF-1, and LMO1. CASP8, when activated through FASLG, results in activation of the extrinsic pathway of apoptosis in the target cells⁵⁵. On the other hand, SH3BP2⁵⁶ and IGF-1⁵⁷ have shown to play a decisive role in NK cell development and cytotoxicity. A strong differential expression signal was observed for several genes that are largely undocumented for their NK cytotoxicity links. These include the overexpression of LMO1 (overexpressed in T lymphocytes from lymphoblastic leukemia)⁵⁸, CNPY1 and ACSBG2 in the NK-killing group, which has not been described in NK cells till today.

Currently, it is difficult to co-incubate single cells in a controlled environment that simultaneously allows study of physical interactions between live cells and its effect on gene expression⁵⁹. Here we report a cell-cell interaction study using an automated environment that precisely controls temperature, humidity, gas composition and media exchange to continuously monitor and measure the distance between cells. By processing the doublets on-IFC (in the same chambers) for cell lysis, reverse transcription, and cDNA amplification using a microfluidic multi-step chemistry, proximity measurements can be directly linked to downstream transcriptomics changes.

The use of this microfluidic system enabled the identification of unique features of NK cells’ anti-cancer activities. Highly coordinated gene expression profiles were identified as a result of dynamic changes in the physical distance of interacting NK and cancer cells, reinforcing transcriptional memory as a key regulatory strategy of cells. We could also trace increased coordination among some specific ligand-protein pairs, manifested through gene expression readouts. In the future, this microfluidic workflow could provide new leads when studying immuno-oncology cellular interactions. which may be considered while developing and administering NK cell-based immunotherapies.

Data Availability

All raw and processed sequencing data generated in this study have been submitted to the NCBI Gene Expression Omnibus (GEO; https://www.ncbi.nlm.nih.gov/geo/) under accession number GSE181591 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE181591)

Author contributions

Study Conception & Design: B.C.T.F., C.D.S., M.L., N.R., S.S.J.

Microfluidic System Development and Support: C.D.S., K.H., N.R.

Performed Experiment or Data Collection: B.C.T.F., C.D.S., N.M, P.K.K.

Computation & Statistical Analysis: S.C., N.M., N.R., D.S.

Data Interpretation & Biological Analysis: B.C.T.F., S.C., N.M, N.R., D.S., S.S.J.

Manuscript Writing, Review & Editing: B.C.T.F., S.C., N.M, L.T.D.C., N.R., D.S., S.S.J. with input from all the authors.

Supervision: L.T.D.C., N.R., D.S., S.S.J.

Acknowledgement

N.M. and P.K.K. were supported in part by the John and Marva Warnock Research Fund.

S.S.J. is supported in part by the Stanford Catalyst for Collaborative Solutions.

B.C.T.F acknowledges the support of the Brazilian National Council for Scientific and Technological Development - CNPq; the Coordination for the Improvement of Higher Education Personnel (CAPES) (grant 88881.190291/2018-01) and the Coordination of A.C.Camargo - Antônio Prudente Foundation.

D.S. acknowledges the support of an intramural start up grant from the Indraprastha Institute of Information Technology, Delhi, and is partially supported by the INSPIRE faculty grant [DST/INSPIRE/04/2015/003068] from the Department of Science and Technology (DST), Govt. of India.

Declaration of interests

K.H and N.R. are employees and stockholders of Fluidigm Corporation.

M.L is a former employee and stockholder of Fluidigm Corporation.

D.S. is a partner and equity holder at CareOnco Biotech Pvt. Ltd.

S.S.J. serves as a scientific advisor for Quantumcyte and Ravel Biotechnology.

Standardization of NK cell activation

The activation of NK cells was established, in principle, by the evaluation of cell lines NK-92MI and MDA-MB-231. Activation was assessed after 24 and 48 hours of cell incubation in a proportion of 3 NK cells for each MDA-MB-231 by evaluating the markers (CD25, CD69, and CD314) expressed from cell-cell contact, using flow cytometry (Sony SH800S).

Culturing MDA-MB-231 and NK-92MI

MDA-MB-231 cells (triple-negative breast carcinoma, ATCC® HTB-26 ™) were cultured using the DMEM (high modified) culture medium supplemented with 10% inactivated fetal bovine serum and 1% penicillin/streptomycin, replicating the culture every three days, keeping the culture at a confluence of 40% at the time of passage, at 37°C in a humid atmosphere at 5% CO2. The adherent cells were dissociated with the TrypLE reagent (Gibco) and resuspended in the complete culture medium.

NK-92MI cells (human NK cells genetically modified to produce interleukin 2, ATCC® CRL-2408) were cultured using the AMEM (Alpha Minimum Essential medium) culture medium, plus 0.2 mM inositol, 0.1 mM mercaptoethanol, 0.02 mM folic acid, 12.5% fetal bovine serum and 12.5% horse serum. The cells were homogenized to separate the clusters of NKs prior to the replication of the culture into a new 25 cm² flask (1mL of cells cultured previously + 9mL new culture medium). Both cell lines were cultured separately for the subsequent cell-cell interaction experiment using the Polaris system.

Selection and incubation using the Polaris system (Fluidigm)

The first step is to prime the IFC with beads that will allow the adhesion of the cells to be incubated in the microchambers. After treatment, reagents and cells already labeled using specific markers to differentiate cell types (for NK cells, celltracker far red, and celltracker orange for cancer cells) and viability (calcein AM) were pipetted into the IFC. The selection of the cells to be incubated was performed on the Polaris system. It was configured for NKs and cancer cells to select positive cells for celltracker far red and calcein AM.

Some IFC wells were maintained with only one cell for later comparison of the gene expression of single cells with the gene expression obtained when the cell was incubated in doublets (NK cells + breast cancer cell lineage). After selection, the equipment was programmed to keep the cells in incubation for 16 hours, performing the replacement of culture medium (20% DMEM + 80% AMEM) every 5 hours, and time-lapse images every hour. Images were captured immediately before and after culture medium change.

71 single cancer cells, 77 single NK cells, and 132 cells in doublets (NK + cancer cells) were incubated, followed by the cell lysis, reverse transcription, and cDNA amplification. Subsequently, sequencing of the single-cell RNA, using NextSeq (Illumina) was performed.

Distance measurement between cells

The distance shown in our study is the shortest distance between the membrane of MDA-MB-231 cell and NK cell on IFC within the doublets (NK + cancer cells) group. The distance has been measured for 13 time points. To analyze the videos frame by frame to yield the distance data, we used ImageJ, a public domain Java-based image processing software developed at the National Institutes of Health⁶⁰.

Data preprocessing

In total, we obtained expression profiles of 340 cells. The cells can be classified as follows. 1. Single NK cells (n = 77); 2. single tumor cells (n = 71); CIDs throughout all time points (n =132), and CIDs that were left with NK cells alone at the terminal time point (n=10). Out of the 340 cells, we removed 50 expression profiles. The cells discarded include — 1. 2 bulk RNA-Seq replicates of each of the NK-92MI and 2 MDA-MB-231 cell-lines; 4 empty chambers (no cells could be traced in the chambers from the start to the end); 8 empty chambers that initially contained tumor cells; 10 empty chambers that initially contained NK cells; 2 CIDs that started with single NK cells; 22 tumor cells that started as CIDs. After removal of these transcriptomes, we were left with 290 cells/doublets. Further, we discarded cells having less than 2000 expressed genes (non-zero RNA-seq by Expectation Maximization (RSEM) expression value⁶¹). Next, we retained genes having RSEM expression >5 in at least 10 cells. After these filtering steps, we were left with an expression matrix constituting 290 cells and protein-coding 8907 genes.

Batch correction, clustering, and visualization of the single cells

The matrix obtained after performing the basic pre-processing steps was used as input for the Seurat, single cell analysis R package, as well as for other downstream analyses. We used Seurat’s data integration workflow to process the data with the genes detected in at least 5 cells, further employing the standard routines for log-normalization, variance stabilizing transformation for identification of highly variable genes using NormalizeData() and FindVariableFeatures() with default parameter settings. The cells used in this study originated from two independent runs. For integration and batch correction, we used the FinIntegrationAnchors() with k.filter=100 for identification of anchor cells that represent matching cell pairs across the two datasets in order to project the transcriptomes into a shared space and IntegrateData() function was used for integrating these anchors, which involve Canonical Correlation Analysis (CCA). These steps provided the batch corrected matrix. 2-D map of the cells was created using the RunUMAP() function.

Differential gene expression analysis

Limma-voom⁶² was used to identify differentially expressed genes (DEGs) among the cell-groups. Top DEGs that qualified adjusted P-value cutoff of 0.05 and absolute log2 fold change cutoff of 1 were further analyzed for their biological significance.

Survival analysis based on the upregulated genes governing NK cell antitumor activity

PROGgeneV2³⁶ was used for gene signature based (based on a gene list) overall survival analysis combined gene signature analysis functionality using the TCGA-BRCA dataset, containing survival information for 594 patients. Out of 187 upregulated genes in NK cells showing antitumor activity, 164 genes were found overlapping with the TCGA-BRCA dataset. PROGgeneV2 produced a Kaplan-Meir plot, segregating the patients into high and low risk groups using the median value of the combined gene expression signature as cut-off.

Regulation of cell-cell distance and transcriptional memory

We considered two matrices to track the association between temporally recorded cell-cell distance and terminally profiled gene expression. First, the processed gene expression matrix of dimensions |G| ×|C|, where |G|, and |C| denote the number of genes and the number of CIDs, respectively. Second, the cell-cell distance matrix of dimensions |C| ×|T|, where |T|denotes the time points when the cell-cell distance measurements were recorded. We used these two matrices to compute the third matrix of dimensions |G| ×|T| containing the Pearson’s correlation coefficients 〖ρ〗_ g,tfor each gene-time point pairs {(g,t) | t∈T,g∈G}. Notably, |G|=2000,|C|=102,and |T|=13. We retained 90 genes with 〖|ρ〗_ g,t|≥0.25for at least one time point for further analysis. These 90 genes were subjected to hierarchical clustering based on the computed correlation values. 4 gene modules were retrieved using cutree(). For each of these modules, motif enrichment and transcription factor activity analyses were performed using RcisTarget³⁵ and ShinyGO³⁴, respectively. RcisTarget discerned enriched TF binding motifs and reported a module-specific list of TFs, using the hg19-tss-centered-10kb-7species.mc9nr.feather database containing genome-wide ranking for the motifs. Gene regulatory networks were constructed using igraph⁶³.

Cell-cell signaling among CIDs

We used iTALK R package featuring information on 2,648 cancer-specific ligand-receptor interactions. We applied the rawParse() function to select the top 50% highly expressed genes from our processed gene expression matrix (constituting 290 cells and 8907 genes), using mean as a method of choice. This identified 230 ligand-receptor pairs using the FindLR() function. We computed Pearson’s correlation coefficient between expression values associated with the selected ligand-receptor pairs across CIDs (TU-NK/TU-NK, TU-NK/NK) and qualified ones that exhibited Pearson's correlation coefficient > 0.4. To rule out the possibility that the co-expression stems solely from either NK or Tumor cells alone, we tracked Pearson's correlation coefficient among NK/NK and TU/TU cases as well. With these criteria, 20 ligand-receptor pairs were selected, and three of these were found directly involved in breast cancer signaling.

Kumar S (2018). Natural killer cell cytotoxicity and its regulation by inhibitory receptors. Immunology 154, 383-93. 10.1111/imm.12921.
Wortzel I, Dror S, Kenific CM, Lyden D (2019). Exosome-Mediated Metastasis: Communication from a Distance. Dev Cell 49, 347-60. 10.1016/j.devcel.2019.04.011.
Dongre A, Weinberg RA (2019). New insights into the mechanisms of epithelial-mesenchymal transition and implications for cancer. Nat Rev Mol Cell Biol 20, 69-84. 10.1038/s41580-018-0080-4.
Baccin C, et al. (2020). Combined single-cell and spatial transcriptomics reveal the molecular, cellular and spatial bone marrow niche organization. Nat Cell Biol 22, 38-48. 10.1038/s41556-019-0439-6.
Berglund E, et al. (2018). Spatial maps of prostate cancer transcriptomes reveal an unexplored landscape of heterogeneity. Nat Commun 9, 2419. 10.1038/s41467-018-04724-5.
Bode D, Cull AH, Rubio-Lara JA, Kent DG (2021). Exploiting Single-Cell Tools in Gene and Cell Therapy. Front Immunol 12, 702636. 10.3389/fimmu.2021.702636.
Giladi A, et al. (2020). Dissecting cellular crosstalk by sequencing physically interacting cells. Nat Biotechnol 38, 629-37. 10.1038/s41587-020-0442-2.
Lewis SM, Asselin-Labat ML, Nguyen Q, Berthelet J, Tan X, Wimmer VC, Merino D, Rogers KL, Naik SH (2021). Spatial omics and multiplexed imaging to explore cancer biology. Nat Methods. 10.1038/s41592-021-01203-6.
Abdel-Latif M, Youness RA (2020). Why natural killer cells in triple negative breast cancer? World J Clin Oncol 11, 464-76. 10.5306/wjco.v11.i7.464.
Lyons TG (2019). Targeted Therapies for Triple-Negative Breast Cancer. Curr Treat Options Oncol 20, 82. 10.1007/s11864-019-0682-x.
Vagia E, Mahalingam D, Cristofanilli M (2020). The Landscape of Targeted Therapies in TNBC. Cancers 12, 916.
Emens LA (2018). Breast Cancer Immunotherapy: Facts and Hopes. Clin Cancer Res 24, 511-20. 10.1158/1078-0432.CCR-16-3001.
Pan C, Liu H, Robins E, Song W, Liu D, Li Z, Zheng L (2020). Next-generation immuno-oncology agents: current momentum shifts in cancer immunotherapy. Journal of Hematology & Oncology 13, 29. 10.1186/s13045-020-00862-w.
Allard B, Longhi MS, Robson SC, Stagg J (2017). The ectonucleotidases CD39 and CD73: Novel checkpoint inhibitor targets. Immunol Rev 276, 121-44. 10.1111/imr.12528.
Castro F, Cardoso AP, Goncalves RM, Serre K, Oliveira MJ (2018). Interferon-Gamma at the Crossroads of Tumor Immune Surveillance or Evasion. Front Immunol 9, 847. 10.3389/fimmu.2018.00847.
Jang JE, Hajdu CH, Liot C, Miller G, Dustin ML, Bar-Sagi D (2017). Crosstalk between Regulatory T Cells and Tumor-Associated Dendritic Cells Negates Anti-tumor Immunity in Pancreatic Cancer. Cell Rep 20, 558-71. 10.1016/j.celrep.2017.06.062.
Schmidt L, Eskiocak B, Kohn R, Dang C, Joshi NS, DuPage M, Lee DY, Jacks T (2019). Enhanced adaptive immune responses in lung adenocarcinoma through natural killer cell stimulation. Proc Natl Acad Sci U S A 116, 17460-69. 10.1073/pnas.1904253116.
Pahl J, Cerwenka A (2017). Tricking the balance: NK cells in anti-cancer immunity. Immunobiology 222, 11-20. 10.1016/j.imbio.2015.07.012.
Sanchez-Correa B, et al. (2019). DNAM-1 and the TIGIT/PVRIG/TACTILE Axis: Novel Immune Checkpoints for Natural Killer Cell-Based Cancer Immunotherapy. Cancers (Basel) 11. 10.3390/cancers11060877.
Shimasaki N, Jain A, Campana D (2020). NK cells for cancer immunotherapy. Nat Rev Drug Discov 19, 200-18. 10.1038/s41573-019-0052-1.
Ingegnere T, et al. (2019). Human CAR NK Cells: A New Non-viral Method Allowing High Efficient Transfection and Strong Tumor Cell Killing. Front Immunol 10, 957. 10.3389/fimmu.2019.00957.
Sivori S, Meazza R, Quintarelli C, Carlomagno S, Della Chiesa M, Falco M, Moretta L, Locatelli F, Pende D (2019). NK Cell-Based Immunotherapy for Hematological Malignancies. J Clin Med 8. 10.3390/jcm8101702.
Sanseviero E (2019). NK Cell-Fc Receptors Advance Tumor Immunotherapy. J Clin Med 8. 10.3390/jcm8101667.
Sungur CM, Murphy WJ (2014). Positive and negative regulation by NK cells in cancer. Crit Rev Oncog 19, 57-66. 10.1615/critrevoncog.2014010805.
Yang C, et al. (2019). Heterogeneity of human bone marrow and blood natural killer cells defined by single-cell transcriptome. Nat Commun 10, 3931. 10.1038/s41467-019-11947-7.
Sanada CD, Ooi AT (2019). Single-Cell Dosing and mRNA Sequencing of Suspension and Adherent Cells Using the Polaris(TM) System. Methods Mol Biol 1979, 185-95. 10.1007/978-1-4939-9240-9_12.
Barkal AA, et al. (2019). CD24 signalling through macrophage Siglec-10 is a target for cancer immunotherapy. Nature 572, 392-96. 10.1038/s41586-019-1456-0.
Mayer S, et al. (2019). Multimodal Single-Cell Analysis Reveals Physiological Maturation in the Developing Human Neocortex. Neuron 102, 143-58 e7. 10.1016/j.neuron.2019.01.027.
Ramalingam N, et al. (2017). Fluidic Logic Used in a Systems Approach to Enable Integrated Single-Cell Functional Analysis. Frontiers in Bioengineering and Biotechnology 4. 10.3389/fbioe.2016.00070.
Wills QF, et al. (2017). The nature and nurture of cell heterogeneity: accounting for macrophage gene-environment interactions with single-cell RNA-Seq. BMC Genomics 18, 53. 10.1186/s12864-016-3445-0.
Wu S, Zhang H, Fouladdel S, Li H, Keller E, Wicha MS, Omenn GS, Azizi E, Guan Y (2020). Cellular, transcriptomic and isoform heterogeneity of breast cancer cell line revealed by full-length single-cell RNA sequencing. Comput Struct Biotechnol J 18, 676-85. 10.1016/j.csbj.2020.03.005.
Wolf FA, Angerer P, Theis FJ (2018). SCANPY: large-scale single-cell gene expression data analysis. Genome Biol 19, 15. 10.1186/s13059-017-1382-0.
Stuart T, et al. (2019). Comprehensive Integration of Single-Cell Data. Cell 177, 1888-902 e21. 10.1016/j.cell.2019.05.031.
Ge SX, Jung D, Yao R (2020). ShinyGO: a graphical gene-set enrichment tool for animals and plants. Bioinformatics 36, 2628-29. 10.1093/bioinformatics/btz931.
Aibar S, et al. (2017). SCENIC: single-cell regulatory network inference and clustering. Nat Methods 14, 1083-86. 10.1038/nmeth.4463.
Goswami CP, Nakshatri H (2014). PROGgeneV2: enhancements on the existing database. BMC Cancer 14, 970. 10.1186/1471-2407-14-970.
Gide TN, et al. (2020). Close proximity of immune and tumor cells underlies response to anti-PD-1 based therapies in metastatic melanoma patients. Oncoimmunology 9, 1659093. 10.1080/2162402X.2019.1659093.
Battich N, Beumer J, de Barbanson B, Krenning L, Baron CS, Tanenbaum ME, Clevers H, van Oudenaarden A (2020). Sequencing metabolically labeled transcripts in single cells reveals mRNA turnover strategies. Science 367, 1151-56. 10.1126/science.aax3072.
Beiter T, Niess AM, Moser D (2020). Transcriptional memory in skeletal muscle. Don't forget (to) exercise. J Cell Physiol 235, 5476-89. 10.1002/jcp.29535.
Orange JS (2008). Formation and function of the lytic NK-cell immunological synapse. Nat Rev Immunol 8, 713-25. 10.1038/nri2381.
Wan M, et al. (2012). Yin Yang 1 plays an essential role in breast cancer and negatively regulates p27. Am J Pathol 180, 2120-33. 10.1016/j.ajpath.2012.01.037.
Shinoda K, et al. (2021). The dystonia gene THAP1 controls DNA double-strand break repair choice. Mol Cell. 10.1016/j.molcel.2021.03.034.
Zhang B, et al. (2015). Aberrantly upregulated TRAP1 is required for tumorigenesis of breast cancer. Oncotarget 6, 44495-508. 10.18632/oncotarget.6252.
Li G, Yang M, Zuo L, Wang MX (2018). MELK as a potential target to control cell proliferation in triple-negative breast cancer MDA-MB-231 cells. Oncol Lett 15, 9934-40. 10.3892/ol.2018.8543.
Xu H, Jiao Y, Yi M, Zhao W, Wu K (2019). EYA2 Correlates With Clinico-Pathological Features of Breast Cancer, Promotes Tumor Proliferation, and Predicts Poor Survival. Front Oncol 9, 26. 10.3389/fonc.2019.00026.
Wang Y, et al. iTALK: an R Package to Characterize and Illustrate Intercellular Communication: bioRxiv, 2019.
Boudhraa Z, Bouchon B, Viallard C, D'Incan M, Degoul F (2016). Annexin A1 localization and its relevance to cancer. Clin Sci (Lond) 130, 205-20. 10.1042/CS20150415.
de Graauw M, et al. (2014). Annexin A2 depletion delays EGFR endocytic trafficking via cofilin activation and enhances EGFR signaling and metastasis formation. Oncogene 33, 2610-9. 10.1038/onc.2013.219.
Jaiswal JK, Lauritzen SP, Scheffer L, Sakaguchi M, Bunkenborg J, Simon SM, Kallunki T, Jaattela M, Nylandsted J (2014). S100A11 is required for efficient plasma membrane repair and survival of invasive cancer cells. Nat Commun 5, 3795. 10.1038/ncomms4795.
Raulf N, Lucarelli P, Thavaraj S, Brown S, Vicencio JM, Sauter T, Tavassoli M (2018). Annexin A1 regulates EGFR activity and alters EGFR-containing tumour-derived exosomes in head and neck cancers. Eur J Cancer 102, 52-68. 10.1016/j.ejca.2018.07.123.
Zhang P, Lu X, Tao K, Shi L, Li W, Wang G, Wu K (2015). Siglec-10 is associated with survival and natural killer cell dysfunction in hepatocellular carcinoma. J Surg Res 194, 107-13. 10.1016/j.jss.2014.09.035.
Ahsan A, et al. (2012). Wild-type EGFR is stabilized by direct interaction with HSP90 in cancer cells and tumors. Neoplasia 14, 670-7. 10.1593/neo.12986.
Friedland JC, Smith DL, Sang J, Acquaviva J, He S, Zhang C, Proia DA (2014). Targeted inhibition of Hsp90 by ganetespib is effective across a broad spectrum of breast cancer subtypes. Invest New Drugs 32, 14-24. 10.1007/s10637-013-9971-6.
Tian Y, Wang C, Chen S, Liu J, Fu Y, Luo Y (2019). Extracellular Hsp90alpha and clusterin synergistically promote breast cancer epithelial-to-mesenchymal transition and metastasis via LRP1. J Cell Sci 132. 10.1242/jcs.228213.
Zhu Y, Huang B, Shi J (2016). Fas ligand and lytic granule differentially control cytotoxic dynamics of natural killer cell against cancer target. Oncotarget 7, 47163-72. 10.18632/oncotarget.9980.
Jevremovic D, Billadeau DD, Schoon RA, Dick CJ, Leibson PJ (2001). Regulation of NK cell-mediated cytotoxicity by the adaptor protein 3BP2. J Immunol 166, 7219-28. 10.4049/jimmunol.166.12.7219.
Ni F, Sun R, Fu B, Wang F, Guo C, Tian Z, Wei H (2013). IGF-1 promotes the development and cytotoxic activity of human NK cells. Nat Commun 4, 1479. 10.1038/ncomms2484.
Liu H, Liu M, You H, Li X, Li X (2020). Oncogenic Network and Hub Genes for Natural Killer/T-Cell Lymphoma Utilizing WGCNA. Front Oncol 10, 223. 10.3389/fonc.2020.00223.
Armingol E, Officer A, Harismendy O, Lewis NE (2021). Deciphering cell-cell interactions and communication from gene expression. Nat Rev Genet 22, 71-88. 10.1038/s41576-020-00292-x.
Schneider CA, Rasband WS, Eliceiri KW (2012). NIH Image to ImageJ: 25 years of image analysis. Nat Methods 9, 671-5. 10.1038/nmeth.2089.
Li B, Dewey CN (2011). RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics 12, 323. 10.1186/1471-2105-12-323.
Law CW, Chen Y, Shi W, Smyth GK (2014). voom: Precision weights unlock linear model analysis tools for RNA-seq read counts. Genome Biol 15, R29. 10.1186/gb-2014-15-2-r29.
Csardi G, Nepusz T (2006). The igraph software package for complex network research. InterJournal, complex systems 1695, 1-9.

Yes there is potential Competing Interest. K.H and N.R. are employees and stockholders of Fluidigm Corporation. M.L is a former employee and stockholder of Fluidigm Corporation. D.S. is a partner and equity holder at CareOnco Biotech Pvt. Ltd. S.S.J. serves as a scientific advisor for Quantumcyte and Ravel Biotechnology.

SupplementaryFigures.docx

Download PDF

Version 1

posted

You are reading this latest preprint version

Microfluidic System-based Time-course Tracking of Physical Proximity Between Cells and Its Effect on Gene Expression for Elucidating Live Single Cancer-immune Cell Interactions

Status:

Version 1

Abstract

Figures

Introduction

Results

Distance Tracking of Cancer-Immune Cell Interactions Over Time Shows Transcriptional Memory

Cancer-Immune Cell Killing Events

Discussion

Cell interaction distances and transcriptional memory

Cancer-Immune Cell Killing Events

Declarations

Material & Methods

References

Additional Declarations

Supplementary Files

Status:

Version 1