Antigen presentation plays positive roles in the regenerative response to cardiac injury in zebrafish.

In contrast to adult mammals, adult zebrafish can fully regenerate injured cardiac tissue, and this regeneration process requires an adequate and tightly controlled immune response. However, which components of the immune response are required during regeneration is unclear. Here, we report positive roles for the antigen presentation-adaptive immunity axis during zebrafish cardiac regeneration. We find that following the initial innate immune response, activated endocardial cells (EdCs), as well as immune cells, start expressing antigen presentation genes. We also observe that T helper cells, a.k.a. Cd4+ T cells, lie in close physical proximity to these antigen-presenting EdCs. We targeted Major Histocompatibility Complex (MHC) class II antigen presentation by generating cd74a; cd74b mutants, which display a defective immune response. In these mutants, Cd4+ T cells and activated EdCs fail to efficiently populate the injured tissue and EdC proliferation is significantly decreased. cd74a; cd74b mutants exhibit additional defects in cardiac regeneration including reduced cardiomyocyte dedifferentiation and proliferation. Notably, Cd74 also becomes activated in neonatal mouse EdCs following cardiac injury. Altogether, these findings point to positive roles for antigen presentation during cardiac regeneration, potentially involving interactions between activated EdCs, classical antigen-presenting cells, and Cd4+ T cells.

Proteomic and transcriptomic characterisation of FIA10, a novel murine leukemic cell line that metastasizes into the brain.

Brain metastasis leads to increased mortality and is a major site of relapse for several cancers, yet the molecular mechanisms of brain metastasis are not well understood. In this study, we established and characterized a new leukemic cell line, FIA10, that metastasizes into the central nervous system (CNS) following injection into the tail vein of syngeneic mice. Mice injected with FIA10 cells developed neurological symptoms such as loss of balance, tremor, ataxic gait and seizures, leading to death within 3 months. Histopathology coupled with PCR analysis clearly showed infiltration of leukemic FIA10 cells into the brain parenchyma of diseased mice, with little involvement of bone marrow, peripheral blood and other organs. To define pathways that contribute to CNS metastasis, global transcriptome and proteome analysis was performed on FIA10 cells and compared with that of the parental stem cell line FDCP-Mix and the related FIA18 cells, which give rise to myeloid leukemia without CNS involvement. 188 expressed genes (RNA level) and 189 proteins were upregulated (log2 ratio FIA10/FIA18 ≥ 1) and 120 mRNAs and 177 proteins were downregulated (log2 ratio FIA10/FIA18 ≤ 1) in FIA10 cells compared with FIA18 cells. Major upregulated pathways in FIA10 cells revealed by biofunctional analyses involved immune response components, adhesion molecules and enzymes implicated in extracellular matrix remodeling, opening and crossing the blood-brain barrier (BBB), molecules supporting migration within the brain parenchyma, alterations in metabolism necessary for growth within the brain microenvironment, and regulators for these functions. Downregulated RNA and protein included several tumor suppressors and DNA repair enzymes. In line with the function of FIA10 cells to specifically infiltrate the brain, FIA10 cells have acquired a phenotype that permits crossing the BBB and adapting to the brain microenvironment thereby escaping immune surveillance. These data and our model system FIA10 will be valuable resources to study the occurrence of brain metastases and may help in the development of potential therapies against brain invasion.

HIF1α-AS1 is a DNA:DNA:RNA triplex-forming lncRNA interacting with the HUSH complex.

DNA:DNA:RNA triplexes that are formed through Hoogsteen base-pairing of the RNA in the major groove of the DNA duplex have been observed in vitro, but the extent to which these interactions occur in cells and how they impact cellular functions remains elusive. Using a combination of bioinformatic techniques, RNA/DNA pulldown and biophysical studies, we set out to identify functionally important DNA:DNA:RNA triplex-forming long non-coding RNAs (lncRNA) in human endothelial cells. The lncRNA HIF1α-AS1 was retrieved as a top hit. Endogenous HIF1α-AS1 reduces the expression of numerous genes, including EPH Receptor A2 and Adrenomedullin through DNA:DNA:RNA triplex formation by acting as an adapter for the repressive human silencing hub complex (HUSH). Moreover, the oxygen-sensitive HIF1α-AS1 is down-regulated in pulmonary hypertension and loss-of-function approaches not only result in gene de-repression but also enhance angiogenic capacity. As exemplified here with HIF1α-AS1, DNA:DNA:RNA triplex formation is a functionally important mechanism of trans-acting gene expression control.

i2dash: Creation of Flexible, Interactive, and Web-based Dashboards for Visualization of Omics Data.

