Search

Abstract Background Several tools have been developed to explore and search Gene Ontology (GO) databases allowing efficient GO enrichment analysis and GO tree visualization. Nevertheless, identification of highly specific GO-terms in complex data sets is relatively complicated and the display of GO term assignments and GO enrichment analysis by simple tables or pie charts is not optimal. Valuable information such as the hierarchical position of a single GO term within the GO tree (topological ordering), or enrichment within a complex set of biological experiments is not displayed. Pie charts based on GO tree levels are, themselves, one-dimensional graphs, which cannot properly or efficiently represent the hierarchical specificity for the biological system being studied. Results Here we present a new method, which we name PCA2GO, capable of GO analysis using complex multidimensional experimental settings. We employed principal component analysis (PCA) and developed a new score, which takes into account the relative frequency of certain GO terms and their specificity (hierarchical position) within the GO graph. We evaluated the correlation between our representation score R and a standard measure of enrichment, namely p-values to convey the versatility of our approach to other methods and point out differences between our method and commonly used enrichment analyses. Although p values and the R score formally measure different quantities they should be correlated, because relative frequencies of GO terms occurrences within a dataset are an indirect measure of protein numbers related to this term. Therefore they are also related to enrichment. We showed that our score enables us to identify more specific GO-terms i.e. those positioned further down the GO-graph than other common tools used for this purpose. PCA2GO allows visualization and detection of multidimensional dependencies both within the acyclic graph (GO tree) and the experimental settings. Our method is intended for the analysis of several experimental sets, not for one set, like standard enrichment tools. To demonstrate the usefulness of our approach we performed a PCA2GO analysis of a fractionated cardiomyocyte protein dataset, which was identified by enhanced liquid chromatography-mass spectrometry (GeLC-MS). The analysis enabled us to detect distinct groups of proteins, which accurately reflect properties of biochemical cell fractions. Conclusions We conclude that PCA2GO is an alternative efficient GO analysis tool with unique features for detection and visualization of multidimensional dependencies within the dataset under study. PCA2GO reveals strongly correlated GO terms within the experimental setting (in this case different fractions) by PCA group formation and improves detection of more specific GO terms within experiment dependent GO term groups than standard p value calculations.

Abstract Background The newt Notophthalmus viridescens possesses the remarkable ability to respond to cardiac damage by formation of new myocardial tissue. Surprisingly little is known about changes in gene activities that occur during the course of regeneration. To begin to decipher the molecular processes, that underlie restoration of functional cardiac tissue, we generated an EST database from regenerating newt hearts and compared the transcriptional profile of selected candidates with genes deregulated during zebrafish heart regeneration. Results A cDNA library of 100,000 cDNA clones was generated from newt hearts 14 days after ventricular injury. Sequencing of 11520 cDNA clones resulted in 2894 assembled contigs. BLAST searches revealed 1695 sequences with potential homology to sequences from the NCBI database. BLAST searches to TrEMBL and Swiss-Prot databases assigned 1116 proteins to Gene Ontology terms. We also identified a relatively large set of 174 ORFs, which are likely to be unique for urodele amphibians. Expression analysis of newt-zebrafish homologues confirmed the deregulation of selected genes during heart regeneration. Sequences, BLAST results and GO annotations were visualized in a relational web based database followed by grouping of identified proteins into clusters of GO Terms. Comparison of data from regenerating zebrafish hearts identified biological processes, which were uniformly overrepresented during cardiac regeneration in newt and zebrafish. Conclusion We concluded that heart regeneration in newts and zebrafish led to the activation of similar sets of genes, which suggests that heart regeneration in both species might follow similar principles. The design of the newly established newt EST database allows identification of molecular pathways important for heart regeneration.

