Thermogenic adipocytes possess a therapeutically appealing, energy-expending capacity, which is canonically cold-induced by ligand-dependent activation of β-adrenergic G protein-coupled receptors (GPCRs). Here, we uncover an alternate paradigm of GPCR-mediated adipose thermogenesis through the constitutively active receptor, GPR3. We show that the N terminus of GPR3 confers intrinsic signaling activity, resulting in continuous Gs-coupling and cAMP production without an exogenous ligand. Thus, transcriptional induction of Gpr3 represents the regulatory parallel to ligand-binding of conventional GPCRs. Consequently, increasing Gpr3 expression in thermogenic adipocytes is alone sufficient to drive energy expenditure and counteract metabolic disease in mice. Gpr3 transcription is cold-stimulated by a lipolytic signal, and dietary fat potentiates GPR3-dependent thermogenesis to amplify the response to caloric excess. Moreover, we find GPR3 to be an essential, adrenergic-independent regulator of human brown adipocytes. Taken together, our findings reveal a noncanonical mechanism of GPCR control and thermogenic activation through the lipolysis-induced expression of constitutively active GPR3.

Northern East Asia was inhabited by modern humans as early as 40 thousand years ago (ka), as demonstrated by the Tianyuan individual. Using genome-wide data obtained from 25 individuals dated to 33.6-3.4 ka from the Amur region, we show that Tianyuan-related ancestry was widespread in northern East Asia before the Last Glacial Maximum (LGM). At the close of the LGM stadial, the earliest northern East Asian appeared in the Amur region, and this population is basal to ancient northern East Asians. Human populations in the Amur region have maintained genetic continuity from 14 ka, and these early inhabitants represent the closest East Asian source known for Ancient Paleo-Siberians. We also observed that EDAR V370A was likely to have been elevated to high frequency after the LGM, suggesting the possible timing for its selection. This study provides a deep look into the population dynamics of northern East Asia.

Searching for factors to improve knockin efficiency for therapeutic applications, biotechnology, and generation of non-human primate models of disease, we found that the strand exchange protein RAD51 can significantly increase Cas9-mediated homozygous knockin in mouse embryos through an interhomolog repair (IHR) mechanism. IHR is a hallmark of meiosis but only occurs at low frequencies in somatic cells, and its occurrence in zygotes is controversial. Using multiple approaches, we provide evidence for an endogenous IHR mechanism in the early embryo that can be enhanced by RAD51. This process can be harnessed to generate homozygotes from wild-type zygotes using exogenous donors and to convert heterozygous alleles into homozygous alleles without exogenous templates. Furthermore, we identify additional IHR-promoting factors and describe features of IHR events. Together, our findings show conclusive evidence for IHR in mouse embryos and describe an efficient method for enhanced gene conversion.

We present a global atlas of 4,728 metagenomic samples from mass-transit systems in 60 cities over 3 years, representing the first systematic, worldwide catalog of the urban microbial ecosystem. This atlas provides an annotated, geospatial profile of microbial strains, functional characteristics, antimicrobial resistance (AMR) markers, and genetic elements, including 10,928 viruses, 1,302 bacteria, 2 archaea, and 838,532 CRISPR arrays not found in reference databases. We identified 4,246 known species of urban microorganisms and a consistent set of 31 species found in 97% of samples that were distinct from human commensal organisms. Profiles of AMR genes varied widely in type and density across cities. Cities showed distinct microbial taxonomic signatures that were driven by climate and geographic differences. These results constitute a high-resolution global metagenomic atlas that enables discovery of organisms and genes, highlights potential public health and forensic applications, and provides a culture-independent view of AMR burden in cities.

Long-term subcellular intravital imaging in mammals is vital to study diverse intercellular behaviors and organelle functions during native physiological processes. However, optical heterogeneity, tissue opacity, and phototoxicity pose great challenges. Here, we propose a computational imaging framework, termed digital adaptive optics scanning light-field mutual iterative tomography (DAOSLIMIT), featuring high-speed, high-resolution 3D imaging, tiled wavefront correction, and low phototoxicity with a compact system. By tomographic imaging of the entire volume simultaneously, we obtained volumetric imaging across 225 × 225 × 16 μm3, with a resolution of up to 220 nm laterally and 400 nm axially, at the millisecond scale, over hundreds of thousands of time points. To establish the capabilities, we investigated large-scale cell migration and neural activities in different species and observed various subcellular dynamics in mammals during neutrophil migration and tumor cell circulation.

