Increased epithelial mTORC1 activity in chronic rhinosinusitis with nasal polyps

Background

The airway epithelium plays a central role in the pathogenesis of chronic respiratory diseases such as asthma and chronic rhinosinusitis with nasal polyps (CRSwNP), but the mechanisms by which airway epithelial cells (EpCs) maintain inflammation are poorly understood.

Objective

We hypothesized that transcriptomic assessment of sorted airway EpCs across the spectrum of differentiation would allow us to define mechanisms by which EpCs perpetuate airway inflammation.

Methods

Ethmoid sinus EpCs from adult patients with CRS were sorted into 3 subsets, bulk RNA sequenced, and analyzed for differentially expressed genes and pathways. Single cell RNA-seq (scRNA-seq) datasets from eosinophilic and non-eosinophilic CRSwNP and bulk RNA-seq of EpCs from mild/moderate and severe asthma were assessed. Immunofluorescent staining and ex vivo functional analysis of sinus EpCs were used to validate our findings.

Results

Analysis within and across purified EpC subsets revealed an enrichment in glycolytic programming in CRSwNP vs CRSsNP. Correlation analysis identified mammalian target of rapamycin complex 1 (mTORC1) as a potential regulator of the glycolytic program and identified EpC expression of cytokines and wound healing genes as potential sequelae. mTORC1 activity was upregulated in CRSwNP, and ex vivo inhibition demonstrated that mTOR is critical for EpC generation of CXCL8, IL-33, and CXCL2. Across patient samples, the degree of glycolytic activity was associated with T2 inflammation in CRSwNP, and with both T2 and non-T2 inflammation in severe asthma.

Conclusions

Together, these findings highlight a metabolic axis required to support epithelial generation of cytokines critical to both chronic T2 and non-T2 inflammation in CRSwNP and asthma.

Single cell RNA-sequencing of human precision-cut lung slices: A novel approach to study the effect of vaping and viral infection on lung health

Vaping is an increasing health threat in the US and worldwide. The damaging impact of vaping on the human distal lung has been highlighted by the recent epidemic of electronic cigarette or vaping use-associated lung injury (EVALI). The pathogenesis of EVALI remains incompletely understood, due to a paucity of models that recapitulate the structural and functional complexity of the human distal lung and the still poorly defined culprit exposures to vaping products and respiratory viral infections. Our aim was to establish the feasibility of using single cell RNA-sequencing (scRNA-seq) technology in human precision-cut lung slices (PCLS) as a more physiologically relevant model to better understand how vaping regulates the antiviral and pro-inflammatory response to influenza A virus infection. Normal healthy donor PCLS were treated with vaping extract and influenza A viruses for scRNA-seq analysis. Vaping extract augmented host antiviral and pro-inflammatory responses in structural cells such as lung epithelial cells and fibroblasts, as well as in immune cells such as macrophages and monocytes. Our findings suggest that human distal lung slice model is useful to study the heterogeneous responses of immune and structural cells under EVALI conditions, such as vaping and respiratory viral infection.

Challenges and opportunities to computationally deconvolve heterogeneous tissue with varying cell sizes using single cell RNA-sequencing datasets

Deconvolution of cell mixtures in "bulk" transcriptomic samples from homogenate human tissue is important for understanding the pathologies of diseases. However, several experimental and computational challenges remain in developing and implementing transcriptomics-based deconvolution approaches, especially those using a single cell/nuclei RNA-seq reference atlas, which are becoming rapidly available across many tissues. Notably, deconvolution algorithms are frequently developed using samples from tissues with similar cell sizes. However, brain tissue or immune cell populations have cell types with substantially different cell sizes, total mRNA expression, and transcriptional activity. When existing deconvolution approaches are applied to these tissues, these systematic differences in cell sizes and transcriptomic activity confound accurate cell proportion estimates and instead may quantify total mRNA content. Furthermore, there is a lack of standard reference atlases and computational approaches to facilitate integrative analyses, including not only bulk and single cell/nuclei RNA-seq data, but also new data modalities from spatial -omic or imaging approaches. New multi-assay datasets need to be collected with orthogonal data types generated from the same tissue block and the same individual, to serve as a "gold standard" for evaluating new and existing deconvolution methods. Below, we discuss these key challenges and how they can be addressed with the acquisition of new datasets and approaches to analysis.

Dynamics of Human Hematopoietic Stem and Progenitor Cell Differentiation to the Erythroid Lineage

