62 Applying machine vision to empower preclinical development of cell engager and adoptive cell therapeutics in patient-derived organoid models of solid tumors

Background

Cell engager and adoptive cell therapeutics have emerged as efficacious and durable treatments in patients with B-cell malignancies. Though many analogous strategies are under development in solid tumors, none have received approval. Preclinical development of these therapies requires cell labeling of immortalized cell lines and/or primary expanded T cells to distinguish target and effector cells. However, cell engager and adoptive cell therapies have had limited evidence of reproducibility in primary patient-derived models such as tumor organoid cultures thus far. Here, we build upon our tumor organoid platform1 to measure organoid specific responses to these therapies. Utilizing machine vision coupled with time-lapse-microscopy, we obtain multiparameter kinetic readouts of patient-derived tumor organoid cell killing and allogeneic MHC-matched primary peripheral blood mononuclear cells (PBMCs).

Methods

The patient-derived tumor organoids were co-cultured with PBMCs in the presence of engagers/activators and vital dyes and incubated for 96 hrs. Cell death was measured by quantifying the caspase 3/7 vital dye pixel intensities at different time points using high throughput imaging. As a first step, a fully convolutional neural network was trained to segment out organoids from brightfield images comprised of organoids, immune cells and potential background artifacts. This segmentation mask was then transferred over to registered caspase 3/7 images to quantify tumor cell specific phenotypes in a rapid and automated manner.

Results

The time-lapse imaging assay allowed for both the tracking of the organoid growth over time as well as the quantification of the kinetics of engagers/activators in comparison to controls resulting in accurate and precise technical reproducibility. Further, this assay allowed for the co-localization of the organoids and the immune cells over time, thus, enabling a spatiotemporal summary of dose dependent efficacy of candidate therapeutics.

Conclusions

We demonstrate the scalability and throughput of a machine vision tumor organoid immune co-culture platform across multiple unique patient-derived tumor organoid lines bearing a target of interest, enabling future discovery of biomarkers of therapeutic response and resistance.

Reference

Larsen B, Kannan M, Langer LF, Khan AA, Salahudeen AA, A pan-cancer organoid platform for precision medicine. Cell Reports 2021; 36:109429

Clinical and Genome-wide Analysis of Cisplatin-induced Tinnitus Implicates Novel Ototoxic Mechanisms

PURPOSE

Cisplatin, a commonly used chemotherapeutic, results in tinnitus, the phantom perception of sound. Our purpose was to identify the clinical and genetic determinants of tinnitus among testicular cancer survivors (TCS) following cisplatin-based chemotherapy.

EXPERIMENTAL DESIGN

TCS (n = 762) were dichotomized to cases (moderate/severe tinnitus; n = 154) and controls (none; n = 608). Logistic regression was used to evaluate associations with comorbidities and SNP dosages in genome-wide association study (GWAS) following quality control and imputation (covariates: age, noise exposure, cisplatin dose, genetic principal components). Pathway over-representation tests and functional studies in mouse auditory cells were performed.

RESULTS

Cisplatin-induced tinnitus (CisIT) significantly associated with age at diagnosis (P = 0.007) and cumulative cisplatin dose (P = 0.007). CisIT prevalence was not significantly greater in 400 mg/m²-treated TCS compared with 300 (P = 0.41), but doses >400 mg/m² (median 580, range 402-828) increased risk by 2.61-fold (P < 0.0001). CisIT cases had worse hearing at each frequency (0.25-12 kHz, P < 0.0001), and reported more vertigo (OR = 6.47; P < 0.0001) and problems hearing in a crowd (OR = 8.22; P < 0.0001) than controls. Cases reported poorer health (P < 0.0001) and greater psychotropic medication use (OR = 2.4; P = 0.003). GWAS suggested a variant near OTOS (rs7606353, P = 2 × 10^-6) and OTOS eQTLs were significantly enriched independently of that SNP (P = 0.018). OTOS overexpression in HEI-OC1, a mouse auditory cell line, resulted in resistance to cisplatin-induced cytotoxicity. Pathway analysis implicated potassium ion transport (q = 0.007).

CONCLUSIONS

