Herpesvirus infection initiates a range of perturbations in the host cell, which remain poorly understood at the level of individual cells. Here, we quantify the transcriptome of single human primary fibroblasts during the first hours of lytic infection with HSV-1. By applying a generalizable analysis scheme, we define a precise temporal order of early viral gene expression and propose a set-wise emergence of viral genes. We identify host cell genes and pathways relevant for infection by combining three different computational approaches: gene and pathway overdispersion analysis, prediction of cell-state transition probabilities, as well as future cell states. One transcriptional program, which correlates with increased resistance to infection, implicates the transcription factor NRF2. Consequently, Bardoxolone methyl and Sulforaphane, two known NRF2 agonists, impair virus production, suggesting that NRF2 activation restricts viral infection. Our study provides insights into early stages of HSV-1 infection and serves as a general blueprint for the investigation of heterogeneous cell states in virus infection.

Circular RNAs (circRNAs) are a regulatory RNA class. While cancer-driving functions have been identified for single circRNAs, how they modulate gene expression in cancer is not well understood. We investigate circRNA expression in the pediatric malignancy, neuroblastoma, through deep whole-transcriptome sequencing in 104 primary neuroblastomas covering all risk groups. We demonstrate that MYCN amplification, which defines a subset of high-risk cases, causes globally suppressed circRNA biogenesis directly dependent on the DHX9 RNA helicase. We detect similar mechanisms in shaping circRNA expression in the pediatric cancer medulloblastoma implying a general MYCN effect. Comparisons to other cancers identify 25 circRNAs that are specifically upregulated in neuroblastoma, including circARID1A. Transcribed from the ARID1A tumor suppressor gene, circARID1A promotes cell growth and survival, mediated by direct interaction with the KHSRP RNA-binding protein. Our study highlights the importance of MYCN regulating circRNAs in cancer and identifies molecular mechanisms, which explain their contribution to neuroblastoma pathogenesis.

The extent of translational control of gene expression in mammalian tissues remains largely unknown. Here we perform genome-wide RNA sequencing and ribosome profiling in heart and liver tissues to investigate strain-specific translational regulation in the spontaneously hypertensive rat (SHR/Ola). For the most part, transcriptional variation is equally apparent at the translational level and there is limited evidence of translational buffering. Remarkably, we observe hundreds of strain-specific differences in translation, almost doubling the number of differentially expressed genes. The integration of genetic, transcriptional and translational data sets reveals distinct signatures in 3'UTR variation, RNA-binding protein motifs and miRNA expression associated with translational regulation of gene expression. We show that a large number of genes associated with heart and liver traits in human genome-wide association studies are primarily translationally regulated. Capturing interindividual differences in the translated genome will lead to new insights into the genes and regulatory pathways underlying disease phenotypes.

Leigh syndrome (LS) is a severe manifestation of mitochondrial disease in children and is currently incurable. The lack of effective models hampers our understanding of the mechanisms underlying the neuronal pathology of LS. Using patient-derived induced pluripotent stem cells and CRISPR/Cas9 engineering, we developed a human model of LS caused by mutations in the complex IV assembly gene SURF1. Single-cell RNA-sequencing and multi-omics analysis revealed compromised neuronal morphogenesis in mutant neural cultures and brain organoids. The defects emerged at the level of neural progenitor cells (NPCs), which retained a glycolytic proliferative state that failed to instruct neuronal morphogenesis. LS NPCs carrying mutations in the complex I gene NDUFS4 recapitulated morphogenesis defects. SURF1 gene augmentation and PGC1A induction via bezafibrate treatment supported the metabolic programming of LS NPCs, leading to restored neuronal morphogenesis. Our findings provide mechanistic insights and suggest potential interventional strategies for a rare mitochondrial disease.

Characterizing the complex composition of solid tumors is fundamental for understanding tumor initiation, progression and metastasis. While patient-derived samples provide valuable insight, they are heterogeneous on multiple molecular levels, and often originate from advanced tumor stages. Here, we use single-cell transcriptome and epitope profiling together with pathway and lineage analyses to study tumorigenesis from a developmental perspective in a mouse model of salivary gland squamous cell carcinoma. We provide a comprehensive cell atlas and characterize tumor-specific cells. We find that these cells are connected along a reproducible developmental trajectory: initiated in basal cells exhibiting an epithelial-to-mesenchymal transition signature, tumorigenesis proceeds through Wnt-differential cancer stem cell-like subpopulations before differentiating into luminal-like cells. Our work provides unbiased insights into tumor-specific cellular identities in a whole tissue environment, and emphasizes the power of using defined genetic model systems.

