Abstract The thyroid has a remarkable evolution, transforming from an exocrine constituent of the chordate endostyle to an endocrine gland in basal vertebrates. In mammals, a second endocrine cell type of presumed neural crest origin appears, although recent lineage tracing has firmly established thyroid C-cells are also endoderm-derived. Here, we characterize the global gene expression in both embryonic thyroid lineages at single-cell level and identify lineage-specific transcription factors and their network regulation of target genes implicated in thyroid development into a dual endocrine organ. Merging of the pharyngeal pouch-derived ultimobranchial bodies with the midline thyroid primordium is an ordered process featuring basement membrane dynamics and epithelial-mesenchymal plasticity required for precursor cells to disseminate and properly integrate, thus forming the typical histoarchitecture of thyroid follicles and parafollicular C-cells. Synchronous lineage growth of compound follicles is recapitulated in mixed-type thyroid carcinoma in which only the neuroendocrine tumor cells escape the follicle basement membrane boundaries and become invasive adopting a C-cell precursor-like migratory phenotype.

Resolving the three-dimensional structure of chromatin with haplotype specificity remains a fundamental challenge for understanding genetic and epigenetic contributions to gene regulation in healthy development and in disease. Phasing of chromatin contacts derived by ligation-dependent methods relies on high single-nucleotide variant (SNV) density, which limits its applicability in human genomes where SNVs are sparse. Here, we present CoPhasing, a strategy that leverages the intrinsic haplotype fidelity of Genome Architecture Mapping (GAM), a ligation-free method that inherently captures haplotype-matched genomic neighborhoods in thin nuclear slices. CoPhasing assigns SNV-free reads to their haplotype through proximity to informative reads, enabling efficient and accurate phasing even in SNV-sparse regions. Applying CoPhasing genome-wide reveals extensive, previously inaccessible structural heterogeneity between human homologous chromosomes, demonstrating the power of GAM CoPhasing to investigate allele-specific differences in 3D genome organization and their relevance to gene regulation and disease.

The thymus is the primary site of T cell development, central to the establishment of self-tolerance and adaptive immune function. In mammals, the thymus undergoes age-related involution, resulting in a global decline in immune function. The thymus has some regenerative ability that relies on pre-existing thymic remnants but is insufficient to prevent involution. Here, we show that the juvenile axolotl (Ambystoma mexicanum) is able to regenerate its thymus de novo after complete removal, constituting an exception among vertebrates. Using single-cell transcriptomics and genetic and transplantation approaches, we demonstrate that de novo thymus regeneration results in the restoration of morphology, cell-type diversity, and function. FOXN1, although it has a conserved role in thymus organogenesis, is dispensable for the initiation of thymic regeneration. In contrast, we identify midkine signaling as a possible early driver of de novo thymus regeneration. This study demonstrates an instance of organ-level regeneration of the lymphoid system, which could guide future clinical strategies seeking to promote thymus regrowth.

The Long-read RNA-Seq Genome Annotation Assessment Project (LRGASP) Consortium was formed to evaluate the effectiveness of long-read approaches for transcriptome analysis. The consortium generated over 427 million long-read sequences from cDNA and direct RNA datasets, encompassing human, mouse, and manatee species, using different protocols and sequencing platforms. These data were utilized by developers to address challenges in transcript isoform detection and quantification, as well as de novo transcript isoform identification. The study revealed that libraries with longer, more accurate sequences produce more accurate transcripts than those with increased read depth, whereas greater read depth improved quantification accuracy. In well-annotated genomes, tools based on reference sequences demonstrated the best performance. When aiming to detect rare and novel transcripts or when using reference-free approaches, incorporating additional orthogonal data and replicate samples are advised. This collaborative study offers a benchmark for current practices and provides direction for future method development in transcriptome analysis.

