Heterogeneity in intrahepatic macrophage populations and druggable target expression in patients with steatotic liver disease-related fibrosis

Background & Aims

Clinical trials for reducing fibrosis in steatotic liver disease (SLD) have targeted macrophages with variable results. We evaluated intrahepatic macrophages in patients with SLD to determine if activity scores or fibrosis stages influenced phenotypes and expression of druggable targets, such as CCR2 and galectin-3.

Methods

Liver biopsies from controls or patients with minimal or advanced fibrosis were subject to gene expression analysis using nCounter to determine differences in macrophage-related genes (n = 30). To investigate variability among individual patients, we compared additional biopsies by staining them with multiplex antibody panels (CD68/CD14/CD16/CD163/Mac387 or CD163/CCR2/galectin-3/Mac387) followed by spectral imaging and spatial analysis. Algorithms that utilize deep learning/artificial intelligence were applied to create cell cluster plots, phenotype profile maps, and to determine levels of protein expression (n = 34).

Results

Several genes known to be pro-fibrotic (e.g. CD206, TREM2, CD163, and ARG1) showed either no significant differences or significantly decreased with advanced fibrosis. Although marked variability in gene expression was observed in individual patients with cirrhosis, several druggable targets and their ligands (e.g. CCR2, CCR5, CCL2, CCL5, and LGALS3) were significantly increased when compared to patients with minimal fibrosis. Antibody panels identified populations that were significantly increased (e.g. Mac387+), decreased (e.g. CD14+), or enriched (e.g. interactions of Mac387) in patients that had progression of disease or advanced fibrosis. Despite heterogeneity in patients with SLD, several macrophage phenotypes and druggable targets showed a positive correlation with increasing NAFLD activity scores and fibrosis stages.

Conclusions

Patients with SLD have markedly varied macrophage- and druggable target-related gene and protein expression in their livers. Several patients had relatively high expression, while others were like controls. Overall, patients with more advanced disease had significantly higher expression of CCR2 and galectin-3 at both the gene and protein levels.

Impact and implications

Appreciating individual differences within the hepatic microenvironment of patients with SLD may be paramount to developing effective treatments. These results may explain why such a small percentage of patients have responded to macrophage-targeting therapies and provide additional support for precision medicine-guided treatment of chronic liver diseases.

Patients with fibrosis from non-alcoholic steatohepatitis have heterogeneous intrahepatic macrophages and therapeutic targets

Background and Aims

In clinical trials for reducing fibrosis in NASH patients, therapeutics that target macrophages have had variable results. We evaluated intrahepatic macrophages in patients with non-alcoholic steatohepatitis to determine if fibrosis influenced phenotypes and expression of CCR2 and Galectin-3.

Approach & Results

We used nCounter to analyze liver biopsies from well-matched patients with minimal (n=12) or advanced (n=12) fibrosis to determine which macrophage-related genes would be significantly different. Known therapy targets (e.g., CCR2 and Galectin-3) were significantly increased in patients with cirrhosis.However, several genes (e.g., CD68, CD16, and CD14) did not show significant differences, and CD163, a marker of pro-fibrotic macrophages was significantly decreased with cirrhosis. Next, we analyzed patients with minimal (n=6) or advanced fibrosis (n=5) using approaches that preserved hepatic architecture by multiplex-staining with anti-CD68, Mac387, CD163, CD14, and CD16. Spectral data were analyzed using deep learning/artificial intelligence to determine percentages and spatial relationships. This approach showed patients with advanced fibrosis had increased CD68+, CD16+, Mac387+, CD163+, and CD16+CD163+ populations. Interaction of CD68+ and Mac387+ populations was significantly increased in patients with cirrhosis and enrichment of these same phenotypes in individuals with minimal fibrosis correlated with poor outcomes. Evaluation of a final set of patients (n=4) also showed heterogenous expression of CD163, CCR2, Galectin-3, and Mac387, and significant differences were not dependent on fibrosis stage or NAFLD activity.

Conclusions

Approaches that leave hepatic architecture intact, like multispectral imaging, may be paramount to developing effective treatments for NASH. In addition, understanding individual differences in patients may be required for optimal responses to macrophage-targeting therapies.

Subclonal evolution and expansion of spatially distinct THY1-positive cells is associated with recurrence in glioblastoma

PURPOSE

Glioblastoma(GBM) is a lethal disease characterized by inevitable recurrence. Here we investigate the molecular pathways mediating resistance, with the goal of identifying novel therapeutic opportunities.

EXPERIMENTAL DESIGN