Developmental transcription factors act in networks, but how these networks achieve cell- and tissue specificity is still poorly understood. Here, we explored pre-B cell leukemia homeobox 1 (PBX1) in adult neurogenesis combining genomic, transcriptomic, and proteomic approaches. ChIP-seq analysis uncovered PBX1 binding to numerous genomic sites. Integration of PBX1 ChIP-seq with ATAC-seq data predicted interaction partners, which were subsequently validated by mass spectrometry. Whole transcriptome spatial RNA analysis revealed shared expression dynamics of Pbx1 and interacting factors. Among these were class I bHLH proteins TCF3 and TCF4. RNA-seq following Pbx1, Tcf3 or Tcf4 knockdown identified proliferation- and differentiation associated genes as shared targets, while sphere formation assays following knockdown argued for functional cooperativity of PBX1 and TCF3 in progenitor cell proliferation. Notably, while physiological PBX1-TCF interaction has not yet been described, chromosomal translocation resulting in genomic TCF3::PBX1 fusion characterizes a subtype of acute lymphoblastic leukemia. Introducing Pbx1 into Nalm6 cells, a pre-B cell line expressing TCF3 but lacking PBX1, upregulated the leukemogenic genes BLK and NOTCH3, arguing that functional PBX1-TCF cooperativity likely extends to hematopoiesis. Our study hence uncovers a transcriptional module orchestrating the balance between progenitor cell proliferation and differentiation in adult neurogenesis with potential implications for leukemia etiology., (© The Author(s) 2024. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Postnatal mammalian cardiomyocytes (CMs) rapidly lose proliferative capacity and exit the cell cycle and undergo further differentiation and maturation. Cell cycle activation has been a major strategy to stimulate postnatal CM proliferation, albeit achieving modest effects. One impediment is that postnatal CMs may need to undergo dedifferentiation before proliferation, if not simultaneously. Here, we report that overexpression of Hdac7 in neonatal mouse CMs results in significant CM dedifferentiation and proliferation. Mechanistically, we show that HDAC7-mediated CM proliferation is contingent on dedifferentiation, which is accomplished through suppressing MEF2. Hdac7 overexpression in CM shifts the chromatin state from binding MEF2, which favors the differentiation transcriptional program to AP-1, which favors the proliferative transcriptional program. Further, we found that HDAC7 interacts with minichromosome maintenance complex (MCM) components to initiate cell cycle progression. Our findings reveal that HDAC7 promotes CM proliferation by its dual action on CM dedifferentiation and proliferation, uncovering a potential new strategy for heart regeneration/repair., (© The Author(s) 2024. Published by Oxford University Press on behalf of Institute of Biochemistry and Cell Biology, Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences.)

In higher organisms, individual cells respond to signals and perturbations by epigenetic regulation and transcriptional adaptation. However, in addition to shifting the expression level of individual genes, the adaptive response of cells can also lead to shifts in the proportions of different cell types. Recent methods such as scRNA-seq allow for the interrogation of expression on the single-cell level, and can quantify individual cell type clusters within complex tissue samples. In order to identify clusters showing differential composition between different biological conditions, differential proportion analysis has recently been introduced. However, bioinformatics tools for robust proportion analysis of both replicated and unreplicated single-cell datasets are critically missing. In this manuscript, we present Scanpro, a modular tool for proportion analysis, seamlessly integrating into widely accepted frameworks in the Python environment. Scanpro is fast, accurate, supports datasets without replicates, and is intended to be used by bioinformatics experts and beginners alike., (© 2024. The Author(s).)

Biomechanical forces, and their molecular transducers, including key mechanosensitive transcription factor genes, such as KLF2 , are required for cardiac valve morphogenesis. However, klf2 mutants fail to completely recapitulate the valveless phenotype observed under no-flow conditions. Here, we identify the transcription factor EGR3 as a conserved biomechanical force transducer critical for cardiac valve formation. We first show that egr3 null zebrafish display a complete and highly penetrant loss of valve leaflets, leading to severe blood regurgitation. Using tissue-specific loss- and gain-of-function tools, we find that during cardiac valve formation, Egr3 functions cell-autonomously in endothelial cells, and identify one of its effectors, the nuclear receptor Nr4a2b. We further find that mechanical forces up-regulate egr3 / EGR3 expression in the developing zebrafish heart and in porcine valvular endothelial cells, as well as during human aortic valve remodeling. Altogether, these findings reveal that EGR3 is necessary to transduce the biomechanical cues required for zebrafish cardiac valve morphogenesis, and potentially for pathological aortic valve remodeling in humans.

Failure of proper ventricular trabeculation is often associated with congenital heart disease. Support from endocardial cells, including the secretion of extracellular matrix and growth factors is critical for trabeculation. However, it is poorly understood how the secretion of extracellular matrix and growth factors is initiated and regulated by endocardial cells. We find that genetic knockout of histone deacetylase 3 in the endocardium in mice results in early embryo lethality and ventricular hypotrabeculation. Single cell RNA sequencing identifies significant downregulation of extracellular matrix components in histone deacetylase 3 knockout endocardial cells. Secretome from cultured histone deacetylase 3 knockout mouse cardiac endothelial cells lacks transforming growth factor ß3 and shows significantly reduced capacity in stimulating cultured cardiomyocyte proliferation, which is remarkably rescued by transforming growth factor ß3 supplementation. Mechanistically, we identify that histone deacetylase 3 knockout induces transforming growth factor ß3 expression through repressing microRNA-129-5p. Our findings provide insights into the pathogenesis of congenital heart disease and conceptual strategies to promote myocardial regeneration., (© 2024. The Author(s).)

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., (© 2024. The Author(s).)

Transcription factors (TFs) are crucial epigenetic regulators, which enable cells to dynamically adjust gene expression in response to environmental signals. Computational procedures like digital genomic footprinting on chromatin accessibility assays such as ATACseq can be used to identify bound TFs in a genome-wide scale. This method utilizes short regions of low accessibility signals due to steric hindrance of DNA bound proteins, called footprints (FPs), which are combined with motif databases for TF identification. However, while over 1600 TFs have been described in the human genome, only ~ 700 of these have a known binding motif. Thus, a substantial number of FPs without overlap to a known DNA motif are normally discarded from FP analysis. In addition, the FP method is restricted to organisms with a substantial number of known TF motifs. Here we present DENIS (DE Novo motIf diScovery), a framework to generate and systematically investigate the potential of de novo TF motif discovery from FPs. DENIS includes functionality (1) to isolate FPs without binding motifs, (2) to perform de novo motif generation and (3) to characterize novel motifs. Here, we show that the framework rediscovers artificially removed TF motifs, quantifies de novo motif usage during an early embryonic development example dataset, and is able to analyze and uncover TF activity in organisms lacking canonical motifs. The latter task is exemplified by an investigation of a scATAC-seq dataset in zebrafish which covers different cell types during hematopoiesis., (© 2024. The Author(s).)

Cardiac contractions and hemodynamic forces are essential for organ development and homeostasis. Control over cardiac contractions can be achieved pharmacologically or optogenetically. However, these approaches lack specificity or require direct access to the heart. Here, we compare two genetic approaches to control cardiac contractions by modulating the levels of the essential sarcomeric protein Tnnt2a in zebrafish. We first recombine a newly generated tnnt2a floxed allele using multiple lines expressing Cre under the control of cardiomyocyte-specific promoters, and show that it does not recapitulate the tnnt2a/silent heart mutant phenotype in embryos. We show that this lack of early cardiac contraction defects is due, at least in part, to the long half-life of tnnt2a mRNA, which masks the gene deletion effects until the early larval stages. We then generate an endogenous Tnnt2a-eGFP fusion line that we use together with the zGRAD system to efficiently degrade Tnnt2a in all cardiomyocytes. Using single-cell transcriptomics, we find that Tnnt2a depletion leads to cardiac phenotypes similar to those observed in tnnt2a mutants, with a loss of blood and pericardial flow-dependent cell types. Furthermore, we achieve conditional degradation of Tnnt2a-eGFP by splitting the zGRAD protein into two fragments that, when combined with the cpFRB2-FKBP system, can be reassembled upon rapamycin treatment. Thus, this Tnnt2a degradation line enables non-invasive control of cardiac contractions with high spatial and temporal specificity and will help further understand how they shape organ development and homeostasis., Competing Interests: Competing interests statement:The authors declare no competing interest.

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., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2024 Just et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Background: Advanced age is unequivocally linked to the development of cardiovascular disease; however, the mechanisms resulting in reduced endothelial cell regeneration remain poorly understood. Here, we investigated novel mechanisms involved in endothelial cell senescence that impact endothelial cell transcription and vascular repair after injury., Methods: Native endothelial cells were isolated from young (20±3.4 years) and aged (80±2.3 years) individuals and subjected to molecular analyses to assess global transcriptional and metabolic changes. In vitro studies were conducted using primary human and murine endothelial cells. A murine aortic re-endothelialization model was used to examine endothelial cell regenerative capacity in vivo., Results: RNA sequencing of native endothelial cells revealed that aging resulted in p53-mediated reprogramming to express senescence-associated genes and suppress glycolysis. Reduced glucose uptake and ATP contributed to attenuated assembly of the telomerase complex, which was required for endothelial cell proliferation. Enhanced p53 activity in aging was linked to its acetylation on K120 due to enhanced activity of the acetyltransferase MOZ (monocytic leukemic zinc finger). Mechanistically, p53 acetylation and translocation were, at least partially, attributed to the loss of the vasoprotective enzyme, CSE (cystathionine γ-lyase). CSE physically anchored p53 in the cytosol to prevent its nuclear translocation and CSE absence inhibited AKT (Protein kinase B)-mediated MOZ phosphorylation, which in turn increased MOZ activity and subsequently p53 acetylation. In mice, the endothelial cell-specific deletion of CSE activated p53, induced premature endothelial senescence, and arrested vascular repair after injury. In contrast, the adeno-associated virus 9-mediated re-expression of an active CSE mutant retained p53 in the cytosol, maintained endothelial glucose metabolism and proliferation, and prevented endothelial cell senescence. Adenoviral overexpression of CSE in native endothelial cells from aged individuals maintained low p53 activity and reactivated telomerase to revert endothelial cell senescence., Conclusions: Aging-associated impairment of vascular repair is partly determined by the vasoprotective enzyme CSE., Competing Interests: Disclosures None.

Postnatal maturation of cardiomyocytes is characterized by a metabolic switch from glycolysis to fatty acid oxidation, chromatin reconfiguration and exit from the cell cycle, instating a barrier for adult heart regeneration 1,2 . Here, to explore whether metabolic reprogramming can overcome this barrier and enable heart regeneration, we abrogate fatty acid oxidation in cardiomyocytes by inactivation of Cpt1b. We find that disablement of fatty acid oxidation in cardiomyocytes improves resistance to hypoxia and stimulates cardiomyocyte proliferation, allowing heart regeneration after ischaemia-reperfusion injury. Metabolic studies reveal profound changes in energy metabolism and accumulation of α-ketoglutarate in Cpt1b-mutant cardiomyocytes, leading to activation of the α-ketoglutarate-dependent lysine demethylase KDM5 (ref.  3 ). Activated KDM5 demethylates broad H3K4me3 domains in genes that drive cardiomyocyte maturation, lowering their transcription levels and shifting cardiomyocytes into a less mature state, thereby promoting proliferation. We conclude that metabolic maturation shapes the epigenetic landscape of cardiomyocytes, creating a roadblock for further cell divisions. Reversal of this process allows repair of damaged hearts., (© 2023. The Author(s).)

Background: The ability of the right ventricle (RV) to adapt to an increased pressure afterload determines survival in patients with pulmonary arterial hypertension. At present, there are no specific treatments available to prevent RV failure, except for heart/lung transplantation. The wingless/int-1 (Wnt) signaling pathway plays an important role in the development of the RV and may also be implicated in adult cardiac remodeling., Methods: Molecular, biochemical, and pharmacological approaches were used both in vitro and in vivo to investigate the role of Wnt signaling in RV remodeling., Results: Wnt/β-catenin signaling molecules are upregulated in RV of patients with pulmonary arterial hypertension and animal models of RV overload (pulmonary artery banding-induced and monocrotaline rat models). Activation of Wnt/β-catenin signaling leads to RV remodeling via transcriptional activation of FOSL1 and FOSL2 (FOS proto-oncogene [FOS] like 1/2, AP-1 [activator protein 1] transcription factor subunit). Immunohistochemical analysis of pulmonary artery banding -exposed BAT-Gal (β-catenin-activated transgene driving expression of nuclear β-galactosidase) reporter mice RVs exhibited an increase in β-catenin expression compared with their respective controls. Genetic inhibition of β-catenin, FOSL1/2, or WNT3A stimulation of RV fibroblasts significantly reduced collagen synthesis and other remodeling genes. Importantly, pharmacological inhibition of Wnt signaling using inhibitor of PORCN (porcupine O-acyltransferase), LGKK-974 attenuated fibrosis and cardiac hypertrophy leading to improvement in RV function in both, pulmonary artery banding - and monocrotaline-induced RV overload., Conclusions: Wnt- β-Catenin-FOSL signaling is centrally involved in the hypertrophic RV response to increased afterload, offering novel targets for therapeutic interference with RV failure in pulmonary hypertension., Competing Interests: Disclosures None.

Background & Aims: Epigenetic processes regulating gene expression contribute markedly to epithelial cell plasticity in colorectal carcinogenesis. The lysine methyltransferase SUV420H2 comprises an important regulator of epithelial plasticity and is primarily responsible for trimethylation of H4K20 (H4K20me3). Loss of H4K20me3 has been suggested as a hallmark of human cancer due to its interaction with DNMT1. However, the role of Suv4-20h2 in colorectal cancer is unknown., Methods: We examined the alterations in histone modifications in patient-derived colorectal cancer organoids. Patient-derived colorectal cancer organoids and mouse intestinal organoids were genetically manipulated for functional studies in patient-derived xenograft and orthotopic transplantation. Gene expression profiling, micrococcal nuclease assay, and chromatin immunoprecipitation were performed to understand epigenetic regulation of chromatin states and gene expression in patient-derived and mouse intestinal organoids., Results: We found that reduced H4K20me3 levels occurred predominantly in right-sided patient-derived colorectal cancer organoids, which were associated with increased chromatin accessibility. Re-compaction of chromatin by methylstat, a histone demethylase inhibitor, resulted in reduced growth selectively in subcutaneously grown tumors derived from right-sided cancers. Using mouse intestinal organoids, we confirmed that Suv4-20h2-mediated H4K20me3 is required for maintaining heterochromatin compaction and to prevent R-loop formation. Cross-species comparison of Suv4-20h2-depleted murine organoids with right-sided colorectal cancer organoids revealed a large overlap of gene signatures involved in chromatin silencing, DNA methylation, and stemness/Wnt signaling., Conclusions: Loss of Suv4-20h2-mediated H4K20me3 drives right-sided colorectal tumorigenesis through an epigenetically controlled mechanism of chromatin compaction. Our findings unravel a conceptually novel approach for subtype-specific therapy of this aggressive form of colorectal cancer., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

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., (© 2022. The Author(s).)

The incidence and prevalence of cardiovascular disease is highest among the elderly. There is a need to further understand the mechanisms behind endothelial cell aging in order to achieve vascular rejuvenation and minimize the onset of age-related vascular diseases. Long non-coding RNAs (lncRNAs) have been proposed to regulate numerous processes in the human genome, yet their function in vascular aging and their therapeutic potential remain largely unknown. This is primarily because the majority of studies investigating the impact of aging on lncRNA expression heavily rely on in vitro studies based on replicative senescence. Here, using a unique collection of young and aged endothelial cells isolated from native human arteries, we sought to characterize the age-related alterations in lncRNA expression profiles. We were able to detect a total of 4463 lncRNAs expressed in the human endothelium from which ∼17% (798) were altered in advanced age. One of the most affected lncRNAs in aging was the primate-specific, Prostate Cancer Associated Transcript (PCAT) 14. In our follow up analysis, using single molecule RNA FISH, we showed that PCAT14 is relatively abundant, localized almost exclusively in the nucleus of young endothelial cells, and silenced in the aged endothelium. Functionally, our studies proposed that downregulation of PCAT14 alters endothelial cell transcription profile and cell functions including endothelial cell migration, sprouting and inflammatory responses in vitro . Taken together, our data highlight that endothelial cell aging correlates with altered expression of lncRNAs, which could impair the endothelial regenerative capacity and enhance inflammatory phenotypes., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Drekolia, Talyan, Cordellini Emídio, Boon, Guenther, Looso, Dumbović and Bibli.)