Erythropoiesis, the development of erythrocytes from hematopoietic stem cells, occurs through four phases: erythroid progenitor development, early erythropoiesis, terminal erythroid differentiation (TED), and maturation. According to the classical model that is based on immunophenotypic profiles of cell populations, each of these phases comprises multiple differentiation states that arise in a hierarchical manner. After segregation of lymphoid potential, erythroid priming begins during progenitor development and progresses through progenitor cell types that have multilineage potential. Complete separation of the erythroid lineage is achieved during early erythropoiesis with the formation of unipotent erythroid progenitors: burst forming unit-erythroid and colony forming unit-erythroid. These erythroid committed progenitors undergo TED and maturation, which involves expulsion of the nucleus and remodeling to form functional biconcave, hemoglobin-filled erythrocytes. In the last decade or so, many studies employing advanced techniques such as single cell RNA-sequencing (scRNA-seq) as well as the conventional methods, including colony forming cell assays and immunophenotyping, have revealed heterogeneity within the stem, progenitor, and erythroblast stages, and uncovered alternate paths for segregation of erythroid lineage potential. In this review, we provide an in-depth account of immunophenotypic profiles of all cell types within erythropoiesis, highlight studies that demonstrate heterogeneous erythroid stages, and describe deviations to the classical model of erythropoiesis. Overall, although scRNA-seq approaches have provided new insights, flow cytometry remains relevant and is the primary method for validation of novel immunophenotypes.

The myonuclear domain in adult skeletal muscle fibres: past, present and future

Most cells in the body are mononuclear whereas skeletal muscle fibres are uniquely multinuclear. The nuclei of muscle fibres (myonuclei) are usually situated peripherally which complicates the equitable distribution of gene products. Myonuclear abundance can also change under conditions such as hypertrophy and atrophy. Specialised zones in muscle fibres have different functions and thus distinct synthetic demands from myonuclei. The complex structure and regulatory requirements of multinuclear muscle cells understandably led to the hypothesis that myonuclei govern defined 'domains' to maintain homeostasis and facilitate adaptation. The purpose of this review is to provide historical context for the myonuclear domain and evaluate its veracity with respect to mRNA and protein distribution resulting from myonuclear transcription. We synthesise insights from past and current in vitro and in vivo genetically modified models for studying the myonuclear domain under dynamic conditions. We also cover the most contemporary knowledge on mRNA and protein transport in muscle cells. Insights from emerging technologies such as single myonuclear RNA-sequencing further inform our discussion of the myonuclear domain. We broadly conclude: (1) the myonuclear domain can be flexible during muscle fibre growth and atrophy, (2) the mechanisms and role of myonuclear loss and motility deserve further consideration, (3) mRNA in muscle is actively transported via microtubules and locally restricted, but proteins may travel far from a myonucleus of origin and (4) myonuclear transcriptional specialisation extends beyond the classic neuromuscular and myotendinous populations. A deeper understanding of the myonuclear domain in muscle may promote effective therapies for ageing and disease.

A comprehensive single cell data analysis of lymphoblastoid cells reveals the role of super-enhancers in maintaining EBV latency

We probed the lifecycle of Epstein-Barr virus (EBV) on a cell-by-cell basis using single cell RNA sequencing (scRNA-seq) data from nine publicly available lymphoblastoid cell lines (LCLs). While the majority of LCLs comprised cells containing EBV in the latent phase, two other clusters of cells were clearly evident and were distinguished by distinct expression of host and viral genes. Notably, both were high expressors of EBV LMP1/BNLF2 and BZLF1 compared to another cluster that expressed neither gene. The two novel clusters differed from each other in their expression of EBV lytic genes, including glycoprotein gene GP350. The first cluster, comprising GP350- LMP1hi cells, expressed high levels of HIF1A and was transcriptionally regulated by HIF1-α. Treatment of LCLs with Pevonedistat, a drug that enhances HIF1-α signaling, markedly induced this cluster. The second cluster, containing GP350+ LMP1hi cells, expressed EBV lytic genes. Host genes that are controlled by super-enhancers (SEs), such as transcription factors MYC and IRF4, had the lowest expression in this cluster. Functionally, the expression of genes regulated by MYC and IRF4 in GP350+ LMP1hi cells were lower compared to other cells. Indeed, induction of EBV lytic reactivation in EBV+ AKATA reduced the expression of these SE-regulated genes. Furthermore, CRISPR-mediated perturbation of the MYC or IRF4 SEs in LCLs induced the lytic EBV gene expression, suggesting that host SEs and/or SE target genes are required for maintenance of EBV latency. Collectively, our study revealed EBV-associated heterogeneity among LCLs that may have functional consequence on host and viral biology.

