skip to main content


Title: Genetic program activity delineates risk, relapse, and therapy responsiveness in multiple myeloma
Abstract

Despite recent advancements in the treatment of multiple myeloma (MM), nearly all patients ultimately relapse and many become refractory to multiple lines of therapies. Therefore, we not only need the ability to predict which patients are at high risk for disease progression but also a means to understand the mechanisms underlying their risk. Here, we report a transcriptional regulatory network (TRN) for MM inferred from cross-sectional multi-omics data from 881 patients that predicts how 124 chromosomal abnormalities and somatic mutations causally perturb 392 transcription regulators of 8549 genes to manifest in distinct clinical phenotypes and outcomes. We identified 141 genetic programs whose activity profiles stratify patients into 25 distinct transcriptional states and proved to be more predictive of outcomes than did mutations. The coherence of these programs and accuracy of our network-based risk prediction was validated in two independent datasets. We observed subtype-specific vulnerabilities to interventions with existing drugs and revealed plausible mechanisms for relapse, including the establishment of an immunosuppressive microenvironment. Investigation of the t(4;14) clinical subtype using the TRN revealed that 16% of these patients exhibit an extreme-risk combination of genetic programs (median progression-free survival of 5 months) that create a distinct phenotype with targetable genes and pathways.

 
more » « less
NSF-PAR ID:
10253197
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
npj Precision Oncology
Volume:
5
Issue:
1
ISSN:
2397-768X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Background Adaptive CD19-targeted chimeric antigen receptor (CAR) T-cell transfer has become a promising treatment for leukemia. Although patient responses vary across different clinical trials, reliable methods to dissect and predict patient responses to novel therapies are currently lacking. Recently, the depiction of patient responses has been achieved using in silico computational models, with prediction application being limited. Methods We established a computational model of CAR T-cell therapy to recapitulate key cellular mechanisms and dynamics during treatment with responses of continuous remission (CR), non-response (NR), and CD19-positive (CD19 + ) and CD19-negative (CD19 − ) relapse. Real-time CAR T-cell and tumor burden data of 209 patients were collected from clinical studies and standardized with unified units in bone marrow. Parameter estimation was conducted using the stochastic approximation expectation maximization algorithm for nonlinear mixed-effect modeling. Results We revealed critical determinants related to patient responses at remission, resistance, and relapse. For CR, NR, and CD19 + relapse, the overall functionality of CAR T-cell led to various outcomes, whereas loss of the CD19 + antigen and the bystander killing effect of CAR T-cells may partly explain the progression of CD19 − relapse. Furthermore, we predicted patient responses by combining the peak and accumulated values of CAR T-cells or by inputting early-stage CAR T-cell dynamics. A clinical trial simulation using virtual patient cohorts generated based on real clinical patient datasets was conducted to further validate the prediction. Conclusions Our model dissected the mechanism behind distinct responses of leukemia to CAR T-cell therapy. This patient-based computational immuno-oncology model can predict late responses and may be informative in clinical treatment and management. 
    more » « less
  2. Abstract

    Maintaining the intracellular iron concentration within the homeostatic range is vital to meet cellular metabolic needs and reduce oxidative stress. Previous research revealed that the haloarchaeon Halobacterium salinarum encodes four diphtheria toxin repressor (DtxR) family transcription factors (TFs) that together regulate the iron response through an interconnected transcriptional regulatory network (TRN). However, the conservation of the TRN and the metal specificity of DtxR TFs remained poorly understood. Here we identified and characterized the TRN of Haloferax volcanii for comparison. Genetic analysis demonstrated that Hfx. volcanii relies on three DtxR transcriptional regulators (Idr, SirR, and TroR), with TroR as the primary regulator of iron homeostasis. Bioinformatics and molecular approaches revealed that TroR binds a conserved cis-regulatory motif located ∼100 nt upstream of the start codon of iron-related target genes. Transcriptomics analysis demonstrated that, under conditions of iron sufficiency, TroR repressed iron uptake and induced iron storage mechanisms. TroR repressed the expression of one other DtxR TF, Idr. This reduced DtxR TRN complexity relative to that of Hbt. salinarum appeared correlated with natural variations in iron availability. Based on these data, we hypothesize that variable environmental conditions such as iron availability appear to select for increasing TRN complexity.

     
    more » « less
  3. null (Ed.)
    Hypertrophic cardiomyopathy (HCM) is the most common heritable heart disease. Although the genetic cause of HCM has been linked to mutations in genes encoding sarcomeric proteins, the ability to predict clinical outcomes based on specific mutations in HCM patients is limited. Moreover, how mutations in different sarcomeric proteins can result in highly similar clinical phenotypes remains unknown. Posttranslational modifications (PTMs) and alternative splicing regulate the function of sarcomeric proteins; hence, it is critical to study HCM at the level of proteoforms to gain insights into the mechanisms underlying HCM. Herein, we employed high-resolution mass spectrometry–based top-down proteomics to comprehensively characterize sarcomeric proteoforms in septal myectomy tissues from HCM patients exhibiting severe outflow track obstruction ( n = 16) compared to nonfailing donor hearts ( n = 16). We observed a complex landscape of sarcomeric proteoforms arising from combinatorial PTMs, alternative splicing, and genetic variation in HCM. A coordinated decrease of phosphorylation in important myofilament and Z-disk proteins with a linear correlation suggests PTM cross-talk in the sarcomere and dysregulation of protein kinase A pathways in HCM. Strikingly, we discovered that the sarcomeric proteoform alterations in the myocardium of HCM patients undergoing septal myectomy were remarkably consistent, regardless of the underlying HCM-causing mutations. This study suggests that the manifestation of severe HCM coalesces at the proteoform level despite distinct genotype, which underscores the importance of molecular characterization of HCM phenotype and presents an opportunity to identify broad-spectrum treatments to mitigate the most severe manifestations of this genetically heterogenous disease. 
    more » « less
  4. Abstract

    Breast carcinoma is the most common cancer among women worldwide that consists of a heterogeneous group of subtype diseases. The whole-slide images (WSIs) can capture the cell-level heterogeneity, and are routinely used for cancer diagnosis by pathologists. However, key driver genetic mutations related to targeted therapies are identified by genomic analysis like high-throughput molecular profiling. In this study, we develop a deep-learning model to predict the genetic mutations and biological pathway activities directly from WSIs. Our study offers unique insights into WSI visual interactions between mutation and its related pathway, enabling a head-to-head comparison to reinforce our major findings. Using the histopathology images from the Genomic Data Commons Database, our model can predict the point mutations of six important genes (AUC 0.68–0.85) and copy number alteration of another six genes (AUC 0.69–0.79). Additionally, the trained models can predict the activities of three out of ten canonical pathways (AUC 0.65–0.79). Next, we visualized the weight maps of tumor tiles in WSI to understand the decision-making process of deep-learning models via a self-attention mechanism. We further validated our models on liver and lung cancers that are related to metastatic breast cancer. Our results provide insights into the association between pathological image features, molecular outcomes, and targeted therapies for breast cancer patients.

     
    more » « less
  5. Abstract

    Triple‐negative breast cancer (TNBC) accounts for 15–25% of diagnosed breast cancers, and its lack of a clinically defined therapeutic target has caused patients to suffer from earlier relapse and higher mortality rates than patients with other breast cancer subtypes. MicroRNAs (miRNAs) are small non‐coding RNAs that regulate the expression of multiple genes through RNA interference to maintain normal tissue function. The tumor suppressor miR‐34a is downregulated in TNBC, and its loss‐of‐expression correlates with worse disease outcomes. Therefore, delivering miR‐34a mimics into TNBC cells is a promising strategy to combat disease progression. To achieve this goal, we synthesized layer‐by‐layer assembled nanoparticles (LbL NPs) comprised of spherical poly(lactic‐co‐glycolic acid) cores surrounded by alternating layers of poly‐L‐lysine (PLL) and miR‐34a. TNBC cells internalized these LbL NPs to a greater extent than polyplexes comprised of PLL and miRNA, and confocal microscopy showed that LbL NPs delivered a substantial fraction of miR‐34a cargo into the cytosol. This yielded robust suppression of the miR‐34a target genes CCND‐1, Notch‐1, Bcl‐2, Survivin, and MDR‐1, which reduced TNBC cell proliferation and induced cell cycle arrest. These data validate that miR‐34a delivery can impair TNBC cell function and support continued investigation of this platform for treatment of TNBC.

     
    more » « less