Cooperativity between transcription factors is important to regulate target gene expression. In particular, the binding grammar of TFs in relation to each other, as well as in the context of other genomic elements, is crucial for TF functionality. However, tools to easily uncover co-occurrence between DNA-binding proteins, and investigate the regulatory modules of TFs, are limited. Here we present TF-COMB (Transcription Factor Co-Occurrence using Market Basket analysis) - a tool to investigate co-occurring TFs and binding grammar within regulatory regions. We found that TF-COMB can accurately identify known co-occurring TFs from ChIP-seq data, as well as uncover preferential localization to other genomic elements. With the use of ATAC-seq footprinting and TF motif locations, we found that TFs exhibit both preferred orientation and distance in relation to each other, and that these are biologically significant. Finally, we extended the analysis to not only investigate individual TF pairs, but also TF pairs in the context of networks, which enabled the investigation of TF complexes and TF hubs. In conclusion, TF-COMB is a flexible tool to investigate various aspects of TF binding grammar., Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (© 2022 The Author(s).)

Chronic obstructive pulmonary disease (COPD) is a disease with an inflammatory phenotype with increasing prevalence in the elderly. Expanded population of mutant blood cells carrying somatic mutations is termed clonal hematopoiesis of indeterminate potential (CHIP). The association between CHIP and COPD and its relevant effects on DNA methylation in aging are mainly unknown. Analyzing the deep-targeted amplicon sequencing from 125 COPD patients, we found enhanced incidence of CHIP mutations (~20%) with a predominance of DNMT3A CHIP-mediated hypomethylation of Phospholipase D Family Member 5 ( PLD5 ), which in turn is positively correlated with increased levels of glycerol phosphocholine, pro-inflammatory cytokines, and deteriorating lung function.

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)., (Copyright © 2022 Beijing Institute of Genomics. Published by Elsevier B.V. All rights reserved.)

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.

