RegenDbase: a comparative database of noncoding RNA regulation of tissue regeneration circuits across multiple taxa

Reusable tutorials for using cloud-based computing environments for the analysis of bacterial gene expression data from bulk RNA sequencing

This manuscript describes the development of a resource module that is part of a learning platform named "NIGMS Sandbox for Cloud-based Learning" https://github.com/NIGMS/NIGMS-Sandbox. The overall genesis of the Sandbox is described in the editorial NIGMS Sandbox at the beginning of this Supplement. This module delivers learning materials on RNA sequencing (RNAseq) data analysis in an interactive format that uses appropriate cloud resources for data access and analyses. Biomedical research is increasingly data-driven, and dependent upon data management and analysis methods that facilitate rigorous, robust, and reproducible research. Cloud-based computing resources provide opportunities to broaden the application of bioinformatics and data science in research. Two obstacles for researchers, particularly those at small institutions, are: (i) access to bioinformatics analysis environments tailored to their research; and (ii) training in how to use Cloud-based computing resources. We developed five reusable tutorials for bulk RNAseq data analysis to address these obstacles. Using Jupyter notebooks run on the Google Cloud Platform, the tutorials guide the user through a workflow featuring an RNAseq dataset from a study of prophage altered drug resistance in Mycobacterium chelonae. The first tutorial uses a subset of the data so users can learn analysis steps rapidly, and the second uses the entire dataset. Next, a tutorial demonstrates how to analyze the read count data to generate lists of differentially expressed genes using R/DESeq2. Additional tutorials generate read counts using the Snakemake workflow manager and Nextflow with Google Batch. All tutorials are open-source and can be used as templates for other analysis.

Transcriptomic data meta-analysis reveals common and injury model specific gene expression changes in the regenerating zebrafish heart

Zebrafish have the capacity to fully regenerate the heart after an injury, which lies in sharp contrast to the irreversible loss of cardiomyocytes after a myocardial infarction in humans. Transcriptomics analysis has contributed to dissect underlying signaling pathways and gene regulatory networks in the zebrafish heart regeneration process. This process has been studied in response to different types of injuries namely: ventricular resection, ventricular cryoinjury, and genetic ablation of cardiomyocytes. However, there exists no database to compare injury specific and core cardiac regeneration responses. Here, we present a meta-analysis of transcriptomic data of regenerating zebrafish hearts in response to these three injury models at 7 days post injury (7dpi). We reanalyzed 36 samples and analyzed the differentially expressed genes (DEG) followed by downstream Gene Ontology Biological Processes (GO:BP) analysis. We found that the three injury models share a common core of DEG encompassing genes involved in cell proliferation, the Wnt signaling pathway and genes that are enriched in fibroblasts. We also found injury-specific gene signatures for resection and genetic ablation, and to a lower extent the cryoinjury model. Finally, we present our data in a user-friendly web interface that displays gene expression signatures across different injury types and highlights the importance to consider injury-specific gene regulatory networks when interpreting the results related to cardiac regeneration in the zebrafish. The analysis is freely available at: https://mybinder.org/v2/gh/MercaderLabAnatomy/PUB_Botos_et_al_2022_shinyapp_binder/HEAD?urlpath=shiny/bus-dashboard/ .

Reconstruction of regulatory network predicts transcription factors driving the dynamics of zebrafish heart regeneration

The limited regenerative capacity in mammals has serious implications for cardiac tissue damage. Meanwhile, zebrafish has a high regenerative capacity, but the regulation of the heart healing process has yet to be elucidated. The dynamic nature of cardiac regeneration requires consideration of the inherent temporal dimension of this process. Here, we conducted a systematic review to find genes that define the regenerative cell state of the zebrafish heart. We then performed an in silico temporal gene regulatory network analysis using transcriptomic data from the zebrafish heart regenerative process obtained from databases. In this analysis, the genes found in the systematic review were used to represent the final cell state of the transition process from a non-regenerative cell state to a regenerative state. We found 135 transcription factors driving the cellular state transition process during zebrafish cardiac regeneration, including Hand2, Nkx2.5, Tbx20, Fosl1, Fosb, Junb, Vdr, Wt1, and Tcf21 previously reported for playing a key role in tissue regeneration. Furthermore, we demonstrate that most regulators are activated in the first days post-injury, indicating that the transition from a non-regenerative to a regenerative state occurs promptly.

Modeling Virus-Induced Inflammation in Zebrafish: A Balance Between Infection Control and Excessive Inflammation

