Integrated analysis of multimodal single-cell data.

The simultaneous measurement of multiple modalities represents an exciting frontier for single-cell genomics and necessitates computational methods that can define cellular states based on multimodal data. Here, we introduce "weighted-nearest neighbor" analysis, an unsupervised framework to learn the relative utility of each data type in each cell, enabling an integrative analysis of multiple modalities. We apply our procedure to a CITE-seq dataset of 211,000 human peripheral blood mononuclear cells (PBMCs) with panels extending to 228 antibodies to construct a multimodal reference atlas of the circulating immune system. Multimodal analysis substantially improves our ability to resolve cell states, allowing us to identify and validate previously unreported lymphoid subpopulations. Moreover, we demonstrate how to leverage this reference to rapidly map new datasets and to interpret immune responses to vaccination and coronavirus disease 2019 (COVID-19). Our approach represents a broadly applicable strategy to analyze single-cell multimodal datasets and to look beyond the transcriptome toward a unified and multimodal definition of cellular identity.

Extremely potent human monoclonal antibodies from COVID-19 convalescent patients.

Human monoclonal antibodies are safe, preventive, and therapeutic tools that can be rapidly developed to help restore the massive health and economic disruption caused by the coronavirus disease 2019 (COVID-19) pandemic. By single-cell sorting 4,277 SARS-CoV-2 spike protein-specific memory B cells from 14 COVID-19 survivors, 453 neutralizing antibodies were identified. The most potent neutralizing antibodies recognized the spike protein receptor-binding domain, followed in potency by antibodies that recognize the S1 domain, the spike protein trimer, and the S2 subunit. Only 1.4% of them neutralized the authentic virus with a potency of 1-10 ng/mL. The most potent monoclonal antibody, engineered to reduce the risk of antibody-dependent enhancement and prolong half-life, neutralized the authentic wild-type virus and emerging variants containing D614G, E484K, and N501Y substitutions. Prophylactic and therapeutic efficacy in the hamster model was observed at 0.25 and 4 mg/kg respectively in absence of Fc functions.

Exon-Mediated Activation of Transcription Starts.

The processing of RNA transcripts from mammalian genes occurs in proximity to their transcription. Here, we describe a phenomenon affecting thousands of genes that we call exon-mediated activation of transcription starts (EMATS), in which the splicing of internal exons impacts promoter choice and the expression level of the gene. We observed that evolutionary gain of internal exons is associated with gain of new transcription start sites (TSSs) nearby and increased gene expression. Inhibiting exon splicing reduced transcription from nearby promoters, and creation of new spliced exons activated transcription from cryptic promoters. The strongest effects occurred for weak promoters located proximal and upstream of efficiently spliced exons. Together, our findings support a model in which splicing recruits transcription machinery locally to influence TSS choice and identify exon gain, loss, and regulatory change as major contributors to the evolution of alternative promoters and gene expression in mammals.

Long-Term Culture Captures Injury-Repair Cycles of Colonic Stem Cells.

The colonic epithelium can undergo multiple rounds of damage and repair, often in response to excessive inflammation. The responsive stem cell that mediates this process is unclear, in part because of a lack of in vitro models that recapitulate key epithelial changes that occur in vivo during damage and repair. Here, we identify a Hopx+ colitis-associated regenerative stem cell (CARSC) population that functionally contributes to mucosal repair in mouse models of colitis. Hopx+ CARSCs, enriched for fetal-like markers, transiently arose from hypertrophic crypts known to facilitate regeneration. Importantly, we established a long-term, self-organizing two-dimensional (2D) epithelial monolayer system to model the regenerative properties and responses of Hopx+ CARSCs. This system can reenact the "homeostasis-injury-regeneration" cycles of epithelial alterations that occur in vivo. Using this system, we found that hypoxia and endoplasmic reticulum stress, insults commonly present in inflammatory bowel diseases, mediated the cyclic switch of cellular status in this process.

LAP2 Proteins Chaperone GLI1 Movement between the Lamina and Chromatin to Regulate Transcription.