Lactate metabolism coordinates macrophage response and regeneration in zebrafish

Rationale: Macrophages are multifunctional cells with a pivotal role on tissue development, homeostasis and regeneration. Indeed, in response to tissue injury and the ensuing regeneration process, macrophages are challenged and undergo massive metabolic adaptations and changes. However, the control of this metabolic reprogramming by macrophage microenvironment has never been deciphered in vivo. Methods: In this study, we used zebrafish model and caudal fin resection as a robust regeneration system. We explored specific changes in gene expression after tissue amputation via single-cell RNA sequencing analysis and whole-tissue transcriptomic analysis. Based on the identification of key modifications, we confirmed the role of the lactate pathway in macrophage response and fin regeneration, through the combination of chemical and genetic inhibitors of this pathway. Results: Single cell RNA sequencing revealed the upregulation of different genes associated with glycolysis and lactate metabolism in macrophages, upon fin regeneration. Hence, using chemical inhibitors of the LDH enzyme, we confirmed the role of lactate in macrophage recruitment and polarization, to promote a pro-inflammatory phenotype and enhance fin regeneration. The genetic modulation of monocarboxylate transporters illustrated a complex regulation of lactate levels, based on both intracellular and extracellular supplies. Commonly, the different sources of lactate resulted in macrophage activation with an increased expression level of inflammatory cytokines such as TNFa during the first 24 hours of regeneration. Transcriptomic analyses confirmed that lactate induced a global modification of gene expression in macrophages. Conclusion: Altogether, our findings highlight the crucial role of lactate at the onset of macrophage differentiation toward a pro-inflammatory phenotype. The deep modifications of macrophage phenotype mediated by lactate and downstream effectors play a key role to coordinate inflammatory response and tissue regeneration.

Transcriptomes of MPO-deficient patients with generalized pustular psoriasis reveals expansion of CD4+ cytotoxic T cells and an involvement of the complement system

Generalized pustular psoriasis (GPP) is a severe psoriatic subtype characterized by epidermal neutrophil infiltration. Although variants in IL36RN and MPO have been shown to affect immune cells, a systematic analysis of neutrophils and peripheral blood mononuclear cells (PBMCs) subsets and their differential gene expression dependent on MPO genotypes was not performed yet. We assessed transcriptomes of MPO-deficient patients using single cell RNA-sequencing (scRNAseq) of PBMCs and RNA-sequencing of neutrophils in stable disease state. Cell type annotation by multimodal reference mapping of scRNAseq data was verified by flow cytometry of surface and intracellular markers; proportions of CD4⁺ cytotoxic T-lymphocytes (CTLs) and other CD4⁺ effector cells were increased in GPP, while frequencies of naïve CD4⁺ T cells were significantly lower. The expression of FGFBP2 marking CD4⁺ CTLs and CD8⁺ effector memory T-cells (TEMs) was elevated in GPP patients with disease-contributing variants compared to non-carriers (p=0.0015). In neutrophils, differentially expressed genes (DEGs) were significantly enriched in genes of the classical complement activation pathway. Future studies assessing affected cell-types and pathways will show their contribution to GPP's pathogenesis, and indicate whether findings can be transferred to the acute epidermal situation and whether depletion or inactivation of CD4⁺ CTLs may be a reasonable therapeutic approach.

Transcriptomic and Single-Cell Analysis Reveals Regulatory Networks and Cellular Heterogeneity in Mouse Primary Sjögren's Syndrome Salivary Glands

Sjögren’s Syndrome (SS) is a chronic autoimmune disease of unknown etiology which primarily affects the salivary and lacrimal glands resulting in the loss of secretory function. Treatment options for SS have been hampered due to the lack of a better understanding of the underlying gene regulatory circuitry and the interplay between the myriad pathological cellular states that contribute to salivary gland dysfunction. To better elucidate the molecular nature of SS, we have performed RNA-sequencing analysis of the submandibular glands (SMG) of a well-established primary Sjögren’s Syndrome (pSS) mouse model. Our comprehensive examination of global gene expression and comparative analyses with additional SS mouse models and human datasets, have identified a number of important pathways and regulatory networks that are relevant in SS pathobiology. To complement these studies, we have performed single-cell RNA sequencing to examine and identify the molecular and cellular heterogeneity of the diseased cell populations of the mouse SMG. Interrogation of the single-cell transcriptomes has shed light on the diversity of immune cells that are dysregulated in SS and importantly, revealed an activated state of the salivary gland epithelial cells that contribute to the global immune mediated responses. Overall, our broad studies have not only revealed key pathways, mediators and new biomarkers, but have also uncovered the complex nature of the cellular populations in the SMG that are likely to drive the progression of SS. These newly discovered insights into the underlying molecular mechanisms and cellular states of SS will better inform targeted therapeutic discoveries.