CisIT associated with several neuro-otological symptoms, increased use of psychotropic medication, and poorer health. OTOS, expressed in the cochlear lateral wall, was implicated as protective. Future studies should investigate otoprotective targets in supporting cochlear cells.

Genome-Wide Association Studies of Chemotherapeutic Toxicities: Genomics of Inequality

With an estimated global population of cancer survivors exceeding 32 million and growing, there is a heightened awareness of the long-term toxicities resulting from cancer treatments and their impact on quality of life. Unexplained heterogeneity in the persistence and development of toxicities, as well as an incomplete understanding of their mechanisms, have generated a growing need for the identification of predictive pharmacogenomic markers. Early studies addressing this need used a candidate gene approach; however, over the last decade, unbiased and comprehensive genome-wide association studies (GWAS) have provided markers of phenotypic risk and potential targets to explore the mechanistic and regulatory pathways of biological functions associated with chemotherapeutic toxicity. In this review, we provide the current status of GWAS of chemotherapeutic toxicities with an emphasis on examining the ancestral diversity of the representative cohorts within these studies. Persistent calls to incorporate both ancestrally diverse and/or admixed populations into genomic efforts resulted in a recent rise in the number of studies utilizing cohorts of East Asian descent; however, few pharmacogenomic studies to date include cohorts of African, Indigenous American, Southwest Asian, and admixed populations. Through comprehensively evaluating sample size, composition by ancestry, genome-wide significant variants, and population-specific minor allele frequencies as reported by HapMap/dbSNP using NCBI PubMed and the NHGRI-EBI GWAS Catalog, we illustrate how allele frequencies and effect sizes tend to vary among individuals of differing ancestries. In an era of personalized medicine, the lack of diversity in genome-wide studies of anticancer agent toxicity may contribute to the health disparity gap. Clin Cancer Res; 23(15); 4010-9. ©2017 AACR.

Colonic transcriptional response to 1α,25(OH)2 vitamin D3 in African- and European-Americans

Colorectal cancer (CRC) is a significant health burden especially among African Americans (AA). Epidemiological studies have correlated low serum vitamin D with CRC risk, and, while hypovitaminosis D is more common and more severe in AA, the mechanisms by which vitamin D modulates CRC risk and how these differ by race are not well understood. Active vitamin D (1α,25(OH)₂D₃) has chemoprotective effects primarily through transcriptional regulation of target genes in the colon. We hypothesized that transcriptional response to 1α,25(OH)₂D₃ differs between AA and European Americans (EA) irrespective of serum vitamin D and that regulatory variants could impact transcriptional response. We treated ex vivo colon cultures from 34 healthy subjects (16 AA and 18 EA) with 0.1μM 1α,25(OH)₂D₃ or vehicle control for 6h and performed genome-wide transcriptional profiling. We found 8 genes with significant differences in transcriptional response to 1α,25(OH)₂D₃ between AA and EA with definitive replication of inter-ethnic differences for uridine phosphorylase 1 (UPP1) and zinc finger-SWIM containing 4 (ZSWIM4). We performed expression quantitative trait loci (eQTL) mapping and identified response cis-eQTLs for ZSWIM4 as well as for histone deacetylase 3 (HDAC3), the latter of which showed a trend toward significant inter-ethnic differences in transcriptional response. Allele frequency differences of eQTLs for ZSWIM4 and HDAC3 accounted for observed transcriptional differences between populations. Taken together, our results demonstrate that transcriptional response to 1α,25(OH)₂D₃ differs between AA and EA independent of serum 25(OH)D levels. We provide evidence in support of a genetic regulatory mechanism underlying transcriptional differences between populations for ZSWIM4 and HDAC3. Further work is needed to elucidate how response eQTLs modify vitamin D response and whether genotype and/or transcriptional response correlate with chemopreventive effects. Relevant biomarkers, such as tissue-specific 1α,25(OH)₂D₃ transcriptional response, could identify individuals likely to benefit from vitamin D for CRC prevention as well as elucidate basic mechanisms underlying CRC disparities.

Variants in WFS1 and Other Mendelian Deafness Genes Are Associated with Cisplatin-Associated Ototoxicity