Background: In vascular endothelial cells, cysteine metabolism by the cystathionine γ lyase (CSE), generates hydrogen sulfide-related sulfane sulfur compounds (H 2 S n ), that exert their biological actions via cysteine S -sulfhydration of target proteins. This study set out to map the " S -sulfhydrome" (ie, the spectrum of proteins targeted by H 2 S n ) in human endothelial cells., Methods: Liquid chromatography with tandem mass spectrometry was used to identify S -sulfhydrated cysteines in endothelial cell proteins and β3 integrin intraprotein disulfide bond rearrangement. Functional studies included endothelial cell adhesion, shear stress-induced cell alignment, blood pressure measurements, and flow-induced vasodilatation in endothelial cell-specific CSE knockout mice and in a small collective of patients with endothelial dysfunction., Results: Three paired sample sets were compared: (1) native human endothelial cells isolated from plaque-free mesenteric arteries (CSE activity high) and plaque-containing carotid arteries (CSE activity low); (2) cultured human endothelial cells kept under static conditions or exposed to fluid shear stress to decrease CSE expression; and (3) cultured endothelial cells exposed to shear stress to decrease CSE expression and treated with solvent or the slow-releasing H 2 S n donor, SG1002. The endothelial cell " S -sulfhydrome" consisted of 3446 individual cysteine residues in 1591 proteins. The most altered family of proteins were the integrins and focusing on β3 integrin in detail we found that S -sulfhydration affected intraprotein disulfide bond formation and was required for the maintenance of an extended-open conformation of the β leg. β3 integrin S -sulfhydration was required for endothelial cell mechanotransduction in vitro as well as flow-induced dilatation in murine mesenteric arteries. In cultured cells, the loss of S -sulfhydration impaired interactions between β3 integrin and Gα13 (guanine nucleotide-binding protein subunit α 13), resulting in the constitutive activation of RhoA (ras homolog family member A) and impaired flow-induced endothelial cell realignment. In humans with atherosclerosis, endothelial function correlated with low H 2 S n generation, impaired flow-induced dilatation, and failure to detect β3 integrin S -sulfhydration, all of which were rescued after the administration of an H 2 S n supplement., Conclusions: Vascular disease is associated with marked changes in the S -sulfhydration of endothelial cell proteins involved in mediating responses to flow. Short-term H 2 S n supplementation improved vascular reactivity in humans highlighting the potential of interfering with this pathway to treat vascular disease.

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.

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., (© 2020 Bhagwat et al.)