Single-Cell RNA Sequencing Reveals Heterogeneity and Functional Diversity of Lymphatic Endothelial Cells

Lymphatic endothelial cells (LECs) line the lymphatic vasculature and play a central role in the immune response. LECs have abilities to regulate immune transport, to promote immune cell survival, and to cross present antigens to dendritic cells. Single-cell RNA sequencing (scRNA) technology has accelerated new discoveries in the field of lymphatic vascular biology. This review will summarize these new findings in regard to embryonic development, LEC heterogeneity with associated functional diversity, and interactions with other cells. Depending on the organ, location in the lymphatic vascular tree, and micro-environmental conditions, LECs feature unique properties and tasks. Furthermore, adjacent stromal cells need the support of LECs for fulfilling their tasks in the immune response, such as immune cell transport and antigen presentation. Although aberrant lymphatic vasculature has been observed in a number of chronic inflammatory diseases, the knowledge on LEC heterogeneity and functional diversity in these diseases is limited. Combining scRNA sequencing data with imaging and more in-depth functional experiments will advance our knowledge of LECs in health and disease. Building the case, the LEC could be put forward as a new therapeutic target in chronic inflammatory diseases, counterweighting the current immune-cell focused therapies.

Integrative analysis of the human brain mural cell transcriptome

Brain mural cells, including pericytes and vascular smooth muscle cells, are important for vascular development, blood-brain barrier function, and neurovascular coupling, but the molecular characteristics of human brain mural cells are incompletely characterized. Single cell RNA-sequencing (scRNA-seq) is increasingly being applied to assess cellular diversity in the human brain, but the scarcity of mural cells in whole brain samples has limited their molecular profiling. Here, we leverage the combined power of multiple independent human brain scRNA-seq datasets to build a transcriptomic database of human brain mural cells. We use this combined dataset to determine human-mouse species differences in mural cell transcriptomes, culture-induced dedifferentiation of human brain pericytes, and human mural cell organotypicity, with several key findings validated by RNA fluorescence in situ hybridization. Together, this work improves knowledge regarding the molecular constituents of human brain mural cells, serves as a resource for hypothesis generation in understanding brain mural cell function, and will facilitate comparative evaluation of animal and in vitro models.

Ribosome Proteins Represented by RPL27A Mark the Development and Metastasis of Triple-Negative Breast Cancer in Mouse and Human

Triple-negative breast cancer (TNBC) is known to have a poor prognosis and limited treatment options. The lack of targeted therapies and poor prognosis of patients with TNBC have made it urgent to discover novel critical diagnosis and therapeutic targets in the TNBC field. Here, in the current study, we integrated the single-cell RNA-sequencing (scRNA-seq) data from four normal mouse mammary tissues and four mouse breast tumors. Comparative analysis was conducted to identify the gene profiles of normal epithelial cells and cancer cells at different models. Surprisingly, two ribosomal protein genes, Rpl27a and Rpl15, were significantly upregulated in the cancer cells in all the TNBC models. Next, we accessed the scRNA-seq data from human primary and metastatic TNBC tissues, and comparative analysis revealed gene profiles of human primary and metastatic TNBC cancer cells. Ribosomal protein genes, represented by RPL27A and RPL15, showed significantly upregulated expression in metastatic TNBC cancer cells. Pathway analysis on the upregulated genes of the metastatic TNBC cancer cells identified the key regulators and signaling pathways that were driving the metastasis of the TNBC cancer cells. Specifically, EIF2 signaling was significantly activated, and major member genes of this signaling pathway were upregulated. In vitro study revealed that targeting RPL27A or EIF2 signaling in a TNBC cell line, MDA-MB-231, significantly reduced cell migration and invasion. Altogether, these data suggested that the RPL27A gene is conducting critical functions in TNBC cancer development and metastasis and is a potential therapeutic target for TNBC.

Marrow adipogenic lineage precursor: A new cellular component of marrow adipose tissue

Bone marrow mesenchymal stromal cells are a highly heterogenic cell population containing mesenchymal stem cells as well as other cell types. With the advance of single cell transcriptome analysis, several recent reports identified a prominent subpopulation of mesenchymal stromal cells that specifically express adipocyte markers but do not contain lipid droplets. We name this cell type marrow adipogenic lineage precursor, MALP, and consider it as a major cellular component of marrow adipose tissue. Here, we review the discovery of MALPs and summarize their unique features and regulatory roles in bone. We further discuss how these findings advance our understanding of bone remodeling, mesenchymal niche regulation of hematopoiesis, and marrow vasculature maintenance.

