The abundance and morphology of human large intestinal goblet and tuft cells during chronic schistosomiasis

Schistosomiasis affects nearly 240 million people in predominately low- and middle-income countries and ranks second in the number of cases and socio-economic burden among all parasitic diseases. Despite the enormous burden posed by schistosomes, our understanding of how schistosomiasis impacts infected human tissues remains limited. Intestinal schistosomiasis in animal models leads to goblet cell hyperplasia, likely increasing mucus production and reflecting an intestinal type 2 immune response. However, it is unknown whether these same changes occur in schistosome-infected humans. Using immunofluorescence and light microscopy, we compared the abundance and morphology of goblet cells in patients diagnosed with schistosomiasis to uninfected controls. The mucin-containing vesicles in goblet cells from schistosome-infected patients were significantly larger (hypertrophic) than uninfected individuals, although goblet cell hyperplasia was absent in chronic human schistosomiasis. In addition, we examined tuft cells in the large intestinal epithelium of control and schistosome-infected patients. Tuft cell numbers expand during helminth infection in mice, but these cells have not been characterized in human parasite infections. We found no evidence of tuft cell hyperplasia during human schistosome infection. Thus, our study provides novel insight into schistosome-associated changes to the intestinal epithelium in humans, suggesting an increase in mucus production by large intestinal goblet cells but relatively minor effects on tuft cell numbers.

DNA methylation-based classification of sinonasal tumors

The diagnosis of sinonasal tumors is challenging due to a heterogeneous spectrum of various differential diagnoses as well as poorly defined, disputed entities such as sinonasal undifferentiated carcinomas (SNUCs). In this study, we apply a machine learning algorithm based on DNA methylation patterns to classify sinonasal tumors with clinical-grade reliability. We further show that sinonasal tumors with SNUC morphology are not as undifferentiated as their current terminology suggests but rather reassigned to four distinct molecular classes defined by epigenetic, mutational and proteomic profiles. This includes two classes with neuroendocrine differentiation, characterized by IDH2 or SMARCA4/ARID1A mutations with an overall favorable clinical course, one class composed of highly aggressive SMARCB1-deficient carcinomas and another class with tumors that represent potentially previously misclassified adenoid cystic carcinomas. Our findings can aid in improving the diagnostic classification of sinonasal tumors and could help to change the current perception of SNUCs.

A real-time GPU-accelerated parallelized image processor for large-scale multiplexed fluorescence microscopy data

Highly multiplexed, single-cell imaging has revolutionized our understanding of spatial cellular interactions associated with health and disease. With ever-increasing numbers of antigens, region sizes, and sample sizes, multiplexed fluorescence imaging experiments routinely produce terabytes of data. Fast and accurate processing of these large-scale, high-dimensional imaging data is essential to ensure reliable segmentation and identification of cell types and for characterization of cellular neighborhoods and inference of mechanistic insights. Here, we describe RAPID, a Real-time, GPU-Accelerated Parallelized Image processing software for large-scale multiplexed fluorescence microscopy Data. RAPID deconvolves large-scale, high-dimensional fluorescence imaging data, stitches and registers images with axial and lateral drift correction, and minimizes tissue autofluorescence such as that introduced by erythrocytes. Incorporation of an open source CUDA-driven, GPU-assisted deconvolution produced results similar to fee-based commercial software. RAPID reduces data processing time and artifacts and improves image contrast and signal-to-noise compared to our previous image processing pipeline, thus providing a useful tool for accurate and robust analysis of large-scale, multiplexed, fluorescence imaging data.

The Tabula Sapiens: A multiple-organ, single-cell transcriptomic atlas of humans

Molecular characterization of cell types using single-cell transcriptome sequencing is revolutionizing cell biology and enabling new insights into the physiology of human organs. We created a human reference atlas comprising nearly 500,000 cells from 24 different tissues and organs, many from the same donor. This atlas enabled molecular characterization of more than 400 cell types, their distribution across tissues, and tissue-specific variation in gene expression. Using multiple tissues from a single donor enabled identification of the clonal distribution of T cells between tissues, identification of the tissue-specific mutation rate in B cells, and analysis of the cell cycle state and proliferative potential of shared cell types across tissues. Cell type-specific RNA splicing was discovered and analyzed across tissues within an individual.

MUC4 is expressed in alveolar rhabdomyosarcoma

AIMS

