Search

A cell's fate involves transitions among its various states, each defined by a distinct gene expression profile governed by the topology of gene regulatory networks, which are affected by 3D genome organization. Here, we develop thermodynamic models to determine the fate of a malignant cell as governed by the tumor suppressor p53 signaling network, taking into account long-range chromatin interactions in the mean-field approximation. The tumor suppressor p53 responds to stress by selectively triggering one of the potential transcription programs that influence many layers of cell signaling. These range from p53 phosphorylation to modulation of its DNA binding affinity, phase separation phenomena, and internal connectivity among cell fate genes. We use the minimum free energy of the system as a fundamental property of biological networks that influences the connection between the gene network topology and the state of the cell. We constructed models based on network topology and equilibrium thermodynamics. Our modeling shows that the binding of phosphorylated p53 to promoters of target genes can have properties of a first order phase transition. We apply our model to cancer cell lines ranging from breast cancer (MCF-7), colon cancer (HCT116), and leukemia (K562), with each one characterized by a specific network topology that determines the cell fate. Our results clarify the biological relevance of these mechanisms and suggest that they represent flexible network designs for switching between developmental decisions. [ABSTRACT FROM AUTHOR]

Spatial transcriptomics technology enables the mapping of gene expression within tissues, allowing researchers to visualize the spatial distribution of RNA molecules and gain insights into cellular organization, interactions, and functions in their native environments. A variety of spatial technologies are now commercially available, each offering distinct technical parameters such as cellular resolution, detection sensitivity, gene coverage, and throughput. This wide range of options can make it challenges or create confusion for researchers to select the most appropriate platform for their specific research objectives. In this paper, we will analyze and compare seven major commercially available spatial platforms to guide researchers in choosing the most suitable option for their needs.

Seeds are complex structures composed of three regions, embryo, endosperm, and seed coat, with each further divided into subregions that consist of tissues, cell layers, and cell types. Although the seed is well characterized anatomically, much less is known about the genetic circuitry that dictates its spatial complexity. To address this issue, we profiled mRNAs from anatomically distinct seed subregions at several developmental stages. Analyses of these profiles showed that all subregions express similar diverse gene numbers and that the small gene numbers expressed subregion specifically provide information about the biological processes that occur in these seed compartments. In parallel, we profiled RNAs in individual nuclei and identified nuclei clusters representing distinct cell identities. Integrating single-nucleus RNA and subregion mRNA transcriptomes allowed most cell identities to be assigned to specific subregions and cell types and/or cell states. The number of cell identities exceeds the number of anatomically distinguishable cell types, emphasizing the spatial complexity of seeds. We defined gene coexpression networks that underlie distinct biological processes during seed development. We showed that network distribution among subregions and cell identities is highly variable. Some networks operate in single subregions and/or cell identities, and many coexpression networks operate in multiple subregions and/or cell identities. We also showed that single subregions and cell identities possess several networks. Together, our studies provide unique insights into the biological processes and genetic circuitry that underlie the spatial landscape of the seed.

Coronavirus disease 2019 (COVID-19) poses significant risks for solid organ transplant recipients, who have atypical but poorly characterized immune responses to infection. We aim to understand the host immunologic and microbial features of COVID-19 in transplant recipients by leveraging a prospective multicenter cohort of 86 transplant recipients age- and sex-matched with 172 non-transplant controls. We find that transplant recipients have higher nasal SARS-CoV-2 viral abundance and impaired viral clearance, and lower anti-spike IgG levels. In addition, transplant recipients exhibit decreased plasmablasts and transitional B cells, and increased senescent T cells. Blood and nasal transcriptional profiling demonstrate unexpected upregulation of innate immune signaling pathways and increased levels of several proinflammatory serum chemokines. Severe disease in transplant recipients, however, is characterized by a less robust induction of pro-inflammatory genes and chemokines. Together, our study reveals distinct immune features and altered viral dynamics in solid organ transplant recipients.

Vertebrates differ over 100,000-fold in responses to pro-inflammatory agonists such as bacterial lipopolysaccharide (LPS), complicating use of animal models to study human sepsis or inflammatory disorders. We compared transcriptomes of resting and LPS-exposed blood from six LPS-sensitive species (rabbit, pig, sheep, cow, chimpanzee, human) and four LPS-resilient species (mice, rats, baboon, rhesus), as well as plasma proteomes and lipidomes. Unexpectedly, at baseline, sensitive species already had enhanced expression of LPS-responsive genes relative to resilient species. After LPS stimulation, maximally different genes in resilient species included genes that detoxify LPS, diminish bacterial growth, discriminate sepsis from SIRS, and play roles in autophagy and apoptosis. The findings reveal the molecular landscape of species differences in inflammation. This may inform better selection of species for pre-clinical models and could lead to new therapeutic strategies that mimic mechanisms in inflammation-resilient species to limit inflammation without causing immunosuppression.