Communication between individuals via molecules, termed chemosignaling, is widespread among animal and plant species. However, we lack knowledge on the specific functions of the substances involved for most systems. The femoral gland is an organ that secretes a waxy substance involved in chemical communication in lizards. Although the lipids and volatile substances secreted by the femoral glands have been investigated in several biochemical studies, the protein composition and functions of secretions remain completely unknown. Applying a proteomic approach, we provide the first attempt to comprehensively characterize the protein composition of femoral gland secretions from the Galápagos marine iguana. Using samples from several organs, the marine iguana proteome was assembled by next-generation sequencing and MS, resulting in 7513 proteins. Of these, 4305 proteins were present in the femoral gland, including keratins, small serum proteins, and fatty acid-binding proteins. Surprisingly, no proteins with discernible roles in partner recognition or inter-species communication could be identified. However, we did find several proteins with direct associations to the innate immune system, including lysozyme C, antileukoproteinase (ALP), pulmonary surfactant protein (SFTPD), and galectin (LGALS1) suggesting that the femoral glands function as an important barrier to infection. Furthermore, we report several novel anti-microbial peptides from the femoral glands that show similar action against Escherichia coli and Bacillus subtilis such as oncocin, a peptide known for its effectiveness against Gram-negative pathogens. This proteomics data set is a valuable resource for future functional protein analysis and demonstrates that femoral gland secretions also perform functions of the innate immune system., Competing Interests: Conflict of interest—Authors declare no competing interests., (© 2020 Tellkamp et al.)

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.

