Spatially Segregated Macrophage Populations Predict Distinct Outcomes In Colon Cancer.

Tumor-associated macrophages are transcriptionally heterogeneous, but the spatial distribution and cell interactions that shape macrophage tissue roles remain poorly characterized. Here, we spatially resolve five distinct human macrophage populations in normal and malignant human breast and colon tissue and reveal their cellular associations. This spatial map reveals that distinct macrophage populations reside in spatially segregated micro-environmental niches with conserved cellular compositions that are repeated across healthy and diseased tissue. We show that IL4I1+ macrophages phagocytose dying cells in areas with high cell turnover and predict good outcome in colon cancer. In contrast, SPP1+ macrophages are enriched in hypoxic and necrotic tumor regions and portend worse outcome in colon cancer. A subset of FOLR2+ macrophages is embedded in plasma cell niches. NLRP3+ macrophages co-localize with neutrophils and activate an inflammasome in tumors. Our findings indicate that a limited number of unique human macrophage niches function as fundamental building blocks in tissue.

Cost-effective DNA methylation profiling by FML-seq.

Current methods for profiling DNA methylation require costly reagents, sequencing, and labor time. We introduce fragmentation at methylated loci and sequencing (FML-seq), a sequencing library protocol that greatly reduces all these costs. Relative to other techniques tested on the same human cell lines, FML-seq produces similar measurements of absolute and differential cytosine methylation at a fraction of the price. FML-seq enables inexpensive, high-throughput experimental designs for large-scale epigenetics research projects.

Unraveling the Relevance of Tissue-Specific Decellularized Extracellular Matrix Hydrogels for Vocal Fold Regenerative Biomaterials: A Comprehensive Proteomic and In Vitro Study.

Decellularized extracellular matrix (dECM) is a promising material for tissue engineering applications. Tissue-specific dECM is often seen as a favorable material that recapitulates a native-like microenvironment for cellular remodeling. However, the minute quantity of dECM derivable from small organs like the vocal fold (VF) hampers manufacturing scalability. Small intestinal submucosa (SIS), a commercial product with proven regenerative capacity, may be a viable option for VF applications. This study aims to compare dECM hydrogels derived from SIS or VF tissue with respect to protein content and functionality using mass spectrometry-based proteomics and in vitro studies. Proteomic analysis reveals that VF and SIS dECM share 75% of core matrisome proteins. Although VF dECM proteins have greater overlap with native VF, SIS dECM shows less cross-sample variability. Following decellularization, significant reductions of soluble collagen (61%), elastin (81%), and hyaluronan (44%) are noted in VF dECM. SIS dECM contains comparable elastin and hyaluronan but 67% greater soluble collagen than VF dECM. Cells deposit more neo-collagen on SIS than VF-dECM hydrogels, whereas neo-elastin (~50 µg/scaffold) and neo-hyaluronan (~ 6 µg/scaffold) are comparable between the two hydrogels. Overall, SIS dECM possesses reasonably similar proteomic profile and regenerative capacity to VF dECM. SIS dECM is considered a promising alternative for dECM-derived biomaterials for VF regeneration.

A multiplexed in vivo approach to identify driver genes in small cell lung cancer.

Small cell lung cancer (SCLC) is a lethal form of lung cancer. Here, we develop a quantitative multiplexed approach on the basis of lentiviral barcoding with somatic CRISPR-Cas9-mediated genome editing to functionally investigate candidate regulators of tumor initiation and growth in genetically engineered mouse models of SCLC. We found that naphthalene pre-treatment enhances lentiviral vector-mediated SCLC initiation, enabling high multiplicity of tumor clones for analysis through high-throughput sequencing methods. Candidate drivers of SCLC identified from a meta-analysis across multiple human SCLC genomic datasets were tested using this approach, which defines both positive and detrimental impacts of inactivating 40 genes across candidate pathways on SCLC development. This analysis and subsequent validation in human SCLC cells establish TSC1 in the PI3K-AKT-mTOR pathway as a robust tumor suppressor in SCLC. This approach should illuminate drivers of SCLC, facilitate the development of precision therapies for defined SCLC genotypes, and identify therapeutic targets.