The inflammatory response to viral infection in humans is a dynamic process with complex cell interactions that are governed by the immune system and influenced by both host and viral factors. Due to this complexity, the relative contributions of the virus and host factors are best studied in vivo using animal models. In this review, we describe how the zebrafish (Danio rerio) has been used as a powerful model to study host-virus interactions and inflammation by combining robust forward and reverse genetic tools with in vivo imaging of transparent embryos and larvae. The innate immune system has an essential role in the initial inflammatory response to viral infection. Focused studies of the innate immune response to viral infection are possible using the zebrafish model as there is a 4-6 week timeframe during development where they have a functional innate immune system dominated by neutrophils and macrophages. During this timeframe, zebrafish lack a functional adaptive immune system, so it is possible to study the innate immune response in isolation. Sequencing of the zebrafish genome has revealed significant genetic conservation with the human genome, and multiple studies have revealed both functional conservation of genes, including those critical to host cell infection and host cell inflammatory response. In addition to studying several fish viruses, zebrafish infection models have been developed for several human viruses, including influenza A, noroviruses, chikungunya, Zika, dengue, herpes simplex virus type 1, Sindbis, and hepatitis C virus. The development of these diverse viral infection models, coupled with the inherent strengths of the zebrafish model, particularly as it relates to our understanding of macrophage and neutrophil biology, offers opportunities for far more intensive studies aimed at understanding conserved host responses to viral infection. In this context, we review aspects relating to the evolution of innate immunity, including the evolution of viral pattern recognition receptors, interferons and interferon receptors, and non-coding RNAs.

Towards comparative analyses of salamander limb regeneration

Among tetrapods, only salamanders can regenerate their limbs and tails throughout life. This amazing regenerative ability has attracted the attention of scientists for hundreds of years. Now that large, salamander genomes are beginning to be sequenced for the first time, omics tools and approaches can be used to integrate new perspectives into the study of tissue regeneration. Here we argue the need to move beyond the primary salamander models to investigate regeneration in other species. Salamanders at first glance come across as a phylogenetically conservative group that has not diverged greatly from their ancestors. While salamanders do present ancestral characteristics of basal tetrapods, including the ability to regenerate limbs, data from fossils and data from studies that have tested for species differences suggest there may be considerable variation in how salamanders develop and regenerate their limbs. We review the case for expanded studies of salamander tissue regeneration and identify questions and approaches that are most likely to reveal commonalities and differences in regeneration among species. We also address challenges that confront such an initiative, some of which are regulatory and not scientific. The time is right to gain evolutionary perspective about mechanisms of tissue regeneration from comparative studies of salamander species.

Exploring the role of microRNAs in axolotl regeneration

The axolotl, Ambystoma mexicanum, is used extensively for research in developmental biology, particularly for its ability to regenerate and restore lost organs, including in the nervous system, to full functionality. Regeneration in mammals typically depends on the healing process and scar formation with limited replacement of lost tissue. Other organisms, such as spiny mice (Acomys cahirinus), salamanders, and zebrafish, are able to regenerate some damaged body components. Blastema is a tissue that is formed after tissue injury in such organisms and is composed of progenitor cells or dedifferentiated cells that differentiate into various cell types during regeneration. Thus, identifying the molecules responsible for initiation of blastema formation is an important aspect for understanding regeneration. Introns, a major source of noncoding RNAs (ncRNAs), have characteristic sizes in the axolotl, particularly in genes associated with development. These ncRNAs, particularly microRNAs (miRNAs), exhibit dynamic regulation during regeneration. These miRNAs play an essential role in timing and control of gene expression to order and organize processes necessary for blastema creation. Master keys or molecules that underlie the remarkable regenerative abilities of the axolotl remain to be fully explored and exploited. Further and ongoing research on regeneration promises new knowledge that may allow improved repair and renewal of human tissues.

The power and limitations of gene expression pathway analyses toward predicting population response to environmental stressors

Rapid environmental changes impact the global distribution and abundance of species, highlighting the urgency to understand and predict how populations will respond. The analysis of differentially expressed genes has elucidated areas of the genome involved in adaptive divergence to past and present environmental change. Such studies however have been hampered by large numbers of differentially expressed genes and limited knowledge of how these genes work in conjunction with each other. Recent methods (broadly termed "pathway analyses") have emerged that aim to group genes that behave in a coordinated fashion to a factor of interest. These methods aid in functional annotation and uncovering biological pathways, thereby collapsing complex datasets into more manageable units, providing more nuanced understandings of both the organism-level effects of modified gene expression, and the targets of adaptive divergence. Here, we reanalyze a dataset that investigated temperature-induced changes in gene expression in marine-adapted and freshwater-adapted threespine stickleback (Gasterosteus aculeatus), using Weighted Gene Co-expression Network Analysis (WGCNA) with PANTHER Gene Ontology (GO)-Slim overrepresentation and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. Six modules exhibited a conserved response and six a divergent response between marine and freshwater stickleback when acclimated to 7°C or 22°C. One divergent module showed freshwater-specific response to temperature, and the remaining divergent modules showed differences in height of reaction norms. PPARAa, a transcription factor that regulates fatty acid metabolism and has been implicated in adaptive divergence, was located in a module that had higher expression at 7°C and in freshwater stickleback. This updated methodology revealed patterns that were not found in the original publication. Although such methods hold promise toward predicting population response to environmental stressors, many limitations remain, particularly with regard to module expression representation, database resources, and cross-database integration.