Gasdermin B (GSDMB) belongs to a large family of pore-forming cytolysins that execute inflammatory cell death programs. While genetic studies have linked GSDMB polymorphisms to human disease, its function in the immunological response to pathogens remains poorly understood. Here, we report a dynamic host-pathogen conflict between GSDMB and the IpaH7.8 effector protein secreted by enteroinvasive Shigella flexneri. We show that IpaH7.8 ubiquitinates and targets GSDMB for 26S proteasome destruction. This virulence strategy protects Shigella from the bacteriocidic activity of natural killer cells by suppressing granzyme-A-mediated activation of GSDMB. In contrast to the canonical function of most gasdermin family members, GSDMB does not inhibit Shigella by lysing host cells. Rather, it exhibits direct microbiocidal activity through recognition of phospholipids found on Gram-negative bacterial membranes. These findings place GSDMB as a central executioner of intracellular bacterial killing and reveal a mechanism employed by pathogens to counteract this host defense system.

Organs are composed of diverse cell types that traverse transient states during organogenesis. To interrogate this diversity during human development, we generate a single-cell transcriptome atlas from multiple developing endodermal organs of the respiratory and gastrointestinal tract. We illuminate cell states, transcription factors, and organ-specific epithelial stem cell and mesenchyme interactions across lineages. We implement the atlas as a high-dimensional search space to benchmark human pluripotent stem cell (hPSC)-derived intestinal organoids (HIOs) under multiple culture conditions. We show that HIOs recapitulate reference cell states and use HIOs to reconstruct the molecular dynamics of intestinal epithelium and mesenchyme emergence. We show that the mesenchyme-derived niche cue NRG1 enhances intestinal stem cell maturation in vitro and that the homeobox transcription factor CDX2 is required for regionalization of intestinal epithelium and mesenchyme in humans. This work combines cell atlases and organoid technologies to understand how human organ development is orchestrated.

Natural antibodies (Abs) can target host glycans on the surface of pathogens. We studied the evolution of glycan-reactive B cells of rhesus macaques and humans using glycosylated HIV-1 envelope (Env) as a model antigen. 2G12 is a broadly neutralizing Ab (bnAb) that targets a conserved glycan patch on Env of geographically diverse HIV-1 strains using a unique heavy-chain (VH) domain-swapped architecture that results in fragment antigen-binding (Fab) dimerization. Here, we describe HIV-1 Env Fab-dimerized glycan (FDG)-reactive bnAbs without VH-swapped domains from simian-human immunodeficiency virus (SHIV)-infected macaques. FDG Abs also recognized cell-surface glycans on diverse pathogens, including yeast and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike. FDG precursors were expanded by glycan-bearing immunogens in macaques and were abundant in HIV-1-naive humans. Moreover, FDG precursors were predominately mutated IgM+IgD+CD27+, thus suggesting that they originated from a pool of antigen-experienced IgM+ or marginal zone B cells.

Organoids capable of forming tissue-like structures have transformed our ability to model human development and disease. With the notable exception of the human heart, lineage-specific self-organizing organoids have been reported for all major organs. Here, we established self-organizing cardioids from human pluripotent stem cells that intrinsically specify, pattern, and morph into chamber-like structures containing a cavity. Cardioid complexity can be controlled by signaling that instructs the separation of cardiomyocyte and endothelial layers and by directing epicardial spreading, inward migration, and differentiation. We find that cavity morphogenesis is governed by a mesodermal WNT-BMP signaling axis and requires its target HAND1, a transcription factor linked to developmental heart chamber defects. Upon cryoinjury, cardioids initiated a cell-type-dependent accumulation of extracellular matrix, an early hallmark of both regeneration and heart disease. Thus, human cardioids represent a powerful platform to mechanistically dissect self-organization, congenital heart defects and serve as a foundation for future translational research.

