Interobserver reproducibility of cervical histology interpretation with and without p16 immunohistochemistry

OBJECTIVES

Histopathological diagnosis of colposcopically identified cervical lesions is a critical step for the recognition of cervical cancer precursors requiring treatment. Although there have been efforts to standardize the histologic diagnosis of cervical biopsy specimens, in terms of terminology and use of biomarkers, there is no uniform approach in the pathology community. Adjunctive p16 immunohistochemistry (IHC) can highlight precancer diagnoses, with use recommendations outlined by the Lower Anogenital Squamous Terminology project.

METHODS

We assessed the diagnostic reproducibility of cervical histopathological biopsy specimens with and without p16 staining among 2 expert pathologists.

RESULTS

Interpretation of p16 IHC as positive vs negative was highly reproducible (92.5% agreement, κ = 0.85); greater variation was seen in the choice of which biopsy specimens required adjunctive p16 staining (78.0% agreement, κ = 0.43). Adjunctive p16 IHC did not significantly increase diagnostic agreement under multitiered grading systems (benign vs cervical intraepithelial neoplasia [CIN] 1/low-grade squamous intraepithelial lesion vs atypical squamous metaplasia vs CIN2/high-grade squamous intraepithelial lesion [HSIL] vs CIN3/HSIL-CIN3 vs cancer) (65.5% agreement, κ = 0.56 without p16; 70.0% agreement, κ = 0.58 with p16). However, when dichotomizing diagnoses based on clinical management (less than HSIL vs HSIL+), diagnostic agreement increased with p16 IHC (90.5% agreement, κ = 0.79 without p16; 92.0% agreement, κ = 0.84 with p16). For biopsy specimens taken from women positive for human papillomavirus (HPV) type 16, agreement was similar with or without adjunctive p16 (κ = 0.80 without p16; κ = 0.78-0.80 with p16). In contrast, p16 IHC substantially improved diagnostic agreement for cervical biopsy specimens taken from women positive for other high-risk HPV strains, producing improvements in κ from 0.03 to 0.24.

CONCLUSIONS

Adjunctive p16 immunostaining provides useful information in the evaluation of cervical biopsies for precancer. In our study, we have demonstrated that it is highly reproducible between 2 pathologists, although the decision of which biopsies warrant its use is less so. Furthermore, although p16 IHC showed a limited increase in diagnostic reproducibility for all biopsies included in our study, it did demonstrate a more sizable gain in biopsies negative for HPV 16 but positive for other high-risk genotypes. Further studies are needed to clarify the role of p16 IHC and how it can be optimized for the detection of cervical precancer, particularly in HPV-vaccinated populations where types other than HPV 16 are relatively more important.

CNN stability training improves robustness to scanner and IHC-based image variability for epithelium segmentation in cervical histology

Background

In digital pathology, image properties such as color, brightness, contrast and blurriness may vary based on the scanner and sample preparation. Convolutional Neural Networks (CNNs) are sensitive to these variations and may underperform on images from a different domain than the one used for training. Robustness to these image property variations is required to enable the use of deep learning in clinical practice and large scale clinical research.

Aims

CNN Stability Training (CST) is proposed and evaluated as a method to increase CNN robustness to scanner and Immunohistochemistry (IHC)-based image variability.

Methods

CST was applied to segment epithelium in immunohistological cervical Whole Slide Images (WSIs). CST randomly distorts input tiles and factors the difference between the CNN prediction for the original and distorted inputs within the loss function. CNNs were trained using 114 p16-stained WSIs from the same scanner, and evaluated on 6 WSI test sets, each with 23 to 24 WSIs of the same tissue but different scanner/IHC combinations. Relative robustness (rAUC) was measured as the difference between the AUC on the training domain test set (i.e., baseline test set) and the remaining test sets.

Results

Across all test sets, The AUC of CST models outperformed "No CST" models (AUC: 0.940-0.989 vs. 0.905-0.986, p < 1e - 8), and obtained an improved robustness (rAUC: [-0.038, -0.003] vs. [-0.081, -0.002]). At a WSI level, CST models showed an increase in performance in 124 of the 142 WSIs. CST models also outperformed models trained with random on-the-fly data augmentation (DA) in all test sets ([0.002, 0.021], p < 1e-6).

Conclusion

CST offers a path to improve CNN performance without the need for more data and allows customizing distortions to specific use cases. A python implementation of CST is publicly available at https://github.com/TIGACenter/CST_v1.

Automated evaluation of p16/Ki-67 dual stain cytology as a biomarker for detection of anal precancer in MSM living with HIV

BACKGROUND

Human papillomavirus (HPV)-related biomarkers such as p16/Ki-67 "dual stain" (DS) cytology have shown promising clinical performance for anal cancer screening. Here, we assessed the performance of automated evaluation of DS cytology to detect anal precancer in men who have sex with men (MSM) living with human immunodeficiency virus (HIV).

METHODS

We conducted a cross-sectional analysis of 320 MSM with HIV undergoing anal cancer screening and high-resolution anoscopy (HRA) between 2009-2010. We evaluated the performance of automated DS based on a deep-learning classifier compared to manual DS evaluation to detect anal intraepithelial neoplasia (AIN) grades 2 and 3 (AIN2+) and AIN3. We evaluated different DS-positive cell thresholds quantified by the automated approach and modeled its performance compared to other screening strategies in a hypothetical population of MSM with HIV.

RESULTS

Automated DS had significantly higher specificity (50.9% vs. 42.2%, p=0.0004) and similar sensitivity (93.2% vs. 92.1%) for detection of AIN2+ compared to manual DS cytology. HPV testing with automated DS triage was significantly more specific (56.5% vs. 50.9%, p=0.0003) with the same sensitivity (93.2%) compared to automated DS alone. In a modeled analysis assuming a 20% AIN2+ prevalence, automated DS detected the most precancers compared to manual DS and anal cytology (186,184, and 162, respectively) and had the lowest HRA referral per AIN2+ detected (3.1,3.5, and 3.3, respectively).

CONCLUSIONS

Compared to manual evaluation, automated DS cytology detects the same number of precancers with lower HRA referral.

Co-expression of YAP and TAZ associates with chromosomal instability in human cholangiocarcinoma

Background

Activation of the oncogene yes-associated protein (YAP) is frequently detected in intrahepatic cholangiocarcinoma (iCCA); however, the expression pattern and the functional impact of its paralogue WW domain-containing transcription regulator 1 (WWTR1; synonym: TAZ) are not well described in different CCA subtypes.

Methods

Immunohistochemical analysis of YAP and TAZ in iCCA and extrahepatic CCA (eCCA) cohorts was performed. YAP/TAZ shuttling and their functional impact on CCA cell lines were investigated. Target genes expression after combined YAP/TAZ inhibition was analyzed.

Results

Immunohistochemical analysis of iCCA and eCCA revealed YAP or TAZ positivity in up to 49.2%; however, oncogene co-expression was less frequent (up to 23%). In contrast, both proteins were jointly detectable in most CCA cell lines and showed nuclear/cytoplasmic shuttling in a cell density-dependent manner. Next to the pro-proliferative function of YAP/TAZ, both transcriptional co-activators cooperated in the regulation of a gene signature that indicated the presence of chromosomal instability (CIN). A correlation between YAP and the CIN marker phospho-H2A histone family member X (pH2AX) was particularly observed in tissues from iCCA and distal CCA (dCCA). The presence of the CIN genes in about 25% of iCCA was statistically associated with worse prognosis.

Conclusions

YAP and TAZ activation is not uncoupled from cell density in CCA cells and both factors cooperatively contribute to proliferation and expression of CIN-associated genes. The corresponding group of CCA patients is characterized by CIN and may benefit from YAP/TAZ-directed therapies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-021-08794-5.

Therapeutic Effects of Inhibition of Sphingosine-1-Phosphate Signaling in HIF-2α Inhibitor-Resistant Clear Cell Renal Cell Carcinoma

Simple Summary

Clear cell renal cell carcinoma is a common malignancy that represents 80% of all kidney tumors. Most tumors harbor an inactivation of the VHL gene, leading to the accumulation of HIF-1α and HIF-2α. Promising clinical results of specific HIF-2α inhibitors will soon lead to new treatment options for advanced cancer patients, although primary and acquired resistance to these agents are common. We here show that Acriflavine, which inhibits both HIF-1α and HIF-2α, and Fingolimod (FTY720), which inhibits sphingosine-1-phosphate signaling, show therapeutic activities in several experimental ccRCC models that are resistant to HIF-2α-inhibitor treatment. Additionally, we show that specific HIF-2α-inhibition suppresses the tumor immune microenvironment, which will be important to consider for future combination studies with immune checkpoint inhibitors.

Abstract

