Stromal cell states are altered in active Crohn's disease (CD), but their origin and phenotypic stability are unknown. Using single-cell spatial transcriptomics, RNA sequencing and ATAC sequencing of ~2,500,000 cells, we map 18 distinct stromal cell states within their cellular and cytokine signaling environments in human full-thickness CD bowel. Inflammatory fibroblasts (IFs) reside in immune cell-rich mucosal ulcers and are induced through combinatorial cytokine exposure suppressing submucosal universal fibroblast programs. The IF state is stabilized through transcription factor (TF) activity of GLI3, TWIST1, ETV4, PRDM1 and RELB, and does not spontaneously revert; however, histone deacetylase inhibition destabilizes the IF state, preventing IF secretome-induced epithelial transmigration and activation of neutrophils. The IF open chromatin configuration is distinct from that of fibrotic contractile stroma, which populate adjoining immune-depleted submucosal fibrotic niches. These findings show that microenvironments in pathological tissue niches shape open chromatin configuration of stromal states that are amenable to modulation by epigenetic modifiers.

In the same way that the mRNA-binding specificities of transfer RNAs define the genetic code, the DNA-binding specificities of transcription factors (TFs) form the molecular basis of the gene regulatory code1,2. The human gene regulatory code is much more complex than the genetic code, in particular because there are more than 1,600 TFs that commonly interact with each other. TF-TF interactions are required for specifying cell fate and executing cell-type-specific transcriptional programs. Despite this, the landscape of interactions between DNA-bound TFs is poorly defined. Here we map the biochemical interactions between DNA-bound TFs using CAP-SELEX, a method that can simultaneously identify individual TF binding preferences, TF-TF interactions and the DNA sequences that are bound by the interacting complexes. A screen of more than 58,000 TF-TF pairs identified 2,198 interacting TF pairs, 1,329 of which preferentially bound to their motifs arranged in a distinct spacing and/or orientation. We also discovered 1,131 TF-TF composite motifs that were markedly different from the motifs of the individual TFs. In total, we estimate that the screen identified between 18% and 47% of all human TF-TF motifs. The novel composite motifs we found were enriched in cell-type-specific elements, active in vivo and more likely to be formed between developmentally co-expressed TFs. Furthermore, TFs that define embryonic axes commonly interacted with different TFs and bound to distinct motifs, explaining how TFs with a similar specificity can define distinct cell types along developmental axes.

Combinations of transcription factors (TFs) regulate gene expression and determine cell fate. Much effort has been devoted to understanding TF activity in different tissues and how tissue-specificity is achieved. However, ultimately gene regulation occurs at the single cell level and the recent explosion in the availability of single cell gene expression data now makes it possible to understand TF activity at this granular level of resolution. Here, we leverage a large collection of Human Cell Atlas (HCA) single cell data to explore TF activity by examining cell-type and tissue-specific sets of target genes, or regulons. We compile a regulon atlas, CellRegulon, and map the activity of TFs in an extensive set of healthy adult and foetal tissues spanning hundreds of cell types. Using CellRegulon, we describe dynamic patterns of co-regulation, associate TF-modules with different cellular functions and characterise the distribution of active TFs and TF families across cell types. We show that CellRegulon can link disease gene expression signatures to cell types and TFs relevant to the disease. Finally, using a newly generated multiome dataset of the adult lung, we show how CellRegulon can be extended into an enhancer-gene regulatory network (eGRN) to improve cell-type associations with genetic risk loci for diseases, such as childhood onset asthma, COPD and IPF, and to identify high risk gene modules. Our database for easy download and interactive exploration allows researchers to understand key gene modules activated at cell type transitions and will therefore be valuable for tasks such as cell type engineering (). ### Competing Interest Statement R.E. is a co-founder and holds equity in Ensocell. In the past 3 years, M.C.N. has received remuneration for scientific advisory board membership from GlaxoSmithKline and AstraZeneca paid to the institution, and compensation for travel and accommodation from Roche, GlaxoSmithKline and AstraZeneca paid to the institution. MCN has received unrestricted research grants from GlaxoSmithKline, AstraZeneca, Sanofi, and Roche paid to the institution. In the past 3 years, S.A.T. has received remuneration for scientific advisory board membership from Sanofi, GlaxoSmithKline, Foresite Labs and Qiagen. S.A.T. is a co-founder and holds equity in Transition Bio and Ensocell. From 8 January 2024, S.A.T. has been a part-time employee of GlaxoSmithKline. K.B.M. is an employee at GlaxoSmithKline. The other authors declare no competing interests.