The thyroid has a remarkable evolution, appearing first as an integral presumably exocrine constituent of the chordate endostyle that is transformed into an endocrine gland during metamorphosis in basal vertebrates. In mammals, the thyroid acquires a second endocrine cell type, calcitonin-producing C-cells, which for long were inferred a neural crest origin, shuttled to the embryonic thyroid by the ultimobranchial bodies, although recent lineage tracing experiments firmly establish these neuroendocrine cells also derive from foregut endoderm. A key question remains unanswered, how thyroid primordia independently develop and, unlike in all non-mammalian vertebrates in which thyroid hormone and calcitonin are secreted from anatomically and functionally unconnected glands, merge into a dual endocrine organ. Here we characterize, leveraging of a single-cell transcriptome atlas derived from mouse pharyngeal endoderm and its subsequent cell fates, the global gene expression profile of thyroid- and ultimobranchial-derived progenitor cells, and identify comprehensive gene regulatory networks of lineage-specific transcription factors and novel target genes predicted to differentially regulate cell proliferation, plasticity and differentiation during development. Spatiotemporal analyses reveal C-cell precursors are triggered to undergo epithelial-mesenchymal transition (EMT) and cell-autonomously down-regulate collagen IV and degrade laminin that encloses the ultimobranchial body. However, the EMT program is not fully deployed until both cell lineages are mixed and propagate conjointly thus forming the typical thyroid histoarchitecture of stationary follicular epithelial cells and parafollicular C-cells, every follicle/C-cell unit being enveloped by a renewed basement membrane. Mixed-type thyroid carcinoma recapitulates a synchronous lineage growth pattern but only the neuroendocrine tumor cells are able escaping the compound follicle boundaries and become invasive adopting a C-cell precursor migratory phenotype.

A detailed understanding of the innate immune mechanisms involved in restricting SARS-CoV-2 infection and how the virus disrupts these processes could reveal new strategies to boost antiviral mechanisms and develop therapeutics for COVID-19. Here, we identify cellular nucleic acid-binding protein (CNBP) as a key host factor controlling SARS-CoV-2 infection. In response to RNA-sensing pathways, CNBP is phosphorylated and translocates from the cytosol to the nucleus where it binds to the interferon-β enhancer to initiate transcription. Because SARS-CoV-2 evades immune detection by the host’s RNA-sensing pathways, CNBP is largely retained in the cytosol where it restricts SARS-CoV-2 directly, leading to a battle between the host and SARS-CoV-2 that extends beyond antiviral immune signaling pathways. We further demonstrated that CNBP binds SARS-CoV-2 viral RNA directly and competes with the viral nucleocapsid protein to prevent viral RNA and nucleocapsid protein from forming liquid–liquid phase separation (LLPS) condensates critical for viral replication. Consequently, cells and animals lacking CNBP have higher viral loads, and CNBP-deficient mice succumb rapidly to infection. Altogether, these findings identify CNBP as a key antiviral factor for SARS-CoV-2, functioning both as a regulator of antiviral IFN gene expression and a cell-intrinsic restriction factor that disrupts LLPS to limit viral replication and spread. In addition, our studies also highlight viral condensates as important targets and strategies for the development of drugs to combat COVID-19.

The continuous evolution of SARS-CoV-2 variants complicates efforts to combat the ongoing pandemic, underscoring the need for a dynamic platform for the rapid development of pan-viral variant therapeutics. Oligonucleotide therapeutics are enhancing the treatment of numerous diseases with unprecedented potency, duration of effect, and safety. Through the systematic screening of hundreds of oligonucleotide sequences, we identified fully chemically stabilized siRNAs and ASOs that target regions of the SARS-CoV-2 genome conserved in all variants of concern, including delta and omicron. We successively evaluated candidates in cellular reporter assays, followed by viral inhibition in cell culture, with eventual testing of leads for in vivo antiviral activity in the lung. Previous attempts to deliver therapeutic oligonucleotides to the lung have met with only modest success. Here, we report the development of a platform for identifying and generating potent, chemically modified multimeric siRNAs bioavailable in the lung after local intranasal and intratracheal delivery. The optimized divalent siRNAs showed robust antiviral activity in human cells and mouse models of SARS-CoV-2 infection and represent a new paradigm for antiviral therapeutic development for current and future pandemics.