Rationale: The adult human heart is an organ with low regenerative potential. Heart failure following acute myocardial infarction is a leading cause of death due to the inability of cardiomyocytes to proliferate and replenish lost cardiac muscle. While the zebrafish has emerged as a powerful model to study endogenous cardiac regeneration, the molecular mechanisms by which cardiomyocytes respond to damage by disassembling sarcomeres, proliferating, and repopulating the injured area remain unclear. Furthermore, we are far from understanding the regulation of the chromatin landscape and epigenetic barriers that must be overcome for cardiac regeneration to occur., Objective: To identify transcription factor regulators of the chromatin landscape, which promote cardiomyocyte regeneration in zebrafish, and investigate their function., Methods and Results: Using the Assay for Transposase-Accessible Chromatin coupled to high-throughput sequencing (ATAC-Seq), we first find that the regenerating cardiomyocyte chromatin accessibility landscape undergoes extensive changes following cryoinjury, and that activator protein-1 (AP-1) binding sites are the most highly enriched motifs in regions that gain accessibility during cardiac regeneration. Furthermore, using bioinformatic and gene expression analyses, we find that the AP-1 response in regenerating adult zebrafish cardiomyocytes is largely different from the response in adult mammalian cardiomyocytes. Using a cardiomyocyte-specific dominant negative approach, we show that blocking AP-1 function leads to defects in cardiomyocyte proliferation as well as decreased chromatin accessibility at the fbxl22 and ilk loci, which regulate sarcomere disassembly and cardiomyocyte protrusion into the injured area, respectively. We further show that overexpression of the AP-1 family members Junb and Fosl1 can promote changes in mammalian cardiomyocyte behavior in vitro., Conclusions: AP-1 transcription factors play an essential role in the cardiomyocyte response to injury by regulating chromatin accessibility changes, thereby allowing the activation of gene expression programs that promote cardiomyocyte dedifferentiation, proliferation, and protrusion into the injured area.