We developed a longitudinal in vivo recurrence model utilizing patient-derived explants to produce paired specimens(pre- and post-recurrence) following temozolomide(TMZ) and radiation(IR). These specimens were evaluated for treatment response and to identify gene expression pathways driving treatment resistance. Findings were clinically validated using spatial transcriptomics of human GBMs.

RESULTS

These studies reveal in replicate cohorts, a gene expression profile characterized by upregulation of mesenchymal and stem-like genes at recurrence. Analyses of clinical databases revealed significant association of this transcriptional profile with worse overall survival and upregulation at recurrence. Notably, gene expression analyses identified upregulation of TGFβ signaling, and more than one-hundred-fold increase in THY1 levels at recurrence. Furthermore, THY1-positive cells represented <10% of cells in treatment-naïve tumors, compared to 75-96% in recurrent tumors. We then isolated THY1-positive cells from treatment-naïve patient samples and determined that they were inherently resistant to chemoradiation in orthotopic models. Additionally, using image-guided biopsies from treatment-naïve human GBM, we conducted spatial transcriptomic analyses. This revealed rare THY1+ regions characterized by mesenchymal/stem-like gene expression, analogous to our recurrent mouse model, which co-localized with macrophages within the perivascular niche. We then inhibited TGFBRI activity in vivo which decreased mesenchymal/stem-like protein levels, including THY1, and restored sensitivity to TMZ/IR in recurrent tumors.

CONCLUSIONS

These findings reveal that GBM recurrence may result from tumor repopulation by pre-existing, therapy-resistant, THY1-positive, mesenchymal cells within the perivascular niche.

Abstract LB013: Subclonal evolution and expansion of spatially distinct THY1-positive cells is associated with recurrence in glioblastoma

Purpose: Glioblastoma (GBM) is a lethal disease characterized by inevitable recurrence. Here we investigate the molecular pathways mediating resistance, with the goal of identifying therapeutic opportunities to target this tumor.

Experimental Design: We developed a longitudinal in vivo recurrence model utilizing patient-derived explants to produce paired specimens (pre- and post-recurrence) following temozolomide(TMZ) and radiation(IR). These specimens were evaluated for treatment response and to identify gene expression pathways driving treatment resistance. Findings were clinically validated using spatial transcriptomics of human GBMs.

Results: These studies reveal in replicate cohorts, a gene expression profile characterized by upregulation of mesenchymal and stem-like genes at recurrence. Analyses of clinical databases revealed increased expression of this transcriptional profile to be significantly associated with worse median overall survival (248 days vs 430 days, p=0.0004), and upregulation of this profile at recurrence. Most notably, we identified upregulation of TGFβ signaling, and more than one-hundred-fold increase in THY1 levels at recurrence. Utilizing cell sorting, we observed that THY1-positive cells represented &lt;10% of cells in the treatment-naïve tumors and 75-96% in the recurrent tumors. We then isolated THY1-positive cells from treatment-naïve patient samples and determined that they were inherently resistant to chemoradiation in orthotopic models. Additionally, using image-guided biopsies from treatment-naïve human GBM, we conducted spatial transcriptomic analyses. This revealed rare THY1+ regions characterized by mesenchymal/stem-like gene expression, analogous to our recurrent mouse model samples, which co-localized with macrophages within the perivascular niche. Since TGFβ signaling contributes to a mesenchymal/stem-like phenotype, we inhibited TGFβRI activity in vivo which resulted in decreased mesenchymal/stem-like protein levels, including THY1, and restored sensitivity to TMZ/IR in recurrent tumors.

Conclusions: These findings reveal that GBM recurrence may result from tumor repopulation by pre-existing, therapy-resistant, THY1-positive, mesenchymal/stem-like cells within the perivascular niche. Furthermore, our data demonstrate the promise of targeting upregulated pathways in resistant subclones as a novel mechanism to achieve therapeutic response, and specifically that THY1 expression may represent a biomarker of response to TGFβ inhibition.

Citation Format: Wajd N. Al-Holou, Hanxiao Wang, Visweswaran Ravikumar, Morgan Oneka, Roel G. Verhaak, Hoon Kim, Drew Pratt, Sandra Camelo-Piragua, Corey Speers, Daniel Wahl, Sunita Shankar, Todd Hollon, Oren Sagher, Jason Heth, Karin M. Muraszko, Theodore S. Lawrence, Ana C. deCarvalho, Tom Mikkelsen, Arvind Rao, Alnawaz Rehemtulla. Subclonal evolution and expansion of spatially distinct THY1-positive cells is associated with recurrence in glioblastoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr LB013.

Stress response gene family expansions correlate with invasive potential in teleost fish