Purpose: Cisplatin is one of the most commonly used chemotherapy drugs worldwide and one of the most ototoxic. We sought to identify genetic variants that modulate cisplatin-associated ototoxicity (CAO).Experimental Design: We performed a genome-wide association study (GWAS) of CAO using quantitative audiometry (4-12 kHz) in 511 testicular cancer survivors of European genetic ancestry. We performed polygenic modeling and functional analyses using a variety of publicly available databases. We used an electronic health record cohort to replicate our top mechanistic finding.Results: One SNP, rs62283056, in the first intron of Mendelian deafness gene WFS1 (wolframin ER transmembrane glycoprotein) and an expression quantitative trait locus (eQTL) for WFS1 met genome-wide significance for association with CAO (P = 1.4 × 10^-8). A significant interaction between cumulative cisplatin dose and rs62283056 genotype was evident, indicating that higher cisplatin doses exacerbate hearing loss in patients with the minor allele (P = 0.035). The association between decreased WFS1 expression and hearing loss was replicated in an independent BioVU cohort (n = 18,620 patients, Bonferroni adjusted P < 0.05). Beyond this top signal, we show CAO is a polygenic trait and that SNPs in and near 84 known Mendelian deafness genes are significantly enriched for low P values in the GWAS (P = 0.048).Conclusions: We show for the first time the role of WFS1 in CAO and document a statistically significant interaction between increasing cumulative cisplatin dose and rs62283056 genotype. Our clinical translational results demonstrate that pretherapy patient genotyping to minimize ototoxicity could be useful when deciding between cisplatin-based chemotherapy regimens of comparable efficacy with different cumulative doses. Clin Cancer Res; 23(13); 3325-33. ©2016 AACR.

Abstract B19: Colonic transcriptional response to 1,25(OH)2 vitamin D in African and European Americans

Background: 1,25(OH)2 vitamin D (1,25vitD) is a steroid hormone with anti-cancer effects that functions through direct transcriptional mechanisms. In colorectal cancer (CRC), 1,25vitD is thought to play a protective role; however, it is not known to what extent inter–ethnic differences in transcriptional response could underlie differences in disease susceptibility in African-Americans (AA) whose incidence of CRC is the highest in the US. To study differences in short-term transcriptional response in the colon, we have optimized an ex vivo culture system using colonic biopsies obtained during colonoscopy, thereby eliminating a number of confounders. We use this system to study transcriptional responses between ethnicities.

Methods: Healthy subjects (n=69; 34 in discovery cohort & 35 in replication cohort) undergoing colonoscopy were consented and colonic biopsies were obtained at the rectosigmoid junction. Two colonic biopsies each were treated with 0.1μM 1,25vitD or vehicle control (EtOH) and incubated for 6 hours at 37°C. For the discovery cohort, RNA from two treatment replicates was pooled and hybridized to Illumina Human HT-12 expression chips. Transcriptional response was analyzed using the packages LUMI and LIMMA in R adjusting for covariates and multiple testing. For the replication cohort, expression for select genes was performed using quantitative real-time PCR. In silico analysis was done to identify vitamin D receptor (VDR) binding peaks using ChIP-sequence data from LS180 cell lines.

Results: In the discovery cohort, 16 AA and 18 European Americans (EA) were included. As proof of concept, the top two significantly up-regulated genes were CD14 (p=1.15x10-25) and CYP24A1 (p=7.65x10-23). There were 10 genes that showed significantly different responses to 1,25vitD between AA and EA including UPP1, UCKL1, ZSWIM4, ERCC1, MFSDA2, EPHA2, CLRN3, KIAA1324, TRIP8 and WDR55. In the replication cohort, 20 AA and 15 EA were included. Of 5 genes tested to date, 3 genes (UPP1, ERCC1 & ZSWIM4) showed evidence of replication with similar direction and effect sizes as the discovery cohort. Two genes (UCKL1 & UPP1) have VDR binding peaks within an intron (UCKL1) or 21kb downstream (UPP1). Within these peaks, there are SNPs with allele frequency differences between Africans and Europeans. In particular, the UCKL1 SNP rs58190754 that is located within the DR3 binding motif and is predicted to be functionally significant for VDR binding has an allele frequency of 10% in YRI & AA and &lt;1% in Europeans.