Aims: To assess the functional relevance and therapeutic potential of the pro-angiogenic long non-coding RNA MANTIS in vascular disease development., Methods and Results: RNA sequencing, CRISPR activation, overexpression, and RNAi demonstrated that MANTIS, especially its Alu-element, limits endothelial ICAM-1 expression in different types of endothelial cells. Loss of MANTIS increased endothelial monocyte adhesion in an ICAM-1-dependent manner. MANTIS reduced the binding of the SWI/SNF chromatin remodelling factor BRG1 at the ICAM-1 promoter. The expression of MANTIS was induced by laminar flow and HMG-CoA-reductase inhibitors (statins) through mechanisms involving epigenetic rearrangements and the transcription factors KLF2 and KLF4. Mutation of the KLF binding motifs in the MANTIS promoter blocked the flow-induced MANTIS expression. Importantly, the expression of MANTIS in human carotid artery endarterectomy material was lower compared with healthy vessels and this effect was prevented by statin therapy. Interestingly, the protective effects of statins were mediated in part through MANTIS, which was required to facilitate the atorvastatin-induced changes in endothelial gene expression. Moreover, the beneficial endothelial effects of statins in culture models (spheroid outgrowth, proliferation, telomerase activity, and vascular organ culture) were lost upon knockdown of MANTIS., Conclusion: MANTIS is tightly regulated by the transcription factors KLF2 and KLF4 and limits the ICAM-1 mediated monocyte adhesion to endothelial cells and thus potentially atherosclerosis development in humans. The beneficial effects of statin treatment and laminar flow are dependent on MANTIS., (Published on behalf of the European Society of Cardiology. All rights reserved. © The Author(s) 2019. For permissions, please email: journals.permissions@oup.com.)

