Search

Analyses of transcriptional bursting from single-cell RNA-sequencing data have revealed patterns of variation and regulation in the kinetic parameters that could be inferred. Here we profiled newly transcribed (4-thiouridine-labelled) RNA across 10,000 individual primary mouse fibroblasts to more broadly infer bursting kinetics and coordination. We demonstrate that inference from new RNA profiles could separate the kinetic parameters that together specify the burst size, and that the synthesis rate (and not the transcriptional off rate) controls the burst size. Importantly, transcriptome-wide inference of transcriptional on and off rates provided conclusive evidence that RNA polymerase II transcribes genes in bursts. Recent reports identified examples of transcriptional co-bursting, yet no global analyses have been performed. The deep new RNA profiles we generated with allelic resolution demonstrated that co-bursting rarely appears more frequently than expected by chance, except for certain gene pairs, notably paralogues located in close genomic proximity. Altogether, new RNA single-cell profiling critically improves the inference of transcriptional bursting and provides strong evidence for independent transcriptional bursting of mammalian genes., (© 2024. The Author(s).)

Recent work employed Pearson residuals from Poisson or negative binomial models to normalize UMI data. To extend this approach to non-UMI data, we model the additional amplification step with a compound distribution: we assume that sequenced RNA molecules follow a negative binomial distribution, and are then replicated following an amplification distribution. We show how this model leads to compound Pearson residuals, which yield meaningful gene selection and embeddings of Smart-seq2 datasets. Further, we suggest that amplification distributions across several sequencing protocols can be described by a broken power law. The resulting compound model captures previously unexplained overdispersion and zero-inflation patterns in non-UMI data.

Rare multipotent stem cells replenish millions of blood cells per second through a time-consuming process, passing through multiple stages of increasingly lineage-restricted progenitors. Although insults to the blood-forming system highlight the need for more rapid blood replenishment from stem cells, established models of hematopoiesis implicate only one mandatory differentiation pathway for each blood cell lineage. Here, we establish a nonhierarchical relationship between distinct stem cells that replenish all blood cell lineages and stem cells that replenish almost exclusively platelets, a lineage essential for hemostasis and with important roles in both the innate and adaptive immune systems. These distinct stem cells use cellularly, molecularly and functionally separate pathways for the replenishment of molecularly distinct megakaryocyte-restricted progenitors: a slower steady-state multipotent pathway and a fast-track emergency-activated platelet-restricted pathway. These findings provide a framework for enhancing platelet replenishment in settings in which slow recovery of platelets remains a major clinical challenge., (© 2024. The Author(s).)

Cell types can be classified according to shared patterns of transcription. Non-genetic variability among individual cells of the same type has been ascribed to stochastic transcriptional bursting and transient cell states. Using high-coverage single-cell RNA profiling, we asked whether long-term, heritable differences in gene expression can impart diversity within cells of the same type. Studying clonal human lymphocytes and mouse brain cells, we uncovered a vast diversity of heritable gene expression patterns among different clones of cells of the same type in vivo. We combined chromatin accessibility and RNA profiling on different lymphocyte clones to reveal thousands of regulatory regions exhibiting interclonal variation, which could be directly linked to interclonal variation in gene expression. Our findings identify a source of cellular diversity, which may have important implications for how cellular populations are shaped by selective processes in development, aging, and disease. A record of this paper's transparent peer review process is included in the supplemental information., Competing Interests: Declaration of interests J.E.M., C.-J.E., J. Lundeberg, and J.F. receive compensation for consultation from 10× Genomics Inc on topics unrelated to this study. J.E.M. also receives consulting fees from Applied Molecular Transport Inc. on topics unrelated to this study., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

Myelodysplastic syndromes with ring sideroblasts (MDS-RS) commonly develop from hematopoietic stem cells (HSC) bearing mutations in the splicing factor SF3B1 (SF3B1mt). Direct studies into MDS-RS pathobiology have been limited by a lack of model systems that fully recapitulate erythroid biology and RS development and the inability to isolate viable human RS. Here, we combined successful direct RS isolation from patient samples, high-throughput multiomics analysis of cells encompassing the SF3B1mt stem-erythroid continuum, and functional assays to investigate the impact of SF3B1mt on erythropoiesis and RS accumulation. The isolated RS differentiated, egressed into the blood, escaped traditional nonsense-mediated decay (NMD) mechanisms, and leveraged stress-survival pathways that hinder wild-type hematopoiesis through pathogenic GDF15 overexpression. Importantly, RS constituted a contaminant of magnetically enriched CD34+ cells, skewing bulk transcriptomic data. Mis-splicing in SF3B1mt cells was intensified by erythroid differentiation through accelerated RNA splicing and decreased NMD activity, and SF3B1mt led to truncations in several MDS-implicated genes. Finally, RNA mis-splicing induced an uncoupling of RNA and protein expression, leading to critical abnormalities in proapoptotic p53 pathway genes. Overall, this characterization of erythropoiesis in SF3B1mt RS provides a resource for studying MDS-RS and uncovers insights into the unexpectedly active biology of the "dead-end" RS., Significance: Ring sideroblast isolation combined with state-of-the-art multiomics identifies survival mechanisms underlying SF3B1-mutant erythropoiesis and establishes an active role for erythroid differentiation and ring sideroblasts themselves in SF3B1-mutant myelodysplastic syndrome pathogenesis., (©2023 The Authors; Published by the American Association for Cancer Research.)

