Tumor-intrinsic and Cancer-associated Fibroblast Subtypes Independently Predict Outcomes in Pancreatic Cancer.

OBJECTIVE: To assess the utility of tumor-intrinsic and cancer-associated fibroblast (CAF) subtypes of pancreatic ductal adenocarcinoma (PDAC) in predicting response to neoadjuvant therapy (NAT) and overall survival. BACKGROUND: PDAC remains a deadly disease with limited treatment options, and both the tumor as well as the microenvironment play an important role in pathogenesis. Gene expression-based tumor-intrinsic subtypes (classical and basal-like) have been shown to predict outcomes, but tumor microenvironment subtypes are still evolving. METHODS: RNA-sequencing was performed on 114 deidentified resected PDAC tumors. Clinical data were collected by retrospective chart review. Single sample classifiers (SSCs) were used to determine classical and basal-like subtypes as well as tumor-permissive permCAF and tumor-restraining restCAF subtypes. Survival was analyzed using log-rank test. RESULTS: Patients who received NAT had an increase in overall survival (OS), with median survival of 27.9 months compared to 20.1 months for those who did not receive NAT, but the difference did not reach statistical significance (HR 0.64, P=0.076). Either tumor-intrinsic or CAF subtypes alone were associated with OS regardless of NAT or no NAT, and patients with classical or restCAF subtype had the best outcomes. When evaluated together, patients with classical-restCAF subtype had the best OS and basal-permCAF the worst OS (P<0.0001). NAT patients with classical-restCAF subtype demonstrated the longest OS compared to the other groups (P=0.00041). CONCLUSIONS: CAF subtypes have an additive effect over tumor-intrinsic subtypes in predicting survival with or without neoadjuvant FOLFIRINOX in PDAC. Molecular subtyping of both tumor and CAF compartments of PDAC may be important steps in selecting first-line systemic therapy.

Efficient computation of high-dimensional penalized generalized linear mixed models by latent factor modeling of the random effects.

Modern biomedical datasets are increasingly high-dimensional and exhibit complex correlation structures. Generalized linear mixed models (GLMMs) have long been employed to account for such dependencies. However, proper specification of the fixed and random effects in GLMMs is increasingly difficult in high dimensions, and computational complexity grows with increasing dimension of the random effects. We present a novel reformulation of the GLMM using a factor model decomposition of the random effects, enabling scalable computation of GLMMs in high dimensions by reducing the latent space from a large number of random effects to a smaller set of latent factors. We also extend our prior work to estimate model parameters using a modified Monte Carlo Expectation Conditional Minimization algorithm, allowing us to perform variable selection on both the fixed and random effects simultaneously. We show through simulation that through this factor model decomposition, our method can fit high-dimensional penalized GLMMs faster than comparable methods and more easily scale to larger dimensions not previously seen in existing approaches.

Determination of permissive and restraining cancer-associated fibroblast (DeCAF) subtypes.

Cancer-associated fibroblast (CAF) subpopulations in pancreatic ductal adenocarcinoma (PDAC) have been identified using single-cell RNA sequencing (scRNAseq) with divergent characteristics, but their clinical relevance remains unclear. We translate scRNAseq-derived CAF cell-subpopulation-specific marker genes to bulk RNAseq data, and develop a single- sample classifier, DeCAF, for the classification of clinically rest raining and perm issive CAF subtypes. We validate DeCAF in 19 independent bulk transcriptomic datasets across four tumor types (PDAC, mesothelioma, bladder and renal cell carcinoma). DeCAF subtypes have distinct histology features, immune landscapes, and are prognostic and predict response to therapy across cancer types. We demonstrate that DeCAF is clinically replicable and robust for the classification of CAF subtypes in patients for multiple tumor types, providing a better framework for the future development and translation of therapies against permissive CAF subtypes and preservation of restraining CAF subtypes. Significance: We introduce a replicable and robust classifier, DeCAF, that delineates the significance of the role of permissive and restraining CAF subtypes in cancer patients. DeCAF is clinically tractable, prognostic and predictive of treatment response in multiple cancer types and lays the translational groundwork for the preclinical and clinical development of CAF subtype specific therapies.

Tumor Intrinsic Subtypes and Gene Expression Signatures in Early-Stage ERBB2/HER2-Positive Breast Cancer: A Pooled Analysis of CALGB 40601, NeoALTTO, and NSABP B-41 Trials.

