Integrative genomics identifies SHPRH as a tumor suppressor gene in lung adenocarcinoma that regulates DNA damage response.

BACKGROUND: Identification of driver mutations and development of targeted therapies has considerably improved outcomes for lung cancer patients. However, significant limitations remain with the lack of identified drivers in a large subset of patients. Here, we aimed to assess the genomic landscape of lung adenocarcinomas (LUADs) from individuals without a history of tobacco use to reveal new genetic drivers of lung cancer. METHODS: Integrative genomic analyses combining whole-exome sequencing, copy number, and mutational information for 83 LUAD tumors was performed and validated using external datasets to identify genetic variants with a predicted functional consequence and assess association with clinical outcomes. LUAD cell lines with alteration of identified candidates were used to functionally characterize tumor suppressive potential using a conditional expression system both in vitro and in vivo. RESULTS: We identified 21 genes with evidence of positive selection, including 12 novel candidates that have yet to be characterized in LUAD. In particular, SNF2 Histone Linker PHD RING Helicase (SHPRH) was identified due to its frequency of biallelic disruption and location within the familial susceptibility locus on chromosome arm 6q. We found that low SHPRH mRNA expression is associated with poor survival outcomes in LUAD patients. Furthermore, we showed that re-expression of SHPRH in LUAD cell lines with inactivating alterations for SHPRH reduces their in vitro colony formation and tumor burden in vivo. Finally, we explored the biological pathways associated SHPRH inactivation and found an association with the tolerance of LUAD cells to DNA damage. CONCLUSIONS: These data suggest that SHPRH is a tumor suppressor gene in LUAD, whereby its expression is associated with more favorable patient outcomes, reduced tumor and mutational burden, and may serve as a predictor of response to DNA damage. Thus, further exploration into the role of SHPRH in LUAD development may make it a valuable biomarker for predicting LUAD risk and prognosis.

An improved algorithm using a Health Canada-approved DNA-image cytometry system for non-invasive screening of high-grade oral lesions.

BACKGROUND: DNA-image cytometry (DNA-ICM) is able to detect gross alterations of cellular DNA-content representing aneuploidy, a biomarker of malignancy. A Health Canada-approved DNA-ICM system, ClearCyte® in combination with a cytopathologist's review, has demonstrated high sensitivity (89%) and specificity (97%) in identifying high-grade oral lesions. The study objective was to create an improved automated algorithm (iClearcyte) and test its robustness in differentiating high grade from benign reactive oral lesions without a cytopathologist's input. METHODS: A set of 214 oral brushing samples of oral cancer (n = 92), severe dysplasia (n = 20), reactive lesions (n = 52), and normal samples (n = 50) were spun down onto slides and stained using Feulgen-Thionin reaction. Following ClearCyte® scan, nuclear features were calculated, and nuclei categorized into "diploid," "hyperdiploid," "tetraploid," and "aneuploid" DNA ploidy groups by the ClearCyte® software. The samples were randomized into training and test sets (70:30) based on patient's age, sex, tobacco use, and lesion site risk. The training set was used to create a new algorithm which was then validated using the remaining samples in the test set, where sensitivity, specificity, positive predictive value (PPV), and negative predictive value (NPV) were calculated. RESULTS: The proposed iClearCyte algorithm (>1 "aneuploid" cell or ≥ 1.7% combined "hyperdiploid" and "tetraploid" nuclei frequency) identified high-grade samples with sensitivity, specificity, PPV, and NPV of 100.0%, 86.7%, 89.7%, and 100.0%, respectively, in the test set. CONCLUSION: The iClearCyte test has potential to serve as a robust non-invasive automated oral cancer screening tool promoting early oral cancer detection and decreasing the number of unnecessary invasive biopsies.

Dual-beam manually actuated distortion-corrected imaging (DMDI): two dimensional scanning with a single-axis galvanometer.

We recently demonstrated a new two-dimensional imaging paradigm called dual-beam manually actuated distortion-corrected imaging (DMDI). This technique uses a single mechanical scanner and two spatially separated beams to determine relative sample velocity and simultaneously corrects image distortions due to manual actuation. DMDI was first demonstrated using a rotating dual-beam micromotor catheter. Here, we present a new implementation of DMDI using a single axis galvanometer to scan a pair of beams in approximately parallel lines onto a sample. Furthermore, we present a method for automated distortion correction based on frame co-registration between images acquired by the two beams. Distortion correction is possible for manually actuated motion both perpendicular and parallel to the galvanometer-scanned lines. Using en face OCT as the imaging modality, we demonstrate DMDI and the automated distortion correction algorithm for imaging a printed paper phantom, a dragon fruit, and a fingerprint.

