- 1Cancer Immunology Group, School of Medicine, University of Nottingham Biodiscovery Institute, Nottingham, United Kingdom
- 2Faculty of Pharmacy, Al-Zaytoonah University of Jordan, Amman, Jordan
- 3Academic Unit of Translational Medical Sciences, School of Medicine, Queens Medical Centre, University of Nottingham, Nottingham, United Kingdom
- 4School of Veterinary Medicine and Sciences, University of Nottingham Biodiscovery Institute, Nottingham, United Kingdom
- 5Department of Pharmacology, Weill Cornell Medicine, New York, NY, United States
- 6Birmingham Tissue Analytics, College of Medical and Dental Sciences, University of Birmingham, Birmingham, United Kingdom
- 7School of Veterinary Medicine and Science, University of Nottingham, Sutton Bonington, United Kingdom
- 8Faculty of Medicine, Institute of Biomedical Sciences, Vilnius University, Vilnius, Lithuania
- 9Institute of Immunology and Immunotherapy, University of Birmingham, Birmingham, United Kingdom
- 10Host-Tumour Interactions Group, School of Medicine, University of Nottingham Biodiscovery Institute, Nottingham, United Kingdom
Introduction: Characterization of the tumour immune infiltrate (notably CD8+ T-cells) has strong predictive survival value for cancer patients. Quantification of CD8 T-cells alone cannot determine antigenic experience, as not all infiltrating T-cells recognize tumour antigens. Activated tumour-specific tissue resident memory CD8 T-cells (TRM) can be defined by the co-express of CD103, CD39 and CD8. We investigated the hypothesis that the abundance and localization of TRM provides a higher-resolution route to patient stratification.
Methods: A comprehensive series of 1000 colorectal cancer (CRC) were arrayed on a tissue microarray, with representative cores from three tumour locations and the adjacent normal mucosa. Using multiplex immunohistochemistry we quantified and determined the localization of TRM.
Results: Across all patients, activated TRM were an independent predictor of survival, and superior to CD8 alone. Patients with the best survival had immune-hot tumours heavily infiltrated throughout with activated TRM. Interestingly, differences between right- and left-sided tumours were apparent. In left-sided CRC, only the presence of activated TRM (and not CD8 alone) was prognostically significant. Patients with low numbers of activated TRM cells had a poor prognosis even with high CD8 T-cell infiltration. In contrast, in right-sided CRC, high CD8 T-cell infiltration with low numbers of activated TRM was a good prognosis.
Conclusion: The presence of high intra-tumoural CD8 T-cells alone is not a predictor of survival in left-sided CRC and potentially risks under treatment of patients. Measuring both high tumour-associated TRM and total CD8 T-cells in left-sided disease has the potential to minimize current under-treatment of patients. The challenge will be to design immunotherapies, for left-sided CRC patients with high CD8 T-cells and low activate TRM,that result in effective immune responses and thereby improve patient survival.
Introduction
Colorectal cancers (CRC) have traditionally been grouped together as one disease due to the anatomic continuity of the colon into the rectum. However, right and left colon differ in terms of their embryonic origin, vascular and nervous supplies, and gut flora (1). Increasing evidence has pointed to the location of the tumour affecting cancer pathology, progression and prognosis (2, 3), and ultimately patient’s response to different cancer treatments (4, 5). This in part may be due to the molecular variations and difference in mutational profiles between right and left-sided colon cancer (3, 6, 7). There is, however, increasing recognition in the literature that the differences in CRC may go beyond currently defined molecular subtypes (7) and that for cancer treatment right and left-sided colon cancer should be treated as separate diseases (3).
It is generally accepted that tumour progression may be influenced by non-malignant cells found in the tumour environment, especially the immune cells (8). T-cell infiltration has been shown to be superior to tumour staging as a prognostic factor (9, 10). Furthermore, responses to immune checkpoint blockade have been shown to be associated with patients having a pre-existing anti-tumour T-cell repertoire (11, 12). There is a strong correlation between the mutational load of a cancer type, the presence of neoepitopes and the response to immune check point blockade (13–15). Currently, only CRC patients with microsatellite instability (MSI), arising through loss of DNA mismatch repair function, are treated with checkpoint inhibitors (14). MSI is predominately found in right-sided colon cancers and tumours with MSI have been shown to have high T-cell infiltration (16, 17). Reflecting this, higher CD8 gene expression has also been demonstrated in right-sided colon cancer compared to the left-side (18).
Based on immune infiltration, tumours have been classified as: “immune-hot” (heavily infiltrated by lymphocytes); “immune-cold” (low levels of lymphocyte infiltration) (19). Since “immune-hot” tumours generally have a better prognosis than immune-cold tumours, this would imply that these immune cells – especially the CD8 T-cells – are activated and have anti-tumour activity. However, it has been shown that some tumour infiltrating T-cells may be inactive bystander cells which recognize viral antigens rather than cancer antigens (20).
Discriminating tumour-specific T-cells from bystander T-cells using markers of activation could aid in defining their role in the tumour microenvironment and patient outcome. One possible marker of activated CD8 T-cells is the co-expression of CD39 and CD103 (21). CD39 is an ectonucleotidase present on activated CD4 and CD8 T-cells following TCR stimulation (22). CD8+CD103+ have been defined as a tissue resident memory T-cell (TRM) signature (23). TRM have been described as T-cells that no-longer circulate but develop and reside in the peripheral tissues as part of the memory response (24). CD103 (integrin alpha E), binds to E-cadherin and is important in retaining TRM in the peripheral tissues (25).
The presence of CD103+ TIL in ovarian (26, 27), breast (26) lung (23, 28) and head and neck cancers (21) has shown a stronger correlation with survival than canonical T-cell markers such as CD3 or CD8. However, in a study in head and neck cancer, Duhen et al. (21) demonstrated that not all CD8+CD103+ T-cells were tumour-specific, and that some were the result of bystander recruitment. An increase in 4-1BB and Ki-67 expression on CD8+CD103+CD39+ T-cells lead the authors to hypothesize that these cells had recently encountered cognate antigen and were proliferating in the tumour. Moreover, these cells expressed granzyme and were cytotoxic. These tumour specific T-cells had an oligoclonal T-cell receptor (TCR) repertoire (21), and were associated with improved survival. Similarly, T-cells recognizing neoepitopes isolated from CRC patients expressed CD103 and CD39 (29), whereas bystander TRM that were specific for viral antigens and not tumour antigens had diverse phenotypes that lacked CD39 expression (20). There is therefore substantial evidence that co-expression of CD103/CD39 can be used as a marker of activated TRM in the tumour environment. Though these studies investigated tumour infiltrating lymphocytes they did not analyze the cells in situ. We hypothesized that intra-tumoural activated TRM confer a survival advantage (by virtue of their tumour-specificity) in patients with CRC. Using multiplex immunohistochemistry, on a tissue microarray (TMA) containing cores from 1000 CRC patients, we quantified the activated and non-activated intra-tumoural CD8 T-cells in 891 cases. This study demonstrated that while high numbers of activated TRM is a good predictor of survival in the overall cohort, there is a difference in the size of the effect depending on tumour location. In left-sided colon cancer high numbers of activated TRM was the sole marker of enhanced survival irrespective of total CD8 T-cell infiltration. In contrast, in the right-sided colon cancer, both high numbers of activated TRM and high total CD8 T-cell infiltration (even with low numbers of activated TRM) predicted good survival. This suggests that there is a fundamental difference in biology between right and left-sided colon cancers.
Materials and methods
Patient cohorts
Cohort 1 consisted of a consecutive series of 1000 CRC patients presenting at Nottingham University Hospitals NHS trust between 2008 and 2012 (53.6 months mean follow up). Cancer specific survival was measured from the date of primary surgical treatment to time of death due to cancer. Ethical approval was obtained for the study (reference no. 05/Q1605/66). The median age was 69 (range 16-94) with TMN stages: 16% stage I, 40% stage II, 32% stage III and 12% stage IV. For clinicopathological details refer to Table 1; Supplementary Table 1. Cohort 2 was obtained from the TCGA database containing mRNA sequencing data and clinicopathological characteristics (Table 1) of 515 colon adenocarcinomas (TCGA-COAD) (Cancer genome Atlas).
Tissue microarrays and immunohistochemistry
All tumours were formalin-fixed and processed into paraffin blocks. The histology of all cases was reviewed, appropriate donor blocks selected and marked to enable sampling of tumours from the luminal surface, center, advancing edge and adjacent normal mucosa. Tissue cores (0.6mm) were obtained from representative tumour regions of each donor block and arrayed into new recipient paraffin block using a tissue microarrayer (Beecher Instruments). Primary antibodies concentrations were optimized using chromogenic immunohistochemistry (30).
Multiplex immunohistochemistry
Multiplex immunohistochemistry was performed using the Opal 4-colour manual IHC kit (Akoya Bioscience) for simultaneous detection of CD8, CD103 and CD39 on a single TMA section. Extensive optimization achieved fluorescence staining with optimal intensity and no bleed-through into other channels. Slides were deparaffinized in xylene and rehydrated in three baths of 100%, 90%, and 70% ethanol. Antigen retrieval was performed with EDTA (pH:9, 35 minutes, 100°C). TMA sections were incubated overnight (4°C) with anti-CD39 [Abcam, EPR20627, 1:300 in antibody diluent/blocking solution], washed with TBST buffer, incubated with anti-mouse-HRP conjugate, washed and incubated with Opal 570 (1:100) for 10 minutes. Antibody stripping was performed by microwaving (100W microwave: 100% for 45 seconds, 20% for 14 minutes) in antigen retrieval buffer (AR9, Akoya Bioscience) (31). Sides were cooled, washed and the procedure repeated with anti-CD103 antibody [Abcam: EPR4166(2) 1 hour, room temperature] paired with Opal 690 (1:100) and then with anti-CD8 [Agilent DAKO clone: C8/144b (M7103) 1 hour, room temperature] paired with Opal 520 (1:100). Slides were counterstained with Spectra DAPI, and coverslips were mounted with ProLong® Diamond Antifade Mountant (Thermo Fisher Scientific). To obtain a spectral library anti-CD8 antibody was paired with each Opal fluorophore (without DAPI counterstain) as a singleplex.
Image acquisition and scoring
TMA slides were scanned (Vectra 3, Akoya Biosciences) at 20 times resolution. The quality of the spectral library was evaluated by reviewing the unmixed images to confirm the absence of spectral overlap or bleed-through between channels. Machine learning (Inform® software, Akoya Biosciences) enabled quantification of cells in the stroma and epithelium. This consisted of training the software to recognize: stroma and epithelial cells; individual cells and phenotype different cells based on fluorescence (Figure 1; Supplementary Figure 1). Tissue segmentation involved training the program to detect differences between the stroma and epithelial cells based on tissue morphologies. DAPI staining allowed nuclei detection which was used to determine the cell segmentation. One schema was used to phenotype all the different subsets. They were trained as following: CD8+CD103+CD39+; CD8+CD103+CD39-; CD8+CD103-CD39+; CD8+CD103-CD39-; CD8-CD103+CD39-; CD8-CD103-CD39+ and other (any other cell). This allowed training to ensure that all double and triple positive cells were as phenotyped and that membrane staining was around the whole cell and not fluorescence from a neighbouring cell (Supplementary Figure 2). A separate schema was used to identify total CD8+ T-cells. All cores were visually inspected to determine agreement with machine learning: any that did not pass initial quality control were trained again (Supplementary Figure 1). The main reason for rejection was tissue segmentation (Supplementary Figure 3 shows examples of failed and retrained images). Out of 4000 TMA cores, 3708 cores were accepted with the rest rejected due to loss of tissue or inaccurately trained after second training. The data from Inform® was merged and processed using phenoptrReports (Akoya Biosciences)
Figure 1 inform training stages for automated detection of tissues, cells and phenotyping individual cells (A) representive multiplex view of tumour cores with multi-stained cells, (B) Tissue segmentation of tumour epithelium (pink) and stroma (green) using pattern recognition (C) Cell segmentation was used to identify individual cell types (D) Training and phenotyping multi-stained and single stained cells was performed using a coloring- code method on infrom software. Black dots represent other cells.
Statistical analysis
Optimal cut-off points were determined using X-tile bioinformatics software (version 3.4.7). Statistical analysis was performed using the SPSS 21.0 software (SPSS, USA). Univariant and multivariant analyses were determined by chi-squared, log rank and Cox regression analysis, respectively. A p-value<0.05 was considered significant. GraphPad PRISM version 8 was used to compare the densities of different CD8+ T-cells. Statistical differences were assessed by two tailed paired T-test; One way Anova with Turkey’s multiple comparison test or Kruskal-Wallis with Dunn’s multiple comparison test.
TCGA analysis
Using the medical reports from the TCGA-COAD database, patients were sorted according to primary tumour location: right-sided (228 patients) or left-sided (131 patients). HTSeq data were used as the gene expression count. Differentially expressed genes were determined in each sample group using Deseq2 (32). These groups were defined in terms of CD8 and CD103 expression: high (4th Quartile) or low (1st Quartile). Inclusion criteria for differentially expressed genes (DEG) was an FDR < 0.05, Log2 FC < -1 and Log2 FC < +1. FDR were analyzed with ClueGO (33), Cytoscape (34) and Webgestalt (http://www.webgestalt.org/) (35).
Results
Multiplex immunohistochemistry
Multiplex IHC was performed on a TMA containing tumour cores from 1000 colorectal cancer patients (cohort 1), each of which was represented by four cores from different tumour regions (see Supplementary Figure 4): luminal side; center of the tumour; invasive margin and adjacent pathologically normal. Representative examples of fluorescently labelled cells are shown in Figure 2. The following subsets were classified: CD8+CD103+CD39+ (activated TRM) CD8+CD103+CD39-, CD8+CD103-CD39+, CD8+CD103-CD39-, and total CD8+ T-cells. Cores from 891 patients remained available to analysis due to the loss of cores during processing and quality control for machine learning.
Figure 2 Immunolocalization of different CD8+ phenotypes in colorectal cancer tissues. (A) representive monoplex expression of CD8, CD103 and CD39 markers, and the tissue nuclear counterstain (DAPI) using fluorescence setting (Inserts x40). (B) Multiplex staining of different CRC TMA cores presenting triple, double and single stained cells with CD8 (green), CD103 (red) and CD39 (yellow) biomarkers. Training and phenotyping multi-stained and single stained cells were performed using coloring- code method on inform software.
Activated TRM T-cells are an independent predictor of survival in CRC
Across the whole cohort, Kaplan Meier analysis indicated high numbers of activated TRM in the epithelium (Figure 3A, p=0.004) and stroma (Supplementary Figure 5, p=0.006) were associated with increased survival. Patients with high numbers of activated TRM in their tumours had a 10-year survival of 92% compared with 79% for total CD8+ T-cells (Figures 3A, B).
Figure 3 Prognostic impact of activated TRM T-cells in CRC patients. Kaplan-Meier plots represent the probability of disease-specific survival for (A) Activated TRM (CD8+CD103+CD39+) T-cells and (B) total CD8+ T-cells in the tumour epithelium. Cut-off points used to stratify CRC patients into high and low-density groups were determined using X-tile. The cut-off point for total CD8+ T cells was 222 cell/mm2 and for activated TRM cells was 74 cell/mm2. The log-rank test was used to compare curves, and p values <0.05 were considered statistically significant. (C) Multivariate analysis (Cox regression) of Intraepithelial total CD8 T cells with vascular invasion, TNM stage, tumour grade, microsatellite status and colorectal cancer-specific survival (D) Multivariate analysis (Cox regression) of Intraepithelial activated TRM with vascular invasion, TNM stage, tumour grade, microsatellite status and colorectal cancer-specific survival.
The abundant presence of high numbers of CD8 T-cells of any phenotype in the intra-epithelial region was significantly associated with disease-free survival (Supplementary Figure 6 CD8+CD103+CD39- p<0.001; CD8+CD103-CD39+ p=0.017; CD8+CD103-CD39- p=0.006). In the stroma only high CD8 T-cells, activated TRM and high CD8+CD103+CD39- T-cells were significantly associated with disease free survival (Supplementary Figure 7). Furthermore, all CD8 T-cell subsets were significantly correlated (p<0.05) with TNM stage, metastases, vascular invasion and microsatellite instability. High infiltration of activated TRM was significantly associated with primary tumour location (see Supplementary Table 2). However, on their own neither MSI nor primary location predicted survival (Supplementary Figure 8). High numbers of intraepithelial activated TRM was an independent prognostic factor (Figure 3C p=0.05). In agreement with the literature (36), high numbers of intraepithelial CD8+ T-cells was also an independent prognostic factor (Figure 3D p=0.05).
Right-sided CRC patients with high CD8 and high CD103 gene expression have increased expression of genes involved in immune pathways
We utilized the data available in the cancer genome atlas (https://portal.gdc.cancer.gov/) to determine the differentially expressed genes (DEG) associated with high (4th Quartile) and low (1st Quartile) expression of CD8 and/or CD103 and the primary tumour location: right-sided (228 patients) or left-sided (131 patients). All groups were compared with CD8 low CD103 low from the same tumour side. Differential analysis identified 7272 DEG in the right-sided CD8 high CD103 high group. This contrast with only 877 DEG identified in the left-sided CD8 high CD103 high. The right- and left-sided CRC patients with CD8 high CD103 low expression had 260 DEG and 201 DEG respectively (Supplementary Figures 9, 10 and Supplementary Table 3). Gene pathway comparison showed the right-sided CRC CD8 high CD103 high group had predominately more gene pathways than the comparable left-sided group (Figures 4A, B). Most of the pathways were associated with immunity. Pathways that are present in right but not left-sided CRC include: protein processing in ER; autophagy; ubiquitin mediated proteolysis; chemokine signaling pathways; cytokine-cytokine interaction and micro RNAs in cancer. There were limited pathways (Supplementary Figure 10) identified in the other groups and this included the high CD8 low CD103 expression group.
Figure 4 Overexpression of pathways predominately in right-sided CRC patients who have high CD8 and CD103 expression: (A) Analysis of the RNAseq data from the TCGA database of right- and left-sided CRC patients: 7245 differentially expressed genes (DEG) for right-sided CD8 high (4th quartile) CD103 high (4th quartile) patients compared with 877 DEG from left-sided CD8 high CD103 high CRC patients (p<0.05 and a logfold at least =1). There were predominately more significant enriched KEGG pathways on the right (shown in red) as opposed to left (where none were significant). The grey nodes represent overlapping pathways. The network of pathways was created using Cytoscape and ClueGo (B) Statistically enriched KEGG pathways with the larger the nodes representing more genes and the darker the red the more significant the pathway.
Only activated TRM T-cells predict survival in left-sided colon cancer
The TCGA analysis indicated the prevalence of immune pathways and DEG associated with high CD8 and high CD103 infiltration in right- versus left-sided CRC. We therefore analysed survival in our patient cohort in terms of tumour location. Although, high numbers of either total CD8+ T-cells or activated TRM predicted survival this was not maintained in left-sided tumours. Indeed, although high infiltration of all immune phenotypes predicted survival in right-sided tumours, activated TRM was the only phenotype that significantly predicted survival in left-sided disease (p=0.008) (Figure 5).
Figure 5 High density of activated TRM T-cells predicts better survival in Left-sided colon cancer patients. Kaplan-Meier plots represent the probability of disease-specific survival for total CD8+ T-cells (A, B), Activated TRM (CD8+CD103+CD39+) T-cells (C, D), CD8+CD103+CD39- T-cells (E, F), CD8+CD103-CD39+ T cells (G, H) and CD8+CD103-CD39- T-cells (I, J) in the tumour epithelium of right-sided colon cancer and Left-sided colorectal cancer patients. Cut-off points used to stratify CRC patients into high and low-density groups were determined using X-tile. The cut-off points for total CD8+ T cells, activated TRM, CD8+CD103+CD39- T-cells, CD8+CD103-CD39+ T-cells and CD8+CD103-CD39- T-cells were 222, 74, 21, 19 and 7 cell/mm2, respectively. The log-rank test was used to compare curves, and p values <0.05 were considered statistically significant.
Disparity between CD8 “immune-hot tumours” and TRM T-cell infiltration
To assess if all “immune-hot” tumours contained elevated numbers of activated TRM, heat maps were employed for all phenotypes. While the majority of patients with increased activated TRM had corresponding elevations of total CD8 T-cells (Supplementary Figure 11), not all patients with high total CD8 T-cells had high levels of activated TRM. Neither the tumour location or MSI status was associated with any of the groups (Supplementary Figure 9).
The absence of activated TRM infiltration associates with poor prognosis irrespective of total CD8 infiltration
Survival was studied in patients with high (“immune-hot”) or low (“immune-cold”) total CD8 T-cells and high or low numbers of activated TRM in the tumour epithelium resulting in 4 groups of patients. Group-1 patients (high CD8 T-cell infiltration and high activated TRM) had the best survival overall Figures 6A(i), B(i) p=0.003). The major CD8+ subset in the tumours were activated TRM. Furthermore, there was greater infiltration of these cells in the epithelial than stroma (p<0.0001) Figure 6B(i): such tumours could be conventionally classified as “immune-hot”.
Figure 6 Tumours with high total CD8 T-cells and low activated TRM cells infiltration predict disease prognosis differently between RCRC and LCRC. (A) CRC patients were classified into four groups depending on total CD8 TILs, and activated TRM infiltration into the TME followed by the assessment of their correlation to survival using Kaplan Meier analysis for the entire (i) CRC cohort, (ii)Right sided and left sided colon cancer patients. (B) The violin plots represent the average densities of the respective CD8+ TIL subsets in the TME for each of the four groups using the Wilcoxon Test: Group 1((i) and (ii); Group 2 ((iii) and (iv)); Group 3 ((v) and (vi)); Group 4 ((vii) and (viii), The left column depicts intraepithelial CD8+ TIL ((i), (iii), (v) and (vii)) and the right column depicts stromal CD8+ TIL ((ii), (iv), (vii) and (viii) (TP: triple positive (CD8+CD103+CD39+), DP 103+ : double positive CD103(CD8+CD103+CD39-), DP 39+: double positive CD39+ (CD8+CD103-CD39+), SP: single positive (CD8+CD103-CD39+)) (C) The pie charts illustrate the numerical proportion of the four groups within RCRC (right-sided colon cancer patient) (n=415) and LCRC (left sided colon cancer patients) (n=326). (D) Comparison of average densities of CD39+ non-lymphoid cells in the stroma of group 1, group 2, group 3 and group 4 tumours. One way Anova/mixed effect analysis with Turkey correction for multiple comparison was used for statistical comparisons in B, while a Kruskal-wallis test with Dunn’s multiple comparison test was used in (D) p-values <0.05 were considered statistically significant. (ns if p > 0.05, *p ≤0.05, **p ≤ 0.01, ***p ≤ 0.001, ****p ≤ 0.0001).
Group-2 tumours lacked pronounced CD8+ T-cell infiltration but nevertheless exhibited high numbers of activated TRM Figures 6A(i), B(ii). These tumours displayed significantly greater CD8 infiltration than Group-4 (low total CD8 low activated TRM) (p<0.0001) Figure 6B(iv). Survival in Group-2 patients was not significantly different to Group-1 indicating that the presence of activated TRM (as opposed to total CD8) are important for survival.
Group-3 had high total CD8 T-cell infiltration but low numbers of activated TRM. Most cells were single positive (CD8+) with few CD8+CD103+ or CD8+CD39+ cells Figure 6B(iii). They had a similar pattern of infiltration in the intraepithelial and stroma. However, despite having high total CD8+ T cells, patients in this group had poor survival. Lastly, group-4 displayed both low total CD8+ infiltration and low numbers of activated TRM with corresponding poor survival Figure 6B. This group were conventional “immune-cold” tumours Figure 6B(iv). Therefore, only patients with high numbers of CD8 T-cells and/or high numbers of activated TRM had improved survival (truly hot tumours). Importantly, patients whose tumours exhibited high total CD8 T-cells in the absence of activated TRM had poor survival Figure 6A(i).
Activated TRM are required for better survival prognosis in left-sided but not right-sided colon cancer
Sub-classification of patients into right and left-sided disease revealed differential survival Figures 6A(ii), (iii). For patients with left-sided tumours survival depended on infiltration by activated TRM. In contrast, for right-sided disease a group of patients exhibited good survival with high total CD8+ T-cells but low activated TRM infiltration. There was no difference in overall number of infiltrating cells between right- and left-sided tumours in any group (data not shown). However, survival was poor in patients with left-sided disease (Figure 6C, Group 3).
The presence of CD39+ on non-lymphocyte stromal cells was associated with poor survival
CD39 and CD73 ectonucleotidases operate in concert to produce immunosuppressive metabolite adenosine in tumour environments (37). We hypothesized that high levels of CD39+ cells in the non-lymphoid stroma cells (e.g., fibroblast) associated with poor survival. As shown in Figure 6D higher numbers of CD39+ cells were observed in group-3 and group-4 patients: those that lacked abundant activated TRM.
Infiltration of total CD8 and activated TRM T-cells consistent across the cores
Colorectal tumours have been shown to exhibit higher levels of immune infiltrate at the invasive margins (38). However, in this study infiltration of both total CD8 and activated TRM was consistent throughout the tumours (Figure 7). Group-1,” immune-hot tumours”, had high numbers of CD8+ T-cells throughout the tumour with lower numbers in adjacent normal tissue. These tumours also had corresponding elevated levels of activated TRM even in the adjacent pathologically normal tissue. Group-3 tumours were characterized by high levels of CD8 T-cells in all cores with no significant difference in any location and limited activated TRM irrespective of region. In contrast, Group 4 “immune-cold” had scant CD8 T-cells throughout, but in adjacent normal tissue these cells were as abundant as those observed in Group-1 and Group-3. This group also had limited infiltration with activated TRM.
Figure 7 Consistent Infiltration of recently activated TRM cells across different tumour regions (LS, CT, IM and ADJ) within the four groups of CRC cohort. (A) The violin plots represent the median density of total CD8 TILs in the respective tumour region matched to the density of total CD8 TILs infiltrating the same tumour region in each patient mapped one below each other in the heat maps (B). Group 1 and group 3 tumours had no significant difference in total CD8 TILs infiltration between LS, CT and IM. The normal adjacent tissue had significantly higher total CD8 TILs infiltration compared to other tumour regions in group 4 tumours (Cold tumours). (C) The violin plots represent the median density of TP TRM cells in the respective tumour region matched to the density of TILs infiltrating the same tumour region in each patient in the heat maps (D). Group 1 tumours (hot tumours) had significantly higher TRM cell infiltration in the LS, CT and IM as compared to ADJ. Meanwhile, group 3 and group 4 tumours had significantly higher activated TRM cell infiltration in the AJD as compared to other regions of the tumour. Mixed effect analysis with Turkeys multiple comparison correction test was used to test for significance. P values <0.05 were considered statistically significant. (*p ≤0.05, **p ≤ 0.01,***p ≤ 0.001, ****p ≤ 0.0001) (LS, luminal side; CT, center of tumour; IM, invasive margin and ADJ, adjacent normal).
Discussion
Our study of TRM cells within the tumour environment is the largest and most detailed study of its kind, having determined the in-situ localization of CD8 T-cell subsets in tissue from nearly 900 CRC patients. We addressed whether the presence of all CD8 T-cells within the tumour are beneficial for survival or whether activated TRM are a better predictor of survival. Using multiplex IHC we established the prognostic significance of activated TRM in CRC. The presence of activated TRM independently predicted survival. Furthermore, the presence of high numbers of activated TRM was a better overall predictor of survival than total CD8+ T-cells. Tumours with high numbers of activated TRM were predominately “immune-hot” (tumours with high total CD8+ T-cell counts), and were highly infiltrated throughout (luminal side, center of the tumour and invasive margins) with both total CD8+ T-cells and activated TRM. While CD8+ T-cells alone did not predict survival in left-sided tumours, activated TRM were uniquely able to predict survival. This contrasted with right-sided tumours where patients exhibiting high total CD8+ T-cells (even with low numbers of activated TRM) had good overall survival.
It has previously been shown that CD8+CD103+CD39+ T-cells identifies tumour reactive T-cells in both head and neck cancer (21) and breast cancer (39), and that their presence associates with survival. While those studies showed the functionality of the activated TRM, we focused on their survival advantage and location in tumours. Our findings that activated TRM are an independent predictor would support the concept that these T-cells have been active in the tumour environment.
The immunoscore uses the presence of CD3 and CD8 T-cell in the invasive margins and center of the tumour to predict survival in CRC (40–42). As far as we are aware we are the first to show that CD8 alone does not predict survival in left-sided colon cancer and that TRM uniquely predict survival. Other studies have shown an increase in CD8 gene expression in right-sided CRC and suggested that these patients have a better immune response (18). Indeed, TCGA analysis of right-sided CRC patients with high CD8 and high CD103 compared to the low CD8 low CD103 had significantly more immune pathway involvement than patients with left-sided tumours, suggesting greater immune responses in right-sided disease. This is the first time both high CD8 and high CD103 expression have been analysed together using data from the TCGA database for CRC and the results are consistent with our data showing that there are differences in immune responses between right and left-sided CRC.
The presence of high CD8 T-cells irrespective of subtype predicts survival in right-sided CRC. Intriguingly, “immune-hot” tumours with low numbers of activated TRM had good survival, suggesting that although activated TRM are important as predictors of survival there may be other markers of antigen specificity and activation that need to be defined for these patients. Why infiltration of high total CD8+ T-cells without TRM is a good prognosis on right-sided but not left-sided colon cancer is unclear, however, other leukocytes within the tumours may be contributing factors. These could include T-regs and tumour associated macrophages. In a study in breast cancer, the presence of T-regs was responsible for the loss of prognostic significance of high numbers of TRM (39).
The difference in survival relating to T-cell infiltration in the right versus the left CRC could be due to differences in somatic mutations and neoantigen generation. The mutations rates in BRAF, POLE, POLD1 and PIK3CA genes associate with right-sided CRC. Patients with these mutations (43) have elevated expression of helper T-cells, class II-related genes, chemokines and inhibitory molecules. In contrast, RAS mutant (KRAS, NRAS) tumours (which are more frequent on the left) associate with poor immune infiltration, low inhibitory molecule expression (43) and recruitment of suppressive myeloid cells. Nevertheless, this does not explain why CRC patients with tumours heavily infiltrated with CD8 T-cells but lacking activated TRM have good prognosis if their tumour is right-sided. It is not that the TRM have not upregulated CD39 in these tumours, but that there are low levels of TRM: most cells are single positive CD8 T-cells. The literature would support the theory that CD8+CD103+CD39- T-cells are bystander without specificity for tumour antigens (20) while activated TRM are tumour specific. Left-sided CRC generally exhibit lower mutational burden, and therefore may present less antigenic opportunity for antigen specific TRM. The association between the absence of activated TRM with high levels of CD39+ stromal cells (e.g. fibroblasts) may reflect the role of CD39 in ATP depletion and suppression of T-cell activation and may consequently account for the lack of TRM in these patients. Quantifying activated TRM in CRC patients may therefore permit stratification of tumours with strong immunity genuinely focused on the tumour.
Bindea et al. (38) in a study of 107 patients, showed elevated total CD8+ cells in the invasive margins than the rest of the tumour. In contrast, we observed in 9-fold more patients that high total CD8+ T-cells were homogeneously distributed throughout the “immune-hot” tumours. Whereas in “immune-cold” tumours there were higher numbers of CD8+ T-cells present in pathological adjacent normal tissue than in the tumours. The difference could be due to their analysis of the whole cohort of CRC patients together compared to our analysis of patients with high and low total CD8 T-cell separately.
CD39 and CD103 upregulation on tumour infiltrating CD8 T-cells is due to TGFβ and chronic TCR stimulation (21). Once a cell is activated it displays markers of exhaustion. The expression of CD39 and TIM3 have been shown to discriminate exhausted cells from memory/effector T-cells (44). Expression of other markers of exhaustion such as PD-1 associate with terminal exhaustion (45). It is of note that in our study the “immune-hot” tumours had high levels of all CD8 phenotypes including both CD8+CD103+ and CD8+CD39+. This could suggest an active immune response in the tumour environment of these patients and therefore the reason activated TRM herald survival is that they have previously responded to the tumour. It can be argued that patients who have previously mounted an immune response could continue to mount new T-cell responses.
Clinical use of immune checkpoint inhibitors has revolutionized the treatment of solid tumours. However, future work is needed to identify CRC patients suitable for immunotherapy and characterize repertoires of T-cells to be targeted, a point illustrated by studies showing the presence of intra-tumoural CD103+CD8+ T cells predicts responses to PDL-1 blockade (46, 47). We have already argued that activated TRM have markers of exhaustion, and that the presence of TRM in the tumour is indicative of active immunity. While there is evidence that exhausted T-cells can respond to checkpoint blockade (48, 49), there is also evidence that these cells are not readily reactivated (50). Indeed, recent studies show that new responses to checkpoint inhibitors involve CD8 T-cells actively replenished from outside the tumour (50) and not from the pre-existing repertoire of exhausted T-cells (51). Patients who have previously mounted an immune response could continue to mount new responses with appropriate immunotherapy. In this paper we have shown over 90% survival in the group of patients with high numbers of total CD8 and high numbers of activated TRM. Therefore, patients who may benefit from future immunotherapies may not be the ones with activated TRM who already have generated a good immune prognosis but those patients who have shown the ability to recruit CD8 T-cells but have so far failed to convert this to a survival advantage: the left-sided colon cancer patients with high CD8 T-cells and low recently activated TRM. The challenge will be to generate activated tumour specific TRM in these patients.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by Oxford Rec B, NHS (reference 05/Q1605/66). The patients/participants provided their written informed consent to participate in this study.
Author contributions
Data acquisition: ST, NM, KH, SM, AN, AL. Analysis and interpretation of data: JR, ST, NM, LF, NK, IS, AJ, MI, AM-G, TD. Building the patient cohort: WF, MI, Conception/design of work: JR, AJ, ST. Oversight of entire project: JR. Manuscript preparation and review: JR, ST, AJ, MI, IS, BW. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the Medical Research Council, UK grant to MI (grant code: MR/N005953/1). Multiplex imaging was supported by a Wellcome Trust Multi-User Equipment Grant awarded to BW. This research was funded in whole, or in part, by the Wellcome Trust (Grant number 218287/Z/19/Z). For the purpose of open access, the author has applied a CC BY public copyright license to any Author Accepted Manuscript version arising from this submission.
Acknowledgments
The authors would like to acknowledge Akoya Biosciences and particularly Roslyn Lloyd for advice and training with the inform software and Matt's Fund (in memory of Matthew Boulton, raising awareness of Malignant Melanoma and funding related research) for their donation to support the computer analysis.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2023.1057292/full#supplementary-material
References
1. Flemer B, Lynch DB, Brown JM, Jeffery IB, Ryan FJ, Claesson MJ, et al. Tumour-associated and non-tumour-associated microbiota in colorectal cancer. Gut (2017) 66:633–43. doi: 10.1136/gutjnl-2015-309595
2. Mukund K, Syulyukina N, Ramamoorthy S, Subramaniam S. Right and left-sided colon cancers - specificity of molecular mechanisms in tumorigenesis and progression. BMC Cancer (2020) 20:317. doi: 10.1186/s12885-020-06784-7
3. Narayanan S, Gabriel E, Attwood K, Boland P, Nurkin S. Association of clinicopathologic and molecular markers on stage-specific survival of right versus left colon cancer. Clin Colorectal Cancer (2018) 17:e671–e8. doi: 10.1016/j.clcc.2018.07.001
4. Arnold D, Lueza B, Douillard JY, Peeters M, Lenz HJ, Venook A, et al. Prognostic and predictive value of primary tumour side in patients with RAS wild-type metastatic colorectal cancer treated with chemotherapy and EGFR directed antibodies in six randomized trials. Ann Oncol (2017) 28:1713–29. doi: 10.1093/annonc/mdx175
5. Baran B, Mert Ozupek N, Yerli Tetik N, Acar E, Bekcioglu O, Baskin Y. Difference between left-sided and right-sided colorectal cancer: a focused review of literature. Gastroenterol Res (2018) 11:264–73. doi: 10.14740/gr1062w
6. Salem ME, Weinberg BA, Xiu J, El-Deiry WS, Hwang JJ, Gatalica Z, et al. Comparative molecular analyses of left-sided colon, right-sided colon, and rectal cancers. Oncotarget (2017) 8:86356–68. doi: 10.18632/oncotarget.21169
7. Cancer Genome Atlas N. Comprehensive molecular characterization of human colon and rectal cancer. Nature (2012) 487:330–7. doi: 10.1038/nature11252
8. McAllister SS, Weinberg RA. The tumour-induced systemic environment as a critical regulator of cancer progression and metastasis. Nat Cell Biol (2014) 16:717–27. doi: 10.1038/ncb3015
9. Galon J, Costes A, Sanchez-Cabo F, Kirilovsky A, Mlecnik B, Lagorce-Pages C, et al. Type, density, and location of immune cells within human colorectal tumors predict clinical outcome. Science (2006) 313:1960–4. doi: 10.1126/science.1129139
10. Simpson JA, Al-Attar A, Watson NF, Scholefield JH, Ilyas M, Durrant LG, et al. Intratumoral T cell infiltration and STAT1 as biomarkers of good prognosis in colorectal cancer. Gut (2010) 59:926–33. doi: 10.1136/gut.2009.194472
11. Herbst RS, Soria JC, Kowanetz M, Fine GD, Hamid O, Gordon MS, et al. Predictive correlates of response to the anti-PD-L1 antibody MPDL3280A in cancer patients. Nature (2014) 515:563–7. doi: 10.1038/nature14011
12. Tumeh PC, Harview CL, Yearley JH, Shintaku IP, Taylor EJ, Robert L, et al. PD-1 blockade induces responses by inhibiting adaptive immune resistance. Nature (2014) 515:568–71. doi: 10.1038/nature13954
13. Alexandrov LB, Nik-Zainal S, Wedge DC, Aparicio SA, Behjati S, Biankin AV, et al. Signatures of mutational processes in human cancer. Nature (2013) 500:415–21. doi: 10.1038/nature12477
14. Le DT, Durham JN, Smith KN, Wang H, Bartlett BR, Aulakh LK, et al. Mismatch repair deficiency predicts response of solid tumors to PD-1 blockade. Science (2017) 357:409–13. doi: 10.1126/science.aan6733
15. Van Allen EM, Miao D, Schilling B, Shukla SA, Blank C, Zimmer L, et al. Genomic correlates of response to CTLA-4 blockade in metastatic melanoma. Science (2015) 350:207–11. doi: 10.1126/science.aad0095
16. Maby P, Tougeron D, Hamieh M, Mlecnik B, Kora H, Bindea G, et al. Correlation between density of CD8+ T-cell infiltrate in microsatellite unstable colorectal cancers and frameshift mutations: a rationale for personalized immunotherapy. Cancer Res (2015) 75:3446–55. doi: 10.1158/0008-5472.CAN-14-3051
17. Mlecnik B, Bindea G, Angell HK, Maby P, Angelova M, Tougeron D, et al. Integrative analyses of colorectal cancer show immunoscore is a stronger predictor of patient survival than microsatellite instability. Immunity (2016) 44:698–711. doi: 10.1016/j.immuni.2016.02.025
18. Zhang L, Zhao Y, Dai Y, Cheng JN, Gong Z, Feng Y, et al. Immune landscape of colorectal cancer tumor microenvironment from different primary tumor location. Front Immunol (2018) 9:1578. doi: 10.3389/fimmu.2018.01578
19. Galon J, Bruni D. Approaches to treat immune hot, altered and cold tumours with combination immunotherapies. Nat Rev Drug Discovery (2019) 18:197–218. doi: 10.1038/s41573-018-0007-y
20. Simoni Y, Becht E, Fehlings M, Loh CY, Koo SL, Teng KWW, et al. Bystander CD8(+) T cells are abundant and phenotypically distinct in human tumour infiltrates. Nature (2018) 557:575–9. doi: 10.1038/s41586-018-0130-2
21. Duhen T, Duhen R, Montler R, Moses J, Moudgil T, de Miranda NF, et al. Co-Expression of CD39 and CD103 identifies tumor-reactive CD8 T cells in human solid tumors. Nat Commun (2018) 9:2724. doi: 10.1038/s41467-018-05072-0
22. Raczkowski F, Rissiek A, Ricklefs I, Heiss K, Schumacher V, Wundenberg K, et al. CD39 is upregulated during activation of mouse and human T cells and attenuates the immune response to listeria monocytogenes. PloS One (2018) 13:e0197151. doi: 10.1371/journal.pone.0197151
23. Ganesan AP, Clarke J, Wood O, Garrido-Martin EM, Chee SJ, Mellows T, et al. Tissue-resident memory features are linked to the magnitude of cytotoxic T cell responses in human lung cancer. Nat Immunol (2017) 18:940–50. doi: 10.1038/ni.3775
24. Schenkel JM, Masopust D. Tissue-resident memory T cells. Immunity (2014) 41:886–97. doi: 10.1016/j.immuni.2014.12.007
25. Okla K, Farber DL, Zou W. Tissue-resident memory T cells in tumor immunity and immunotherapy. J Exp Med (2021) 218(4): e20201605. doi: 10.1084/jem.20201605
26. Wang ZQ, Milne K, Derocher H, Webb JR, Nelson BH, Watson PH. CD103 and intratumoral immune response in breast cancer. Clin Cancer Res (2016) 22:6290–7. doi: 10.1158/1078-0432.CCR-16-0732
27. Webb JR, Milne K, Watson P, Deleeuw RJ, Nelson BH. Tumor-infiltrating lymphocytes expressing the tissue resident memory marker CD103 are associated with increased survival in high-grade serous ovarian cancer. Clin Cancer Res (2014) 20:434–44. doi: 10.1158/1078-0432.CCR-13-1877
28. Djenidi F, Adam J, Goubar A, Durgeau A, Meurice G, de Montpreville V, et al. CD8+CD103+ tumor-infiltrating lymphocytes are tumor-specific tissue-resident memory T cells and a prognostic factor for survival in lung cancer patients. J Immunol (2015) 194:3475–86. doi: 10.4049/jimmunol.1402711
29. van den Bulk J, Verdegaal EME, Ruano D, Ijsselsteijn ME, Visser M, van der Breggen R, et al. Neoantigen-specific immunity in low mutation burden colorectal cancers of the consensus molecular subtype 4. Genome Med (2019) 11:87. doi: 10.1186/s13073-019-0697-8
30. Popple A, Durrant LG, Spendlove I, Rolland P, Scott IV, Deen S, et al. The chemokine, CXCL12, is an independent predictor of poor survival in ovarian cancer. Br J Cancer (2012) 106:1306–13. doi: 10.1038/bjc.2012.49
31. Mori H, Bolen J, Schuetter L, Massion P, Hoyt CC, VandenBerg S, et al. Characterizing the tumor immune microenvironment with tyramide-based multiplex immunofluorescence. J Mammary Gland Biol Neoplasia (2020) 25:417–32. doi: 10.1007/s10911-021-09479-2
32. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol (2014) 15:550. doi: 10.1186/s13059-014-0550-8
33. Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, et al. ClueGO: a cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics (2009) 25:1091–3. doi: 10.1093/bioinformatics/btp101
34. Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res (2003) 13:2498–504. doi: 10.1101/gr.1239303
35. Liao Y, Wang J, Jaehnig EJ, Shi Z, Zhang B. WebGestalt 2019: gene set analysis toolkit with revamped UIs and APIs. Nucleic Acids Res (2019) 47:W199–205. doi: 10.1093/nar/gkz401
36. Bruni D, Angell HK, Galon J. The immune contexture and immunoscore in cancer prognosis and therapeutic efficacy. Nat Rev Cancer (2020) 20:662–80. doi: 10.1038/s41568-020-0285-7
37. Ghiringhelli F, Bruchard M, Chalmin F, Rebe C. Production of adenosine by ectonucleotidases: a key factor in tumor immunoescape. J BioMed Biotechnol (2012) 2012:473712. doi: 10.1155/2012/473712
38. Bindea G, Mlecnik B, Tosolini M, Kirilovsky A, Waldner M, Obenauf AC, et al. Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer. Immunity (2013) 39:782–95. doi: 10.1016/j.immuni.2013.10.003
39. Losurdo A, Scirgolea C, Alvisi G, Brummelman J, Errico V, Di Tommaso L, et al. Single-cell profiling defines the prognostic benefit of CD39(high) tissue resident memory CD8+ T cells in luminal-like breast cancer. Commun Biol (2021) 4:1117. doi: 10.1038/s42003-021-02595-z
40. Pages F, Mlecnik B, Marliot F, Bindea G, Ou FS, Bifulco C, et al. International validation of the consensus immunoscore for the classification of colon cancer: a prognostic and accuracy study. Lancet (2018) 391:2128–39. doi: 10.1016/S0140-6736(18)30789-X
41. Galon J, Pages F, Marincola FM, Angell HK, Thurin M, Lugli A, et al. Cancer classification using the immunoscore: a worldwide task force. J Transl Med (2012) 10:205. doi: 10.1186/1479-5876-10-205
42. Angell HK, Bruni D, Barrett JC, Herbst R, Galon J. The immunoscore: colon cancer and beyond. Clin Cancer Res (2020) 26:332–9. doi: 10.1158/1078-0432.CCR-18-1851
43. Lal N, Beggs AD, Willcox BE, Middleton GW. An immunogenomic stratification of colorectal cancer: implications for development of targeted immunotherapy. Oncoimmunology (2015) 4:e976052. doi: 10.4161/2162402X.2014.976052
44. Sade-Feldman M, Yizhak K, Bjorgaard SL, Ray JP, de Boer CG, Jenkins RW, et al. Defining T cell states associated with response to checkpoint immunotherapy in melanoma. Cell (2018) 175:998–1013.e20. doi: 10.1016/j.cell.2018.10.038
45. Gupta PK, Godec J, Wolski D, Adland E, Yates K, Pauken KE, et al. CD39 expression identifies terminally exhausted CD8+ T cells. PloS Pathog (2015) 11:e1005177. doi: 10.1371/journal.ppat.1005177
46. Banchereau R, Chitre AS, Scherl A, Wu TD, Patil NS, de Almeida P, et al. Intratumoral CD103+ CD8+ T cells predict response to PD-L1 blockade. J Immunother Cancer (2021) 9:e002231. doi: 10.1136/jitc-2020-002231
47. Hsu CL, Ou DL, Bai LY, Chen CW, Lin L, Huang SF, et al. Exploring markers of exhausted CD8 T cells to predict response to immune checkpoint inhibitor therapy for hepatocellular carcinoma. Liver Cancer (2021) 10:346–59. doi: 10.1159/000515305
48. Nelson CE, Mills LJ, McCurtain JL, Thompson EA, Seelig DM, Bhela S, et al. Reprogramming responsiveness to checkpoint blockade in dysfunctional CD8 T cells. Proc Natl Acad Sci U.S.A. (2019) 116:2640–5. doi: 10.1073/pnas.1810326116
49. Kallies A, Zehn D, Utzschneider DT. Precursor exhausted T cells: key to successful immunotherapy? Nat Rev Immunol (2020) 20:128–36. doi: 10.1038/s41577-019-0223-7
50. Yost KE, Satpathy AT, Wells DK, Qi Y, Wang C, Kageyama R, et al. Clonal replacement of tumor-specific T cells following PD-1 blockade. Nat Med (2019) 25:1251–9. doi: 10.1038/s41591-019-0522-3
Keywords: colorectal cancer, T-cells, multiplex IHC/IF, tissue resident T cells, immune microenvironment, CD8 T-cells, cancer
Citation: Talhouni S, Fadhil W, Mongan NP, Field L, Hunter K, Makhsous S, Maciel-Guerra A, Kaur N, Nestarenkaite A, Laurinavicius A, Willcox BE, Dottorini T, Spendlove I, Jackson AM, Ilyas M and Ramage JM (2023) Activated tissue resident memory T-cells (CD8+CD103+CD39+) uniquely predict survival in left sided “immune-hot” colorectal cancers. Front. Immunol. 14:1057292. doi: 10.3389/fimmu.2023.1057292
Received: 29 September 2022; Accepted: 21 April 2023;
Published: 11 May 2023.
Edited by:
Catherine Sautes-Fridman, INSERM U1138 Centre de Recherche des Cordeliers (CRC), FranceReviewed by:
Datian Chen, Haimen City People’s Hospital, ChinaKevin D. Pavelko, Mayo Clinic, United States
Copyright © 2023 Talhouni, Fadhil, Mongan, Field, Hunter, Makhsous, Maciel-Guerra, Kaur, Nestarenkaite, Laurinavicius, Willcox, Dottorini, Spendlove, Jackson, Ilyas and Ramage. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Judith M. Ramage, SnVkaXRoLnJhbWFnZUBub3R0aW5naGFtLmFjLnVr