Tissue-resident stem and progenitor cells are present in many adult organs, where they are important for organ homeostasis and repair in response to injury. However, the signals that activate these cells and the mechanisms governing how these cells renew or differentiate are highly context-dependent and incompletely understood, particularly in non-hematopoietic tissues. In the skin, melanocyte stem and progenitor cells are responsible for replenishing mature pigmented melanocytes. In mammals, these cells reside in the hair follicle bulge and bulb niches where they are activated during homeostatic hair follicle turnover and following melanocyte destruction, as occurs in vitiligo and other skin hypopigmentation disorders. Recently, we identified melanocyte progenitors in adult zebrafish skin. To elucidate mechanisms governing melanocyte progenitor renewal and differentiation we analyzed individual transcriptomes from thousands of melanocyte lineage cells during the regeneration process. We identified transcriptional signatures for progenitors, deciphered transcriptional changes and intermediate cell states during regeneration, and analyzed cell–cell signaling changes to discover mechanisms governing melanocyte regeneration. We identified KIT signaling via the RAS/MAPK pathway as a regulator of melanocyte progenitor direct differentiation and asymmetric division. Our findings show how activation of different subpopulations of mitfa -positive cells underlies cellular transitions required to properly reconstitute the melanocyte pigmentary system following injury.

Approaches to study human pharyngeal foregut endoderm-a developmental intermediate that is linked to various human syndromes involving pharynx development and organogenesis of tissues such as thymus, parathyroid, and thyroid-have been hampered by scarcity of tissue access and cellular models. We present an efficient stepwise differentiation method to generate human pharyngeal foregut endoderm from pluripotent stem cells. We determine dose and temporal requirements of signaling pathway engagement for optimized differentiation and characterize the differentiation products on cellular and integrated molecular level. We present a computational classification tool, "CellMatch," and transcriptomic classification of differentiation products on an integrated mouse scRNA-seq developmental roadmap confirms cellular maturation. Integrated transcriptomic and chromatin analyses infer differentiation stage-specific gene regulatory networks. Our work provides the method and integrated multiomic resource for the investigation of disease-relevant loci and gene regulatory networks and their role in developmental defects affecting the pharyngeal endoderm and its derivatives.

Long non-coding RNAs (lncRNAs) have emerged as fundamental regulators in various biological processes, including embryonic development and cellular differentiation. Despite much progress over the past decade, the genome-wide annotation of lncRNAs remains incomplete and many known non-coding loci are still poorly characterized. Here, we report the discovery of a previously unannotated lncRNA that is transcribed 230 kb upstream of the SOX17 gene and located within the same topologically associating domain. We termed it T-REX17 ( T ranscript R egulating E ndoderm and activated by so X17 ) and show that it is induced following SOX17 activation but its expression is more tightly restricted to early definitive endoderm. Loss of T-REX17 affects crucial functions independent of SOX17 and leads to an aberrant endodermal transcriptome, signaling pathway deregulation and epithelial to mesenchymal transition defects. Consequently, cells lacking the lncRNA cannot further differentiate into more mature endodermal cell types. Taken together, our study identified and characterized T-REX17 as a transiently expressed and essential non-coding regulator in early human endoderm differentiation.

Maldevelopment of the pharyngeal endoderm, an embryonic tissue critical for patterning of the pharyngeal region and ensuing organogenesis, ultimately contributes to several classes of human developmental syndromes and disorders. Such syndromes are characterized by a spectrum of phenotypes that currently cannot be fully explained by known mutations or genetic variants due to gaps in characterization of critical drivers of normal and dysfunctional development. Despite the disease-relevance of pharyngeal endoderm, we still lack a comprehensive and integrative view of the molecular basis and gene regulatory networks driving pharyngeal endoderm development. To close this gap, we apply transcriptomic and chromatin accessibility single-cell sequencing technologies to generate a multi-omic developmental resource spanning pharyngeal endoderm patterning to the emergence of organ-specific epithelia in the developing mouse embryo. We identify cell-type specific gene regulation, distill GRN models that define developing organ domains, and characterize the role of an immunodeficiency-associated forkhead box transcription factor.