Quantification of large scale DNA organization for predicting prostate cancer recurrence.

This study investigates whether Genomic Organization at Large Scales (which we propose to call GOALS) as quantified via nuclear phenotype characteristics and cell sociology features (describing cell organization within tissue) collected from prostate tissue microarrays (TMAs) can separate biochemical failure from biochemical nonevidence of disease (BNED) after radical prostatectomy (RP). Of the 78 prostate cancer tissue cores collected from patients treated with RP, 16 who developed biochemical relapse (failure group) and 16 who were BNED patients (nonfailure group) were included in the analyses (36 cores from 32 patients). A section from this TMA was stained stoichiometrically for DNA using the Feulgen-Thionin methodology, and scanned with a Pannoramic MIDI scanner. Approximately 110 nuclear phenotypic features, predominately quantifying large scale DNA organization (GOALS), were extracted from each segmented nuclei. In addition, the centers of these segmented nuclei defined a Voronoi tessellation and subsequent architectural analysis. Prostate TMA core classification as biochemical failure or BNED after RP using GOALS features was conducted (a) based on cell type and cell position within the epithelium (all cells, all epithelial cells, epithelial >2 cell layers away from basement membrane) from all cores, and (b) based on epithelial cells more than two cell layers from the basement membrane using a Classifier trained on Gleason 6, 8, 9 (16 cores) only and applied to a Test set consisting of the Gleason 7 cores (20 cores). Successful core classification as biochemical failure or BNED after RP by a linear classifier was 75% using all cells, 83% using all epithelial cells, and 86% using epithelial >2 layers. Overall success of predicted classification by the linear Classifier of (b) was 87.5% using the Training Set and 80% using the Test Set. Overall success of predicted progression using Gleason score alone was 75% for Gleason >7 as failures and 69% for Gleason >6 as failures. © 2017 International Society for Advancement of Cytometry.

Dual-beam manually-actuated distortion-corrected imaging (DMDI) with micromotor catheters.

We present a new paradigm for performing two-dimensional scanning called dual-beam manually-actuated distortion-corrected imaging (DMDI). DMDI operates by imaging the same object with two spatially-separated beams that are being mechanically scanned rapidly in one dimension with slower manual actuation along a second dimension. Registration of common features between the two imaging channels allows remapping of the images to correct for distortions due to manual actuation. We demonstrate DMDI using a 4.7 mm OD rotationally scanning dual-beam micromotor catheter (DBMC). The DBMC requires a simple, one-time calibration of the beam paths by imaging a patterned phantom. DMDI allows for distortion correction of non-uniform axial speed and rotational motion of the DBMC. We show the utility of this technique by demonstrating en face OCT image distortion correction of a manually-scanned checkerboard phantom and fingerprint scan.

Developmental transcription factor NFIB is a putative target of oncofetal miRNAs and is associated with tumour aggressiveness in lung adenocarcinoma.

Genes involved in fetal lung development are thought to play crucial roles in the malignant transformation of adult lung cells. Consequently, the study of lung tumour biology in the context of lung development has the potential to reveal key developmentally relevant genes that play critical roles in lung cancer initiation/progression. Here, we describe for the first time a comprehensive characterization of miRNA expression in human fetal lung tissue, with subsequent identification of 37 miRNAs in non-small cell lung cancer (NSCLC) that recapitulate their fetal expression patterns. Nuclear factor I/B (NFIB), a transcription factor essential for lung development, was identified as a potential frequent target for these 'oncofetal' miRNAs. Concordantly, analysis of NFIB expression in multiple NSCLC independent cohorts revealed its recurrent underexpression (in ∼40-70% of tumours). Interrogation of NFIB copy number, methylation, and mutation status revealed that DNA level disruption of this gene is rare, and further supports the notion that oncofetal miRNAs are likely the primary mechanism responsible for NFIB underexpression in NSCLC. Reflecting its functional role in regulating lung differentiation, low expression of NFIB was significantly associated with biologically more aggressive subtypes and, ultimately, poorer survival in lung adenocarcinoma patients. Copyright © 2016 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Endoscopic high-resolution autofluorescence imaging and OCT of pulmonary vascular networks.

