Small nucleolar RNAs (snoRNAs) are non-coding RNAs known for guiding RNA modifications, including 2'-O-methylation (Nm) and pseudouridine (Ψ). While snoRNAs may also interact with other RNAs, such as mRNA, the full repertoire of RNAs targeted by snoRNA remains elusive due to the lack of effective technologies that identify snoRNA targets transcriptome wide. Here, we develop a chemical crosslinking-based approach that comprehensively detects cellular RNA targets of snoRNAs, yielding thousands of previously unrecognized snoRNA-mRNA interactions in human cells and mouse brain tissues. Many interactions occur outside of snoRNA-guided RNA modification sites, hinting at non-canonical functions beyond RNA modification. We find that one of these snoRNAs, SNORA73, targets mRNAs that encode secretory proteins and membrane proteins. SNORA73 also interacts with 7SL RNA, part of the signal recognition particle (SRP) required for protein secretion. The mRNA-SNORA73-7SL RNA interactions enhance the association of the SNORA73-target mRNAs with SRP, thereby facilitating the secretion of encoded proteins.

Brown seaweeds are keystone species of coastal ecosystems, often forming extensive underwater forests, and are under considerable threat from climate change. In this study, analysis of multiple genomes has provided insights across the entire evolutionary history of this lineage, from initial emergence, through later diversification of the brown algal orders, down to microevolutionary events at the genus level. Emergence of the brown algal lineage was associated with a marked gain of new orthologous gene families, enhanced protein domain rearrangement, increased horizontal gene transfer events, and the acquisition of novel signaling molecules and key metabolic pathways, the latter notably related to biosynthesis of the alginate-based extracellular matrix, and halogen and phlorotannin biosynthesis. We show that brown algal genome diversification is tightly linked to phenotypic divergence, including changes in life cycle strategy and zoid flagellar structure. The study also showed that integration of large viral genomes has had a significant impact on brown algal genome content throughout the emergence of the lineage.

Stringent control of T cell activity in the tumor microenvironment is essential for the generation of protective antitumor immunity. However, the identity, differentiation, and functions of the cells that create critical fibroblastic niches promoting tumor-infiltrating T cells remain elusive. Here, we show that CCL19-expressing fibroblastic reticular cells (FRCs) generate interconnected T cell environments (TEs) in human non-small cell lung cancer, including tertiary lymphoid structures and T cell tracks. Analysis of the FRC-T cell interactome in TEs indicated molecular networks regulating niche-specific differentiation of CCL19-expressing fibroblasts and T cell activation pathways. Single-cell transcriptomics and cell fate-mapping analyses in mice confirmed that FRCs in TEs originate from mural and adventitial progenitors. Ablation of intratumoral FRC precursors decreased antitumor T cell activity, resulting in reduced tumor control during coronavirus vector-based immunotherapy. In summary, specialized FRC niches in the tumor microenvironment govern the quality and extent of antitumor T cell immunity.

We present replication-aware single-molecule accessibility mapping (RASAM), a method to nondestructively measure replication status and protein-DNA interactions on chromatin genome-wide. Using RASAM, we uncover a genome-wide state of single-molecule "hyperaccessibility" post-replication that resolves over several hours. Combining RASAM with cellular models for rapid protein degradation, we demonstrate that histone chaperone CAF-1 reduces nascent chromatin accessibility by filling single-molecular "gaps" and generating closely spaced dinucleosomes on replicated DNA. At cis-regulatory elements, we observe unique modes by which nascent chromatin hyperaccessibility resolves: at CCCTC-binding factor (CTCF)-binding sites, CTCF and nucleosomes compete, reducing CTCF occupancy and motif accessibility post-replication; at active transcription start sites, high chromatin accessibility is maintained, implying rapid re-establishment of nucleosome-free regions. Our study introduces a new paradigm for studying replicated chromatin fiber organization. More broadly, we uncover a unique organization of newly replicated chromatin that must be reset by active processes, providing a substrate for epigenetic reprogramming.

