Kinome state is predictive of cell viability in pancreatic cancer tumor and cancer-associated fibroblast cell lines.

Numerous aspects of cellular signaling are regulated by the kinome-the network of over 500 protein kinases that guides and modulates information transfer throughout the cell. The key role played by both individual kinases and assemblies of kinases organized into functional subnetworks leads to kinome dysregulation driving many diseases, particularly cancer. In the case of pancreatic ductal adenocarcinoma (PDAC), a variety of kinases and associated signaling pathways have been identified for their key role in the establishment of disease as well as its progression. However, the identification of additional relevant therapeutic targets has been slow and is further confounded by interactions between the tumor and the surrounding tumor microenvironment. In this work, we attempt to link the state of the human kinome, or kinotype, with cell viability in treated, patient-derived PDAC tumor and cancer-associated fibroblast cell lines. We applied classification models to independent kinome perturbation and kinase inhibitor cell screen data, and found that the inferred kinotype of a cell has a significant and predictive relationship with cell viability. We further find that models are able to identify a set of kinases whose behavior in response to perturbation drive the majority of viability responses in these cell lines, including the understudied kinases CSNK2A1/3, CAMKK2, and PIP4K2C. We next utilized these models to predict the response of new, clinical kinase inhibitors that were not present in the initial dataset for model devlopment and conducted a validation screen that confirmed the accuracy of the models. These results suggest that characterizing the perturbed state of the human protein kinome provides significant opportunity for better understanding of signaling behavior and downstream cell phenotypes, as well as providing insight into the broader design of potential therapeutic strategies for PDAC.

Purity Independent Subtyping of Tumors Pancreatic Cancer Classifier: Analytic Validation of a 16-RNA Expression Signature Distinguishing Basal and Classical Subtypes.

The two major molecular subtypes of pancreatic adenocarcinoma reportedly have differential response to FOLFIRINOX-based therapy. To promote rapid assignment of basal versus classical subtypes, an array-based single-sample classifier assay was developed and applied to 74 formalin-fixed, paraffin-embedded biopsy or resection specimens of known subtype based on transcriptomics. The Purity Independent Subtyping of Tumors algorithm assigns subtype based on relative expression of 16 RNAs counted by RNA sequencing (RNAseq) versus more practical array-based NanoString nCounter Elements XT technology. Subtype calls were largely concordant between RNAseq and array methods (72/74, 97% agreement). Compared with the lengthy RNAseq protocol, the array-based assay takes just 3 working days to analyze, permitting rapid reporting of tumor subtype. In conclusion, the Purity Independent Subtyping of Tumors pancreatic cancer classifier has robust performance to classify pancreatic adenocarcinoma into basal versus classical subtypes. Clinical validation studies are underway to evaluate outcome in patients whose standard-of-care chemotherapy regimen is selected on the basis of rapid subtype assignment (NCT04683315).

Deeply-Learned Generalized Linear Models with Missing Data.

Deep Learning (DL) methods have dramatically increased in popularity in recent years, with significant growth in their application to various supervised learning problems. However, the greater prevalence and complexity of missing data in such datasets present significant challenges for DL methods. Here, we provide a formal treatment of missing data in the context of deeply learned generalized linear models, a supervised DL architecture for regression and classification problems. We propose a new architecture, dlglm, that is one of the first to be able to flexibly account for both ignorable and non-ignorable patterns of missingness in input features and response at training time. We demonstrate through statistical simulation that our method outperforms existing approaches for supervised learning tasks in the presence of missing not at random (MNAR) missingness. We conclude with a case study of the Bank Marketing dataset from the UCI Machine Learning Repository, in which we predict whether clients subscribed to a product based on phone survey data. Supplementary materials for this article are available online.

Effect of hydrothermal pretreatment on leachate fed Scenedesmus sp. biomass solubilization and biogas production.