High-resolution imaging from within airways may allow new methods for studying lung disease. In this work, we report an endoscopic imaging system capable of high-resolution autofluorescence imaging (AFI) and optical coherence tomography (OCT) in peripheral airways using a 0.9 mm diameter double-clad fiber (DCF) catheter. In this system, AFI excitation light is coupled into the core of the DCF, enabling tightly focused excitation light while maintaining efficient collection of autofluorescence emission through the large diameter inner cladding of the DCF. We demonstrate the ability of this imaging system to visualize pulmonary vasculature as small as 12 µm in vivo.

Quantification of confocal fluorescence microscopy for the detection of cervical intraepithelial neoplasia.

BACKGROUND: Cervical cancer remains a major health problem, especially in developing countries. Colposcopic examination is used to detect high-grade lesions in patients with a history of abnormal pap smears. New technologies are needed to improve the sensitivity and specificity of this technique. We propose to test the potential of fluorescence confocal microscopy to identify high-grade lesions. METHODS: We examined the quantification of ex vivo confocal fluorescence microscopy to differentiate among normal cervical tissue, low-grade Cervical Intraepithelial Neoplasia (CIN), and high-grade CIN. We sought to (1) quantify nuclear morphology and tissue architecture features by analyzing images of cervical biopsies; and (2) determine the accuracy of high-grade CIN detection via confocal microscopy relative to the accuracy of detection by colposcopic impression. Forty-six biopsies obtained from colposcopically normal and abnormal cervical sites were evaluated. Confocal images were acquired at different depths from the epithelial surface and histological images were analyzed using in-house software. RESULTS: The features calculated from the confocal images compared well with those features obtained from the histological images and histopathological reviews of the specimens (obtained by a gynecologic pathologist). The correlations between two of these features (the nuclear-cytoplasmic ratio and the average of three nearest Delaunay-neighbors distance) and the grade of dysplasia were higher than that of colposcopic impression. The sensitivity of detecting high-grade dysplasia by analysing images collected at the surface of the epithelium, and at 15 and 30 µm below the epithelial surface were respectively 100, 100, and 92 %. CONCLUSIONS: Quantitative analysis of confocal fluorescence images showed its capacity for discriminating high-grade CIN lesions vs. low-grade CIN lesions and normal tissues, at different depth of imaging. This approach could be used to help clinicians identify high-grade CIN in clinical settings.

DNA methylation is globally disrupted and associated with expression changes in chronic obstructive pulmonary disease small airways.

DNA methylation is an epigenetic modification that is highly disrupted in response to cigarette smoke and involved in a wide spectrum of malignant and nonmalignant diseases, but surprisingly not previously assessed in small airways of patients with chronic obstructive pulmonary disease (COPD). Small airways are the primary sites of airflow obstruction in COPD. We sought to determine whether DNA methylation patterns are disrupted in small airway epithelia of patients with COPD, and evaluate whether changes in gene expression are associated with these disruptions. Genome-wide methylation and gene expression analysis were performed on small airway epithelial DNA and RNA obtained from the same patient during bronchoscopy, using Illumina's Infinium HM27 and Affymetrix's Genechip Human Gene 1.0 ST arrays. To control for known effects of cigarette smoking on DNA methylation, methylation and gene expression profiles were compared between former smokers with and without COPD matched for age, pack-years, and years of smoking cessation. Our results indicate that aberrant DNA methylation is (1) a genome-wide phenomenon in small airways of patients with COPD, and (2) associated with altered expression of genes and pathways important to COPD, such as the NF-E2-related factor 2 oxidative response pathway. DNA methylation is likely an important mechanism contributing to modulation of genes important to COPD pathology. Because these methylation events may underlie disease-specific gene expression changes, their characterization is a critical first step toward the development of epigenetic markers and an opportunity for developing novel epigenetic therapeutic interventions for COPD.

A statistical model for removing inter-device differences in spectroscopy.

We are investigating spectroscopic devices designed to make in vivo cervical tissue measurements to detect pre-cancerous and cancerous lesions. All devices have the same design and ideally should record identical measurements. However, we observed consistent differences among them. An experiment was designed to study the sources of variation in the measurements recorded. Here we present a log additive statistical model that incorporates the sources of variability we identified. Based on this model, we estimated correction factors from the experimental data needed to eliminate the inter-device variability and other sources of variation. These correction factors are intended to improve the accuracy and repeatability of such devices when making future measurements on patient tissue.

Fiber-optic polarization diversity detection for rotary probe optical coherence tomography.