Underlying variation in height are regulatory changes to chondrocytes, cartilage cells comprising long-bone growth plates. Currently, we lack knowledge on epigenetic regulation and gene expression of chondrocytes sampled across the human skeleton, and therefore we cannot understand basic regulatory mechanisms controlling height biology. We first rectify this issue by generating extensive epigenetic and transcriptomic maps from chondrocytes sampled from different growth plates across developing human skeletons, discovering novel regulatory networks shaping human bone/joint development. Next, using these maps in tandem with height genome-wide association study (GWAS) signals, we disentangle the regulatory impacts that skeletal element-specific versus global-acting variants have on skeletal growth, revealing the prime importance of regulatory pleiotropy in controlling height variation. Finally, as height is highly heritable, and thus often the test case for complex-trait genetics, we leverage these datasets within a testable omnigenic model framework to discover novel chondrocyte developmental modules and peripheral-acting factors shaping height biology and skeletal growth.

Transcription-coupled DNA repair (TCR) removes bulky DNA lesions impeding RNA polymerase II (RNAPII) transcription. Recent studies have outlined the stepwise assembly of TCR factors CSB, CSA, UVSSA, and transcription factor IIH (TFIIH) around lesion-stalled RNAPII. However, the mechanism and factors required for the transition to downstream repair steps, including RNAPII removal to provide repair proteins access to the DNA lesion, remain unclear. Here, we identify STK19 as a TCR factor facilitating this transition. Loss of STK19 does not impact initial TCR complex assembly or RNAPII ubiquitylation but delays lesion-stalled RNAPII clearance, thereby interfering with the downstream repair reaction. Cryoelectron microscopy (cryo-EM) and mutational analysis reveal that STK19 associates with the TCR complex, positioning itself between RNAPII, UVSSA, and CSA. The structural insights and molecular modeling suggest that STK19 positions the ATPase subunits of TFIIH onto DNA in front of RNAPII. Together, these findings provide new insights into the factors and mechanisms required for TCR.

In transcription-coupled nucleotide excision repair (TC-NER), stalled RNA polymerase II (RNA Pol II) binds CSB and CRL4CSA, which cooperate with UVSSA and ELOF1 to recruit TFIIH. To explore the mechanism of TC-NER, we recapitulated this reaction in vitro. When a plasmid containing a site-specific lesion is transcribed in frog egg extract, error-free repair is observed that depends on CSB, CRL4CSA, UVSSA, and ELOF1. Repair also requires STK19, a factor previously implicated in transcription recovery after UV exposure. A 1.9-Å cryo-electron microscopy structure shows that STK19 binds the TC-NER complex through CSA and the RPB1 subunit of RNA Pol II. Furthermore, AlphaFold predicts that STK19 interacts with the XPD subunit of TFIIH, and disrupting this interface impairs cell-free repair. Molecular modeling suggests that STK19 positions TFIIH ahead of RNA Pol II for lesion verification. Our analysis of cell-free TC-NER suggests that STK19 couples RNA Pol II stalling to downstream repair events.

The identification of individual protein-protein interactions (PPIs) began more than 40 years ago, using protein affinity chromatography and antibody co-immunoprecipitation. As new technologies emerged, analysis of PPIs increased to a genome-wide scale with the introduction of intracellular tagging methods, affinity purification (AP) followed by mass spectrometry (MS), and co-fractionation MS (CF-MS). Now, combining the resulting catalogs of interactions with complementary methods, including crosslinking MS (XL-MS) and cryogenic electron microscopy (cryo-EM), helps distinguish direct interactions from indirect ones within the same or between different protein complexes. These powerful approaches and the promise of artificial intelligence applications like AlphaFold herald a future where PPIs and protein complexes, including energy-driven protein machines, will be understood in exquisite detail, unlocking new insights in the contexts of both basic biology and disease.

