Multiplex single-cell chemical genomics reveals the kinase dependence of the response to targeted therapy.

Chemical genetic screens are a powerful tool for exploring how cancer cells' response to drugs is shaped by their mutations, yet they lack a molecular view of the contribution of individual genes to the response to exposure. Here, we present sci-Plex-Gene-by-Environment (sci-Plex-GxE), a platform for combined single-cell genetic and chemical screening at scale. We highlight the advantages of large-scale, unbiased screening by defining the contribution of each of 522 human kinases to the response of glioblastoma to different drugs designed to abrogate signaling from the receptor tyrosine kinase pathway. In total, we probed 14,121 gene-by-environment combinations across 1,052,205 single-cell transcriptomes. We identify an expression signature characteristic of compensatory adaptive signaling regulated in a MEK/MAPK-dependent manner. Further analyses aimed at preventing adaptation revealed promising combination therapies, including dual MEK and CDC7/CDK9 or nuclear factor κB (NF-κB) inhibitors, as potent means of preventing transcriptional adaptation of glioblastoma to targeted therapy.

Dystrophin deficiency impairs cell junction formation during embryonic myogenesis.

Mutations in the DMD gene lead to Duchenne muscular dystrophy, a severe X-linked neuromuscular disorder that manifests itself as young boys acquire motor functions. DMD is typically diagnosed at 2 to 4 years of age, but the absence of dystrophin negatively impacts muscle structure and function before overt symptoms appear in patients, which poses a serious challenge in the optimization of standards of care. In this report, we investigated the early consequences of dystrophin deficiency during skeletal muscle development. We used single-cell transcriptome profiling to characterize the myogenic trajectory of human pluripotent stem cells and showed that DMD cells bifurcate to an alternative branch when they reach the somite stage. Here, dystrophin deficiency was linked to marked dysregulations of cell junction protein families involved in the cell state transitions characteristic of embryonic somitogenesis. Altogether, this work demonstrates that in vitro, dystrophin deficiency has deleterious effects on cell-cell communication during myogenic development, which should be considered in future therapeutic strategies for DMD.

Proteostasis governs differential temperature sensitivity across embryonic cell types.

Embryonic development is remarkably robust, but temperature stress can degrade its ability to generate animals with invariant anatomy. Phenotypes associated with environmental stress suggest that some cell types are more sensitive to stress than others, but the basis of this sensitivity is unknown. Here, we characterize hundreds of individual zebrafish embryos under temperature stress using whole-animal single-cell RNA sequencing (RNA-seq) to identify cell types and molecular programs driving phenotypic variability. We find that temperature perturbs the normal proportions and gene expression programs of numerous cell types and also introduces asynchrony in developmental timing. The notochord is particularly sensitive to temperature, which we map to a specialized cell type: sheath cells. These cells accumulate misfolded protein at elevated temperature, leading to a cascading structural failure of the notochord and anatomic defects. Our study demonstrates that whole-animal single-cell RNA-seq can identify mechanisms for developmental robustness and pinpoint cell types that constitute key failure points.

Multiplex single-cell chemical genomics reveals the kinase dependence of the response to targeted therapy.

Chemical genetic screens are a powerful tool for exploring how cancer cells' response to drugs is shaped by their mutations, yet they lack a molecular view of the contribution of individual genes to the response to exposure. Here, we present sci-Plex-Gene-by-Environment (sci-Plex-GxE), a platform for combined single-cell genetic and chemical screening at scale. We highlight the advantages of large-scale, unbiased screening by defining the contribution of each of 522 human kinases to the response of glioblastoma to different drugs designed to abrogate signaling from the receptor tyrosine kinase pathway. In total, we probed 14,121 gene-by-environment combinations across 1,052,205 single-cell transcriptomes. We identify an expression signature characteristic of compensatory adaptive signaling regulated in a MEK/MAPK-dependent manner. Further analyses aimed at preventing adaptation revealed promising combination therapies, including dual MEK and CDC7/CDK9 or NF-kB inhibitors, as potent means of preventing transcriptional adaptation of glioblastoma to targeted therapy.

Differentiation latency and dormancy signatures define fetal liver HSCs at single cell resolution.

Decoding the gene regulatory mechanisms mediating self-renewal of hematopoietic stem cells (HSCs) during their amplification in the fetal liver (FL) is relevant for advancing therapeutic applications aiming to expand transplantable HSCs, a long-standing challenge. Here, to explore intrinsic and extrinsic regulation of self-renewal in FL-HSCs at the single cell level, we engineered a culture platform designed to recapitulate the FL endothelial niche, which supports the amplification of serially engraftable HSCs ex vivo. Leveraging this platform in combination with single cell index flow cytometry, serial transplantation assays, and single cell RNA-sequencing, we elucidated previously unrecognized heterogeneity in immunophenotypically defined FL-HSCs and demonstrated that differentiation latency and transcriptional signatures of biosynthetic dormancy are distinguishing properties of self-renewing FL-HSCs with capacity for serial, long-term multilineage hematopoietic reconstitution. Altogether, our findings provide key insights into HSC expansion and generate a novel resource for future exploration of the intrinsic and niche-derived signaling pathways that support FL-HSC self-renewal.