Mucin 4 (MUC4) is a transmembrane glycoprotein normally expressed by several human epithelial surfaces, including those of the colon, vagina, and respiratory tract. Although MUC4 overexpression is seen in various carcinomas, its expression among mesenchymal neoplasms is fairly specific to low-grade fibromyxoid sarcoma and sclerosing epithelioid fibrosarcoma. Having observed unanticipated anti-MUC4 immunoreactivity in rhabdomyosarcoma, we aimed to further characterise its expression.

METHODS AND RESULTS

Expression of MUC4 was assessed by immunohistochemistry in a total of 97 rhabdomyosarcomas using formalin-fixed paraffin-embedded tissue sections. MUC4 was expressed by 21 of 26 PAX3/7-FOXO1 fusion-positive cases, wherein immunoreactivity, varying from weak to strong, was present in 20-100% of neoplastic cells. With the exception of one sclerosing rhabdomyosarcoma showing immunoreactivity in 20% of cells, MUC4 was not expressed by embryonal (n = 28), sclerosing (n = 20), or pleomorphic (n = 23) rhabdomyosarcomas. Analysing published gene expression microarray data from a separate cohort of 33 fusion-positive and 25 fusion-negative rhabdomyosarcomas, we found on average a 11.4-fold increased expression in fusion-positive tumours (P = 0.0004).

CONCLUSIONS

MUC4 is expressed to a variable extent in the majority of PAX3/7-FOXO1 fusion-positive (alveolar) rhabdomyosarcomas, while expression in other rhabdomyosarcoma subtypes is rare.

Morphological, immunophenotypical and molecular features of hypermutation in colorectal carcinomas with mutations in DNA polymerase ε (POLE)

AIMS

Colorectal carcinomas (CRC) with mismatch repair (MMR) deficiency have increased tumour mutation burden and respond to immune check-point inhibitor therapy. The Cancer Genome Atlas identified hypermutated CRCs with somatic mutations in DNA polymerase ε (POLE) with mutation burdens exceeding that of MMR-deficient CRCs.

METHODS AND RESULTS

To identify the morphological, immunophenotypical and molecular features of POLE-mutated CRCs, 63 consecutive MMR-intact CRCs were evaluated by Sanger sequencing for POLE exonuclease domain mutations in exons 9, 11, 13 and 14 and confirmed by next-generation sequencing. Tumour immune microenvironment and IMMUNOSCORE®¹ were assessed in POLE-mutated CRCs using immunohistochemistry to detect CD3⁺ /CD8⁺ tumour-infiltrating lymphocytes and compared to 59 non-POLE mutated MMR-intact CRC, 10 non-POLE mutated MMR-deficient CRCs and 223 normal colonic mucosa.

CONCLUSIONS

A total of 4.8% CRC (four MMR-intact primary and one MMR-intact metastasis) harboured POLE mutations in amino acid 286 in exon 9 (p.P286R) or exon 13 (p.V411L). POLE-mutated CRCs arose in the transverse colon and rectum, were male-predominant, younger and showed increased tumour-infiltrating lymphocytes and immune cells at the tumour-stromal interface. The patient with metastatic POLE-mutated CRC was placed on PD-1 inhibitor treatment with marked and sustained response. These data indicate that POLE-mutated CRCs have hypermutated phenotypes despite MMR-intact status, with mutation burdens higher than that in microsatellite-unstable CRCs. Given the recent approval for treatment of microsatellite-unstable cancer with immune check-point inhibitors, assessment of POLE status may help to guide therapeutic decisions for hypermutated tumours with intact MMR that would otherwise be missed by routine testing.

Detection of Circulating Tumor DNA in Patients With Uterine Leiomyomas

PURPOSE

The preoperative distinction between uterine leiomyoma (LM) and leiomyosarcoma (LMS) is difficult, which may result in dissemination of an unexpected malignancy during surgery for a presumed benign lesion. An assay based on circulating tumor DNA (ctDNA) could help in the preoperative distinction between LM and LMS. This study addresses the feasibility of applying the two most frequently used approaches for detection of ctDNA: profiling of copy number alterations (CNAs) and point mutations in the plasma of patients with LM.

PATIENTS AND METHODS

By shallow whole-genome sequencing, we prospectively examined whether LM-derived ctDNA could be detected in plasma specimens of 12 patients. Plasma levels of lactate dehydrogenase, a marker suggested for the distinction between LM and LMS by prior studies, were also determined. We also profiled 36 LM tumor specimens by exome sequencing to develop a panel for targeted detection of point mutations in ctDNA of patients with LM.

