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.

Zebrafish display a distinct ability to regenerate their heart following injury. However, this ability is not shared by another teleost, the medaka. In order to identify cellular and molecular bases for this difference, we performed comparative transcriptomic analyses following cardiac cryoinjury. This comparison points to major differences in immune cell dynamics between these models. Upon closer examination, we observed delayed and reduced macrophage recruitment in medaka, along with delayed neutrophil clearance. To investigate the role of immune responses in cardiac regeneration, we delayed macrophage recruitment in zebrafish and observed compromised neovascularization, neutrophil clearance, cardiomyocyte proliferation and scar resolution. In contrast, stimulating Toll-like receptor signaling in medaka enhanced immune cell dynamics and promoted neovascularization, neutrophil clearance, cardiomyocyte proliferation and scar resolution. Altogether, these data provide further insight into the complex role of the immune response during regeneration, and serve as a platform to identify and test additional regulators of cardiac repair.DOI: http://dx.doi.org/10.7554/eLife.25605.001

Cell-type-specific splicing generates numerous alternatively spliced transcripts playing important roles for organ development and homeostasis, but only a few tissue-specific splicing factors have been identified. We found that RBM24 governs a large number of muscle-specific splicing events that are critically involved in cardiac and skeletal muscle development and disease. Targeted inactivation of RBM24 in mice disrupted cardiac development and impaired sarcomerogenesis in striated muscles. In vitro splicing assays revealed that recombinant RBM24 is sufficient to promote muscle-specific exon inclusion in nuclear extracts of nonmuscle cells. Furthermore, we demonstrate that binding of RBM24 to an intronic splicing enhancer (ISE) is essential and sufficient to overcome repression of exon inclusion by an exonic splicing silencer (ESS) containing PTB and hnRNP A1/A2 binding sites. Introduction of ESS and ISE converted a constitutive exon into an RMB24-dependent alternative exon. We reason that RBM24 is a major regulator of alternative splicing in striated muscles.

Rationale: Activated cardiac fibroblasts (CF) are crucial players in the cardiac damage response; excess fibrosis, however, may result in myocardial stiffening and heart failure development. Inhibition of activated CF has been suggested as a therapeutic strategy in cardiac disease, but whether this truly improves cardiac function is unclear. Objective: To study the effect of CF ablation on cardiac remodeling. Methods and Results: We characterized subgroups of murine CF by single-cell expression analysis and identified periostin as the marker showing the highest correlation to an activated CF phenotype. We generated bacterial artificial chromosome–transgenic mice allowing tamoxifen-inducible Cre expression in periostin-positive cells as well as their diphtheria toxin-mediated ablation. In the healthy heart, periostin expression was restricted to valvular fibroblasts; ablation of this population did not affect cardiac function. After chronic angiotensin II exposure, ablation of activated CF resulted in significantly reduced cardiac fibrosis and improved cardiac function. After myocardial infarction, ablation of periostin-expressing CF resulted in reduced fibrosis without compromising scar stability, and cardiac function was significantly improved. Single-cell transcriptional analysis revealed reduced CF activation but increased expression of prohypertrophic factors in cardiac macrophages and cardiomyocytes, resulting in localized cardiomyocyte hypertrophy. Conclusions: Modulation of the activated CF population is a promising approach to prevent adverse cardiac remodeling in response to angiotensin II and after myocardial infarction.

Supplemental Digital Content is available in the text.

Endothelial cells play a critical role in the adaptation of tissues to injury. Tissue ischemia induced by infarction leads to profound changes in endothelial cell functions and can induce transition to a mesenchymal state. Here we explore the kinetics and individual cellular responses of endothelial cells after myocardial infarction by using single cell RNA sequencing. This study demonstrates a time dependent switch in endothelial cell proliferation and inflammation associated with transient changes in metabolic gene signatures. Trajectory analysis reveals that the majority of endothelial cells 3 to 7 days after myocardial infarction acquire a transient state, characterized by mesenchymal gene expression, which returns to baseline 14 days after injury. Lineage tracing, using the Cdh5-CreERT2;mT/mG mice followed by single cell RNA sequencing, confirms the transient mesenchymal transition and reveals additional hypoxic and inflammatory signatures of endothelial cells during early and late states after injury. These data suggest that endothelial cells undergo a transient mes-enchymal activation concomitant with a metabolic adaptation within the first days after myocardial infarction but do not acquire a long-term mesenchymal fate. This mesenchymal activation may facilitate endothelial cell migration and clonal expansion to regenerate the vascular network.

The cnidarian freshwater polyp Hydra sp. exhibits an unparalleled regeneration capacity in the animal kingdom. Using an integrative transcriptomic and stable isotope labeling by amino acids in cell culture proteomic/phosphoproteomic approach, we studied stem cell-based regeneration in Hydra polyps. As major contributors to head regeneration, we identified diverse signaling pathways adopted for the regeneration response as well as enriched novel genes. Our global analysis reveals two distinct molecular cascades: an early injury response and a subsequent, signaling driven patterning of the regenerating tissue. A key factor of the initial injury response is a general stabilization of proteins and a net upregulation of transcripts, which is followed by a subsequent activation cascade of signaling molecules including Wnts and transforming growth factor (TGF) beta-related factors. We observed moderate overlap between the factors contributing to proteomic and transcriptomic responses suggesting a decoupled regulation between the transcriptional and translational levels. Our data also indicate that interstitial stem cells and their derivatives (e.g., neurons) have no major role in Hydra head regeneration. Remarkably, we found an enrichment of evolutionarily more recent genes in the early regeneration response, whereas conserved genes are more enriched in the late phase. In addition, genes specific to the early injury response were enriched in transposon insertions. Genetic dynamicity and taxon-specific factors might therefore play a hitherto underestimated role in Hydra regeneration.

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.