Transcriptional Adaptation Research Articles

SelectRepair Knockout: Efficient PTC-Free Gene Knockout Through Selectable Homology-Directed DNA Repair.

Generating nonessential gene knockouts using CRISPR/Cas9 technology is becoming increasingly common in biological research. In a typical workflow, the Cas9 endonuclease is used to induce a DNA double-strand break that relies on nonhomologous end-joining (NHEJ) to introduce a premature termination codon (PTC) in the target gene. The goal is to isolate clones in which the gene produces PTC-containing mRNA transcripts that are degraded via nonsense-mediated mRNA decay (NMD) to cause loss of gene function. Unfortunately, this approach is laborious, and not all PTCs trigger NMD. More importantly, mounting evidence suggest that PTC mutations can also result in a transcriptional adaptation response that can mask the effects of a PTC-mediated gene knockout. In this chapter, we present a PTC-free gene knockout strategy that implements homology-directed DNA repair (HDR) with selectable markers to substantially reduce the complexity of the screening and validation of genome edits in cells containing more than one gene copy as in the case of the commonly used hypotriploid HEK293 cell line. We describe how to obtain a complete knockout of the Ligase IV protein (LIG4) and provide considerations for the application of this SelectRepair Knockout method to other genes.

Diversity of transcriptional regulatory adaptation in E. coli.

The Transcriptional Regulatory Network (TRN) in bacteria is thought to rapidly evolve in response to selection pressures, modulating transcription factor (TF) activities and interactions. In order to probe the limits and mechanisms surrounding the short-term adaptability of the TRN, we generated, evolved, and characterized knockout (KO) strains in E. coli for 11 regulators selected based on measured growth impact on glucose minimal media. All but one knockout strain (Δlrp) were able to recover growth and did so requiring few convergent mutations. We found that the TF knockout adaptations could be divided into four categories: 1) Strains (ΔargR, ΔbasR, Δlon, ΔzntR, Δzur) that recovered growth without any regulator-specific adaptations, likely due to minimal activity of the regulator on the growth condition, 2) Strains (ΔcytR, ΔmlrA, ΔybaO) that recovered growth without TF-specific mutations but with differential expression of regulators with overlapping regulons to the KO'ed TF, 3) Strains (Δcrp, Δfur) that recovered growth using convergent mutations within their regulatory networks, including regulated promoters and connected regulators, and 4) Strains (Δlrp) that were unable to fully recover growth, seemingly due to the broad connectivity of the TF within the TRN. Analyzing growth capabilities in evolved and unevolved strains indicated that growth adaptation can restore fitness to diverse substrates often despite a lack of TF-specific mutations. This work reveals the breadth of TRN adaptive mechanisms and suggests these mechanisms can be anticipated based on the network and functional context of the perturbed TFs.

TMIC-67. EMERGING ROLE OF EXTRACELLULAR MATRIX STIFFNESS AS A HALLMARK OF GLIOBLASTOMA: INSIGHT TO MECHANICAL MEMORY FORMATION AND PROLONGED INVASION IN THE TUMOR PERIPHERY

Abstract The extracellular matrix (ECM) plays a critical role in cancer progression, with its biophysical properties influencing cancer behavior. ECM stiffness can induce long-lasting adaptations in cancer cells, promoting sustained invasion through a process known as “mechanical memory” (MM). Glioblastoma (GBM) exhibits spatial heterogeneity in stiffness, with stiffer ECM found at the tumor periphery where invasive cells reside and promote recurrence. We hypothesize that increased stiffness at the tumor periphery primes GBM cells for sustained invasion through MM. METHODS Cell lines were established from image-guided human GBM samples from tumor margin and core (n=10 patients, 20 cell lines). Cells from random locations were conditioned in physiologically relevant extracellular matrices for two weeks. Human organotypic assay (n=3), Migration Assays, Atomic Force Microscopy (AFM), Immunocytochemistry, Western Blot, and RNA-sequencing (n=6) were performed. RESULTS GBM cells at higher stiff surfaces displayed increased migration (p<0.01). Magnetic Resonance Elastography (n=9) confirmed GBM heterogeneity, with softer cores (0.2-0.5kPa) transitioning to stiffer margins (>2kPa). Organotypic scaffolds mimicking the tumor margin promoted enhanced migration of GFP-labeled GBM cells. Margin cells and cells conditioned on stiffer matrices demonstrated higher migratory potential compared to their core and soft-conditioned counterparts (p<0.05). These trends persist over time and across multiple passages. ECM stiffness exposure causes nanomechanical, cytoskeletal, and transcriptional adaptations. AFM shows that cells from the margin (n=6) are softer (p<0.01) and more deformable (p<0.01). RNAseq identified upregulation of MM-specific genes in GBM margin cells. The mechanosensitive transcription factor YAP1 emerged as a key driver of MM. Notably, YAP1 is targetable by the FDA-approved drug Verteporfin. YAP1 expression was elevated in margin and stiff-conditioned cells (p<0.05). Stiff-conditioned cells and margin cells displayed increased sensitivity to Verteporfin treatment (p<0.05), with treatment significantly reducing margin migration compared to the core (p<0.001). CONCLUSION Stiffer GBM margins promote the formation of aggressive phenotypes through MM. Targeting YAP1 with Verteporfin, a readily available drug, presents a promising strategy to halt GBM invasion and improve patient outcomes.