Specific inhibitors of HIF-2α have recently been approved for the treatment of ccRCC in VHL disease patients and have shown encouraging results in clinical trials for metastatic sporadic ccRCC. However, not all patients respond to therapy and pre-clinical and clinical studies indicate that intrinsic as well as acquired resistance mechanisms to HIF-2α inhibitors are likely to represent upcoming clinical challenges. It would be desirable to have additional therapeutic options for the treatment of HIF-2α inhibitor resistant ccRCCs. Here we investigated the effects on tumor growth and on the tumor microenvironment of three different direct and indirect HIF-α inhibitors, namely the HIF-2α-specific inhibitor PT2399, the dual HIF-1α/HIF-2α inhibitor Acriflavine, and the S1P signaling pathway inhibitor FTY720, in the autochthonous Vhl/Trp53/Rb1 mutant ccRCC mouse model and validated these findings in human ccRCC cell culture models. We show that FTY720 and Acriflavine exhibit therapeutic activity in several different settings of HIF-2α inhibitor resistance. We also identify that HIF-2α inhibition strongly suppresses T cell activation in ccRCC. These findings suggest prioritization of sphingosine pathway inhibitors for clinical testing in ccRCC patients and also suggest that HIF-2α inhibitors may inhibit anti-tumor immunity and might therefore be contraindicated for combination therapies with immune checkpoint inhibitors.

Accuracy and Efficiency of Deep-Learning-Based Automation of Dual Stain Cytology in Cervical Cancer Screening

Background

With the advent of primary human papillomavirus testing followed by cytology for cervical cancer screening, visual interpretation of cytology slides remains the last subjective analysis step and suffers from low sensitivity and reproducibility.

Methods

We developed a cloud-based whole-slide imaging platform with a deep-learning classifier for p16/Ki-67 dual-stained (DS) slides trained on biopsy-based gold standards. We compared it with conventional Pap and manual DS in 3 epidemiological studies of cervical and anal precancers from Kaiser Permanente Northern California and the University of Oklahoma comprising 4253 patients. All statistical tests were 2-sided.

Results

In independent validation at Kaiser Permanente Northern California, artificial intelligence (AI)-based DS had lower positivity than cytology (P < .001) and manual DS (P < .001) with equal sensitivity and substantially higher specificity compared with both Pap (P < .001) and manual DS (P < .001), respectively. Compared with Pap, AI-based DS reduced referral to colposcopy by one-third (41.9% vs 60.1%, P < .001). At a higher cutoff, AI-based DS had similar performance to high-grade squamous intraepithelial lesions cytology, indicating a risk high enough to allow for immediate treatment. The classifier was robust, showing comparable performance in 2 cytology systems and in anal cytology.

Conclusions

Automated DS evaluation removes the remaining subjective component from cervical cancer screening and delivers consistent quality for providers and patients. Moving from Pap to automated DS substantially reduces the number of colposcopies and also achieves excellent performance in a simulated fully vaccinated population. Through cloud-based implementation, this approach is globally accessible. Our results demonstrate that AI not only provides automation and objectivity but also delivers a substantial benefit for women by reduction of unnecessary colposcopies.

Common clonal origin of conventional T cells and induced regulatory T cells in breast cancer patients

Regulatory CD4+ T cells (Treg) prevent tumor clearance by conventional T cells (Tconv) comprising a major obstacle of cancer immune-surveillance. Hitherto, the mechanisms of Treg repertoire formation in human cancers remain largely unclear. Here, we analyze Treg clonal origin in breast cancer patients using T-Cell Receptor and single-cell transcriptome sequencing. While Treg in peripheral blood and breast tumors are clonally distinct, Tconv clones, including tumor-antigen reactive effectors (Teff), are detected in both compartments. Tumor-infiltrating CD4+ cells accumulate into distinct transcriptome clusters, including early activated Tconv, uncommitted Teff, Th1 Teff, suppressive Treg and pro-tumorigenic Treg. Trajectory analysis suggests early activated Tconv differentiation either into Th1 Teff or into suppressive and pro-tumorigenic Treg. Importantly, Tconv, activated Tconv and Treg share highly-expanded clones contributing up to 65% of intratumoral Treg. Here we show that Treg in human breast cancer may considerably stem from antigen-experienced Tconv converting into secondary induced Treg through intratumoral activation.

The local immunological microenvironment in colorectal cancer as a prognostic factor for treatment decisions in the clinic: The way ahead

Analysis of the local immunological microenvironment in colorectal cancer lesions yielded prognostic markers. Harnessing these insights for clinical application however requires the use of sophisticated technology and algorithms, especially the robust and reproducible quantification of immune cells. These technologies are available and will allow individualized treatment decisions beyond the current standard.

Sequential metastases of colorectal cancer: Immunophenotypes and spatial distributions of infiltrating immune cells in relation to time and treatments

The role of the immune system in the course of colorectal cancer has been elucidated in the last decade. While quantification of immune cell infiltrates within the resected specimen at diagnosis has a clear power to estimate the prognosis of the patient, the role of infiltrating immune cells within the metastatic situation and especially within the metastatic lesion itself requires further detailed analyses. Recent analyses of infiltrates in colorectal cancer liver metastases revealed a role for the infiltrate density not only for prognosis but also in the prediction of treatment response. This not only broadens the view on these infiltrates and indicates a systematic role of the local immunological microenvironment, but also raises the question how these infiltrates change during repeated courses of treatment (i.e., resection, chemotherapy, etc.). To address this question, sequential lung or sequential liver metastases of colorectal cancer patients were analyzed using whole slide image quantification after immunohistochemical staining against CD3, CD8, FOXP3, CD68 and Granzyme B. The clinical data and interventions were associated with each individual patient and the metastatic lesions. The resulting cell densities reveal a heterogeneous profile: after successful treatment of a metastatic lesion, the recurrent lesion can still have the same immunophenotype with similar cell distributions. In a situation of a favorable immune cell profile, this profile can return and apparently convey a similar favorable course throughout the disease. But also the opposite was found: the recurrent metastatic lesion could have a different profile with alterations in specific immune cell subsets over time. Further analyses are required to elucidate the different patterns and their associations to the treatment, the tumor cell phenotype and other dynamic factors. However, it is clear from this data however, that there is an immune cell plasticity that needs to be analyzed for individual patients.

Precision medicine: myocardial fibrosis burden and genotype predict outcome in non-ischemic dilated cardiomyopathy (DCM)

Introduction

Myocardial fibrosis occurs during pathological remodeling of the heart and can be associated with worse outcome in affected patients. Genetic background is also known to affect patients' survival. In this study we sought to estimate patients' overall fibrosis burden by combining independent modalities including left ventricular endomyocardial biopsy (LV-EMB), circulating biomarkers, and cMRI. We also aimed to use patients' genetics information to predict outcome. Furthermore, we investigated the correlation between cardiac fibrosis and genetic variations to detect novel susceptibility loci for fibrosis in DCM.

Methods

A total number of 542 DCM patients were included. Collagen volume fraction (CVF) was automatically estimated from biopsies. 13 circulating fibrosis biomarkers were measured using Human Magnetic Luminex Screening Assays, and the cMRIs were screened for presence of LGE. Whole exome sequencing (WES) was performed in 410 patients of the cohort using illumina HiSeq 2000. Common (MAF ≥0.05 in the study population OR gnomAD NFE AF ≥0.01) and non-common missense variants (MAF &lt;0.05 in the study population AND gnomAD NFE AF &lt;0.01) in 42 DCM genes were tested for associations with end points using single variants and burden analyses respectively. Analyses were adjusted for age and sex and performed using R and SKAT. End points were all-cause mortality and a composite of heart failure (HF) associated events.

Results

The median follow-up time was 43.2 months (2084 patient-years). Sixty-two patients reached the composite end point and 55 died. LV-EMB proved to be a safe procedure with a total complication rate of 2.3%. Machine learning based characterization of biopsies was highly accurate. Although the 3 different modalities did not significantly correlate with one another, the extent of CVF, levels of MMP-2, TIMP-1, OPN, and GDF-15, and presence of LGE were each significantly associated with worse outcome. Four possible susceptibility loci for cardiac fibrosis in DCM were introduced and underwent eQTL analyses. Rare missense variants in a list of 11 DCM-related genes showed to be associated with the 2 outcome measures or fibrosis burden.

Conclusions

LV-EMB, fibrosis biomarkers, and cMRI likely capture different aspects of a detrimental fibrosis process and may be combined to estimate patients' prognosis and monitor therapeutic success. Phenotype-genotype association studies help elucidate novel disease pathomechanisms and individualize patients' treatment.

Funding Acknowledgement

Type of funding source: Other. Main funding source(s): German Centre for Cardiovascular Research: DZHK

The impact of SAMHD1 expression and mutation status in mantle cell lymphoma: An analysis of the MCL Younger and Elderly trial

The sterile alpha motif and histidine-aspartic domain-containing protein 1 (SAMHD1) has been demonstrated to predict the response to high-dose cytarabine consolidation treatment in acute myeloid leukemia patients. Here, we evaluated SAMHD1 as potential biomarker for the response to high-dose cytarabine in mantle cell lymphoma (MCL) patients. We quantified SAMHD1 protein expression and determined the mutation status in patients of the MCL Younger and Elderly trials (n = 189), who had received high-dose cytarabine- or fludarabine-based polychemotherapy. Additionally, we quantified SAMHD1 expression in B cell lymphoma cell lines and exposed them to cytarabine, fludarabine, and clinically relevant combinations. Across both trials investigated, SAMHD1 mutations had a frequency of 7.1% (n = 13) and did not significantly affect the failure-free survival (FFS, P = .47). In patients treated with high-dose cytarabine- or fludarabine-containing regimes, SAMHD1 expression was not significantly associated with FFS or complete remission rate. SAMHD1 expression in B cell lymphoma cell lines, however, inversely correlated with their in vitro response to cytarabine as single agent (R = .65, P = .0065). This correlation could be reversed by combining cytarabine with other chemotherapeutics, such as oxaliplatin and vincristine, similar to the treatment regime of the MCL Younger trial. We conclude that this might explain why we did not observe a significant association between SAMHD1 protein expression and the outcome of MCL patients upon cytarabine-based treatment.