Osteoarthritis is the third most rapidly growing health condition associated with disability, after dementia and diabetes1. By 2050, the total number of patients with osteoarthritis is estimated to reach 1 billion worldwide2. As no disease-modifying treatments exist for osteoarthritis, a better understanding of disease aetiopathology is urgently needed. Here we perform a genome-wide association study meta-analyses across up to 489,975 cases and 1,472,094 controls, establishing 962 independent associations, 513 of which have not been previously reported. Using single-cell multiomics data, we identify signal enrichment in embryonic skeletal development pathways. We integrate orthogonal lines of evidence, including transcriptome, proteome and epigenome profiles of primary joint tissues, and implicate 700 effector genes. Within these, we find rare coding-variant burden associations with effect sizes that are consistently higher than common frequency variant associations. We highlight eight biological processes in which we find convergent involvement of multiple effector genes, including the circadian clock, glial-cell-related processes and pathways with an established role in osteoarthritis (TGFβ, FGF, WNT, BMP and retinoic acid signalling, and extracellular matrix organization). We find that 10% of the effector genes express a protein that is the target of approved drugs, offering repurposing opportunities, which can accelerate translation.

In humans, the bone marrow becomes the primary site for B lymphopoiesis during the second trimester of pregnancy and continues throughout life. Prenatal and adult B cell progenitors play distinct roles in the aetiology and pathology of paediatric and adult hematopoietic malignancies, though the molecular drivers of these differences remain unclear. Here, we created a comprehensive multiomics atlas of over 500k cells covering immune and stromal compartment from prenatal and adult bone marrow, enabling high-resolution analysis of the cell-intrinsic and cell-extrinsic processes that modulate prenatal and adult B lymphopoiesis. Even though B cells follow broadly similar developmental trajectories, we identify a novel postnatal ‘lateCLP’ subset equivalent to prenatal ‘preProB’ cells, and uncover prenatal cells carry several signatures characteristic of leukemias, including enhanced proliferation, higher RAG1/RAG2 activity, extrinsic B cell signals such as IL7 , and lower retention signals in the bone marrow. We also developed and characterised, at both cellular and molecular levels, a new experimental framework for generating B cell precursors from human induced pluripotent stem cells (hiPSCs), and show that it faithfully recapitulates key stages of B cell differentiation. Together, our single-cell multiomics atlas of B lymphopoiesis in vivo and in vitro offers detailed insights into the unique molecular features of prenatal and adult B cell lymphopoiesis, and serves as a powerful resource for investigating the early events that contribute to haematological disorders. ### Competing Interest Statement In the past 3 years, S.A.T. has received remuneration for scientific advisory board membership from Sanofi, GlaxoSmithKline, Foresite Labs and Qiagen. S.A.T. is a co-founder and holds equity in Transition Bio and Ensocell. From 8 January 2024, S.A.T. has been a part-time employee of GlaxoSmithKline. Y.S.M. and P.W.Z. are named inventors on patents for T cell differentiation technology. P.W.Z. is a co-founder of Notch Therapeutics. P.W.Z. and Y.S.M. are co-founders of Apiary Therapeutics. P.W.Z. consult for cell therapy companies. This work is currently under consideration for patent application. Wellcome Trust, 220540/Z/20/A Wellcome Leap HOPE Program EMBO long-term, ALTF 737-2021 UKRI Postdoctoral Fellowship awarded based on a selected proposal under the European Union’s Horizon 2021 research and innovation programme, within the Marie Skłodowska-Curie Actions Spanish Health Institute Carlos III, CP22/00127 Deutsche Forschungsgemeinschaft, 468499998

Human embryonic bone and joint formation is determined by coordinated differentiation of progenitors in the nascent skeleton. The cell states, epigenetic processes and key regulatory factors that underlie lineage commitment of these cells remain elusive. Here we applied paired transcriptional and epigenetic profiling of approximately 336,000 nucleus droplets and spatial transcriptomics to establish a multi-omic atlas of human embryonic joint and cranium development between 5 and 11 weeks after conception. Using combined modelling of transcriptional and epigenetic data, we characterized regionally distinct limb and cranial osteoprogenitor trajectories across the embryonic skeleton and further described regulatory networks that govern intramembranous and endochondral ossification. Spatial localization of cell clusters in our in situ sequencing data using a new tool, ISS-Patcher, revealed mechanisms of progenitor zonation during bone and joint formation. Through trajectory analysis, we predicted potential non-canonical cellular origins for human chondrocytes from Schwann cells. We also introduce SNP2Cell, a tool to link cell-type-specific regulatory networks to polygenic traits such as osteoarthritis. Using osteolineage trajectories characterized here, we simulated in silico perturbations of genes that cause monogenic craniosynostosis and implicate potential cell states and disease mechanisms. This work forms a detailed and dynamic regulatory atlas of bone and cartilage maturation and advances our fundamental understanding of cell-fate determination in human skeletal development.