How to Get Started with Single Cell RNA Sequencing Data Analysis

Over the last 5 years, single cell methods have enabled the monitoring of gene and protein expression, genetic, and epigenetic changes in thousands of individual cells in a single experiment. With the improved measurement and the decreasing cost of the reactions and sequencing, the size of these datasets is increasing rapidly. The critical bottleneck remains the analysis of the wealth of information generated by single cell experiments. In this review, we give a simplified overview of the analysis pipelines, as they are typically used in the field today. We aim to enable researchers starting out in single cell analysis to gain an overview of challenges and the most commonly used analytical tools. In addition, we hope to empower others to gain an understanding of how typical readouts from single cell datasets are presented in the published literature.

Single Cell Transcriptomic Analyses Reveal the Impact of bHLH Factors on Human Retinal Organoid Development

The developing retina expresses multiple bHLH transcription factors. Their precise functions and interactions in uncommitted retinal progenitors remain to be fully elucidated. Here, we investigate the roles of bHLH factors ATOH7 and Neurog2 in human ES cell-derived retinal organoids. Single cell transcriptome analyses identify three states of proliferating retinal progenitors: pre-neurogenic, neurogenic, and cell cycle-exiting progenitors. Each shows different expression profile of bHLH factors. The cell cycle-exiting progenitors feed into a postmitotic heterozygous neuroblast pool that gives rise to early born neuronal lineages. Elevating ATOH7 or Neurog2 expression accelerates the transition from the pre-neurogenic to the neurogenic state, and expands the exiting progenitor and neuroblast populations. In addition, ATOH7 and Neurog2 significantly, yet differentially, enhance retinal ganglion cell and cone photoreceptor production. Moreover, single cell transcriptome analyses reveal that ATOH7 and Neurog2 each assert positive autoregulation, and both suppress key bHLH factors associated with the pre-neurogenic and states and elevate bHLH factors expressed by exiting progenitors and differentiating neuroblasts. This study thus provides novel insight regarding how ATOH7 and Neurog2 impact human retinal progenitor behaviors and neuroblast fate choices.

Single-Cell Profiling and SCOPE-Seq Reveal Lineage Dynamics of Adult Ventricular-Subventricular Zone Neurogenesis and NOTUM as a Key Regulator

In the adult ventricular-subventricular zone (V-SVZ), neural stem cells (NSCs) generate new olfactory bulb (OB) neurons and glia throughout life. To map adult neuronal lineage progression, we profiled >56,000 V-SVZ and OB cells by single-cell RNA sequencing (scRNA-seq). Our analyses reveal the molecular diversity of OB neurons, including fate-mapped neurons, lineage progression dynamics, and an NSC intermediate enriched for Notum, which encodes a secreted WNT antagonist. SCOPE-seq technology, which links live-cell imaging with scRNA-seq, uncovers cell-size transitions during NSC differentiation and preferential NOTUM binding to proliferating neuronal precursors. Consistently, application of NOTUM protein in slice cultures and pharmacological inhibition of NOTUM in slice cultures and in vivo demonstrated that NOTUM negatively regulates V-SVZ proliferation. Timely, context-dependent neurogenesis demands adaptive signaling among neighboring progenitors. Our findings highlight a critical regulatory state during NSC activation marked by NOTUM, which attenuates WNT-stimulated proliferation in NSC progeny.

In Pursuit of Adult Progenitors of Thymic Epithelial Cells

Peripheral T cells capable of discriminating between self and non-self antigens are major components of a robust adaptive immune system. The development of self-tolerant T cells is orchestrated by thymic epithelial cells (TECs), which are localized in the thymic cortex (cortical TECs, cTECs) and medulla (medullary TECs, mTECs). cTECs and mTECs are essential for differentiation, proliferation, and positive and negative selection of thymocytes. Recent advances in single-cell RNA-sequencing technology have revealed a previously unknown degree of TEC heterogeneity, but we still lack a clear picture of the identity of TEC progenitors in the adult thymus. In this review, we describe both earlier and recent findings that shed light on features of these elusive adult progenitors in the context of tissue homeostasis, as well as recovery from stress-induced thymic atrophy.

Isolation of Plant Root Nuclei for Single Cell RNA Sequencing