Background: Formation of tissue-specific transcriptional programs underlies multicellular development, including dorsoventral (DV) patterning of the Drosophila embryo. This involves interactions between transcriptional enhancers and promoters in a chromatin context, but how the chromatin landscape influences transcription is not fully understood., Results: Here we comprehensively resolve differential transcriptional and chromatin states during Drosophila DV patterning. We find that RNA Polymerase II pausing is established at DV promoters prior to zygotic genome activation (ZGA), that pausing persists irrespective of cell fate, but that release into productive elongation is tightly regulated and accompanied by tissue-specific P-TEFb recruitment. DV enhancers acquire distinct tissue-specific chromatin states through CBP-mediated histone acetylation that predict the transcriptional output of target genes, whereas promoter states are more tissue-invariant. Transcriptome-wide inference of burst kinetics in different cell types revealed that while DV genes are generally characterized by a high burst size, either burst size or frequency can differ between tissues., Conclusions: The data suggest that pausing is established by pioneer transcription factors prior to ZGA and that release from pausing is imparted by enhancer chromatin state to regulate bursting in a tissue-specific manner in the early embryo. Our results uncover how developmental patterning is orchestrated by tissue-specific bursts of transcription from Pol II primed promoters in response to enhancer regulatory cues., (© 2023. The Author(s).)

Tissues within an organism and even cell types within a tissue can age with different velocities. However, it is unclear whether cells of one type experience different aging trajectories within a tissue depending on their spatial location. Here, we used spatial transcriptomics in combination with single-cell ATAC-seq and RNA-seq, lipidomics and functional assays to address how cells in the male murine liver are affected by age-related changes in the microenvironment. Integration of the datasets revealed zonation-specific and age-related changes in metabolic states, the epigenome and transcriptome. The epigenome changed in a zonation-dependent manner and functionally, periportal hepatocytes were characterized by decreased mitochondrial fitness, whereas pericentral hepatocytes accumulated large lipid droplets. Together, we provide evidence that changing microenvironments within a tissue exert strong influences on their resident cells that can shape epigenetic, metabolic and phenotypic outputs., (© 2023. The Author(s).)

Current single-cell RNA sequencing (scRNA-seq) methods with high cellular throughputs sacrifice full-transcript coverage and often sensitivity. Here we describe Smart-seq3xpress, which miniaturizes and streamlines the Smart-seq3 protocol to substantially reduce reagent use and increase cellular throughput. Smart-seq3xpress analysis of peripheral blood mononuclear cells resulted in a granular atlas complete with common and rare cell types. Compared with droplet-based single-cell RNA sequencing that sequences RNA ends, the additional full-transcript coverage revealed cell-type-associated isoform variation., (© 2022. The Author(s).)

The variable clinical course of follicular lymphoma (FL) is determined by the molecular heterogeneity of tumor cells and complex interactions within the tumor microenvironment (TME). IL-4 producing follicular helper T cells (T FH ) are critical components of the FL TME. Binding of IL-4 to IL-4R on FL cells activates JAK/STAT signaling. We identified STAT6 mutations (STAT6 MUT ) in 13% of FL (N = 33/258), all clustered within the DNA binding domain. Gene expression data and immunohistochemistry showed upregulation of IL-4/STAT6 target genes in STAT6 MUT FL, including CCL17, CCL22, and FCER2 (CD23). Functionally, STAT6 MUT was gain-of-function by serial replating phenotype in pre-B CFU assays. Expression of STAT6 MUT enhanced IL-4 induced FCER2/CD23, CCL17 and CCL22 expression and was associated with nuclear accumulation of pSTAT6. RNA sequencing identified PARP14 -a transcriptional switch and co-activator of STAT6- among the top differentially upregulated genes in IL-4 stimulated STAT6 MUT lymphoma cells and in STAT6 MUT primary FL cells. Quantitative chromatin immunoprecipitation (qChIP) demonstrated binding of STAT6 MUT but not STAT6 WT to the PARP14 promotor. Reporter assays showed increased IL-4 induced transactivation activity of STAT6 MUT at the PARP14 promotor, suggesting a self-reinforcing regulatory circuit. Knock-down of PARP14 or PARP-inhibition abrogated the STAT6 MUT gain-of-function phenotype. Thus, our results identify PARP14 as a novel therapeutic target in STAT6 MUT FL., (© 2022. The Author(s).)

Plasmacytoid and conventional dendritic cells (pDC and cDC) are generated from progenitor cells in the bone marrow and commitment to pDCs or cDC subtypes may occur in earlier and later progenitor stages. Cells within the CD11c + MHCII -/lo Siglec-H + CCR9 lo DC precursor fraction of the mouse bone marrow generate both pDCs and cDCs. Here we investigate the heterogeneity and commitment of subsets in this compartment by single-cell transcriptomics and high-dimensional flow cytometry combined with cell fate analysis: Within the CD11c + MHCII -/lo Siglec-H + CCR9 lo DC precursor pool cells expressing high levels of Ly6D and lacking expression of transcription factor Zbtb46 contain CCR9 lo B220 hi immediate pDC precursors and CCR9 lo B220 lo (lo-lo) cells which still generate pDCs and cDCs in vitro and in vivo under steady state conditions. cDC-primed cells within the Ly6D hi Zbtb46 - lo-lo precursors rapidly upregulate Zbtb46 and pass through a Zbtb46 + Ly6D + intermediate stage before acquiring cDC phenotype after cell division. Type I IFN stimulation limits cDC and promotes pDC output from this precursor fraction by arresting cDC-primed cells in the Zbtb46 + Ly6D + stage preventing their expansion and differentiation into cDCs. Modulation of pDC versus cDC output from precursors by external factors may allow for adaptation of DC subset composition at later differentiation stages., (© 2022. The Author(s).)