Every cell must solve the problem of how to fold its genome. We describe how the folded state of chromosomes is the result of the combined activity of multiple conserved mechanisms. Homotypic affinity-driven interactions lead to spatial partitioning of active and inactive loci. Molecular motors fold chromosomes through loop extrusion. Topological features such as supercoiling and entanglements contribute to chromosome folding and its dynamics, and tethering loci to sub-nuclear structures adds additional constraints. Dramatically diverse chromosome conformations observed throughout the cell cycle and across the tree of life can be explained through differential regulation and implementation of these basic mechanisms. We propose that the first functions of chromosome folding are to mediate genome replication, compaction, and segregation and that mechanisms of folding have subsequently been co-opted for other roles, including long-range gene regulation, in different conditions, cell types, and species.

The microbiota in individual habitats differ in both relative composition and absolute abundance. While sequencing approaches determine the relative abundances of taxa and genes, they do not provide information on their absolute abundances. Here, we developed a machine-learning approach to predict fecal microbial loads (microbial cells per gram) solely from relative abundance data. Applying our prediction model to a large-scale metagenomic dataset (n = 34,539), we demonstrated that microbial load is the major determinant of gut microbiome variation and is associated with numerous host factors, including age, diet, and medication. We further found that for several diseases, changes in microbial load, rather than the disease condition itself, more strongly explained alterations in patients' gut microbiome. Adjusting for this effect substantially reduced the statistical significance of the majority of disease-associated species. Our analysis reveals that the fecal microbial load is a major confounder in microbiome studies, highlighting its importance for understanding microbiome variation in health and disease.

Kernel dehydration rate (KDR) is a crucial production trait that affects mechanized harvesting and kernel quality in maize; however, the underlying mechanisms remain unclear. Here, we identified a quantitative trait locus (QTL), qKDR1, as a non-coding sequence that regulates the expression of qKDR1 REGULATED PEPTIDE GENE (RPG). RPG encodes a 31 amino acid micropeptide, microRPG1, which controls KDR by precisely modulating the expression of two genes, ZmETHYLENE-INSENSITIVE3-like 1 and 3, in the ethylene signaling pathway in the kernels after filling. microRPG1 is a Zea genus-specific micropeptide and originated de novo from a non-coding sequence. Knockouts of microRPG1 result in faster KDR in maize. By contrast, overexpression or exogenous application of the micropeptide shows the opposite effect both in maize and Arabidopsis. Our findings reveal the molecular mechanism of microRPG1 in kernel dehydration and provide an important tool for future crop breeding.

β-Hydroxybutyrate (BHB) is an abundant ketone body. To date, all known pathways of BHB metabolism involve the interconversion of BHB and primary energy intermediates. Here, we identify a previously undescribed BHB secondary metabolic pathway via CNDP2-dependent enzymatic conjugation of BHB and free amino acids. This BHB shunt pathway generates a family of anti-obesity ketone metabolites, the BHB-amino acids. Genetic ablation of CNDP2 in mice eliminates tissue amino acid BHB-ylation activity and reduces BHB-amino acid levels. The most abundant BHB-amino acid, BHB-Phe, is a ketosis-inducible congener of Lac-Phe that activates hypothalamic and brainstem neurons and suppresses feeding. Conversely, CNDP2-KO mice exhibit increased food intake and body weight following exogenous ketone ester supplementation or a ketogenic diet. CNDP2-dependent amino acid BHB-ylation and BHB-amino acid metabolites are also conserved in humans. Therefore, enzymatic amino acid BHB-ylation defines a ketone shunt pathway and bioactive ketone metabolites linked to energy balance.