TMIC-22. MYELOID CELLS INDUCE IMMUNOSUPPRESSION AND HYPORESPONSIVENESS OF CHIMERIC ANTIGEN RECEPTOR T CELLS IN THE NEURONAL MICROENVIRONMENT OF GLIOBLASTOMA

Abstract Chimeric antigen receptor T (CART)-cell therapy has shown little success in glioblastoma and other solid tumors, unlike in hematologic malignancies. The immunosuppressive milieu and intratumoral heterogeneity are some obstacles limiting the effectiveness of CART-cell therapy in glioblastoma. Our goal in this study is to investigate how CART-cell transcriptional adaptations and cellular interactions in the glioblastoma microenvironment drive T-cell hyporesponsiveness. Using a human neocortical slice model, we performed physically interacting cell sequencing (PIC-seq) to investigate the transcriptional diversity of NKG2D CAR-T cells and control T-cells in the glioblastoma ecosystem. Through the integration of spatially resolved T-cell receptors sequencing and PIC-seq, we investigated the cellular and receptor-ligand interactions in CART-treated glioblastoma cells. Gene regulatory networks were reconstructed to identify key transcriptional regulators of CD8 CART-cell dysfunction. We used in-silico perturbation to understand transcriptional drivers that could potentially enhance the cytotoxic response of CAR-Ts. Primary patient-derived glioma cells were inoculated into human neocortical slices from surgical access tissue specimens. After 3 days of tumor growth, slices were treated with NKG2D CARTs or control T-cells. We observed that after 2 days of anti-tumor T-cell response in the CART group, tumor growth did not differ between conditions. After 7 days of treatment, PIC-seq was used to profile CARTs (tomato+) and tumor cells (GFP+) as well as physically connected cells (+/+). We identified a transcriptional shift towards T-cell exhaustion in CART cells compared to T-cells lacking the NKG2D construct. This transcriptional phenotype was driven by myeloid-CART and myeloid-tumor interaction. We identified that tumor-associated macrophages with enhanced phagocytic function in proximity to mesenchymal-like tumor cells located within hypoxic niches are the main drivers of T-cell dysfunction. Gene regulatory network analysis demonstrated that MAF and BACH2 are key transcriptional regulators of CART-cell function. Myeloid-tumor and myeloid-T-cell interactions promote an immunosuppressive microenvironment and hyporesponsive T-cells in glioblastoma. Our data has the potential to catalyze the re-engineering and optimization of CART-cell therapies with sustained cytotoxic effect against glioblastoma.

Impact of the elderly lung mucosa on Mycobacterium tuberculosis transcriptional adaptation during infection of alveolar epithelial cells.