Surfactant Expression Defines an Inflamed Subtype of Lung Adenocarcinoma Brain Metastases that Correlates with Prolonged Survival

PURPOSE

To provide a better understanding of the interplay between the immune system and brain metastases to advance therapeutic options for this life-threatening disease.

EXPERIMENTAL DESIGN

Tumor-infiltrating lymphocytes (TIL) were quantified by semiautomated whole-slide analysis in brain metastases from 81 lung adenocarcinomas. Multi-color staining enabled phenotyping of TILs (CD3, CD8, and FOXP3) on a single-cell resolution. Molecular determinants of the extent of TILs in brain metastases were analyzed by transcriptomics in a subset of 63 patients. Findings in lung adenocarcinoma brain metastases were related to published multi-omic primary lung adenocarcinoma The Cancer Genome Atlas data (n = 230) and single-cell RNA-sequencing (scRNA-seq) data (n = 52,698).

RESULTS

TIL numbers within tumor islands was an independent prognostic marker in patients with lung adenocarcinoma brain metastases. Comparative transcriptomics revealed that expression of three surfactant metabolism-related genes (SFTPA1, SFTPB, and NAPSA) was closely associated with TIL numbers. Their expression was not only prognostic in brain metastasis but also in primary lung adenocarcinoma. Correlation with scRNA-seq data revealed that brain metastases with high expression of surfactant genes might originate from tumor cells resembling alveolar type 2 cells. Methylome-based estimation of immune cell fractions in primary lung adenocarcinoma confirmed a positive association between lymphocyte infiltration and surfactant expression. Tumors with a high surfactant expression displayed a transcriptomic profile of an inflammatory microenvironment.

CONCLUSIONS

The expression of surfactant metabolism-related genes (SFTPA1, SFTPB, and NAPSA) defines an inflamed subtype of lung adenocarcinoma brain metastases characterized by high abundance of TILs in close vicinity to tumor cells, a prolonged survival, and a tumor microenvironment which might be more accessible to immunotherapeutic approaches.

High endothelial venules are associated with microsatellite instability, hereditary background and immune evasion in colorectal cancer

BACKGROUND

Microsatellite-unstable (MSI) tumours show a high load of mutational neoantigens, as a consequence of DNA mismatch repair deficiency. Consequently, MSI tumours commonly present with dense immune infiltration and develop immune evasion mechanisms. Whether improved lymphocyte recruitment contributes to the pronounced immune infiltration in MSI tumours is unknown. We analysed the density of high endothelial venules (HEV) and postcapillary blood vessels specialised for lymphocyte trafficking, in MSI colorectal cancers (CRC).

METHODS

HEV density was determined by immunohistochemical staining of FFPE tissue sections from MSI (n = 48) and microsatellite-stable (MSS, n = 35) CRCs. Associations with clinical and pathological variables were analysed.

RESULTS

We found elevated HEV densities in MSI compared with MSS CRCs (median 0.049 vs 0.000 counts/mm2, respectively, p = 0.0002), with the highest densities in Lynch syndrome MSI CRCs. Dramatically elevated HEV densities were observed in B2M-mutant Lynch syndrome CRCs, pointing towards a link between lymphocyte recruitment and immune evasion (median 0.485 vs 0.0885 counts/mm2 in B2M-wild-type tumours, p = 0.0237).

CONCLUSIONS

Our findings for the first time indicate a significant contribution of lymphocyte trafficking in immune responses against MSI CRC, particularly in the context of Lynch syndrome. High HEV densities in B2M-mutant tumours underline the significance of immunoediting during tumour evolution.

Robust and accurate quantification of biomarkers of immune cells in lung cancer micro-environment using deep convolutional neural networks

Recent years have seen a growing awareness of the role the immune system plays in successful cancer treatment, especially in novel therapies like immunotherapy. The characterization of the immunological composition of tumors and their micro-environment is thus becoming a necessity. In this paper we introduce a deep learning-based immune cell detection and quantification method, which is based on supervised learning, i.e., the input data for training comprises labeled images. Our approach objectively deals with staining variation and staining artifacts in immunohistochemically stained lung cancer tissue and is as precise as humans. This is evidenced by the low cell count difference to humans of 0.033 cells on average. This method, which is based on convolutional neural networks, has the potential to provide a new quantitative basis for research on immunotherapy.

Large-scale in-silico identification of a tumor-specific antigen pool for targeted immunotherapy in triple-negative breast cancer

Since the advent of cetuximab, clinical cancer treatment has evolved from the standard, relatively nonspecific chemo- and radiotherapy with significant cytotoxic side effects towards immunotherapeutic approaches with selective, target-mechanism-based effects. Antibody therapies as the most successful form of cancer immunotherapy led to approved treatments for specific cancer types with increased patient survival. Thus, the identification of tumor antigens with high immunogenicity is in central focus now. In this study, we applied computational methods to comprehensively discover overexpressed molecular targets with high therapeutic relevance for clinical, immunotherapeutic cancer treatment in triple-negative breast cancer (TNBC). By actively modeling potential negative side effects utilizing expression data of 29 different, normal human tissues, we were able to develop a highly-specific coverage of TNBC patients with RNA targets. We identified here more than 400 potential tumor-specific antigens suitable for targeted therapy, including several already identified as potential targets for TNBC and other solid tumors. A specific cocktail of MAGEB4, CT83, TLX3, ACTL8, PRDM13 achieved almost 94% patient coverage in TNBC. Overall, these results show that our approach can identify and prioritize TNBC targets suitable for targeted therapy. Therefore, our method has the potential to lead to new and more effective immunotherapeutic cancer treatment.

Retinal stem cells modulate proliferative parameters to coordinate post-embryonic morphogenesis in the eye of fish

Combining clonal analysis with a computational agent based model, we investigate how tissue-specific stem cells for neural retina (NR) and retinal pigmented epithelium (RPE) of the teleost medaka (Oryzias latipes) coordinate their growth rates. NR cell division timing is less variable, consistent with an upstream role as growth inducer. RPE cells divide with greater variability, consistent with a downstream role responding to inductive signals. Strikingly, the arrangement of the retinal ciliary marginal zone niche results in a spatially biased random lineage loss, where stem- and progenitor cell domains emerge spontaneously. Further, our data indicate that NR cells orient division axes to regulate organ shape and retinal topology. We highlight an unappreciated mechanism for growth coordination, where one tissue integrates cues to synchronize growth of nearby tissues. This strategy may enable evolution to modulate cell proliferation parameters in one tissue to adapt whole-organ morphogenesis in a complex vertebrate organ.

By the time babies reach adulthood, they have grown many times larger than they were at birth. This development is driven by an increase in the number and size of cells in the body. In particular, special types of cells, called stem cells, act as a reservoir for tissues: they divide to create new cells that will mature into various specialized structures.

The retina is the light-sensitive part of the eye. It consists of the neural retina, a tissue that contains light-detecting cells, which is supported by the retinal pigment epithelium or RPE. In fish, the RPE and neural retina are replenished by distinct groups of stem cells that do not mix, despite the tissues being close together.

Unlike humans, fish grow throughout adulthood, and their eyes must then keep pace with the body. This means that the different tissues in the retina must somehow coordinate to expand at the same rate: otherwise, the retina would get wrinkled and not work properly. Tsingos et al. therefore wanted to determine how stem cells in the neural retina and RPE co-operated to produce the right number of new cells at the right time.

First, stem cells in the eyes of newly hatched fish were labelled with a visible marker so that their divisions could be tracked over time to build cell family trees. This showed that stem cells behaved differently in the neural retina and the RPE. Computer simulations of the growing retina explained this behavior: stem cells in the neural retina were telling the RPE stem cells when it was time to divide. Combining results from the simulations with data from the experiments revealed that a stem cell decided to keep up dividing partly because of its position in the tissue, and partly because of random chance.

To be healthy, the body needs to fine-tune the number of cells it produces: creating too few cells may make it difficult to heal after injury, but making too many could lead to diseases such as cancer. Understanding how tissues normally agree to grow together could therefore open new avenues of treatment for these conditions.

Comparison of the Reverse-Remodeling Effect of Pharmacological Soluble Guanylate Cyclase Activation With Pressure Unloading in Pathological Myocardial Left Ventricular Hypertrophy

Background: Pressure unloading induces the regression of left ventricular myocardial hypertrophy (LVH). Recent findings indicate that pharmacological activation of the soluble guanylate cyclase (sGC) – cyclic guanosine monophosphate (cGMP) pathway may also exert reverse-remodeling properties in the myocardium. Therefore, we aimed to investigate the effects of the sGC activator cinaciguat in a rat model of LVH and compare it to the “gold standard” pressure unloading therapy.