The aim of the present study was to assess the effect of hydrothermal pretreatment on the solubilization and anaerobic digestion (AD) of Scenedesmus sp. biomass. At first, the microalgae was cultivated in 5% fresh leachate (FL) to recover nutrients such as nitrogen and phosphorus. Scenedesmus sp. grown in 5% FL obtained 100%, 77% and 97% removal efficiency of ammonium nitrogen (NH4+ - N), total Kjeldahl nitrogen (TKN) and phosphate phosphorous (PO43- -P), respectively. In the following step, the hydrothermal pretreatment of Scenedesmus sp. biomass was carried out at 120, 150 and 170 °C and retention time of 0, 30 and 60 min to evaluate its solubilization and biogas production through AD in batch test. Soluble chemical oxygen demand (sCOD) increased by 260% compared to untreated microalgae at 170 °C for 60 min. In comparison to untreated microalgae, the highest increase in biogas (70%) and methane yield (100%) was observed for 150 °C and 60 min pretreated microalgae as a consequence of hydrothermal pretreatment. Hydrothermal pretreatment has shown effectiveness in enhancing biomass solubilization and increasing biogas yield. Nevertheless, further research at the pilot scale is necessary to thoroughly evaluate the potential and feasibility of hydrothermal pretreatment for full-scale implementation.

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.

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.

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.

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.

Evaluation of cell disruption methods for protein and coenzyme Q10 quantification in purple non-sulfur bacteria.

A recent focus has been on the recovery of single-cell protein and other nutritionally valuable bioproducts, such as Coenzyme Q10 (CoQ10) from purple non-sulfur bacteria (PNSB) biomass following wastewater treatment. However, due to PNSB's peculiar cell envelope (e.g., increased membrane cross-section for energy transduction) and relatively smaller cell size compared to well-studied microbial protein sources like yeast and microalgae, the effectiveness of common cell disruption methods for protein quantification from PNSB may differ. Thus, this study examines the efficiency of selected chemical (NaOH and EDTA), mechanical (homogenization and bead milling), physical (thermal and bath/probe sonication), and combined chemical-mechanical/physical treatment techniques on the PNSB cell lysis. PNSB biomass was recovered from the treatment of gas-to-liquid process water. Biomass protein and CoQ10 contents were quantified based on extraction efficiency. Considering single-treatment techniques, bead milling resulted in the best protein yields (p < 0.001), with the other techniques resulting in poor yields. However, the NaOH-assisted sonication (combined chemical/physical treatment technique) resulted in similar protein recovery (p = 1.00) with bead milling, with the former having a better amino acid profile. For example, close to 50% of the amino acids, such as sensitive ones like tryptophan, threonine, cystine, and methionine, were detected in higher concentrations in NaOH-assisted sonication (>10% relative difference) compared to bead-milling due to its less disruptive nature and improved solubility of amino acids in alkaline conditions. Overall, PNSB required more intensive protein extraction techniques than were reported to be effective on other single-cell organisms. NaOH was the preferred chemical for chemical-aided mechanical/physical extraction as EDTA was observed to interfere with the Lowry protein kit, resulting in significantly lower concentrations. However, EDTA was the preferred chemical agent for CoQ10 extraction and quantification. CoQ10 extraction efficiency was also suspected to be adversely influenced by pH and temperature.

Differential transcript usage analysis incorporating quantification uncertainty via compositional measurement error regression modeling.

Differential transcript usage (DTU) occurs when the relative expression of multiple transcripts arising from the same gene changes between different conditions. Existing approaches to detect DTU often rely on computational procedures that can have speed and scalability issues as the number of samples increases. Here we propose a new method, CompDTU, that uses compositional regression to model the relative abundance proportions of each transcript that are of interest in DTU analyses. This procedure leverages fast matrix-based computations that make it ideally suited for DTU analysis with larger sample sizes. This method also allows for the testing of and adjustment for multiple categorical or continuous covariates. Additionally, many existing approaches for DTU ignore quantification uncertainty in the expression estimates for each transcript in RNA-seq data. We extend our CompDTU method to incorporate quantification uncertainty leveraging common output from RNA-seq expression quantification tool in a novel method CompDTUme. Through several power analyses, we show that CompDTU has excellent sensitivity and reduces false positive results relative to existing methods. Additionally, CompDTUme results in further improvements in performance over CompDTU with sufficient sample size for genes with high levels of quantification uncertainty, while also maintaining favorable speed and scalability. We motivate our methods using data from the Cancer Genome Atlas Breast Invasive Carcinoma data set, specifically using RNA-seq data from primary tumors for 740 patients with breast cancer. We show greatly reduced computation time from our new methods as well as the ability to detect several novel genes with significant DTU across different breast cancer subtypes.

