Interstitial macrophages are a focus of viral takeover and inflammation in COVID-19 initiation in human lung.

Early stages of deadly respiratory diseases including COVID-19 are challenging to elucidate in humans. Here, we define cellular tropism and transcriptomic effects of SARS-CoV-2 virus by productively infecting healthy human lung tissue and using scRNA-seq to reconstruct the transcriptional program in "infection pseudotime" for individual lung cell types. SARS-CoV-2 predominantly infected activated interstitial macrophages (IMs), which can accumulate thousands of viral RNA molecules, taking over 60% of the cell transcriptome and forming dense viral RNA bodies while inducing host profibrotic (TGFB1, SPP1) and inflammatory (early interferon response, CCL2/7/8/13, CXCL10, and IL6/10) programs and destroying host cell architecture. Infected alveolar macrophages (AMs) showed none of these extreme responses. Spike-dependent viral entry into AMs used ACE2 and Sialoadhesin/CD169, whereas IM entry used DC-SIGN/CD209. These results identify activated IMs as a prominent site of viral takeover, the focus of inflammation and fibrosis, and suggest targeting CD209 to prevent early pathology in COVID-19 pneumonia. This approach can be generalized to any human lung infection and to evaluate therapeutics.

Single-cell transcriptomic analysis of corneal organoids during development.

Corneal organoids are useful tools for disease modeling and tissue transplantation; however, they have not yet been well studied during maturation. We characterized human iPSC-derived corneal organoids at 1, 2, 3, and 4 months of development using single-cell RNA sequencing to determine the cellular heterogeneity at each stage. We found pluripotent cell clusters committed to epithelial cell lineage at 1 month; early corneal epithelial, endothelial, and stromal cell markers at 2 months; keratocytes as the largest cell population at 3 months; and a large epithelial cell population at 4 months. We compared organoid to fetal corneal development at different stages and found that 4-month organoids closely resemble the corneal cellular complexity of the fetal (16 post conception week) and adult cornea. Using RNA velocity trajectory analysis, we found that less differentiated cells appear to give rise to corneal epithelial cells during development.

KRAS(G12D) drives lepidic adenocarcinoma through stem-cell reprogramming.

Many cancers originate from stem or progenitor cells hijacked by somatic mutations that drive replication, exemplified by adenomatous transformation of pulmonary alveolar epithelial type II (AT2) cells1. Here we demonstrate a different scenario: expression of KRAS(G12D) in differentiated AT1 cells reprograms them slowly and asynchronously back into AT2 stem cells that go on to generate indolent tumours. Like human lepidic adenocarcinoma, the tumour cells slowly spread along alveolar walls in a non-destructive manner and have low ERK activity. We find that AT1 and AT2 cells act as distinct cells of origin and manifest divergent responses to concomitant WNT activation and KRAS(G12D) induction, which accelerates AT2-derived but inhibits AT1-derived adenoma proliferation. Augmentation of ERK activity in KRAS(G12D)-induced AT1 cells increases transformation efficiency, proliferation and progression from lepidic to mixed tumour histology. Overall, we have identified a new cell of origin for lung adenocarcinoma, the AT1 cell, which recapitulates features of human lepidic cancer. In so doing, we also uncover a capacity for oncogenic KRAS to reprogram a differentiated and quiescent cell back into its parent stem cell en route to adenomatous transformation. Our work further reveals that irrespective of a given cancer's current molecular profile and driver oncogene, the cell of origin exerts a pervasive and perduring influence on its subsequent behaviour.

Purification and functional characterization of novel human skeletal stem cell lineages.

