Pan-cancer single-cell dissection reveals phenotypically distinct B cell subtypes.

Characterizing the compositional and phenotypic characteristics of tumor-infiltrating B cells (TIBs) is important for advancing our understanding of their role in cancer development. Here, we establish a comprehensive resource of human B cells by integrating single-cell RNA sequencing data of B cells from 649 patients across 19 major cancer types. We demonstrate substantial heterogeneity in their total abundance and subtype composition and observe immunoglobulin G (IgG)-skewness of antibody-secreting cell isotypes. Moreover, we identify stress-response memory B cells and tumor-associated atypical B cells (TAABs), two tumor-enriched subpopulations with prognostic potential, shared in a pan-cancer manner. In particular, TAABs, characterized by a high clonal expansion level and proliferative capacity as well as by close interactions with activated CD4 T cells in tumors, are predictive of immunotherapy response. Our integrative resource depicts distinct clinically relevant TIB subsets, laying a foundation for further exploration of functional commonality and diversity of B cells in cancer.

A pan-cancer single-cell transcriptional atlas of tumor infiltrating myeloid cells.

Tumor-infiltrating myeloid cells (TIMs) are key regulators in tumor progression, but the similarity and distinction of their fundamental properties across different tumors remain elusive. Here, by performing a pan-cancer analysis of single myeloid cells from 210 patients across 15 human cancer types, we identified distinct features of TIMs across cancer types. Mast cells in nasopharyngeal cancer were found to be associated with better prognosis and exhibited an anti-tumor phenotype with a high ratio of TNF+/VEGFA+ cells. Systematic comparison between cDC1- and cDC2-derived LAMP3+ cDCs revealed their differences in transcription factors and external stimulus. Additionally, pro-angiogenic tumor-associated macrophages (TAMs) were characterized with diverse markers across different cancer types, and the composition of TIMs appeared to be associated with certain features of somatic mutations and gene expressions. Our results provide a systematic view of the highly heterogeneous TIMs and suggest future avenues for rational, targeted immunotherapies.

COVID-19 immune features revealed by a large-scale single-cell transcriptome atlas.

A dysfunctional immune response in coronavirus disease 2019 (COVID-19) patients is a recurrent theme impacting symptoms and mortality, yet a detailed understanding of pertinent immune cells is not complete. We applied single-cell RNA sequencing to 284 samples from 196 COVID-19 patients and controls and created a comprehensive immune landscape with 1.46 million cells. The large dataset enabled us to identify that different peripheral immune subtype changes are associated with distinct clinical features, including age, sex, severity, and disease stages of COVID-19. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) RNA was found in diverse epithelial and immune cell types, accompanied by dramatic transcriptomic changes within virus-positive cells. Systemic upregulation of S100A8/A9, mainly by megakaryocytes and monocytes in the peripheral blood, may contribute to the cytokine storms frequently observed in severe patients. Our data provide a rich resource for understanding the pathogenesis of and developing effective therapeutic strategies for COVID-19.

Redox State Controls Phase Separation of the Yeast Ataxin-2 Protein via Reversible Oxidation of Its Methionine-Rich Low-Complexity Domain.

Yeast ataxin-2, also known as Pbp1, senses the activity state of mitochondria in order to regulate TORC1. A domain of Pbp1 required to adapt cells to mitochondrial activity is of low sequence complexity. The low-complexity (LC) domain of Pbp1 forms labile, cross-ß polymers that facilitate phase transition of the protein into liquid-like or gel-like states. Phase transition for other LC domains is reliant upon widely distributed aromatic amino acids. In place of tyrosine or phenylalanine residues prototypically used for phase separation, Pbp1 contains 24 similarly disposed methionine residues. Here, we show that the Pbp1 methionine residues are sensitive to hydrogen peroxide (H2O2)-mediated oxidation in vitro and in living cells. Methionine oxidation melts Pbp1 liquid-like droplets in a manner reversed by methionine sulfoxide reductase enzymes. These observations explain how reversible formation of labile polymers by the Pbp1 LC domain enables the protein to function as a sensor of cellular redox state.

Yeast Ataxin-2 Forms an Intracellular Condensate Required for the Inhibition of TORC1 Signaling during Respiratory Growth.