Conclusions: Using short-term ex vivo culture of primary colonic tissue, we found significantly up- and down-regulated genes in response to 1,25vitD including a number of established 1,25vitD targets. There were 10 genes with significantly different transcriptional responses to 1,25vitD treatment between AA and EA of which 3 genes have been replicated in preliminary work. Two genes (UPP1 & UCKL1) harbor VDR binding peaks with candidate SNPs that could be functionally relevant for gene expression differences between populations. These genes are involved in pyrimidine salvage and could have a role in cancer susceptibility. Future work will genotype SNPs and determine whether genotype predicts expression. In addition, functional assays of pyrimidine salvage are being tested. Finally, replication of additional genes is ongoing. Further studies of these genes may reveal important pathways underlying CRC disparities.

Citation Format: Sonia Kupfer, Brandon Mapes, Shigeki Nakagome, David Witonsky, Anna Di Rienzo. Colonic transcriptional response to 1,25(OH)2 vitamin D in African and European Americans. [abstract]. In: Proceedings of the Eighth AACR Conference on The Science of Health Disparities in Racial/Ethnic Minorities and the Medically Underserved; Nov 13-16, 2015; Atlanta, GA. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2016;25(3 Suppl):Abstract nr B19.

Structure-Function Analysis of the Non-Muscle Myosin Light Chain Kinase (nmMLCK) Isoform by NMR Spectroscopy and Molecular Modeling: Influence of MYLK Variants

The MYLK gene encodes the multifunctional enzyme, myosin light chain kinase (MLCK), involved in isoform-specific non-muscle and smooth muscle contraction and regulation of vascular permeability during inflammation. Three MYLK SNPs (P21H, S147P, V261A) alter the N-terminal amino acid sequence of the non-muscle isoform of MLCK (nmMLCK) and are highly associated with susceptibility to acute lung injury (ALI) and asthma, especially in individuals of African descent. To understand the functional effects of SNP associations, we examined the N-terminal segments of nmMLCK by ¹H-¹⁵N heteronuclear single quantum correlation (HSQC) spectroscopy, a 2-D NMR technique, and by in silico molecular modeling. Both NMR analysis and molecular modeling indicated SNP localization to loops that connect the immunoglobulin-like domains of nmMLCK, consistent with minimal structural changes evoked by these SNPs. Molecular modeling analysis identified protein-protein interaction motifs adversely affected by these MYLK SNPs including binding by the scaffold protein 14-3-3, results confirmed by immunoprecipitation and western blot studies. These structure-function studies suggest novel mechanisms for nmMLCK regulation, which may confirm MYLK as a candidate gene in inflammatory lung disease and advance knowledge of the genetic underpinning of lung-related health disparities.

The NAMPT promoter is regulated by mechanical stress, signal transducer and activator of transcription 5, and acute respiratory distress syndrome-associated genetic variants

Increased nicotinamide phosphoribosyltransferase (NAMPT) transcription is mechanistically linked to ventilator-induced inflammatory lung injury (VILI), with VILI severity attenuated by reduced NAMPT bioavailability. The molecular mechanisms of NAMPT promoter regulation in response to excessive mechanical stress remain poorly understood. The objective of this study was to define the contribution of specific transcription factors, acute respiratory distress syndrome (ARDS)-associated single nucleotide polymorphisms (SNPs), and promoter demethylation to NAMPT transcriptional regulation in response to mechanical stress. In vivo NAMPT protein expression levels were examined in mice exposed to high tidal volume mechanical ventilation. In vitro NAMPT expression levels were examined in human pulmonary artery endothelial cells exposed to 5 or 18% cyclic stretch (CS), with NAMPT promoter activity assessed using NAMPT promoter luciferase reporter constructs with a series of nested deletions. In vitro NAMPT transcriptional regulation was further characterized by measuring luciferase activity, DNA demethylation, and chromatin immunoprecipitation. VILI-challenged mice exhibited significantly increased NAMPT expression in bronchoalveolar lavage leukocytes and in lung endothelium. A mechanical stress-inducible region (MSIR) was identified in the NAMPT promoter from -2,428 to -2,128 bp. This MSIR regulates NAMPT promoter activity, mRNA expression, and signal transducer and activator of transcription 5 (STAT5) binding, which is significantly increased by 18% CS. In addition, NAMPT promoter activity was increased by pharmacologic promoter demethylation and inhibited by STAT5 silencing. ARDS-associated NAMPT promoter SNPs rs59744560 (-948G/T) and rs7789066 (-2,422A/G) each significantly elevated NAMPT promoter activity in response to 18% CS in a STAT5-dependent manner. Our results show that NAMPT is a key novel ARDS therapeutic target and candidate gene with genetic/epigenetic transcriptional regulation in response to excessive mechanical stress.