Tuberculosis is one of the leading causes of death due to a single infectious agent. Upon infection, Mycobacterium tuberculosis (M.tb) is deposited in the alveoli and encounters the lung mucosa or alveolar lining fluid (ALF). We previously showed that, as we age, ALF presents a higher degree of oxidation and inflammatory mediators, which favors M.tb replication in human macrophages and alveolar epithelial cells (ATs). Here, we define the transcriptional profile of M.tb when exposed to healthy ALF from adult (A-ALF) or elderly (E-ALF) humans before and during infection of ATs. Prior to infection, M.tb exposure to E-ALF upregulated genes essential for bacterial host adaptation directly involved in M.tb pathogenesis. During infection of ATs, E-ALF exposed M.tb further upregulated genes involved in its ability to escape into the AT cytosol bypassing critical host defense mechanisms, as well as genes associated with defense against oxidative stress. These findings demonstrate how alterations in human ALF during the aging process can impact the metabolic status of M.tb, potentially enabling a greater adaptation and survival within host cells. Importantly, we present the first transcriptomic analysis on the impact of the elderly lung mucosa on M.tb pathogenesis during intracellular replication in ATs.IMPORTANCETuberculosis is one of the leading causes of death due to a single infectious agent. Upon infection, Mycobacterium tuberculosis (M.tb) is deposited in the alveoli and comes in contact with the alveolar lining fluid (ALF). We previously showed that elderly ALF favors M.tb replication in human macrophages and alveolar epithelial cells (ATs). Here we define the transcriptional profile of when exposed to healthy ALF from adult (A-ALF) or elderly (E-ALF) humans before and during infection of ATs. Prior to infection, exposure to E-ALF upregulates genes essential for bacterial host adaptation and pathogenesis. During infection of ATs, E-ALF further upregulates M.tb genes involved in its ability to escape into the AT cytosol, as well as genes for defense against oxidative stress. These findings demonstrate how alterations in human ALF during the aging process can impact the metabolic status of M.tb, potentially enabling a greater adaptation and survival within host cells.

Cocaine taking and craving produce distinct transcriptional profiles in dopamine neurons.

Dopamine (DA) signaling plays an essential role in reward valence attribution and in encoding the reinforcing properties of natural and artificial rewards. The adaptive responses from midbrain dopamine neurons to artificial rewards such as drugs of abuse are therefore important for understanding the development of substance use disorders. Drug-induced changes in gene expression are one such adaptation that can determine the activity of dopamine signaling in projection regions of the brain reward system. One of the major challenges to obtaining this understanding involves the complex cellular makeup of the brain, where each neuron population can be defined by a distinct transcriptional profile. To bridge this gap, we have adapted a virus-based method for labeling and capture of dopamine nuclei, coupled with nuclear RNA-sequencing, to study the transcriptional adaptations, specifically, of dopamine neurons in the ventral tegmental area (VTA) during cocaine taking and cocaine craving, using a mouse model of cocaine intravenous self-administration (IVSA). Our results show significant changes in gene expression across non-drug operant training, cocaine taking, and cocaine craving, highlighted by an enrichment of repressive epigenetic modifying enzyme gene expression during cocaine craving. Immunohistochemical validation further revealed an increase of H3K9me3 deposition in DA neurons during cocaine craving. These results demonstrate that cocaine-induced transcriptional adaptations in dopamine neurons vary by phase of self-administration and underscore the utility of this approach for identifying relevant phase-specific molecular targets to study the behavioral course of substance use disorders.

Neutrophils in nasal polyps exhibit transcriptional adaptation and proinflammatory role depend on local polyp milieu.

Chronic rhinosinusitis with nasal polyps (CRSwNP) is an inflammatory upper airway disease, divided into eosinophilic CRSwNP (eCRSwNP) and noneosinophilic CRSwNP (neCRSwNP) according to eosinophilic levels. Neutrophils are major effector cells in CRSwNP. but their role in different inflammatory environments remain largely unclear. We performed an integrated transcriptome analysis of polyp-infiltrating neutrophils from CRSwNP patients, using healthy donor blood as a control. Flow cytometry and in vitro studies showed that neutrophils are activated in both CRSwNP type. The scRNA-sequencing analysis demonstrated that neutrophils were classified into five functional subsets, with GBP5+ neutrophils occurring mainly in neCRSwNPs and a high proportion of CXCL8+ neutrophils in both subendotypes. GBP5+ neutrophils exhibited significant IFN-I pathway activity in neCRSwNPs. CXCL8+ neutrophils displayed increased neutrophil activation scores and mainly secrete Oncostatin M (OSM), which facilitates communication with other cells. In vitro experiments revealed that OSM could enhance IL-13- or IL-17-mediated immune responses in nasal epithelial cells and fibroblasts. Our findings revealed that neutrophils exhibited transcriptional plasticity and activation when exposed to polyp tissue and exert their proinflammatory role in the pathogenesis of CRSwNP by releasing OSM to interact with epithelial cells and fibroblasts in a manner dependent on the inflammatory milieu.

Immunoresponsive Gene 1 Facilitates TLR4-agonist-Induced Augmentation of Innate Antimicrobial Immunity.