Methods: Abdominal aortic banding was performed for 6 or 12 weeks. Sham operated animals served as controls. Pressure unloading was induced by removing the aortic constriction after week 6. The animals were treated from week 7 to 12, with 10 mg/kg/day cinaciguat or with placebo p.o., respectively. Cardiac function and morphology were assessed by left ventricular pressure-volume analysis and echocardiography. Additionally, key markers of myocardial hypertrophy, fibrosis, nitro-oxidative stress, apoptosis and cGMP signaling were analyzed.

Results: Pressure unloading effectively reversed LVH, decreased collagen accumulation and provided protection against oxidative stress and apoptosis. Regression of LVH was also associated with a full recovery of cardiac function. In contrast, chronic activation of the sGC enzyme by cinaciguat at sustained pressure overload only slightly influenced pre-established hypertrophy. However, it led to increased PKG activity and had a significant impact on interstitial fibrosis, nitro-oxidative stress and apoptosis. Amelioration of the pathological structural alterations prevented the deterioration of LV systolic function (contractility and ejection fraction) and improved myocardial stiffness.

Conclusion: Our results indicate that both cinaciguat treatment and pressure unloading evoked anti-remodeling effects and improved LV function, however in a differing manners.

Correlation between genomic index lesions and mpMRI and 68Ga-PSMA-PET/CT imaging features in primary prostate cancer

Magnetic resonance imaging (MRI) and prostate specific membrane antigen (PSMA)- positron emission tomography (PET)/computed tomography (CT)-imaging of prostate cancer (PCa) are emerging techniques to assess the presence of significant disease and tumor progression. It is not known, however, whether and to what extent lesions detected by these imaging techniques correlate with genomic features of PCa. The aim of this study was therefore to define a genomic index lesion based on chromosomal copy number alterations (CNAs) as marker for tumor aggressiveness in prostate biopsies in direct correlation to multiparametric (mp) MRI and ⁶⁸Ga-PSMA-PET/CT imaging features. CNA profiles of 46 biopsies from five consecutive patients with clinically high-risk PCa were obtained from radiologically suspicious and unsuspicious areas. All patients underwent mpMRI, MRI/TRUS-fusion biopsy, ⁶⁸Ga-PSMA-PET/CT and a radical prostatectomy. CNAs were directly correlated to imaging features and radiogenomic analyses were performed. Highly significant CNAs (≥10 Mbp) were found in 22 of 46 biopsies. Chromosome 8p, 13q and 5q losses were the most common findings. There was an strong correspondence between the radiologic and the genomic index lesions. The radiogenomic analyses suggest the feasibility of developing radiologic signatures that can distinguish between genomically more or less aggressive lesions. In conclusion, imaging features of mpMRI and ⁶⁸Ga-PSMA-PET/CT can guide to the genomically most aggressive lesion of a PCa. Radiogenomics may help to better differentiate between indolent and aggressive PCa in the future.

Teledermatology: Comparison of Store-and-Forward Versus Live Interactive Video Conferencing

A decreasing number of dermatologists and an increasing number of patients in Western countries have led to a relative lack of clinicians providing expert dermatologic care. This, in turn, has prolonged wait times for patients to be examined, putting them at risk. Store-and-forward teledermatology improves patient access to dermatologists through asynchronous consultations, reducing wait times to obtain a consultation. However, live video conferencing as a synchronous service is also frequently used by practitioners because it allows immediate interaction between patient and physician. This raises the question of which of the two approaches is superior in terms of quality of care and convenience. There are pros and cons for each in terms of technical requirements and features. This viewpoint compares the two techniques based on a literature review and a clinical perspective to help dermatologists assess the value of teledermatology and determine which techniques would be valuable in their practice.

Skin Cancer Classification Using Convolutional Neural Networks: Systematic Review

Background

State-of-the-art classifiers based on convolutional neural networks (CNNs) were shown to classify images of skin cancer on par with dermatologists and could enable lifesaving and fast diagnoses, even outside the hospital via installation of apps on mobile devices. To our knowledge, at present there is no review of the current work in this research area.

Objective

This study presents the first systematic review of the state-of-the-art research on classifying skin lesions with CNNs. We limit our review to skin lesion classifiers. In particular, methods that apply a CNN only for segmentation or for the classification of dermoscopic patterns are not considered here. Furthermore, this study discusses why the comparability of the presented procedures is very difficult and which challenges must be addressed in the future.

Methods

We searched the Google Scholar, PubMed, Medline, ScienceDirect, and Web of Science databases for systematic reviews and original research articles published in English. Only papers that reported sufficient scientific proceedings are included in this review.

Results

We found 13 papers that classified skin lesions using CNNs. In principle, classification methods can be differentiated according to three principles. Approaches that use a CNN already trained by means of another large dataset and then optimize its parameters to the classification of skin lesions are the most common ones used and they display the best performance with the currently available limited datasets.

Conclusions

CNNs display a high performance as state-of-the-art skin lesion classifiers. Unfortunately, it is difficult to compare different classification methods because some approaches use nonpublic datasets for training and/or testing, thereby making reproducibility difficult. Future publications should use publicly available benchmarks and fully disclose methods used for training to allow comparability.

Altered density, composition and microanatomical distribution of infiltrating immune cells in cutaneous squamous cell carcinoma of organ transplant recipients

BACKGROUND

The inflammatory tumour microenvironment is crucial for effective tumour control, and long-term immunosuppression has been identified as a major risk factor for skin carcinogenesis. In solid organ transplant recipients (OTRs) undergoing long-term pharmacological immunosuppression, an increased incidence of cutaneous squamous cell carcinoma (SCC) and more aggressive tumour growth compared with immunocompetent patients has been reported.

OBJECTIVES

To determine the density and phenotype of immune cells infiltrating SCC and surrounding skin in OTRs, and to characterize the microanatomical distribution patterns in comparison with immunocompetent patients.

METHODS

We analysed immune cell infiltrates within SCC and at defined regions of interest (ROIs) of tumour-surrounding skin in formalin-fixed paraffin-embedded tissue of 20 renal transplant patients and 18 carefully matched immunocompetent patients by high-resolution semiautomated microscopy on complete tissue sections stained for CD4, CD8, CD20 and CD68.

RESULTS

The overall immune cell density of SCC arising in OTRs was significantly reduced compared with immunocompetent patients. Particularly CD4⁺ infiltrates at the directly invasive margin and tumour vicinity, intratumoral CD8⁺ T-cell densities and the overall density of CD20⁺ tumour-infiltrating B cells were significantly reduced in the tissue of OTRs.

CONCLUSIONS

Immune cell infiltrates within SCC and at defined ROIs of tumour-surrounding skin in OTRs differ markedly in their composition and microanatomical distribution compared with tumours arising in immunocompetent patients. Our findings substantially broaden the understanding of how long-term systemic immunosuppression modulates the local inflammatory microenvironment in the skin and at the site of invasive SCC.

Identification of CRKII, CFL1, CNTN1, NME2, and TKT as Novel and Frequent T-Cell Targets in Human IDH-Mutant Glioma

Purpose: Successful immunotherapies for IDH^mut gliomas require better knowledge of T-cell target antigens. Here, we elucidated their antigen repertoire recognized by spontaneous T-cell responses using an unbiased proteomic approach.Experimental Design: Protein fractionations of tissue lysates from IDH^mut gliomas (n = 4) were performed. Fractions were tested by IFNγ ELISpot assay for recognition through patients' T cells. Proteins of immunogenic fractions were identified by mass spectrometry and validated by in silico-predicted synthetic long peptides in patients of origin, additional IDH^mut glioma patients (n = 16), and healthy donors (n = 13). mRNA and protein expression of immunogenic antigens was analyzed in tumor tissues and IDH^mut glioma stem-like cells (GSC). HLA-A*02-restricted T-cell epitopes were functionally determined by short peptides and numbers of antigen-specific T cells by HLA-peptide tetramer analysis.Results: A total of 2,897 proteins were identified in immunogenic tumor fractions. Based on a thorough filter process, 79 proteins were selected as potential T-cell antigens. Twenty-six of these were recognized by the patients' T cells, and five of them (CRKII, CFL1, CNTN1, NME2, and TKT) in up to 56% unrelated IDH^mut glioma patients. Most immunogenic tumor-associated antigens (TAA) were expressed in IDH^mut gliomas and GSCs, while being almost absent in normal brain tissues. Finally, we identified HLA-A*02-restricted epitopes for CRKII, NME2, and TKT that were recognized by up to 2.82% of antigen-specific peripheral cytotoxic T cells in IDH^mut glioma patients.Conclusions: By analyzing the repertoire of T-cell target antigens in IDH^mut glioma patients, we identified five novel immunogenic TAAs and confirmed their expression on IDH^mut tumors and GSCs. Clin Cancer Res; 24(12); 2951-62. ©2018 AACR.

Transcriptomic analysis of aggressive meningiomas identifies PTTG1 and LEPR as prognostic biomarkers independent of WHO grade