We report a polarization diversity detection scheme for optical coherence tomography with a new, custom, miniaturized fiber coupler with single mode (SM) fiber inputs and polarization maintaining (PM) fiber outputs. The SM fiber inputs obviate matching the optical lengths of the X and Y OCT polarization channels prior to interference and the PM fiber outputs ensure defined X and Y axes after interference. Advantages for this scheme include easier alignment, lower cost, and easier miniaturization compared to designs with free-space bulk optical components. We demonstrate the utility of the detection system to mitigate the effects of rapidly changing polarization states when imaging with rotating fiber optic probes in Intralipid suspension and during in vivo imaging of human airways.

Smoking status impacts microRNA mediated prognosis and lung adenocarcinoma biology.

BACKGROUND: Cigarette smoke is associated with the majority of lung cancers: however, 25% of lung cancer patients are non-smokers, and half of all newly diagnosed lung cancer patients are former smokers. Lung tumors exhibit distinct epidemiological, clinical, pathological, and molecular features depending on smoking status, suggesting divergent mechanisms underlie tumorigenesis in smokers and non-smokers. MicroRNAs (miRNAs) are integral contributors to tumorigenesis and mediate biological responses to smoking. Based on the hypothesis that smoking-specific miRNA differences in lung adenocarcinomas reflect distinct tumorigenic processes selected by different smoking and non-smoking environments, we investigated the contribution of miRNA disruption to lung tumor biology and patient outcome in the context of smoking status. METHODS: We applied a whole transcriptome sequencing based approach to interrogate miRNA levels in 94 patient-matched lung adenocarcinoma and non-malignant lung parenchymal tissue pairs from current, former and never smokers. RESULTS: We discovered novel and distinct smoking status-specific patterns of miRNA and miRNA-mediated gene networks, and identified miRNAs that were prognostically significant in a smoking dependent manner. CONCLUSIONS: We conclude that miRNAs disrupted in a smoking status-dependent manner affect distinct cellular pathways and differentially influence lung cancer patient prognosis in current, former and never smokers. Our findings may represent promising biologically relevant markers for lung cancer prognosis or therapeutic intervention.

Molecular fixative enables expression microarray analysis of microdissected clinical cervical specimens.

Formalin-fixed tissue has been a mainstay of clinical pathology laboratories, but formalin alters many biomolecules, including nucleic acids and proteins. Meanwhile, frozen tissues contain better-preserved biomolecules, but tissue morphology is affected, limiting their diagnostic utility. Molecular fixatives promise to bridge this gap by simultaneously preserving morphology and biomolecules, enabling clinical diagnosis and molecular analyses on the same specimen. While previous reports have broadly evaluated the use of molecular fixative in various human tissues, we present here the first detailed assessment of the applicability of molecular fixative to both routine histopathological diagnosis and molecular analysis of cervical tissues. Ten specimens excised via the loop electrosurgical excision procedure, which removes conical tissue samples from the cervix, were cut into alternating pieces preserved in either formalin or molecular fixative. Cervical specimens preserved in molecular fixative were easily interpretable, despite featuring more eosinophilic cytoplasm and more recognizable chromatin texture than formalin-fixed specimens. Immunohistochemical staining patterns of p16 and Ki-67 were similar between fixatives, although Ki-67 staining was stronger in the molecular fixative specimens. The RNA of molecular fixative specimens from seven cases representing various dysplasia grades was assessed for utility in expression microarray analysis. Cluster analysis and scatter plots of duplicate samples suggest that data of sufficient quality can be obtained from as little as 50ng of RNA from molecular fixative samples. Taken together, our results show that molecular fixative may be a more versatile substitute for formalin, simultaneously preserving tissue morphology for clinical diagnosis and biomolecules for immunohistochemistry and gene expression analysis.

Influence of fluorescence on screening decisions for oral mucosal lesions in community dental practices.