Yeast ataxin-2, also known as Pbp1 (polyA binding protein-binding protein 1), is an intrinsically disordered protein implicated in stress granule formation, RNA biology, and neurodegenerative disease. To understand the endogenous function of this protein, we identify Pbp1 as a dedicated regulator of TORC1 signaling and autophagy under conditions that require mitochondrial respiration. Pbp1 binds to TORC1 specifically during respiratory growth, but utilizes an additional methionine-rich, low complexity (LC) region to inhibit TORC1. This LC region causes phase separation, forms reversible fibrils, and enables self-association into assemblies required for TORC1 inhibition. Mutants that weaken phase separation in vitro exhibit reduced capacity to inhibit TORC1 and induce autophagy. Loss of Pbp1 leads to mitochondrial dysfunction and reduced fitness during nutritional stress. Thus, Pbp1 forms a condensate in response to respiratory status to regulate TORC1 signaling.

Transcription factors TCF-1 and GATA3 are key factors for the epigenetic priming of early innate lymphoid progenitors toward distinct cell fates.

The differentiation of innate lymphoid cells (ILCs) from hematopoietic stem cells needs to go through several multipotent progenitor stages. However, it remains unclear whether the fates of multipotent progenitors are predefined by epigenetic states. Here, we report the identification of distinct accessible chromatin regions in all lymphoid progenitors (ALPs), EILPs, and ILC precursors (ILCPs). Single-cell MNase-seq analyses revealed that EILPs contained distinct subpopulations epigenetically primed toward either dendritic cell lineages or ILC lineages. We found that TCF-1 and GATA3 co-bound to the lineage-defining sites for ILCs (LDS-Is), whereas PU.1 binding was enriched in the LDSs for alternative dendritic cells (LDS-As). TCF-1 and GATA3 were indispensable for the epigenetic priming of LDSs at the EILP stage. Our results suggest that the multipotency of progenitor cells is defined by the existence of a heterogeneous population of cells epigenetically primed for distinct downstream lineages, which are regulated by key transcription factors.

Widespread macromolecular interaction perturbations in human genetic disorders.

How disease-associated mutations impair protein activities in the context of biological networks remains mostly undetermined. Although a few renowned alleles are well characterized, functional information is missing for over 100,000 disease-associated variants. Here we functionally profile several thousand missense mutations across a spectrum of Mendelian disorders using various interaction assays. The majority of disease-associated alleles exhibit wild-type chaperone binding profiles, suggesting they preserve protein folding or stability. While common variants from healthy individuals rarely affect interactions, two-thirds of disease-associated alleles perturb protein-protein interactions, with half corresponding to "edgetic" alleles affecting only a subset of interactions while leaving most other interactions unperturbed. With transcription factors, many alleles that leave protein-protein interactions intact affect DNA binding. Different mutations in the same gene leading to different interaction profiles often result in distinct disease phenotypes. Thus disease-associated alleles that perturb distinct protein activities rather than grossly affecting folding and stability are relatively widespread.

A lanthanide-rich kilonova in the aftermath of a long gamma-ray burst.

Observationally, kilonovae are astrophysical transients powered by the radioactive decay of nuclei heavier than iron, thought to be synthesized in the merger of two compact objects1-4. Over the first few days, the kilonova evolution is dominated by a large number of radioactive isotopes contributing to the heating rate2,5. On timescales of weeks to months, its behaviour is predicted to differ depending on the ejecta composition and the merger remnant6-8. Previous work has shown that the kilonova associated with gamma-ray burst 230307A is similar to kilonova AT2017gfo (ref. 9), and mid-infrared spectra revealed an emission line at 2.15 micrometres that was attributed to tellurium. Here we report a multi-wavelength analysis, including publicly available James Webb Space Telescope data9 and our own Hubble Space Telescope data, for the same gamma-ray burst. We model its evolution up to two months after the burst and show that, at these late times, the recession of the photospheric radius and the rapidly decaying bolometric luminosity (Lbol ∝ t-2.7±0.4, where t is time) support the recombination of lanthanide-rich ejecta as they cool.

Tumour-selective activity of RAS-GTP inhibition in pancreatic cancer.

