- Research article
- Open Access
Collagen density regulates the activity of tumor-infiltrating T cells
Journal for ImmunoTherapy of Cancer volume 7, Article number: 68 (2019)
Tumor progression is accompanied by dramatic remodeling of the surrounding extracellular matrix leading to the formation of a tumor-specific ECM, which is often more collagen-rich and of increased stiffness. The altered ECM of the tumor supports cancer growth and metastasis, but it is unknown if this effect involves modulation of T cell activity. To investigate if a high-density tumor-specific ECM could influence the ability of T cells to kill cancer cells, we here studied how T cells respond to 3D culture in different collagen densities.
T cells cultured in 3D conditions surrounded by a high or low collagen density were imaged using confocal fluorescent microscopy. The effects of the different collagen densities on T cell proliferation, survival, and differentiation were examined using flow cytometry. Cancer cell proliferation in similar 3D conditions was also measured. Triple-negative breast cancer specimens were analyzed for the number of infiltrating CD8+ T cells and for the collagen density. Whole-transcriptome analyses were applied to investigate in detail the effects of collagen density on T cells. Computational analyses were used to identify transcription factors involved in the collagen density-induced gene regulation. Observed changes were confirmed by qRT-PCR analysis.
T cell proliferation was significantly reduced in a high-density matrix compared to a low-density matrix and prolonged culture in a high-density matrix led to a higher ratio of CD4+ to CD8+ T cells. The proliferation of cancer cells was unaffected by the surrounding collagen-density. Consistently, we observed a reduction in the number of infiltrating CD8+ T-cells in mammary tumors with high collagen-density indicating that collagen-density has a role in regulating T cell abundance in human breast cancer.
Whole-transcriptome analysis of 3D-cultured T cells revealed that a high-density matrix induces downregulation of cytotoxic activity markers and upregulation of regulatory T cell markers. These transcriptional changes were predicted to involve autocrine TGF-β signaling and they were accompanied by an impaired ability of tumor-infiltrating T cells to kill autologous cancer cells.
Our study identifies a new immune modulatory mechanism, which could be essential for suppression of T cell activity in the tumor microenvironment.
Solid tumors consist of cancer cells interacting with the tumor microenvironment, which includes stromal cells, immune cells, and the extracellular matrix. Infiltration of tumors by lymphocytes, and in particular CD8+ cytotoxic T cells, is known to predict good prognosis in many types of cancer . T cell infiltration into solid tumors is also associated with increased clinical efficacy of immunotherapies . The beneficial effect of a high abundance of tumor-infiltrating T cells reflects the ability of the immune cells to mount a response against cancer cells . As an important way for the cancer cells to evade immune destruction, tumors can develop a strongly immunosuppressive tumor microenvironment [4,5,6]. This includes the accumulation of cell types with immunosuppressive activity, such as tumor-associated macrophages (TAMs), myeloid-derived suppressor cells (MDSCs), and regulatory T cells (Tregs) . Upregulation of PD-L1 by cells in the tumor microenvironment constitutes an important and well-studied immune escape mechanism in which the interaction with its receptor PD-1 on T cells lead to inactivation of the T cells . Antibody-mediated blockade of the PD1-PD-L1 interaction has demonstrated remarkable clinical efficacy for many cancer patients and stimulated the research aiming at identifying additional targetable immunosuppressive mechanisms in the tumor microenvironment .
In addition to non-malignant stromal cells and immune cells, the tumor microenvironment consists of the extracellular matrix (ECM). Degradation of the ECM surrounding a tumor is an essential part of invasive cancer growth and a main reason for the destruction of the normal tissue [8,9,10]. Of central importance, the degradation of the ECM is accompanied by the deposition of a different tumor-specific ECM . This new ECM is often of increased density and stiffness, and contains components that are not typically present in the original ECM [12, 13]. Within the last decade, a number of important discoveries have emphasized how the ECM can affect cancer biology [14, 15]. A strong correlation between the density of collagen type I, which is the most abundant component of the tumor ECM, and poor prognosis of breast cancer, gastric cancer and oral cancer has been demonstrated [16,17,18], and in vitro studies have shown that a high-density and stiff ECM can induce a process in epithelial cells resembling malignant transformation [19,20,21]. Other cell types such as fibroblasts and mesenchymal stem cells have also been demonstrated to respond to the mechanical properties of the surrounding ECM by a process termed cellular mechanosensing [22, 23]. Through the effects on cancer cells and stromal cells, the ECM can augment many of the hallmarks of cancer, such as the induction of angiogenesis  and the activation of invasion and metastasis [15, 25]. It remains quite speculative if the ECM can also modulate the immunosuppressive tumor microenvironment and thereby support the cancer’s evasion of immune destruction [15, 26]. It should, however, be noted that the presentation of antigens linked to a stiff surface has been demonstrated to impair TCR-mediated T cell activation, suggesting that T cells possess mechanosensing abilities  and others have confirmed that the TCR is affected by mechanical force . Furthermore, tumor-associated remodeling of the ECM, can lead to the deposition of ECM components such as osteopontin, SPARC, versican, and tenascin C, which have been suggested to possess immunosuppressive properties [29,30,31].
Although cell culture in 3D environments is widely used in the field of cancer biology , 3D culture of T cells is less common and has mainly been used for the study of cell migration [33, 34]. The interaction between collagen and tumor-infiltrating T cells was, however, studied by Salmon et al. through the use of elegant ex vivo culture of lung tumor slices combined with real-time imaging . Here, collagen fibers were suggested to prevent the migration of T cells from the stromal compartment into the tumor islets. Collagen-density mediated inhibition of directional T cell migration has also been suggested as the reason for intratumoral T cell exclusion in pancreatic cancers . It was not addressed in these studies if the collagen also influences the activity of the T cells. In this study, we employed 3D culture assays to investigate if the collagen-density can directly impact the activity of T cells.
Materials and methods
T cell isolation and culture
Human peripheral blood mononuclear cells (PBMCs) were isolated from healthy donors by gradient centrifugation using Lymphoprep (Alere AS) separation and frozen in fetal calf serum (Sigma Aldrich) with 10% DMSO (Sigma Aldrich). For RNAseq experiments, the PBMCs were enriched for T cells by allowing the cells to adhere overnight and collecting only the non-adherent and loosely adherent cells. For confocal microscopy, T cells were isolated from healthy donors at the Combined Technical Research Core facilities at the National Institute of Dental and Craniofacial Research (NIDCR), NIH using elutriation. For all other experiments, T cells were isolated from PBMCs using magnetic anti-CD3 microbeads (Miltenyi Biotec) according to the manufacturer’s instruction. Cells were cultured at 37 °C in a humidified 5% CO2 environment in X-vivo media (Sartorius) with 5% human serum (Sigma Aldrich).
3D culture in collagen gels
Type I collagen gels were prepared using a modified protocol from Artym and Matsumoto . Briefly, cells were resuspended in a mix of rat tail collagen type I (Corning), 0.02 N acetic acid, 10x DMEM with Phenol Red (Sigma Aldrich) and 10x reconstitution buffer (0 .2M Hepes (Gibco) and 0.2 62 M NaHCO3). To neutralize the pH, 2 N NaOH was added to the reconstitution buffer before use. For all cell lines, low density (LD) and high density (HD) gels contained 1 mg/ml and 4 mg/ml collagen type I respectively. First, 300 to 400 μl of the collagen solution was plated per well of a non-tissue culture treated 24-wells plate (Corning) and allowed to polymerize at 37 °C for at least 30 min. Afterwards a second layer of collagen solution containing the cells was seeded on top of the first gel and allowed to polymerize for 1–2 h at 37 °C, after which 600 μl of culture media was added. For culturing cells in 2D conditions, the same number of cells were seeded in the wells of a tissue culture treated 24-well plate.
T cell proliferation assay
T cells from healthy donors were transiently stimulated with 10 nM PMA (Sigma Aldrich) and 175 nM ionomycin (Sigma Aldrich) and labelled with CellTrace Violet dye (Thermo Fisher Scientific). 8 × 105 T cells were seeded in each well within LD or HD collagen gels or on tissue culture plastic. After 5 days of culture in 3D collagen gels or on plastic (2D), cells were treated with 3 mg/ml collagenase (Worthington) solution for 45–60 min at 37 °C to extract the cells from the gels, washed once with media and once with DPBS and resuspended in FACS buffer containing the following: Live/Dead Fixable Near-IR Dead cell stain (Thermo Fisher Scientific), anti-CD3-APC (cl. SK7), anti-CD4-PE (cl. SK3), and anti-CD8-FITC (cl. HIT8a) (all BD Biosciences). Cells were incubated with the antibodies in the dark at 4 °C for 30 min, washed twice with FACS buffer, resuspended in FACS buffer and acquired using a BD FACSCanto II flow cytometer (BD Biosciences). Analysis was performed with FlowJo V10 software. Experiments were repeated three times using T cells isolated from different donors.
T cell proliferation was also measured using the APC BrdU Flow Kit (BD Biosciences) according to the manufacturer’s instructions. Briefly, 8 × 105 T cells were seeded in each well within LD or HD collagen gels or on tissue culture plastic for 3 days. For labelling, BrdU was added to the medium of each well (final concentration 10 μM) and cells were incubated for 90 min at 37 °C in 5% CO2. Afterwards, media of the wells with collagen gels was aspirated and replaced with 3 mg/ml collagenase solution (Worthington). Collagenase was also added to the media of 2D cultured cells in concentration similar to which cells in 3D gels were exposed. After complete digestion of the collagen gels, cells were collected and washed once with DPBS (Lonza) and stained with Live/Dead Fixable Near-IR Dead cell stain (Thermo Fisher Scientific), anti-CD3-FITC (cl- UCHT1), anti-CD4-PE (cl. SK3), and anti-CD8-BV421 (cl. RPA-T8) (all BD Biosciences). Next, cells were fixed, permeabilized, and stained with APC anti-BrdU antibody according to the manufacturer’s instructions. Cells were resuspended in FACS buffer and acquired using a BD FACSCanto II or LSR II flow cytometer (BD Biosciences). Analysis was performed with FlowJo V10 software.
Cancer cell proliferation assay
Cell proliferation was determined using the APC BrdU Flow Kit (BD Biosciences). 50.000 cells were cultured in the different 3D or 2D culture conditions for 5 days. For labelling, BrdU was added to the medium of each well (final concentration 10 μM) and cells were incubated for 45 min at 37 °C in 5% CO2. After digestion of collagen gels and washing, cells were stained with Zombie Aqua Fixable Viability Dye (BioLegend) to determine cell viability and analyzed for the fraction of BrdU positive cells as described above.
RNA extraction, cDNA synthesis and quantitative real-time-PCR
For RNA isolation, 8 × 105 PBMCs enriched for T cells, purified T cells (isolated with anti-CD3 microbeads, Miltenyi Biotec), or purified CD8 positive cells (isolated with CD8 MicroBeads, Miltenyi Biotec) were transiently stimulated with PMA and ionomycin and seeded within LD or HD collagen gels or on tissue culture plastic. After 2 days, total RNA from cell cultures was purified using RNeasy kit (Qiagen) according to the manufacturer’s instructions. Quality of samples was measured using an Agilent 2100 BioAnalyzer (Agilent Genomics). Afterwards 500 ng – 1 μg RNA per sample was reverse transcribed using iScript cDNA Synthesis Kit (Bio-Rad). The synthesized cDNA was used as template in a real-time quantitative PCR reaction using Brilliant III Ultra-Fast SYBR Green (Agilent Technologies) according to the manufacturer’s standard protocol. Equal amounts of cDNA were applied in each reaction mixture. As a control for the specificity of the quantitative real-time PCR, a sample without template was included. The real-time cycler conditions were as follows: initial activation step at 95 °C for 3 min, 40 cycles of denaturing at 95 °C for 5 s, and annealing/extension at 60 °C for 20 s, followed by a melting curve analysis of 65–95 °C with 0 .5°C increment, 5 s per step.
All measurements were based on triplicates or quadruplicates of each cell culture condition measured in duplicates and normalized to the internal control gene, ACTB. Four independent experiments were performed. The comparative cycle threshold (ΔΔCT) method was used to calculate the relative fold changes.
Primers were designed using the Primer-BLAST tool (NCBI, NIH). All primers spanned exon-exon junction, and maximum product length was 250 bp. Primer efficiencies were measured for all primer sets and found to be between 85 and 103%. Primers are listed in Additional file 1: Table S2.
T cell cultures from three different healthy donors were established and transduced with high affinity MAGE-A3a3a, MAGE-A3a3c, or MART-1 TCR, respectively, as previously described . For MART-1 TCR transduction, the TCRα and β sequences of the HLA-A2-restricted MART-1-specific (MART_27–36: ELAGIGILTV) A42 T cell clone were used .
1 × 106 TCR transduced T cells were embedded in collagen of high (4 mg/ml) or low (1 mg/ml) density together with 1 × 105 FM82 melanoma cells  and cultured for 2 days, after which total RNA was purified. Extensive T cell-mediated melanoma cytotoxicity during the co-culture period resulted in the presence of even fewer melanoma cells as compared to the starting point.
The quantity and purity of isolated RNA were assessed using an Agilent 2100 BioAnalyzer (Agilent Genomics). Total RNA (500 ng) was prepared for sequencing using polydT-mediated cDNA synthesis in accordance with the manufacturer’s (Illumina) instructions. Libraries were made with a NEBNext RNA Library Preparation Kit for Illumina. Library quality was assessed using Fragment Analyzer (AATI), followed by library quantification (Illumina Library Quantification Kit). Sequencing was done on a HiSeq1500 platform (Illumina) with a read length of 50 bp. Sequenced reads were aligned to the human genome assembly hg19 using STAR . Uniquely aligned reads were quantified at exons of annotated genes and normalized to sequence depth and gene length using HOMER . Sequencing depth and alignment information is in Additional file 1: Table S3. The number of reads per kilobase per million mapped (RPKM) for all RefSeq annotated genes can be found in Additional file 2: Table S1. The analysis of differential expression was performed using DESeq2 package in R . Principal component analysis was performed using R (prcomp package). Heatmaps were generated from z-score normalized RPKM values using R (pheatmap package) on selected sets of genes. MA and Volcano plots were generated using R.
A 51Cr-release assay for T cell-mediated cytotoxicity was used to assess the cytotoxicity of tumor infiltrating T cells after 3D culture in collagen matrices of high or low density. Autologous melanoma cells MM33 were used as target cells . Effector cells (T cells) were pre-cultured on plastic or in 3D collagen cultures for 3 days, after which cells were treated with 3 mg/ml collagenase solution for 45–60 min to extract them from the gels, and washed twice with media. Typically, 5 × 105 target cells in 150 μl RPMI were labelled with 20 μl 51Cr (Perkin Elmer) in a 15 ml falcon tube at 37 °C for 1 h. After washing, 5 × 103 target cells per well were plated out in a 96-well plate (Corning) and T cells were added at various effector to target cell (E:T) ratios. Cells were incubated at 37 °C for 4 h. Next, the level of 51Cr in the supernatant was measured using a Wallac Wizard 1470 automatic γ counter (Perkin Elmer). The maximum 51Cr release was determined by addition of 100 μl 10% Triton X-100, and minimum release was determined by addition of 100 μl of RPMI to target cells. Specific lysis was calculated using the following formula: ((cpm sample – cpm minimum release)/(cpm maximum release – cpm minimum release)) × 100%.
All individual experiments were performed at least three times with at least three replicates per condition and results are presented as mean ± standard error of the mean (SEM) unless otherwise specified. For two-group comparisons between cells from the same donors, paired two-tailed Student’s t-tests were performed. For multi-group comparisons, one-way analysis of variance (ANOVA) was used followed by paired two-tailed Student’s t-tests. All statistical analyses were performed using GraphPad Prism. A p-value < 0.05 was considered statistically significant. Correlation between PSR-positive area and CD8+ cell abundance was analyzed by Pearson correlation.
Correlation of expression levels of collagen and T cell activation markers in oral cancer
To study the correlation between collagen type I (COL1A1) levels and the expression level of 7 markers of T cell activation, we used an already published single-cell RNA sequencing dataset available in the Gene Expression Omnibus (GSE103322). In the study by Puram et al. , head and neck tumors from 18 patients and 5 matching lymph nodes were disaggregated, sorted into single cells and profiled using Smart-seq2. Here, we used the already processed data discarding lymph node samples (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE103322), which contain expression levels of 23,686 genes in 4541 single-cells (including fibroblasts, B cells, T cells, dendritic cells, endothelial cells, macrophages, mast cells, myocytes, and cancer and non-cancer cells). To investigate correlation between collagen and T cell activation markers, we first calculated average expression levels of COL1A1 across all cell types and expression of each T cell marker across T cells for each patient. Then, we calculated pairwise Pearson correlation between average expression levels of COL1A1 and each T cell activation marker. All data processing and correlation analysis were performed using Pythons SciPy and Pandas .
Additional materials and methods
Detailed information about cancer cell culture, confocal microscopy, flow cytometry analysis of T cell subsets, histology, and ELISA can be found in the Additional file 3.
3D culture of T cells in different collagen densities impairs proliferation without compromising viability
To investigate if 3D culture in collagen matrices of different collagen concentrations affected the viability of T cells, we isolated T cells from healthy donors and transiently stimulated the cells with PMA and ionomycin. This type of stimulation bypasses T cell receptor activation but acts on several of the same downstream signaling pathways including Protein Kinase C . The T cells were embedded in collagen matrices of high (4 mg/ml) or low (1 mg/ml) collagen concentration, or seeded on regular tissue culture plastic (2D culture) and cultured for 5 days. The selected collagen concentration of 1 mg/ml is representative of healthy normal tissue such as lung or mammary gland whereas 4 mg/ml collagen gels mimic the tissue stiffening occurring in solid tumors [19, 47]. To completely avoid cellular contact with the plastic surface of the wells, the 3D culture was established on top of a pre-generated collagen matrix without cells (Fig. 1a). To evaluate if viability of the T cells was affected by the different culture conditions, cells were extracted from the collagen matrices by a brief collagenase-treatment, stained with a live/dead cell marker and analyzed by flow cytometry (Fig. 1b). A high viability of more than 95% was observed in both 2D culture and in 3D culture in different collagen densities. To visualize the 3D culture of T cells in collagen matrices of different collagen concentrations, purified T cells were embedded in collagen matrices and imaged by confocal fluorescent microscopy (Fig. 1c-e). As expected, 3D cultured T cells were completely surrounded by collagen (Fig. 1c). The structure and density of collagen fibers were clearly different in the matrices of different collagen concentrations but no apparent morphological changes were observed between T cells in 3D culture of high or low collagen density (Fig. 1d-e). To examine if T cell proliferation was affected by collagen density, T cells were transiently PMA/ionomycin stimulated, CellTrace Violet (CTV) labeled, and embedded in a high- or low-density collagen matrix or cultured on regular tissue culture plastic. Flow cytometry-based analysis of CTV dilution in CD3+ cells showed a clear reduction in proliferation when cells were cultured in 3D compared to 2D. Interestingly, we also observed a smaller but still significant reduction in proliferation when T cells were cultured in a high-density collagen matrix compared to a low-density collagen matrix (Fig. 1f-g). This result was confirmed using a BrdU-based proliferation assay (Fig. 1h). In consistence with other reports [21, 48] we did not observe that proliferation of cancer cell lines was similarly impaired in a high-density collagen matrix (Fig. 1i).
The effects were similar for CD4+ and CD8+ T cells (Additional file 1: Figure S1) and consistently the different culture conditions did not change the ratio of CD4+ cells to CD8+ cells after 2 days in culture (Fig. 1j). However, after prolonged 3D culture of T cells for 5 days, the ratio of CD4+ cells to CD8+ cells was higher in a high-density matrix compared to a low-density matrix (Fig. 1k), indicating that proliferation and/or survival is slightly more impaired for CD8+ cells than for CD4+ cells in a high-density collagen matrix. Analysis of the distribution of the different T cell differentiation subsets after 5 days of 3D culture in a high- or low-density collagen matrix or regular 2D culture (Additional file 1: Figure S2) suggested an increase in effector memory T cells and a decrease in central memory T cells after 3D culture in a high-density collagen matrix compared to a low-density collagen matrix.
To study if the reduced T cell proliferation in a high-density collagen matrix could be reflected in the differential abundance of tumor-infiltrating T cells in breast tumors, we examined 20 samples of resected triple-negative breast cancers, which were all of histological grade 3 and had a diameter between 10 and 20 mm. The samples were immunostained for CD8 and picrosirius red-stained for fibrillar collagen (Fig. 2a-d). Consistent with our 3D culture data, samples that contained a high collagen density often had fewer infiltrating CD8+ T cells, although this negative correlation failed to reach statistically significance with this limited material (Fig. 2e). Analysis of the number of CD8+ cells located in intratumoral regions of either high or low collagen density within each tumor sample (Additional file 1: Figure S3) also revealed a tendency of fewer CD8+ T cells in regions of high collagen density (Fig. 2f).
The gene expression profile of T cells is regulated by surrounding extracellular matrix
To further elucidate the response of T cells to 3D culture in different collagen densities, T cells were 3D cultured for 2 days in low- or high-density collagen matrices or cultured on regular tissue culture plastic and subjected to RNA sequencing. A principal component analysis shows that the gene expression profile of cells cultured in 3D (low and high density) separate very clearly from the 2D cultured cells (cultured on plastic) and additionally that cells cultured in a high-density collagen matrix cluster separately from cells cultured in a low-density collagen matrix (Fig. 3a). The full gene expression dataset is in Additional file 2: Table S1. The clear difference between 2D and 3D culture was reflected in 683 differentially regulated genes (FDR < 0.01 and fold change > 1.5) for cells cultured in a low-density collagen matrix compared to regular 2D culture on tissue culture plastic (Fig. 3b and Additional file 1: Figure S4A) and 1928 differentially regulated genes (FDR < 0.01 and fold change > 1.5) between culture in high-density collagen compared to regular 2D culture (Additional file 1: Figure S4B-C). In consistence with the reduced proliferation (Fig. 1f-h), downregulated genes were involved in cell cycle processes (Fig. 3c and Additional file 1: Figure S4E). To investigate if T cell activity was affected by 3D culture, we examined the expression levels of a panel of T cell activity markers (Fig. 3d), regulatory T cell (Treg) markers (Fig. 3e), and exhaustion markers (Fig. 3f). The heatmaps did not show any clear 3D culture-induced changes in T cell activity although some of the genes were indeed significantly regulated, suggesting that 3D culture could influence T cell biology. To identify putative transcription factors (TFs) and TF families that may be responsible for the differential gene expression observed after 3D culture of T cells, we used the computational method ISMARA (Integrated Motif Activity Response Analysis, https://ismara.unibas.ch/mara/) to model transcription factor activity . Among the top-ranked identified downregulated TF motifs, several known pro-proliferative factors such as Myb and members of the E2F family were included [50, 51] (Additional file 1: Figure S5). This is consistent with the gene ontology analysis (Fig. 3c) and the observed reduction in cellular proliferation after 3D culture of T cells (Fig. 1f-h).
T cells cultured in a high-density collagen matrix downregulate markers of cytotoxic activity
To investigate the gene regulation specifically induced by collagen density, we compared T cells cultured in high- vs. low-density collagen matrices and found that 351 genes were differentially expressed (FDR < 0.01 and fold change > 1.5) (Fig. 4a and Additional file 1: Figure S4D). In alignment with the observed collagen density-induced reduction in proliferation (Fig. 1f-h), downregulated genes were involved in cell cycle processes (Fig. 4b). Upregulated genes were involved in processes such as chemotaxis and cell migration (Fig. 4b). Of the 10 most significantly up- and down-regulated genes (Table 1), several have been suggested to impact T cell activity. CD101, which is upregulated in a high-density matrix, has been suggested to be involved in the negative regulation of T cell activity  and the expression levels on Tregs correlate with their immunosuppressive potency . CIP2A, which is downregulated in a high-density matrix, has been suggested to promote T cell activation .
To investigate the potential regulation of T cell activity by collagen density, we visualized the regulation of a panel of cytotoxic activity markers, exhaustion markers, and markers of Tregs (Fig. 4c-e). Transcripts encoding six out of seven T cell activity markers were significantly downregulated by culture in high-density collagen compared to low-density collagen (Fig. 4c). Only TNF did not follow the same trend and was instead significantly upregulated. At the same time Treg markers were slightly upregulated (Fig. 4d). The markers of T cell exhaustion were not systematically regulated by collagen density (Fig. 4e). To identify potential transcriptional regulators of these effects, we again used ISMARA to model TF activity. Among the top-ranked downregulated TF motifs, we identified several E2F motifs (Fig. 4f), which is consistent with the observed changes in cell proliferation. Among the top-ranked upregulated TF motifs, we identified SMAD4 and FOXO1 motifs (Fig. 4f), which are important mediators of TGF-β induced Treg differentiation [55, 56]. In addition, upregulated motifs include PML and FOXM1, which have been suggested to be important for sustaining TGF-β signaling [57, 58].
Altogether, the striking gene expression differences induced by culture in high-density collagen compared to low-density collagen suggests that the T cells acquire a less cytotoxic and more regulatory phenotype.
All of the RNA sequencing experiments were performed using PBMCs isolated from a single healthy donor. To examine if the collagen-density induced gene regulations were donor-specific or reflected a general regulatory mechanism, we examined if similar gene-regulation was observed in cells isolated from four additional healthy donors. A panel of six transcripts, which were significantly regulated by collagen-density in the RNA sequencing experiments (two upregulated and four downregulated), were analyzed by qRT-PCR and compared to the RNA sequencing data (Fig. 5a-b). The gene-regulations followed the same pattern as observed in the RNA sequencing experiments, with four of the six genes showing a statistically significant regulation. This experiment as well as the RNA sequencing was performed with PBMCs enriched for T cells by excluding the strongly adherent fraction of the PBMCs, and consequently approximately 70% of the cells were T cells (data not shown). To investigate if the remaining 30% cells, which include myeloid cells, could be critical for mediating this effect of collagen-density we purified T cells from PBMCs from four healthy donors using magnetic anti-CD3 microbeads, resulting in a T cell purity of more than 96% (data not shown). The purified T cells were again transiently PMA/ionomycin-stimulated and cultured in high-density collagen or low-density collagen. Analysis by qRT-PCR of the six-gene panel showed a similar pattern of gene-regulations, indicating that T cells are directly affected by the surrounding collagen density (Fig. 5c). The only gene that did not follow the same pattern in the purified T cell was TNFRSF9, potentially indicating the myeloid cells were involved in mediating the previously observed upregulation of TNFRSF9. Analysis of purified CD8+ T cells after culture in collagen of high or low density showed that cytotoxic CD8-positive T cells responded to the surrounding collagen matrix in a comparable manner as observed for total CD3+ T cells (Fig. 5d). To investigate if TCR-stimulated T cells responded similarly to the surrounding collagen density, we co-cultured MAGE-A3 or MART-1 TCR-transduced T cells together with the melanoma cell line FM82  within collagen matrices (Fig. 5e). Analysis of the same six-gene panel showed a similar response as observed for the purified T cells stimulated with PMA/ionomycin (Fig. 5c-d).
Collagen density modulates the tumoricidal activity of tumor-infiltrating T cells
To investigate if the collagen density-induced transcriptional regulation of T cells was also reflected in an altered cytotoxic activity, we used a matched set of cultured T cells and melanoma cells isolated from the same tumor fragment. The T cell culture from this patient (MM33) has previously been shown to contain tumor-reactive T cells and to have the ability to lyze autologous melanoma cells . The MM33 T cells were transiently PMA/ionomycin-stimulated and 3D cultured for 2 days in a low-density collagen matrix or a high-density collagen matrix or on regular tissue culture plastic. Subsequently, the cells were extracted from the matrices and assayed for their ability to lyze melanoma cells in a standard four-hour Cr-51 release assay (Fig. 6a). The incubation of melanoma cells with increasing numbers of T cells led to increased cell lysis, but the cytotoxic activity of the T cells was impaired after 3D culture compared to 2D culture. Strikingly, cytotoxicity was particularly low for T cells cultured in high-density collagen matrix (Fig. 6a-b). The collagen density-induced regulation of T cell cytotoxicity was also reflected in a reduced level of secreted IFN-γ for the MM33 T cells cultured in a high-density collagen matrix (Fig. 6c). In these experiments, stimulation with PMA/ionomycin led to noticeable cell death for a large fraction of the cells, which is probably due to the preceding expansion of the cells in high-dose IL-2 containing media. Therefore, we also 3D cultured the T cell for 2 days without any stimulation followed by extraction of cells for cytotoxicity measurements (Fig. 6d). In this situation we also observed a reduced cytotoxic activity of the T cells cultured in a high-density collagen matrix compared to the T cells cultured in a low-density collagen matrix (Fig. 6d-e). The reduced cytotoxicity was again accompanied by reduced secretion of IFN-γ, although it should be noted that IFN-γ levels were much lower than for the PMA/ionomycin treated T cells (Fig. 6f).
Taken together, our study reveals a novel immunosuppressive mechanism, which could be important for cancer progression and for cancer immunotherapy efficacy.
Collagen density within tumors and their T cell infiltration comprise strong prognostic indicators of poor and good prognosis, respectively. It has, however, not been investigated if these two parameters are independent or if they are interrelated. In this study, we used 3D culture to address whether T cells respond to the surrounding collagen matrix and if the collagen-density alters their activity.
A first important observation made was that T cells indeed respond to their extracellular matrix environment. This was reflected in a dramatically altered transcriptional profile and reduced proliferation for 3D cultured cells compared to regular 2D culture on tissue culture plastic. Gene ontology analysis to identify biological process categories that were statistically enriched confirmed that cell proliferation was significantly affected in 3D. This suggests that 3D culture models, which more accurately mimic tissue environments, could be highly relevant for studies of T cell biology. Although 3D culture of T cells, compared to regular 2D culture, led to a substantial change in the transcriptional profile, inspection of genes encoding markers of activation, Tregs, or exhaustion did not reveal any obvious pattern in the effects on T cell activity.
Another important observation in this study was that T cells clearly respond to the density of the surrounding collagen matrix. Culture of T cells within collagen-matrices of either high or low density resulted in fewer differentially regulated genes compared to the culture of T cells in 2D vs. 3D but nevertheless resulted in striking functional differences. A high collagen-density reduced T cell proliferation and 5 days of culture in this matrix favored CD4+ T cells over CD8+ T cells. A similar reduction in cell proliferation was not observed for three different breast cancer cell lines, suggesting that within a tumor of high collagen density, T cell proliferation (especially CD8+ T cell proliferation) is impaired whereas cancer cells are unaffected. In line with these observations, we observed in a limited number of triple-negative breast cancer samples, that collagen density seemed to negatively impact CD8+ T cell abundance. The effect of collagen density on T cell proliferation could constitute a new immunosuppressive mechanism within the tumor microenvironment and provide an explanation for the correlation between collagen density in tumors and cancer patient prognosis.
In alignment with the effect of the surrounding collagen density on T cell proliferation, inspection of differentially regulated genes confirmed that cell cycle processes were affected. In addition, markers of cytotoxic T cell activity were clearly downregulated by a high-density collagen matrix compared to a low-density matrix whereas markers of Tregs were upregulated. This striking observation suggests that collagen density, in addition to reducing T cell proliferation, impairs cytotoxic activity. Analysis of putative transcription factor motifs, which could mediate the observed transcriptional changes, identify decreased activity of E2F motifs and increased activity of transcription factor motifs downstream of TGF-β signaling. These findings suggest that E2F transcription factors could be involved in the collagen density-induced inhibition of T cell proliferation and that autocrine TGF-β signaling could be centrally engaged in the modulation of T cell activity.
Using a unique T cell culture and melanoma cell culture established from the same resected melanoma, we could also demonstrate that the tumor-infiltrating T cells were indeed less efficient at killing the melanoma cells after culture in a high-density collagen matrix compared to a low-density collagen matrix. This observation suggests that the collagen-density of the tumor microenvironment can support the cancer cells’ escape from immune destruction by reducing T cell proliferation and by modulating the cytotoxic activity of tumor-infiltrating T cells. This identified immunosuppressive mechanism could be of relevance during tumor progression but also have importance for the efficacy of cancer immunotherapy.
Although this is the first study to directly assess the response of T cells to the surrounding collagen density, the potential of collagen to modulate immune activity is supported by a study of tissue regeneration, in which collagen implantation in wounded muscles of mice was shown to induce an immunosuppressive microenvironment . This effect involved M2-polarization of macrophages, which led to Th2-polarization of T cells. In our study we focused on the ability of the surrounding collagen to directly regulate T cell activity, but it is possible that collagen-induced M2-polarization of macrophages could further augment this modulation of T cells.
In addition to the effects of collagen-density on T cell activity observed in this study, others have suggested that stromal collagen could limit the migration of T cells into the tumor islets and thereby impede their contact with cancer cells [35, 36]. The impaired T cell migration into tumor islets was suggested to be caused by reduced motility in collagen-dense region combined with inappropriate guidance of T cells by the collagen fibers aligned in parallel to tumor islets. These studies further underscore that collagen within the tumor microenvironment could be an important regulator of anti-cancer immunity.
By using 3D culture of T cells, we have identified collagen density as a novel regulator of anti-cancer T cell activity. This immunosuppressive mechanism could be of central importance for the cancer cells’ evasion of immune-destruction and could constitute a novel therapeutic target for enhancing immunotherapy efficacy.
Cell trace violet
Myeloid-derived suppressor cells
Peripheral blood mononuclear cells
Programmed cell death protein 1
Programmed death-ligand 1
Phorbol 12-myristate 13-acetate
T cell receptor
Transforming growth factor β
Regulatory T cells
Fridman WH, Zitvogel L, Sautès–Fridman C, Kroemer G. The immune contexture in cancer prognosis and treatment. Nat Rev Clin Oncol. 2017;14:717–34.
Tumeh PC, Harview CL, Yearley JH, Shintaku IP, Taylor EJM, Robert L, et al. PD-1 blockade induces responses by inhibiting adaptive immune resistance. Nature. 2014;515:568–71.
Dunn GP, Bruce AT, Ikeda H, Old LJ, Schreiber RD. Cancer immunoediting: from immuno-surveillance to tumor escape. Nat Immunol. 2002;3:991–8.
Munn DH, Bronte V. Immune suppressive mechanisms in the tumor microenvironment. Curr Opin Immunol. 2016;39:1–6.
Joyce JA, Fearon DT. T cell exclusion, immune privilege, and the tumor microenvironment. Science. 2015;348:74–80.
Kumar V, Patel S, Tcyganov E, Gabrilovich DI. The nature of myeloid-derived suppressor cells in the tumor microenvironment. Trends Immunol. 2016;37:208–20.
Sharma P, Allison JP. The future of immune checkpoint therapy. Science. 2015;348:56–61.
Madsen DH, Bugge TH. The source of matrix-degrading enzymes in human cancer: problems of research reproducibility and possible solutions. J Cell Biol. 2015;209:195–8.
Madsen DH, Jürgensen HJ, Siersbæk MS, Kuczek DE, Grey Cloud L, Liu S, et al. Tumor-associated macrophages derived from circulating inflammatory monocytes degrade collagen through cellular uptake. Cell Rep. 2017;21:3662–71.
Kessenbrock K, Plaks V, Werb Z. Matrix metalloproteinases: regulators of the tumor microenvironment. Cell. 2010;141:52–67.
Schedin P, Keely PJ. Mammary gland ECM remodeling, stiffness, and mechanosignaling in normal development and tumor progression. Cold Spring Harb Perspect Biol. 2011;3:a003228.
Cox TR, Erler JT. Molecular pathways: connecting fibrosis and solid tumor metastasis. Clin Cancer Res. 2014;20:3637–43.
Naba A, Clauser KR, Lamar JM, Carr SA, Hynes RO. Extracellular matrix signatures of human mammary carcinoma identify novel metastasis promoters. Elife. 2014;3:e01308.
Cox TR, Bird D, Baker A-M, Barker HE, Ho MW-Y, Lang G, et al. LOX-mediated collagen crosslinking is responsible for fibrosis-enhanced metastasis. Cancer Res. 2013;73:1721–32.
Pickup MW, Mouw JK, Weaver VM. The extracellular matrix modulates the hallmarks of cancer. EMBO Rep. 2014;15:1243–53.
Conklin MW, Eickhoff JC, Riching KM, Pehlke CA, Eliceiri KW, Provenzano PP, et al. Aligned collagen is a prognostic signature for survival in human breast carcinoma. Am J Pathol. 2011;178:1221–32.
Ohno S, Tachibana M, Fujii T, Ueda S, Kubota H, Nagasue N. Role of stromal collagen in immunomodulation and prognosis of advanced gastric carcinoma. Int J Cancer. 2002;97:770–4.
Li H-X, Zheng J-H, Fan H-X, Li H-P, Gao Z-X, Chen D. Expression of αvβ6 integrin and collagen fibre in oral squamous cell carcinoma: association with clinical outcomes and prognostic implications. J Oral Pathol Med. 2013;42:547–56.
Paszek MJ, Zahir N, Johnson KR, Lakins JN, Rozenberg GI, Gefen A, et al. Tensional homeostasis and the malignant phenotype. Cancer Cell. 2005;8:241–54.
Levental KR, Yu H, Kass L, Lakins JN, Egeblad M, Erler JT, et al. Matrix crosslinking forces tumor progression by enhancing integrin signaling. Cell. 2009;139:891–906.
Provenzano PP, Inman DR, Eliceiri KW, Knittel JG, Yan L, Rueden CT, et al. Collagen density promotes mammary tumor initiation and progression. BMC Med. 2008;6:11.
Puig M, Lugo R, Gabasa M, Giménez A, Velásquez A, Galgoczy R, et al. Matrix stiffening and β1 integrin drive subtype-specific fibroblast accumulation in lung cancer. Mol Cancer Res. 2015;13:161–73.
Engler AJ, Sen S, Sweeney HL, Discher DE. Matrix elasticity directs stem cell lineage specification. Cell. 2006;126:677–89.
Mammoto A, Connor KM, Mammoto T, Yung CW, Huh D, Aderman CM, et al. A mechanosensitive transcriptional mechanism that controls angiogenesis. Nature. 2009;457:1103–8.
Leight JL, Wozniak MA, Chen S, Lynch ML, Chen CS. Matrix rigidity regulates a switch between TGF-β1-induced apoptosis and epithelial-mesenchymal transition. Mol Biol Cell. 2012;23:781–91.
Mushtaq MU, Papadas A, Pagenkopf A, Flietner E, Morrow Z, Chaudhary SG, et al. Tumor matrix remodeling and novel immunotherapies: the promise of matrix-derived immune biomarkers. J Immunother Cancer. 2018;6:65.
O’Connor RS, Hao X, Shen K, Bashour K, Akimova T, Hancock WW, et al. Substrate rigidity regulates human T cell activation and proliferation. J Immunol. 2012;189:1330–9.
Feng Y, Reinherz EL, Lang MJ. αβ T cell receptor Mechanosensing forces out serial engagement. Trends Immunol. 2018;39:596–609.
McMahon M, Ye S, Izzard L, Dlugolenski D, Tripp RA, Bean AGD, et al. ADAMTS5 is a critical regulator of virus-specific T cell immunity. PLoS Biol. 2016;14:e1002580.
Jachetti E, Caputo S, Mazzoleni S, Brambillasca CS, Parigi SM, Grioni M, et al. Tenascin-C protects cancer stem-like cells from immune surveillance by arresting T-cell activation. Cancer Res. 2015;75:2095–108.
Sangaletti S, Chiodoni C, Tripodo C, Colombo MP. Common extracellular matrix regulation of myeloid cell activity in the bone marrow and tumor microenvironments. Cancer Immunol Immunother. 2017;66:1059–67.
Yamada KM, Cukierman E. Modeling tissue morphogenesis and cancer in 3D. Cell. 2007;130:601–10.
Klemke M, Kramer E, Konstandin MH, Wabnitz GH, Samstag Y. An MEK-cofilin signalling module controls migration of human T cells in 3D but not 2D environments. EMBO J. 2010;29:2915–29.
Wolf K, te Lindert M, Krause M, Alexander S, te Riet J, Willis AL, et al. Physical limits of cell migration: control by ECM space and nuclear deformation and tuning by proteolysis and traction force. J Cell Biol. 2013;201:1069–84.
Salmon H, Franciszkiewicz K, Damotte D, Dieu-Nosjean M-C, Validire P, Trautmann A, et al. Matrix architecture defines the preferential localization and migration of T cells into the stroma of human lung tumors. J Clin Invest. 2012;122:899–910.
Hartmann N, Giese NA, Giese T, Poschke I, Offringa R, Werner J, et al. Prevailing role of contact guidance in intrastromal T-cell trapping in human pancreatic cancer. Clin Cancer Res. 2014;20:3422–33.
Artym VV, Matsumoto K. Imaging cells in three-dimensional collagen matrix. Curr Protoc Cell Biol. 2010;48:10.18.1–10.18.20.
Idorn M, Olsen M, Halldórsdóttir HR, Skadborg SK, Pedersen M, Høgdall C, et al. Improved migration of tumor ascites lymphocytes to ovarian cancer microenvironment by CXCR2 transduction. Oncoimmunology. 2018;7:e1412029.
Cole DJ, Weil DP, Shamamian P, Rivoltini L, Kawakami Y, Topalian S, et al. Identification of MART-1-specific T-cell receptors: T cells utilizing distinct T-cell receptor variable and joining regions recognize the same tumor epitope. Cancer Res. 1994;54:5265–8.
Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29:15–21.
Heinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P, et al. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol Cell. 2010;38:576–89.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.
Jespersen H, Lindberg MF, Donia M, Söderberg EMV, Andersen R, Keller U, et al. Clinical responses to adoptive T-cell transfer can be modeled in an autologous immune-humanized mouse model. Nat Commun. 2017;8:707.
Puram SV, Tirosh I, Parikh AS, Patel AP, Yizhak K, Gillespie S, et al. Single-cell transcriptomic analysis of primary and metastatic tumor ecosystems in head and neck Cancer. Cell. 2017;171:1611–1624.e24.
Oliphant TE. Python for scientific computing. Comput Sci Eng. 2007;9:10–20.
Isakov N, Altman A. Protein kinase Cθ in t cell activation. Annu Rev Immunol. 2002;20:761–94.
Wullkopf L, West A-KV, Leijnse N, Cox TR, Madsen CD, Oddershede LB, et al. Cancer cell ability to mechanically adjust to extracellular matrix stiffness correlates with their invasive potential. Mol Biol Cell. 2018;29(20):2378–85 mbc. E18-05-0319.
Provenzano PP, Inman DR, Eliceiri KW, Keely PJ. Matrix density-induced mechanoregulation of breast cell phenotype, signaling and gene expression through a FAK–ERK linkage. Oncogene. 2009;28:4326–43.
Balwierz PJ, Pachkov M, Arnold P, Gruber AJ, Zavolan M, van Nimwegen E. ISMARA: automated modeling of genomic signals as a democracy of regulatory motifs. Genome Res. 2014;24:869–84.
Oh I-H, Reddy EP. The myb gene family in cell growth, differentiation and apoptosis. Oncogene. 1999;18:3017–33.
Bertoli C, Skotheim JM, de Bruin RAM. Control of cell cycle transcription during G1 and S phases. Nat Rev Mol Cell Biol. 2013;14:518–28.
Rivas A, Ruegg CL, Zeitung J, Laus R, Warnke R, Benike C, et al. V7, a novel leukocyte surface protein that participates in T cell activation. I. Tissue distribution and functional studies. J Immunol. 1995;154:4423–33.
Fernandez I, Zeiser R, Karsunky H, Kambham N, Beilhack A, Soderstrom K, et al. CD101 surface expression discriminates potency among murine FoxP3+ regulatory T cells. J Immunol. 2007;179:2808–14.
Côme C, Cvrljevic A, Khan MM, Treise I, Adler T, Aguilar-Pimentel JA, et al. CIP2A promotes T-cell activation and immune response to listeria monocytogenes infection. PLoS One. 2016;11:e0152996.
Shi Y, Massagué J. Mechanisms of TGF-beta signaling from cell membrane to the nucleus. Cell. 2003;113:685–700.
Ouyang W, Beckett O, Ma Q, Paik J, DePinho RA, Li MO. Foxo proteins cooperatively control the differentiation of Foxp3+ regulatory T cells. Nat Immunol. 2010;11:618–27.
Lin H-K, Bergmann S, Pandolfi PP. Cytoplasmic PML function in TGF-β signalling. Nature. 2004;431:205–11.
Xue J, Lin X, Chiu W-T, Chen Y-H, Yu G, Liu M, et al. Sustained activation of SMAD3/SMAD4 by FOXM1 promotes TGF-β–dependent cancer metastasis. J Clin Invest. 2014;124:564–79.
Sadtler K, Estrellas K, Allen BW, Wolf MT, Fan H, Tam AJ, et al. Developing a pro-regenerative biomaterial scaffold microenvironment requires T helper 2 cells. Science (80- ). 2016;352:366–70.
We thank Drs. Thomas H. Bugge and Marie Kveiborg for critically reviewing this manuscript. We thank Dr. Alberto Santos for his help with bioinformatics analyses. We would like to acknowledge Dr. Andrew Gerry and Dr. Bent Jakobsen of Adaptimmune Ltd. (Oxfordshire, UK) for supplying us with lentiviral vectors for high affinity MAGE-A3 TCR transduction.
This study was supported by the Danish Cancer Society (D.H.M.), Knæk Cancer (D.H.M.), Novo Nordisk Foundation (D.H.M.), Dagmar Marshalls Foundation (D.E.K., D.H.M.), Dansk Kræftforskningsfond (D.E.K.), Einar Willumsen Foundation (D.E.K.), and Immutrain Project, EU’s H2020 under Marie-Sklodowska-Curie grant agreement 641549 (A.M.C.S.).
Availability of data and materials
The datasets obtained in the current study are available from the corresponding author or available from the indicated sources.
Ethics approval and consent to participate
Healthy donor PBLs were obtained from buffy coats available from the central blood bank of the capital region of Copenhagen and informed consent was obtained from all donors. Breast cancer samples were anonymized and used upon approval by the Scientific Ethics Committee for The Capital Region of Denmark.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Proliferation of CDLA and CD8A T cells. Figure S2. Fraction of T cell subsets after HD culture or %D culture in different collagen densities. Figure S3. Analysis of CD8+ T cell abundance in areas of high and low collagen density in 20 triple negative breast cancer samples. Figure S4. RNAseq data. Figure S5. Regulated transcription factors after 2D culture or 3D culture in low density collagen. Table S2. Sequences of primers used for RT-qPCR. All primers were designed using the NCBI gene database and the primer-BLAST tool. Table S3. RNA sequencing depth and alignment info. (PDF 1041 kb)
The number of reads per kilobase per million mapped (RPKM) for all RefSeq annotated genes. (XLSX 5563 kb)
Supplementary materials and methods. (PDF 91 kb)