Understanding transcription factor navigation through the nucleus remains critical for developing targeted therapeutics. The GLI1 transcription factor must maintain maximal Hedgehog pathway output in basal cell carcinomas (BCCs), and we have previously shown that resistant BCCs increase GLI1 deacetylation through atypical protein kinase Cι/λ (aPKC) and HDAC1. Here we identify a lamina-associated polypeptide 2 (LAP2) isoform-dependent nuclear chaperoning system that regulates GLI1 movement between the nuclear lamina and nucleoplasm to achieve maximal activation. LAP2ß forms a two-site interaction with the GLI1 zinc-finger domain and acetylation site, stabilizing an acetylation-dependent reserve on the inner nuclear membrane (INM). By contrast, the nucleoplasmic LAP2α competes with LAP2ß for GLI1 while scaffolding HDAC1 to deacetylate the secondary binding site. aPKC functions to promote GLI1 association with LAP2α, promoting egress off the INM. GLI1 intranuclear trafficking by LAP2 isoforms represents a powerful signal amplifier in BCCs with implications for zinc finger-based signal transduction and therapeutics.

Mapping the Mouse Cell Atlas by Microwell-Seq.

Single-cell RNA sequencing (scRNA-seq) technologies are poised to reshape the current cell-type classification system. However, a transcriptome-based single-cell atlas has not been achieved for complex mammalian systems. Here, we developed Microwell-seq, a high-throughput and low-cost scRNA-seq platform using simple, inexpensive devices. Using Microwell-seq, we analyzed more than 400,000 single cells covering all of the major mouse organs and constructed a basic scheme for a mouse cell atlas (MCA). We reveal a single-cell hierarchy for many tissues that have not been well characterized previously. We built a web-based "single-cell MCA analysis" pipeline that accurately defines cell types based on single-cell digital expression. Our study demonstrates the wide applicability of the Microwell-seq technology and MCA resource.

Aster Proteins Facilitate Nonvesicular Plasma Membrane to ER Cholesterol Transport in Mammalian Cells.

The mechanisms underlying sterol transport in mammalian cells are poorly understood. In particular, how cholesterol internalized from HDL is made available to the cell for storage or modification is unknown. Here, we describe three ER-resident proteins (Aster-A, -B, -C) that bind cholesterol and facilitate its removal from the plasma membrane. The crystal structure of the central domain of Aster-A broadly resembles the sterol-binding fold of mammalian StARD proteins, but sequence differences in the Aster pocket result in a distinct mode of ligand binding. The Aster N-terminal GRAM domain binds phosphatidylserine and mediates Aster recruitment to plasma membrane-ER contact sites in response to cholesterol accumulation in the plasma membrane. Mice lacking Aster-B are deficient in adrenal cholesterol ester storage and steroidogenesis because of an inability to transport cholesterol from SR-BI to the ER. These findings identify a nonvesicular pathway for plasma membrane to ER sterol trafficking in mammals.

Noise facilitates transcriptional control under dynamic inputs.

Cells must respond sensitively to time-varying inputs in complex signaling environments. To understand how signaling networks process dynamic inputs into gene expression outputs and the role of noise in cellular information processing, we studied the immune pathway NF-κB under periodic cytokine inputs using microfluidic single-cell measurements and stochastic modeling. We find that NF-κB dynamics in fibroblasts synchronize with oscillating TNF signal and become entrained, leading to significantly increased NF-κB oscillation amplitude and mRNA output compared to non-entrained response. Simulations show that intrinsic biochemical noise in individual cells improves NF-κB oscillation and entrainment, whereas cell-to-cell variability in NF-κB natural frequency creates population robustness, together enabling entrainment over a wider range of dynamic inputs. This wide range is confirmed by experiments where entrained cells were measured under all input periods. These results indicate that synergy between oscillation and noise allows cells to achieve efficient gene expression in dynamically changing signaling environments.

Excessive Polyamine Generation in Keratinocytes Promotes Self-RNA Sensing by Dendritic Cells in Psoriasis.

Psoriasis is a chronic inflammatory disease whose etiology is multifactorial. The contributions of cellular metabolism to psoriasis are unclear. Here, we report that interleukin-17 (IL-17) downregulated Protein Phosphatase 6 (PP6) in psoriatic keratinocytes, causing phosphorylation and activation of the transcription factor C/EBP-ß and subsequent generation of arginase-1. Mice lacking Pp6 in keratinocytes were predisposed to psoriasis-like skin inflammation. Accumulation of arginase-1 in Pp6-deficient keratinocytes drove polyamine production from the urea cycle. Polyamines protected self-RNA released by psoriatic keratinocytes from degradation and facilitated the endocytosis of self-RNA by myeloid dendritic cells to promote toll-like receptor-7 (TLR7)-dependent RNA sensing and IL-6 production. An arginase inhibitor improved skin inflammation in murine and non-human primate models of psoriasis. Our findings suggest that urea cycle hyperreactivity and excessive polyamine generation in psoriatic keratinocytes promote self-RNA sensation and PP6 deregulation in keratinocytes is a pivotal event that amplifies the inflammatory circuits in psoriasis.