Single-cell sequencing methods rely on molecule-counting strategies to account for amplification biases, yet no experimental strategy to evaluate counting performance exists. Here, we introduce molecular spikes-RNA spike-ins containing built-in unique molecular identifiers (UMIs) that we use to identify critical experimental and computational conditions for accurate RNA counting in single-cell RNA-sequencing (scRNA-seq). Using molecular spikes, we uncovered impaired RNA counting in methods that were not informative for cellular RNA abundances due to inflated UMI counts. We further leverage molecular spikes to improve estimates of total endogenous RNA amounts in cells, and introduce a strategy to correct experiments with impaired RNA counting. The molecular spikes and the accompanying R package UMIcountR ( https://github.com/cziegenhain/UMIcountR ) will improve the validation of new methods, better estimate and adjust for cellular mRNA amounts and enable more indepth characterization of RNA counting in scRNA-seq., (© 2022. The Author(s).)

Cost-efficient library generation by early barcoding has been central in propelling single-cell RNA sequencing. Here, we optimize and validate prime-seq, an early barcoding bulk RNA-seq method. We show that it performs equivalently to TruSeq, a standard bulk RNA-seq method, but is fourfold more cost-efficient due to almost 50-fold cheaper library costs. We also validate a direct RNA isolation step, show that intronic reads are derived from RNA, and compare cost-efficiencies of available protocols. We conclude that prime-seq is currently one of the best options to set up an early barcoding bulk RNA-seq protocol from which many labs would profit., (© 2022. The Author(s).)

An increasing number of long noncoding RNAs (lncRNAs) have experimentally confirmed functions, yet little is known about their transcriptional dynamics and it is challenging to determine their regulatory effects. Here, we used allele-sensitive single-cell RNA sequencing to demonstrate that, compared to messenger RNAs, lncRNAs have twice as long duration between two transcriptional bursts. Additionally, we observed increased cell-to-cell variability in lncRNA expression due to lower frequency bursting producing larger numbers of RNA molecules. Exploiting heterogeneity in asynchronously growing cells, we identified and experimentally validated lncRNAs with cell state-specific functions involved in cell cycle progression and apoptosis. Finally, we identified cis-functioning lncRNAs and showed that knockdown of these lncRNAs modulated the nearby protein-coding gene's transcriptional burst frequency or size. In summary, we identified distinct transcriptional regulation of lncRNAs and demonstrated a role for lncRNAs in the regulation of mRNA transcriptional bursting., (© 2022. The Author(s).)

The risks associated with re-identification of human genetic data are severely limiting open data sharing in life sciences, even in studies where donor-related genetic variant information is not of primary interest. Here, we developed BAMboozle, a versatile tool to eliminate critical types of sensitive genetic information in human sequence data by reverting aligned reads to the genome reference sequence. Applying BAMboozle to functional genomics data, such as single-cell RNA-seq (scRNA-seq) and scATAC-seq datasets, confirmed the removal of donor-related single nucleotide polymorphisms (SNPs) and indels in a manner that did not disclose the altered positions. Importantly, BAMboozle only removes the genetic sequence variants of the sample (i.e., donor) while preserving other important aspects of the raw sequence data. For example, BAMboozled scRNA-seq data contained accurate cell-type associated gene expression signatures, splice kinetic information, and can be used for methods benchmarking. Altogether, BAMboozle efficiently removes genetic variation in aligned sequence data, which represents a step forward towards open data sharing in many areas of genomics where the genetic variant information is not of primary interest., (© 2021. The Author(s).)

Transcriptional bursts render substantial biological noise in cellular transcriptomes. Here, we investigated the theoretical extent of allelic expression resulting from transcriptional bursting and how it compared to the amount biallelic, monoallelic and allele-biased expression observed in single-cell RNA-sequencing (scRNA-seq) data. We found that transcriptional bursting can explain the allelic expression patterns observed in single cells, including the frequent observations of autosomal monoallelic gene expression. Importantly, we identified that the burst frequency largely determined the fraction of cells with monoallelic expression, whereas the burst size had little effect on monoallelic observations. The high consistency between the bursting model predictions and scRNA-seq observations made it possible to assess the heterogeneity of a group of cells as their deviation in allelic observations from the expected. Finally, both burst frequency and size contributed to allelic imbalance observations and reinforced that studies of allelic imbalance can be confounded from the inherent noise in transcriptional bursting. Altogether, we demonstrate that allele-level transcriptional bursting renders widespread, although predictable, amounts of monoallelic and biallelic expression in single cells and cell populations., Competing Interests: The authors have declared that no competing interests exist.

A detailed understanding of intestinal stem cell (ISC) self-renewal and differentiation is required to treat chronic intestinal diseases. However, the different models of ISC lineage hierarchy 1-6 and segregation 7-12 are subject to debate. Here, we have discovered non-canonical Wnt/planar cell polarity (PCP)-activated ISCs that are primed towards the enteroendocrine or Paneth cell lineage. Strikingly, integration of time-resolved lineage labelling with single-cell gene expression analysis revealed that both lineages are directly recruited from ISCs via unipotent transition states, challenging the existence of formerly predicted bi- or multipotent secretory progenitors 7-12 . Transitory cells that mature into Paneth cells are quiescent and express both stem cell and secretory lineage genes, indicating that these cells are the previously described Lgr5 + label-retaining cells 7 . Finally, Wnt/PCP-activated Lgr5 + ISCs are molecularly indistinguishable from Wnt/β-catenin-activated Lgr5 + ISCs, suggesting that lineage priming and cell-cycle exit is triggered at the post-transcriptional level by polarity cues and a switch from canonical to non-canonical Wnt/PCP signalling. Taken together, we redefine the mechanisms underlying ISC lineage hierarchy and identify the Wnt/PCP pathway as a new niche signal preceding lateral inhibition in ISC lineage priming and segregation.

Genome-wide DNA demethylation is a unique feature of mammalian development and naïve pluripotent stem cells. Here, we describe a recently evolved pathway in which global hypomethylation is achieved by the coupling of active and passive demethylation. TET activity is required, albeit indirectly, for global demethylation, which mostly occurs at sites devoid of TET binding. Instead, TET-mediated active demethylation is locus-specific and necessary for activating a subset of genes, including the naïve pluripotency and germline marker Dppa3 (Stella, Pgc7). DPPA3 in turn drives large-scale passive demethylation by directly binding and displacing UHRF1 from chromatin, thereby inhibiting maintenance DNA methylation. Although unique to mammals, we show that DPPA3 alone is capable of inducing global DNA demethylation in non-mammalian species (Xenopus and medaka) despite their evolutionary divergence from mammals more than 300 million years ago. Our findings suggest that the evolution of Dppa3 facilitated the emergence of global DNA demethylation in mammals.

Large-scale sequencing of RNA from individual cells can reveal patterns of gene, isoform and allelic expression across cell types and states 1 . However, current short-read single-cell RNA-sequencing methods have limited ability to count RNAs at allele and isoform resolution, and long-read sequencing techniques lack the depth required for large-scale applications across cells 2,3 . Here we introduce Smart-seq3, which combines full-length transcriptome coverage with a 5' unique molecular identifier RNA counting strategy that enables in silico reconstruction of thousands of RNA molecules per cell. Of the counted and reconstructed molecules, 60% could be directly assigned to allelic origin and 30-50% to specific isoforms, and we identified substantial differences in isoform usage in different mouse strains and human cell types. Smart-seq3 greatly increased sensitivity compared to Smart-seq2, typically detecting thousands more transcripts per cell. We expect that Smart-seq3 will enable large-scale characterization of cell types and states across tissues and organisms.

Single-cell RNA sequencing (scRNA-seq) is the leading technique for characterizing the transcriptomes of individual cells in a sample. The latest protocols are scalable to thousands of cells and are being used to compile cell atlases of tissues, organs and organisms. However, the protocols differ substantially with respect to their RNA capture efficiency, bias, scale and costs, and their relative advantages for different applications are unclear. In the present study, we generated benchmark datasets to systematically evaluate protocols in terms of their power to comprehensively describe cell types and states. We performed a multicenter study comparing 13 commonly used scRNA-seq and single-nucleus RNA-seq protocols applied to a heterogeneous reference sample resource. Comparative analysis revealed marked differences in protocol performance. The protocols differed in library complexity and their ability to detect cell-type markers, impacting their predictive value and suitability for integration into reference cell atlases. These results provide guidance both for individual researchers and for consortium projects such as the Human Cell Atlas.

Single-nucleotide polymorphisms and locus amplification link the NF-κB transcription factor c-Rel to human autoimmune diseases and B cell lymphomas, respectively. However, the functional consequences of enhanced c-Rel levels remain enigmatic. Here, we overexpressed c-Rel specifically in mouse B cells from BAC-transgenic gene loci and demonstrate that c-Rel protein levels linearly dictated expansion of germinal center B (GCB) cells and isotype-switched plasma cells. c-Rel expression in B cells of otherwise c-Rel-deficient mice fully rescued terminal B cell differentiation, underscoring its critical B cell-intrinsic roles. Unexpectedly, in GCB cells transcription-independent regulation produced the highest c-Rel protein levels among B cell subsets. In c-Rel-overexpressing GCB cells this caused enhanced nuclear translocation, a profoundly altered transcriptional program, and increased proliferation. Finally, we provide a link between c-Rel gain and autoimmunity by showing that c-Rel overexpression in B cells caused autoantibody production and renal immune complex deposition.

The recent rapid spread of single cell RNA sequencing (scRNA-seq) methods has created a large variety of experimental and computational pipelines for which best practices have not yet been established. Here, we use simulations based on five scRNA-seq library protocols in combination with nine realistic differential expression (DE) setups to systematically evaluate three mapping, four imputation, seven normalisation and four differential expression testing approaches resulting in ~3000 pipelines, allowing us to also assess interactions among pipeline steps. We find that choices of normalisation and library preparation protocols have the biggest impact on scRNA-seq analyses. Specifically, we find that library preparation determines the ability to detect symmetric expression differences, while normalisation dominates pipeline performance in asymmetric DE-setups. Finally, we illustrate the importance of informed choices by showing that a good scRNA-seq pipeline can have the same impact on detecting a biological signal as quadrupling the sample size.

Actin is a protein of central importance for many cellular key processes. It is regulated by local interactions with a large number of actin binding proteins (ABPs). Various compounds are known to either increase or decrease the polymerization dynamics of actin. However, no actin binding compound has been developed for clinical applications yet because of selectivity issues. We provide a crystal structure of the natural product chivosazole A (ChivoA) bound to actin and show that-in addition to inhibiting nucleation, polymerization, and severing of F-actin filaments-it selectively modulates binding of ABPs to G-actin: Although unphysiological actin dimers are induced by ChivoA, interaction with gelsolin, profilin, cofilin, and thymosin-β4 is inhibited. Moreover, ChivoA causes transcriptional effects differing from latrunculin B, an actin binder with a different binding site. Our data show that ChivoA and related compounds could serve as scaffolds for the development of actin binding molecules selectively targeting specific actin functions.

Actin binding compounds are widely used tools in cell biology. We compare the biological and biochemical effects of miuraenamide A and jasplakinolide, a structurally related prototypic actin stabilizer. Though both compounds have similar effects on cytoskeletal morphology and proliferation, they affect migration and transcription in a distinctive manner, as shown by a transcriptome approach in endothelial cells. In vitro, miuraenamide A acts as an actin nucleating, F-actin polymerizing and stabilizing compound, just like described for jasplakinolide. However, in contrast to jasplakinolide, miuraenamide A competes with cofilin, but not gelsolin or Arp2/3 for binding to F-actin. We propose a binding mode of miuraenamide A, explaining both its similarities and its differences to jasplakinolide. Molecular dynamics simulations suggest that the bromophenol group of miurenamide A interacts with residues Tyr133, Tyr143, and Phe352 of actin. This shifts the D-loop of the neighboring actin, creating tighter packing of the monomers, and occluding the binding site of cofilin. Since relatively small changes in the molecular structure give rise to this selectivity, actin binding compounds surprisingly are promising scaffolds for creating actin binders with specific functionality instead of just "stabilizers".

Genetic disruptions of the forkhead box transcription factor FOXP2 in humans cause an autosomal-dominant speech and language disorder. While FOXP2 expression pattern are highly conserved, its role in specific brain areas for mammalian social behaviors remains largely unknown. Here we studied mice carrying a homozygous cortical Foxp2 deletion. The postnatal development and gross morphological architecture of mutant mice was indistinguishable from wildtype (WT) littermates. Unbiased behavioral profiling of adult mice revealed abnormalities in approach behavior towards conspecifics as well as in the reciprocal responses of WT interaction partners. Furthermore mutant mice showed alterations in acoustical parameters of ultrasonic vocalizations, which also differed in function of the social context. Cell type-specific gene expression profiling of cortical pyramidal neurons revealed aberrant regulation of genes involved in social behavior. In particular Foxp2 mutants showed the downregulation of Mint2 (Apba2), a gene involved in approach behavior in mice and autism spectrum disorder in humans. Taken together these data demonstrate that cortical Foxp2 is required for normal social behaviors in mice., (© The Author(s) 2018. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Bacterial opportunistic pathogens are feared for their difficult-to-treat nosocomial infections and for causing morbidity in immunocompromised patients. Here, we study how such a versatile opportunist, Pseudomonas aeruginosa, adapts to conditions inside and outside its model host Caenorhabditis elegans, and use phenotypic and genotypic screens to identify the mechanistic basis of virulence evolution. We found that virulence significantly dropped in unstructured environments both in the presence and absence of the host, but remained unchanged in spatially structured environments. Reduction of virulence was either driven by a substantial decline in the production of siderophores (in treatments without hosts) or toxins and proteases (in treatments with hosts). Whole-genome sequencing of evolved clones revealed positive selection and parallel evolution across replicates, and showed an accumulation of mutations in regulator genes controlling virulence factor expression. Our study identifies the spatial structure of the non-host environment as a key driver of virulence evolution in an opportunistic pathogen.

Actin has emerged as a versatile regulator of gene transcription. Cytoplasmatic actin regulates mechanosensitive-signaling pathways such as MRTF-SRF and Hippo-YAP/TAZ. In the nucleus, both polymerized and monomeric actin directly interfere with transcription-associated molecular machineries. Natural actin-binding compounds are frequently used tools to study actin-related processes in cell biology. However, their influence on transcriptional regulation and intranuclear actin polymerization is poorly understood to date. Here, we analyze the effects of two representative actin-binding compounds, Miuraenamide A (polymerizing properties) and Latrunculin B (depolymerizing properties), on transcriptional regulation in primary cells. We find that actin stabilizing and destabilizing compounds inversely shift nuclear actin levels without a direct influence on polymerization state and intranuclear aspects of transcriptional regulation. Furthermore, we identify Miuraenamide A as a potent inducer of G-actin-dependent SRF target gene expression. In contrast, the F-actin-regulated Hippo-YAP/TAZ axis remains largely unaffected by compound-induced actin aggregation. This is due to the inability of AMOTp130 to bind to the amorphous actin aggregates resulting from treatment with miuraenamide. We conclude that actin-binding compounds predominantly regulate transcription via their influence on cytoplasmatic G-actin levels, while transcriptional processes relying on intranuclear actin polymerization or functional F-actin networks are not targeted by these compounds at tolerable doses.

Single-cell RNA sequencing (scRNA-seq) has emerged as a central genome-wide method to characterize cellular identities and processes. Consequently, improving its sensitivity, flexibility, and cost-efficiency can advance many research questions. Among the flexible plate-based methods, single-cell RNA barcoding and sequencing (SCRB-seq) is highly sensitive and efficient. Here, we systematically evaluate experimental conditions of this protocol and find that adding polyethylene glycol considerably increases sensitivity by enhancing cDNA synthesis. Furthermore, using Terra polymerase increases efficiency due to a more even cDNA amplification that requires less sequencing of libraries. We combined these and other improvements to develop a scRNA-seq library protocol we call molecular crowding SCRB-seq (mcSCRB-seq), which we show to be one of the most sensitive, efficient, and flexible scRNA-seq methods to date.

Single-cell RNA sequencing (scRNA-seq) is currently transforming our understanding of biology, as it is a powerful tool to resolve cellular heterogeneity and molecular networks. Over 50 protocols have been developed in recent years and also data processing and analyzes tools are evolving fast. Here, we review the basic principles underlying the different experimental protocols and how to benchmark them. We also review and compare the essential methods to process scRNA-seq data from mapping, filtering, normalization and batch corrections to basic differential expression analysis. We hope that this helps to choose appropriate experimental and computational methods for the research question at hand.

Background: Single-cell RNA-sequencing (scRNA-seq) experiments typically analyze hundreds or thousands of cells after amplification of the cDNA. The high throughput is made possible by the early introduction of sample-specific bar codes (BCs), and the amplification bias is alleviated by unique molecular identifiers (UMIs). Thus, the ideal analysis pipeline for scRNA-seq data needs to efficiently tabulate reads according to both BC and UMI., Findings: zUMIs is a pipeline that can handle both known and random BCs and also efficiently collapse UMIs, either just for exon mapping reads or for both exon and intron mapping reads. If BC annotation is missing, zUMIs can accurately detect intact cells from the distribution of sequencing reads. Another unique feature of zUMIs is the adaptive downsampling function that facilitates dealing with hugely varying library sizes but also allows the user to evaluate whether the library has been sequenced to saturation. To illustrate the utility of zUMIs, we analyzed a single-nucleus RNA-seq dataset and show that more than 35% of all reads map to introns. Also, we show that these intronic reads are informative about expression levels, significantly increasing the number of detected genes and improving the cluster resolution., Conclusions: zUMIs flexibility makes if possible to accommodate data generated with any of the major scRNA-seq protocols that use BCs and UMIs and is the most feature-rich, fast, and user-friendly pipeline to process such scRNA-seq data.

Developmental processes, such as angiogenesis, are associated with a constant remodeling of the actin cytoskeleton in response to different mechanical stimuli. The mechanosensitive transcription factors MRTF-A (MKL1) and YAP (also known as YAP1) are important mediators of this challenging adaptation process. However, it is as yet unknown whether both pathways respond in an identical or in a divergent manner to a given microenvironmental guidance cue. Here, we use a micropatterning approach to dissect single aspects of cellular behavior in a spatiotemporally controllable setting. Using the exemplary process of angiogenesis, we show that cell-cell contacts and adhesive surface area are shared regulatory parameters of MRTF and YAP on rigid 2D surfaces. By analyzing MRTF and YAP under laminar flow conditions and during cell migration on dumbbell-shaped microstructures, we demonstrate that they exhibit different translocation kinetics. In conclusion, our work promotes the application of micropatterning techniques as a cell biological tool to study mechanosensitive signaling in the context of angiogenesis., Competing Interests: Competing interestsThe authors declare no competing or financial interests., (© 2018. Published by The Company of Biologists Ltd.)

The poor correlation of mutational landscapes with phenotypes limits our understanding of the pathogenesis and metastasis of pancreatic ductal adenocarcinoma (PDAC). Here we show that oncogenic dosage-variation has a critical role in PDAC biology and phenotypic diversification. We find an increase in gene dosage of mutant KRAS in human PDAC precursors, which drives both early tumorigenesis and metastasis and thus rationalizes early PDAC dissemination. To overcome the limitations posed to gene dosage studies by the stromal richness of PDAC, we have developed large cell culture resources of metastatic mouse PDAC. Integration of cell culture genomes, transcriptomes and tumour phenotypes with functional studies and human data reveals additional widespread effects of oncogenic dosage variation on cell morphology and plasticity, histopathology and clinical outcome, with the highest Kras MUT levels underlying aggressive undifferentiated phenotypes. We also identify alternative oncogenic gains (Myc, Yap1 or Nfkb2), which collaborate with heterozygous Kras MUT in driving tumorigenesis, but have lower metastatic potential. Mechanistically, different oncogenic gains and dosages evolve along distinct evolutionary routes, licensed by defined allelic states and/or combinations of hallmark tumour suppressor alterations (Cdkn2a, Trp53, Tgfβ-pathway). Thus, evolutionary constraints and contingencies direct oncogenic dosage gain and variation along defined routes to drive the early progression of PDAC and shape its downstream biology. Our study uncovers universal principles of Ras-driven oncogenesis that have potential relevance beyond pancreatic cancer.

To elucidate the molecular basis of BMP4-induced differentiation of human pluripotent stem cells (PSCs) toward progeny with trophectoderm characteristics, we produced transcriptome, epigenome H3K4me3, H3K27me3, and CpG methylation maps of trophoblast progenitors, purified using the surface marker APA. We combined them with the temporally resolved transcriptome of the preprogenitor phase and of single APA+ cells. This revealed a circuit of bivalent TFAP2A, TFAP2C, GATA2, and GATA3 transcription factors, coined collectively the "trophectoderm four" (TEtra), which are also present in human trophectoderm in vivo. At the onset of differentiation, the TEtra factors occupy multiple sites in epigenetically inactive placental genes and in OCT4 Functional manipulation of GATA3 and TFAP2A indicated that they directly couple trophoblast-specific gene induction with suppression of pluripotency. In accordance, knocking down GATA3 in primate embryos resulted in a failure to form trophectoderm. The discovery of the TEtra circuit indicates how trophectoderm commitment is regulated in human embryogenesis., Competing Interests: The authors declare no conflict of interest., (Copyright © 2017 the Author(s). Published by PNAS.)

Summary: Power analysis is essential to optimize the design of RNA-seq experiments and to assess and compare the power to detect differentially expressed genes in RNA-seq data. PowsimR is a flexible tool to simulate and evaluate differential expression from bulk and especially single-cell RNA-seq data making it suitable for a priori and posterior power analyses., Availability and Implementation: The R package and associated tutorial are freely available at https://github.com/bvieth/powsimR., Contact: vieth@bio.lmu.de or hellmann@bio.lmu.de., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author 2017. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com)

Effectively targeting leukemia-initiating cells (LIC) in FLT3 -ITD-mutated acute myeloid leukemia (AML) is crucial for cure. Tyrosine kinase inhibitors (TKI) have limited impact as single agents, failing to eradicate LIC in the bone marrow. Using primary AML samples and a patient-derived xenograft model, we investigated whether combining the FLT3-selective TKI crenolanib with the hypomethylating agent azacitidine (AZA) eliminates FLT3 -ITD LIC and whether efficacy of this combination depends on co-existing mutations. Using multiparameter flow cytometry, we show FLT3 -ITD occurs within the most primitive Lin - /CD33 (+) /CD45 dim /CD34 + CD38 - LIC compartment. Crenolanib alone could not target FLT3 -ITD LIC in contact with niche cells while addition of AZA overcame stromal protection resulting in dramatically reduced clonogenic capacity of LIC in vitro and severely impaired engraftment in NSG mice. Strikingly, FLT3 -mutated samples harboring TET2 mutations were completely resistant to crenolanib whereas neither NPM1 nor DNMT3A mutations influenced response. Conversely, primary AML LIC harboring either TET2, DNMT3A or NPM1 mutations did not show increased sensitivity to AZA. In summary, resistance of FLT3 -ITD LIC to TKI depends on co-existing epigenetic mutations. However, AZA + crenolanib effectively abrogates stromal protection and inhibits survival of FLT3 -ITD LIC irrespective of mutations, providing evidence for this combination as a means to suppress residual LIC., Competing Interests: CONFLICTS OF INTEREST KSG received honoraria from Celgene. All other authors declare no potential conflicts of interest.

We identify SMARCD2 (SWI/SNF-related, matrix-associated, actin-dependent regulator of chromatin, subfamily D, member 2), also known as BAF60b (BRG1/Brahma-associated factor 60b), as a critical regulator of myeloid differentiation in humans, mice, and zebrafish. Studying patients from three unrelated pedigrees characterized by neutropenia, specific granule deficiency, myelodysplasia with excess of blast cells, and various developmental aberrations, we identified three homozygous loss-of-function mutations in SMARCD2. Using mice and zebrafish as model systems, we showed that SMARCD2 controls early steps in the differentiation of myeloid-erythroid progenitor cells. In vitro, SMARCD2 interacts with the transcription factor CEBPɛ and controls expression of neutrophil proteins stored in specific granules. Defective expression of SMARCD2 leads to transcriptional and chromatin changes in acute myeloid leukemia (AML) human promyelocytic cells. In summary, SMARCD2 is a key factor controlling myelopoiesis and is a potential tumor suppressor in leukemia.

Single-cell RNA sequencing (scRNA-seq) offers new possibilities to address biological and medical questions. However, systematic comparisons of the performance of diverse scRNA-seq protocols are lacking. We generated data from 583 mouse embryonic stem cells to evaluate six prominent scRNA-seq methods: CEL-seq2, Drop-seq, MARS-seq, SCRB-seq, Smart-seq, and Smart-seq2. While Smart-seq2 detected the most genes per cell and across cells, CEL-seq2, Drop-seq, MARS-seq, and SCRB-seq quantified mRNA levels with less amplification noise due to the use of unique molecular identifiers (UMIs). Power simulations at different sequencing depths showed that Drop-seq is more cost-efficient for transcriptome quantification of large numbers of cells, while MARS-seq, SCRB-seq, and Smart-seq2 are more efficient when analyzing fewer cells. Our quantitative comparison offers the basis for an informed choice among six prominent scRNA-seq methods, and it provides a framework for benchmarking further improvements of scRNA-seq protocols., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Here, we show that the Wnt5a-haploinsufficient niche regenerates dysfunctional HSCs, which do not successfully engraft in secondary recipients. RNA sequencing of the regenerated donor Lin - SCA-1 + KIT + (LSK) cells shows dysregulated expression of ZEB1-associated genes involved in the small GTPase-dependent actin polymerization pathway. Misexpression of DOCK2, WAVE2, and activation of CDC42 results in apolar F-actin localization, leading to defects in adhesion, migration and homing of HSCs regenerated in a Wnt5a-haploinsufficient microenvironment. Moreover, these cells show increased differentiation in vitro, with rapid loss of HSC-enriched LSK cells. Our study further shows that the Wnt5a-haploinsufficient environment similarly affects BCR-ABL p185 leukemia-initiating cells, which fail to generate leukemia in 42% of the studied recipients, or to transfer leukemia to secondary hosts. Thus, we show that WNT5A in the bone marrow niche is required to regenerate HSCs and leukemic cells with functional ability to rearrange the actin cytoskeleton and engraft successfully., (© 2017 Schreck et al.)

Tumor relapse is associated with dismal prognosis, but responsible biological principles remain incompletely understood. To isolate and characterize relapse-inducing cells, we used genetic engineering and proliferation-sensitive dyes in patient-derived xenografts of acute lymphoblastic leukemia (ALL). We identified a rare subpopulation that resembled relapse-inducing cells with combined properties of long-term dormancy, treatment resistance, and stemness. Single-cell and bulk expression profiling revealed their similarity to primary ALL cells isolated from pediatric and adult patients at minimal residual disease (MRD). Therapeutically adverse characteristics were reversible, as resistant, dormant cells became sensitive to treatment and started proliferating when dissociated from the in vivo environment. Our data suggest that ALL patients might profit from therapeutic strategies that release MRD cells from the niche., (Copyright © 2016 The Authors. Published by Elsevier Inc. All rights reserved.)

Currently, quantitative RNA-seq methods are pushed to work with increasingly small starting amounts of RNA that require amplification. However, it is unclear how much noise or bias amplification introduces and how this affects precision and accuracy of RNA quantification. To assess the effects of amplification, reads that originated from the same RNA molecule (PCR-duplicates) need to be identified. Computationally, read duplicates are defined by their mapping position, which does not distinguish PCR- from natural duplicates and hence it is unclear how to treat duplicated reads. Here, we generate and analyse RNA-seq data sets prepared using three different protocols (Smart-Seq, TruSeq and UMI-seq). We find that a large fraction of computationally identified read duplicates are not PCR duplicates and can be explained by sampling and fragmentation bias. Consequently, the computational removal of duplicates does improve neither accuracy nor precision and can actually worsen the power and the False Discovery Rate (FDR) for differential gene expression. Even when duplicates are experimentally identified by unique molecular identifiers (UMIs), power and FDR are only mildly improved. However, the pooling of samples as made possible by the early barcoding of the UMI-protocol leads to an appreciable increase in the power to detect differentially expressed genes.

Objective: Macroglial cells like astrocytes are key targets for the formation of HIV-1 reservoirs in the brain. The 'shock-and-kill' HIV-1 cure strategy proposes eradication of reservoirs by clinical treatment with latency reversing agents (LRAs). However, virus activation may endanger the brain, due to limited cell turnover, viral neurotoxicity and poor penetration of antiretroviral drugs. Since the brain is not accessible to clinical sampling, we established an experimental model to investigate the LRA effects on HIV-1 latency in macroglial reservoirs., Design: Human neural stem cells (HNSC.100) were used to generate a system that models HIV-1 transcriptional latency in proliferating progenitor, as well as differentiated macroglial cell populations and latency-modulating effects of LRAs and compounds targeting HIV-1 transcription were analysed., Methods: HNSCs were infected with pseudotyped Env-defective HIV-1 viruses. HIV-1 DNA and RNA levels were quantified by qPCR. Expression of latent GFP-reporter viruses was analysed by confocal microscopy and flow cytometry. NF-κB signalling was investigated by confocal microscopy and chromatin immunoprecipitation., Results: Two of the eight well known LRAs (tumour necrosis factor-alpha, suberoylanilide hydroxamic acid) reactivated HIV-1 in latently infected HNSCs. Tumour necrosis factor-alpha reactivated HIV-1 in progenitor and differentiated populations, whereas suberoylanilide hydroxamic acid was more potent in progenitors. Pre-treatment with inhibitors of key HIV-1 transcription factors (NF-κB, Cdk9) suppressed HIV-1 reactivation., Conclusion: We conclude that latent HIV-1 in macroglial reservoirs can be activated by selected LRAs. Identification of small molecules that suppress HIV-1 reactivation supports functional cure strategies. We propose using the HNSC model to develop novel strategies to enforce provirus quiescence in the brain.