Mechanical stress induces pre-B-cell colony-enhancing factor/NAMPT expression via epigenetic regulation by miR-374a and miR-568 in human lung endothelium

Increased lung vascular permeability and alveolar edema are cardinal features of inflammatory conditions such as acute respiratory distress syndrome (ARDS) and ventilator-induced lung injury (VILI). We previously demonstrated that pre-B-cell colony-enhancing factor (PBEF)/NAMPT, the proinflammatory cytokine encoded by NAMPT, participates in ARDS and VILI inflammatory syndromes. The present study evaluated posttranscriptional regulation of PBEF/NAMPT gene expression in human lung endothelium via 3'-untranslated region (UTR) microRNA (miRNA) binding. In silico analysis identified hsa-miR-374a and hsa-miR-568 as potential miRNA candidates. Increased PBEF/NAMPT transcription (by RT-PCR) and expression (by Western blotting) induced by 18% cyclic stretch (CS) (2 h: 3.4 ± 0.06 mRNA fold increase (FI); 10 h: 1.5 ± 0.06 protein FI) and by LPS (4 h: 3.8 ± 0.2 mRNA FI; 48 h: 2.6 ± 0.2 protein FI) were significantly attenuated by transfection with mimics of hsa-miR-374a or hsa-miR-568 (40-60% reductions each). LPS and 18% CS increased the activity of a PBEF/NAMPT 3'-UTR luciferase reporter (2.4-3.25 FI) with induction reduced by mimics of each miRNA (44-60% reduction). Specific miRNA inhibitors (antagomirs) for each PBEF/NAMPT miRNA significantly increased the endogenous PBEF/NAMPT mRNA (1.4-3.4 ± 0.1 FI) and protein levels (1.2-1.4 ± 0.1 FI) and 3'-UTR luciferase activity (1.4-1.7 ± 0.1 FI) compared with negative antagomir controls. Collectively, these data demonstrate that increased PBEF/NAMPT expression induced by bioactive agonists (i.e., excessive mechanical stress, LPS) involves epigenetic regulation with hsa-miR-374a and hsa-miR-568, representing novel therapeutic strategies to reduce inflammatory lung injury.

Ex vivo culture of primary human colonic tissue for studying transcriptional responses to 1α,25(OH)2 and 25(OH) vitamin D

1α,25-Dihydroxyvitamin D3 [1α,25(OH)2D3] is a steroid hormone derived from circulating 25(OH) vitamin D [25(OH)D] with chemopreventive effects in colorectal cancer. 1α,25(OH)2D3 acts through transcriptional mechanisms; however, our understanding of vitamin D transcriptional responses in the colon is derived from studies in transformed cancer cell lines which may not represent responses in normal healthy tissue. Here, we describe the optimization of an ex vivo culture model using primary colonic biopsy samples for studying short-term transcriptional response induced by 1α,25(OH)2D3 and 25(OH)D treatment. Colon biopsy samples from healthy subjects were maintained in primary culture and treated in parallel with 100 nM 1α,25(OH)2D3 or 62.5 nM 25(OH)D and vehicle control (ethanol). Viability was assessed using histology and enzymatic assays. Genome-wide transcriptional responses to 1α,25(OH)2D3 were assessed and expression of 25(OH)D targets CYP27B1 and CYP24A1 were measured by real time PCR. We show that ex vivo culture of colonic tissue remains viable for up to 8 h. The largest number of differentially expressed genes in response to 1α,25(OH)2D3 was noted after 6 h (n = 120). As proof of concept, the top upregulated gene was CYP24A1, a well-established vitamin D-responsive gene. With 25(OH)D treatment, mRNA expression of CYP27B1 was significantly increased after 1 h, while expression of CYP24A1 was greatest at 8 h. Ex vivo culture can be used to assess short-term transcriptional responses to 1α,25(OH)2D3 and 25(OH)D in primary tissue from human colon. Future studies will address interindividual differences in transcriptional responses.