The characterization of the transcriptional similarities and differences existing between plant cells and cell types is important to better understand the biology of each cell composing the plant, to reveal new molecular mechanisms controlling gene activity, and to ultimately implement meaningful strategies to enhance plant cell biology. To gain a deeper understanding of the regulation of plant gene activity, the individual transcriptome of each plant cell needs to be established. Until recently, single cell approaches were mostly limited to bulk transcriptomic studies on selected cell types. Accessing specific cell types required the development of labor-intensive strategies. Recently, single cell sequencing strategies were successfully applied on isolated Arabidopsis thaliana root protoplasts. However, this strategy relies on the successful isolation of viable protoplasts upon the optimization of the enzymatic cocktails required to digest the cell wall and on the compatibility of fragile plant protoplasts with the use of microfluidic systems to generate single cell transcriptomic libraries. To overcome these difficulties, we present a simple and fast alternative strategy: the isolation and use of plant nuclei to access meaningful transcriptomic information from plant cells. This protocol was specifically developed to enable the use of the plant nuclei with 10× Genomics' Chromium technology partitions technology. Briefly, the plant nuclei are released from the root by chopping into a nuclei isolation buffer before purification by filtration then nuclei sorting. Upon sorting, the nuclei are resuspended in a low divalent ion buffer compatible with the Chromium technology in order to create single nuclei ribonucleic acid-sequencing libraries (sNucRNA-seq). © 2020 Wiley Periodicals LLC. Basic Protocol 1: Arabidopsis seed sterilization and planting Basic Protocol 2: Nuclei isolation from Arabidopsis roots Basic Protocol 3: Fluorescent-activated nuclei sorting (FANS) purification Support Protocol: Estimation of nuclei density using Countess II automated cell counter Alternate Protocol 1: Proper growth conditions for Medicago truncatula and Sorghum bicolor Alternate Protocol 2: Estimation of nuclei density using sNucRNA-seq technology.

FEATS: feature selection-based clustering of single-cell RNA-seq data

MOTIVATION

Advances in next-generation sequencing have made it possible to carry out transcriptomic studies at single-cell resolution and generate vast amounts of single-cell RNA sequencing (RNA-seq) data rapidly. Thus, tools to analyze this data need to evolve as well as to improve accuracy and efficiency.

RESULTS

We present FEATS, a Python software package, that performs clustering on single-cell RNA-seq data. FEATS is capable of performing multiple tasks such as estimating the number of clusters, conducting outlier detection and integrating data from various experiments. We develop a univariate feature selection-based approach for clustering, which involves the selection of top informative features to improve clustering performance. This is motivated by the fact that cell types are often manually determined using the expression of only a few known marker genes. On a variety of single-cell RNA-seq datasets, FEATS gives superior performance compared with the current tools, in terms of adjusted Rand index and estimating the number of clusters. It achieves a 22% improvement in clustering and more accurately estimates the number of clusters when compared with other tools. In addition to cluster estimation, FEATS also performs outlier detection and data integration while giving an excellent computational performance. Thus, FEATS is a comprehensive clustering tool capable of addressing the challenges during the clustering of single-cell RNA-seq data.

AVAILABILITY

The installation instructions and documentation of FEATS is available at https://edwinv87.github.io/feats/.

SUPPLEMENTARY DATA

Supplementary data are available online at https://academic.oup.com/bib.

Human Primordial Germ Cells Are Specified from Lineage-Primed Progenitors

In vitro gametogenesis is the process of making germline cells from human pluripotent stem cells. The foundation of this model is the quality of the first progenitors called primordial germ cells (PGCs), which in vivo are specified during the peri-implantation window of human development. Here, we show that human PGC (hPGC) specification begins at day 12 post-fertilization. Using single-cell RNA sequencing of hPGC-like cells (hPGCLCs) differentiated from pluripotent stem cells, we discovered that hPGCLC specification involves resetting pluripotency toward a transitional state with shared characteristics between naive and primed pluripotency, followed by differentiation into lineage-primed TFAP2A⁺ progenitors. Applying the germline trajectory to TFAP2C mutants reveals that TFAP2C functions in the TFAP2A⁺ progenitors upstream of PRDM1 to regulate the expression of SOX17. This serves to protect hPGCLCs from crossing the Weismann’s barrier to adopt somatic cell fates and, therefore, is an essential mechanism for successfully initiating in vitro gametogenesis.

Using genetics, genomics, and single-cell RNA-seq, Chen et al. characterize the human germline trajectory, revealing two pluripotent cell transitions during primordial germ cell specification. They reveal the identity of primordial germ cell progenitors and show that TFAP2C prevents gastrulation and amnion-like fate at the point of primordial germ cell specification.