Data visualization and interactive data exploration are important aspects of illustrating complex concepts and results from analyses of omics data. A suitable visualization has to be intuitive and accessible. Web-based dashboards have become popular tools for the arrangement, consolidation, and display of such visualizations. However, the combination of automated data processing pipelines handling omics data and dynamically generated, interactive dashboards is poorly solved. Here, we present i2dash, an R package intended to encapsulate functionality for the programmatic creation of customized dashboards. It supports interactive and responsive (linked) visualizations across a set of predefined graphical layouts. i2dash addresses the needs of data analysts/software developers for a tool that is compatible and attachable to any R-based analysis pipeline, thereby fostering the separation of data visualization on one hand and data analysis tasks on the other hand. In addition, the generic design of i2dash enables the development of modular extensions for specific needs. As a proof of principle, we provide an extension of i2dash optimized for single-cell RNA sequencing analysis, supporting the creation of dashboards for the visualization needs of such experiments. Equipped with these features, i2dash is suitable for extensive use in large-scale sequencing/bioinformatics facilities. Along this line, we provide i2dash as a containerized solution, enabling a straightforward large-scale deployment and sharing of dashboards using cloud services. i2dash is freely available via the R package archive CRAN (https://CRAN.R-project.org/package=i2dash).

Interleukin-11 signaling promotes cellular reprogramming and limits fibrotic scarring during tissue regeneration.

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.

multicrispr: gRNA design for prime editing and parallel targeting of thousands of targets.

Targeting the coding genome to introduce nucleotide deletions/insertions via the CRISPR/Cas9 technology has become a standard procedure. It has quickly spawned a multitude of methods such as prime editing, APEX proximity labeling, or homology directed repair, for which supporting bioinformatics tools are, however, lagging behind. New CRISPR/Cas9 applications often require specific gRNA design functionality, and a generic tool is critically missing. Here, we introduce multicrispr, an R/bioconductor tool, intended to design individual gRNAs and complex gRNA libraries. The package is easy to use; detects, scores, and filters gRNAs on both efficiency and specificity; visualizes and aggregates results per target or CRISPR/Cas9 sequence; and finally returns both genomic ranges and sequences of gRNAs. To be generic, multicrispr defines and implements a genomic arithmetic framework as a basis for facile adaptation to techniques recently introduced such as prime editing or yet to arise. Its performance and design concepts such as target set-specific filtering render multicrispr a tool of choice when dealing with screening-like approaches.

ATAC-seq footprinting unravels kinetics of transcription factor binding during zygotic genome activation.

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-of-concept, 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.

Mono- and multi-nucleated ventricular cardiomyocytes constitute a transcriptionally homogenous cell population.

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.

Connect-four: genomic analyses of regenerating stem cells identifies zygotic Dux factors as tumor initiators.

How, if and in which cell types embryonic gene expression programs are elicited to induce tumor formation remains poorly understood. Through genomic analyses of regenerating, p53 deficient muscle stem cells we identified various oncogenomic amplifications, including but not limited to, the zygotic transcription factor Duxbl/DUXB to initiate tumorigenic transformation.

A whole organism small molecule screen identifies novel regulators of pancreatic endocrine development.

An early step in pancreas development is marked by the expression of the transcription factor Pdx1 within the pancreatic endoderm, where it is required for the specification of all endocrine cell types. Subsequently, Pdx1 expression becomes restricted to the ß-cell lineage, where it plays a central role in ß-cell function. This pivotal role of Pdx1 at various stages of pancreas development makes it an attractive target to enhance pancreatic ß-cell differentiation and increase ß-cell function. In this study, we used a newly generated zebrafish reporter to screen over 8000 small molecules for modulators of pdx1 expression. We found four hit compounds and validated their efficacy at different stages of pancreas development. Notably, valproic acid treatment increased pancreatic endoderm formation, while inhibition of TGFß signaling led to α-cell to ß-cell transdifferentiation. HC toxin, another HDAC inhibitor, enhances ß-cell function in primary mouse and human islets. Thus, using a whole organism screening strategy, this study identified new pdx1 expression modulators that can be used to influence different steps in pancreas and ß-cell development.

The histone demethylase PHF8 facilitates alternative splicing of the histocompatibility antigen HLA-G.

Histone3-lysine9 (H3K9) residues not only control gene expression, but also contribute to RNA splicing. Here, the H3K9 histone demethylase PHF8 was investigated in endothelial cells for its involvement in alternative splicing. An angiogenic sprouting assay shows the importance of PHF8 for endothelial cells. Immunoprecipitation reveals that PHF8 interacts with U1 spliceosomal proteins, such as SRPK1 and snRNP70. We identify the histocompatibility antigen HLA-G as a target of PHF8. The inclusion of HLA-G intron 4, with concomitant RNA Polymerase II accumulation at this intron is controlled by PHF8 and H3K9. Soluble HLA-G is generated after PHF8 knockdown, which leads to reduced T-cell proliferation. Collectively, PHF8 knockdown generates the immunosuppressive alternative splice product soluble HLA-G, which is secreted by endothelial cells to elicit a potential inhibitory effect on inflammation.

WIlsON: Web-based Interactive Omics VisualizatioN.

MOTIVATION: High throughput (HT) screens in the omics field are typically analyzed by automated pipelines that generate static visualizations and comprehensive spreadsheet data for scientists. However, exploratory and hypothesis driven data analysis are key aspects of the understanding of biological systems, both generating extensive need for customized and dynamic visualization. RESULTS: Here we describe WIlsON, an interactive workbench for analysis and visualization of multi-omics data. It is primarily intended to empower screening platforms to offer access to pre-calculated HT screen results to the non-computational scientist. Facilitated by an open file format, WIlsON supports all types of omics screens, serves results via a web-based dashboard, and enables end users to perform analyses and generate publication-ready plots. AVAILABILITY AND IMPLEMENTATION: We implemented WIlsON in R with a focus on extensibility using the modular Shiny and Plotly frameworks. A demo of the interactive workbench without limitations may be accessed at http://loosolab.mpi-bn.mpg.de. A standalone Docker container as well as the source code of WIlsON are freely available from our Docker hub https://hub.docker. com/r/loosolab/wilson, CRAN https://cran.r-project.org/web/packages/wilson/, and GitHub repository https://github.molgen.mpg.de/loosolab/wilson-apps, respectively.

Oncogenic Amplification of Zygotic Dux Factors in Regenerating p53-Deficient Muscle Stem Cells Defines a Molecular Cancer Subtype.

The identity of tumor-initiating cells in many cancer types is unknown. Tumors often express genes associated with embryonic development, although the contributions of zygotic programs to tumor initiation and formation are poorly understood. Here, we show that regeneration-induced loss of quiescence in p53-deficient muscle stem cells (MuSCs) results in rhabdomyosarcoma formation with 100% penetrance. Genomic analyses of purified tumor cells revealed spontaneous and discrete oncogenic amplifications in MuSCs that drive tumorigenesis, including, but not limited to, the amplification of the cleavage-stage Dux transcription factor (TF) Duxbl. We further found that Dux factors drive an early embryonic gene signature that defines a molecular subtype across a broad range of human cancers. Duxbl initiates tumorigenesis by enforcing a mesenchymal-to-epithelial transition, and targeted inactivation of Duxbl specifically in Duxbl-expressing tumor cells abolishes their expansion. These findings reveal how regeneration and genomic instability can interact to activate zygotic genes that drive tumor initiation and growth.

Myh10 deficiency leads to defective extracellular matrix remodeling and pulmonary disease.

Impaired alveolar formation and maintenance are features of many pulmonary diseases that are associated with significant morbidity and mortality. In a forward genetic screen for modulators of mouse lung development, we identified the non-muscle myosin II heavy chain gene, Myh10. Myh10 mutant pups exhibit cyanosis and respiratory distress, and die shortly after birth from differentiation defects in alveolar epithelium and mesenchyme. From omics analyses and follow up studies, we find decreased Thrombospondin expression accompanied with increased matrix metalloproteinase activity in both mutant lungs and cultured mutant fibroblasts, as well as disrupted extracellular matrix (ECM) remodeling. Loss of Myh10 specifically in mesenchymal cells results in ECM deposition defects and alveolar simplification. Notably, MYH10 expression is downregulated in the lung of emphysema patients. Altogether, our findings reveal critical roles for Myh10 in alveologenesis at least in part via the regulation of ECM remodeling, which may contribute to the pathogenesis of emphysema.

Single cell RNA-seq and ATAC-seq analysis of cardiac progenitor cell transition states and lineage settlement.

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.

The potassium channel KCNJ13 is essential for smooth muscle cytoskeletal organization during mouse tracheal tubulogenesis.

Tubulogenesis is essential for the formation and function of internal organs. One such organ is the trachea, which allows gas exchange between the external environment and the lungs. However, the cellular and molecular mechanisms underlying tracheal tube development remain poorly understood. Here, we show that the potassium channel KCNJ13 is a critical modulator of tracheal tubulogenesis. We identify Kcnj13 in an ethylnitrosourea forward genetic screen for regulators of mouse respiratory organ development. Kcnj13 mutants exhibit a shorter trachea as well as defective smooth muscle (SM) cell alignment and polarity. KCNJ13 is essential to maintain ion homeostasis in tracheal SM cells, which is required for actin polymerization. This process appears to be mediated, at least in part, through activation of the actin regulator AKT, as pharmacological increase of AKT phosphorylation ameliorates the Kcnj13-mutant trachea phenotypes. These results provide insight into the role of ion homeostasis in cytoskeletal organization during tubulogenesis.

20-HETE promotes glucose-stimulated insulin secretion in an autocrine manner through FFAR1.

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.

Pulmonary endothelial cell DNA methylation signature in pulmonary arterial hypertension.

Pulmonary arterial hypertension (PAH) is a severe and incurable pulmonary vascular disease. One of the primary origins of PAH is pulmonary endothelial dysfunction leading to vasoconstriction, aberrant angiogenesis and smooth muscle cell proliferation, endothelial-to-mesenchymal transition, thrombosis and inflammation. Our objective was to study the epigenetic variations in pulmonary endothelial cells (PEC) through a specific pattern of DNA methylation. DNA was extracted from cultured PEC from idiopathic PAH (n = 11), heritable PAH (n = 10) and controls (n = 18). DNA methylation was assessed using the Illumina HumanMethylation450 Assay. After normalization, samples and probes were clustered according to their methylation profile. Differential clusters were functionally analyzed using bioinformatics tools. Unsupervised hierarchical clustering allowed the identification of two clusters of probes that discriminates controls and PAH patients. Among 147 differential methylated promoters, 46 promoters coding for proteins or miRNAs were related to lipid metabolism. Top 10 up and down-regulated genes were involved in lipid transport including ABCA1, ABCB4, ADIPOQ, miR-26A, BCL2L11. NextBio meta-analysis suggested a contribution of ABCA1 in PAH. We confirmed ABCA1 mRNA and protein downregulation specifically in PAH PEC by qPCR and immunohistochemistry and made the proof-of-concept in an experimental model of the disease that its targeting may offer novel therapeutic options. In conclusion, DNA methylation analysis identifies a set of genes mainly involved in lipid transport pathway which could be relevant to PAH pathophysiology.

UROPA: a tool for Universal RObust Peak Annotation.

The annotation of genomic ranges of interest represents a recurring task for bioinformatics analyses. These ranges can originate from various sources, including peaks called for transcription factor binding sites (TFBS) or histone modification ChIP-seq experiments, chromatin structure and accessibility experiments (such as ATAC-seq), but also from other types of predictions that result in genomic ranges. While peak annotation primarily driven by ChiP-seq was extensively explored, many approaches remain simplistic ("most closely located TSS"), rely on fixed pre-built references, or require complex scripting tasks on behalf of the user. An adaptable, fast, and universal tool, capable to annotate genomic ranges in the respective biological context is critically missing. UROPA (Universal RObust Peak Annotator) is a command line based tool, intended for universal genomic range annotation. Based on a configuration file, different target features can be prioritized with multiple integrated queries. These can be sensitive for feature type, distance, strand specificity, feature attributes (e.g. protein_coding) or anchor position relative to the feature. UROPA can incorporate reference annotation files (GTF) from different sources (Gencode, Ensembl, RefSeq), as well as custom reference annotation files. Statistics and plots transparently summarize the annotation process. UROPA is implemented in Python and R.

Promoter hypermethylation as a mechanism for Lamin A/C silencing in a subset of neuroblastoma cells.

Nuclear lamins support the nuclear envelope and provide anchorage sites for chromatin. They are involved in DNA synthesis, transcription, and replication. It has previously been reported that the lack of Lamin A/C expression in lymphoma and leukaemia is due to CpG island promoter hypermethylation. Here, we provide evidence that Lamin A/C is silenced via this mechanism in a subset of neuroblastoma cells. Moreover, Lamin A/C expression can be restored with a demethylating agent. Importantly, Lamin A/C reintroduction reduced cell growth kinetics and impaired migration, invasion, and anchorage-independent cell growth. Cytoskeletal restructuring was also induced. In addition, the introduction of lamin Δ50, known as Progerin, caused senescence in these neuroblastoma cells. These cells were stiffer and developed a cytoskeletal structure that differed from that observed upon Lamin A/C introduction. Of relevance, short hairpin RNA Lamin A/C depletion in unmethylated neuroblastoma cells enhanced the aforementioned tumour properties. A cytoskeletal structure similar to that observed in methylated cells was induced. Furthermore, atomic force microscopy revealed that Lamin A/C knockdown decreased cellular stiffness in the lamellar region. Finally, the bioinformatic analysis of a set of methylation arrays of neuroblastoma primary tumours showed that a group of patients (around 3%) gives a methylation signal in some of the CpG sites located within the Lamin A/C promoter region analysed by bisulphite sequencing PCR. These findings highlight the importance of Lamin A/C epigenetic inactivation for a subset of neuroblastomas, leading to enhanced tumour properties and cytoskeletal changes. Additionally, these findings may have treatment implications because tumour cells lacking Lamin A/C exhibit more aggressive behaviour.