In the context of publishing data sets acquired by mass spectrometry or works based on such molecular screens, metadata documenting the instrument settings are of central importance to the evaluation and reproduction of results. A single experiment may be linked to hundreds of data acquisitions, which are frequently stored in proprietary file formats. Together with community-, repository-, as well as publisher-specific reporting standards, this state of affairs frequently leads to manual -and thus error prone-metadata extraction and formatting. Data extracted from a single file also often stand in for an entire file set, implying a risk for unreported parameter divergence. To support quality control and data reporting, the C# application MARMoSET extracts and reduces publication relevant metadata from Thermo Fischer Scientific RAW files. It is integrated with an R package for easy reporting. The tool is expected to be particularly useful to high throughput environments such as service facilities with large project numbers and/or sizes., (© 2019 Kiweler et al.)

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., Competing Interests: Competing interestsThe authors declare no competing or financial interests., (© 2019. Published by The Company of Biologists Ltd.)

Background: The potential of next-generation sequencing (NGS) for hypothesis-free pathogen diagnosis from (poly-)microbially contaminated, formalin-fixed, paraffin embedded tissue samples from patients with invasive fungal infections and amebiasis was investigated. Samples from patients with chromoblastomycosis (n = 3), coccidioidomycosis (n = 2), histoplasmosis (n = 4), histoplasmosis or cryptococcosis with poor histological discriminability (n = 1), mucormycosis (n = 2), mycetoma (n = 3), rhinosporidiosis (n = 2), and invasive Entamoeba histolytica infections (n = 6) were analyzed by NGS (each one Illumina v3 run per sample). To discriminate contamination from putative infections in NGS analysis, mean and standard deviation of the number of specific sequence fragments (paired reads) were determined and compared in all samples examined for the pathogens in question., Results: For matches between NGS results and histological diagnoses, a percentage of species-specific reads greater than the 4th standard deviation above the mean value of all 23 assessed sample materials was required. Potentially etiologically relevant pathogens could be identified by NGS in 5 out of 17 samples of patients with invasive mycoses and in 1 out of 6 samples of patients with amebiasis., Conclusions: The use of NGS for hypothesis-free pathogen diagnosis from contamination-prone formalin-fixed, paraffin-embedded tissue requires further standardization.

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., (© The Author(s) 2018. Published by Oxford University Press.)

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., (© 2019 Federation of European Biochemical Societies.)

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.

Strict regulation of proliferation is vital for development, whereas unregulated cell proliferation is a fundamental characteristic of cancer. The polarity protein atypical protein kinase C lambda/iota (aPKCλ) is associated with cell proliferation through unknown mechanisms. In endothelial cells, suppression of aPKCλ impairs proliferation despite hyperactivated mitogenic signaling. Here we show that aPKCλ phosphorylates the DNA binding domain of forkhead box O1 (FoxO1) transcription factor, a gatekeeper of endothelial growth. Although mitogenic signaling excludes FoxO1 from the nucleus, consequently increasing c-Myc abundance and proliferation, aPKCλ controls c-Myc expression via FoxO1/miR-34c signaling without affecting its localization. We find this pathway is strongly activated in the malignant vascular sarcoma, angiosarcoma, and aPKC inhibition reduces c-Myc expression and proliferation of angiosarcoma cells. Moreover, FoxO1 phosphorylation at Ser218 and aPKC expression correlates with poor patient prognosis. Our findings may provide a potential therapeutic strategy for treatment of malignant cancers, like angiosarcoma.

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., (Copyright © 2018 Elsevier Inc. All rights reserved.)