Broad-spectrum RAS inhibition has the potential to benefit roughly a quarter of human patients with cancer whose tumours are driven by RAS mutations1,2. RMC-7977 is a highly selective inhibitor of the active GTP-bound forms of KRAS, HRAS and NRAS, with affinity for both mutant and wild-type variants3. More than 90% of cases of human pancreatic ductal adenocarcinoma (PDAC) are driven by activating mutations in KRAS4. Here we assessed the therapeutic potential of RMC-7977 in a comprehensive range of PDAC models. We observed broad and pronounced anti-tumour activity across models following direct RAS inhibition at exposures that were well-tolerated in vivo. Pharmacological analyses revealed divergent responses to RMC-7977 in tumour versus normal tissues. Treated tumours exhibited waves of apoptosis along with sustained proliferative arrest, whereas normal tissues underwent only transient decreases in proliferation, with no evidence of apoptosis. In the autochthonous KPC mouse model, RMC-7977 treatment resulted in a profound extension of survival followed by on-treatment relapse. Analysis of relapsed tumours identified Myc copy number gain as a prevalent candidate resistance mechanism, which could be overcome by combinatorial TEAD inhibition in vitro. Together, these data establish a strong preclinical rationale for the use of broad-spectrum RAS-GTP inhibition in the setting of PDAC and identify a promising candidate combination therapeutic regimen to overcome monotherapy resistance.

A long-duration gamma-ray burst with a peculiar origin.

It is generally believed that long-duration gamma-ray bursts (GRBs) are associated with massive star core collapse1, whereas short-duration GRBs are associated with mergers of compact star binaries2. However, growing observations3-6 have suggested that oddball GRBs do exist, and several criteria (prompt emission properties, supernova/kilonova associations and host galaxy properties) rather than burst duration only are needed to classify GRBs physically7. A previously reported long-duration burst, GRB 060614 (ref. 3), could be viewed as a short GRB with extended emission if it were observed at a larger distance8 and was associated with a kilonova-like feature9. As a result, it belongs to the type I (compact star merger) GRB category and is probably of binary neutron star (NS) merger origin. Here we report a peculiar long-duration burst, GRB 211211A, whose prompt emission properties in many aspects differ from all known type I GRBs, yet its multiband observations suggest a non-massive-star origin. In particular, substantial excess emission in both optical and near-infrared wavelengths has been discovered (see also ref. 10), which resembles kilonova emission, as observed in some type I GRBs. These observations point towards a new progenitor type of GRBs. A scenario invoking a white dwarf (WD)-NS merger with a post-merger magnetar engine provides a self-consistent interpretation for all the observations, including prompt gamma rays, early X-ray afterglow, as well as the engine-fed11,12 kilonova emission.

Structured community transitions explain the switching capacity of microbial systems.

Microbial systems appear to exhibit a relatively high switching capacity of moving back and forth among few dominant communities (taxon memberships). While this switching behavior has been mainly attributed to random environmental factors, it remains unclear the extent to which internal community dynamics affect the switching capacity of microbial systems. Here, we integrate ecological theory and empirical data to demonstrate that structured community transitions increase the dependency of future communities on the current taxon membership, enhancing the switching capacity of microbial systems. Following a structuralist approach, we propose that each community is feasible within a unique domain in environmental parameter space. Then, structured transitions between any two communities can happen with probability proportional to the size of their feasibility domains and inversely proportional to their distance in environmental parameter space-which can be treated as a special case of the gravity model. We detect two broad classes of systems with structured transitions: one class where switching capacity is high across a wide range of community sizes and another class where switching capacity is high only inside a narrow size range. We corroborate our theory using temporal data of gut and oral microbiota (belonging to class 1) as well as vaginal and ocean microbiota (belonging to class 2). These results reveal that the topology of feasibility domains in environmental parameter space is a relevant property to understand the changing behavior of microbial systems. This knowledge can be potentially used to understand the relevant community size at which internal dynamics can be operating in microbial systems.

An improved tetracycline-inducible expression system for fission yeast.

The ability to manipulate gene expression is valuable for elucidating gene function. In the fission yeast Schizosaccharomyces pombe, the most widely used regulatable expression system is the nmt1 promoter and its two attenuated variants. However, these promoters have limitations, including a long lag, incompatibility with rich media, and unsuitability for non-dividing cells. Here, we present a tetracycline-inducible system free of these shortcomings. Our system features the enotetS promoter, which achieves a similar induced level and a higher induction ratio compared to the nmt1 promoter, without exhibiting a lag. Additionally, our system includes four weakened enotetS variants, offering an expression range similar to the nmt1 series promoters but with more intermediate levels. To enhance usability, each promoter is combined with a Tet-repressor-expressing cassette in an integration plasmid. Importantly, our system can be used in non-dividing cells, enabling the development of a synchronous meiosis induction method with high spore viability. Moreover, our system allows for the shutdown of gene expression and the generation of conditional loss-of-function mutants. This system provides a versatile and powerful tool for manipulating gene expression in fission yeast.

Pbp1 associates with Puf3 and promotes translation of its target mRNAs involved in mitochondrial biogenesis.