RESULTS

We identified tumor-derived CNAs in the plasma DNA of 50% (six of 12) of patients with LM. The lactate dehydrogenase levels did not allow for an accurate distinction between patients with LM and patients with LMS. We identified only two recurrently mutated genes in LM tumors (MED12 and ACLY).

CONCLUSION

Our results show that LMs do shed DNA into the circulation, which provides an opportunity for the development of ctDNA-based testing to distinguish LM from LMS. Although we could not design an LM-specific panel for ctDNA profiling, we propose that the detection of CNAs or point mutations in selected tumor suppressor genes in ctDNA may favor a diagnosis of LMS, since these genes are not affected in LM.

Combination Approach for Detecting Different Types of Alterations in Circulating Tumor DNA in Leiomyosarcoma

Purpose: The clinical utility of circulating tumor DNA (ctDNA) monitoring has been shown in tumors that harbor highly recurrent mutations. Leiomyosarcoma represents a type of tumor with a wide spectrum of heterogeneous genomic abnormalities; thus, targeting hotspot mutations or a narrow genomic region for ctDNA detection may not be practical. Here, we demonstrate a combinatorial approach that integrates different sequencing protocols for the orthogonal detection of single-nucleotide variants (SNV), small indels, and copy-number alterations (CNA) in ctDNA.Experimental Design: We employed Cancer Personalized Profiling by deep Sequencing (CAPP-Seq) for the analysis of SNVs and indels, together with a genome-wide interrogation of CNAs by Genome Representation Profiling (GRP). We profiled 28 longitudinal plasma samples and 25 tumor specimens from 7 patients with leiomyosarcoma.Results: We detected ctDNA in 6 of 7 of these patients with >98% specificity for mutant allele fractions down to a level of 0.01%. We show that results from CAPP-Seq and GRP are highly concordant, and the combination of these methods allows for more comprehensive monitoring of ctDNA by profiling a wide spectrum of tumor-specific markers. By analyzing multiple tumor specimens in individual patients obtained from different sites and at different times during treatment, we observed clonal evolution of these tumors that was reflected by ctDNA profiles.Conclusions: Our strategy allows for the comprehensive monitoring of a broad spectrum of tumor-specific markers in plasma. Our approach may be clinically useful not only in leiomyosarcoma but also in other tumor types that lack recurrent genomic alterations. Clin Cancer Res; 24(11); 2688-99. ©2018 AACR.

Abstract A05: Circulating tumor DNA levels correlate with response to treatment in LMS patients

Circulating tumor DNA (ctDNA) has significant potential for several clinical applications, including assessment of treatment response and monitoring of recurrent/residual disease. We performed a pilot study to explore the feasibility of ctDNA monitoring in patients with leiomyosarcoma (LMS).

We profiled matching plasma and FFPE tumor specimens from 9 LMS patients. We analyzed between 2 to 6 longitudinal plasma samples (median of 5) and between 1 to 7 tumor specimens (median of 2) per patient. ctDNA analysis was performed on plasma samples collected pre-/post-surgery, throughout chemo-/radiotherapy and during follow-up. We used two separate approaches in our study: 1) targeted deep sequencing of ctDNA, tumor DNA and germline DNA to detect single nucleotide variants and indels using Cancer Personalized Profiling by deep Sequencing with integrated digital error suppression (CAPP-Seq; with a median deduplicated depth of sequencing of 2,136x); 2) copy number variant analysis in ctDNA by genome representation profiling (GRP; median coverage across the whole genome 0.23x) and in the matched tumors by SNP arrays. One patient was excluded from the analysis due to inadequate sequencing coverage in tumor specimen.