Importance: Biologic features may affect pathologic complete response (pCR) and event-free survival (EFS) after neoadjuvant chemotherapy plus ERBB2/HER2 blockade in ERBB2/HER2-positive early breast cancer (EBC). Objective: To define the quantitative association between pCR and EFS by intrinsic subtype and by other gene expression signatures in a pooled analysis of 3 phase 3 trials: CALGB 40601, NeoALTTO, and NSABP B-41. Design, Setting, and Participants: In this retrospective pooled analysis, 1289 patients with EBC received chemotherapy plus either trastuzumab, lapatinib, or the combination, with a combined median follow-up of 5.5 years. Gene expression profiling by RNA sequencing was obtained from 758 samples, and intrinsic subtypes and 618 gene expression signatures were calculated. Data analyses were performed from June 1, 2020, to January 1, 2023. Main Outcomes and Measures: The association of clinical variables and gene expression biomarkers with pCR and EFS were studied by logistic regression and Cox analyses. Results: In the pooled analysis, of 758 women, median age was 49 years, 12% were Asian, 6% Black, and 75% were White. Overall, pCR results were associated with EFS in the ERBB2-enriched (hazard ratio [HR], 0.45; 95% CI, 0.29-0.70; P < .001) and basal-like (HR, 0.19; 95% CI, 0.04-0.86; P = .03) subtypes but not in luminal A or B tumors. Dual trastuzumab plus lapatinib blockade over trastuzumab alone had a trend toward EFS benefit in the intention-to-treat population; however, in the ERBB2-enriched subtype there was a significant and independent EFS benefit of trastuzumab plus lapatinib vs trastuzumab alone (HR, 0.47; 95% CI, 0.27-0.83; P = .009). Overall, 275 of 618 gene expression signatures (44.5%) were significantly associated with pCR and 9 of 618 (1.5%) with EFS. The ERBB2/HER2 amplicon and multiple immune signatures were significantly associated with pCR. Luminal-related signatures were associated with lower pCR rates but better EFS, especially among patients with residual disease and independent of hormone receptor status. There was significant adjusted HR for pCR ranging from 0.45 to 0.81 (higher pCR) and 1.21-1.94 (lower pCR rate); significant adjusted HR for EFS ranged from 0.71 to 0.94. Conclusions and relevance: In patients with ERBB2/HER2-positive EBC, the association between pCR and EFS differed by tumor intrinsic subtype, and the benefit of dual ERBB2/HER2 blockade was limited to ERBB2-enriched tumors. Immune-activated signatures were concordantly associated with higher pCR rates and better EFS, whereas luminal signatures were associated with lower pCR rates.

glmmPen: High Dimensional Penalized Generalized Linear Mixed Models.

Generalized linear mixed models (GLMMs) are widely used in research for their ability to model correlated outcomes with non-Gaussian conditional distributions. The proper selection of fixed and random effects is a critical part of the modeling process, where model misspecification may lead to significant bias. However, the joint selection of fixed and random effects has historically been limited to lower dimensional GLMMs, largely due to the use of criterion-based model selection strategies. Here we present the R package glmmPen, one of the first to select fixed and random effects in higher dimension using a penalized GLMM modeling framework. Model parameters are estimated using a Monte Carlo expectation conditional minimization (MCECM) algorithm, which leverages Stan and RcppArmadillo for increased computational efficiency. Our package supports the Binomial, Gaussian, and Poisson families and multiple penalty functions. In this manuscript we discuss the modeling procedure, estimation scheme, and software implementation through application to a pancreatic cancer subtyping study. Simulation results show our method has good performance in selecting both the fixed and random effects in high dimensional GLMMs.

Association of PIK3CA Mutation With Pathologic Complete Response and Outcome by Hormone Receptor Status and Intrinsic Subtype in Early-Stage ERBB2/HER2-Positive Breast Cancer.