BACKGROUND: Quality of oral screening examinations is dependent upon the experience of the clinician and can vary widely. Deciding when a patient needs to be referred is a critical and difficult decision for general practice clinicians. A device to aid in this decision would be beneficial. The objective of this study was to to examine the utility of direct fluorescence visualization (FV) by dental practitioners as an aid in decision-making during screening for cancer and other oral lesions. METHODS: Dentists were trained to use a stepwise protocol for evaluation of the oral mucosa: medical history, head, neck and oral exam, and fluorescent visualization exam. They were asked to use clinical features to categorize lesions as low (LR), intermediate (IR), or high (HR) risk and then to determine FV status of these lesions. Clinicians made the decision of which lesions to reassess in 3 weeks and based on this reassessment, to refer forward. RESULTS: Of 2404 patients screened over 11 months, 357 initially had lesions with 325 (15%) identified as LR, 16 (4.5%) IR, and 16 (4.5%) HR. Lesions assessed initially as IR and HR had a 2.7-fold increased risk of FV loss persisting to the reassessment appointment versus the LR lesions. The most predictive model for lesion persistence included both FV status and lesion risk assessment. CONCLUSION: A protocol for screening (assess risk, reassess, and refer) is recommended for the screening of abnormal intraoral lesions. Integrating FV into a process of assessing and reassessing lesions significantly improved this model.

Pre-profiling factors influencing serum microRNA levels.

BACKGROUND: MicroRNAs (miRNAs) are non-coding RNAs that negatively regulate gene expression by preventing the translation of specific mRNA transcripts. Recent studies have shown that miRNAs are stably expressed in human serum samples, making them good candidates for the non-invasive detection of disease. However, before circulating miRNAs can be used reliably as biomarkers of disease, the pre-measurement variables that may affect serum miRNA levels must be assessed. METHODS: In this study we used quantitative RT-PCR to examine the effect of hemolysis, fasting, and smoking on the levels of 742 miRNAs in the serum of healthy individuals. We also compared serum miRNA profiles of samples taken from healthy individuals over different time periods to assess normal serum miRNA fluctuations. RESULTS: We have found that mechanical hemolysis of blood samples can significantly alter serum miRNA quantification and have identified 162 miRNAs that are significantly up-regulated in hemolysed serum samples. Conversely, fasting and smoking were demonstrated to not have a significant effect on the overall serum miRNA profiles of healthy individuals. The serum miRNA profiles of matched samples taken from individuals over varying time periods showed a high correlation and no miRNAs were significantly differentially expressed in these samples further suggesting the utility of serum miRNAs as biomarkers of disease. Taking the above results into consideration, we have identified miR-99a-5p and miR-139-5p as novel endogenous controls for serum miRNA studies due to their consistency across all sample sets. CONCLUSION: These results identify important pre-profiling factors that should be taken into consideration when identifying endogenous controls and candidate biomarkers for circulating miRNA studies.

Three-Dimension Epithelial Segmentation in Optical Coherence Tomography of the Oral Cavity Using Deep Learning.

This paper aims to simplify the application of optical coherence tomography (OCT) for the examination of subsurface morphology in the oral cavity and reduce barriers towards the adoption of OCT as a biopsy guidance device. The aim of this work was to develop automated software tools for the simplified analysis of the large volume of data collected during OCT. Imaging and corresponding histopathology were acquired in-clinic using a wide-field endoscopic OCT system. An annotated dataset (n = 294 images) from 60 patients (34 male and 26 female) was assembled to train four unique neural networks. A deep learning pipeline was built using convolutional and modified u-net models to detect the imaging field of view (network 1), detect artifacts (network 2), identify the tissue surface (network 3), and identify the presence and location of the epithelial-stromal boundary (network 4). The area under the curve of the image and artifact detection networks was 1.00 and 0.94, respectively. The Dice similarity score for the surface and epithelial-stromal boundary segmentation networks was 0.98 and 0.83, respectively. Deep learning (DL) techniques can identify the location and variations in the epithelial surface and epithelial-stromal boundary in OCT images of the oral mucosa. Segmentation results can be synthesized into accessible en face maps to allow easier visualization of changes.

DNA Ploidy as a Potential Adjunct Prognostic Marker of Low-Risk Prostate Cancer Progression after Radical Prostatectomy.