Transcriptomics is a powerful approach for functional genomics and systems biology, yet it can also be used for genetic part discovery. Here, we derive constitutive and light-regulated promoters directly from transcriptomics data of the basidiomycete red yeast Xanthophyllomyces dendrorhous CBS 6938 (anamorph Phaffia rhodozyma) and use these promoters with other genetic elements to create a modular synthetic biology parts collection for this organism. X. dendrorhous is currently the sole biotechnologically relevant yeast in the Tremellomycete class-it produces large amounts of astaxanthin, especially under oxidative stress and exposure to light. Thus, we performed transcriptomics on X. dendrorhous under different wavelengths of light (red, green, blue, and ultraviolet) and oxidative stress. Differential gene expression analysis (DGE) revealed that terpenoid biosynthesis was primarily upregulated by light through crtI, while oxidative stress upregulated several genes in the pathway. Further gene ontology (GO) analysis revealed a complex survival response to ultraviolet (UV) where X. dendrorhous upregulates aromatic amino acid and tetraterpenoid biosynthesis and downregulates central carbon metabolism and respiration. The DGE data was also used to identify 26 constitutive and regulated genes, and then, putative promoters for each of the 26 genes were derived from the genome. Simultaneously, a modular cloning system for X. dendrorhous was developed, including integration sites, terminators, selection markers, and reporters. Each of the 26 putative promoters were integrated into the genome and characterized by luciferase assay in the dark and under UV light. The putative constitutive promoters were constitutive in the synthetic genetic context, but so were many of the putative regulated promoters. Notably, one putative promoter, derived from a hypothetical gene, showed ninefold activation upon UV exposure. Thus, this study reveals metabolic pathway regulation and develops a genetic parts collection for X. dendrorhous from transcriptomic data. Therefore, this study demonstrates that combining systems biology and synthetic biology into an omics-to-parts workflow can simultaneously provide useful biological insight and genetic tools for nonconventional microbes, particularly those without a related model organism. This approach can enhance current efforts to engineer diverse microbes. KEY POINTS: • Transcriptomics revealed further insights into the photobiology of X. dendrorhous, specifically metabolic nodes that are transcriptionally regulated by light. • A modular genetic part collection was developed, including 26 constitutive and regulated promoters derived from the transcriptomics of X. dendrorhous. • Omics-to-parts can be applied to nonconventional microbes for rapid onboarding.

Background and objectiveWe characterized tumor prostate-specific membrane antigen (PSMA) levels as a reflection of cancer biology and treatment sensitivities for treatment-naïve prostate cancer.MethodsWe first correlated PSMA positron emission tomography (PET) maximum standardized uptake values (SUVmax) in primary prostate cancer with tumor FOLH1 (PSMA RNA abundance) to establish RNA as a proxy (n = 55). We then discovered and validated molecular pathways associated with PSMA RNA levels in two large primary tumor cohorts. We validated those associations in independent cohorts (18 total; 5684 tumor samples) to characterize the pathways and treatment responses associated with PSMA.Key findings and limitationsPSMA RNA abundance correlates moderately with SUVmax (ρ = 0.41). In independent cohorts, androgen receptor signaling is more active in tumors with high PSMA. Accordingly, patients with high PSMA tumors experienced longer cancer-specific survival when managed with androgen deprivation therapy for biochemical recurrence (adjusted hazard ratio [AHR] 0.54 [0.34-0.87]; n = 174). PSMA low tumors possess molecular markers of resistance to radiotherapy. Consistent with this, patients with high PSMA tumors experience longer time to recurrence following primary radiotherapy (AHR 0.50 [0.28-0.90]; n = 248). In the SAKK09/10 trial (n = 224), patients with high PSMA tumors who were managed with salvage radiotherapy experienced longer time to progression in the 64-Gy arm (restricted mean survival time [RMST] +7.60 [0.05-15.16]), but this effect was mitigated in the 70-Gy arm (RMST 3.52 [-3.30 to 10.33]). Limitations include using PSMA RNA as a surrogate for PET SUVmax.Conclusions and clinical implicationsPSMA levels in treatment-naïve prostate cancer differentiate tumor biology and treatment susceptibilities. These results warrant validation using PET metrics to substantiate management decisions based on imaging.

The pathogenesis of Alzheimer's disease (AD) depends on environmental and heritable factors, with its molecular etiology still unclear. Here we present a spatial transcriptomic (ST) and single-nucleus transcriptomic survey of late-onset sporadic AD and AD in Down syndrome (DSAD). Studying DSAD provides an opportunity to enhance our understanding of the AD transcriptome, potentially bridging the gap between genetic mouse models and sporadic AD. We identified transcriptomic changes that may underlie cortical layer-preferential pathology accumulation. Spatial co-expression network analyses revealed transient and regionally restricted disease processes, including a glial inflammatory program dysregulated in upper cortical layers and implicated in AD genetic risk and amyloid-associated processes. Cell-cell communication analysis further contextualized this gene program in dysregulated signaling networks. Finally, we generated ST data from an amyloid AD mouse model to identify cross-species amyloid-proximal transcriptomic changes with conformational context.

BACKGROUND: Despite nearly 100% 5-year survival for localised prostate cancer, the survival rate for metastatic prostate cancer significantly declines to 32%. Thus, it is crucial to identify molecular indicators that reflect the progression from localised disease to metastatic prostate cancer. METHODS: To search for molecular indicators associated with prostate cancer metastasis, we performed proteomic analysis of rapid autopsy tissue samples from metastatic prostate cancer (N = 8) and localised prostate cancer (N = 2). Then, we utilised multiple independent, publicly available prostate cancer patient datasets to select candidates that also correlate with worse prostate cancer clinical prognosis. RESULTS: We identified 154 proteins with increased expressions in metastases relative to localised prostate cancer through proteomic analysis. From the subset of these candidates that correlate with prostate cancer recurrence (N = 28) and shorter disease-free survival (N = 37), we identified a 5-gene signature panel with improved performance in predicting worse clinical prognosis relative to individual candidates. CONCLUSIONS: Our study presents a new 5-gene signature panel that is associated with worse clinical prognosis and is elevated in prostate cancer metastasis on both protein and mRNA levels. Our 5-gene signature panel represents a potential modality for the prediction of prostate cancer progression towards the onset of metastasis.