Human skeletal stem cells (hSSCs) hold tremendous therapeutic potential for developing new clinical strategies to effectively combat congenital and age-related musculoskeletal disorders. Unfortunately, refined methodologies for the proper isolation of bona fide hSSCs and the development of functional assays that accurately recapitulate their physiology within the skeleton have been lacking. Bone marrow-derived mesenchymal stromal cells (BMSCs), commonly used to describe the source of precursors for osteoblasts, chondrocytes, adipocytes and stroma, have held great promise as the basis of various approaches for cell therapy. However, the reproducibility and clinical efficacy of these attempts have been obscured by the heterogeneous nature of BMSCs due to their isolation by plastic adherence techniques. To address these limitations, our group has refined the purity of individual progenitor populations that are encompassed by BMSCs by identifying defined populations of bona fide hSSCs and their downstream progenitors that strictly give rise to skeletally restricted cell lineages. Here, we describe an advanced flow cytometric approach that utilizes an extensive panel of eight cell surface markers to define hSSCs; bone, cartilage and stromal progenitors; and more differentiated unipotent subtypes, including an osteogenic subset and three chondroprogenitors. We provide detailed instructions for the FACS-based isolation of hSSCs from various tissue sources, in vitro and in vivo skeletogenic functional assays, human xenograft mouse models and single-cell RNA sequencing analysis. This application of hSSC isolation can be performed by any researcher with basic skills in biology and flow cytometry within 1-2 days. The downstream functional assays can be performed within a range of 1-2 months.

Purification and characterization of human neural stem and progenitor cells.

The human brain undergoes rapid development at mid-gestation from a pool of neural stem and progenitor cells (NSPCs) that give rise to the neurons, oligodendrocytes, and astrocytes of the mature brain. Functional study of these cell types has been hampered by a lack of precise purification methods. We describe a method for prospectively isolating ten distinct NSPC types from the developing human brain using cell-surface markers. CD24-THY1-/lo cells were enriched for radial glia, which robustly engrafted and differentiated into all three neural lineages in the mouse brain. THY1hi cells marked unipotent oligodendrocyte precursors committed to an oligodendroglial fate, and CD24+THY1-/lo cells marked committed excitatory and inhibitory neuronal lineages. Notably, we identify and functionally characterize a transcriptomically distinct THY1hiEGFRhiPDGFRA- bipotent glial progenitor cell (GPC), which is lineage-restricted to astrocytes and oligodendrocytes, but not to neurons. Our study provides a framework for the functional study of distinct cell types in human neurodevelopment.

Intermittent fasting induces rapid hepatocyte proliferation to restore the hepatostat in the mouse liver.

Nutrient availability fluctuates in most natural populations, forcing organisms to undergo periods of fasting and re-feeding. It is unknown how dietary changes influence liver homeostasis. Here, we show that a switch from ad libitum feeding to intermittent fasting (IF) promotes rapid hepatocyte proliferation. Mechanistically, IF-induced hepatocyte proliferation is driven by the combined action of systemic FGF15 and localized WNT signaling. Hepatocyte proliferation during periods of fasting and re-feeding re-establishes a constant liver-to-body mass ratio, thus maintaining the hepatostat. This study provides the first example of dietary influence on adult hepatocyte proliferation and challenges the widely held view that liver tissue is mostly quiescent unless chemically or mechanically injured.

Deconvoluting complex correlates of COVID-19 severity with a multi-omic pandemic tracking strategy.

The SARS-CoV-2 pandemic has differentially impacted populations across race and ethnicity. A multi-omic approach represents a powerful tool to examine risk across multi-ancestry genomes. We leverage a pandemic tracking strategy in which we sequence viral and host genomes and transcriptomes from nasopharyngeal swabs of 1049 individuals (736 SARS-CoV-2 positive and 313 SARS-CoV-2 negative) and integrate them with digital phenotypes from electronic health records from a diverse catchment area in Northern California. Genome-wide association disaggregated by admixture mapping reveals novel COVID-19-severity-associated regions containing previously reported markers of neurologic, pulmonary and viral disease susceptibility. Phylodynamic tracking of consensus viral genomes reveals no association with disease severity or inferred ancestry. Summary data from multiomic investigation reveals metagenomic and HLA associations with severe COVID-19. The wealth of data available from residual nasopharyngeal swabs in combination with clinical data abstracted automatically at scale highlights a powerful strategy for pandemic tracking, and reveals distinct epidemiologic, genetic, and biological associations for those at the highest risk.