Pbp1 (poly(A)-binding protein-binding protein 1) is a cytoplasmic stress granule marker that is capable of forming condensates that function in the negative regulation of TORC1 signaling under respiratory conditions. Polyglutamine expansions in its mammalian ortholog ataxin-2 lead to spinocerebellar dysfunction due to toxic protein aggregation. Here, we show that loss of Pbp1 in S. cerevisiae leads to decreased amounts of mRNAs and mitochondrial proteins which are targets of Puf3, a member of the PUF (Pumilio and FBF) family of RNA-binding proteins. We found that Pbp1 supports the translation of Puf3-target mRNAs in respiratory conditions, such as those involved in the assembly of cytochrome c oxidase and subunits of mitochondrial ribosomes. We further show that Pbp1 and Puf3 interact through their respective low complexity domains, which is required for Puf3-target mRNA translation. Our findings reveal a key role for Pbp1-containing assemblies in enabling the translation of mRNAs critical for mitochondrial biogenesis and respiration. They may further explain prior associations of Pbp1/ataxin-2 with RNA, stress granule biology, mitochondrial function, and neuronal health.

Broader functions of TIR domains in Arabidopsis immunity.

TIR domains are NAD-degrading enzymes that function during immune signaling in prokaryotes, plants, and animals. In plants, most TIR domains are incorporated into intracellular immune receptors termed TNLs. In Arabidopsis, TIR-derived small molecules bind and activate EDS1 heterodimers, which in turn activate RNLs, a class of cation channel-forming immune receptors. RNL activation drives cytoplasmic Ca2+ influx, transcriptional reprogramming, pathogen resistance, and host cell death. We screened for mutants that suppress an RNL activation mimic allele and identified a TNL, SADR1. Despite being required for the function of an autoactivated RNL, SADR1 is not required for defense signaling triggered by other tested TNLs. SADR1 is required for defense signaling initiated by some transmembrane pattern recognition receptors and contributes to the unbridled spread of cell death in lesion simulating disease 1. Together with RNLs, SADR1 regulates defense gene expression at infection site borders, likely in a non-cell autonomous manner. RNL mutants that cannot sustain this pattern of gene expression are unable to prevent disease spread beyond localized infection sites, suggesting that this pattern corresponds to a pathogen containment mechanism. SADR1 potentiates RNL-driven immune signaling not only through the activation of EDS1 but also partially independently of EDS1. We studied EDS1-independent TIR function using nicotinamide, an NADase inhibitor. Nicotinamide decreased defense induction from transmembrane pattern recognition receptors and decreased calcium influx, pathogen growth restriction, and host cell death following intracellular immune receptor activation. We demonstrate that TIR domains can potentiate calcium influx and defense and are thus broadly required for Arabidopsis immunity.

α-Amino-ß-Carboxymuconate-ε-Semialdehyde Decarboxylase Catalyzes Enol/Keto Tautomerization of Oxaloacetate.

ACMSD (α-amino-ß-carboxymuconate-ε-semialdehyde decarboxylase) is a key metalloenzyme critical for regulating de novo endogenous NAD+/NADH biosynthesis through the tryptophan-kynurenine pathway. This decarboxylase is a recognized target implicated in mitochondrial diseases and neurodegenerative disorders. However, unraveling its enzyme-substrate complex has been challenging due to its high catalytic efficiency. Here, we present a combined biochemical and structural study wherein we determined the crystal structure of ACMSD in complex with malonate. Our analysis revealed significant rearrangements in the active site, particularly in residues crucial for ACMS decarboxylation, including Arg51, Arg239* (a residue from an adjacent subunit), His228, and Trp194. Docking modeling studies proposed a putative ACMS binding mode. Additionally, we found that ACMSD catalyzes oxaloacetic acid (OAA) tautomerization at a rate of 6.51 ± 0.42 s-1 but not decarboxylation. The isomerase activity of ACMSD on OAA warrants further investigation in future biological studies. Subsequent mutagenesis studies and crystallographic analysis of W194A variant shed light on the roles of specific second-coordination sphere residues. Our findings indicate that Arg51 and Arg239* are crucial for OAA tautomerization. Moreover, our comparative analysis with related isomerase superfamily members underscores a general strategy employing arginine residues to promote OAA isomerization. Given the observed isomerase activity of ACMSD on OAA and its structural similarity to ACMS, we propose that ACMSD may facilitate isomerization to ensure ACMS is in the optimal tautomeric form for subsequent decarboxylation initiated by the zinc-bound hydroxide ion. Overall, these findings deepen the understanding of the structure and function of ACMSD, offering insights into potential therapeutic interventions.