Importance: PIK3CA mutations may be associated with outcomes of patients with ERBB2/HER2-positive early breast cancer (EBC). Objectives: To assess if PIK3CA mutations among patients with ERBB2/HER2-positive EBC are associated with treatment response and outcome, and if these associations vary by hormone receptor (HR) status or intrinsic molecular subtype (IMS). Design, Setting, and Participants: This cohort study derived data on 184 patients from the phase 3 neoadjuvant Cancer and Leukemia Group B (CALGB) 40601 trial that enrolled patients with ERBB2/HER2-positive EBC in North America between January 1, 2008, and December 31, 2012. Participants received neoadjuvant paclitaxel with trastuzumab, lapatinib, or both. Statistical analysis was performed from March 23, 2022, to March 9, 2023. Exposures: Gene expression profiling by RNA sequencing with Prediction Analysis of Microarray 50-determined IMS and PIK3CA mutations from whole-exome sequencing were obtained from pretreatment biopsies from 184 of 305 trial participants. Main Outcomes and Measures: The primary end point was pathologic complete response (pCR) and the secondary end point of event-free survival (EFS). The association of PIK3CA mutations with pCR and EFS by HR status and IMS was estimated using logistic and Cox proportional hazards regression models. Results: All 184 participants were women, with a median age of 49 years (range 24-75 years). A total of 121 participants (66%) had clinical stage II tumors; 32 (17%) had PIK3CA mutations, most frequently H1047R (38% [12 of 32]) and E545K (22% [7 of 32]). PIK3CA mutations were present in 20 of 102 cases of HR-positive EBC (20%) and 12 of 82 cases HR-negative EBC (15%) and varied by IMS (luminal B, 9 of 25 [36%]; luminal A, 2 of 21 [10%]; and ERBB2/HER2-enriched tumors, 19 of 102 [19%]). Pathologic complete response rates were lower in PIK3CA mutated than PIK3CA wild type in the overall population (34% [11 of 32] vs 49% [74 of 152]; P = .14) and were significantly different among those receiving trastuzumab (30% [7 of 23] vs 54% [63 of 117]; P = .045). At a median follow-up of 9 years, PIK3CA mutations were significantly associated with worse EFS in the overall cohort (hazard ratio, 2.58 [95% CI, 1.24-5.35]; P = .01), which persisted in a multivariable model including pCR, HR status, stage, and IMS (hazard ratio, 2.52 [95% CI, 1.16-5.47]; P = .02). The negative association of PIK3CA mutation was significant in HR-positive (hazard ratio, 3.60 [95% CI, 1.45-8.96]; P = .006) and luminal subtypes (hazard ratio, 4.84 [95% CI, 1.08-21.70]; P = .04), but not in nonluminal and HR-negative tumors. Conclusions and Relevance: In ERBB2/HER2-positive EBC, PIK3CA mutations were associated with lower pCR rates and independently associated with worse long-term EFS. These findings appear to be associated with PIK3CA mutations in HR-positive and luminal EBC.

RHOAL57V drives the development of diffuse gastric cancer through IGF1R-PAK1-YAP1 signaling.

Cancer-associated mutations in the guanosine triphosphatase (GTPase) RHOA are found at different locations from the mutational hotspots in the structurally and biochemically related RAS. Tyr42-to-Cys (Y42C) and Leu57-to-Val (L57V) substitutions are the two most prevalent RHOA mutations in diffuse gastric cancer (DGC). RHOAY42C exhibits a gain-of-function phenotype and is an oncogenic driver in DGC. Here, we determined how RHOAL57V promotes DGC growth. In mouse gastric organoids with deletion of Cdh1, which encodes the cell adhesion protein E-cadherin, the expression of RHOAL57V, but not of wild-type RHOA, induced an abnormal morphology similar to that of patient-derived DGC organoids. RHOAL57V also exhibited a gain-of-function phenotype and promoted F-actin stress fiber formation and cell migration. RHOAL57V retained interaction with effectors but exhibited impaired RHOA-intrinsic and GAP-catalyzed GTP hydrolysis, which favored formation of the active GTP-bound state. Introduction of missense mutations at KRAS residues analogous to Tyr42 and Leu57 in RHOA did not activate KRAS oncogenic potential, indicating distinct functional effects in otherwise highly related GTPases. Both RHOA mutants stimulated the transcriptional co-activator YAP1 through actin dynamics to promote DGC progression; however, RHOAL57V additionally did so by activating the kinases IGF1R and PAK1, distinct from the FAK-mediated mechanism induced by RHOAY42C. Our results reveal that RHOAL57V and RHOAY42C drive the development of DGC through distinct biochemical and signaling mechanisms.