Machine-learning-optimized Cas12a barcoding enables the recovery of single-cell lineages and transcriptional profiles.

The development of CRISPR-based barcoding methods creates an exciting opportunity to understand cellular phylogenies. We present a compact, tunable, high-capacity Cas12a barcoding system called dual acting inverted site array (DAISY). We combined high-throughput screening and machine learning to predict and optimize the 60-bp DAISY barcode sequences. After optimization, top-performing barcodes had ∼10-fold increased capacity relative to the best random-screened designs and performed reliably across diverse cell types. DAISY barcode arrays generated ∼12 bits of entropy and ∼66,000 unique barcodes. Thus, DAISY barcodes-at a fraction of the size of Cas9 barcodes-achieved high-capacity barcoding. We coupled DAISY barcoding with single-cell RNA-seq to recover lineages and gene expression profiles from ∼47,000 human melanoma cells. A single DAISY barcode recovered up to ∼700 lineages from one parental cell. This analysis revealed heritable single-cell gene expression and potential epigenetic modulation of memory gene transcription. Overall, Cas12a DAISY barcoding is an efficient tool for investigating cell-state dynamics.

The Tabula Sapiens: A multiple-organ, single-cell transcriptomic atlas of humans.

Molecular characterization of cell types using single-cell transcriptome sequencing is revolutionizing cell biology and enabling new insights into the physiology of human organs. We created a human reference atlas comprising nearly 500,000 cells from 24 different tissues and organs, many from the same donor. This atlas enabled molecular characterization of more than 400 cell types, their distribution across tissues, and tissue-specific variation in gene expression. Using multiple tissues from a single donor enabled identification of the clonal distribution of T cells between tissues, identification of the tissue-specific mutation rate in B cells, and analysis of the cell cycle state and proliferative potential of shared cell types across tissues. Cell type-specific RNA splicing was discovered and analyzed across tissues within an individual.

uPIC-M: Efficient and Scalable Preparation of Clonal Single Mutant Libraries for High-Throughput Protein Biochemistry.

New high-throughput biochemistry techniques complement selection-based approaches and provide quantitative kinetic and thermodynamic data for thousands of protein variants in parallel. With these advances, library generation rather than data collection has become rate-limiting. Unlike pooled selection approaches, high-throughput biochemistry requires mutant libraries in which individual sequences are rationally designed, efficiently recovered, sequence-validated, and separated from one another, but current strategies are unable to produce these libraries at the needed scale and specificity at reasonable cost. Here, we present a scalable, rapid, and inexpensive approach for creating User-designed Physically Isolated Clonal-Mutant (uPIC-M) libraries that utilizes recent advances in oligo synthesis, high-throughput sample preparation, and next-generation sequencing. To demonstrate uPIC-M, we created a scanning mutant library of SpAP, a 541 amino acid alkaline phosphatase, and recovered 94% of desired mutants in a single iteration. uPIC-M uses commonly available equipment and freely downloadable custom software and can produce a 5000 mutant library at 1/3 the cost and 1/5 the time of traditional techniques.

Aged skeletal stem cells generate an inflammatory degenerative niche.