A statistical physics approach for disease module detection.

Extensive evidence indicates that the pathobiological processes of a complex disease are associated with perturbation in specific neighborhoods of the human protein-protein interaction (PPI) network (also known as the interactome), often referred to as the disease module. Many computational methods have been developed to integrate the interactome and omics profiles to extract context-dependent disease modules. Yet, existing methods all have fundamental limitations in terms of rigor and/or efficiency. Here, we developed a statistical physics approach based on the random-field Ising model (RFIM) for disease module detection, which is both mathematically rigorous and computationally efficient. We applied our RFIM approach to genome-wide association studies (GWAS) of ten complex diseases to examine its performance for disease module detection. We found that our RFIM approach outperforms existing methods in terms of computational efficiency, connectivity of disease modules, and robustness to the interactome incompleteness.

MPCLCDA: predicting circRNA-disease associations by using automatically selected meta-path and contrastive learning.

Circular RNA (circRNA) is closely associated with human diseases. Accordingly, identifying the associations between human diseases and circRNA can help in disease prevention, diagnosis and treatment. Traditional methods are time consuming and laborious. Meanwhile, computational models can effectively predict potential circRNA-disease associations (CDAs), but are restricted by limited data, resulting in data with high dimension and imbalance. In this study, we propose a model based on automatically selected meta-path and contrastive learning, called the MPCLCDA model. First, the model constructs a new heterogeneous network based on circRNA similarity, disease similarity and known association, via automatically selected meta-path and obtains the low-dimensional fusion features of nodes via graph convolutional networks. Then, contrastive learning is used to optimize the fusion features further, and obtain the node features that make the distinction between positive and negative samples more evident. Finally, circRNA-disease scores are predicted through a multilayer perceptron. The proposed method is compared with advanced methods on four datasets. The average area under the receiver operating characteristic curve, area under the precision-recall curve and F1 score under 5-fold cross-validation reached 0.9752, 0.9831 and 0.9745, respectively. Simultaneously, case studies on human diseases further prove the predictive ability and application value of this method.

NSRGRN: a network structure refinement method for gene regulatory network inference.

The elucidation of gene regulatory networks (GRNs) is one of the central challenges of systems biology, which is crucial for understanding pathogenesis and curing diseases. Various computational methods have been developed for GRN inference, but identifying redundant regulation remains a fundamental problem. Although considering topological properties and edge importance measures simultaneously can identify and reduce redundant regulations, how to address their respective weaknesses whilst leveraging their strengths is a critical problem faced by researchers. Here, we propose a network structure refinement method for GRN (NSRGRN) that effectively combines the topological properties and edge importance measures during GRN inference. NSRGRN has two major parts. The first part constructs a preliminary ranking list of gene regulations to avoid starting the GRN inference from a directed complete graph. The second part develops a novel network structure refinement (NSR) algorithm to refine the network structure from local and global topology perspectives. Specifically, the Conditional Mutual Information with Directionality and network motifs are applied to optimise the local topology, and the lower and upper networks are used to balance the bilateral relationship between the local topology's optimisation and the global topology's maintenance. NSRGRN is compared with six state-of-the-art methods on three datasets (26 networks in total), and it shows the best all-round performance. Furthermore, when acting as a post-processing step, the NSR algorithm can improve the results of other methods in most datasets.

The role of ncRNA regulatory mechanisms in diseases-case on gestational diabetes.

Non-coding RNAs (ncRNAs) are a class of RNA molecules that do not have the potential to encode proteins. Meanwhile, they can occupy a significant portion of the human genome and participate in gene expression regulation through various mechanisms. Gestational diabetes mellitus (GDM) is a pathologic condition of carbohydrate intolerance that begins or is first detected during pregnancy, making it one of the most common pregnancy complications. Although the exact pathogenesis of GDM remains unclear, several recent studies have shown that ncRNAs play a crucial regulatory role in GDM. Herein, we present a comprehensive review on the multiple mechanisms of ncRNAs in GDM along with their potential role as biomarkers. In addition, we investigate the contribution of deep learning-based models in discovering disease-specific ncRNA biomarkers and elucidate the underlying mechanisms of ncRNA. This might assist community-wide efforts to obtain insights into the regulatory mechanisms of ncRNAs in disease and guide a novel approach for early diagnosis and treatment of disease.