DNA Mutational Profiling in Patients With Colorectal Cancer Treated With Standard of Care Reveals Differences in Outcome and Racial Distribution of Mutations.

PURPOSE: CALGB (Alliance)/SWOG 80405 was a randomized phase III trial that in first-line patients with metastatic colorectal cancer (mCRC) treated with bevacizumab or cetuximab with chemotherapy. We aimed to discover novel mutated genes associated with prognosis and differential response to therapy with the biologics. METHODS: Primary tumor DNA from 548 patients was sequenced using FoundationOne. The effect of mutated genes and mutations on overall survival (OS) was tested adjusting for microsatellite instability status, BRAF V600E, all RAS mutations, arm, sex, and age. RESULTS: The median number (lower-upper quartile) of mutated genes was 5 (3-7), 5 (3-6) in microsatellite stable and 12.5 (4.5-32) in microsatellite instability-high tumors. Mutated KRAS and APC were more frequent in Black (53% and 85%) than White (27% and 65%, respectively) patients while BRAF V600E was less frequent in Black (5%) than White (14%) patients. The median OS in patients with BRAF non-V600E (2.2% of patients) was 31.9 months (95% CI, 15.1 to not applicable [NA]) similar to that of BRAF wild-type (WT) patients (31.2 months [95% CI, 29.0 to 33.9]). Mutated LRP1B (10.7% of patients) was associated with improved OS compared with WT LRP1B (hazard ratio, 0.57 [95% CI, 0.40 to 0.80]). RNF43 (5.6% of patients) interacted with treatment arms as, in the cetuximab arm, patients with mutated RNF43 had a median OS of 11.5 (95% CI, 10.8 to NA) months compared with 30.1 (95% CI, 24.9 to 35.3) months in patients with WT RNF43, whereas in the bevacizumab arm, patients with mutated RNF43 had a median OS of 25.0 (95% CI, 14.2 to NA) months compared with 31.3 (95% CI, 29.0 to 34.3) months in patients with WT RNF43. CONCLUSION: These results can provide new tools to predict patient outcome and improve therapeutic decisions and trial participation in patient minorities. The molecular alterations identified in this study may direct biomarker-driven studies.

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.

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.

Sub-Cluster Identification through Semi-Supervised Optimization of Rare-Cell Silhouettes (SCISSORS) in single-cell RNA-sequencing.

MOTIVATION: Single-cell RNA-sequencing (scRNA-seq) has enabled the molecular profiling of thousands to millions of cells simultaneously in biologically heterogenous samples. Currently, the common practice in scRNA-seq is to determine cell type labels through unsupervised clustering and the examination of cluster-specific genes. However, even small differences in analysis and parameter choosing can greatly alter clustering results and thus impose great influence on which cell types are identified. Existing methods largely focus on determining the optimal number of robust clusters, which can be problematic for identifying cells of extremely low abundance due to their subtle contributions toward overall patterns of gene expression. RESULTS: Here, we present a carefully designed framework, SCISSORS, which accurately profiles subclusters within broad cluster(s) for the identification of rare cell types in scRNA-seq data. SCISSORS employs silhouette scoring for the estimation of heterogeneity of clusters and reveals rare cells in heterogenous clusters by a multi-step semi-supervised reclustering process. Additionally, SCISSORS provides a method for the identification of marker genes of high specificity to the cell type. SCISSORS is wrapped around the popular Seurat R package and can be easily integrated into existing Seurat pipelines. AVAILABILITY AND IMPLEMENTATION: SCISSORS, including source code and vignettes, are freely available at https://github.com/jr-leary7/SCISSORS.

Open-source curation of a pancreatic ductal adenocarcinoma gene expression analysis platform (pdacR) supports a two-subtype model.