An Atlas of Vagal Sensory Neurons and Their Molecular Specialization

Sensory functions of the vagus nerve are critical for conscious perceptions and for monitoring visceral functions in the cardio-pulmonary and gastrointestinal systems. Here, we present a comprehensive identification, classification, and validation of the neuron types in the neural crest (jugular) and placode (nodose) derived vagal ganglia by single-cell RNA sequencing (scRNA-seq) transcriptomic analysis. Our results reveal major differences between neurons derived from different embryonic origins. Jugular neurons exhibit fundamental similarities to the somatosensory spinal neurons, including major types, such as C-low threshold mechanoreceptors (C-LTMRs), A-LTMRs, Aδ-nociceptors, and cold-, and mechano-heat C-nociceptors. In contrast, the nodose ganglion contains 18 distinct types dedicated to surveying the physiological state of the internal body. Our results reveal a vast diversity of vagal neuron types, including many previously unanticipated types, as well as proposed types that are consistent with chemoreceptors, nutrient detectors, baroreceptors, and stretch and volume mechanoreceptors of the respiratory, gastrointestinal, and cardiovascular systems.

•

A comprehensive molecular identification of neuronal types in vagal ganglion complex

•

Prdm12⁺ jugular ganglion neurons share features with spinal somatosensory neurons

•

Phox2b⁺ viscerosensory nodose neurons are molecularly versatile and highly specialized

•

Nodose neuron types are consistent with chemo-, baro-, stretch-, tension-, and volume-sensors

Integration and reanalysis of transcriptomics and methylomics data derived from blood and testis tissue of men with 47,XXY Klinefelter syndrome indicates the primary involvement of Sertoli cells in the testicular pathogenesis

Klinefelter syndrome (KS; 47,XXY) is the most common sex chromosomal anomaly and causes a multitude of symptoms. Often the most noticeable symptom is infertility caused by azoospermia with testicular histology showing hyalinization of tubules, germ cells loss, and Leydig cell hyperplasia. The germ cell loss begins early in life leading to partial hyalinization of the testis at puberty, but the mechanistic drivers behind this remain poorly understood. In this systematic review, we summarize the current knowledge on developmental changes in the cellularity of KS gonads supplemented by a comparative analysis of the fetal and adult gonadal transcriptome, and blood transcriptome and methylome of men with KS. We identified a high fraction of upregulated genes that escape X-chromosome inactivation, thus supporting previous hypotheses that these are the main drivers of the testicular phenotype in KS. Enrichment analysis showed overrepresentation of genes from the X- and Y-chromosome and testicular transcription factors. Furthermore, by re-evaluation of recent single cell RNA-sequencing data originating from adult KS testis, we found novel evidence that the Sertoli cell is the most affected cell type. Our results are consistent with disturbed cross-talk between somatic and germ cells in the KS testis, and with X-escapee genes acting as mediators.

A Human Skeletal Muscle Atlas Identifies the Trajectories of Stem and Progenitor Cells across Development and from Human Pluripotent Stem Cells

The developmental trajectory of human skeletal myogenesis and the transition between progenitor and stem cell states are unclear. We used single-cell RNA sequencing to profile human skeletal muscle tissues from embryonic, fetal, and postnatal stages. In silico, we identified myogenic as well as other cell types and constructed a "roadmap" of human skeletal muscle ontogeny across development. In a similar fashion, we also profiled the heterogeneous cell cultures generated from multiple human pluripotent stem cell (hPSC) myogenic differentiation protocols and mapped hPSC-derived myogenic progenitors to an embryonic-to-fetal transition period. We found differentially enriched biological processes and discovered co-regulated gene networks and transcription factors present at distinct myogenic stages. This work serves as a resource for advancing our knowledge of human myogenesis. It also provides a tool for a better understanding of hPSC-derived myogenic progenitors for translational applications in skeletal muscle-based regenerative medicine.

Transcriptional profiling of lung cell populations in idiopathic pulmonary arterial hypertension