Cost-effective DNA methylation profiling by FML-seq.

Current methods for profiling DNA methylation require costly reagents, sequencing, or labor time. We introduce FML-seq, a sequencing library protocol that greatly reduces all these costs. Relative to other techniques tested on the same human cell lines, FML-seq produces similar measurements of absolute and differential cytosine methylation at a fraction of the price. FML-seq enables inexpensive, high-throughput experimental designs for large-scale epigenetics research projects.

A spatial map of human macrophage niches reveals context-dependent macrophage functions in colon and breast cancer.

Tumor-associated macrophages (TAMs) display heterogeneous phenotypes. Yet the exact tissue cues that shape macrophage functional diversity are incompletely understood. Here we discriminate, spatially resolve and reveal the function of five distinct macrophage niches within malignant and benign breast and colon tissue. We found that SPP1 TAMs reside in hypoxic and necrotic tumor regions, and a novel subset of FOLR2 tissue resident macrophages (TRMs) supports the plasma cell tissue niche. We discover that IL4I1 macrophages populate niches with high cell turnover where they phagocytose dying cells. Significantly, IL4I1 TAMs abundance correlates with anti-PD1 treatment response in breast cancer. Furthermore, NLRP3 inflammasome activation in NLRP3 TAMs correlates with neutrophil infiltration in the tumors and is associated with poor outcome in breast cancer patients. This suggests the NLRP3 inflammasome as a novel cancer immunetherapy target. Our work uncovers context-dependent roles of macrophage subsets, and suggests novel predictive markers and macrophage subset-specific therapy targets.

Molecular classification and biomarkers of clinical outcome in breast ductal carcinoma in situ: Analysis of TBCRC 038 and RAHBT cohorts.

Ductal carcinoma in situ (DCIS) is the most common precursor of invasive breast cancer (IBC), with variable propensity for progression. We perform multiscale, integrated molecular profiling of DCIS with clinical outcomes by analyzing 774 DCIS samples from 542 patients with 7.3 years median follow-up from the Translational Breast Cancer Research Consortium 038 study and the Resource of Archival Breast Tissue cohorts. We identify 812 genes associated with ipsilateral recurrence within 5 years from treatment and develop a classifier that predicts DCIS or IBC recurrence in both cohorts. Pathways associated with recurrence include proliferation, immune response, and metabolism. Distinct stromal expression patterns and immune cell compositions are identified. Our multiscale approach employed in situ methods to generate a spatially resolved atlas of breast precancers, where complementary modalities can be directly compared and correlated with conventional pathology findings, disease states, and clinical outcome.

The microdissected gene expression landscape of nasopharyngeal cancer reveals vulnerabilities in FGF and noncanonical NF-κB signaling.

Nasopharyngeal cancer (NPC) is an Epstein-Barr virus (EBV)-positive epithelial malignancy with an extensive inflammatory infiltrate. Traditional RNA-sequencing techniques uncovered only microenvironment signatures, while the gene expression of the tumor epithelial compartment has remained a mystery. Here, we use Smart-3SEQ to prepare transcriptome-wide gene expression profiles from microdissected NPC tumors, dysplasia, and normal controls. We describe changes in biological pathways across the normal to tumor spectrum and show that fibroblast growth factor (FGF) ligands are overexpressed in NPC tumors, while negative regulators of FGF signaling, including SPRY1, SPRY2, and LGALS3, are down-regulated early in carcinogenesis. Within the NF-κB signaling pathway, the critical noncanonical transcription factors, RELB and NFKB2, are enriched in the majority of NPC tumors. We confirm the responsiveness of EBV-positive NPC cell lines to targeted inhibition of these pathways, reflecting the heterogeneity in NPC patient tumors. Our data comprehensively describe the gene expression landscape of NPC and unravel the mysteries of receptor tyrosine kinase and NF-κB pathways in NPC.

Transcriptome and genome evolution during HER2-amplified breast neoplasia.

