While footprinting analysis of ATAC-seq data can theoretically enable investigation of transcription factor (TF) binding, the lack of a computational tool able to conduct different levels of footprinting analysis has so-far hindered the widespread application of this method. Here we present TOBIAS, a comprehensive, accurate, and fast footprinting framework enabling genome-wide investigation of TF binding dynamics for hundreds of TFs simultaneously. We validate TOBIAS using paired ATAC-seq and ChIP-seq data, and find that TOBIAS outperforms existing methods for bias correction and footprinting. As a proof-ofconcept, we illustrate how TOBIAS can unveil complex TF dynamics during zygotic genome activation in both humans and mice, and propose how zygotic Dux activates cascades of TFs, binds to repeat elements and induces expression of novel genetic elements.

Supplemental Digital Content is available in the text.

Formation and segregation of cell lineages forming the heart have been studied extensively but the underlying gene regulatory networks and epigenetic changes driving cell fate transitions during early cardiogenesis are still only partially understood. Here, we comprehensively characterize mouse cardiac progenitor cells (CPCs) marked by Nkx2-5 and Isl1 expression from E7.5 to E9.5 using single-cell RNA sequencing and transposase-accessible chromatin profiling (ATAC-seq). By leveraging on cell-to-cell transcriptome and chromatin accessibility heterogeneity, we identify different previously unknown cardiac subpopulations. Reconstruction of developmental trajectories reveal that multipotent Isl1+ CPC pass through an attractor state before separating into different developmental branches, whereas extended expression of Nkx2-5 commits CPC to an unidirectional cardiomyocyte fate. Furthermore, we show that CPC fate transitions are associated with distinct open chromatin states critically depending on Isl1 and Nkx2-5. Our data provide a model of transcriptional and epigenetic regulations during cardiac progenitor cell fate decisions at single-cell resolution.

BackgroundNotophthalmus viridescens, an urodelian amphibian, represents an excellent model organism to study regenerative processes, but mechanistic insights into molecular processes driving regeneration have been hindered by a paucity and poor annotation of coding nucleotide sequences. The enormous genome size and the lack of a closely related reference genome have so far prevented assembly of the urodelian genome.ResultsWe describe the de novo assembly of the transcriptome of the newt Notophthalmus viridescens and its experimental validation. RNA pools covering embryonic and larval development, different stages of heart, appendage and lens regeneration, as well as a collection of different undamaged tissues were used to generate sequencing datasets on Sanger, Illumina and 454 platforms. Through a sequential de novo assembly strategy, hybrid datasets were converged into one comprehensive transcriptome comprising 120,922 non-redundant transcripts with a N50 of 975. From this, 38,384 putative transcripts were annotated and around 15,000 transcripts were experimentally validated as protein coding by mass spectrometry-based proteomics. Bioinformatical analysis of coding transcripts identified 826 proteins specific for urodeles. Several newly identified proteins establish novel protein families based on the presence of new sequence motifs without counterparts in public databases, while others containing known protein domains extend already existing families and also constitute new ones.ConclusionsWe demonstrate that our multistep assembly approach allows de novo assembly of the newt transcriptome with an annotation grade comparable to well characterized organisms. Our data provide the groundwork for mechanistic experiments to answer the question whether urodeles utilize proprietary sets of genes for tissue regeneration.

The long-chain fatty acid receptor FFAR1 is highly expressed in pancreatic β-cells. Synthetic FFAR1 agonists can be used as antidiabetic drugs to promote glucose-stimulated insulin secretion (GSIS). However, the physiological role of FFAR1 in β-cells remains poorly understood. Here we show that 20-HETE activates FFAR1 and promotes GSIS via FFAR1 with higher potency and efficacy than dietary fatty acids such as palmitic, linoleic, and α-linolenic acid. Murine and human β-cells produce 20-HETE, and the ω-hydroxylase-mediated formation and release of 20-HETE is strongly stimulated by glucose. Pharmacological inhibition of 20-HETE formation and blockade of FFAR1 in islets inhibits GSIS. In islets from type-2 diabetic humans and mice, glucose-stimulated 20-HETE formation and 20-HETE-dependent stimulation of GSIS are strongly reduced. We show that 20-HETE is an FFAR1 agonist, which functions as an autocrine positive feed-forward regulator of GSIS, and that a reduced glucose-induced 20-HETE formation contributes to inefficient GSIS in type-2 diabetes.