Engineering a niche supporting hematopoietic stem cell development using integrated single-cell transcriptomics.

Hematopoietic stem cells (HSCs) develop from hemogenic endothelium within embryonic arterial vessels such as the aorta of the aorta-gonad-mesonephros region (AGM). To identify the signals responsible for HSC formation, here we use single cell RNA-sequencing to simultaneously analyze the transcriptional profiles of AGM-derived cells transitioning from hemogenic endothelium to HSCs, and AGM-derived endothelial cells which provide signals sufficient to support HSC maturation and self-renewal. Pseudotemporal ordering reveals dynamics of gene expression during the hemogenic endothelium to HSC transition, identifying surface receptors specifically expressed on developing HSCs. Transcriptional profiling of niche endothelial cells identifies corresponding ligands, including those signaling to Notch receptors, VLA-4 integrin, and CXCR4, which, when integrated in an engineered platform, are sufficient to support the generation of engrafting HSCs. These studies provide a transcriptional map of the signaling interactions necessary for the development of HSCs and advance the goal of engineering HSCs for therapeutic applications.

Multipotent progenitors and hematopoietic stem cells arise independently from hemogenic endothelium in the mouse embryo.

During embryogenesis, waves of hematopoietic progenitors develop from hemogenic endothelium (HE) prior to the emergence of self-renewing hematopoietic stem cells (HSCs). Although previous studies have shown that yolk-sac-derived erythromyeloid progenitors and HSCs emerge from distinct populations of HE, it remains unknown whether the earliest lymphoid-competent progenitors, multipotent progenitors, and HSCs originate from common HE. In this study, we demonstrate by clonal assays and single-cell transcriptomics that rare HE with functional HSC potential in the early murine embryo are distinct from more abundant HE with multilineage hematopoietic potential that fail to generate HSCs. Specifically, HSC-competent HE are characterized by expression of CXCR4 surface marker and by higher expression of genes tied to arterial programs regulating HSC dormancy and self-renewal. Taken together, these findings suggest a revised model of developmental hematopoiesis in which the initial populations of multipotent progenitors and HSCs arise independently from HE with distinct phenotypic and transcriptional properties.

Single-cell landscape of nuclear configuration and gene expression during stem cell differentiation and X inactivation.

BACKGROUND: Mammalian development is associated with extensive changes in gene expression, chromatin accessibility, and nuclear structure. Here, we follow such changes associated with mouse embryonic stem cell differentiation and X inactivation by integrating, for the first time, allele-specific data from these three modalities obtained by high-throughput single-cell RNA-seq, ATAC-seq, and Hi-C. RESULTS: Allele-specific contact decay profiles obtained by single-cell Hi-C clearly show that the inactive X chromosome has a unique profile in differentiated cells that have undergone X inactivation. Loss of this inactive X-specific structure at mitosis is followed by its reappearance during the cell cycle, suggesting a "bookmark" mechanism. Differentiation of embryonic stem cells to follow the onset of X inactivation is associated with changes in contact decay profiles that occur in parallel on both the X chromosomes and autosomes. Single-cell RNA-seq and ATAC-seq show evidence of a delay in female versus male cells, due to the presence of two active X chromosomes at early stages of differentiation. The onset of the inactive X-specific structure in single cells occurs later than gene silencing, consistent with the idea that chromatin compaction is a late event of X inactivation. Single-cell Hi-C highlights evidence of discrete changes in nuclear structure characterized by the acquisition of very long-range contacts throughout the nucleus. Novel computational approaches allow for the effective alignment of single-cell gene expression, chromatin accessibility, and 3D chromosome structure. CONCLUSIONS: Based on trajectory analyses, three distinct nuclear structure states are detected reflecting discrete and profound simultaneous changes not only to the structure of the X chromosomes, but also to that of autosomes during differentiation. Our study reveals that long-range structural changes to chromosomes appear as discrete events, unlike progressive changes in gene expression and chromatin accessibility.

Trajectory analysis quantifies transcriptional plasticity during macrophage polarization.