The bluegill sunfish Lepomis macrochirus and the closely related redear sunfish Lepomis microlophus have important ecological and recreational value and are widely used for research and aquaculture. While both species have been introduced outside of their native ranges, only the bluegill is considered invasive. Here, we report de novo transcriptome assemblies for these fish as a resource for sunfish biology. Comparative analyses of the transcriptomes revealed an unexpected, bluegill-specific expansion in the HSP70 and HSP90 molecular chaperone gene families. These expansions were not unique to the bluegill as expansions in HSP70s and HSP90s were identified in the genomes of other teleost fish using the NCBI RefSeq database. To determine whether gene family expansions are specific for thermal stress responses, GST and SOD gene families that are associated with oxidative stress responses were also analyzed. Species-specific expansions were also observed for these gene families in distinct fish species. Validating our approach, previously described expansions in the MHC gene family were also identified. Intriguingly, the number of HSP70 paralogs was positively correlated with thermotolerance range for each species, suggesting that these expansions can impact organismal physiology. Furthermore, fish that are considered invasive contained a higher average number of HSP70 paralogs than non-invasive fish. Invasive fish also had higher average numbers of HSP90, MHC and GST paralogs, but not SOD paralogs. Taken together, we propose that expansions in key cellular stress response gene families represent novel genetic signatures that correlate with invasive potential.

Developing and deploying an integrated workshop curriculum teaching computational skills for reproducible research

Inspired by well-established material and pedagogy provided by The Carpentries (Wilson, 2016), we developed a two-day workshop curriculum that teaches introductory R programming for managing, analyzing, plotting and reporting data using packages from the tidyverse (Wickham et al., 2019), the Unix shell, version control with git, and GitHub. While the official Software Carpentry curriculum is comprehensive, we found that it contains too much content for a two-day workshop. We also felt that the independent nature of the lessons left learners confused about how to integrate the newly acquired programming skills in their own work. Thus, we developed a new curriculum that aims to teach novices how to implement reproducible research principles in their own data analysis. The curriculum integrates live coding lessons with individual-level and group-based practice exercises, and also serves as a succinct resource that learners can reference both during and after the workshop. Moreover, it lowers the entry barrier for new instructors as they do not have to develop their own teaching materials or sift through extensive content. We developed this curriculum during a two-day sprint, successfully used it to host a two-day virtual workshop with almost 40 participants, and updated the material based on instructor and learner feedback. We hope that our new curriculum will prove useful to future instructors interested in teaching workshops with similar learning objectives.

CGAT: Cell Graph ATtention Network for Grading of Pancreatic Disease Histology Images

Early detection of Pancreatic Ductal Adenocarcinoma (PDAC), one of the most aggressive malignancies of the pancreas, is crucial to avoid metastatic spread to other body regions. Detection of pancreatic cancer is typically carried out by assessing the distribution and arrangement of tumor and immune cells in histology images. This is further complicated due to morphological similarities with chronic pancreatitis (CP), and the co-occurrence of precursor lesions in the same tissue. Most of the current automated methods for grading pancreatic cancers rely on extensive feature engineering involving accurate identification of cell features or utilising single number spatially informed indices for grading purposes. Moreover, sophisticated methods involving black-box approaches, such as neural networks, do not offer insights into the model's ability to accurately identify the correct disease grade. In this paper, we develop a novel cell-graph based Cell-Graph Attention (CGAT) network for the precise classification of pancreatic cancer and its precursors from multiplexed immunofluorescence histology images into the six different types of pancreatic diseases. The issue of class imbalance is addressed through bootstrapping multiple CGAT-nets, while the self-attention mechanism facilitates visualization of cell-cell features that are likely responsible for the predictive capabilities of the model. It is also shown that the model significantly outperforms the decision tree classifiers built using spatially informed metric, such as the Morisita-Horn (MH) indices.

GWAS of thyroid stimulating hormone highlights pleiotropic effects and inverse association with thyroid cancer

Thyroid stimulating hormone (TSH) is critical for normal development and metabolism. To better understand the genetic contribution to TSH levels, we conduct a GWAS meta-analysis at 22.4 million genetic markers in up to 119,715 individuals and identify 74 genome-wide significant loci for TSH, of which 28 are previously unreported. Functional experiments show that the thyroglobulin protein-altering variants P118L and G67S impact thyroglobulin secretion. Phenome-wide association analysis in the UK Biobank demonstrates the pleiotropic effects of TSH-associated variants and a polygenic score for higher TSH levels is associated with a reduced risk of thyroid cancer in the UK Biobank and three other independent studies. Two-sample Mendelian randomization using TSH index variants as instrumental variables suggests a protective effect of higher TSH levels (indicating lower thyroid function) on risk of thyroid cancer and goiter. Our findings highlight the pleiotropic effects of TSH-associated variants on thyroid function and growth of malignant and benign thyroid tumors.