BACKGROUND: The acquisition of oncogenic drivers is a critical feature of cancer progression. For some carcinomas, it is clear that certain genetic drivers occur early in neoplasia and others late. Why these drivers are selected and how these changes alter the neoplasia's fitness is less understood. METHODS: Here we use spatially oriented genomic approaches to identify transcriptomic and genetic changes at the single-duct level within precursor neoplasia associated with invasive breast cancer. We study HER2 amplification in ductal carcinoma in situ (DCIS) as an event that can be both quantified and spatially located via fluorescence in situ hybridization (FISH) and immunohistochemistry on fixed paraffin-embedded tissue. RESULTS: By combining the HER2-FISH with the laser capture microdissection (LCM) Smart-3SEQ method, we found that HER2 amplification in DCIS alters the transcriptomic profiles and increases diversity of copy number variations (CNVs). Particularly, interferon signaling pathway is activated by HER2 amplification in DCIS, which may provide a prolonged interferon signaling activation in HER2-positive breast cancer. Multiple subclones of HER2-amplified DCIS with distinct CNV profiles are observed, suggesting that multiple events occurred for the acquisition of HER2 amplification. Notably, DCIS acquires key transcriptomic changes and CNV events prior to HER2 amplification, suggesting that pre-amplified DCIS may create a cellular state primed to gain HER2 amplification for growth advantage. CONCLUSION: By using genomic methods that are spatially oriented, this study identifies several features that appear to generate insights into neoplastic progression in precancer lesions at a single-duct level.

Gene Expression Profiling of Head and Neck Tumors Identifies FOXP1 and SOX10 Expression as Useful for Distinguishing Ameloblastoma From Basaloid Salivary Gland Tumors.

Odontogenic tumors show considerable morphologic heterogeneity and at times the diagnosis can be challenging. Ameloblastoma, the most common odontogenic tumor, can have morphologic similarity to some salivary gland tumors and therefore we sought to identify biomarkers that might aid in the diagnosis by performing transcriptome wide gene expression profiling of 80 odontogenic and salivary gland neoplasms. These data identified the FOXP1/SOX10 expression profile as characteristic of many odontogenic tumors including ameloblastoma but largely absent in salivary gland tumors. We then assessed 173 salivary gland tumors and 108 odontogenic tumors by immunohistochemistry for FOXP1 and SOX10 expression and found that 34/35 (97%) cases of ameloblastomas were diffusely positive for FOXP1 but completely negative for SOX10. None of the basaloid salivary neoplasms (basal cell adenoma, adenoid cystic carcinoma, polymorphous adenocarcinoma, and myoepitheloma) demonstrated FOXP1/SOX10 expression pattern. Taken together, the results of this study suggest that the FOXP1/SOX10 immunophenotype is common in odontogenic tumors including ameloblastoma and might be useful distinguishing these from similar appearing basaloid salivary gland tumors.

Gene expression profiling of single cells from archival tissue with laser-capture microdissection and Smart-3SEQ.

RNA sequencing (RNA-seq) is a sensitive and accurate method for quantifying gene expression. Small samples or those whose RNA is degraded, such as formalin-fixed paraffin-embedded (FFPE) tissue, remain challenging to study with nonspecialized RNA-seq protocols. Here, we present a new method, Smart-3SEQ, that accurately quantifies transcript abundance even with small amounts of total RNA and effectively characterizes small samples extracted by laser-capture microdissection (LCM) from FFPE tissue. We also obtain distinct biological profiles from FFPE single cells, which have been impossible to study with previous RNA-seq protocols, and we use these data to identify possible new macrophage phenotypes associated with the tumor microenvironment. We propose Smart-3SEQ as a highly cost-effective method to enable large gene expression profiling experiments unconstrained by sample size and tissue availability. In particular, Smart-3SEQ's compatibility with FFPE tissue unlocks an enormous number of archived clinical samples; combined with LCM it allows unprecedented studies of small cell populations and single cells isolated by their in situ context.

Genomic analysis of benign prostatic hyperplasia implicates cellular re-landscaping in disease pathogenesis.