For CAPP-Seq analysis, we designed a custom 184kb capture panel targeting 89 genes that are recurrently mutated in LMS. Using strict variant calling criteria (requiring that variants be present on each strand of the original DNA “duplex” molecule) our panel identified a median of one nonsynonymous coding/splicing variant per tumor. We detected the same variants in TP53, RB1 and ATRX genes in ctDNA of 6/8 patients (with a baseline sensitivity of 87.5% and overall specificity of 98.96% calculated using plasma from 24 healthy donors). These six patients presented with advanced disease at the time of the first blood collection and were progressing throughout multiple lines of therapy. Two patients who did not have any variants detectable by CAPP-Seq in plasma had localized disease at the time of the first blood collection and/or responded well to the therapy. We found that changes in ctDNA levels appear to correspond with the extent of disease and response to treatment. Specifically, ctDNA levels decreased in a subset of patients after surgery or at the time of temporary response to chemo- and/or radiotherapy. Congruently, increases in ctDNA levels correlated with progression in most of the patients. There was a high correlation between ctDNA levels detected by CAPP-Seq (quantified as mutant molecules/mL plasma) and GRP (quantified as percent of genome showing copy number aberrations) across all plasma samples (Pearson's r= 0.88, p < 0.0001), but in a few samples ctDNA was detected by only one of the two assays.

Our results suggest that serial analysis of ctDNA is a promising approach for evaluation of treatment response in LMS patients. Validation of these findings in a prospective study on a larger group of patients will be required to determine the use of this approach in a clinical setting.

References:

CAPP-Seq: PMIDs 24705333, 27018799

GRP: PMIDs 25585704, 26687610

Citation Format: Joanna Przybyl, Jacob J. Chabon, Lien Spans, Kristen Ganjoo, Sujay Vennam, Aaron M. Newman, Erna Forgó, Sushama Varma, Shirley Zhu, Maria Debiec-Rychter, Ash Alizadeh, Maximilian Diehn, Matt van de Rijn. Circulating tumor DNA levels correlate with response to treatment in LMS patients [abstract]. In: Proceedings of the AACR Conference on Advances in Sarcomas: From Basic Science to Clinical Translation; May 16-19, 2017; Philadelphia, PA. Philadelphia (PA): AACR; Clin Cancer Res 2018;24(2_Suppl):Abstract nr A05.

CDX2 as a Prognostic Biomarker in Stage II and Stage III Colon Cancer

Background The identification of high-risk stage II colon cancers is key to the selection of patients who require adjuvant treatment after surgery. Microarray-based multigene-expression signatures derived from stem cells and progenitor cells hold promise, but they are difficult to use in clinical practice. Methods We used a new bioinformatics approach to search for biomarkers of colon epithelial differentiation across gene-expression arrays and then ranked candidate genes according to the availability of clinical-grade diagnostic assays. With the use of subgroup analysis involving independent and retrospective cohorts of patients with stage II or stage III colon cancer, the top candidate gene was tested for its association with disease-free survival and a benefit from adjuvant chemotherapy. Results The transcription factor CDX2 ranked first in our screening test. A group of 87 of 2115 tumor samples (4.1%) lacked CDX2 expression. In the discovery data set, which included 466 patients, the rate of 5-year disease-free survival was lower among the 32 patients (6.9%) with CDX2-negative colon cancers than among the 434 (93.1%) with CDX2-positive colon cancers (hazard ratio for disease recurrence, 3.44; 95% confidence interval [CI], 1.60 to 7.38; P=0.002). In the validation data set, which included 314 patients, the rate of 5-year disease-free survival was lower among the 38 patients (12.1%) with CDX2 protein-negative colon cancers than among the 276 (87.9%) with CDX2 protein-positive colon cancers (hazard ratio, 2.42; 95% CI, 1.36 to 4.29; P=0.003). In both these groups, these findings were independent of the patient's age, sex, and tumor stage and grade. Among patients with stage II cancer, the difference in 5-year disease-free survival was significant both in the discovery data set (49% among 15 patients with CDX2-negative tumors vs. 87% among 191 patients with CDX2-positive tumors, P=0.003) and in the validation data set (51% among 15 patients with CDX2-negative tumors vs. 80% among 106 patients with CDX2-positive tumors, P=0.004). In a pooled database of all patient cohorts, the rate of 5-year disease-free survival was higher among 23 patients with stage II CDX2-negative tumors who were treated with adjuvant chemotherapy than among 25 who were not treated with adjuvant chemotherapy (91% vs. 56%, P=0.006). Conclusions Lack of CDX2 expression identified a subgroup of patients with high-risk stage II colon cancer who appeared to benefit from adjuvant chemotherapy. (Funded by the National Comprehensive Cancer Network, the National Institutes of Health, and others.).

Molecular subtyping of leiomyosarcoma with 3' end RNA sequencing

Leiomyosarcoma (LMS) is a malignant neoplasm with smooth muscle differentiation. Little is known about its molecular heterogeneity and no targeted therapy currently exists for LMS. We performed expression profiling on 99 cases of LMS with 3'end RNA sequencing (3SEQ) and demonstrated the existence of 3 molecular subtypes in this cohort. We consequently showed that these molecular subtypes are reproducible using an independent cohort of 82 LMS cases from TCGA. Two new formalin-fixed, paraffin-embedded (FFPE) tissue-compatible diagnostic immunohistochemical markers were identified for two of the three subtypes: LMOD1 for subtype I LMS and ARL4C for subtype II LMS. Subtype I and subtype II LMS were associated with good and poor prognosis, respectively. Here, we describe the details of LMS diagnosis, RNA isolation, 3SEQ library construction, 3SEQ sequencing data analysis and molecular subtype determination. The 3SEQ data produced in this study was deposited into Gene Expression Omnibus (GEO) under GSE45510.

Chromosomal copy number alterations for associations of ductal carcinoma in situ with invasive breast cancer

INTRODUCTION

Screening mammography has contributed to a significant increase in the diagnosis of ductal carcinoma in situ (DCIS), raising concerns about overdiagnosis and overtreatment. Building on prior observations from lineage evolution analysis, we examined whether measuring genomic features of DCIS would predict association with invasive breast carcinoma (IBC). The long-term goal is to enhance standard clinicopathologic measures of low- versus high-risk DCIS and to enable risk-appropriate treatment.

METHODS

We studied three common chromosomal copy number alterations (CNA) in IBC and designed fluorescence in situ hybridization-based assay to measure copy number at these loci in DCIS samples. Clinicopathologic data were extracted from the electronic medical records of Stanford Cancer Institute and linked to demographic data from the population-based California Cancer Registry; results were integrated with data from tissue microarrays of specimens containing DCIS that did not develop IBC versus DCIS with concurrent IBC. Multivariable logistic regression analysis was performed to describe associations of CNAs with these two groups of DCIS.

RESULTS

We examined 271 patients with DCIS (120 that did not develop IBC and 151 with concurrent IBC) for the presence of 1q, 8q24 and 11q13 copy number gains. Compared to DCIS-only patients, patients with concurrent IBC had higher frequencies of CNAs in their DCIS samples. On multivariable analysis with conventional clinicopathologic features, the copy number gains were significantly associated with concurrent IBC. The state of two of the three copy number gains in DCIS was associated with a risk of IBC that was 9.07 times that of no copy number gains, and the presence of gains at all three genomic loci in DCIS was associated with a more than 17-fold risk (P = 0.0013).

CONCLUSIONS

CNAs have the potential to improve the identification of high-risk DCIS, defined by presence of concurrent IBC. Expanding and validating this approach in both additional cross-sectional and longitudinal cohorts may enable improved risk stratification and risk-appropriate treatment in DCIS.

Clinically Relevant Molecular Subtypes in Leiomyosarcoma

PURPOSE

Leiomyosarcoma is a malignant neoplasm with smooth muscle differentiation. Little is known about its molecular heterogeneity and no targeted therapy currently exists for leiomyosarcoma. Recognition of different molecular subtypes is necessary to evaluate novel therapeutic options. In a previous study on 51 leiomyosarcomas, we identified three molecular subtypes in leiomyosarcoma. The current study was performed to determine whether the existence of these subtypes could be confirmed in independent cohorts.

EXPERIMENTAL DESIGN

Ninety-nine cases of leiomyosarcoma were expression profiled with 3'end RNA-Sequencing (3SEQ). Consensus clustering was conducted to determine the optimal number of subtypes.

RESULTS

We identified 3 leiomyosarcoma molecular subtypes and confirmed this finding by analyzing publically available data on 82 leiomyosarcoma from The Cancer Genome Atlas (TCGA). We identified two new formalin-fixed, paraffin-embedded tissue-compatible diagnostic immunohistochemical markers; LMOD1 for subtype I leiomyosarcoma and ARL4C for subtype II leiomyosarcoma. A leiomyosarcoma tissue microarray with known clinical outcome was used to show that subtype I leiomyosarcoma is associated with good outcome in extrauterine leiomyosarcoma while subtype II leiomyosarcoma is associated with poor prognosis in both uterine and extrauterine leiomyosarcoma. The leiomyosarcoma subtypes showed significant differences in expression levels for genes for which novel targeted therapies are being developed, suggesting that leiomyosarcoma subtypes may respond differentially to these targeted therapies.

CONCLUSIONS

We confirm the existence of 3 molecular subtypes in leiomyosarcoma using two independent datasets and show that the different molecular subtypes are associated with distinct clinical outcomes. The findings offer an opportunity for treating leiomyosarcoma in a subtype-specific targeted approach.