Treatment with the toll-like receptor (TLR) 4 agonist monophosphoryl lipid A (MPLA) conditions innate immunocytes to respond robustly to subsequent infection, a phenotype termed innate immune memory. Our published studies show that metabolic reprogramming of macrophages is a prominent feature of the memory phenotype. We undertook studies to define the functional contributions of tricarboxylic acid (TCA) cycle reprogramming to innate immune memory. We observed that priming of wild type (WT) mice with MPLA potently facilitated accumulation of the TCA cycle metabolite itaconate at sites of infection and enhanced microbial clearance. Augmentation of itaconate accumulation and microbial clearance was ablated in immuneresponsive gene 1 (Irg1) -deficient mice. We further observed that MPLA potently induces expression of Irg1 and accumulation of itaconate in macrophages. Compared to WT macrophages, the ability of Irg1-deficient macrophages to kill Pseudomonas aeruginosa was impaired. We further observed that itaconate is directly antimicrobial against P. aeruginosa at pH 5, which is characteristic of the phagolysosome, and is facilitated by reactive oxygen species. MPLA-induced augmentation of glycolysis, oxidative phosphorylation and accumulation of the TCA cycle metabolites succinate and malate was decreased in Irg1 KO macrophages compared to WT controls. RNA sequencing revealed suppressed transcription of genes associated with phagolysosome function and increased expression of genes associated with cytokine production and chemotaxis in Irg1 deficient macrophages. This study identifies a contribution of itaconate to MPLA-induced augmentation of innate antimicrobial immunity via facilitation of microbial killing as well as impact on metabolic and transcriptional adaptations.

Abstract C018: Targeting KRAS inhibitor resistance by ferroptosis induction in PDAC