BackgroundPhysical activity is an adjunctive therapy that improves symptoms in people living with systemic lupus erythematosus (SLE), yet the mechanisms underlying this benefit remain unclear.MethodsWe carried out a cohort study of 123 patients with SLE enrolled in the California Lupus Epidemiology Study (CLUES). The primary predictor variable was self-reported physical activity, which was measured using a previously validated instrument. We analyzed peripheral blood mononuclear cell (PBMC) single-cell RNA sequencing (scRNA-seq) data available from the cohort. From the scRNA-seq data, we compared immune cell frequencies, cell-specific gene expression, biological signalling pathways, and upstream cytokine activation states between physically active and inactive patients, adjusting for age, sex and race.FindingsWe found that physical activity influenced immune cell frequencies, with sedentary patients most notably demonstrating greater CD4+ T cell lymphopenia (Padj = 0.028). Differential gene expression analysis identified a transcriptional signature of physical inactivity across five cell types. In CD4+ and CD8+ T cells, this signature was characterized by 686 and 445 differentially expressed genes (Padj

Plant-microbe symbioses require intense interaction and genetic coordination to successfully establish in specific cell types of the host and symbiont. Traditional RNA-seq methodologies lack the cellular resolution to fully capture these complexities, but single-cell and spatial transcriptomics (ST) are now allowing scientists to probe symbiotic interactions at an unprecedented level of detail. Here, we discuss the advantages that novel spatial and single-cell transcriptomic technologies provide in studying plant-microbe endosymbioses and highlight key recent studies. Finally, we consider the remaining limitations of applying these approaches to symbiosis research, which are mainly related to the simultaneous capture of both plant and microbial transcripts within the same cells.

Resistance training is used to combat skeletal muscle function decline in older adults. Few studies have been designed specific for females, resulting in very limited treatment options for skeletal muscle atrophy in aging women. Here, we analyzed the gene expression profiles of skeletal muscle samples from sedentary young women, sedentary older women, and resistance-trained older women, using microarray data from public database. A total of 45 genes that were differentially expressed during female muscle aging and reversed by resistance training were identified. Functional and pathway enrichment analysis, protein–protein interaction network analysis, and receiver operating characteristic analysis were performed to reveal the key genes and pathways involved in the effects of resistance training on female muscle aging. The collagen genes COL1A1, COL3A1, and COL4A1 were identified important regulators of female muscle aging and resistance training, by modulating multiple signaling pathways, such as PI3 kinase-Akt signaling, focal adhesions, extracellular matrix-receptor interactions, and relaxin signaling. Interestingly, the expression of CDKN1A and TP63 were increased during aging, and further upregulated by resistance training in older women, suggesting they may negatively affect resistance training outcomes. Our findings provide novel insights into the molecular mechanisms of resistance training on female muscle aging and identify potential biomarkers and targets for clinical intervention. [ABSTRACT FROM AUTHOR]

Cancer transcriptional patterns reflect both unique features and shared hallmarks across diverse cancer types, but whether differences in these patterns are sufficient to characterize the full breadth of tumor phenotype heterogeneity remains an open question. We hypothesized that these shared transcriptomic signatures reflect repurposed versions of functional tasks performed by normal tissues. Starting with normal tissue transcriptomic profiles, we use non-negative matrix factorization to derive six distinct transcriptomic phenotypes, called archetypes, which combine to describe both normal tissue patterns and variations across a broad spectrum of malignancies. We show that differential enrichment of these signatures correlates with key tumor characteristics, including overall patient survival and drug sensitivity, independent of clinically actionable DNA alterations. Additionally, we show that in HR+/HER2- breast cancers, metastatic tumors adopt transcriptomic signatures consistent with the invaded tissue. Broadly, our findings suggest that cancer often arrogates normal tissue transcriptomic characteristics as a component of both malignant progression and drug response. This quantitative framework provides a strategy for connecting the diversity of cancer phenotypes and could potentially help manage individual patients.

There are myriad types of biomedical data-molecular, clinical images, and others. When a group of patients with the same underlying disease exhibits similarities across multiple types of data, this is called a subtype. Existing subtyping approaches struggle to handle diverse data types with missing information. To improve subtype discovery, we exploited changes in the correlation-structure between different data types to create iSubGen, an algorithm for integrative subtype generation. iSubGen can accommodate any feature that can be compared with a similarity metric to create subtypes versatilely. It can combine arbitrary data types for subtype discovery, such as merging genetic, transcriptomic, proteomic, and pathway data. iSubGen recapitulates known subtypes across multiple cancers even with substantial missing data and identifies subtypes with distinct clinical behaviors. It performs equally with or superior to other subtyping methods, offering greater stability and robustness to missing data and flexibility to new data types. It is available at https://cran.r-project.org/web/packages/iSubGen.

Effects from aging in single cells are heterogenous, whereas at the organ- and tissue-levels aging phenotypes tend to appear as stereotypical changes. The mammary epithelium is a bilayer of two major phenotypically and functionally distinct cell lineages: luminal epithelial and myoepithelial cells. Mammary luminal epithelia exhibit substantial stereotypical changes with age that merit attention because these cells are the putative cells-of-origin for breast cancers. We hypothesize that effects from aging that impinge upon maintenance of lineage fidelity increase susceptibility to cancer initiation. We generated and analyzed transcriptomes from primary luminal epithelial and myoepithelial cells from younger 55 y women. In addition to age-dependent directional changes in gene expression, we observed increased transcriptional variance with age that contributed to genome-wide loss of lineage fidelity. Age-dependent variant responses were common to both lineages, whereas directional changes were almost exclusively detected in luminal epithelia and involved altered regulation of chromatin and genome organizers such as SATB1. Epithelial expression variance of gap junction protein GJB6 increased with age, and modulation of GJB6 expression in heterochronous co-cultures revealed that it provided a communication conduit from myoepithelial cells that drove directional change in luminal cells. Age-dependent luminal transcriptomes comprised a prominent signal that could be detected in bulk tissue during aging and transition into cancers. A machine learning classifier based on luminal-specific aging distinguished normal from cancer tissue and was highly predictive of breast cancer subtype. We speculate that luminal epithelia are the ultimate site of integration of the variant responses to aging in their surrounding tissue, and that their emergent phenotype both endows cells with the ability to become cancer-cells-of-origin and represents a biosensor that presages cancer susceptibility.

Two correspondences raised concerns or comments about our analyses regarding exaggerated false positives found by differential expression (DE) methods. Here, we discuss the points they raise and explain why we agree or disagree with these points. We add new analysis to confirm that the Wilcoxon rank-sum test remains the most robust method compared to the other five DE methods (DESeq2, edgeR, limma-voom, dearseq, and NOISeq) in two-condition DE analyses after considering normalization and winsorization, the data preprocessing steps discussed in the two correspondences.

The dog, Canis lupus familiaris, is an important model for studying human diseases. Unlike many model organisms, the dog genome has a comparatively poor functional annotation, which hampers gene discovery for development, morphology, disease, and behavior. To fill this gap, we established a comprehensive tissue biobank for both the dog and wolf samples. The biobank consists of 5485 samples representing 132 tissues from 13 dogs, 12 dog embryos, and 24 wolves. In a subset of 100 tissues from nine dogs and 12 embryos, we characterized gene expression activity for each promoter, including alternative and novel, i.e., previously not annotated, promoter regions, using the 5 targeting RNA sequencing technology STRT2-seq. We identified over 100,000 promoter region candidates in the recent canine genome assembly, CanFam4, including over 45,000 highly reproducible sites with gene expression and respective tissue enrichment levels. We provide a promoter and gene expression atlas with interactive, open data resources, including a data coordination center and genome browser track hubs. We demonstrated the applicability of Dog Genome Annotation (DoGA) data and resources using multiple examples spanning canine embryonic development, morphology and behavior, and diseases across species.

Torpor encompasses diverse adaptations to extreme environmental stressors such as hibernation, aestivation, brumation, and daily torpor. Here we introduce StrokeofGenus, an analytic pipeline that identifies distinct transcriptomic states and shared gene expression patterns across studies, tissues, and species. We use StrokeofGenus to study multiple and diverse forms of torpor from publicly-available RNA-seq datasets that span eight species and two classes. We identify three transcriptionally distinct states during the cycle of heterothermia: euthermia, torpor, and interbout arousal. We also identify torpor-specific gene expression patterns that are shared both across tissues and between species with over three hundred million years of evolutionary divergence. We further demonstrate the general sharing of gene expression patterns in multiple forms of torpor, implying a common evolutionary origin for this process. Although here we apply StrokeofGenus to analysis of torpor, it can be used to interrogate any other complex physiological processes defined by transient transcriptomic states.

UNLABELLED: Microorganisms with simplified genomes represent interesting cell chassis for systems and synthetic biology. However, genome reduction can lead to undesired traits, such as decreased growth rate and metabolic imbalances. To investigate the impact of genome reduction on Escherichia coli strain DGF-298, a strain in which ~ 36% of the genome has been removed, we reconstructed a strain-specific metabolic model (iAC1061), investigated the regulation of gene expression using iModulon-based transcriptome analysis, and performed adaptive laboratory evolution to let the strain correct potential imbalances that arose during its simplification. The model notably predicted that the removal of all three key pathways for glycolaldehyde disposal in this microorganism would lead to a metabolic bottleneck through folate starvation. Glycolaldehyde is also known to cause self-generation of reactive oxygen species, as evidenced by the increased expression of oxidative stress resistance genes in the SoxS iModulon. The reintroduction of the aldA gene, responsible for one native glycolaldehyde disposal route, alleviated the constitutive oxidative stress response. Our results suggest that systems-level approaches and adaptive laboratory evolution have additive benefits when trying to repair and optimize genome-engineered strains. IMPORTANCE: Genomic streamlining can be employed in model organisms to reduce complexity and enhance strain predictability. One of the most striking examples is the bacterial strain Escherichia coli DGF-298, notable for having over one-third of its genome deleted. However, such extensive genome modifications raise the question of how similar this simplified cell remains when compared with its parent, and what are the possible unintended consequences of this simplification. In this study, we used metabolic modeling along with iModulon-based transcriptomic analysis in different growth conditions to assess the impact of genome reduction on metabolism and gene regulation. We observed little impact of genomic reduction on the regulatory network of E. coli DGF-298 and identified a potential metabolic bottleneck leading to the constitutive activity of the SoxS iModulon. We then leveraged the models predictions to successfully restore SoxS activity to the basal level.

From single-cell RNA-sequencing (scRNA-seq) and spatial transcriptomics (ST), one can extract high-dimensional gene expression patterns that can be described by intercellular communication networks or decoupled gene modules. These two descriptions of information flow are often assumed to occur independently. However, intercellular communication drives directed flows of information that are mediated by intracellular gene modules, in turn triggering outflows of other signals. Methodologies to describe such intercellular flows are lacking. We present FlowSig, a method that infers communication-driven intercellular flows from scRNA-seq or ST data using graphical causal modeling and conditional independence. We benchmark FlowSig using newly generated experimental cortical organoid data and synthetic data generated from mathematical modeling. We demonstrate FlowSig's utility by applying it to various studies, showing that FlowSig can capture stimulation-induced changes to paracrine signaling in pancreatic islets, demonstrate shifts in intercellular flows due to increasing COVID-19 severity and reconstruct morphogen-driven activator-inhibitor patterns in mouse embryogenesis.

This article focuses on screening the major secreted proteins by the ischemia-challenged cardiac stromal fibroblasts (CF), the assessment of their expression status and functional role in the post-ischemic left ventricle (LV) and in the ischemia-challenged CF culture and to phenotype CF at single cell resolution based on the positivity of the identified mediators. The expression level of CRSP2, HSP27, IL-8, Cofilin-1, and HSP90 in the LV tissues following coronary artery bypass graft (CABG) and myocardial infarction (MI) and CF cells followed the screening profile derived from the MS/MS findings. The histology data unveiled ECM disorganization, inflammation and fibrosis reflecting the ischemic pathology. CRSP2, HSP27, and HSP90 were significantly upregulated in the LV-CABG tissues with a concomitant reduction ion LV-MI whereas Cofilin-1, IL8, Nrf2, and Troponin I were downregulated in LV-CABG and increased in LV-MI. Similar trends were exhibited by ischemic CF. Single cell transcriptomics revealed multiple sub-phenotypes of CF based on their respective upregulation of CRSP2, HSP27, IL-8, Cofilin-1, HSP90, Troponin I and Nrf2 unveiling pathological and pro-healing phenotypes. Further investigations regarding the underlying signaling mechanisms and validation of sub-populations would offer novel translational avenues for the management of cardiac diseases.

Despite recent work linking mixed phenotype acute leukemia (MPAL) to certain genetic lesions, specific driver mutations remain undefined for a significant proportion of patients and no genetic subtype is predictive of clinical outcomes. Moreover, therapeutic strategy for MPAL remains unclear, and prognosis is overall poor. We performed multiomic single cell profiling of 14 newly diagnosed adult MPAL patients to characterize the inter- and intra-tumoral transcriptional, immunophenotypic, and genetic landscapes of MPAL. We show that neither genetic profile nor transcriptome reliably correlate with specific MPAL immunophenotypes. Despite this, we find that MPAL blasts express a shared stem cell-like transcriptional profile indicative of high differentiation potential. Patients with the highest differentiation potential demonstrate inferior survival in our dataset. A gene set score, MPAL95, derived from genes highly enriched in the most stem-like MPAL cells, is applicable to bulk RNA sequencing data and is predictive of survival in an independent patient cohort, suggesting a potential strategy for clinical risk stratification.

The rapid growth of large-scale spatial gene expression data demands efficient and reliable computational tools to extract major trends of gene expression in their native spatial context. Here, we used stability-driven unsupervised learning (i.e., staNMF) to identify principal patterns (PPs) of 3D gene expression profiles and understand spatial gene distribution and anatomical localization at the whole mouse brain level. Our subsequent spatial correlation analysis systematically compared the PPs to known anatomical regions and ontology from the Allen Mouse Brain Atlas using spatial neighborhoods. We demonstrate that our stable and spatially coherent PPs, whose linear combinations accurately approximate the spatial gene data, are highly correlated with combinations of expert-annotated brain regions. These PPs yield a brain ontology based purely on spatial gene expression. Our PP identification approach outperforms principal component analysis and typical clustering algorithms on the same task. Moreover, we show that the stable PPs reveal marked regional imbalance of brainwide genetic architecture, leading to region-specific marker genes and gene coexpression networks. Our findings highlight the advantages of stability-driven machine learning for plausible biological discovery from dense spatial gene expression data, streamlining tasks that are infeasible by conventional manual approaches.

Sterile inflammation after myocardial infarction is classically credited to myeloid cells interacting with dead cell debris in the infarct zone1,2. Here we show that cardiomyocytes are the dominant initiators of a previously undescribed type I interferon response in the infarct borderzone. Using spatial transcriptomics analysis in mice and humans, we find that myocardial infarction induces colonies of interferon-induced cells (IFNICs) expressing interferon-stimulated genes decorating the borderzone, where cardiomyocytes experience mechanical stress, nuclear rupture and escape of chromosomal DNA. Cardiomyocyte-selective deletion of Irf3 abrogated IFNIC colonies, whereas mice lacking Irf3 in fibroblasts, macrophages, neutrophils or endothelial cells, Ccr2-deficient mice or plasmacytoid-dendritic-cell-depleted mice did not. Interferons blunted the protective matricellular programs and contractile function of borderzone fibroblasts, and increased vulnerability to pathological remodelling. In mice that died after myocardial infarction, IFNIC colonies were immediately adjacent to sites of ventricular rupture, while mice lacking IFNICs were protected from rupture and exhibited improved survival3. Together, these results reveal a pathological borderzone niche characterized by a cardiomyocyte-initiated innate immune response. We suggest that selective inhibition of IRF3 activation in non-immune cells could limit ischaemic cardiomyopathy while avoiding broad immunosuppression.

Zinc (Zn) is a major soil contaminant and high Zn levels can disrupt growth, survival, and reproduction of fungi. Some fungal species evolved Zn tolerance through cell processes mitigating Zn toxicity, although the genes and detailed mechanisms underlying mycorrhizal fungal Zn tolerance remain unexplored. To fill this gap in knowledge, we investigated the gene expression of Zn tolerance in the ectomycorrhizal fungus Suillus luteus. We found that Zn tolerance in this species is mainly a constitutive trait that can also be environmentally dependent. Zinc tolerance in S. luteus is associated with differences in the expression of genes involved in metal exclusion and immobilization, as well as recognition and mitigation of metal-induced oxidative stress. Differentially expressed genes were predicted to be involved in transmembrane transport, metal chelation, oxidoreductase activity, and signal transduction. Some of these genes were previously reported as candidates for S. luteus Zn tolerance, while others are reported here for the first time. Our results contribute to understanding the mechanisms of fungal metal tolerance and pave the way for further research on the role of fungal metal tolerance in mycorrhizal associations.

Understanding cell fate decision-making during complex biological processes is an open challenge that is now aided by high-resolution single-cell sequencing technologies. Specifically, it remains challenging to identify and characterize transition states corresponding to "tipping points" whereby cells commit to new cell states. Here, we present a computational method that takes advantage of single-cell transcriptomics data to infer the stability and gene regulatory networks (GRNs) along cell lineages. Our method uses the unspliced and spliced counts from single-cell RNA sequencing data and cell ordering along lineage trajectories to train an RNA splicing multivariate model, from which cell-state stability along the lineage is inferred based on spectral analysis of the model's Jacobian matrix. Moreover, the model infers the RNA cross-species interactions resulting in GRNs and their variation along the cell lineage. When applied to epithelial-mesenchymal transition in ovarian and lung cancer-derived cell lines, our model predicts a saddle-node transition between the epithelial and mesenchymal states passing through an unstable, intermediate cell state. Furthermore, we show that the underlying GRN controlling epithelial-mesenchymal transition rearranges during the transition, resulting in denser and less modular networks in the intermediate state. Overall, our method represents a flexible tool to study cell lineages with a combination of theory-driven modeling and single-cell transcriptomics data.

Inferring the driving regulatory programs from comparative analysis of gene expression data is a cornerstone of systems biology. Many computational frameworks were developed to address this problem, including our iPAGE (information-theoretic Pathway Analysis of Gene Expression) toolset that uses information theory to detect non-random patterns of expression associated with given pathways or regulons. Our recent observations, however, indicate that existing approaches are susceptible to the technical biases that are inherent to most real world annotations. To address this, we have extended our information-theoretic framework to account for specific biases and artifacts in biological networks using the concept of conditional information. To showcase pyPAGE, we performed a comprehensive analysis of regulatory perturbations that underlie the molecular etiology of Alzheimers disease (AD). pyPAGE successfully recapitulated several known AD-associated gene expression programs. We also discovered several additional regulons whose differential activity is significantly associated with AD. We further explored how these regulators relate to pathological processes in AD through cell-type specific analysis of single cell and spatial gene expression datasets. Our findings showcase the utility of pyPAGE as a precise and reliable biomarker discovery in complex diseases such as Alzheimers disease.

Multiple myeloma is a treatable, but currently incurable, hematological malignancy of plasma cells characterized by diverse and complex tumor genetics for which precision medicine approaches to treatment are lacking. The Multiple Myeloma Research Foundations Relating Clinical Outcomes in Multiple Myeloma to Personal Assessment of Genetic Profile study ( NCT01454297 ) is a longitudinal, observational clinical study of newly diagnosed patients with multiple myeloma (n = 1,143) where tumor samples are characterized using whole-genome sequencing, whole-exome sequencing and RNA sequencing at diagnosis and progression, and clinical data are collected every 3 months. Analyses of the baseline cohort identified genes that are the target of recurrent gain-of-function and loss-of-function events. Consensus clustering identified 8 and 12 unique copy number and expression subtypes of myeloma, respectively, identifying high-risk genetic subtypes and elucidating many of the molecular underpinnings of these unique biological groups. Analysis of serial samples showed that 25.5% of patients transition to a high-risk expression subtype at progression. We observed robust expression of immunotherapy targets in this subtype, suggesting a potential therapeutic option.

BACKGROUND: RNA sequencing (RNA-Seq) offers profound insights into the complex transcriptomes of diverse biological systems. However, standard differential expression analysis pipelines based on DESeq2 and edgeR encounter challenges when applied to the immediate early transcriptomes of Chlamydia spp., obligate intracellular bacteria. These challenges arise from their reliance on assumptions that do not hold in scenarios characterized by extensive transcriptomic activation and limited repression. RESULTS: Standard analyses using unique chlamydial RNA-Seq reads alone identify nearly 300 upregulated and about 300 downregulated genes, significantly deviating from actual RNA-Seq read trends. By incorporating both chlamydial and host reads or adjusting for total sequencing depth, the revised normalization methods each detected over 700 upregulated genes and 30 or fewer downregulated genes, closely aligned with observed RNA-Seq data. Further validation through qRT-PCR analysis confirmed the effectiveness of these adjusted approaches in capturing the true extent of transcriptomic activation during the immediate early phase of chlamydial infection. CONCLUSIONS: This study highlights the limitations of standard RNA-Seq analysis tools in scenarios with extensive transcriptomic activation, such as in Chlamydia spp. during early infection. Our revised normalization methods, incorporating host reads or total sequencing depth, provide a more accurate representation of gene expression dynamics. These approaches may inform similar adjustments in other systems with unbalanced gene expression dynamics, enhancing the accuracy of transcriptomic analysis.

BACKGROUND: Two important questions regarding the genetics of pancreatic adenocarcinoma (PDAC) are 1. Which germline genetic variants influence the incidence of this cancer; and 2. Whether PDAC causally predisposes to associated non-malignant phenotypes, such as type 2 diabetes (T2D) and venous thromboembolism (VTE). METHODS: In this study of 8803 patients with PDAC and 67,523 controls, we first performed a large-scale transcriptome-wide association study to investigate the association between genetically determined gene expression in normal pancreas tissue and PDAC risk. Secondly, we used Mendelian Randomization (MR) to analyse the causal relationships among PDAC, T2D (74,124 cases and 824,006 controls) and VTE (30,234 cases and 172,122 controls). FINDINGS: Sixteen genes showed an association with PDAC risk (FDR

Merging diverse single-cell RNA sequencing (scRNA-seq) data from numerous experiments, laboratories and technologies can uncover important biological insights. Nonetheless, integrating scRNA-seq data encounters special challenges when the datasets are composed of diverse cell type compositions. Scanorama offers a robust solution for improving the quality and interpretation of heterogeneous scRNA-seq data by effectively merging information from diverse sources. Scanorama is designed to address the technical variation introduced by differences in sample preparation, sequencing depth and experimental batches that can confound the analysis of multiple scRNA-seq datasets. Here we provide a detailed protocol for using Scanorama within a Scanpy-based single-cell analysis workflow coupled with Google Colaboratory, a cloud-based free Jupyter notebook environment service. The protocol involves Scanorama integration, a process that typically spans 0.5-3 h. Scanorama integration requires a basic understanding of cellular biology, transcriptomic technologies and bioinformatics. Our protocol and new Scanorama-Colaboratory resource should make scRNA-seq integration more widely accessible to researchers.

Microgravity has been shown to lead to both muscle atrophy and impaired muscle regeneration. The purpose was to study the efficacy of microgravity to model impaired muscle regeneration in an engineered muscle platform and then to demonstrate the feasibility of performing drug screening in this model. Engineered human muscle was launched to the International Space Station National Laboratory, where the effect of microgravity exposure for 7 days was examined by transcriptomics and proteomics approaches. Gene set enrichment analysis of engineered muscle cultured in microgravity, compared to normal gravity conditions, highlighted a metabolic shift toward lipid and fatty acid metabolism, along with increased apoptotic gene expression. The addition of pro-regenerative drugs, insulin-like growth factor-1 (IGF-1) and a 15-hydroxyprostaglandin dehydrogenase inhibitor (15-PGDH-i), partially inhibited the effects of microgravity. In summary, microgravity mimics aspects of impaired myogenesis, and the addition of these drugs could partially inhibit the effects induced by microgravity.

BackgroundDeeper insights into ERBB2-driven cancers are essential to develop new treatment approaches for ERBB2+ breast cancers (BCs). We employed the Collaborative Cross (CC) mouse model to unearth genetic factors underpinning Erbb2-driven mammary tumour development and metastasis.Methods732 F1 hybrid female mice between FVB/N MMTV-Erbb2 and 30 CC strains were monitored for mammary tumour phenotypes. GWAS pinpointed SNPs that influence various tumour phenotypes. Multivariate analyses and models were used to construct the polygenic score and to develop a mouse tumour susceptibility gene signature (mTSGS), where the corresponding human ortholog was identified and designated as hTSGS. The importance and clinical value of hTSGS in human BC was evaluated using public datasets, encompassing TCGA, METABRIC, GSE96058, and I-SPY2 cohorts. The predictive power of mTSGS for response to chemotherapy was validated in vivo using genetically diverse MMTV-Erbb2 mice.FindingsDistinct variances in tumour onset, multiplicity, and metastatic patterns were observed in F1-hybrid female mice between FVB/N MMTV-Erbb2 and 30 CC strains. Besides lung metastasis, liver and kidney metastases emerged in specific CC strains. GWAS identified specific SNPs significantly associated with tumour onset, multiplicity, lung metastasis, and liver metastasis. Multivariate analyses flagged SNPs in 20 genes (Stx6, Ramp1, Traf3ip1, Nckap5, Pfkfb2, Trmt1l, Rprd1b, Rer1, Sepsecs, Rhobtb1, Tsen15, Abcc3, Arid5b, Tnr, Dock2, Tti1, Fam81a, Oxr1, Plxna2, and Tbc1d31) independently tied to various tumour characteristics, designated as a mTSGS. hTSGS scores (hTSGSS) based on their transcriptional level showed prognostic values, superseding clinical factors and PAM50 subtype across multiple human BC cohorts, and predicted pathological complete response independent of and superior to MammaPrint score in I-SPY2 study. The power of mTSGS score for predicting chemotherapy response was further validated in an in vivo mouse MMTV-Erbb2 model, showing that, like findings in human patients, mouse tumours with low mTSGS scores were most likely to respond to treatment.InterpretationOur investigation has unveiled many new genes predisposing individuals to ERBB2-driven cancer. Translational findings indicate that hTSGS holds promise as a biomarker for refining treatment strategies for patients with BC.FundingThe U.S. Department of Defense (DoD) Breast Cancer Research Program (BCRP) (BC190820), United States; MCIN/AEI/10.13039/501100011039 (PID2020-118527RB-I00, PDC2021-121735-I00), the "European Union Next Generation EU/PRTR," the Regional Government of Castile and León (CSI144P20), European Union.

Natural killer (NK) cells are innate lymphoid cells (ILCs) contributing to immune responses to microbes and tumors. Historically, their classification hinged on a limited array of surface protein markers. Here, we used single-cell RNA sequencing (scRNA-seq) and cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) to dissect the heterogeneity of NK cells. We identified three prominent NK cell subsets in healthy human blood: NK1, NK2 and NK3, further differentiated into six distinct subgroups. Our findings delineate the molecular characteristics, key transcription factors, biological functions, metabolic traits and cytokine responses of each subgroup. These data also suggest two separate ontogenetic origins for NK cells, leading to divergent transcriptional trajectories. Furthermore, we analyzed the distribution of NK cell subsets in the lung, tonsils and intraepithelial lymphocytes isolated from healthy individuals and in 22 tumor types. This standardized terminology aims at fostering clarity and consistency in future research, thereby improving cross-study comparisons.

ObjectivePeople with HIV (PWH) are at greater risk of developing lung diseases even when they are antiretroviral therapy (ART)-adherent and virally suppressed. The most common pulmonary function abnormality in PWH is that of impaired diffusing capacity of the lungs for carbon monoxide (DL CO ), which is an independent risk factor for increased mortality in PWH. Earlier work has identified several plasma biomarkers of inflammation and immune activation to be associated with decreased DL CO . However, the underpinning molecular mechanisms of HIV-associated impaired DL CO are largely unknown.DesignCross-sectional pilot study with PWH with normal DL CO (values greater than or equal to the lower limit of normal, DL CO  ≥ LLN, N = 9) or abnormal DL CO (DL CO

Systemic lupus erythematosus (SLE) is prototypical autoimmune disease driven by pathological T cell-B cell interactions1,2. Expansion of T follicular helper (TFH) and T peripheral helper (TPH) cells, two T cell populations that provide help to B cells, is a prominent feature of SLE3,4. Human TFH and TPH cells characteristically produce high levels of the B cell chemoattractant CXCL13 (refs. 5,6), yet regulation of T cell CXCL13 production and the relationship between CXCL13+ T cells and other T cell states remains unclear. Here, we identify an imbalance in CD4+ T cell phenotypes in patients with SLE, with expansion of PD-1+/ICOS+ CXCL13+ T cells and reduction of CD96hi IL-22+ T cells. Using CRISPR screens, we identify the aryl hydrocarbon receptor (AHR) as a potent negative regulator of CXCL13 production by human CD4+ T cells. Transcriptomic, epigenetic and functional studies demonstrate that AHR coordinates with AP-1 family member JUN to prevent CXCL13+ TPH/TFH cell differentiation and promote an IL-22+ phenotype. Type I interferon, a pathogenic driver of SLE7, opposes AHR and JUN to promote T cell production of CXCL13. These results place CXCL13+ TPH/TFH cells on a polarization axis opposite from T helper 22 (TH22) cells and reveal AHR, JUN and interferon as key regulators of these divergent T cell states.

Deciphering the genetic basis of prostate-specific antigen (PSA) levels may improve their utility for prostate cancer (PCa) screening. Using genome-wide association study (GWAS) summary statistics from 95,768 PCa-free men, we conducted a transcriptome-wide association study (TWAS) to examine impacts of genetically predicted gene expression on PSA. Analyses identified 41 statistically significant (p 0.5: CCNA2 and HIST1H2BN. Six of the 20 identified genes are not known to impact PCa risk. Fine-mapping based on whole blood and prostate tissue revealed five protein-coding genes with evidence of causal relationships with PSA levels. Of these five genes, four exhibited evidence of colocalization and one was conditionally independent of previous GWAS findings. These results yield hypotheses that should be further explored to improve understanding of genetic factors underlying PSA levels.

Melanoma clinical outcomes emerge from incompletely understood genetic mechanisms operating within the tumor and its microenvironment. Here, we used single-cell RNA-based spatial molecular imaging (RNA-SMI) in patient-derived archival tumors to reveal clinically relevant markers of malignancy progression and prognosis. We examined spatial gene expression of 203,472 cells inside benign and malignant melanocytic neoplasms, including melanocytic nevi and primary invasive and metastatic melanomas. Algorithmic cell clustering paired with intratumoral comparative two-dimensional analyses visualized synergistic, spatial gene signatures linking cellular proliferation, metabolism, and malignancy, validated by protein expression. Metastatic niches included up-regulation of CDK2 and FABP5, which independently predicted poor clinical outcome in 473 patients with melanoma via Cox regression analysis. More generally, our work demonstrates a framework for applying single-cell RNA-SMI technology toward identifying gene regulatory landscapes pertinent to cancer progression and patient survival.

Bone marrow aspirate concentrate (BMAC) and adipose-derived stromal vascular fraction (ADSVF) are the most marketed stem cell therapies to treat a variety of conditions in the general population and elite athletes. Both tissues have been used interchangeably clinically even though their detailed composition, heterogeneity, and mechanisms of action have neither been rigorously inventoried nor compared. This lack of information has prevented investigations into ideal dosages and has facilitated anecdata and misinformation. Here, we analyzed single-cell transcriptomes, proteomes, and flow cytometry profiles from paired clinical-grade BMAC and ADSVF. This comparative transcriptional atlas challenges the prevalent notion that there is one therapeutic cell type present in both tissues. We also provide data of surface markers that may enable isolation and investigation of cell (sub)populations. Furthermore, the proteome atlas highlights intertissue and interpatient heterogeneity of injected proteins with potentially regenerative or immunomodulatory capacities. An interactive webtool is available online.

The Long-read RNA-Seq Genome Annotation Assessment Project Consortium was formed to evaluate the effectiveness of long-read approaches for transcriptome analysis. Using different protocols and sequencing platforms, the consortium generated over 427 million long-read sequences from complementary DNA and direct RNA datasets, encompassing human, mouse and manatee species. Developers utilized these data to address challenges in transcript isoform detection, quantification and de novo transcript detection. The study revealed that libraries with longer, more accurate sequences produce more accurate transcripts than those with increased read depth, whereas greater read depth improved quantification accuracy. In well-annotated genomes, tools based on reference sequences demonstrated the best performance. Incorporating additional orthogonal data and replicate samples is advised when aiming to detect rare and novel transcripts or using reference-free approaches. This collaborative study offers a benchmark for current practices and provides direction for future method development in transcriptome analysis.

Cardiac allograft vasculopathy (CAV) causes late graft failure and mortality after heart transplantation. Donor-specific antibodies (DSAs) lead to chronic endothelial cell injury, inflammation, and arterial intimal thickening. In this study, GeoMx digital spatial profiling was used to analyze arterial areas of interest (AOIs) from CAV+DSA+ rejected cardiac allografts (N = 3; 22 AOIs total). AOIs were categorized based on CAV neointimal thickening and underwent whole transcriptome and protein profiling. By comparing our transcriptomic data with that of healthy control vessels of rapid autopsy myocardial tissue, we pinpointed specific pathways and transcripts indicative of heightened inflammatory profiles in CAV lesions. Moreover, we identified protein and transcriptomic signatures distinguishing CAV lesions exhibiting low and high neointimal lesions. AOIs with low neointima showed increased markers for activated inflammatory infiltrates, endothelial cell activation transcripts, and gene modules involved in metalloproteinase activation and TP53 regulation of caspases. Inflammatory and apoptotic proteins correlated with inflammatory modules in low neointima AOIs. High neointima AOIs exhibited elevated TGFβ-regulated transcripts and modules enriched for platelet activation/aggregation. Proteins associated with growth factors/survival correlated with modules enriched for proliferation/repair in high neointima AOIs. Our findings reveal novel insight into immunological mechanisms mediating CAV pathogenesis.