Abstract The lung is a complex organ with various cell types having distinct roles. Antisense oligonucleotides (ASOs) have been studied in the lung, but it has been challenging to determine their effectiveness in each cell type due to the lack of appropriate analytical methods. We employed three distinct approaches to study silencing efficacy within different cell types. First, we used lineage markers to identify cell types in flow cytometry, and simultaneously measured ASO-induced silencing of cell-surface proteins CD47 or CD98. Second, we applied single-cell RNA sequencing (scRNA-seq) to measure silencing efficacy in distinct cell types; to the best of our knowledge, this is the first time scRNA-seq has been applied to measure the efficacy of oligonucleotide therapeutics. In both approaches, fibroblasts were the most susceptible to locally delivered ASOs, with significant silencing also in endothelial cells. Third, we confirmed that the robust silencing in fibroblasts is broadly applicable by silencing two targets expressed mainly in fibroblasts, Mfap4 and Adam33. Across independent approaches, we demonstrate that intratracheally administered LNA gapmer ASOs robustly induce gene silencing in lung fibroblasts. ASO-induced gene silencing in fibroblasts was durable, lasting 4–8 weeks after a single dose. Thus, lung fibroblasts are well aligned with ASOs as therapeutics.

ABSTRACT Long non-coding RNAs (lncRNAs) have emerged as fundamental regulators in various biological processes, including embryonic development and cellular differentiation. Despite much progress over the past decade, the genome-wide annotation of lncRNAs remains incomplete and many known non-coding loci are still poorly characterized. Here, we report the discovery of a previously not annotated lncRNA that is transcribed upstream of the SOX17 gene and located within the same topologically associating domain. We termed it LNCSOX17 and show that it is induced following SOX17 activation but its expression is more tightly restricted to early definitive endoderm. Loss of LNCSOX17 affects crucial functions independent of SOX17 and leads to an aberrant endodermal transcriptome, signaling pathway deregulation and epithelial to mesenchymal transition defects. Consequently, cells lacking the lncRNA cannot further differentiate into more mature endodermal cell types. We identified and characterized LNCSOX17 as an essential new actor in early human endoderm, thereby further expanding the list of functionally important non-coding regulators.

Novel pathogenic severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) continues to pose an imminent global threat since its initial outbreak in December 2019. A simple in vitro model system using cell lines highly susceptible to SARS-CoV-2 infection are critical to facilitate the study of the virus cycle and to discover effective antivirals against the virus. Human lung alveolar A549 cells are regarded as a useful and valuable model for respiratory virus infection. However, SARS-CoV-2 uses the ACE2 as receptor for viral entry and the TMPRSS2 to prime the Spike protein, both of which are negligibly expressed in A549 cells. Here, we report the generation of a robust human lung epithelial cell-based model by transducing ACE2 and TMPRSS2 into A549 cells and show that the ACE2 enriched A549ACE2/TMPRSS2 cells (ACE2plus) and its single-cell-derived subclone (ACE2plusC3) are highly susceptible to SARS-CoV-2 infection. These engineered ACE2plus showed higher ACE2 and TMPRSS2 mRNA expression levels than currently used Calu3 and commercial A549ACE2/TMPRSS2 cells. ACE2 and TMPRSS2 proteins were also highly and ubiquitously expressed in ACE2plusC3 cells. Additionally, antiviral drugs like Camostat mesylate, EIDD-1931, and Remdesivir strongly inhibited SARS-CoV-2 replication. Notably, multinucleated syncytia, a clinical feature commonly observed in severe COVID-19 patients was induced in ACE2plusC3 cells either by virus infection or by overexpressing the Spike proteins of different variants of SARS-CoV-2. Syncytial process was effectively blocked by the furin protease inhibitor, Decanoyl-RVKR-CMK. Taken together, we have developed a robust human A549 lung epithelial cell-based model that can be applied to probe SARS-CoV-2 replication and to facilitate the discovery of SARS-CoV-2 inhibitors.