Excessive mechanical stress increases HMGB1 expression in human lung microvascular endothelial cells via STAT3

Ventilator-induced lung injury (VILI) occurs when the lung parenchyma and vasculature are exposed to repetitive and excessive mechanical stress via mechanical ventilation utilized as supportive care for the adult respiratory distress syndrome (ARDS). VILI induces gene expression and systemic release of inflammatory mediators that contribute to the multi-organ dysfunction and morbidity and mortality of ARDS. HMGB1, an intracellular transcription factor with cytokine properties, is a late mediator in sepsis and ARDS pathobiology, however, the role of HMGB1 in VILI remains poorly described. We now report HMGB1 expression in human lung microvessel endothelial cells (ECs) exposed to excessive, equibiaxial mechanical stress, an in vitro correlate of VILI. We determined that high amplitude cyclic stretch (18% CS) increased HMGB1 expression (2-4-fold) via a signaling pathway with critical involvement of the transcription factor, STAT3. Concomitant exposure to 18% CS and oxidative stress (H₂O₂) augmented HMGB1 expression (~13 fold increase) whereas lipopolysaccharide (LPS) challenge increased HMGB1 expression in static EC, but not in 18% CS-challenged EC. In contrast, physiologic, low amplitude cyclic stretch (5% CS) attenuated both oxidative H₂O₂- and LPS-induced increases in HMGB1 expression, suggesting that physiologic mechanical stress is protective. These results indicate that HMGB1 gene expression is markedly responsive to VILI-mediated mechanical stress, an effect that is augmented by oxidative stress. We speculate that VILI-induced HMGB1 expression acts locally to increase vascular permeability and alveolar flooding, thereby exacerbating systemic inflammatory responses and increasing the likelihood of multi-organ dysfunction.

MicroRNA regulation of nonmuscle myosin light chain kinase expression in human lung endothelium

Increased lung vascular permeability, the consequence of endothelial cell (EC) barrier dysfunction, is a cardinal feature of inflammatory conditions such as acute lung injury and sepsis and leads to lethal physiological dysfunction characterized by alveolar flooding, hypoxemia, and pulmonary edema. We previously demonstrated that the nonmuscle myosin light chain kinase isoform (nmMLCK) plays a key role in agonist-induced pulmonary EC barrier regulation. The present study evaluated posttranscriptional regulation of MYLK expression, the gene encoding nmMLCK, via 3' untranslated region (UTR) binding by microRNAs (miRNAs) with in silico analysis identifying hsa-miR-374a, hsa-miR-374b, hsa-miR-520c-3p, and hsa-miR-1290 as miRNA candidates. We identified increased MYLK gene transcription induced by TNF-α (24 h; 4.7 ± 0.45 fold increase [FI]), LPS (4 h; 2.85 ± 0.15 [FI]), and 18% cyclic stretch (24 h; 4.6 ± 0.24 FI) that was attenuated by transfection of human lung ECs with mimics of hsa-miR-374a, hsa-miR-374b, hsa-miR-520c-3p, or hsa-miR-1290 (20-80% reductions by each miRNA). TNF-α, LPS, and 18% cyclic stretch each increased the activity of a MYLK 3'UTR luciferase reporter (2.5-7.0 FI) with induction reduced by mimics of each miRNA (30-60% reduction). MiRNA inhibitors (antagomirs) for each MYLK miRNA significantly increased 3'UTR luciferase activity (1.2-2.3 FI) and rescued the decreased MLCK-3'UTR reporter activity produced by miRNA mimics (70-110% increases for each miRNA; P < 0.05). These data demonstrate that increased human lung EC expression of MYLK by bioactive agonists (excessive mechanical stress, LPS, TNF-α) is regulated in part by specific miRNAs (hsa-miR-374a, hsa-miR-374b, hsa-miR-520c-3p, and hsa-miR-1290), representing a novel therapeutic strategy for reducing inflammatory lung injury.