The cellular organization of gastrointestinal crypts is orchestrated by different cells of the stromal niche but available in vitro models fail to fully recapitulate the interplay between epithelium and stroma. Here, we establish a colon assembloid system comprising the epithelium and diverse stromal cell subtypes. These assembloids recapitulate the development of mature crypts resembling in vivo cellular diversity and organization, including maintenance of a stem/progenitor cell compartment in the base and their maturation into secretory/absorptive cell types. This process is supported by self-organizing stromal cells around the crypts that resemble in vivo organization, with cell types that support stem cell turnover adjacent to the stem cell compartment. Assembloids that lack BMP receptors either in epithelial or stromal cells fail to undergo proper crypt formation. Our data highlight the crucial role of bidirectional signaling between epithelium and stroma, with BMP as a central determinant of compartmentalization along the crypt axis.

Skeletal muscle stem cells, called satellite cells and defined by the transcription factor PAX7, are responsible for postnatal muscle growth, homeostasis and regeneration. Attempts to utilize the regenerative potential of muscle stem cells for therapeutic purposes so far failed. We previously established the existence of human PAX7-positive cell colonies with high regenerative potential. We now identified PAX7-negative human muscle-derived cell colonies also positive for the myogenic markers desmin and MYF5. These include cells from a patient with a homozygous PAX7 c.86-1G > A mutation (PAX7null). Single cell and bulk transcriptome analysis show high intra- and inter-donor heterogeneity and reveal the endothelial cell marker CLEC14A to be highly expressed in PAX7null cells. All PAX7-negative cell populations, including PAX7null, form myofibers after transplantation into mice, and regenerate muscle after reinjury. Transplanted PAX7neg cells repopulate the satellite cell niche where they re-express PAX7, or, strikingly, CLEC14A. In conclusion, transplanted human cells do not depend on PAX7 for muscle regeneration.

Protein subcellular localization is fundamental to the establishment of the body axis, cell migration, synaptic plasticity, and a vast range of other biological processes. Protein localization occurs through three mechanisms: protein transport, mRNA localization, and local translation. However, the relative contribution of each process to neuronal polarity remains unknown. Using neurons differentiated from mouse embryonic stem cells, we analyze protein and RNA expression and translation rates in isolated cell bodies and neurites genome-wide. We quantify 7323 proteins and the entire transcriptome, and identify hundreds of neurite-localized proteins and locally translated mRNAs. Our results demonstrate that mRNA localization is the primary mechanism for protein localization in neurites that may account for half of the neurite-localized proteome. Moreover, we identify multiple neurite-targeted non-coding RNAs and RNA-binding proteins with potential regulatory roles. These results provide further insight into the mechanisms underlying the establishment of neuronal polarity.Subcellular localization of RNAs and proteins is important for polarized cells such as neurons. Here the authors differentiate mouse embryonic stem cells into neurons, and analyze the local transcriptome, proteome, and translated transcriptome in their cell bodies and neurites, providing a unique resource for future studies on neuronal polarity.

Understanding regulation of an mRNA requires knowledge of its regulators. However, methods for reliable de-novo identification of proteins binding to a particular RNA are scarce and were thus far only successfully applied to abundant noncoding RNAs in cell culture. Here, we present vIPR, an RNA-protein crosslink, RNA pulldown, and shotgun proteomics approach to identify proteins bound to selected mRNAs in C. elegans. Applying vIPR to the germline-specific transcript gld-1 led to enrichment of known and novel interactors. By comparing enrichment upon gld-1 and lin-41 pulldown, we demonstrate that vIPR recovers both common and specific RNA-binding proteins, and we validate DAZ-1 as a specific gld-1 regulator. Finally, combining vIPR with small RNA sequencing, we recover known and biologically important transcript-specific miRNA interactions, and we identify miR-84 as a specific interactor of the gld-1 transcript. We envision that vIPR will provide a platform for investigating RNA in vivo regulation in diverse biological systems.

Octopuses are mollusks that have evolved intricate neural systems comparable with vertebrates in terms of cell number, complexity and size. The brain cell types that control their sophisticated behavioral repertoire are still unknown. Here, we profile the cell diversity of the paralarval Octopus vulgaris brain to build a cell type atlas that comprises mostly neural cells, but also multiple glial subtypes, endothelial cells and fibroblasts. We spatially map cell types to the vertical, subesophageal and optic lobes. Investigation of cell type conservation reveals a shared gene signature between glial cells of mouse, fly and octopus. Genes related to learning and memory are enriched in vertical lobe cells, which show molecular similarities with Kenyon cells in Drosophila. We construct a cell type taxonomy revealing transcriptionally related cell types, which tend to appear in the same brain region. Together, our data sheds light on cell type diversity and evolution in the octopus brain.