Dorsal root ganglia (DRGs) play a crucial role in processing sensory information, making it essential to understand their development. Here, we construct a single-cell spatiotemporal transcriptomic atlas of human embryonic DRG. This atlas reveals the diversity of cell types and highlights the extrinsic signaling cascades and intrinsic regulatory hierarchies that guide cell fate decisions, including neuronal/glial lineage restriction, sensory neuron differentiation and specification, and the formation of neuron-satellite glial cell (SGC) units. Additionally, we identify a human-enriched NTRK3+/DCC+ nociceptor subtype, which is involved in multimodal nociceptive processing. Mimicking the programmed activation of signaling pathways in vivo, we successfully establish functional human DRG organoids and underscore the critical roles of transcriptional regulators in the fate commitment of unspecialized sensory neurons (uSNs). Overall, our research elucidates the multilevel signaling pathways and transcription factor (TF) regulatory hierarchies that underpin the diversity of somatosensory neurons, emphasizing the phenotypic distinctions in human nociceptor subtypes.

Quantifying spatiotemporal dynamics during embryogenesis is crucial for understanding congenital diseases. We developed Spateo (https://github.com/aristoteleo/spateo-release), a 3D spatiotemporal modeling framework, and applied it to a 3D mouse embryogenesis atlas at E9.5 and E11.5, capturing eight million cells. Spateo enables scalable, partial, non-rigid alignment, multi-slice refinement, and mesh correction to create molecular holograms of whole embryos. It introduces digitization methods to uncover multi-level biology from subcellular to whole organ, identifying expression gradients along orthogonal axes of emergent 3D structures, e.g., secondary organizers such as midbrain-hindbrain boundary (MHB). Spateo further jointly models intercellular and intracellular interaction to dissect signaling landscapes in 3D structures, including the zona limitans intrathalamica (ZLI). Lastly, Spateo introduces "morphometric vector fields" of cell migration and integrates spatial differential geometry to unveil molecular programs underlying asymmetrical murine heart organogenesis and others, bridging macroscopic changes with molecular dynamics. Thus, Spateo enables the study of organ ecology at a molecular level in 3D space over time.

While apolipoprotein E (APOE) is the strongest genetic modifier for late-onset Alzheimer's disease (LOAD), the molecular mechanisms underlying isoform-dependent risk and the relevance of ApoE-associated lipids remain elusive. Here, we report that impaired low-density lipoprotein (LDL) receptor (LDLR) binding of lipidated ApoE2 (lipApoE2) avoids LDLR recycling defects observed with lipApoE3/E4 and decreases the uptake of cholesteryl esters (CEs), which are lipids linked to neurodegeneration. In human neurons, the addition of ApoE carrying polyunsaturated fatty acids (PUFAs)-CE revealed an allelic series (ApoE4 > ApoE3 > ApoE2) associated with lipofuscinosis, an age-related lysosomal pathology resulting from lipid peroxidation. Lipofuscin increased lysosomal accumulation of tau fibrils and was elevated in the APOE4 mouse brain with exacerbation by tau pathology. Intrahippocampal injection of PUFA-CE-lipApoE4 was sufficient to induce lipofuscinosis in wild-type mice. Finally, the protective Christchurch mutation also reduced LDLR binding and phenocopied ApoE2. Collectively, our data strongly suggest decreased lipApoE-LDLR interactions minimize LOAD risk by reducing the deleterious effects of endolysosomal targeting of ApoE and associated pathogenic lipids.

Proximity labeling technologies are limited to indexing localized protein residents. Such data-although valuable-cannot inform on small-molecule responsivity of local residents. We here bridge this gap by demonstrating in live C. elegans how electrophile-sensing propensity in specific organs can be quantitatively mapped and ranked. Using this method, >70% of tissue-specific responders exhibit electrophile responsivity, independent of tissue-specific abundance. One responder, cyp-33e1-for which both human and worm orthologs are electrophile responsive-marshals stress-dependent gut functions, despite manifesting uniform abundance across all tissues studied. Cyp-33e1's localized electrophile responsivity operates site specifically, triggering multifaceted responses: electrophile sensing through the catalytic-site cysteine results in partitioning between enzyme inhibition and localized production of a critical metabolite that governs global lipid availability, whereas rapid dual-cysteine site-specific sensing modulates gut homeostasis. Beyond pinpointing chemical actionability within local proteomes, organ-specific electrophile responsivity mapping illuminates otherwise intractable locale-specific metabolite signaling and stress response programs influencing organ-specific decision-making.

Ubiquitin-dependent proteolysis regulates diverse cellular functions with high substrate specificity, which hinges on the ability of ubiquitin E3 ligases to decode the targets' degradation signals, i.e., degrons. Here, we show that BACH1, a transcription repressor of antioxidant response genes, features two distinct unconventional degrons encrypted in the quaternary structure of its homodimeric BTB domain. These two degrons are both functionalized by oxidative stress and are deciphered by two complementary E3s. FBXO22 recognizes a degron constructed by the BACH1 BTB domain dimer interface, which is unmasked from transcriptional co-repressors after oxidative stress releases BACH1 from chromatin. When this degron is impaired by oxidation, a second BACH1 degron manifested by its destabilized BTB dimer is probed by a pair of FBXL17 proteins that remodels the substrate into E3-bound monomers for ubiquitination. Our findings highlight the multidimensionality of protein degradation signals and the functional complementarity of different ubiquitin ligases targeting the same substrate.

We envision "AI scientists" as systems capable of skeptical learning and reasoning that empower biomedical research through collaborative agents that integrate AI models and biomedical tools with experimental platforms. Rather than taking humans out of the discovery process, biomedical AI agents combine human creativity and expertise with AI's ability to analyze large datasets, navigate hypothesis spaces, and execute repetitive tasks. AI agents are poised to be proficient in various tasks, planning discovery workflows and performing self-assessment to identify and mitigate gaps in their knowledge. These agents use large language models and generative models to feature structured memory for continual learning and use machine learning tools to incorporate scientific knowledge, biological principles, and theories. AI agents can impact areas ranging from virtual cell simulation, programmable control of phenotypes, and the design of cellular circuits to developing new therapies.

The autophagy-lysosome system directs the degradation of a wide variety of cargo and is also involved in tumor progression. Here, we show that the immunity-related GTPase family Q protein (IRGQ), an uncharacterized protein to date, acts in the quality control of major histocompatibility complex class I (MHC class I) molecules. IRGQ directs misfolded MHC class I toward lysosomal degradation through its binding mode to GABARAPL2 and LC3B. In the absence of IRGQ, free MHC class I heavy chains do not only accumulate in the cell but are also transported to the cell surface, thereby promoting an immune response. Mice and human patients suffering from hepatocellular carcinoma show improved survival rates with reduced IRGQ levels due to increased reactivity of CD8+ T cells toward IRGQ knockout tumor cells. Thus, we reveal IRGQ as a regulator of MHC class I quality control, mediating tumor immune evasion.

Genetic variation in host immunity impacts the disproportionate burden of infectious diseases that can be experienced by First Nations peoples. Polymorphic human leukocyte antigen (HLA) class I and killer cell immunoglobulin-like receptors (KIRs) are key regulators of natural killer (NK) cells, which mediate early infection control. How this variation impacts their responses across populations is unclear. We show that HLA-A∗24:02 became the dominant ligand for inhibitory KIR3DL1 in First Nations peoples across Oceania, through positive natural selection. We identify KIR3DL1∗114, widespread across and unique to Oceania, as an allele lineage derived from archaic humans. KIR3DL1∗114+NK cells from First Nations Australian donors are inhibited through binding HLA-A∗24:02. The KIR3DL1∗114 lineage is defined by phenylalanine at residue 166. Structural and binding studies show phenylalanine 166 forms multiple unique contacts with HLA-peptide complexes, increasing both affinity and specificity. Accordingly, assessing immunogenetic variation and the functional implications for immunity are fundamental toward understanding population-based disease associations.