Skeletal muscle aging is a key contributor to age-related frailty and sarcopenia with substantial implications for global health. Here we profiled 90,902 single cells and 92,259 single nuclei from 17 donors to map the aging process in the adult human intercostal muscle, identifying cellular changes in each muscle compartment. We found that distinct subsets of muscle stem cells exhibit decreased ribosome biogenesis genes and increased CCL2 expression, causing different aging phenotypes. Our atlas also highlights an expansion of nuclei associated with the neuromuscular junction, which may reflect re-innervation, and outlines how the loss of fast-twitch myofibers is mitigated through regeneration and upregulation of fast-type markers in slow-twitch myofibers with age. Furthermore, we document the function of aging muscle microenvironment in immune cell attraction. Overall, we present a comprehensive human skeletal muscle aging resource ( https://www.muscleageingcellatlas.org/ ) together with an in-house mouse muscle atlas to study common features of muscle aging across species.

Bone and joint formation in the developing skeleton rely on co-ordinated differentiation of progenitors in the nascent developing limbs and joints. The cell states, epigenetic processes and key regulatory factors underlying their lineage commitment to osteogenic and other mesenchymal populations during ossification and joint formation remain poorly understood and are largely unexplored in human studies. Here, we apply paired single-nuclei transcriptional and epigenetic profiling of 336,000 droplets, in addition to spatial transcriptomics, to construct a comprehensive atlas of human bone, cartilage and joint development in the shoulder, hip, knee and cranium from 5 to 11 post-conception weeks. Spatial mapping of cell clusters to our highly multiplexed in situ sequencing (ISS) data using our newly developed tool ISS-Patcher revealed new cellular mechanisms of zonation during bone and joint formation. Combined modelling of chromatin accessibility and RNA expression allowed the identification of the transcriptional and epigenetic regulatory landscapes that drive differentiation of mesenchymal lineages including osteogenic and chondrogenic lineages, and novel chondrocyte cell states. In particular, we define regionally distinct limb and cranial osteoprogenitor populations and trajectories across the fetal skeleton and characterise differential regulatory networks that govern intramembranous and endochondral ossification. Through somatic mutation analysis, we predict two new potential cell origins for human chondrocyte development. We also introduce SNP2Cell, a tool to link cell-type specific regulatory networks to numerous polygenic traits such as osteoarthritis. We also conduct in silico perturbations of genes that cause monogenic craniosynostosis and implicate potential pathogenic cell states and disease mechanisms involved. This work forms a detailed and dynamic regulatory atlas of human fetal skeletal maturation and advances our fundamental understanding of cell fate determination in human skeletal development. ### Competing Interest Statement C.D.B is a founder of Mestag Therapeutics. In the past 3 years, S.A.T. has received remuneration for scientific advisory board membership from Sanofi, GlaxoSmithKline, Foresite Labs and Qiagen. S.A.T. is a co-founder and holds equity in Transition Bio and Ensocell. From 8 January 2024, S.A.T. is a part-time employee of GlaxoSmithKline. The remaining authors declare no competing interests.

Studies of human lung development have focused on epithelial and mesenchymal cell types and function, but much less is known about the developing lung immune cells, even though the airways are a major site of mucosal immunity after birth. An unanswered question is whether tissue-resident immune cells play a role in shaping the tissue as it develops in utero. Here, we profiled human embryonic and fetal lung immune cells using scRNA-seq, smFISH, and immunohistochemistry. At the embryonic stage, we observed an early wave of innate immune cells, including innate lymphoid cells, natural killer cells, myeloid cells, and lineage progenitors. By the canalicular stage, we detected naive T lymphocytes expressing high levels of cytotoxicity genes and the presence of mature B lymphocytes, including B-1 cells. Our analysis suggests that fetal lungs provide a niche for full B cell maturation. Given the presence and diversity of immune cells during development, we also investigated their possible effect on epithelial maturation. We found that IL-1β drives epithelial progenitor exit from self-renewal and differentiation to basal cells in vitro. In vivo, IL-1β–producing myeloid cells were found throughout the lung and adjacent to epithelial tips, suggesting that immune cells may direct human lung epithelial development.

Developmental dynamics encompass both the specification of cell types and their spatial organisation into multicellular niches. Here we harness the power of single-cell and spatial multiomics to unravel embryonic and foetal cardiac tissue niches, which lead to the development of a new tool, TissueTypist. We infer developmental cell trajectories, including evidence for lineage relationships based on shared somatic mutations, within first- and second-trimester human hearts. We reveal that cardiac-resident macrophages likely originate from the yolk sac, forming heterogeneous subsets. CX3CR1+ macrophages with a microglia-like profile localise in the sinoatrial node, which may contribute to axon guidance for the innervating autonomic neurons. Foetal pacemaker cells exhibit distinct characteristics compared to their adult counterparts, including the expression of genes that are known to promote parasympathetic innervation. By comparing somatic mutation profiles of cardiomyocytes, we identify an early branching point where pacemaker cells diverge from working cardiomyocytes. We highlight the enhancer-mediated gene regulatory networks governing atrial and ventricular cardiomyocyte specification. The maturation of atrial cardiomyocytes into distinct left and right phenotypes, driven by transcription factors linked to atrial septal defect genes, underscores the significance of this process for healthy heart development. In the ventricle, cellular and transcriptional gradients along both pseudotime and the transmural axis provide a new molecular understanding of myocardial compaction. Finally, generating data from Trisomy 21 hearts and comparing this with the euploid atlas, we reveal a reduced abundance of specific cell types including compact cardiomyocytes. Overall, this extensive dataset and our precomputed models will form a valuable resource for the field. ### Competing Interest Statement S.A.T. is a scientific advisory board member of ForeSite Labs, Qiagen and Element Biosciences, and a co-founder and equity holder of TransitionBio and EnsoCell Therapeutics, and a part-time employee of GlaxoSmithKline. S.S. is a co-founder and equity holder of ABS Biotechnologies. The remaining authors declare no competing interests.

Childhood-onset asthma is characterized by Type 2-inflammation and airway wall remodeling, but mechanisms of asthma development in the first years of life remain unclear. Here, we investigate transcriptional changes in airway wall biopsies of 22 symptomatic one year old children and relate these to asthma at school age. We demonstrate that pre-asthmatic children (n = 10) overexpressed a gene signature characteristic for an airway epithelial differentiation trajectory via hillock cells towards squamous cells (adjusted p-value 8.06e-16), whilst there was no association with gene signatures of Type 2-inflammation or eosinophil activation. Genes expressed along this trajectory are linked to an altered epithelial barrier function, innate immune activation and extracellular matrix remodeling. Functional GWAS analysis supports a causal link between childhood-onset, but not adult-onset asthma, and the hillock-squamous cell differentiation trajectory. Next, we confirmed the presence of hillock-like cells at the RNA and protein level in pediatric upper and lower airway samples. These findings identify a novel mechanism by which an aberrant airway epithelial differentiation trajectory may contribute to a pre-asthmatic state, highlighting the difference between the early origins of childhood-onset asthma and adult asthma, and point to possible new targets for the early diagnosis and treatment of asthma in the first two years of life. One Sentence Summary RNA sequencing in bronchial biopsies from wheezing infants and children < 2 years shows evidence for an airway epithelial hillock-to-squamous differentiation pathway that marks the development of asthma.

Abstract Introduction Cardiogenesis requires the differentiation and coordinated function of a large variety of cell types. Congenital heart disease (CHD), affecting 1% of live births[1], and several adult-onset heart diseases result when these processes go awry[2–4]. To understand the cellular and molecular mechanisms at play we must define the repertoire of cell types in the developing heart as well as the cellular niches within which they interact. Single-cell RNA sequencing is a powerful means of defining cell populations, however existing studies of the human foetal heart[5,6] have been hampered by low cell numbers which limits detection of rare but important cell types. To create an improved atlas we increased the number of cells, increased the age range of foetal hearts sampled, and generated spatial transcriptomic data to map the defined cells to multicellular niches of co-located and interacting cells. To better infer the transcription factors (TFs) governing cell identity using gene regulatory network (GRN) inference we also generated paired single-cell resolution RNA and chromatin accessibility profiles using Assay for Transposase-Accessible Chromatin (ATAC) sequencing. Methods We performed paired single-nuclei RNA and ATAC sequencing (10X Multiome) and single-cell RNAseq (10X 3’ RNAseq) on 23 human foetal hearts ranging from 4 to 20 post-conception weeks. In addition we generated spatial transcriptomic data (10X Visium) from these hearts. We validated cell types defined in single-cell data using protein-level immunofluorescence imaging. Results Our healthy reference atlas includes 300k cells and nuclei. Using their gene expression profile we resolved 46 fine-grained cell types (FIGURE 1). We used cell2location[7], a machine learning classifier trained on the atlas to map the defined cell types to the spatial data (FIGURE 2) and detected niches of co-locating cells, such as a macrophage population in the walls of the great vessels predicted to interact with mural cells. We were able to define location-specific transcriptional signatures such as cardiomyocytes of each chamber (as well as the compact and trabeculated myocardium) and smooth muscle cells specific to great vessels vs coronary arteries. Using a combination of RNA and ATAC data we inferred cell-specific GRNs, which implicated novel TFs such as MSX2 in sinoatrial node pacemaker cells and the enrichment of CHD genes amongst cardiomyocyte-specific GRNs. In the spatial data we showed that gene modules associated with different forms of CHD enrich in specific anatomical regions, such as the outflow tract. Discussion: Together these data constitute the most comprehensive cardiac developmental cell atlas to date. It sheds new light on the diversity of cell states, their niches, and their vulnerability to the genetic causes of CHD.Figure 1:UMAP showing atlasFigure 2:Spatial transcriptomic mapping

Scientists employing omics in life science studies face challenges such as the modeling of multiassay studies, recording of all relevant parameters, and managing many samples with their metadata. They must manage many large files that are the results of the assays or subsequent computation. Users with diverse backgrounds, ranging from computational scientists to wet-lab scientists, have dissimilar needs when it comes to data access, with programmatic interfaces being favored by the former and graphical ones by the latter. We introduce SODAR, the system for omics data access and retrieval. SODAR is a software package that addresses these challenges by providing a web-based graphical user interface for managing multiassay studies and describing them using the ISA (Investigation, Study, Assay) data model and the ISA-Tab file format. Data storage is handled using the iRODS data management system, which handles large quantities of files and substantial amounts of data. SODAR also offers programmable APIs and command-line access for metadata and file storage. SODAR supports complex omics integration studies and can be easily installed. The software is written in Python 3 and freely available at https://github.com/bihealth/sodar-server under the MIT license.

The function of a cell is defined by its intrinsic characteristics and its niche: the tissue microenvironment in which it dwells. Here we combine single-cell and spatial transcriptomics data to discover cellular niches within eight regions of the human heart. We map cells to microanatomical locations and integrate knowledge-based and unsupervised structural annotations. We also profile the cells of the human cardiac conduction system1. The results revealed their distinctive repertoire of ion channels, G-protein-coupled receptors (GPCRs) and regulatory networks, and implicated FOXP2 in the pacemaker phenotype. We show that the sinoatrial node is compartmentalized, with a core of pacemaker cells, fibroblasts and glial cells supporting glutamatergic signalling. Using a custom CellPhoneDB.org module, we identify trans-synaptic pacemaker cell interactions with glia. We introduce a druggable target prediction tool, drug2cell, which leverages single-cell profiles and drug-target interactions to provide mechanistic insights into the chronotropic effects of drugs, including GLP-1 analogues. In the epicardium, we show enrichment of both IgG+ and IgA+ plasma cells forming immune niches that may contribute to infection defence. Overall, we provide new clarity to cardiac electro-anatomy and immunology, and our suite of computational approaches can be applied to other tissues and organs.

Summary Cdr1as is a conserved circular RNA (circRNA) enriched in the CNS and important for maintaining brain homeostasis. The loss of Cdr1as results in aberrant synaptic transmission and deregulation of stress response and circadian clock genes. However, it is not known whether the expression of Cdr1as or circRNAs, in general, follows a circadian pattern in different tissues. Here, using newly generated and public RNA-Seq data, we monitor circRNA expression throughout circadian rhythm in various mouse brain regions. We demonstrate that Cdr1as , despite its stable character, has a highly dynamic expression during the circadian cycle in the mouse suprachiasmatic nucleus (SCN). Cdr1as is one of the highest expressed RNAs in a cluster associated with light-induced synaptic transmission and phase shift in the SCN. Further, we identified that another brain enriched circRNA, mbl , is also substantially deregulated upon light induction in the fly head. Our study highlights the potential impact of abundant and conserved circRNAs on maintaining a healthy circadian cycle across species.

We present a multiomic cell atlas of human lung development that combines single-cell RNA and ATAC sequencing, high-throughput spatial transcriptomics, and single-cell imaging. Coupling single-cell methods with spatial analysis has allowed a comprehensive cellular survey of the epithelial, mesenchymal, endothelial, and erythrocyte/leukocyte compartments from 5-22 post-conception weeks. We identify previously unchar-acterized cell states in all compartments. These include developmental-specific secretory progenitors and a subtype of neuroendocrine cell related to human small cell lung cancer. Our datasets are available through our web interface (https://lungcellatlas.org). To illustrate its general utility, we use our cell atlas to generate predictions about cell-cell signaling and transcription factor hierarchies which we rigorously test using orga-noid models.

Single-cell transcriptomics has allowed unprecedented resolution of cell types/states in the human lung, but their spatial context is less well defined. To (re)define tissue architecture of lung and airways, we profiled five proximal-to-distal locations of healthy human lungs in depth using multi-omic single cell/nuclei and spatial transcriptomics (queryable at lungcellatlas.org ). Using computational data integration and analysis, we extend beyond the suspension cell paradigm and discover macro and micro-anatomical tissue compartments including previously unannotated cell types in the epithelial, vascular, stromal and nerve bundle micro-environments. We identify and implicate peribronchial fibroblasts in lung disease. Importantly, we discover and validate a survival niche for IgA plasma cells in the airway submucosal glands (SMG). We show that gland epithelial cells recruit B cells and IgA plasma cells, and promote longevity and antibody secretion locally through expression of CCL28, APRIL and IL-6. This new 'gland-associated immune niche' has implications for respiratory health.

Circular RNAs (circRNAs) are a large class of relatively stable RNA molecules that are highly expressed in animal brains. Many circRNAs have been associated with CNS disorders accompanied by an aberrant wake-sleep cycle. However, the regulation of circRNAs in brain homeostasis over daily light-dark (LD) cycles has not been characterized. Here, we aim to quantify the daily expression changes of circRNAs in physiological conditions in healthy adult animals. Using newly generated and public RNA-Seq data, we monitored circRNA expression throughout the 12:12 h LD cycle in various mouse brain regions. We identified that Cdr1as, a conserved circRNA that regulates synaptic transmission, is highly expressed in the suprachiasmatic nucleus (SCN), the master circadian pacemaker. Despite its high stability, Cdr1as has a very dynamic expression in the SCN throughout the LD cycle, as well as a significant regulation in the hippocampus following the entry into the dark phase. Computational integration of different public datasets predicted that Cdr1as is important for regulating light entrainment in the SCN. We hypothesize that the expression changes of Cdr1as in the SCN, particularly during the dark phase, are associated with light-induced phase shifts. Importantly, our work revises the current beliefs about natural circRNA stability and suggests that the time component must be considered when studying circRNA regulation.

Lung carcinoid tumors, also referred to as pulmonary neuroendocrine tumors or lung carcinoids, are rare neoplasms of the lung with a more favorable prognosis than other subtypes of lung cancer. Still, some patients suffer from relapsed disease and metastatic spread while no consensus treatment exists for metastasized carcinoids. Several recent single-cell studies have provided detailed insights into the cellular heterogeneity of more common lung cancers, such as adeno- and squamous cell carcinoma. However, the characteristics of lung carcinoids on the single-cell level are yet completely unknown. To study the cellular composition and single-cell gene expression profiles in lung carcinoids, we applied single-cell RNA sequencing to three lung carcinoid tumor samples and normal lung tissue. The single-cell transcriptomes of carcinoid tumor cells reflected intertumoral heterogeneity associated with clinicopathological features, such as tumor necrosis and proliferation index. The immune microenvironment was specifically enriched in noninflammatory monocyte-derived myeloid cells. Tumor-associated endothelial cells were characterized by distinct gene expression profiles. A spectrum of vascular smooth muscle cells and pericytes predominated the stromal microenvironment. We found a small proportion of myofibroblasts exhibiting features reminiscent of cancer-associated fibroblasts. Stromal and immune cells exhibited potential paracrine interactions which may shape the microenvironment via NOTCH, VEGF, TGFβ and JAK/STAT signaling. Moreover, single-cell gene signatures of pericytes and myofibroblasts demonstrated prognostic value in bulk gene expression data. Here, we provide first comprehensive insights into the cellular composition and single-cell gene expression profiles in lung carcinoids, demonstrating the non-inflammatory and vessel-rich nature of their tumor microenvironment, and outlining relevant intercellular interactions which could serve as future therapeutic targets.