Abstract KRAS inhibitors (KRASi) have shown promising antitumor activity in PDAC patients; however, intrinsic and acquired resistance will likely preclude durable monotherapy responses in the clinic. Genetic and pharmacologic targeting of the KRAS pathway is associated with transcriptional and metabolic adaptations, including epithelial-to-mesenchymal transition (EMT) and dysregulated iron homeostasis that contribute to resistance. Therapy resistant mesenchymal cancer cells have a high reliance on iron that increases their susceptibility to ferroptosis, an iron-dependent form of cell death characterized by generation of lipid reactive oxygen species (ROS). We hypothesize that KRASi resistance characterized by EMT can be overcome by ferroptosis induction in PDAC. Here, we use in vitro and ex vivo PDAC models and transcriptomic and proteomic profiling to define critical ferroptotic liabilities in KRASi-treated PDAC. First, we evaluated the effect of inhibiting different ferroptosis defense mechanisms (e.g. SLC7A11, GPX4, FSP1) using newly developed ferroptosis inducers (FINs) in a panel of epithelial and mesenchymal PDAC cells. Mesenchymal cells were more sensitive to ferroptosis induction compared to epithelial cells. Since ferroptotic death is regulated by multiple ferroptosis defense pathways, and in general PDAC are relatively resistant to ferroptosis induction, it is unlikely that a single ferroptosis inducer will be successful. Therefore, targeting multiple nodes with FINs in a combinatorial approach may lead to a more pronounced effect. Combination GPX4 and SLC7A11 inhibition and GPX4 and FSP1 inhibition were the most effective combinations compared to respective monotherapies. Next, we evaluated if combination of KRASi and ferroptosis induction is synergistic. Transcriptomic analysis of KRAS-mutant PDAC cell lines sensitive or with a priori resistance to KRASi revealed that EMT-related gene sets were upregulated in KRASi resistant cell lines. Quantitative proteomics revealed that PDAC cells treated long-term with KRASi undergo EMT and likewise long-term KRASi in in vivo PDAC models induces EMT. A KRASi-anchored CRISPR/Cas9 screen showed that sgRNAs targeting GPX4 dropped out in multiple cell lines indicating that ferroptosis induction may be synergistic upon KRAS inhibition. Remarkably, KRASi induced an increase of lipid peroxidation levels in PDAC cells which was rescued with antioxidant ferrostatin-1. We also observed that previously generated KRASi resistant cells were more sensitive to GPX4 inhibition. Combination of GPX4 inhibitor with KRASi led to a synergistic effect on cell viability and a significantly higher level of lipid peroxidation compared to monotherapy. These findings highlight the importance of investigating ferroptosis induction as a novel strategy to overcome therapeutic resistance in KRAS-mutant PDAC. We aim to further evaluate these combinations in vivo to determine the best combinatorial strategies in PDAC to prevent or circumvent KRASi resistance. Citation Format: Aparna Padhye, Qijia Yu, Nicole Sindoni, Huan Zhang, Julien Dilly, Hanrong Feng, Andrew Aguirre, Joseph Mancias. Targeting KRAS inhibitor resistance by ferroptosis induction in PDAC [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: Advances in Pancreatic Cancer Research; 2024 Sep 15-18; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2024;84(17 Suppl_2):Abstract nr C018.

Malonate is relevant to the lung environment and induces genome-wide stress responses in Pseudomonas aeruginosa.

Versatility in carbon source utilization is a major contributor to niche adaptation in Pseudomonas aeruginosa. Malonate is among the abundant carbon sources in the lung airways, yet it is understudied. Recently, we characterized how malonate impacts quorum sensing regulation, antibiotic resistance, and virulence factor production in P. aeruginosa. Herein, we show that malonate as a carbon source supports more robust growth in comparison to glycerol in several cystic fibrosis isolates of P. aeruginosa. Furthermore, we show phenotypic responses to malonate were conserved among clinical strains, i.e., formation of biomineralized biofilm-like aggregates, increased tolerance to kanamycin, and increased susceptibility to norfloxacin. Moreover, we explored transcriptional adaptations of P. aeruginosa UCBPP-PA14 (PA14) in response to malonate versus glycerol as a sole carbon source using transcriptomics. Malonate utilization activated glyoxylate and methylcitrate cycles and induced several stress responses, including oxidative, anaerobic, and metal stress responses associated with increases in intracellular aluminum and strontium. We identified several genes that were required for optimal growth of P. aeruginosa in malonate. Our findings reveal important remodeling of P. aeruginosa gene expression during its growth on malonate as a sole carbon source that is accompanied by several important phenotypic changes. These findings add to the accumulating literature highlighting the role of different carbon sources in the physiology of P. aeruginosa and its niche adaptation.

Sex-Specific Skeletal Muscle Gene Expression Responses to Exercise Reveal Novel Direct Mediators of Insulin Sensitivity Change.

Understanding the causal pathways, systems, and mechanisms through which exercise impacts human health is complex. This study explores molecular signaling related to whole-body insulin sensitivity (Si) by examining changes in skeletal muscle gene expression. The analysis considers differences by biological sex, exercise amount, and exercise intensity to identify potential molecular targets for developing pharmacologic agents that replicate the health benefits of exercise. The study involved 53 participants from the STRRIDE I and II trials who completed eight months of aerobic training. Skeletal muscle gene expression was measured using Affymetrix and Illumina technologies, while pre- and post-training Si was assessed via an intravenous glucose tolerance test. A novel gene discovery protocol, integrating three literature-derived and data-driven modeling strategies, was employed to identify causal pathways and direct causal factors based on differentially expressed transcripts associated with exercise intensity and amount. In women, the transcription factor targets identified were primarily influenced by exercise amount and were generally inhibitory. In contrast, in men, these targets were driven by exercise intensity and were generally activating. Transcription factors such as ATF1, CEBPA, BACH2, and STAT1 were commonly activating in both sexes. Specific transcriptional targets related to exercise-induced Si improvements included TACR3 and TMC7 for intensity-driven effects, and GRIN3B and EIF3B for amount-driven effects. Two key signaling pathways mediating aerobic exercise-induced Si improvements were identified: one centered on estrogen signaling and the other on phorbol ester (PKC) signaling, both converging on the epidermal growth factor receptor (EGFR) and other relevant targets. The signaling pathways mediating Si improvements from aerobic exercise differed by sex and were further distinguished by exercise intensity and amount. Transcriptional adaptations in skeletal muscle related to Si improvements appear to be causally linked to estrogen and PKC signaling, with EGFR and other identified targets emerging as potential skeletal muscle-specific drug targets to mimic the beneficial effects of exercise on Si.

Human Pluripotent Stem Cell-Derived Astrocyte Functionality Compares Favorably with Primary Rat Astrocytes.

Astrocytes are essential for the formation and maintenance of neural networks. However, a major technical challenge for investigating astrocyte function and disease-related pathophysiology has been the limited ability to obtain functional human astrocytes. Despite recent advances in human pluripotent stem cell (hPSC) techniques, primary rodent astrocytes remain the gold standard in coculture with human neurons. We demonstrate that a combination of leukemia inhibitory factor (LIF) and bone morphogenetic protein-4 (BMP4) directs hPSC-derived neural precursor cells to a highly pure population of astroglia in 28 d. Using single-cell RNA sequencing, we confirm the astroglial identity of these cells and highlight profound transcriptional adaptations in cocultured hPSC-derived astrocytes and neurons, consistent with their further maturation. In coculture with human neurons, multielectrode array recordings revealed robust network activity of human neurons in a coculture with hPSC-derived or rat astrocytes [3.63 ± 0.44 min-1 (hPSC-derived), 2.86 ± 0.64 min-1 (rat); p = 0.19]. In comparison, we found increased spike frequency within network bursts of human neurons cocultured with hPSC-derived astrocytes [56.31 ± 8.56 Hz (hPSC-derived), 24.77 ± 4.04 Hz (rat); p < 0.01], and whole-cell patch-clamp recordings revealed an increase of postsynaptic currents [2.76 ± 0.39 Hz (hPSC-derived), 1.07 ± 0.14 Hz (rat); p < 0.001], consistent with a corresponding increase in synapse density [14.90 ± 1.27/100 μm2 (hPSC-derived), 8.39 ± 0.63/100 μm2 (rat); p < 0.001]. Taken together, we show that hPSC-derived astrocytes compare favorably with rat astrocytes in supporting human neural network activity and maturation, providing a fully human platform for investigating astrocyte function and neuronal-glial interactions.

At the crossroads: strigolactones mediate changes in cytokinin synthesis and signalling in response to nitrogen limitation.

Strigolactones (SLs) are key regulators of shoot growth and responses to environmental stimuli. Numerous studies have indicated that nitrogen (N) limitation induces SL biosynthesis, suggesting that SLs may play a pivotal role in coordinating systemic responses to N availability, but this idea has not been clearly demonstrated. Here, we generated triple knockout mutants in the SL synthesis gene TaDWARF17 (TaD17) in bread wheat and investigated their phenotypic and transcriptional responses under N limitation, aiming to elucidate the role of SLs in the adaptation to N limitation. Tad17 mutants display typical SL mutant phenotypes, and fail to adapt their shoot growth appropriately to N. Despite exhibiting an increased tillering phenotype, Tad17 mutants continued to respond to N limitation by reducing tiller number, suggesting that SLs are not the sole regulators of tillering in response to N availability. RNA-seq analysis of basal nodes revealed that the loss of D17 significantly altered the transcriptional response of N-responsive genes, including changes in the expression profiles of key N response master regulators. Crucially, our findings suggest that SLs are required for the transcriptional downregulation of cytokinin (CK) synthesis and signalling in response to N limitation. Collectively, our results suggest that SLs are essential for the appropriate morphological and transcriptional adaptation to N limitation in wheat, and that the repressive effect of SLs on shoot growth is partly mediated by their repression of CK synthesis.

Prevalence of and gene regulatory constraints on transcriptional adaptation in single cells

BackgroundCells and tissues have a remarkable ability to adapt to genetic perturbations via a variety of molecular mechanisms. Nonsense-induced transcriptional compensation, a form of transcriptional adaptation, has recently emerged as one such mechanism, in which nonsense mutations in a gene trigger upregulation of related genes, possibly conferring robustness at cellular and organismal levels. However, beyond a handful of developmental contexts and curated sets of genes, no comprehensive genome-wide investigation of this behavior has been undertaken for mammalian cell types and conditions. How the regulatory-level effects of inherently stochastic compensatory gene networks contribute to phenotypic penetrance in single cells remains unclear.ResultsWe analyze existing bulk and single-cell transcriptomic datasets to uncover the prevalence of transcriptional adaptation in mammalian systems across diverse contexts and cell types. We perform regulon gene expression analyses of transcription factor target sets in both bulk and pooled single-cell genetic perturbation datasets. Our results reveal greater robustness in expression of regulons of transcription factors exhibiting transcriptional adaptation compared to those of transcription factors that do not. Stochastic mathematical modeling of minimal compensatory gene networks qualitatively recapitulates several aspects of transcriptional adaptation, including paralog upregulation and robustness to mutation. Combined with machine learning analysis of network features of interest, our framework offers potential explanations for which regulatory steps are most important for transcriptional adaptation.ConclusionsOur integrative approach identifies several putative hits—genes demonstrating possible transcriptional adaptation—to follow-up on experimentally and provides a formal quantitative framework to test and refine models of transcriptional adaptation.

Transitions in chromatin conformation shaped by fatty acids and the circadian clock underlie hepatic transcriptional reorganization in obese mice

The circadian clock system coordinates metabolic, physiological, and behavioral functions across a 24-h cycle, crucial for adapting to environmental changes. Disruptions in circadian rhythms contribute to major metabolic pathologies like obesity and Type 2 diabetes. Understanding the regulatory mechanisms governing circadian control is vital for identifying therapeutic targets. It is well characterized that chromatin remodeling and 3D structure at genome regulatory elements contributes to circadian transcriptional cycles; yet the impact of rhythmic chromatin topology in metabolic disease is largely unexplored. In this study, we explore how the spatial configuration of the genome adapts to diet, rewiring circadian transcription and contributing to dysfunctional metabolism. We describe daily fluctuations in chromatin contacts between distal regulatory elements of metabolic control genes in livers from lean and obese mice and identify specific lipid-responsive regions recruiting the clock molecular machinery. Interestingly, under high-fat feeding, a distinct interactome for the clock-controlled gene Dbp strategically promotes the expression of distal metabolic genes including Fgf21. Alongside, new chromatin loops between regulatory elements from genes involved in lipid metabolism control contribute to their transcriptional activation. These enhancers are responsive to lipids through CEBPβ, counteracting the circadian repressor REVERBa. Our findings highlight the intricate coupling of circadian gene expression to a dynamic nuclear environment under high-fat feeding, supporting a temporally regulated program of gene expression and transcriptional adaptation to diet.

3D Brain Vascular Niche Model Captures Invasive Behavior and Gene Signatures of Glioblastoma.

Glioblastoma (GBM) is a lethal brain cancer with no effective treatment; understanding how GBM cells respond to tumor microenvironment remains challenging as conventional cell cultures lack proper cytoarchitecture while in vivo animal models present complexity all at once. Developing a culture system to bridge the gap is thus crucial. Here, we employed a multicellular approach using human glia and vascular cells to optimize a 3-dimensional (3D) brain vascular niche model that enabled not only long-term culture of patient derived GBM cells but also recapitulation of key features of GBM heterogeneity, in particular invasion behavior and vascular association. Comparative transcriptomics of identical patient derived GBM cells in 3D and in vivo xenotransplants models revealed that glia-vascular contact induced genes concerning neural/glia development, synaptic regulation, as well as immune suppression. This gene signature displayed region specific enrichment in the leading edge and microvascular proliferation zones in human GBM and predicted poor prognosis. Gene variance analysis also uncovered histone demethylation and xylosyltransferase activity as main themes for gene adaption of GBM cells in vivo . Furthermore, our 3D model also demonstrated the capacity to provide a quiescence and a protective niche against chemotherapy. In summary, an advanced 3D brain vascular model can bridge the gap between 2D cultures and in vivo models in capturing key features of GBM heterogeneity and unveil previously unrecognized influence of glia-vascular contact for transcriptional adaption in GBM cells featuring neural/synaptic interaction and immunosuppression.

Scanpro is a tool for robust proportion analysis of single-cell resolution data

In higher organisms, individual cells respond to signals and perturbations by epigenetic regulation and transcriptional adaptation. However, in addition to shifting the expression level of individual genes, the adaptive response of cells can also lead to shifts in the proportions of different cell types. Recent methods such as scRNA-seq allow for the interrogation of expression on the single-cell level, and can quantify individual cell type clusters within complex tissue samples. In order to identify clusters showing differential composition between different biological conditions, differential proportion analysis has recently been introduced. However, bioinformatics tools for robust proportion analysis of both replicated and unreplicated single-cell datasets are critically missing. In this manuscript, we present Scanpro, a modular tool for proportion analysis, seamlessly integrating into widely accepted frameworks in the Python environment. Scanpro is fast, accurate, supports datasets without replicates, and is intended to be used by bioinformatics experts and beginners alike.

Skeletal muscle BMAL1 is necessary for transcriptional adaptation of local and peripheral tissues in response to endurance exercise training

ObjectiveIn this investigation, we addressed the contribution of the core circadian clock factor, BMAL1, in skeletal muscle to both acute transcriptional responses to exercise and transcriptional remodeling in response to exercise training. Additionally, we adopted a systems biology approach to investigate how loss of skeletal muscle BMAL1 altered peripheral tissue homeostasis as well as exercise training adaptations in iWAT, liver, heart, and lung of male mice. MethodsCombining inducible skeletal muscle specific BMAL1 knockout mice, physiological testing and standardized exercise protocols, we performed a multi-omic analysis (transcriptomics, chromatin accessibility and metabolomics) to explore loss of muscle BMAL1 on muscle and peripheral tissue responses to exercise. ResultsMuscle-specific BMAL1 knockout mice demonstrated a blunted transcriptional response to acute exercise, characterized by the lack of upregulation of well-established exercise responsive transcription factors including Nr4a3 and Ppargc1a. Six weeks of exercise training in muscle-specific BMAL1 knockout mice induced significantly greater and divergent transcriptomic and metabolomic changes in muscle. Surprisingly, liver, lung, inguinal white adipose and heart showed divergent exercise training transcriptomes with less than 5% of ‘exercise-training’ responsive genes shared for each tissue between genotypes. ConclusionsOur investigation has uncovered the critical role that BMAL1 plays in skeletal muscle as a key regulator of gene expression programs for both acute exercise and training adaptations. In addition, our work has uncovered the significant impact that altered exercise response in muscle and its likely impact on the system plays in the peripheral tissue adaptations to exercise training. Our work also demonstrates that if the muscle adaptations diverge to a more maladaptive state this is linked to increased gene expression signatures of inflammation across many tissues. Understanding the molecular targets and pathways contributing to health vs. maladaptive exercise adaptations will be critical for the next stage of therapeutic design for exercise mimetics.

Bacterial single-cell RNA sequencing captures biofilm transcriptional heterogeneity and differential responses to immune pressure.

Biofilm formation is an important mechanism of survival and persistence for many bacterial pathogens. These multicellular communities contain subpopulations of cells that display vast metabolic and transcriptional diversity along with high recalcitrance to antibiotics and host immune defenses. Investigating the complex heterogeneity within biofilm has been hindered by the lack of a sensitive and high-throughput method to assess stochastic transcriptional activity and regulation between bacterial subpopulations, which requires single-cell resolution. We have developed an optimized bacterial single-cell RNA sequencing method, BaSSSh-seq, to study Staphylococcus aureus diversity during biofilm growth and transcriptional adaptations following immune cell exposure. We validated the ability of BaSSSh-seq to capture extensive transcriptional heterogeneity during biofilm compared to planktonic growth. Application of new computational tools revealed transcriptional regulatory networks across the heterogeneous biofilm subpopulations and identification of gene sets that were associated with a trajectory from planktonic to biofilm growth. BaSSSh-seq also detected alterations in biofilm metabolism, stress response, and virulence that were tailored to distinct immune cell populations. This work provides an innovative platform to explore biofilm dynamics at single-cell resolution, unlocking the potential for identifying biofilm adaptations to environmental signals and immune pressure.

Identification of Therapy-Induced Clonal Evolution and Resistance Pathways in Minimal Residual Clones in Multiple Myeloma through Single-Cell Sequencing.

In multiple myeloma (MM), therapy-induced clonal evolution is associated with treatment resistance and is one of the most important hindrances toward a cure for MM. To further understand the molecular mechanisms controlling the clonal evolution of MM, we applied single-cell RNA sequencing (scRNA-seq) to paired diagnostic and posttreatment bone marrow (BM) samples. scRNA-seq was performed on 38 BM samples from patients with monoclonal gammopathy of undetermined significance (n = 1), MM patients at diagnosis (n = 19), MM posttreatment (n = 17), and one healthy donor (HD). The single-cell transcriptome data of malignant plasma cells (PC) and the surrounding immune microenvironment were analyzed. Profiling by scRNA-seq data revealed three primary trajectories of transcriptional evolution after treatment: clonal elimination in patients with undetectable minimal residual disease (MRD-) and clonal stabilization and clonal selection in detectable MRD (MRD+) patients. We noted a metabolic shift toward fatty acid oxidation in cycling-resistant PCs, whereas selective PCs favored the NF-κB pathway. Intriguingly, when comparing the genetic and transcriptional dynamics, we found a significant correlation between genetic and nongenetic factors in driving the clonal evolution. Furthermore, we identified variations in cellular interactions between malignant PCs and the tumor microenvironment. Selective PCs showed the most robust cellular interactions with the tumor microenvironment. These data suggest that MM cells could rapidly adapt to induction treatment through transcriptional adaptation, metabolic adaptation, and specialized immune evasion. Targeting therapy-induced resistance mechanisms may help to avert refractory disease in MM.