In recent years, macrophages have been shown to be tremendously plastic in both in vitro and in vivo settings; however, it remains unclear whether macrophages retain any persistent memory of past polarization states which may then impact their future repolarization to new states. Here, we perform deep transcriptomic profiling at high temporal resolution as macrophages are polarized with cytokines that drive them into "M1" and "M2" molecular states. We find through trajectory analysis of their global transcriptomic profiles that macrophages which are first polarized to M1 or M2 and then subsequently repolarized demonstrate little to no memory of their polarization history. We observe complete repolarization both from M1 to M2 and vice versa, and we find that macrophage transcriptional phenotypes are defined by the current cell microenvironment, rather than an amalgamation of past and present states.

Single-Cell Transcriptomic Comparison of Human Fetal Retina, hPSC-Derived Retinal Organoids, and Long-Term Retinal Cultures.

To study the development of the human retina, we use single-cell RNA sequencing (RNA-seq) at key fetal stages and follow the development of the major cell types as well as populations of transitional cells. We also analyze stem cell (hPSC)-derived retinal organoids; although organoids have a very similar cellular composition at equivalent ages as the fetal retina, there are some differences in gene expression of particular cell types. Moreover, the inner retinal lamination is disrupted at more advanced stages of organoids compared with fetal retina. To determine whether the disorganization in the inner retina is due to the culture conditions, we analyze retinal development in fetal retina maintained under similar conditions. These retinospheres develop for at least 6 months, displaying better inner retinal lamination than retinal organoids. Our single-cell RNA sequencing (scRNA-seq) comparisons of fetal retina, retinal organoids, and retinospheres provide a resource for developing better in vitro models for retinal disease.

Massively multiplex chemical transcriptomics at single-cell resolution.

High-throughput chemical screens typically use coarse assays such as cell survival, limiting what can be learned about mechanisms of action, off-target effects, and heterogeneous responses. Here, we introduce "sci-Plex," which uses "nuclear hashing" to quantify global transcriptional responses to thousands of independent perturbations at single-cell resolution. As a proof of concept, we applied sci-Plex to screen three cancer cell lines exposed to 188 compounds. In total, we profiled ~650,000 single-cell transcriptomes across ~5000 independent samples in one experiment. Our results reveal substantial intercellular heterogeneity in response to specific compounds, commonalities in response to families of compounds, and insight into differential properties within families. In particular, our results with histone deacetylase inhibitors support the view that chromatin acts as an important reservoir of acetate in cancer cells.

A Seminomadic Population in Bangladesh with Extensive Exposure to Macaques Does Not Exhibit High Levels of Zoonotic Simian Foamy Virus Infection.

Simian foamy viruses (SVF) are ubiquitous in nonhuman primates (NHP). SFV can be zoonotically transmitted to humans who either work with or live commensally with NHP. We analyzed the blood of 45 Bangladeshi performing monkey owners (an ethnic group called the Bedey) for SFV infection. Surprisingly, a PCR assay failed to detect SFV infection in any of these participants. This is in contrast to our previously reported infection rate of about 5% among Bangladeshi villagers.

New World simian foamy virus infections in vivo and in vitro.

Foamy viruses (FV) are complex retroviruses that naturally infect all nonhuman primates (NHP) studied to date. Zoonotic transmission of Old World NHP simian foamy viruses (SFV) has been documented, leading to nonpathogenic persistent infections. To date, there have been no reports concerning zoonotic transmission of New World monkey (NWM) SFV to humans and resulting infection. In this study, we developed a Western blot assay to detect antibodies to NWM SFV, a nested PCR assay to detect NWM SFV DNA, and a ß-galactosidase-containing indicator cell line to assay replication of NWM SFV. Using these tools, we analyzed the plasma and blood of 116 primatologists, of whom 69 had reported exposures to NWM. While 8 of the primatologists tested were seropositive for SFV from a NWM, the spider monkey, none had detectable levels of viral DNA in their blood. We found that SFV isolated from three different species of NWM replicated in some, but not all, human cell lines. From our data, we conclude that while humans exposed to NWM SFV produce antibodies, there is no evidence for long-term viral persistence.

Population dynamics of rhesus macaques and associated foamy virus in Bangladesh.

Foamy viruses are complex retroviruses that have been shown to be transmitted from nonhuman primates to humans. In Bangladesh, infection with simian foamy virus (SFV) is ubiquitous among rhesus macaques, which come into contact with humans in diverse locations and contexts throughout the country. We analyzed microsatellite DNA from 126 macaques at six sites in Bangladesh in order to characterize geographic patterns of macaque population structure. We also included in this study 38 macaques owned by nomadic people who train them to perform for audiences. PCR was used to analyze a portion of the proviral gag gene from all SFV-positive macaques, and multiple clones were sequenced. Phylogenetic analysis was used to infer long-term patterns of viral transmission. Analyses of SFV gag gene sequences indicated that macaque populations from different areas harbor genetically distinct strains of SFV, suggesting that geographic features such as forest cover play a role in determining the dispersal of macaques and SFV. We also found evidence suggesting that humans traveling the region with performing macaques likely play a role in the translocation of macaques and SFV. Our studies found that individual animals can harbor more than one strain of SFV and that presence of more than one SFV strain is more common among older animals. Some macaques are infected with SFV that appears to be recombinant. These findings paint a more detailed picture of how geographic and sociocultural factors influence the spectrum of simian-borne retroviruses.