Skeletal muscle myoblasts (iMyoblasts) were generated from human induced pluripotent stem cells (iPSCs) using an efficient and reliable transgene-free induction and stem cell selection protocol. Immunofluorescence, flow cytometry, qPCR, digital RNA expression profiling, and scRNA-Seq studies identify iMyoblasts as a PAX3+/MYOD1+ skeletal myogenic lineage with a fetal-like transcriptome signature, distinct from adult muscle biopsy myoblasts (bMyoblasts) and iPSC-induced muscle progenitors. iMyoblasts can be stably propagated for >12 passages or 30 population doublings while retaining their dual commitment for myotube differentiation and regeneration of reserve cells. iMyoblasts also efficiently xenoengrafted into irradiated and injured mouse muscle where they undergo differentiation and fetal-adult MYH isoform switching, demonstrating their regulatory plasticity for adult muscle maturation in response to signals in the host muscle. Xenograft muscle retains PAX3+ muscle progenitors and can regenerate human muscle in response to secondary injury. As models of disease, iMyoblasts from individuals with Facioscapulohumeral Muscular Dystrophy revealed a previously unknown epigenetic regulatory mechanism controlling developmental expression of the pathological DUX4 gene. iMyoblasts from Limb-Girdle Muscular Dystrophy R7 and R9 and Walker Warburg Syndrome patients modeled their molecular disease pathologies and were responsive to small molecule and gene editing therapeutics. These findings establish the utility of iMyoblasts for ex vivo and in vivo investigations of human myogenesis and disease pathogenesis and for the development of muscle stem cell therapeutics.

Chromosome conformation capture (3C) assays are used to map chromatin interactions genome-wide. Chromatin interaction maps provide insights into the spatial organization of chromosomes and the mechanisms by which they fold. Hi-C and Micro-C are widely used 3C protocols that differ in key experimental parameters including cross-linking chemistry and chromatin fragmentation strategy. To understand how the choice of experimental protocol determines the ability to detect and quantify aspects of chromosome folding we have performed a systematic evaluation of 3C experimental parameters. We identified optimal protocol variants for either loop or compartment detection, optimizing fragment size and cross-linking chemistry. We used this knowledge to develop a greatly improved Hi-C protocol (Hi-C 3.0) that can detect both loops and compartments relatively effectively. In addition to providing benchmarked protocols, this work produced ultra-deep chromatin interaction maps using Micro-C, conventional Hi-C and Hi-C 3.0 for key cell lines used by the 4D Nucleome project.

Coronaviruses are a family of RNA viruses that cause acute and chronic diseases of the upper and lower respiratory tract in humans and other animals. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a recently emerged coronavirus that has led to a global pandemic, causing a severe respiratory disease known as coronavirus disease 2019 (COVID-19) with substantial morbidity and mortality worldwide. The development of antiviral therapeutics is urgently needed while vaccine programs roll out worldwide. Here, we describe a diamidobenzimidazole compound, diABZI-4, that activates stimulator of interferon genes (STING) and is highly effective in limiting SARS-CoV-2 replication in cells and animals. diABZI-4 inhibited SARS-CoV-2 replication in lung epithelial cells. Administration of diABZI-4 intranasally before or even after virus infection conferred complete protection from severe respiratory disease in K18-ACE2 transgenic mice infected with SARS-CoV-2. Intranasal delivery of diABZI-4 induced a rapid short-lived activation of STING, leading to transient proinflammatory cytokine production and lymphocyte activation in the lung associated with inhibition of viral replication. Our study supports the use of diABZI-4 as a host-directed therapy that mobilizes antiviral defenses for the treatment and prevention of COVID-19.

