ABSTRACTType 2 diabetes mellitus (T2DM), characterized by insulin resistance, is closely associated with Alzheimer's disease (AD). Cerebrovascular dysfunction is manifested in both T2DM and AD, and is often considered as a pathological link between the two diseases. Insulin signaling regulates critical functions of the blood–brain barrier (BBB), and endothelial insulin resistance could lead to BBB dysfunction, aggravating AD pathology. However, insulin signaling is intrinsically dynamic and involves interactions among numerous molecular mediators. Hence, a mechanistic systems biology model is needed to understand how insulin regulates BBB physiology and the consequences of its impairment in T2DM and AD. In this study, we investigated the pharmacodynamic effect of insulin on the expression of vascular cell adhesion molecule 1 (VCAM1), a marker of cerebrovascular inflammation. Intriguingly, normal insulin concentrations selectively activated the PI3K‐AKT pathway, leading to decreased VCAM1 expression. However, exposure to supraphysiological insulin levels, which is present in insulin resistance, activated both PI3K‐AKT and MEK–ERK pathways, and increased VCAM1 expression. We developed a mathematical model that adequately described the dynamics of various insulin signaling nodes and VCAM1 expression. Further, the model was integrated with in vitro proteomics and transcriptomics data from AD patients to simulate VCAM1 expression under pathological conditions. This approach allowed us to establish a quantitative systems pharmacology framework to investigate BBB dysfunction in AD and metabolic syndrome, thereby offering opportunities to identify specific disruptions in molecular networks that will enable us to identify novel therapeutic targets.
OmicsFootPrint: a framework to integrate and interpret multi-omics data using circular images and deep neural networks
Abstract The OmicsFootPrint framework addresses the need for advanced multi-omics data analysis methodologies by transforming data into intuitive two-dimensional circular images and facilitating the interpretation of complex diseases. Utilizing deep neural networks and incorporating the SHapley Additive exPlanations algorithm, the framework enhances model interpretability. Tested with The Cancer Genome Atlas data, OmicsFootPrint effectively classified lung and breast cancer subtypes, achieving high area under the curve (AUC) scores—0.98 ± 0.02 for lung cancer subtype differentiation and 0.83 ± 0.07 for breast cancer PAM50 subtypes, and successfully distinguished between invasive lobular and ductal carcinomas in breast cancer, showcasing its robustness. It also demonstrated notable performance in predicting drug responses in cancer cell lines, with a median AUC of 0.74, surpassing nine existing methods. Furthermore, its effectiveness persists even with reduced training sample sizes. OmicsFootPrint marks an enhancement in multi-omics research, offering a novel, efficient and interpretable approach that contributes to a deeper understanding of disease mechanisms.
Luminal androgen receptor breast cancer subtype and investigation of the microenvironment and neoadjuvant chemotherapy response
Abstract Triple-negative breast cancer (TNBC) is the most aggressive breast cancer subtype with low overall survival rates and high molecular heterogeneity; therefore, few targeted therapies are available. The luminal androgen receptor (LAR) is the most consistently identified TNBC subtype, but the clinical utility has yet to be established. Here, we constructed a novel genomic classifier, LAR-Sig, that distinguishes the LAR subtype from other TNBC subtypes and provide evidence that it is a clinically distinct disease. A meta-analysis of seven TNBC datasets (n = 1086 samples) from neoadjuvant clinical trials demonstrated that LAR patients have significantly reduced response (pCR) rates than non-LAR TNBC patients (odds ratio = 2.11, 95% CI: 1.33, 2.89). Moreover, deconvolution of the tumor microenvironment confirmed an enrichment of luminal epithelium corresponding with a decrease in basal and myoepithelium in LAR TNBC tumors. Increased immunosuppression in LAR patients may lead to a decreased presence of cycling T-cells and plasma cells. While, an increased presence of myofibroblast-like cancer-associated cells may impede drug delivery and treatment. In summary, the lower levels of tumor infiltrating lymphocytes (TILs), reduced immune activity in the micro-environment, and lower pCR rates after NAC, suggest that new therapeutic strategies for the LAR TNBC subtype need to be developed.
Characteristics and Spatially Defined Immune (micro)landscapes of Early-stage PD-L1–positive Triple-negative Breast Cancer
Abstract Purpose: Programmed death ligand 1 [PD-(L)1]-targeted therapies have shown modest survival benefit in triple-negative breast cancer (TNBC). PD-L1+ microenvironments in TNBC are not well characterized and may inform combinatorial immune therapies. Herein, we characterized clinicopathologic features, RNA-based immune signatures, and spatially defined protein-based tumor–immune microenvironments (TIME) in early-stage PD-L1+ and PD-L1− TNBC. Experimental Design: From a large cohort of chemotherapy-naïve TNBC, clinicopathologic features, deconvoluted RNA immune signatures, and intraepithelial and stromal TIME (Nanostring GeoMX) were identified in subsets of PD-L1+ and PD-L1− TNBC, as defined by FDA-approved PD-L1 companion assays. Results: 228 of 499 (46%) TNBC were PD-L1+ (SP142: ≥1% immune cells-positive). Using PD-L1 22C3, 46% had combined positive score (CPS) ≥ 1 and 16% had CPS ≥10. PD-L1+ TNBC were higher grade with higher tumor-infiltrating lymphocytes (TIL; P < 0.05). PD-L1 was not associated with improved survival following adjustment for TILs and other variables. RNA profiles of PD-L1+ TNBC had increased dendritic cell, macrophage, and T/B cell subset features; and decreased myeloid-derived suppressor cells. PD-L1+ stromal and intraepithelial TIMEs were highly enriched in IDO-1, HLA-DR, CD40, and CD163 compared with PD-L1-TIME, with spatially specific alterations in CTLA-4, Stimulator of Interferon Genes (STING), and fibronectin. Macrophage- and antigen presentation–related proteins correlated most strongly with PD-L1 protein. Conclusions: In this early-stage TNBC cohort, nearly 50% were PD-L1+ (SP142 companion assay) while 16% were PD-L1+ with the 22C3 companion assay. PD-L1+ TNBC had specific myeloid-derived and lymphoid features. Spatially defined PD-L1+ TIME were enriched in several clinically actionable immune proteins. These data may inform future studies on combinatorial immunotherapies for patients with PD-L1+ TNBC. See related commentary by Symmans, p. 5446
Circular RNAs and their associations with breast cancer subtypes
Abstract Rapid development of next generation sequencing technology has enabled the identification of genomic alterations from short sequencing reads. There are a number of software pipelines available for calling single nucleotide variants from genomic DNA but, no comprehensive pipelines to identify, annotate and prioritize expressed SNVs (eSNVs) from non-directional paired-end RNA-Seq data. We have developed the eSNV-Detect, a novel computational system, which utilizes data from multiple aligners to call, even at low read depths, and rank variants from RNA-Seq. Multi-platform comparisons with the eSNV-Detect variant candidates were performed. The method was first applied to RNA-Seq from a lymphoblastoid cell-line, achieving 99.7% precision and 91.0% sensitivity in the expressed SNPs for the matching HumanOmni2.5 BeadChip data. Comparison of RNA-Seq eSNV candidates from 25 ER+ breast tumors from The Cancer Genome Atlas (TCGA) project with whole exome coding data showed 90.6–96.8% precision and 91.6–95.7% sensitivity. Contrasting single-cell mRNA-Seq variants with matching traditional multicellular RNA-Seq data for the MD-MB231 breast cancer cell-line delineated variant heterogeneity among the single-cells. Further, Sanger sequencing validation was performed for an ER+ breast tumor with paired normal adjacent tissue validating 29 out of 31 candidate eSNVs. The source code and user manuals of the eSNV-Detect pipeline for Sun Grid Engine and virtual machine are available at http://bioinformaticstools.mayo.edu/research/esnv-detect/.
AT‐101 downregulates BCL2 and MCL1 and potentiates the cytotoxic effects of lenalidomide and dexamethasone in preclinical models of multiple myeloma and Waldenström macroglobulinaemia
SummaryMultiple myeloma, the second most common haematological malignancy in the U.S., is currently incurable. Disruption of the intrinsic apoptotic pathway by BCL2 and MCL1 upregulation is observed in >80% of myeloma cases and is associated with an aggressive clinical course. Remarkably, there is no approved drug with the ability to target BCL2 or MCL1. Thus, we investigated the anti‐tumour effects of a pan‐BCL2 inhibitor, AT‐101, which has high binding specificity for BCL2 and MCL1 in preclinical models of plasma cell cancers (Multiple myeloma and Waldenström macroglobulinaemia). Gene expression and immunoblot analysis of six plasma cell cancer models showed upregulation of BCL2 family members. AT‐101 was able to downregulate BCL2 and MCL1 in all plasma cell cancer models and induced apoptotic cell death in a caspase‐dependent manner by altering mitochondrial membrane permeability. This cytotoxic effect and BCL2 downregulation were further potentiated when AT‐101 was combined with lenalidomide/dexamethasone (LDA). NanoString nCounter mRNA quantification and Ingenuity Pathways Analysis revealed differential changes in the CCNA2, FRZB, FYN, IRF1, PTPN11 genes in LDA‐treated cells. In summary, we describe for the first time the cellular and molecular events associated with the use of AT‐101 in combination with lenalidomide/dexamethasone in preclinical models of plasma cell malignancy.