Zoonotic simian foamy virus in Bangladesh reflects diverse patterns of transmission and co-infection.

Simian foamy viruses (SFVs) are ubiquitous in non-human primates (NHPs). As in all retroviruses, reverse transcription of SFV leads to recombination and mutation. Because more humans have been shown to be infected with SFV than with any other simian borne virus, SFV is a potentially powerful model for studying the virology and epidemiology of viruses at the human/NHP interface. In Asia, SFV is likely transmitted to humans through macaque bites and scratches that occur in the context of everyday life. We analyzed multiple proviral sequences from the SFV gag gene from both humans and macaques in order to characterize retroviral transmission at the human/NHP interface in Bangladesh. Here we report evidence that humans can be concurrently infected with multiple SFV strains, with some individuals infected by both an autochthonous SFV strain as well as a strain similar to SFV found in macaques from another geographic area. These data, combined with previous results, suggest that both human-facilitated movement of macaques leading to the introduction of non-resident strains of SFV and retroviral recombination in macaques contribute to SFV diversity among humans in Bangladesh.

Expression of prototype foamy virus pol as a Gag-Pol fusion protein does not change the timing of reverse transcription.

Foamy viruses are retroviruses whose Pol protein is synthesized without Gag from a spliced mRNA. Unlike orthoretroviruses, reverse transcription occurs during viral assembly, leading to DNA-containing virions. When prototype foamy virus Pol is expressed as an orthoretroviral-like Gag-Pol fusion protein, reverse transcription also occurs late in viral replication, as measured by the timing of reverse transcriptase sensitivity to the inhibitor 3'-azido-3'deoxythymidine (AZT). Thus, timing of reverse transcription is intrinsic to Pol itself.

Foamy virus Pol protein expressed as a Gag-Pol fusion retains enzymatic activities, allowing for infectious virus production.

Foamy viruses (FV) synthesize Pol from a spliced pol mRNA independently of Gag, unlike orthoretroviruses, which synthesize Pol as a Gag-Pol protein that coassembles with Gag. We found that prototype FV (PFV) mutants expressing Gag and Pol only as a Gag-Pol protein without the spliced Pol contain protease activity equivalent to that of wild-type (WT) Pol. Regardless of the presence or absence of the spliced Pol, the PFV Gag-Pol proteins can assemble into virus-like particles (VLPs), in contrast to the orthoretroviral Gag-Pol proteins, which cannot form VLPs. However, the PFV Gag-Pol VLPs have aberrant morphologies and are not infectious. In the absence of the spliced Pol, coexpression of a PFV Gag-Pol protein with Gag can produce infectious virions. Our results suggest that enzymes encoded by PFV pol (protease, reverse transcriptase, and integrase) are enzymatically active if they are synthesized as part of a Gag-Pol protein.

The DEAD-box RNA helicase DDX6 is required for efficient encapsidation of a retroviral genome.

Viruses have to encapsidate their own genomes during the assembly process. For most RNA viruses, there are sequences within the viral RNA and virion proteins needed for high efficiency of genome encapsidation. However, the roles of host proteins in this process are not understood. Here we find that the cellular DEAD-box RNA helicase DDX6 is required for efficient genome packaging of foamy virus, a spumaretrovirus. After infection, a significant amount of DDX6, normally concentrated in P bodies and stress granules, re-localizes to the pericentriolar site where viral RNAs and Gag capsid proteins are concentrated and capsids are assembled. Knockdown of DDX6 by siRNA leads to a decreased level of viral nucleic acids in extracellular particles, although viral protein expression, capsid assembly and release, and accumulation of viral RNA and Gag protein at the assembly site are little affected. DDX6 does not interact stably with Gag proteins nor is it incorporated into particles. However, we find that the ATPase/helicase motif of DDX6 is essential for viral replication. This suggests that the ATP hydrolysis and/or the RNA unwinding activities of DDX6 function in moderating the viral RNA conformation and/or viral RNA-Gag ribonucleoprotein complex in a transient manner to facilitate incorporation of the viral RNA into particles. These results reveal a unique role for a highly conserved cellular protein of RNA metabolism in specifically re-locating to the site of viral assembly for its function as a catalyst in retroviral RNA packaging.