Damage-induced fibrotic scarring limits tissue regeneration in mammals and is a leading cause of morbidity. In contrast, species like zebrafish can regenerate damaged tissues without excessive fibrosis. However, whether specific signaling pathways can both limit fibrosis and promote regeneration is unclear. Here, we show that interleukin-11 (Il-11)/Stat3 signaling has such a dual function. Zebrafish lacking Il-11 receptor function display severely compromised heart, fin, and scale regeneration. Deep phenotyping and transcriptional analysis of adult hearts and fins show that Il-11 signaling drives cellular reprogramming to orchestrate global and tissue-specific regenerative programs and broadly antagonizes hallmarks of adult mammalian scarring. Mechanistically, our data indicate that IL-11 signaling in endothelial cells antagonizes profibrotic transforming growth factor- signaling and endothelial-to-mesenchymal transition, limiting scarring and promoting cardiomyocyte repopulation, after injury. Overall, our findings position damage-induced Il-11/Stat3 signaling in a key role limiting fibrosis and promoting regeneration, revealing novel targets for regenerative therapies.

G-protein-coupled receptor (GPCR) expression is extensively studied in bulk cDNA, but heterogeneity and functional patterning of GPCR expression in individual vascular cells is poorly understood. Here, we perform a microfluidic-based single-cell GPCR expression analysis in primary smooth muscle cells (SMC) and endothelial cells (EC). GPCR expression is highly heterogeneous in all cell types, which is confirmed in reporter mice, on the protein level and in human cells. Inflammatory activation in murine models of sepsis or atherosclerosis results in characteristic changes in the GPCR repertoire, and we identify functionally relevant subgroups of cells that are characterized by specific GPCR patterns. We further show that dedifferentiating SMC upregulate GPCRs such as Gpr39, Gprc5b, Gprc5c or Gpr124, and that selective targeting of Gprc5b modulates their differentiation state. Taken together, single-cell profiling identifies receptors expressed on pathologically relevant subpopulations and provides a basis for the development of new therapeutic strategies in vascular diseases.

Individual adult ventricular cardiomyocytes are either mono- or multi-nucleated and undergo morphological changes during cardiac hypertrophy. However, corresponding transcriptional signatures, reflecting potentially different functions or the ability for cell-cycle entry, are not known. The aim of this study was to determine the transcriptional profile of mono- and multi-nucleated adult cardiomyocytes by single-cell RNA-sequencing (scRNA-seq) and to investigate heterogeneity among cardiomyocytes under baseline conditions and in pressure-induced cardiac hypertrophy. We developed an array-based approach for scRNA-seq of rod-shaped multi-nucleated cardiomyocytes from both healthy and hypertrophic hearts. Single-cell transcriptomes of mono- or multi-nucleated cardiomyocytes were highly similar, although a certain degree of variation was noted across both populations. Non-image-based quality control allowing inclusion of damaged cardiomyocytes generated artificial cell clusters demonstrating the need for strict exclusion criteria. In contrast, cardiomyocytes isolated from hypertrophic heart after transverse aortic constriction showed heterogeneous transcriptional signatures, characteristic for hypoxia-induced responses. Immunofluorescence analysis revealed an inverse correlation between HIF1α + cells and CD31-stained vessels, suggesting that imbalanced vascular growth in the hypertrophied heart induces cellular heterogeneity. Our study demonstrates that individual mono- and multi-nucleated cardiomyocytes express nearly identical sets of genes. Homogeneity among cardiomyocytes was lost after induction of hypertrophy due to differential HIF1α-dependent responses most likely caused by none-homogenous vessel growth. Electronic supplementary material The online version of this article (10.1007/s00395-019-0744-z) contains supplementary material, which is available to authorized users.