Despite recent improvements in management of idiopathic pulmonary arterial hypertension, mortality remains high. Understanding the alterations in the transcriptome–phenotype of the key lung cells involved could provide insight into the drivers of pathogenesis. In this study, we examined differential gene expression of cell types implicated in idiopathic pulmonary arterial hypertension from lung explants of patients with idiopathic pulmonary arterial hypertension compared to control lungs. After tissue digestion, we analyzed all cells from three idiopathic pulmonary arterial hypertension and six control lungs using droplet-based single cell RNA-sequencing. After dimensional reduction by t-stochastic neighbor embedding, we compared the transcriptomes of endothelial cells, pericyte/smooth muscle cells, fibroblasts, and macrophage clusters, examining differential gene expression and pathways implicated by analysis of Gene Ontology Enrichment. We found that endothelial cells and pericyte/smooth muscle cells had the most differentially expressed gene profile compared to other cell types. Top differentially upregulated genes in endothelial cells included novel genes: ROBO4, APCDD1, NDST1, MMRN2, NOTCH4, and DOCK6, as well as previously reported genes: ENG, ORAI2, TFDP1, KDR, AMOTL2, PDGFB, FGFR1, EDN1, and NOTCH1. Several transcription factors were also found to be upregulated in idiopathic pulmonary arterial hypertension endothelial cells including SOX18, STRA13, LYL1, and ELK, which have known roles in regulating endothelial cell phenotype. In particular, SOX18 was implicated through bioinformatics analyses in regulating the idiopathic pulmonary arterial hypertension endothelial cell transcriptome. Furthermore, idiopathic pulmonary arterial hypertension endothelial cells upregulated expression of FAM60A and HDAC7, potentially affecting epigenetic changes in idiopathic pulmonary arterial hypertension endothelial cells. Pericyte/smooth muscle cells expressed genes implicated in regulation of cellular apoptosis and extracellular matrix organization, and several ligands for genes showing increased expression in endothelial cells. In conclusion, our study represents the first detailed look at the transcriptomic landscape across idiopathic pulmonary arterial hypertension lung cells and provides robust insight into alterations that occur in vivo in idiopathic pulmonary arterial hypertension lungs.

Identification of the Human Skeletal Stem Cell

Stem cell regulation and hierarchical organization of human skeletal progenitors remain largely unexplored. Here, we report the isolation of a self-renewing and multipotent human skeletal stem cell (hSSC) that generates progenitors of bone, cartilage, and stroma, but not fat. Self-renewing and multipotent hSSCs are present in fetal and adult bones and can also be derived from BMP2-treated human adipose stroma (B-HAS) and induced pluripotent stem cells (iPSCs). Gene expression analysis of individual hSSCs reveals overall similarity between hSSCs obtained from different sources and partially explains skewed differentiation toward cartilage in fetal and iPSC-derived hSSCs. hSSCs undergo local expansion in response to acute skeletal injury. In addition, hSSC-derived stroma can maintain human hematopoietic stem cells (hHSCs) in serum-free culture conditions. Finally, we combine gene expression and epigenetic data of mouse skeletal stem cells (mSSCs) and hSSCs to identify evolutionarily conserved and divergent pathways driving SSC-mediated skeletogenesis. VIDEO ABSTRACT.

Podocyte development, disease, and stem cell research

The glomerular podocyte is one of the major targets of kidney research. Recent establishment of kidney organoids from pluripotent stem cells has enabled the detailed analysis of human podocytes in both development and disease. The podocytes in organoids express slit diaphragm-related genes and proteins and exhibit characteristic morphology, especially upon experimental transplantation. Organoid technology is now used to reproduce hereditary podocyte diseases, and selective podocyte induction methods have also been reported. Moreover, single-cell RNA-sequencing of human fetal and adult kidneys has revealed the detailed molecular features of this cell lineage, as well as serving as references for kidney organoids in which podocytes are still immature. Here, we discuss the recent progress and limitations of podocyte research from the viewpoint of developmental biology and kidney organoids.

Stromal Microenvironment Shapes the Intratumoral Architecture of Pancreatic Cancer

Single-cell technologies have described heterogeneity across tissues, but the spatial distribution and forces that drive single-cell phenotypes have not been well defined. Combining single-cell RNA and protein analytics in studying the role of stromal cancer-associated fibroblasts (CAFs) in modulating heterogeneity in pancreatic cancer (pancreatic ductal adenocarcinoma [PDAC]) model systems, we have identified significant single-cell population shifts toward invasive epithelial-to-mesenchymal transition (EMT) and proliferative (PRO) phenotypes linked with mitogen-activated protein kinase (MAPK) and signal transducer and activator of transcription 3 (STAT3) signaling. Using high-content digital imaging of RNA in situ hybridization in 195 PDAC tumors, we quantified these EMT and PRO subpopulations in 319,626 individual cancer cells that can be classified within the context of distinct tumor gland "units." Tumor gland typing provided an additional layer of intratumoral heterogeneity that was associated with differences in stromal abundance and clinical outcomes. This demonstrates the impact of the stroma in shaping tumor architecture by altering inherent patterns of tumor glands in human PDAC.