Loss of skeletal integrity during ageing and disease is associated with an imbalance in the opposing actions of osteoblasts and osteoclasts1. Here we show that intrinsic ageing of skeletal stem cells (SSCs)2 in mice alters signalling in the bone marrow niche and skews the differentiation of bone and blood lineages, leading to fragile bones that regenerate poorly. Functionally, aged SSCs have a decreased bone- and cartilage-forming potential but produce more stromal lineages that express high levels of pro-inflammatory and pro-resorptive cytokines. Single-cell RNA-sequencing studies link the functional loss to a diminished transcriptomic diversity of SSCs in aged mice, which thereby contributes to the transformation of the bone marrow niche. Exposure to a youthful circulation through heterochronic parabiosis or systemic reconstitution with young haematopoietic stem cells did not reverse the diminished osteochondrogenic activity of aged SSCs, or improve bone mass or skeletal healing parameters in aged mice. Conversely, the aged SSC lineage promoted osteoclastic activity and myeloid skewing by haematopoietic stem and progenitor cells, suggesting that the ageing of SSCs is a driver of haematopoietic ageing. Deficient bone regeneration in aged mice could only be returned to youthful levels by applying a combinatorial treatment of BMP2 and a CSF1 antagonist locally to fractures, which reactivated aged SSCs and simultaneously ablated the inflammatory, pro-osteoclastic milieu. Our findings provide mechanistic insights into the complex, multifactorial mechanisms that underlie skeletal ageing and offer prospects for rejuvenating the aged skeletal system.

Distinct skeletal stem cell types orchestrate long bone skeletogenesis.

Skeletal stem and progenitor cell populations are crucial for bone physiology. Characterization of these cell types remains restricted to heterogenous bulk populations with limited information on whether they are unique or overlap with previously characterized cell types. Here we show, through comprehensive functional and single-cell transcriptomic analyses, that postnatal long bones of mice contain at least two types of bone progenitors with bona fide skeletal stem cell (SSC) characteristics. An early osteochondral SSC (ocSSC) facilitates long bone growth and repair, while a second type, a perivascular SSC (pvSSC), co-emerges with long bone marrow and contributes to shape the hematopoietic stem cell niche and regenerative demand. We establish that pvSSCs, but not ocSSCs, are the origin of bone marrow adipose tissue. Lastly, we also provide insight into residual SSC heterogeneity as well as potential crosstalk between the two spatially distinct cell populations. These findings comprehensively address previously unappreciated shortcomings of SSC research.

Light-activated Frizzled7 reveals a permissive role of non-canonical wnt signaling in mesendoderm cell migration.

Non-canonical Wnt signaling plays a central role for coordinated cell polarization and directed migration in metazoan development. While spatiotemporally restricted activation of non-canonical Wnt-signaling drives cell polarization in epithelial tissues, it remains unclear whether such instructive activity is also critical for directed mesenchymal cell migration. Here, we developed a light-activated version of the non-canonical Wnt receptor Frizzled 7 (Fz7) to analyze how restricted activation of non-canonical Wnt signaling affects directed anterior axial mesendoderm (prechordal plate, ppl) cell migration within the zebrafish gastrula. We found that Fz7 signaling is required for ppl cell protrusion formation and migration and that spatiotemporally restricted ectopic activation is capable of redirecting their migration. Finally, we show that uniform activation of Fz7 signaling in ppl cells fully rescues defective directed cell migration in fz7 mutant embryos. Together, our findings reveal that in contrast to the situation in epithelial cells, non-canonical Wnt signaling functions permissively rather than instructively in directed mesenchymal cell migration during gastrulation.

Optical functionalization of human Class A orphan G-protein-coupled receptors.

G-protein-coupled receptors (GPCRs) form the largest receptor family, relay environmental stimuli to changes in cell behavior and represent prime drug targets. Many GPCRs are classified as orphan receptors because of the limited knowledge on their ligands and coupling to cellular signaling machineries. Here, we engineer a library of 63 chimeric receptors that contain the signaling domains of human orphan and understudied GPCRs functionally linked to the light-sensing domain of rhodopsin. Upon stimulation with visible light, we identify activation of canonical cell signaling pathways, including cAMP-, Ca2+-, MAPK/ERK-, and Rho-dependent pathways, downstream of the engineered receptors. For the human pseudogene GPR33, we resurrect a signaling function that supports its hypothesized role as a pathogen entry site. These results demonstrate that substituting unknown chemical activators with a light switch can reveal information about protein function and provide an optically controlled protein library for exploring the physiology and therapeutic potential of understudied GPCRs.