Pancreatic ductal adenocarcinoma (PDAC) is an aggressive disease for which potent therapies have limited efficacy. Several studies have described the transcriptomic landscape of PDAC tumors to provide insight into potentially actionable gene expression signatures to improve patient outcomes. Despite centralization efforts from multiple organizations and increased transparency requirements from funding agencies and publishers, analysis of public PDAC data remains difficult. Bioinformatic pitfalls litter public transcriptomic data, such as subtle inclusion of low-purity and non-adenocarcinoma cases. These pitfalls can introduce non-specificity to gene signatures without appropriate data curation, which can negatively impact findings. To reduce barriers to analysis, we have created pdacR ( http://pdacR.bmi.stonybrook.edu , github.com/rmoffitt/pdacR), an open-source software package and web-tool with annotated datasets from landmark studies and an interface for user-friendly analysis in clustering, differential expression, survival, and dimensionality reduction. Using this tool, we present a multi-dataset analysis of PDAC transcriptomics that confirms the basal-like/classical model over alternatives.

Prognostic and Predictive Value of Immune-Related Gene Expression Signatures vs Tumor-Infiltrating Lymphocytes in Early-Stage ERBB2/HER2-Positive Breast Cancer: A Correlative Analysis of the CALGB 40601 and PAMELA Trials.

Importance: Both tumor-infiltrating lymphocytes (TILs) assessment and immune-related gene expression signatures by RNA profiling predict higher pathologic complete response (pCR) and improved event-free survival (EFS) in patients with early-stage ERBB2/HER2-positive breast cancer. However, whether these 2 measures of immune activation provide similar or additive prognostic value is not known. Objective: To examine the prognostic ability of TILs and immune-related gene expression signatures, alone and in combination, to predict pCR and EFS in patients with early-stage ERBB2/HER2-positive breast cancer treated in 2 clinical trials. Design, Setting, and Participants: In this prognostic study, a correlative analysis was performed on the Cancer and Leukemia Group B (CALGB) 40601 trial and the PAMELA trial. In the CALGB 40601 trial, 305 patients were randomly assigned to weekly paclitaxel with trastuzumab, lapatinib, or both for 16 weeks. The primary end point was pCR, with a secondary end point of EFS. In the PAMELA trial, 151 patients received neoadjuvant treatment with trastuzumab and lapatinib for 18 weeks. The primary end point was the ability of the HER2-enriched subtype to predict pCR. The studies were conducted from October 2013 to November 2015 (PAMELA) and from December 2008 to February 2012 (CALGB 40601). Data analyses were performed from June 1, 2020, to January 1, 2022. Main Outcomes and Measures: Immune-related gene expression profiling by RNA sequencing and TILs were assessed on 230 CALGB 40601 trial pretreatment tumors and 138 PAMELA trial pretreatment tumors. The association of these biomarkers with pCR (CALGB 40601 and PAMELA) and EFS (CALGB 40601) was studied by logistic regression and Cox analyses. Results: The median age of the patients was 50 years (IQR, 42-50 years), and 305 (100%) were women. Of 202 immune signatures tested, 166 (82.2%) were significantly correlated with TILs. In both trials combined, TILs were significantly associated with pCR (odds ratio, 1.01; 95% CI, 1.01-1.02; P = .02). In addition to TILs, 36 immune signatures were significantly associated with higher pCR rates. Seven of these signatures outperformed TILs for predicting pCR, 6 of which were B-cell related. In a multivariable Cox model adjusted for clinicopathologic factors, including PAM50 intrinsic tumor subtype, the immunoglobulin G signature, but not TILs, was independently associated with EFS (immunoglobulin G signature-adjusted hazard ratio, 0.63; 95% CI, 0.42-0.93; P = .02; TIL-adjusted hazard ratio, 1.00; 95% CI, 0.98-1.02; P = .99). Conclusions and Relevance: Results of this study suggest that multiple B-cell-related signatures were more strongly associated with pCR and EFS than TILs, which largely represent T cells. When both TILs and gene expression are available, the prognostic value of immune-related signatures appears to be superior.

Estimating cell type composition using isoform expression one gene at a time.

Human tissue samples are often mixtures of heterogeneous cell types, which can confound the analyses of gene expression data derived from such tissues. The cell type composition of a tissue sample may itself be of interest and is needed for proper analysis of differential gene expression. A variety of computational methods have been developed to estimate cell type proportions using gene-level expression data. However, RNA isoforms can also be differentially expressed across cell types, and isoform-level expression could be equally or more informative for determining cell type origin than gene-level expression. We propose a new computational method, IsoDeconvMM, which estimates cell type fractions using isoform-level gene expression data. A novel and useful feature of IsoDeconvMM is that it can estimate cell type proportions using only a single gene, though in practice we recommend aggregating estimates of a few dozen genes to obtain more accurate results. We demonstrate the performance of IsoDeconvMM using a unique data set with cell type-specific RNA-seq data across more than 135 individuals. This data set allows us to evaluate different methods given the biological variation of cell type-specific gene expression data across individuals. We further complement this analysis with additional simulations.

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.

Gene signatures derived from transcriptomic-causal networks stratified colorectal cancer patients for effective targeted therapy.

Predictive and prognostic gene signatures derived from interconnectivity among genes can tailor clinical care to patients in cancer treatment. We identified gene interconnectivity as the transcriptomic-causal network by integrating germline genotyping and tumor RNA-seq data from 1,165 patients with metastatic colorectal cancer (CRC). The patients were enrolled in a clinical trial with randomized treatment, either cetuximab or bevacizumab in combination with chemotherapy. We linked the network to overall survival (OS) and detected novel biomarkers by controlling for confounding genes. Our data-driven approach discerned sets of genes, each set collectively stratify patients based on OS. Two signatures under the cetuximab treatment were related to wound healing and macrophages. The signature under the bevacizumab treatment was related to cytotoxicity and we replicated its effect on OS using an external cohort. We also showed that the genes influencing OS within the signatures are downregulated in CRC tumor vs. normal tissue using another external cohort. Furthermore, the corresponding proteins encoded by the genes within the signatures interact each other and are functionally related. In conclusion, this study identified a group of genes that collectively stratified patients based on OS and uncovered promising novel prognostic biomarkers for personalized treatment of CRC using transcriptomic causal networks.

Kinetics of Biodiesel Production from Microalgae Using Microbubble Interfacial Technology.

As an alternative to fossil fuels, biodiesel can be a source of clean and environmentally friendly energy source. However, its commercial application is limited by expensive feedstock and the slow nature of the pretreatment step-acid catalysis. The conventional approach to carry out this reaction uses stirred tank reactors. Recently, the lab-scale experiments using microbubble mediated mass transfer technology have demonstrated its potential use at commercial scale. However, all the studies conducted so far have been at a lab scale~100 mL of feedstock. To analyze the feasibility of microbubble technology, a larger pilot scale study is required. In this context, a kinetic study of microbubble technology at an intermediate scale is conducted (3 L of oil). Owing to the target for industrial application of the process, a commercial feedstock (Spirulina), microalgae oil (MO) and a commercial catalyst para-toluene sulfonic acid (PTSA) are used. Experiments to characterize the kinetics space (response surface, RSM) required for up-scaling are designed to develop a robust model. The model is compared with that developed by the gated recurrent unit (GRU) method. The maximum biodiesel conversion of 99.45 ± 1.3% is achieved by using these conditions: the molar ratio of MO to MeOH of 1:23.73 ratio, time of 60 min, and a catalyst loading of 3.3 wt% MO with an MO volume of 3 L. Furthermore, predicted models of RSM and GRU show proper fits to the experimental result. It was found that GRU produced a more accurate and robust model with correlation coefficient R2 = 0.9999 and root-mean-squared error (RSME) = 0.0515 in comparison with RSM model with R2 = 0.9844 and RMSE = 3.0832, respectively. Although RSM and GRU are fully empirical representations, they can be used for reactor up-scaling horizontally with microbubbles if the liquid layer height is held constant while the microbubble injection replicates along the floor of the reactor vessel-maintaining the tessellation pattern of the smaller vessel. This scaling approach maintains the local mixing profile, which is the major uncontrolled variable in conventional stirred tank reactor up-scaling.