Autonomous membrane IgE signaling prevents IgE-memory formation.

Aberrant production of IgE antibodies can lead to allergic diseases. Normally, IgE(+) B cells rarely differentiate into memory B cells (Bmem) or long-lived plasma cells (LLPCs), as they only transiently participate in the germinal center (GC), but the mechanism behind this remains elusive. We found that membrane IgE (mIgE) autonomously triggered rapid plasma-cell differentiation and apoptosis independently of antigen or cellular context, predominantly through the mutually independent CD19-PI3K-Akt-IRF4 and BLNK-Jnk/p38 pathways, respectively, and we identified the ectodomains of mIgE as being responsible. Accordingly, deregulated GC IgE(+) B cell proliferation and prolonged IgE production with exaggerated anaphylaxis were observed in CD19- and BLNK-deficient mice. Our findings reveal an autonomous mIgE signaling mechanism that normally prevents IgE(+) Bmem and LLPC formation, providing insights into the molecular pathogenesis of allergic diseases.

Enzymatic Activity of HPGD in Treg Cells Suppresses Tconv Cells to Maintain Adipose Tissue Homeostasis and Prevent Metabolic Dysfunction.

Regulatory T cells (Treg cells) are important for preventing autoimmunity and maintaining tissue homeostasis, but whether Treg cells can adopt tissue- or immune-context-specific suppressive mechanisms is unclear. Here, we found that the enzyme hydroxyprostaglandin dehydrogenase (HPGD), which catabolizes prostaglandin E2 (PGE2) into the metabolite 15-keto PGE2, was highly expressed in Treg cells, particularly those in visceral adipose tissue (VAT). Nuclear receptor peroxisome proliferator-activated receptor-γ (PPARγ)-induced HPGD expression in VAT Treg cells, and consequential Treg-cell-mediated generation of 15-keto PGE2 suppressed conventional T cell activation and proliferation. Conditional deletion of Hpgd in mouse Treg cells resulted in the accumulation of functionally impaired Treg cells specifically in VAT, causing local inflammation and systemic insulin resistance. Consistent with this mechanism, humans with type 2 diabetes showed decreased HPGD expression in Treg cells. These data indicate that HPGD-mediated suppression is a tissue- and context-dependent suppressive mechanism used by Treg cells to maintain adipose tissue homeostasis.

The proliferation-quiescence decision is controlled by a bifurcation in CDK2 activity at mitotic exit.

Tissue homeostasis in metazoans is regulated by transitions of cells between quiescence and proliferation. The hallmark of proliferating populations is progression through the cell cycle, which is driven by cyclin-dependent kinase (CDK) activity. Here, we introduce a live-cell sensor for CDK2 activity and unexpectedly found that proliferating cells bifurcate into two populations as they exit mitosis. Many cells immediately commit to the next cell cycle by building up CDK2 activity from an intermediate level, while other cells lack CDK2 activity and enter a transient state of quiescence. This bifurcation is directly controlled by the CDK inhibitor p21 and is regulated by mitogens during a restriction window at the end of the previous cell cycle. Thus, cells decide at the end of mitosis to either start the next cell cycle by immediately building up CDK2 activity or to enter a transient G0-like state by suppressing CDK2 activity.

MicroRNAs Cause Accelerated Decay of Short-Tailed Target mRNAs.

MicroRNAs (miRNAs) specify the recruitment of deadenylases to mRNA targets. Despite this recruitment, we find that miRNAs have almost no effect on steady-state poly(A)-tail lengths of their targets in mouse fibroblasts, which motivates the acquisition of pre-steady-state measurements of the effects of miRNAs on tail lengths, mRNA levels, and translational efficiencies. Effects on translational efficiency are minimal compared to effects on mRNA levels, even for newly transcribed target mRNAs. Effects on target mRNA levels accumulate as the mRNA population approaches steady state, whereas effects on tail lengths peak for recently transcribed target mRNAs and then subside. Computational modeling of this phenomenon reveals that miRNAs cause not only accelerated deadenylation of their targets but also accelerated decay of short-tailed target molecules. This unanticipated effect of miRNAs largely prevents short-tailed target mRNAs from accumulating despite accelerated target deadenylation. The net result is a nearly imperceptible change to the steady-state tail-length distribution of targeted mRNAs.