Benign prostatic hyperplasia (BPH) is the most common cause of lower urinary tract symptoms in men. Current treatments target prostate physiology rather than BPH pathophysiology and are only partially effective. Here, we applied next-generation sequencing to gain new insight into BPH. By RNAseq, we uncovered transcriptional heterogeneity among BPH cases, where a 65-gene BPH stromal signature correlated with symptom severity. Stromal signaling molecules BMP5 and CXCL13 were enriched in BPH while estrogen regulated pathways were depleted. Notably, BMP5 addition to cultured prostatic myofibroblasts altered their expression profile towards a BPH profile that included the BPH stromal signature. RNAseq also suggested an altered cellular milieu in BPH, which we verified by immunohistochemistry and single-cell RNAseq. In particular, BPH tissues exhibited enrichment of myofibroblast subsets, whilst depletion of neuroendocrine cells and an estrogen receptor (ESR1)-positive fibroblast cell type residing near epithelium. By whole-exome sequencing, we uncovered somatic single-nucleotide variants (SNVs) in BPH, of uncertain pathogenic significance but indicative of clonal cell expansions. Thus, genomic characterization of BPH has identified a clinically-relevant stromal signature and new candidate disease pathways (including a likely role for BMP5 signaling), and reveals BPH to be not merely a hyperplasia, but rather a fundamental re-landscaping of cell types.

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.

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.

Cell-lineage heterogeneity and driver mutation recurrence in pre-invasive breast neoplasia.

BACKGROUND: All cells in an individual are related to one another by a bifurcating lineage tree, in which each node is an ancestral cell that divided into two, each branch connects two nodes, and the root is the zygote. When a somatic mutation occurs in an ancestral cell, all its descendants carry the mutation, which can then serve as a lineage marker for the phylogenetic reconstruction of tumor progression. Using this concept, we investigate cell lineage relationships and genetic heterogeneity of pre-invasive neoplasias compared to invasive carcinomas. METHODS: We deeply sequenced over a thousand phylogenetically informative somatic variants in 66 morphologically independent samples from six patients that represent a spectrum of normal, early neoplasia, carcinoma in situ, and invasive carcinoma. For each patient, we obtained a highly resolved lineage tree that establishes the phylogenetic relationships among the pre-invasive lesions and with the invasive carcinoma. RESULTS: The trees reveal lineage heterogeneity of pre-invasive lesions, both within the same lesion, and between histologically similar ones. On the basis of the lineage trees, we identified a large number of independent recurrences of PIK3CA H1047 mutations in separate lesions in four of the six patients, often separate from the diagnostic carcinoma. CONCLUSIONS: Our analyses demonstrate that multi-sample phylogenetic inference provides insights on the origin of driver mutations, lineage heterogeneity of neoplastic proliferations, and the relationship of genomically aberrant neoplasias with the primary tumors. PIK3CA driver mutations may be comparatively benign inducers of cellular proliferation.

A shared transcriptional program in early breast neoplasias despite genetic and clinical distinctions.

BACKGROUND: The earliest recognizable stages of breast neoplasia are lesions that represent a heterogeneous collection of epithelial proliferations currently classified based on morphology. Their role in the development of breast cancer is not well understood but insight into the critical events at this early stage will improve efforts in breast cancer detection and prevention. These microscopic lesions are technically difficult to study so very little is known about their molecular alterations. RESULTS: To characterize the transcriptional changes of early breast neoplasia, we sequenced 3'- end enriched RNAseq libraries from formalin-fixed paraffin-embedded tissue of early neoplasia samples and matched normal breast and carcinoma samples from 25 patients. We find that gene expression patterns within early neoplasias are distinct from both normal and breast cancer patterns and identify a pattern of pro-oncogenic changes, including elevated transcription of ERBB2, FOXA1, and GATA3 at this early stage. We validate these findings on a second independent gene expression profile data set generated by whole transcriptome sequencing. Measurements of protein expression by immunohistochemistry on an independent set of early neoplasias confirms that ER pathway regulators FOXA1 and GATA3, as well as ER itself, are consistently upregulated at this early stage. The early neoplasia samples also demonstrate coordinated changes in long non-coding RNA expression and microenvironment stromal gene expression patterns. CONCLUSIONS: This study is the first examination of global gene expression in early breast neoplasia, and the genes identified here represent candidate participants in the earliest molecular events in the development of breast cancer.