Clear cell renal carcinoma (ccRCC) is a heterogeneous disease with a variable post-surgical course. To assemble a comprehensive ccRCC tumor microenvironment (TME) atlas, we performed single-cell RNA sequencing (scRNA-seq) of hematopoietic and non-hematopoietic subpopulations from tumor and tumor-adjacent tissue of treatment-naive ccRCC resections. We leveraged the VIPER algorithm to quantitate single-cell protein activity and validated this approach by comparison to flow cytometry. The analysis identified key TME subpopulations, as well as their master regulators and candidate cell-cell interactions, revealing clinically relevant populations, undetectable by gene-expression analysis. Specifically, we uncovered a tumor-specific macrophage subpopulation characterized by upregulation of TREM2/APOE/C1Q, validated by spatially resolved, quantitative multispectral immunofluorescence. In a large clinical validation cohort, these markers were significantly enriched in tumors from patients who recurred following surgery. The study thus identifies TREM2/APOE/C1Q-positive macrophage infiltration as a potential prognostic biomarker for ccRCC recurrence, as well as a candidate therapeutic target.

Monoclonal antibodies (mAbs) are a focus in vaccine and therapeutic design to counteract severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and its variants. Here, we combined B cell sorting with single-cell VDJ and RNA sequencing (RNA-seq) and mAb structures to characterize B cell responses against SARS-CoV-2. We show that the SARS-CoV-2-specific B cell repertoire consists of transcriptionally distinct B cell populations with cells producing potently neutralizing antibodies (nAbs) localized in two clusters that resemble memory and activated B cells. Cryo-electron microscopy structures of selected nAbs from these two clusters complexed with SARS-CoV-2 spike trimers show recognition of various receptor-binding domain (RBD) epitopes. One of these mAbs, BG10-19, locks the spike trimer in a closed conformation to potently neutralize SARS-CoV-2, the recently arising mutants B.1.1.7 and B.1.351, and SARS-CoV and cross-reacts with heterologous RBDs. Together, our results characterize transcriptional differences among SARS-CoV-2-specific B cells and uncover cross-neutralizing Ab targets that will inform immunogen and therapeutic design against coronaviruses.

Defining long-term protective immunity to SARS-CoV-2 is one of the most pressing questions of our time and will require a detailed understanding of potential ways this virus can evolve to escape immune protection. Immune protection will most likely be mediated by antibodies that bind to the viral entry protein, spike (S). Here, we used Phage-DMS, an approach that comprehensively interrogates the effect of all possible mutations on binding to a protein of interest, to define the profile of antibody escape to the SARS-CoV-2 S protein using coronavirus disease 2019 (COVID-19) convalescent plasma. Antibody binding was common in two regions, the fusion peptide and the linker region upstream of the heptad repeat region 2. However, escape mutations were variable within these immunodominant regions. There was also individual variation in less commonly targeted epitopes. This study provides a granular view of potential antibody escape pathways and suggests there will be individual variation in antibody-mediated virus evolution.

Plant species have evolved myriads of solutions, including complex cell type development and regulation, to adapt to dynamic environments. To understand this cellular diversity, we profiled tomato root cell type translatomes. Using xylem differentiation in tomato, examples of functional innovation, repurposing, and conservation of transcription factors are described, relative to the model plant Arabidopsis. Repurposing and innovation of genes are further observed within an exodermis regulatory network and illustrate its function. Comparative translatome analyses of rice, tomato, and Arabidopsis cell populations suggest increased expression conservation of root meristems compared with other homologous populations. In addition, the functions of constitutively expressed genes are more conserved than those of cell type/tissue-enriched genes. These observations suggest that higher order properties of cell type and pan-cell type regulation are evolutionarily conserved between plants and animals.

Gene expression by RNA polymerase II (RNAPII) is tightly controlled by cyclin-dependent kinases (CDKs) at discrete checkpoints during the transcription cycle. The pausing checkpoint following transcription initiation is primarily controlled by CDK9. We discovered that CDK9-mediated, RNAPII-driven transcription is functionally opposed by a protein phosphatase 2A (PP2A) complex that is recruited to transcription sites by the Integrator complex subunit INTS6. PP2A dynamically antagonizes phosphorylation of key CDK9 substrates including DSIF and RNAPII-CTD. Loss of INTS6 results in resistance to tumor cell death mediated by CDK9 inhibition, decreased turnover of CDK9 phospho-substrates, and amplification of acute oncogenic transcriptional responses. Pharmacological PP2A activation synergizes with CDK9 inhibition to kill both leukemic and solid tumor cells, providing therapeutic benefit in vivo. These data demonstrate that fine control of gene expression relies on the balance between kinase and phosphatase activity throughout the transcription cycle, a process dysregulated in cancer that can be exploited therapeutically.