The Dynamics of Cytoplasmic mRNA Metabolism.

For all but a few mRNAs, the dynamics of metabolism are unknown. Here, we developed an experimental and analytical framework for examining these dynamics for mRNAs from thousands of genes. mRNAs of mouse fibroblasts exit the nucleus with diverse intragenic and intergenic poly(A)-tail lengths. Once in the cytoplasm, they have a broad (1000-fold) range of deadenylation rate constants, which correspond to cytoplasmic lifetimes. Indeed, with few exceptions, degradation appears to occur primarily through deadenylation-linked mechanisms, with little contribution from either endonucleolytic cleavage or deadenylation-independent decapping. Most mRNA molecules degrade only after their tail lengths fall below 25 nt. Decay rate constants of short-tailed mRNAs vary broadly (1000-fold) and are larger for short-tailed mRNAs that have previously undergone more rapid deadenylation. This coupling helps clear rapidly deadenylated mRNAs, enabling the large range in deadenylation rate constants to impart a similarly large range in stabilities.

The adrenergic-induced ERK3 pathway drives lipolysis and suppresses energy dissipation.

Obesity-induced diabetes affects >400 million people worldwide. Uncontrolled lipolysis (free fatty acid release from adipocytes) can contribute to diabetes and obesity. To identify future therapeutic avenues targeting this pathway, we performed a high-throughput screen and identified the extracellular-regulated kinase 3 (ERK3) as a hit. We demonstrated that ß-adrenergic stimulation stabilizes ERK3, leading to the formation of a complex with the cofactor MAP kinase-activated protein kinase 5 (MK5), thereby driving lipolysis. Mechanistically, we identified a downstream target of the ERK3/MK5 pathway, the transcription factor FOXO1, which promotes the expression of the major lipolytic enzyme ATGL. Finally, we provide evidence that targeted deletion of ERK3 in mouse adipocytes inhibits lipolysis, but elevates energy dissipation, promoting lean phenotype and ameliorating diabetes. Thus, ERK3/MK5 represents a previously unrecognized signaling axis in adipose tissue and an attractive target for future therapies aiming to combat obesity-induced diabetes.

An Antibody Targeting the Fusion Machinery Neutralizes Dual-Tropic Infection and Defines a Site of Vulnerability on Epstein-Barr Virus.

Epstein-Barr virus (EBV) is a causative agent of infectious mononucleosis and is associated with 200,000 new cases of cancer and 140,000 deaths annually. Subunit vaccines against this pathogen have focused on the gp350 glycoprotein and remain unsuccessful. We isolated human antibodies recognizing the EBV fusion machinery (gH/gL and gB) from rare memory B cells. One anti-gH/gL antibody, AMMO1, potently neutralized infection of B cells and epithelial cells, the two major cell types targeted by EBV. We determined a cryo-electron microscopy reconstruction of the gH/gL-gp42-AMMO1 complex and demonstrated that AMMO1 bound to a discontinuous epitope formed by both gH and gL at the Domain-I/Domain-II interface. Integrating structural, biochemical, and infectivity data, we propose that AMMO1 inhibits fusion of the viral and cellular membranes. This work identifies a crucial epitope that may aid in the design of next-generation subunit vaccines against this major public health burden.

Multifunctional Pan-ebolavirus Antibody Recognizes a Site of Broad Vulnerability on the Ebolavirus Glycoprotein.

Ebolaviruses cause severe disease in humans, and identification of monoclonal antibodies (mAbs) that are effective against multiple ebolaviruses are important for therapeutics development. Here we describe a distinct class of broadly neutralizing human mAbs with protective capacity against three ebolaviruses infectious for humans: Ebola (EBOV), Sudan (SUDV), and Bundibugyo (BDBV) viruses. We isolated mAbs from human survivors of ebolavirus disease and identified a potent mAb, EBOV-520, which bound to an epitope in the glycoprotein (GP) base region. EBOV-520 efficiently neutralized EBOV, BDBV, and SUDV and also showed protective capacity in relevant animal models of these infections. EBOV-520 mediated protection principally by direct virus neutralization and exhibited multifunctional properties. This study identified a potent naturally occurring mAb and defined key features of the human antibody response that may contribute to broad protection. This multifunctional mAb and related clones are promising candidates for development as broadly protective pan-ebolavirus therapeutic molecules.