Meningiomas are frequent central nervous system tumors. Although most meningiomas are benign (WHO grade I) and curable by surgery, WHO grade II and III tumors remain therapeutically challenging due to frequent recurrence. Interestingly, relapse also occurs in some WHO grade I meningiomas. Hence, we investigated the transcriptional features defining aggressive (recurrent, malignantly progressing or WHO grade III) meningiomas in 144 cases. Meningiomas were categorized into non-recurrent (NR), recurrent (R), and tumors undergoing malignant progression (M) in addition to their WHO grade. Unsupervised transcriptomic analysis in 62 meningiomas revealed transcriptional profiles lining up according to WHO grade and clinical subgroup. Notably aggressive subgroups (R+M tumors and WHO grade III) shared a large set of differentially expressed genes (n=332; p<0.01, FC>1.25). In an independent multicenter validation set (n=82), differential expression of 10 genes between WHO grades was confirmed. Additionally, among WHO grade I tumors differential expression between NR and aggressive R+M tumors was affirmed for PTTG1, AURKB, ECT2, UBE2C and PRC1, while MN1 and LEPR discriminated between NR and R+M WHO grade II tumors. Univariate survival analysis revealed a significant association with progression-free survival for PTTG1, LEPR, MN1, ECT2, PRC1, COX10, UBE2C expression, while multivariate analysis identified a prediction for PTTG1 and LEPR mRNA expression independent of gender, WHO grade and extent of resection. Finally, stainings of PTTG1 and LEPR confirmed malignancy-associated protein expression changes. In conclusion, based on the so far largest study sample of WHO grade III and recurrent meningiomas we report a comprehensive transcriptional landscape and two prognostic markers.

Tumor Cell Load and Heterogeneity Estimation From Diffusion-Weighted MRI Calibrated With Histological Data: an Example From Lung Cancer

Diffusion-weighted magnetic resonance imaging (DWI) is a key non-invasive imaging technique for cancer diagnosis and tumor treatment assessment, reflecting Brownian movement of water molecules in tissues. Since densely packed cells restrict molecule mobility, tumor tissues produce usually higher signal (a.k.a. less attenuated signal) on isotropic maps compared with normal tissues. However, no general quantitative relation between DWI data and the cell density has been established. In order to link low-resolution clinical cross-sectional data with high-resolution histological information, we developed an image processing and analysis chain, which was used to study the correlation between the diffusion coefficient (D value) estimated from DWI and tumor cellularity from serial histological slides of a resected non-small cell lung cancer tumor. Color deconvolution followed by cell nuclei segmentation was performed on digitized histological images to determine local and cell-type specific 2d (two-dimensional) densities. From these, the 3d cell density was inferred by a model-based sampling technique, which is necessary for the calculation of local and global 3d tumor cell count. Next, DWI sequence information was overlaid with high-resolution CT data and the resected histology using prominent anatomical hallmarks for co-registration of histology tissue blocks and non-invasive imaging modalities' data. The integration of cell numbers information and DWI data derived from different tumor areas revealed a clear negative correlation between cell density and D value. Importantly, spatial tumor cell density can be calculated based on DWI data. In summary, our results demonstrate that tumor cell count and heterogeneity can be predicted from DWI data, which may open new opportunities for personalized diagnosis and therapy optimization.

High numbers of PDCD1 (PD-1)-positive T cells and B2M mutations in microsatellite-unstable colorectal cancer

DNA mismatch repair (MMR)-deficient cancers accumulate high numbers of coding microsatellite mutations, which lead to the generation of highly immunogenic frameshift peptide (FSP) neoantigens. MMR-deficient cells can grow out to clinically manifest cancers either if they evade immune cell attack or if local T-cells get exhausted. Therefore, a subset of MSI cancer patients responds particularly well to treatment with immune checkpoint inhibitors. We analyzed whether immune evasion in MMR-deficient cancer mediated by loss of HLA class I or II antigens is related to local immune cell activation status. Microsatellites located in Beta2-microglobulin (B2M) and the HLA class II-regulatory genes RFX5 and CIITA were analyzed for mutations in MMR-deficient colorectal cancers (n = 53). The results were related to CD3-positive and PDCD1 (PD-1)-positive T-cell infiltration. PDCD1 (PD-1)-positive T-cell counts were significantly higher in B2M-mutant compared to B2M-wild type tumors (median: 22.2 cells per 0.25 mm² vs. 2.0 cells per 0.25 mm², Wilcoxon test p = 0.002). Increasing PDCD1 (PD-1)-positive T-cell infiltration was significantly related to an increased likelihood of B2M mutations (OR = 1.81). HLA class II antigen expression status was significantly associated with enhanced overall T-cell infiltration, but not related to PDCD1 (PD-1)-positive T-cells. These results suggest that immune evasion mediated by B2M mutation-induced loss of HLA class I antigen expression predominantly occurs in an environment of activated PDCD1 (PD-1)-positive T cell infiltration. If B2M mutations interfere with anti-PDCD1 (PD-1)/CD274 (PD-L1) therapy success, we predict that resistance towards anti-PDCD1 (PD-1) therapy may - counterintuitively - be particularly common in patients with MMR-deficient cancers that show high PDCD1 (PD-1)-positive T cell infiltration.

Glyoxalase 1 expression is associated with an unfavorable prognosis of oropharyngeal squamous cell carcinoma

Background

Glyoxalase 1 is a key enzyme in the detoxification of reactive metabolites such as methylglyoxal and induced Glyoxalase 1 expression has been demonstrated for several human malignancies. However, the regulation and clinical relevance of Glyoxalase 1 in the context of head and neck squamous cell carcinoma has not been addressed so far.

Methods

Argpyrimidine modification as a surrogate for methylglyoxal accumulation and Glyoxalase 1 expression in tumor cells was assessed by immunohistochemical staining of tissue microarrays with specimens from oropharyngeal squamous cell carcinoma patients (n = 154). Prognostic values of distinct Glyoxalase 1 staining patterns were demonstrated by Kaplan-Meier, univariate and multivariate Cox proportional hazard model analysis. The impact of exogenous methylglyoxal or a Glyoxalase 1 inhibitor on the viability of two established tumor cell lines was monitored by a colony-forming assay in vitro.

Results

Glyoxalase 1 expression in tumor cells of oropharyngeal squamous cell carcinoma patients was positively correlated with the presence of Argpyrimidine modification and administration of exogenous methylglyoxal induced Glyoxalase 1 protein levels in FaDu and Cal27 cells in vitro. Cal27 cells with lower basal and methylglyoxal-induced Glyoxalase 1 expression were more sensitive to the cytotoxic effect at high methylgyoxal concentrations and both cell lines showed a decrease in colony formation with increasing amounts of a Glyoxalase 1 inhibitor. A high and nuclear Glyoxalase 1 staining was significantly correlated with shorter progression-free and disease-specific survival, and served as an independent risk factor for an unfavorable prognosis of oropharyngeal squamous cell carcinoma patients.

Conclusions

Induced Glyoxalase 1 expression is a common feature in the pathogenesis of oropharyngeal squamous cell carcinoma and most likely represents an adaptive response to the accumulation of cytotoxic metabolites. Oropharyngeal squamous cell carcinoma patients with a high and nuclear Glyoxalase 1 staining pattern have a high risk for treatment failure, but might benefit from pharmacological targeting Glyoxalase 1 activity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-017-3367-5) contains supplementary material, which is available to authorized users.

Quantitative diagnosis of breast tumors by morphometric classification of microenvironmental myoepithelial cells using a machine learning approach

Machine learning systems have recently received increased attention for their broad applications in several fields. In this study, we show for the first time that histological types of breast tumors can be classified using subtle morphological differences of microenvironmental myoepithelial cell nuclei without any direct information about neoplastic tumor cells. We quantitatively measured 11661 nuclei on the four histological types: normal cases, usual ductal hyperplasia and low/high grade ductal carcinoma in situ (DCIS). Using a machine learning system, we succeeded in classifying the four histological types with 90.9% accuracy. Electron microscopy observations suggested that the activity of typical myoepithelial cells in DCIS was lowered. Through these observations as well as meta-analytic database analyses, we developed a paracrine cross-talk-based biological mechanism of DCIS progressing to invasive cancer. Our observations support novel approaches in clinical computational diagnostics as well as in therapy development against progression.

F-18 labelled PSMA-1007: biodistribution, radiation dosimetry and histopathological validation of tumor lesions in prostate cancer patients

PURPOSE

The prostate-specific membrane antigen (PSMA) targeted positron-emitting-tomography (PET) tracer 68Ga-PSMA-11 shows great promise in the detection of prostate cancer. However, 68Ga has several shortcomings as a radiolabel including short half-life and non-ideal energies, and this has motivated consideration of 18F-labelled analogs. 18F-PSMA-1007 was selected among several 18F-PSMA-ligand candidate compounds because it demonstrated high labelling yields, outstanding tumor uptake and fast, non-urinary background clearance. Here, we describe the properties of 18F-PSMA-1007 in human volunteers and patients.

METHODS

Radiation dosimetry of 18F-PSMA-1007 was determined in three healthy volunteers who underwent whole-body PET-scans and concomitant blood and urine sampling. Following this, ten patients with high-risk prostate cancer underwent 18F-PSMA-1007 PET/CT (1 h and 3 h p.i.) and normal organ biodistribution and tumor uptakes were examined. Eight patients underwent prostatectomy with extended pelvic lymphadenectomy. Uptake in intra-prostatic lesions and lymph node metastases were correlated with final histopathology, including PSMA immunostaining.

RESULTS