The isocortex and hippocampal formation (HPF) in the mammalian brain play critical roles in perception, cognition, emotion, and learning. We profiled ∼1.3 million cells covering the entire adult mouse isocortex and HPF and derived a transcriptomic cell-type taxonomy revealing a comprehensive repertoire of glutamatergic and GABAergic neuron types. Contrary to the traditional view of HPF as having a simpler cellular organization, we discover a complete set of glutamatergic types in HPF homologous to all major subclasses found in the six-layered isocortex, suggesting that HPF and the isocortex share a common circuit organization. We also identify large-scale continuous and graded variations of cell types along isocortical depth, across the isocortical sheet, and in multiple dimensions in hippocampus and subiculum. Overall, our study establishes a molecular architecture of the mammalian isocortex and hippocampal formation and begins to shed light on its underlying relationship with the development, evolution, connectivity, and function of these two brain structures.

Glycans modify lipids and proteins to mediate inter- and intramolecular interactions across all domains of life. RNA is not thought to be a major target of glycosylation. Here, we challenge this view with evidence that mammals use RNA as a third scaffold for glycosylation. Using a battery of chemical and biochemical approaches, we found that conserved small noncoding RNAs bear sialylated glycans. These "glycoRNAs" were present in multiple cell types and mammalian species, in cultured cells, and in vivo. GlycoRNA assembly depends on canonical N-glycan biosynthetic machinery and results in structures enriched in sialic acid and fucose. Analysis of living cells revealed that the majority of glycoRNAs were present on the cell surface and can interact with anti-dsRNA antibodies and members of the Siglec receptor family. Collectively, these findings suggest the existence of a direct interface between RNA biology and glycobiology, and an expanded role for RNA in extracellular biology.

Since establishment of the first embryonic stem cells (ESCs), in vitro culture of totipotent cells functionally and molecularly comparable with in vivo blastomeres with embryonic and extraembryonic developmental potential has been a challenge. Here we report that spliceosomal repression in mouse ESCs drives a pluripotent-to-totipotent state transition. Using the splicing inhibitor pladienolide B, we achieve stable in vitro culture of totipotent ESCs comparable at molecular levels with 2- and 4-cell blastomeres, which we call totipotent blastomere-like cells (TBLCs). Mouse chimeric assays combined with single-cell RNA sequencing (scRNA-seq) demonstrate that TBLCs have a robust bidirectional developmental capability to generate multiple embryonic and extraembryonic cell lineages. Mechanically, spliceosomal repression causes widespread splicing inhibition of pluripotent genes, whereas totipotent genes, which contain few short introns, are efficiently spliced and transcriptionally activated. Our study provides a means for capturing and maintaining totipotent stem cells.

We identified an emerging severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant by viral whole-genome sequencing of 2,172 nasal/nasopharyngeal swab samples from 44 counties in California, a state in the western United States. Named B.1.427/B.1.429 to denote its two lineages, the variant emerged in May 2020 and increased from 0% to >50% of sequenced cases from September 2020 to January 2021, showing 18.6%-24% increased transmissibility relative to wild-type circulating strains. The variant carries three mutations in the spike protein, including an L452R substitution. We found 2-fold increased B.1.427/B.1.429 viral shedding in vivo and increased L452R pseudovirus infection of cell cultures and lung organoids, albeit decreased relative to pseudoviruses carrying the N501Y mutation common to variants B.1.1.7, B.1.351, and P.1. Antibody neutralization assays revealed 4.0- to 6.7-fold and 2.0-fold decreases in neutralizing titers from convalescent patients and vaccine recipients, respectively. The increased prevalence of a more transmissible variant in California exhibiting decreased antibody neutralization warrants further investigation.

Microbes can play both pathogenic and commensal roles, and it is common to label them as either detrimental or beneficial. However, the lines between good and bad can be blurred. This graphical summary attempts to illustrate the complexity of host-microbe interactions, with outcomes for human health being highly context specific. To view this SnapShot, open or download the PDF.