Self-Renewal and Toll-like Receptor Signaling Sustain Exhausted Plasmacytoid Dendritic Cells during Chronic Viral Infection.

Although characterization of T cell exhaustion has unlocked powerful immunotherapies, the mechanisms sustaining adaptations of short-lived innate cells to chronic inflammatory settings remain unknown. During murine chronic viral infection, we found that concerted events in bone marrow and spleen mediated by type I interferon (IFN-I) and Toll-like receptor 7 (TLR7) maintained a pool of functionally exhausted plasmacytoid dendritic cells (pDCs). In the bone marrow, IFN-I compromised the number and the developmental capacity of pDC progenitors, which generated dysfunctional pDCs. Concurrently, exhausted pDCs in the periphery were maintained by self-renewal via IFN-I- and TLR7-induced proliferation of CD4- subsets. On the other hand, pDC functional loss was mediated by TLR7, leading to compromised IFN-I production and resistance to secondary infection. These findings unveil the mechanisms sustaining a self-perpetuating pool of functionally exhausted pDCs and provide a framework for deciphering long-term exhaustion of other short-lived innate cells during chronic inflammation.

Brown remodeling of white adipose tissue by SirT1-dependent deacetylation of Pparγ.

Brown adipose tissue (BAT) can disperse stored energy as heat. Promoting BAT-like features in white adipose (WAT) is an attractive, if elusive, therapeutic approach to staunch the current obesity epidemic. Here we report that gain of function of the NAD-dependent deacetylase SirT1 or loss of function of its endogenous inhibitor Deleted in breast cancer-1 (Dbc1) promote "browning" of WAT by deacetylating peroxisome proliferator-activated receptor (Ppar)-γ on Lys268 and Lys293. SirT1-dependent deacetylation of Lys268 and Lys293 is required to recruit the BAT program coactivator Prdm16 to Pparγ, leading to selective induction of BAT genes and repression of visceral WAT genes associated with insulin resistance. An acetylation-defective Pparγ mutant induces a brown phenotype in white adipocytes, whereas an acetylated mimetic fails to induce "brown" genes but retains the ability to activate "white" genes. We propose that SirT1-dependent Pparγ deacetylation is a form of selective Pparγ modulation of potential therapeutic import.

Design of biologically active binary protein 2D materials.

Ordered two-dimensional arrays such as S-layers1,2 and designed analogues3-5 have intrigued bioengineers6,7, but with the exception of a single lattice formed with flexible linkers8, they are constituted from just one protein component. Materials composed of two components have considerable potential advantages for modulating assembly dynamics and incorporating more complex functionality9-12. Here we describe a computational method to generate co-assembling binary layers by designing rigid interfaces between pairs of dihedral protein building blocks, and use it to design a p6m lattice. The designed array components are soluble at millimolar concentrations, but when combined at nanomolar concentrations, they rapidly assemble into nearly crystalline micrometre-scale arrays nearly identical to the computational design model in vitro and in cells without the need for a two-dimensional support. Because the material is designed from the ground up, the components can be readily functionalized and their symmetry reconfigured, enabling formation of ligand arrays with distinguishable surfaces, which we demonstrate can drive extensive receptor clustering, downstream protein recruitment and signalling. Using atomic force microscopy on supported bilayers and quantitative microscopy on living cells, we show that arrays assembled on membranes have component stoichiometry and structure similar to arrays formed in vitro, and that our material can therefore impose order onto fundamentally disordered substrates such as cell membranes. In contrast to previously characterized cell surface receptor binding assemblies such as antibodies and nanocages, which are rapidly endocytosed, we find that large arrays assembled at the cell surface suppress endocytosis in a tunable manner, with potential therapeutic relevance for extending receptor engagement and immune evasion. Our work provides a foundation for a synthetic cell biology in which multi-protein macroscale materials are designed to modulate cell responses and reshape synthetic and living systems.