With an effective dose of approximately 4.4-5.5 mSv per 200-250 MBq examination, 18F-PSMA-1007 behaves similar to other PSMA-PET agents as well as to other 18F-labelled PET-tracers. In comparison to other PSMA-targeting PET-tracers, 18F-PSMA-1007 has reduced urinary clearance enabling excellent assessment of the prostate. Similar to 18F-DCFPyL and with slightly slower clearance kinetics than PSMA-11, favorable tumor-to-background ratios are observed 2-3 h after injection. In eight patients, diagnostic findings were successfully validated by histopathology. 18F-PSMA-1007 PET/CT detected 18 of 19 lymph node metastases in the pelvis, including nodes as small as 1 mm in diameter.

CONCLUSION

18F-PSMA-1007 performs at least comparably to 68Ga-PSMA-11, but its longer half-life combined with its superior energy characteristics and non-urinary excretion overcomes some practical limitations of 68Ga-labelled PSMA-targeted tracers.

A 3D self-organizing multicellular epidermis model of barrier formation and hydration with realistic cell morphology based on EPISIM

The epidermis and the stratum corneum (SC) as its outermost layer have evolved to protect the body from evaporative water loss to the environment. To morphologically represent the extremely flattened cells of the SC - and thereby the epidermal barrier - in a multicellular computational model, we developed a 3D biomechanical model (BM) based on ellipsoid cell shapes. We integrated the BM in the multicellular modelling and simulation platform EPISIM. We created a cell behavioural model (CBM) with EPISIM encompassing regulatory feedback loops between the epidermal barrier, water loss to the environment, and water and calcium flow within the tissue. This CBM allows a small number of stem cells to initiate self-organizing epidermal stratification, yielding the spontaneous emergence of water and calcium gradients comparable to experimental data. We find that the 3D in silico epidermis attains homeostasis most quickly at high ambient humidity, and once in homeostasis the epidermal barrier robustly buffers changes in humidity. Our model yields an in silico epidermis with a previously unattained realistic morphology, whose cell neighbour topology is validated with experimental data obtained from in vivo images. This work paves the way to computationally investigate how an impaired SC barrier precipitates disease.

Spatial intratumoral heterogeneity of proliferation in immunohistochemical images of solid tumors

PURPOSE

The interactions of neoplastic cells with each other and the microenvironment are complex. To understand intratumoral heterogeneity, subtle differences should be quantified. Main factors contributing to heterogeneity include the gradient ischemic level within neoplasms, action of microenvironment, mechanisms of intercellular transfer of genetic information, and differential mechanisms of modifications of genetic material/proteins. This may reflect on the expression of biomarkers in the context of prognosis/stratification. Hence, a rigorous approach for assessing the spatial intratumoral heterogeneity of histological biomarker expression with accuracy and reproducibility is required, since patterns in immunohistochemical images can be challenging to identify and describe.

METHODS

A quantitative method that is useful for characterizing complex irregular structures is lacunarity; it is a multiscale technique that exhaustively samples the image, while the decay of its index as a function of window size follows characteristic patterns for different spatial arrangements. In histological images, lacunarity provides a useful measure for the spatial organization of a biomarker when a sampling scheme is employed and relevant features are computed. The proposed approach quantifies the segmented proliferative cells and not the textural content of the histological slide, thus providing a more realistic measure of heterogeneity within the sample space of the tumor region. The aim is to investigate in whole sections of primary pancreatic neuroendocrine neoplasms (pNENs), using whole-slide imaging and image analysis, the spatial intratumoral heterogeneity of Ki-67 immunostains. Unsupervised learning is employed to verify that the approach can partition the tissue sections according to distributional heterogeneity.

RESULTS

The architectural complexity of histological images has shown that single measurements are often insufficient. Inhomogeneity of distribution depends not only on percentage content of proliferation phase but also on how the phase fills the space. Lacunarity curves demonstrate variations in the sampled image sections. Since the spatial distribution of proliferation in each case is different, the width of the curves changes too. Image sections that have smaller numerical variations in the computed features correspond to neoplasms with spatially homogeneous proliferation, while larger variations correspond to cases where proliferation shows various degrees of clumping. Grade 1 (uniform/nonuniform: 74%/26%) and grade 3 (uniform: 100%) pNENs demonstrate a more homogeneous proliferation with grade 1 neoplasms being more variant, while grade 2 tumor regions render a more diverse landscape (50%/50%). Hence, some cases show an increased degree of spatial heterogeneity comparing to others with similar grade. Whether this is a sign of different tumor biology and an association with a more benign/malignant clinical course needs to be investigated further. The extent and range of spatial heterogeneity has the potential to be evaluated as a prognostic marker.

CONCLUSIONS

The association with tumor grade as well as the rationale that the methodology reflects true tumor architecture supports the technical soundness of the method. This reflects a general approach which is relevant to other solid tumors and biomarkers. Drawing upon the merits of computational biomedicine, the approach uncovers salient features for use in future studies of clinical relevance.

Reduced Microvascular Density in Omental Biopsies of Children with Chronic Kidney Disease

BACKGROUND

Endothelial dysfunction is an early manifestation of cardiovascular disease (CVD) and consistently observed in patients with chronic kidney disease (CKD). We hypothesized that CKD is associated with systemic damage to the microcirculation, preceding macrovascular pathology. To assess the degree of "uremic microangiopathy", we have measured microvascular density in biopsies of the omentum of children with CKD.

PATIENTS AND METHODS

Omental tissue was collected from 32 healthy children (0-18 years) undergoing elective abdominal surgery and from 23 age-matched cases with stage 5 CKD at the time of catheter insertion for initiation of peritoneal dialysis. Biopsies were analyzed by independent observers using either a manual or an automated imaging system for the assessment of microvascular density. Quantitative immunohistochemistry was performed for markers of autophagy and apoptosis, and for the abundance of the angiogenesis-regulating proteins VEGF-A, VEGF-R2, Angpt1 and Angpt2.

RESULTS