Through veiled mirrors: Fish fins giving insight into size regulation

Faithful establishment and maintenance of proportion is seen across biological systems and provides a glimpse at fundamental rules of scaling that underlie development and evolution. Dysregulation of proportion is observed in a range of human diseases and growth disorders, indicating that proper scaling is an essential component of normal anatomy and physiology. However, when viewed through an evolutionary lens, shifts in the regulation of relative proportion are one of the most striking sources of morphological diversity among organisms. To date, the mechanisms via which relative proportion is specified and maintained remain unclear. Through the application of powerful experimental, genetic and molecular approaches, the teleost fin has provided an effective model to investigate the regulation of scaling, size, and relative growth in vertebrate organisms. This article is categorized under: Establishment of Spatial and Temporal Patterns > Regulation of Size, Proportion, and Timing Adult Stem Cells, Tissue Renewal, and Regeneration > Regeneration Comparative Development and Evolution > Regulation of Organ Diversity.

A Temporal Map of Gene Expression Pattern During Zebrafish Liver Regeneration

Zebrafish is increasingly being used to study liver injury and regeneration. However, very little is known about molecular players that respond to injury and those important for liver regeneration. We use a metronidazole nitroreductase (MTZ-nfsb)-based system to selectively ablate hepatocytes in adult zebrafish to create a model for liver injury and regeneration. In this study, we generate a comprehensive temporal map of gene expression changes during regeneration through RNA sequencing of liver samples at various stages of injury and regeneration. Analyzing these data, we find that soon after injury the immediate early transcription factor MYC induces a battery of genes that respond to the MTZ-induced ROS by activating oxido-reductase pathways and apoptosis machinery. Immediately after injury, liver cells downregulate many functional genes, including complement protein synthesis, bile acid, and lipid biosynthesis, in a concerted manner. At 6 days postinjury, we find a dramatic induction of cholesterol biosynthesis and protein folding machinery, with expression levels returning to predamage levels by 8 days, suggesting an important role for these pathways in liver regeneration. This chronological transcriptomic map of liver regeneration in zebrafish would serve as a framework for further studies in understanding, and for screening for compounds that augment liver regeneration.

zfRegeneration: a database for gene expression profiling during regeneration

MOTIVATION

Zebrafish is a model organism with the ability to regenerate many different organs. Although RNA-Seq has been used extensively to study this process, there are no databases that allow easy access to data.

RESULTS

Here we develop the first regeneration database that provides easy access to a large number of RNA-Seq datasets through custom-made plots of expression levels, differential expression analyses, correlations of genes and comparisons of the different datasets. zfRegeneration has a user-friendly web interface designed to enhance regeneration studies and to overcome the barriers between different research groups that study the regeneration of distinct organs. Using several case studies, we demonstrate that zfRegeneration provides a unique platform to analyse and understand gene expression during regeneration.

AVAILABILITY AND IMPLEMENTATION

zfRegeneration is freely available at www.zfregeneration.org.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Emerging Roles for Immune Cells and MicroRNAs in Modulating the Response to Cardiac Injury

Stimulating cardiomyocyte regeneration after an acute injury remains the central goal in cardiovascular regenerative biology. While adult mammals respond to cardiac damage with deposition of rigid scar tissue, adult zebrafish and salamander unleash a regenerative program that culminates in new cardiomyocyte formation, resolution of scar tissue, and recovery of heart function. Recent studies have shown that immune cells are key to regulating pro-inflammatory and pro-regenerative signals that shift the injury microenvironment toward regeneration. Defining the genetic regulators that control the dynamic interplay between immune cells and injured cardiac tissue is crucial to decoding the endogenous mechanism of heart regeneration. In this review, we discuss our current understanding of the extent that macrophage and regulatory T cells influence cardiomyocyte proliferation and how microRNAs (miRNAs) regulate their activity in the injured heart.