Carbon source utilization and inhibitor tolerance of 45 oleaginous yeast species.

Conversion of lignocellulosic hydrolysates to lipids using oleaginous (high lipid) yeasts requires alignment of the hydrolysate composition with the characteristics of the yeast strain, including ability to utilize certain nutrients, ability to grow independently of costly nutrients such as vitamins, and ability to tolerate inhibitors. Some combination of these characteristics may be present in wild strains. In this study, 48 oleaginous yeast strains belonging to 45 species were tested for ability to utilize carbon sources associated with lignocellulosic hydrolysates, tolerate inhibitors, and grow in medium without supplemented vitamins. Some well-studied oleaginous yeast species, as well as some that have not been frequently utilized in research or industrial production, emerged as promising candidates for industrial use due to ability to utilize many carbon sources, including Cryptococcus aureus, Cryptococcus laurentii, Hannaella aff. zeae, Tremella encephala, and Trichosporon coremiiforme. Other species excelled in inhibitor tolerance, including Candida aff. tropicalis, Cyberlindnera jadinii, Metschnikowia pulcherrima, Schwanniomyces occidentalis and Wickerhamomyces ciferrii. No yeast tested could utilize all carbon sources and tolerate all inhibitors tested. These results indicate that yeast strains should be selected based on characteristics compatible with the composition of the targeted hydrolysate. Other factors to consider include the production of valuable co-products such as carotenoids, availability of genetic tools, biosafety level, and flocculation of the yeast strain. The data generated in this study will aid in aligning yeasts with compatible hydrolysates for conversion of carbohydrates to lipids to be used for biofuels and other oleochemicals.

Improved glucose metabolism in vitro and in vivo by an allosteric monoclonal antibody that increases insulin receptor binding affinity.

Previously we reported studies of XMetA, an agonist antibody to the insulin receptor (INSR). We have now utilized phage display to identify XMetS, a novel monoclonal antibody to the INSR. Biophysical studies demonstrated that XMetS bound to the human and mouse INSR with picomolar affinity. Unlike monoclonal antibody XMetA, XMetS alone had little or no agonist effect on the INSR. However, XMetS was a strong positive allosteric modulator of the INSR that increased the binding affinity for insulin nearly 20-fold. XMetS potentiated insulin-stimulated INSR signaling ∼15-fold or greater including; autophosphorylation of the INSR, phosphorylation of Akt, a major enzyme in the metabolic pathway, and phosphorylation of Erk, a major enzyme in the growth pathway. The enhanced signaling effects of XMetS were more pronounced with Akt than with Erk. In cultured cells, XMetS also enhanced insulin-stimulated glucose transport. In contrast to its effects on the INSR, XMetS did not potentiate IGF-1 activation of the IGF-1 receptor. We studied the effect of XMetS treatment in two mouse models of insulin resistance and diabetes. The first was the diet induced obesity mouse, a hyperinsulinemic, insulin resistant animal, and the second was the multi-low dose streptozotocin/high-fat diet mouse, an insulinopenic, insulin resistant animal. In both models, XMetS normalized fasting blood glucose levels and glucose tolerance. In concert with its ability to potentiate insulin action at the INSR, XMetS reduced insulin and C-peptide levels in both mouse models. XMetS improved the response to exogenous insulin without causing hypoglycemia. These data indicate that an allosteric monoclonal antibody can be generated that markedly enhances the binding affinity of insulin to the INSR. These data also suggest that an INSR monoclonal antibody with these characteristics may have the potential to both improve glucose metabolism in insulinopenic type 2 diabetes mellitus and correct compensatory hyperinsulinism in insulin resistant conditions.