Treating premature infants with high oxygen is a routine intervention in the context of neonatal intensive care. Unfortunately, the increase in survival rates is associated with various detrimental sequalae of hyperoxia exposure, most notably bronchopulmonary dysplasia (BPD), a disease of disrupted lung development. The effects of high oxygen exposure on other developing organs of the infant, as well as the possible impact such disrupted development may have on later life remain poorly understood. Using a neonatal mouse model to investigate the effects of hyperoxia on the immature immune system we observed a dramatic involution of the thymic medulla, and this lesion was associated with disrupted FoxP3+ regulatory T cell generation and T cell autoreactivity. Significantly, administration of mesenchymal stromal cell-derived extracellular vesicles (MEx) restored thymic medullary architecture and physiological thymocyte profiles. Using single cell transcriptomics, we further demonstrated preferential impact of MEx treatment on the thymic medullary antigen presentation axis, as evidenced by enrichment of antigen presentation and antioxidative-stress related genes in dendritic cells (DCs) and medullary epithelial cells (mTECs). Our study demonstrates that MEx treatment represents a promising restorative therapeutic approach for oxygen-induced thymic injury, thus promoting normal development of both central tolerance and adaptive immunity.

Rationale: Coronavirus disease 2019 (COVID-19) is characterized by increased incidence of microthrombosis with hyperactive platelets sporadically containing viral RNA. It is unclear if SARS-CoV-2 (severe acute respiratory syndrome, corona virus-2) directly alters platelet activation or if these changes are a reaction to infection-mediated global inflammatory alterations. Importantly, the direct effect of SARS-CoV-2 on platelets has yet to be studied. Objective: To characterize the direct SARS-CoV-2–platelet interactions using in vitro studies with purified infectious virions and samples from infected patients. Methods and Results: Platelet RNA analyzed by ARTIC v3 sequencing for SARS-CoV-2 showed presence of fragmented viral genome in all patients with COVID-19. Immunofluorescent imaging of platelets from patients with COVID-19 confirmed presence of SARS-CoV-2 proteins, whereas there was no detection of viral RNA by real-time quantitative polymerase chain reaction. Transmission electron microscopy of platelets incubated with purified SARS-CoV-2 virions demonstrated rapid internalization and digestion leading to distinct morphological changes and resulted in a release of extracellular vesicles. Interactions between SARS-CoV-2 and platelets occurred with or without ACE (angiotensin-converting enzyme) 2 presence as measured by immunofluorescence. Transmission electron microscopy showed that SARS-CoV-2 virions became internalized when they were attached to microparticles, bypassing the need for ACE2. Enrichment analysis of platelet transcriptome from patients with acute COVID-19, compared with those with clinical thrombosis, suggested upregulation of pathways related to virally mediated cell death, specifically necroptosis and apoptosis. Platelets incubated with infectious virus appeared to undergo cell death in 30 minutes postincubation as assessed by transmission electron microscopy and platelets from patients with COVID-19 showed evidence of increased markers of apoptosis and necroptosis by Western blot. Immunofluorescence confirmed colocalization of SARS-CoV-2 with phospho-MLKL (mixed lineage kinase domain-like pseudokinase) and caspase 3 on nonpermeabilized platelets in vitro and in COVID-19 platelets. Conclusions: Platelets internalize SARS-CoV-2 virions, directly or attached to microparticles, and viral internalization leads to rapid digestion, programmed cell death, and extracellular vesicle release. During COVID-19, platelets mediate a rapid response to SARS-CoV-2 and this response can contribute to dysregulated immunity and thrombosis.

Processed single-cell RNA-seq data and R code for the analyses in Guo et al., eLife 2022;10:e70341 The datasets include a Seurat object with raw UMI counts for all genes in all cells; a Seurat object with some cells filtered out that also includes normalized counts, PCA and UMAP coordinates, and clustering results; and pseudobulk counts for each donor in various cell-types that are defined in part by this clustering.