PURPOSE: Post prostatectomy PSA kinetics and General Grade Groups (GGG) are the strongest prognostic markers of biochemical recurrence (BCR) and prostate cancer (PCa)-specific mortality after radical prostatectomy. Despite having low-risk PCa, some patients will experience BCR, for some, clinically significant BCR. There is a need for an objective prognostic marker at the time of prostatectomy to improve risk stratification within this population. In this study, we investigated the prognostic potential of DNA ploidy. MATERIALS AND METHODS: Prostatectomy samples from 97 patients with GGG1 and GGG2 with a low-risk CAPRA-S score were included in this study. PCa tissue with the worst Gleason pattern underwent tissue disaggregation, cell isolation and staining with a DNA stoichiometric stain. Using image cytometry, DNA ploidy was measured and a Ploidy Score (PS) was generated. RESULTS: Among the 97 patients, 79 had no BCR, 18 experienced BCR, of which 14 had a PSA doubling time (PSA-DT) >1 year (low-risk group) and 4 had a PSA-DT of <1 year (high-risk group). Using Logistic regression analysis, only pathological T stage (pT) and PS independently predicted BCR with PS being the most significant (p = 0.001). The number of aneuploid cells was significantly higher in the high-risk group compared to the other groups (p = 1.7x10-11). PS combined with GGG diagnosis further stratified risk groups of biochemical recurrence free survival within CAPRA-S low-risk cohort. CONCLUSION: DNA ploidy is an independent prognostic marker of BCR in low-risk PCa after radical prostatectomy, which could early on identify potentially aggressive PCa recurrences and introduce a more personalized approach to salvage treatments.

A tiling resolution DNA microarray with complete coverage of the human genome.

We constructed a tiling resolution array consisting of 32,433 overlapping BAC clones covering the entire human genome. This increases our ability to identify genetic alterations and their boundaries throughout the genome in a single comparative genomic hybridization (CGH) experiment. At this tiling resolution, we identified minute DNA alterations not previously reported. These alterations include microamplifications and deletions containing oncogenes, tumor-suppressor genes and new genes that may be associated with multiple tumor types. Our findings show the need to move beyond conventional marker-based genome comparison approaches, that rely on inference of continuity between interval markers. Our submegabase resolution tiling set for array CGH (SMRT array) allows comprehensive assessment of genomic integrity and thereby the identification of new genes associated with disease.

Prediction of Disease Progression to Upfront Pembrolizumab Monotherapy in Advanced Non-Small-Cell Lung Cancer with High PD-L1 Expression Using Baseline CT Disease Quantification and Smoking Pack Years.

Health Canada approved pembrolizumab in the first-line setting for advanced non-small-cell lung cancer with PD-L1 ≥ 50% and no EGFR/ALK aberration. The keynote 024 trial showed 55% of such patients progress with pembrolizumab monotherapy. We propose that the combination of baseline CT and clinical factors can help identify those patients who may progress. In 138 eligible patients from our institution, we retrospectively collected their baseline variables, including baseline CT findings (primary lung tumor size and metastatic site), smoking pack years, performance status, tumor pathology, and demographics. The treatment response was assessed via RECIST 1.1 using the baseline and first follow-up CT. Associations between the baseline variables and progressive disease (PD) were tested by logistic regression analyses. The results showed 46/138 patients had PD. The baseline CT "number of involved organs" by metastasis and smoking pack years were independently associated with PD (p < 0.05), and the ROC analysis showed a good performance of the model that integrated these variables in predicting PD (AUC: 0.79). This pilot study suggests that the combination of baseline CT disease and smoking PY can identify who may progress on pembrolizumab monotherapy and can potentially facilitate decision-making for the optimal first-line treatment in the high PD-L1 cohort.

Delineating spatial cell-cell interactions in the solid tumour microenvironment through the lens of highly multiplexed imaging.

The growth and metastasis of solid tumours is known to be facilitated by the tumour microenvironment (TME), which is composed of a highly diverse collection of cell types that interact and communicate with one another extensively. Many of these interactions involve the immune cell population within the TME, referred to as the tumour immune microenvironment (TIME). These non-cell autonomous interactions exert substantial influence over cell behaviour and contribute to the reprogramming of immune and stromal cells into numerous pro-tumourigenic phenotypes. The study of some of these interactions, such as the PD-1/PD-L1 axis that induces CD8+ T cell exhaustion, has led to the development of breakthrough therapeutic advances. Yet many common analyses of the TME either do not retain the spatial data necessary to assess cell-cell interactions, or interrogate few (<10) markers, limiting the capacity for cell phenotyping. Recently developed digital pathology technologies, together with sophisticated bioimage analysis programs, now enable the high-resolution, highly-multiplexed analysis of diverse immune and stromal cell markers within the TME of clinical specimens. In this article, we review the tumour-promoting non-cell autonomous interactions in the TME and their impact on tumour behaviour. We additionally survey commonly used image analysis programs and highly-multiplexed spatial imaging technologies, and we discuss their relative advantages and limitations. The spatial organization of the TME varies enormously between patients, and so leveraging these technologies in future studies to further characterize how non-cell autonomous interactions impact tumour behaviour may inform the personalization of cancer treatment.​.