Microvascular density was significantly reduced in uremic children compared to healthy controls, both by manual imaging with a digital microscope (median surface area 0.61% vs. 0.95%, p<0.0021 and by automated quantification (total microvascular surface area 0.89% vs. 1.17% p = 0.01). Density measured by manual imaging was significantly associated with age, height, weight and body surface area in CKD patients and healthy controls. In multivariate analysis, age and serum creatinine level were the only independent, significant predictors of microvascular density (r2 = 0.73). There was no immunohistochemical evidence for apoptosis or autophagy. Quantitative staining showed similar expression levels of the angiogenesis regulators VEGF-A, VEGF-receptor 2 and Angpt1 (p = 0.11), but Angpt2 was significantly lower in CKD children (p = 0.01).

CONCLUSIONS

Microvascular density is profoundly reduced in omental biopsies of children with stage 5 CKD and associated with diminished Angpt2 signaling. Microvascular rarefaction could be an early systemic manifestation of CKD-induced cardiovascular disease.

Abstract A085: Immune evasion and PD-1-positive T cell infiltration in DNA mismatch repair-deficient colorectal cancer

Background and Aims: DNA mismatch repair-deficient cancers show microsatellite instability (MSI) and accumulate high numbers of somatic insertion/deletion mutations at repetitive sequence stretches. The high mutational load of MSI cancers leads to the generation of multiple frameshift peptide (FSP) neoantigens, which are highly immunogenic and recognized by the immune system as foreign antigens. MSI cancers respond particularly well to treatment with immune checkpoint inhibitors such as anti-PD-1 antibodies. We systematically analyzed immune evasion phenomena in MSI cancers and aimed to identify factors that are related to and potentially determine PD-1 expression on tumor-infiltrating lymphocytes, focusing on alterations of HLA class I and II antigen presentation pathways in cancer cells.

Methods: Microsatellites located in Beta2-microglobulin (B2M) and the HLA class II-regulatory genes RFX5 and CIITA were analyzed for mutations in MSI colorectal cancer specimens (n = 53). HLA class I and II antigen expression was examined by immunohistochemistry. In addition, tumor-infiltrating lymphocytes (CD3-positive T cells, PD-1-positive T cells) were quantified using a semi-automated system.

Results: We related B2M mutation and HLA class II antigen expression status of MSI colorectal cancer specimens (n = 56) to CD3- and PD-1 positive T cell infiltration in the tumor. PD-1-positive T cell infiltration was significantly higher in B2M-mutant (mt) compared to B2M-wild type (wt) tumors (median: 22.2 cells per 0.25 mm2 in B2M-mt vs. 2.0 cells per 0.25 mm2 in B2M-wt, Wilcoxon's rank sum test p = 0.002). Increasing PD-1-positive T cell infiltration was significantly related to an increased likelihood of B2M mutation and loss of HLA class I antigen expression (OR = 1.81). In contrast, HLA class II antigen expression status was not related to the proportion of PD-1-positive lymphocytes, but significantly associated with enhanced overall T cell infiltration.

Conclusions: These results suggest that immune evasion mediated by B2M mutation-induced loss of HLA class I antigen expression predominantly occurs in an environment of activated PD-1-positive T cell infiltration, supporting the validity of the immunoediting concept in MSI colorectal cancers. Moreover, B2M mutation status may predict therapy resistance against anti-PD-1 therapy.

Note: This abstract was not presented at the conference.

Citation Format: Matthias Kloor, Jonas Janikovits, Julia Krzykalla, Axel Benner, Niels Grabe, Fabian Echterdiek, Meike Mueller, Sandrina Koerner, Aysel Ahadova, Magnus von Knebel Doeberitz. Immune evasion and PD-1-positive T cell infiltration in DNA mismatch repair-deficient colorectal cancer [abstract]. In: Proceedings of the Second CRI-CIMT-EATI-AACR International Cancer Immunotherapy Conference: Translating Science into Survival; 2016 Sept 25-28; New York, NY. Philadelphia (PA): AACR; Cancer Immunol Res 2016;4(11 Suppl):Abstract nr A085.

Spatial transcriptome analysis reveals Notch pathway-associated prognostic markers in IDH1 wild-type glioblastoma involving the subventricular zone

BACKGROUND

The spatial relationship of glioblastoma (GBM) to the subventricular zone (SVZ) is associated with inferior patient survival. However, the underlying molecular phenotype is largely unknown. We interrogated an SVZ-dependent transcriptome and potential location-specific prognostic markers.

METHODS

mRNA microarray data of a discovery set (n = 36 GBMs) were analyzed for SVZ-dependent gene expression and process networks using the MetaCore™ workflow. Differential gene expression was confirmed by qPCR in a validation set of 142 IDH1 wild-type GBMs that was also used for survival analysis.

RESULTS

Microarray analysis revealed a transcriptome distinctive of SVZ+ GBM that was enriched for genes associated with Notch signaling. No overlap was found to The Cancer Genome Atlas's molecular subtypes. Independent validation of SVZ-dependent expression confirmed four genes with simultaneous prognostic impact: overexpression of HES4 (p = 0.034; HR 1.55) and DLL3 (p = 0.017; HR 1.61) predicted inferior, and overexpression of NTRK2 (p = 0.049; HR 0.66) and PIR (p = 0.025; HR 0.62) superior overall survival (OS). Additionally, overexpression of DLL3 was predictive of shorter progression-free survival (PFS) (p = 0.043; HR 1.64). Multivariate analysis revealed overexpression of HES4 to be independently associated with inferior OS (p = 0.033; HR 2.03), and overexpression of DLL3 with inferior PFS (p = 0.046; HR 1.65).

CONCLUSIONS

We identified four genes with SVZ-dependent expression and prognostic significance, among those HES4 and DLL3 as part of Notch signaling, suggesting further evaluation of location-tailored targeted therapies.

Identification of immunotherapeutic targets by genomic profiling of rectal NET metastases

Neuroendocrine tumors (NETs) of the gastrointestinal tract are a rare and heterogeneous group of neoplasms with unique tumor biology and clinical management issues. While surgery is the only curative treatment option in patients with early stage NETs, the optimal management strategy for patients with advanced metastatic NETs is unknown. Based on the tremendous success of immunotherapeutic approaches, we sought to investigate such approaches in a case of metastatic rectal NET. Here, we apply an integrative approach using various computational and experimental methods to explore several aspects of the tumor-host immune interactions for immunotherapeutic options. Sequencing of six different liver metastases revealed a quite homogenous set of mutations, and further analysis of these mutations for immunogenicity revealed few neo-epitopes with pre-existing T cell reactivity, which can be used in therapeutic vaccines. Staining for immunomodulatory proteins and cytokine profiling showed that the immune setting is surprisingly different, when compared to liver metastases of colorectal cancer for instance. Taken together, our results highlight the broad range and complexity of tumor-host immune interaction and underline the value of an integrative approach.

In-depth tissue profiling using multiplexed immunohistochemical consecutive staining on single slide

Despite remarkable recent achievements of immunotherapy strategies in cancer treatment, clinical responses remain limited to subsets of patients. Predictive markers of disease course and response to immunotherapy are urgently needed. Recent results have revealed the potential predictive value of immune cell phenotype and spatial distribution at the tumor site, prompting the need for multidimensional immunohistochemical analyses of tumor tissues. To address this need, we developed a sample-sparing, highly multiplexed immunohistochemistry technique based on iterative cycles of tagging, image scanning, and destaining of chromogenic substrate on a single slide. This assay, in combination with a newly developed automated digital landscaping solution, democratizes access to high-dimensional immunohistochemical analyses by capturing the complexity of the immunome using routine pathology standards. Applications of the method extend beyond cancer to screen and validate comprehensive panels of tissue-based prognostic and predictive markers, perform in-depth in situ monitoring of therapies, and identify targets of disease.

Spatial niche formation but not malignant progression is a driving force for intratumoural heterogeneity

Intratumoural heterogeneity (ITH) is a major cause of cancer-associated lethality. Extensive genomic ITH has previously been reported in clear cell renal cell carcinoma (ccRCC). Here we address the question whether ITH increases with malignant progression and can hence be exploited as a prognostic marker. Unexpectedly, precision quantitative image analysis reveals that the degree of functional ITH is virtually identical between primary ccRCCs of the lowest stage and advanced, metastatic tumours. Functional ITH was found to show a stage-independent topological pattern with peak proliferative and signalling activities almost exclusively in the tumour periphery. Exome sequencing of matching peripheral and central primary tumour specimens reveals various region-specific mutations. However, these mutations cannot directly explain the zonal pattern suggesting a role of microenvironmental factors in shaping functional ITH. In conclusion, our results indicate that ITH is an early and general characteristic of malignant growth rather than a consequence of malignant progression.

It has been increasingly recognised that tumours are not made up of a homogeneous population of cells. Here, the authors show heterogeneous expression of five protein markers in renal cell cancer and demonstrate that the progression of the tumour does not influence the degree of heterogeneity in the tumour.

Oral treatment with a zinc complex of acetylsalicylic acid prevents diabetic cardiomyopathy in a rat model of type-2 diabetes: activation of the Akt pathway

BACKGROUND

Type-2 diabetics have an increased risk of cardiomyopathy, and heart failure is a major cause of death among these patients. Growing evidence indicates that proinflammatory cytokines may induce the development of insulin resistance, and that anti-inflammatory medications may reverse this process. We investigated the effects of the oral administration of zinc and acetylsalicylic acid, in the form of bis(aspirinato)zinc(II)-complex Zn(ASA)2, on different aspects of cardiac damage in Zucker diabetic fatty (ZDF) rats, an experimental model of type-2 diabetic cardiomyopathy.

METHODS

Nondiabetic control (ZL) and ZDF rats were treated orally with vehicle or Zn(ASA)2 for 24 days. At the age of 29-30 weeks, the electrical activities, left-ventricular functional parameters and left-ventricular wall thicknesses were assessed. Nitrotyrosine immunohistochemistry, TUNEL-assay, and hematoxylin-eosin staining were performed. The protein expression of the insulin-receptor and PI3K/AKT pathway were quantified by Western blot.

RESULTS

Zn(ASA)2-treatment significantly decreased plasma glucose concentration in ZDF rats (39.0 ± 3.6 vs 49.4 ± 2.8 mM, P < 0.05) while serum insulin-levels were similar among the groups. Data from cardiac catheterization showed that Zn(ASA)2 normalized the increased left-ventricular diastolic stiffness (end-diastolic pressure-volume relationship: 0.064 ± 0.008 vs 0.084 ± 0.014 mmHg/µl; end-diastolic pressure: 6.5 ± 0.6 vs 7.9 ± 0.7 mmHg, P < 0.05). Furthermore, ECG-recordings revealed a restoration of prolonged QT-intervals (63 ± 3 vs 83 ± 4 ms, P < 0.05) with Zn(ASA)2. Left-ventricular wall thickness, assessed by echocardiography, did not differ among the groups. However histological examination revealed an increase in the cardiomyocytes' transverse cross-section area in ZDF compared to the ZL rats, which was significantly decreased after Zn(ASA)2-treatment. Additionally, a significant fibrotic remodeling was observed in the diabetic rats compared to ZL rats, and Zn(ASA)2-administered ZDF rats showed a similar collagen content as ZL animals. In diabetic hearts Zn(ASA)2 significantly decreased DNA-fragmentation, and nitro-oxidative stress, and up-regulated myocardial phosphorylated-AKT/AKT protein expression. Zn(ASA)2 reduced cardiomyocyte death in a cellular model of oxidative stress. Zn(ASA)2 had no effects on altered myocardial CD36, GLUT-4, and PI3K protein expression.

CONCLUSIONS

We demonstrated that treatment of type-2 diabetic rats with Zn(ASA)2 reduced plasma glucose-levels and prevented diabetic cardiomyopathy. The increased myocardial AKT activation could, in part, help to explain the cardioprotective effects of Zn(ASA)2. The oral administration of Zn(ASA)2 may have therapeutic potential, aiming to prevent/treat cardiac complications in type-2 diabetic patients.

Tumoral Immune Cell Exploitation in Colorectal Cancer Metastases Can Be Targeted Effectively by Anti-CCR5 Therapy in Cancer Patients

The immune response influences the clinical course of colorectal cancer (CRC). Analyzing the invasive margin of human CRC liver metastases, we identified a mechanism of immune cell exploitation by tumor cells. While two distinct subsets of myeloid cells induce an influx of T cells into the invasive margin via CXCL9/CXCL10, CCL5 is produced by these T cells and stimulates pro-tumoral effects via CCR5. CCR5 blockade in patient-derived functional in vitro organotypic culture models showed a macrophage repolarization with anti-tumoral effects. These anti-tumoral effects were then confirmed in a phase I trial with a CCR5 antagonist in patients with liver metastases of advanced refractory CRC. Mitigation of tumor-promoting inflammation within the tumor tissue and objective tumor responses in CRC were observed.

Combined Treatment of ATRA with Epigenetic Drugs Increases Aggressiveness of Glioma Xenografts

BACKGROUND/AIM

Recently, anti-tumourigenic effects of all-trans-retinoic-acid (ATRA) on glioblastoma stem cells were demonstrated. Therefore we investigated if these beneficial effects could be enhanced by co-medication with epigenetic drugs such as the histone deacetylase (HDAC) inhibitor suberoylanilide hydroxamic acid (SAHA) or the DNA-methyltransferase inhibitor 5-aza-2'deoxycytidine (5-AZA).

MATERIALS AND METHODS

Glioma stem cell xenografts were treated for 42 days with ATRA plus SAHA or ATRA plus 5-AZA or the correspondent monotherapies. Tumour sizes, histological features, proliferation and apoptosis rates were assessed.

RESULTS

Neither SAHA nor 5-AZA were able to enhance the anti-tumourigenic effect of ATRA. Instead, tumours became more aggressive. Combination of ATRA plus 5-AZA increased tumour size (p<0.05) and induced more frequent and larger necroses (p<0.05) and tumours were more invasive (p<0.05) in comparison to controls. A similar trend was observed for the combination of ATRA plus SAHA.

CONCLUSION

Combining ATRA with epigenetic drug therapies led to the unwanted opposite effect and increased aggressiveness of glioma xenografts, arguing against future clinical applications of such combinations.

LOC283731 promoter hypermethylation prognosticates survival after radiochemotherapy in IDH1 wild-type glioblastoma patients

MGMT promoter methylation status is currently the only established molecular prognosticator in IDH wild-type glioblastoma multiforme (GBM). Therefore, we aimed to discover novel therapy-associated epigenetic biomarkers. After enrichment for hypermethylated fractions using methyl-CpG-immunoprecipitation (MCIp), we performed global DNA methylation profiling for 14 long-term (LTS; >36 months) and 15 short-term (STS; 6-10 months) surviving GBM patients. Even after exclusion of the G-CIMP phenotype, we observed marked differences between the LTS and STS methylome. A total of 1,247 probes in 706 genes were hypermethylated in LTS and 463 probes in 305 genes were found to be hypermethylated in STS patients (p values < 0.05, log2 fold change ± 0.5). We identified 13 differentially methylated regions (DMRs) with a minimum of four differentially methylated probes per gene. Indeed, we were able to validate a subset of these DMRs through a second, independent method (MassARRAY) in our LTS/STS training set (ADCY1, GPC3, LOC283731/ISLR2). These DMRs were further assessed for their prognostic capability in an independent validation cohort (n = 62) of non-G-CIMP GBMs from the TCGA. Hypermethylation of multiple CpGs mapping to the promoter region of LOC283731 correlated with improved patient outcome (p = 0.03). The prognostic performance of LOC283731 promoter hypermethylation was confirmed in a third independent study cohort (n = 89), and was independent of gender, performance (KPS) and MGMT status (p = 0.0485, HR = 0.63). Intriguingly, the prediction was most pronounced in younger GBM patients (<60 years). In conclusion, we provide compelling evidence that promoter methylation status of this novel gene is a prognostic biomarker in IDH1 wild-type/non-G-CIMP GBMs.

Directed random walks and constraint programming reveal active pathways in hepatocyte growth factor signaling

An effective means to analyze mRNA expression data is to take advantage of established knowledge from pathway databases, using methods such as pathway-enrichment analyses. However, pathway databases are not case-specific and expression data could be used to infer gene-regulation patterns in the context of specific pathways. In addition, canonical pathways may not always describe the signaling mechanisms properly, because interactions can frequently occur between genes in different pathways. Relatively few methods have been proposed to date for generating and analyzing such networks, preserving the causality between gene interactions and reasoning over the qualitative logic of regulatory effects. We present an algorithm (MCWalk) integrated with a logic programming approach, to discover subgraphs in large-scale signaling networks by random walks in a fully automated pipeline. As an exemplary application, we uncover the signal transduction mechanisms in a gene interaction network describing hepatocyte growth factor-stimulated cell migration and proliferation from gene-expression measured with microarray and RT-qPCR using in-house perturbation experiments in a keratinocyte-fibroblast co-culture. The resulting subgraphs illustrate possible associations of hepatocyte growth factor receptor c-Met nodes, differentially expressed genes and cellular states. Using perturbation experiments and Answer Set programming, we are able to select those which are more consistent with the experimental data. We discover key regulator nodes by measuring the frequency with which they are traversed when connecting signaling between receptors and significantly regulated genes and predict their expression-shift consistently with the measured data. The Java implementation of MCWalk is publicly available under the MIT license at: https://bitbucket.org/akittas/biosubg.

Patch-Based Nonlinear Image Registration for Gigapixel Whole Slide Images

OBJECTIVE

Image registration of whole slide histology images allows the fusion of fine-grained information-like different immunohistochemical stains-from neighboring tissue slides. Traditionally, pathologists fuse this information by looking subsequently at one slide at a time. If the slides are digitized and accurately aligned at cell level, automatic analysis can be used to ease the pathologist's work. However, the size of those images exceeds the memory capacity of regular computers.

METHODS

We address the challenge to combine a global motion model that takes the physical cutting process of the tissue into account with image data that is not simultaneously globally available. Typical approaches either reduce the amount of data to be processed or partition the data into smaller chunks to be processed separately. Our novel method first registers the complete images on a low resolution with a nonlinear deformation model and later refines this result on patches by using a second nonlinear registration on each patch. Finally, the deformations computed on all patches are combined by interpolation to form one globally smooth nonlinear deformation. The NGF distance measure is used to handle multistain images.

RESULTS

The method is applied to ten whole slide image pairs of human lung cancer data. The alignment of 85 corresponding structures is measured by comparing manual segmentations from neighboring slides. Their offset improves significantly, by at least 15%, compared to the low-resolution nonlinear registration.

CONCLUSION/SIGNIFICANCE

The proposed method significantly improves the accuracy of multistain registration which allows us to compare different antibodies at cell level.

Low density of FOXP3-positive T cells in normal colonic mucosa is related to the presence of beta2-microglobulin mutations in Lynch syndrome-associated colorectal cancer

Microsatellite instability (MSI-H) is caused by DNA mismatch repair deficiency and occurs in 15% of colorectal cancers. MSI-H cancers generate highly immunogenic frameshift peptide (FSP) antigens, which elicit pronounced local immune responses. A subset of MSI-H colorectal cancers develops in frame of Lynch syndrome, which represents an ideal human model for studying the concept of immunoediting. Immunoediting describes how continuous anti-tumoral immune surveillance of the host eventually leads to the selection of tumor cells that escape immune cell recognition and destruction. Between 30 and 40% of Lynch syndrome-associated colorectal cancers display loss of HLA class I antigen expression as a result of Beta2-microglobulin (B2M) mutations. Whether B2M mutations result from immunoediting has been unknown. To address this question, we related B2M mutation status of Lynch syndrome-associated colorectal cancer specimens (n = 30) to CD3-positive, CD8-positive and FOXP3-positive T cell infiltration in both tumor and normal mucosa. No significant correlation between B2M mutations and immune cell infiltration was observed in tumor tissue. However, FOXP3-positive T cell infiltration was significantly lower in normal mucosa adjacent to B2M-mutant (mt) compared to B2M-wild type (wt) tumors (mean: 0.98% FOXP3-positive area/region of interest (ROI) in B2M-wt vs. 0.52% FOXP3-positive area/ROI in B2M-mt, p = 0.023). Our results suggest that in the absence of immune-suppressive regulatory T cells (Treg), the outgrowth of less immunogenic B2M-mt tumor cells is favored. This finding supports the immunoediting concept in human solid cancer development and indicates a critical role of the immune milieu in normal colonic mucosa for the course of disease.

Molecular profiling of long-term survivors identifies a subgroup of glioblastoma characterized by chromosome 19/20 co-gain

Glioblastoma (GBM) is a devastating tumor and few patients survive beyond 3 years. Defining the molecular determinants underlying long-term survival is essential for insights into tumor biology and biomarker identification. We therefore investigated homogeneously treated, IDH (wt) long-term (LTS, n = 10) and short-term survivors (STS, n = 6) by microarray transcription profiling. While there was no association of clinical parameters and molecular subtypes with long-term survival, STS tumors were characterized by differential polarization of infiltrating microglia with predominance of the M2 phenotype detectable both on the mRNA and protein level. Furthermore, transcriptional signatures of LTS and STS predicted patient outcome in a large, IDH (wt) cohort (n = 468). Interrogation of overlapping genomic alterations identified concurrent gain of chromosomes 19 and 20 as a favorable prognostic marker. The strong association of this co-gain with survival was validated by aCGH in a second, independent cohort (n = 124). Finally, FISH and gene expression data revealed gains to constitute low-amplitude, clonal events with a strong impact on transcription. In conclusion, these findings provide important insights into the manipulation of the innate immune system by particularly aggressive GBM tumors. Furthermore, we genomically characterize a previously unknown, clinically relevant subgroup of glioblastoma, which can easily be identified through modern neuropathological workup.