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.

Dissecting the spatiotemporal diversity of adult neural stem cells

Adult stem cells are important for tissue turnover and regeneration. However, in most adult systems it remains elusive how stem cells assume different functional states and support spatially patterned tissue architecture. Here, we dissected the diversity of neural stem cells in the adult zebrafish brain, an organ that is characterized by pronounced zonation and high regenerative capacity. We combined single-cell transcriptomics of dissected brain regions with massively parallel lineage tracing and in vivo RNA metabolic labeling to analyze the regulation of neural stem cells in space and time. We detected a large diversity of neural stem cells, with some subtypes being restricted to a single brain region, while others were found globally across the brain. Global stem cell states are linked to neurogenic differentiation, with different states being involved in proliferative and non-proliferative differentiation. Our work reveals principles of adult stem cell organization and establishes a resource for the functional manipulation of neural stem cell subtypes.

Inference of differentiation trajectories by transfer learning across biological processes

Stem cells differentiate into distinct fates by transitioning through a series of transcriptional states. Current computational approaches allow reconstruction of differentiation trajectories from single-cell transcriptomics data, but it remains unknown to what degree differentiation can be predicted across biological processes. Here, we use transfer learning to infer differentiation processes and quantify predictability in early embryonic development and adult hematopoiesis. Overall, we find that non-linear methods outperform linear approaches, and we achieved the best predictions with a custom variational autoencoder that explicitly models changes in transcriptional variance. We observed a high accuracy of predictions in embryonic development, but we found somewhat lower agreement with the real data in adult hematopoiesis. We demonstrate that this discrepancy can be explained by a higher degree of concordant transcriptional processes along embryonic differentiation compared with adult homeostasis. In summary, we establish a framework for quantifying and exploiting predictability of cellular differentiation trajectories.

A single-cell RNA labeling strategy for measuring stress response upon tissue dissociation

Tissue dissociation, a crucial step in single-cell sample preparation, can alter the transcriptional state of a sample through the intrinsic cellular stress response. Here we demonstrate a general approach for measuring transcriptional response during sample preparation. In our method, transcripts made during dissociation are labeled for later identification upon sequencing. We found general as well as cell-type-specific dissociation response programs in zebrafish larvae, and we observed sample-to-sample variation in the dissociation response of mouse cardiomyocytes despite well-controlled experimental conditions. Finally, we showed that dissociation of the mouse hippocampus can lead to the artificial activation of microglia. In summary, our approach facilitates experimental optimization of dissociation procedures as well as computational removal of transcriptional perturbation response.

Long-term imaging reveals behavioral plasticity during C. elegans dauer exit

BACKGROUND

During their lifetime, animals must adapt their behavior to survive in changing environments. This ability requires the nervous system to undergo adjustments at distinct temporal scales, from short-term dynamic changes in expression of neurotransmitters and receptors to longer-term growth, spatial and connectivity reorganization, while integrating external stimuli. The nematode Caenorhabditis elegans provides a model of nervous system plasticity, in particular its dauer exit decision. Under unfavorable conditions, larvae will enter the non-feeding and non-reproductive stress-resistant dauer stage and adapt their behavior to cope with the harsh new environment, with active reversal under improved conditions leading to resumption of reproductive development. However, how different environmental stimuli regulate the exit decision mechanism and thereby drive the larva's behavioral change is unknown. To fill this gap and provide insights on behavioral changes over extended periods of time, we developed a new open hardware method for long-term imaging (12h) of C. elegans larvae.

RESULTS

Our WormObserver platform comprises open hardware and software components for video acquisition, automated processing of large image data (> 80k images/experiment) and data analysis. We identified dauer-specific behavioral motifs and characterized the behavioral trajectory of dauer exit in different environments and genetic backgrounds to identify key decision points and stimuli promoting dauer exit. Combining long-term behavioral imaging with transcriptomics data, we find that bacterial ingestion triggers a change in neuropeptide gene expression to establish post-dauer behavior.

CONCLUSIONS

Taken together, we show how a developing nervous system can robustly integrate environmental changes activate a developmental switch and adapt the organism's behavior to a new environment. WormObserver is generally applicable to other research questions within and beyond the C. elegans field, having a modular and customizable character and allowing assessment of behavioral plasticity over longer periods.

Origin and function of activated fibroblast states during zebrafish heart regeneration

The adult zebrafish heart has a high capacity for regeneration following injury. However, the composition of the regenerative niche has remained largely elusive. Here, we dissected the diversity of activated cell states in the regenerating zebrafish heart based on single-cell transcriptomics and spatiotemporal analysis. We observed the emergence of several transient cell states with fibroblast characteristics following injury, and we outlined the proregenerative function of collagen-12-expressing fibroblasts. To understand the cascade of events leading to heart regeneration, we determined the origin of these cell states by high-throughput lineage tracing. We found that activated fibroblasts were derived from two separate sources: the epicardium and the endocardium. Mechanistically, we determined Wnt signalling as a regulator of the endocardial fibroblast response. In summary, our work identifies specialized activated fibroblast cell states that contribute to heart regeneration, thereby opening up possible approaches to modulating the regenerative capacity of the vertebrate heart.

Single-cell RNA sequencing and spatiotemporal analysis of the regenerating zebrafish heart identify transient proregenerative fibroblast-like cells that are derived from the epicardium and the endocardium. Wnt signalling regulates the endocardial fibroblast response.

A single-cell atlas of de novo β-cell regeneration reveals the contribution of hybrid β/δ-cells to diabetes recovery in zebrafish

Regeneration-competent species possess the ability to reverse the progression of severe diseases by restoring the function of the damaged tissue. However, the cellular dynamics underlying this capability remain unexplored. Here, we have used single-cell transcriptomics to map de novo β-cell regeneration during induction and recovery from diabetes in zebrafish. We show that the zebrafish has evolved two distinct types of somatostatin-producing δ-cells, which we term δ1- and δ2-cells. Moreover, we characterize a small population of glucose-responsive islet cells, which share the hormones and fate-determinants of both β- and δ1-cells. The transcriptomic analysis of β-cell regeneration reveals that β/δ hybrid cells provide a prominent source of insulin expression during diabetes recovery. Using in vivo calcium imaging and cell tracking, we further show that the hybrid cells form de novo and acquire glucose-responsiveness in the course of regeneration. The overexpression of dkk3, a gene enriched in hybrid cells, increases their formation in the absence of β-cell injury. Finally, interspecies comparison shows that plastic δ1-cells are partially related to PP cells in the human pancreas. Our work provides an atlas of β-cell regeneration and indicates that the rapid formation of glucose-responsive hybrid cells contributes to the resolution of diabetes in zebrafish.

Single-cell-resolved dynamics of chromatin architecture delineate cell and regulatory states in zebrafish embryos

DNA accessibility of cis-regulatory elements (CREs) dictates transcriptional activity and drives cell differentiation during development. While many genes regulating embryonic development have been identified, the underlying CRE dynamics controlling their expression remain largely uncharacterized. To address this, we produced a multimodal resource and genomic regulatory map for the zebrafish community, which integrates single-cell combinatorial indexing assay for transposase-accessible chromatin with high-throughput sequencing (sci-ATAC-seq) with bulk histone PTMs and Hi-C data to achieve a genome-wide classification of the regulatory architecture determining transcriptional activity in the 24-h post-fertilization (hpf) embryo. We characterized the genome-wide chromatin architecture at bulk and single-cell resolution, applying sci-ATAC-seq on whole 24-hpf stage zebrafish embryos, generating accessibility profiles for ∼23,000 single nuclei. We developed a genome segmentation method, ScregSeg (single-cell regulatory landscape segmentation), for defining regulatory programs, and candidate CREs, specific to one or more cell types. We integrated the ScregSeg output with bulk measurements for histone post-translational modifications and 3D genome organization and identified new regulatory principles between chromatin modalities prevalent during zebrafish development. Sci-ATAC-seq profiling of npas4l/cloche mutant embryos identified novel cellular roles for this hematovascular transcriptional master regulator and suggests an intricate mechanism regulating its expression. Our work defines regulatory architecture and principles in the zebrafish embryo and establishes a resource of cell-type-specific genome-wide regulatory annotations and candidate CREs, providing a valuable open resource for genomics, developmental, molecular, and computational biology.

OUP accepted manuscript

AIMS

Arrhythmogenic cardiomyopathy (ACM) is an inherited cardiac disorder that is characterized by progressive loss of myocardium that is replaced by fibro-fatty cells, arrhythmias, and sudden cardiac death. While myocardial degeneration and fibro-fatty replacement occur in specific locations, the underlying molecular changes remain poorly characterized. Here, we aim to delineate local changes in gene expression to identify new genes and pathways that are relevant for specific remodelling processes occurring during ACM.

METHODS AND RESULTS

Using Tomo-Seq, genome-wide transcriptional profiling with high spatial resolution, we created transmural epicardial-to-endocardial gene expression atlases of explanted ACM hearts to gain molecular insights into disease-driving processes. This enabled us to link gene expression profiles to the different regional remodelling responses and allowed us to identify genes that are potentially relevant for disease progression. In doing so, we identified distinct gene expression profiles marking regions of cardiomyocyte degeneration and fibro-fatty remodelling and revealed Zinc finger and BTB domain-containing protein 11 (ZBTB11) to be specifically enriched at sites of active fibro-fatty replacement of myocardium. Immunohistochemistry indicated ZBTB11 to be induced in cardiomyocytes flanking fibro-fatty areas, which could be confirmed in multiple cardiomyopathy patients. Forced overexpression of ZBTB11 induced autophagy and cell death-related gene programmes in human cardiomyocytes, leading to increased apoptosis.

CONCLUSION

Our study shows the power of Tomo-Seq to unveil new molecular mechanisms in human cardiomyopathy and uncovers ZBTB11 as a novel driver of cardiomyocyte loss.

Single cell biology-a Keystone Symposia report

Single cell biology has the potential to elucidate many critical biological processes and diseases, from development and regeneration to cancer. Single cell analyses are uncovering the molecular diversity of cells, revealing a clearer picture of the variation among and between different cell types. New techniques are beginning to unravel how differences in cell state-transcriptional, epigenetic, and other characteristics-can lead to different cell fates among genetically identical cells, which underlies complex processes such as embryonic development, drug resistance, response to injury, and cellular reprogramming. Single cell technologies also pose significant challenges relating to processing and analyzing vast amounts of data collected. To realize the potential of single cell technologies, new computational approaches are needed. On March 17-19, 2021, experts in single cell biology met virtually for the Keystone eSymposium "Single Cell Biology" to discuss advances both in single cell applications and technologies.

Spatio-temporal mRNA tracking in the early zebrafish embryo

Early stages of embryogenesis depend on subcellular localization and transport of maternal mRNA. However, systematic analysis of these processes is hindered by a lack of spatio-temporal information in single-cell RNA sequencing. Here, we combine spatially-resolved transcriptomics and single-cell RNA labeling to perform a spatio-temporal analysis of the transcriptome during early zebrafish development. We measure spatial localization of mRNA molecules within the one-cell stage embryo, which allows us to identify a class of mRNAs that are specifically localized at an extraembryonic position, the vegetal pole. Furthermore, we establish a method for high-throughput single-cell RNA labeling in early zebrafish embryos, which enables us to follow the fate of individual maternal transcripts until gastrulation. This approach reveals that many localized transcripts are specifically transported to the primordial germ cells. Finally, we acquire spatial transcriptomes of two xenopus species and compare evolutionary conservation of localized genes as well as enriched sequence motifs.

Systematic identification of A-to-I RNA editing in zebrafish development and adult organs

A-to-I RNA editing is a common post transcriptional mechanism, mediated by the Adenosine deaminase that acts on RNA (ADAR) enzymes, that increases transcript and protein diversity. The study of RNA editing is limited by the absence of editing maps for most model organisms, hindering the understanding of its impact on various physiological conditions. Here, we mapped the vertebrate developmental landscape of A-to-I RNA editing, and generated the first comprehensive atlas of editing sites in zebrafish. Tens of thousands unique editing events and 149 coding sites were identified with high-accuracy. Some of these edited sites are conserved between zebrafish and humans. Sequence analysis of RNA over seven developmental stages revealed high levels of editing activity in early stages of embryogenesis, when embryos rely on maternal mRNAs and proteins. In contrast to the other organisms studied so far, the highest levels of editing were detected in the zebrafish ovary and testes. This resource can serve as the basis for understanding of the role of editing during zebrafish development and maturity.

Single-cell genomics to study developmental cell fate decisions in zebrafish

New developments in single-cell genomics have transformed developmental biology in recent years by enabling systematic analysis of embryonic cell types and differentiation trajectories. Ongoing efforts in experimental and computational method development aim to reveal gene-regulatory mechanisms and to provide additional spatio-temporal information about developmental cell fate decisions. Here, we discuss recent technological developments as well as biological applications of single-cell genomics, with a particular focus on analysis of developmental cell fate decisions. Although the approaches described here are generally applicable to a broad range of model systems, we focus our discussion on applications in zebrafish, which has proven to be a particularly powerful model organism for establishing novel methods in single-cell genomics.

Variability of an Early Developmental Cell Population Underlies Stochastic Laterality Defects

Embryonic development seemingly proceeds with almost perfect precision. However, it is largely unknown how much underlying microscopic variability is compatible with normal development. Here, we quantify embryo-to-embryo variability in vertebrate development by studying cell number variation in the zebrafish endoderm. We notice that the size of a sub-population of the endoderm, the dorsal forerunner cells (DFCs, which later form the left-right organizer), exhibits significantly more embryo-to-embryo variation than the rest of the endoderm. We find that, with incubation of the embryos at elevated temperature, the frequency of left-right laterality defects is increased drastically in embryos with a low number of DFCs. Furthermore, we observe that these fluctuations have a large stochastic component among fish of the same genetic background. Hence, a stochastic variation in early development leads to a remarkably strong macroscopic phenotype. These fluctuations appear to be associated with maternal effects in the specification of the DFCs.

•

High embryo-to-embryo variability of dorsal forerunner cell numbers

•

Fluctuations of dorsal forerunner cells have a large stochastic component

•

Embryos with fewer dorsal forerunner cells frequently develop laterality defects

•

Variability of dorsal forerunner cell numbers is associated to maternal effects

The transcriptome dynamics of single cells during the cell cycle

The cell cycle is among the most basic phenomena in biology. Despite advances in single-cell analysis, dynamics and topology of the cell cycle in high-dimensional gene expression space remain largely unknown. We developed a linear analysis of transcriptome data which reveals that cells move along a planar circular trajectory in transcriptome space during the cycle. Non-cycling gene expression adds a third dimension causing helical motion on a cylinder. We find in immortalized cell lines that cell cycle transcriptome dynamics occur largely independently from other cellular processes. We offer a simple method ("Revelio") to order unsynchronized cells in time. Precise removal of cell cycle effects from the data becomes a straightforward operation. The shape of the trajectory implies that each gene is upregulated only once during the cycle, and only two dynamic components represented by groups of genes drive transcriptome dynamics. It indicates that the cell cycle has evolved to minimize changes of transcriptional activity and the related regulatory effort. This design principle of the cell cycle may be of relevance to many other cellular differentiation processes.

Inclusion of temporal information in single cell transcriptomics

Single cell transcriptomics has emerged as a powerful method for dissecting cell type diversity and for understanding mechanisms of cell fate decisions. However, inclusion of temporal information remains challenging, since each cell can be measured only once by sequencing analysis. Here, we discuss recent progress and current efforts towards inclusion of temporal information in single cell transcriptomics. Even from snapshot data, temporal dynamics can be computationally inferred via pseudo-temporal ordering of single cell transcriptomes. Temporal information can also come from analysis of intronic reads or from RNA metabolic labeling, which can provide additional evidence for pseudo-time trajectories and enable more fine-grained analysis of gene regulatory interactions. These approaches measure dynamics on short timescales of hours. Emerging methods for high-throughput lineage tracing now enable information storage over long timescales by using CRISPR/Cas9 to record information in the genome, which can later be read out by sequencing.

RNA Tomography for Spatially Resolved Transcriptomics (Tomo-Seq)

Embryonic development is heavily dependent on temporally and spatially restricted gene expression. Spatially resolved measurements of gene expression are therefore crucial for identifying novel regulators and the understanding of their function. However, in situ methods do not resolve global gene expression, and sequencing-based methods usually do not provide spatial information. Here, we describe tomo-seq, a method that combines classical histological sectioning of embryos or tissues with a highly sensitive RNA-sequencing technique. Application of tomo-seq to zebrafish embryos allows reconstructing the spatial gene expression of thousands of genes.

Simultaneous lineage tracing and cell-type identification using CRISPR-Cas9-induced genetic scars

A key goal of developmental biology is to understand how a single cell transforms into a full-grown organism comprising many different cell types. Single-cell RNA-sequencing (scRNA-seq) is commonly used to identify cell types in a tissue or organ1. However, organizing the resulting taxonomy of cell types into lineage trees to understand developmental origin of cells remains challenging. Here we present LINNAEUS (LINeage tracing by Nuclease-Activated Editing of Ubiquitous Sequences)—a strategy for simultaneous lineage tracing and transcriptome profiling in thousands of single cells. By combining scRNA-seq with computational analysis of lineage barcodes, generated by genome editing of transgenic reporter genes, we reconstruct developmental lineage trees in zebrafish larvae, and in heart, liver, pancreas and telencephalon of adult fish. LINNAEUS provides a systematic approach for tracing the origin of novel cell types, or known cell types under different conditions.

Methods for lineage tracing on the organism-wide level

Determining the lineage origin of cell types is a major goal in developmental biology. Furthermore, lineage tracing is a powerful approach for understanding the origin of developmental defects as well as the origin of diseases such as cancer. There is now a variety of complementary approaches for identifying lineage relationships, ranging from direct observation of cell divisions by light microscopy to genetic labeling of cells using inducible recombinases and fluorescent reporters. A recent development, and the main topic of this review article, is the use of high-throughput sequencing data for lineage analysis. This emerging approach holds the promise of increased multiplexing capacity, allowing lineage analysis of large cell numbers up to the organism-wide level combined with simultaneous transcription profiling by single cell RNA sequencing.

Tomo-Seq Identifies SOX9 as a Key Regulator of Cardiac Fibrosis During Ischemic Injury

Background:

Cardiac ischemic injury induces a pathological remodeling response, which can ultimately lead to heart failure. Detailed mechanistic insights into molecular signaling pathways relevant for different aspects of cardiac remodeling will support the identification of novel therapeutic targets.

Methods:

Although genome-wide transcriptome analysis on diseased tissues has greatly advanced our understanding of the regulatory networks that drive pathological changes in the heart, this approach has been disadvantaged by the fact that the signals are derived from tissue homogenates. Here we used tomo-seq to obtain a genome-wide gene expression signature with high spatial resolution spanning from the infarcted area to the remote to identify new regulators of cardiac remodeling. Cardiac tissue samples from patients suffering from ischemic heart disease were used to validate our findings.

Results:

Tracing transcriptional differences with a high spatial resolution across the infarcted heart enabled us to identify gene clusters that share a comparable expression profile. The spatial distribution patterns indicated a separation of expressional changes for genes involved in specific aspects of cardiac remodeling, such as fibrosis, cardiomyocyte hypertrophy, and calcium handling ( Col1a2 , Nppa , and Serca2 ). Subsequent correlation analysis allowed for the identification of novel factors that share a comparable transcriptional regulation pattern across the infarcted tissue. The strong correlation between the expression levels of these known marker genes and the expression of the coregulated genes could be confirmed in human ischemic cardiac tissue samples. Follow-up analysis identified SOX9 as common transcriptional regulator of a large portion of the fibrosis-related genes that become activated under conditions of ischemic injury. Lineage-tracing experiments indicated that the majority of COL1-positive fibroblasts stem from a pool of SOX9-expressing cells, and in vivo loss of Sox9 blunted the cardiac fibrotic response on ischemic injury. The colocalization between SOX9 and COL1 could also be confirmed in patients suffering from ischemic heart disease.

Conclusions:

Based on the exact local expression cues, tomo-seq can serve to reveal novel genes and key transcription factors involved in specific aspects of cardiac remodeling. Using tomo-seq, we were able to unveil the unknown relevance of SOX9 as a key regulator of cardiac fibrosis, pointing to SOX9 as a potential therapeutic target for cardiac fibrosis.

Patched Receptors Sense, Interpret, and Establish an Epidermal Hedgehog Signaling Gradient

By using the sensitivity of single-molecule fluorescent in situ hybridization, we have precisely quantified the levels and defined the temporal and spatial distribution of Hedgehog signaling activity during embryonic skin development and discovered that there is a Hedgehog signaling gradient along the proximal-distal axis of developing hair follicles. To explore the contribution of Hedgehog receptors Ptch1 and Ptch2 in establishing the epidermal signaling gradient, we quantitated the level of pathway activity generated in Ptch1- and Ptch1;Ptch2-deficient skin and defined the contribution of each receptor to regulation of the levels of Hedgehog signaling identified in wild-type skin. Moreover, we show that both the cellular phenotype and level of pathway activity featured in Ptch1;Ptch2-deficient cells faithfully recapitulates the Peak level of endogenous Hedgehog signaling detected at the base of developing follicles, where the concentration of endogenous Shh is predicted to be highest. Taken together, these data show that both Ptch1 and Ptch2 play a crucial role in sensing the concentration of Hedgehog ligand and regulating the appropriate dose-dependent response.

Tomo-seq: A method to obtain genome-wide expression data with spatial resolution

To improve our understanding of pattern formation during development and disease we heavily rely on the identification of novel regulators and pathways. While RNA sequencing yields genome-wide expression data that suit this purpose, it lacks spatial resolution. Such spatial resolution can be obtained by microscopy-based methods like in situ hybridization, but these fail to provide information on more than a few genes at a time. Here, we describe tomo-seq, a technique that combines the advantages of the above-mentioned approaches and provides genome-wide expression data with spatial information. The tomo-seq technique is based on cryosectioning of an embryo or tissue of interest and performing RNA-seq on individual sections. Using this method, we have generated genome-wide transcriptomics with high spatial resolution of the whole zebrafish embryo at various stages of development (Junker et al., 2014) and of adult zebrafish hearts after injury (Wu et al., 2016).

Spatially Resolved Genome-wide Transcriptional Profiling Identifies BMP Signaling as Essential Regulator of Zebrafish Cardiomyocyte Regeneration

In contrast to mammals, zebrafish regenerate heart injuries via proliferation of cardiomyocytes located near the wound border. To identify regulators of cardiomyocyte proliferation, we used spatially resolved RNA sequencing (tomo-seq) and generated a high-resolution genome-wide atlas of gene expression in the regenerating zebrafish heart. Interestingly, we identified two wound border zones with distinct expression profiles, including the re-expression of embryonic cardiac genes and targets of bone morphogenetic protein (BMP) signaling. Endogenous BMP signaling has been reported to be detrimental to mammalian cardiac repair. In contrast, we find that genetic or chemical inhibition of BMP signaling in zebrafish reduces cardiomyocyte dedifferentiation and proliferation, ultimately compromising myocardial regeneration, while bmp2b overexpression is sufficient to enhance it. Our results provide a resource for further studies on the molecular regulation of cardiac regeneration and reveal intriguing differential cellular responses of cardiomyocytes to a conserved signaling pathway in regenerative versus non-regenerative hearts.

Single-cell transcriptomics enters the age of mass production

Two publications in the current issue of Cell introduce novel methods for high-throughput single-cell transcriptomics by using droplet microfluidics and sophisticated barcoding schemes for transcriptional profiling of thousands of individual cells.

Abstract 373: Spatial Expression Analysis of the Ischemic Heart Identifies Novel Factors Involved in Cardiac Remodeling and Disease

Cardiac ischemic injury is associated with expressional changes in a broad range of coding and non-coding genes. While genome-wide expression techniques like RNA sequencing, give valuable information about the transcriptional activity of genes, spatial clues are lost.

Here, we use a novel method (TomoSEQ) based on traditional histological techniques with low-input RNA sequencing to generate a spatial gene expression atlas for the infarcted heart. Transcriptome profiling of 50 consecutive sections from hearts exposed to ischemia reperfusion or sham surgery, allowed us to trace differential gene regulation throughout the infarct, borderzone and remote area. Fourteen days post-injury, we clearly observed differential expression curves specific for genes known to be relevant for cardiac function and disease, like Nppa (remodeling/hypertrophy), Col3a1 (fibrosis) and Atp2a2 (calcium handling), with the first two being more abundantly expressed in the infarct, and the latter more dominant towards the remote. Interestingly, the exact spatial information on these 3 reference genes enabled us to identify a set of genes that show an identical regulation, potentially linking these transcripts to the aforementioned biological function of the reference genes. The validity of using this method to identify genes that are expressionally linked was underscored by the observation that the correlation between the expression level of the Nppa , Col3a1 , Atp2a2 and the newly identified genes could be confirmed in human cardiac tissue samples from patients suffering from ischemic heart disease. We next used an in silico approach to identify common transcription factors among the co-regulated genes. TomoSEQ analysis revealed a strong spatial correlation between these identified transcription factors and their target genes. This provides further support for our methodology, but also provides an opportunity to identify novel transcription factors relevant for cardiac function and disease.

Our study shows the power of spatial transcriptome profiling throughout the infarcted heart by giving local clues on changes in gene expression and allowing for the identification of novel genes and key transcription factors that could be relevant for cardiac repair upon injury.

DAZL regulates Tet1 translation in murine embryonic stem cells

Embryonic stem cell (ESC) cultures display a heterogeneous gene expression profile, ranging from a pristine naïve pluripotent state to a primed epiblast state. Addition of inhibitors of GSK3β and MEK (so-called 2i conditions) pushes ESC cultures toward a more homogeneous naïve pluripotent state, but the molecular underpinnings of this naïve transition are not completely understood. Here, we demonstrate that DAZL, an RNA-binding protein known to play a key role in germ-cell development, marks a subpopulation of ESCs that is actively transitioning toward naïve pluripotency. Moreover, DAZL plays an essential role in the active reprogramming of cytosine methylation. We demonstrate that DAZL associates with mRNA of Tet1, a catalyst of 5-hydroxylation of methyl-cytosine, and enhances Tet1 mRNA translation. Overexpression of DAZL in heterogeneous ESC cultures results in elevated TET1 protein levels as well as increased global hydroxymethylation. Conversely, null mutation of Dazl severely stunts 2i-mediated TET1 induction and hydroxymethylation. Our results provide insight into the regulation of the acquisition of naïve pluripotency and demonstrate that DAZL enhances TET1-mediated cytosine hydroxymethylation in ESCs that are actively reprogramming to a pluripotent ground state.

Licensing of primordial germ cells for gametogenesis depends on genital ridge signaling

In mouse embryos at mid-gestation, primordial germ cells (PGCs) undergo licensing to become gametogenesis-competent cells (GCCs), gaining the capacity for meiotic initiation and sexual differentiation. GCCs then initiate either oogenesis or spermatogenesis in response to gonadal cues. Germ cell licensing has been considered to be a cell-autonomous and gonad-independent event, based on observations that some PGCs, having migrated not to the gonad but to the adrenal gland, nonetheless enter meiosis in a time frame parallel to ovarian germ cells -- and do so regardless of the sex of the embryo. Here we test the hypothesis that germ cell licensing is cell-autonomous by examining the fate of PGCs in Gata4 conditional mutant (Gata4 cKO) mouse embryos. Gata4, which is expressed only in somatic cells, is known to be required for genital ridge initiation. PGCs in Gata4 cKO mutants migrated to the area where the genital ridge, the precursor of the gonad, would ordinarily be formed. However, these germ cells did not undergo licensing and instead retained characteristics of PGCs. Our results indicate that licensing is not purely cell-autonomous but is induced by the somatic genital ridge.

During embryonic development, stem cell-like primordial germ cells travel across the developing embryo to the genital ridge, which gives rise to the gonad. Around the time of their arrival, the primordial germ cells gain the capacity to undertake sexual specialization and meiosis—a process called germ cell licensing. Based on the observation that meiosis and sexual differentiation can occur when primordial germ cells stray into the area of the adrenal gland, the primordial germ cell has been thought to be responsible for its own licensing. We tested this notion by examining the licensing process in mutant mouse embryos that did not form a genital ridge. We discovered that in the absence of the genital ridge, primordial germ cells migrate across the developing embryo properly, but instead of undergoing licensing, these cells retain their primordial germ cell characteristics. We conclude that licensing of embryonic primordial germ cells for gametogenesis is dependent on signaling from the genital ridge.

Ascl2 acts as an R-spondin/Wnt-responsive switch to control stemness in intestinal crypts

The Wnt signaling pathway controls stem cell identity in the intestinal epithelium and in many other adult organs. The transcription factor Ascl2 (a Wnt target gene) is a master regulator of intestinal stem cell identity. It is unclear how the continuous Wnt gradient along the crypt axis is translated into discrete expression of Ascl2 and discrete specification of stem cells at crypt bottoms. We show that (1) Ascl2 is regulated in a direct autoactivatory loop, leading to a distinct on/off expression pattern, and (2) Wnt/R-spondin can activate this regulatory loop. This mechanism interprets the Wnt levels in the intestinal crypt and translates the continuous Wnt signal into a discrete Ascl2 "on" or "off" decision. In turn, Ascl2, together with β-catenin/Tcf, activates the genes fundamental to the stem cell state. In this manner, Ascl2 forms a transcriptional switch that is both Wnt responsive and Wnt dependent to define stem cell identity.

A predictive model of bifunctional transcription factor signaling during embryonic tissue patterning

Hedgehog signaling controls pattern formation in many vertebrate tissues. The downstream effectors of the pathway are the bifunctional Gli transcription factors, which, depending on hedgehog concentration, act as either transcriptional activators or repressors. Quantitatively understanding the interplay between Gli activator and repressor forms for patterning complex tissues is an open challenge. Here, we describe a reductionist mathematical model for how Gli activators and repressors are integrated in space and time to regulate transcriptional outputs of hedgehog signaling, using the pathway readouts Gli1 and Ptch1 as a model system. Spatially resolved measurements of absolute transcript numbers for these genes allow us to infer spatiotemporal variations of Gli activator and repressor levels. We validate our model by successfully predicting expression changes of Gli1 and Ptch1 in mutants at different developmental stages and in different tissues. Our results provide a starting point for understanding gene regulation by bifunctional transcription factors during mammalian development.

Neural-specific Sox2 input and differential Gli-binding affinity provide context and positional information in Shh-directed neural patterning

In the vertebrate neural tube, regional Sonic hedgehog (Shh) signaling invokes a time- and concentration-dependent induction of six different cell populations mediated through Gli transcriptional regulators. Elsewhere in the embryo, Shh/Gli responses invoke different tissue-appropriate regulatory programs. A genome-scale analysis of DNA binding by Gli1 and Sox2, a pan-neural determinant, identified a set of shared regulatory regions associated with key factors central to cell fate determination and neural tube patterning. Functional analysis in transgenic mice validates core enhancers for each of these factors and demonstrates the dual requirement for Gli1 and Sox2 inputs for neural enhancer activity. Furthermore, through an unbiased determination of Gli-binding site preferences and analysis of binding site variants in the developing mammalian CNS, we demonstrate that differential Gli-binding affinity underlies threshold-level activator responses to Shh input. In summary, our results highlight Sox2 input as a context-specific determinant of the neural-specific Shh response and differential Gli-binding site affinity as an important cis-regulatory property critical for interpreting Shh morphogen action in the mammalian neural tube.

When noisy neighbors are a blessing: analysis of gene expression noise identifies coregulated genes

In this issue of Molecular Cell, Stewart-Ornstein et al. (2012) use systematic pair-wise correlation analysis of expression noise in a large number of yeast genes to identify clusters of functionally related genes and signaling pathways responsible for elevated noise.

Ligand-dependent equilibrium fluctuations of single calmodulin molecules

Single-molecule force spectroscopy allows superb mechanical control of protein conformation. We used a custom-built low-drift atomic force microscope to observe mechanically induced conformational equilibrium fluctuations of single molecules of the eukaryotic calcium-dependent signal transducer calmodulin (CaM). From this data, the ligand dependence of the full energy landscape can be reconstructed. We find that calcium ions affect the folding kinetics of the individual CaM domains, whereas target peptides stabilize the already folded structure. Single-molecule data of full length CaM reveal that a wasp venom peptide binds noncooperatively to CaM with 2:1 stoichiometry, whereas a target enzyme peptide binds cooperatively with 1:1 stoichiometry. If mechanical load is applied directly to the target peptide, real-time binding/unbinding transitions can be observed.

Cysteine engineering of polyproteins for single-molecule force spectroscopy

Single-molecule methods such as force spectroscopy give experimental access to the mechanical properties of protein molecules. So far, owing to the limitations of recombinant construction of polyproteins, experimental access has been limited to mostly the N-to-C terminal direction of force application. This protocol gives a fast and simple alternative to current recombinant strategies for preparing polyproteins. We describe in detail the method to construct polyproteins with precisely controlled linkage topologies, based on the pairwise introduction of cysteines into protein structure and subsequent polymerization in solution. Stretching such constructed polyproteins in an atomic force microscope allows mechanical force application to a single protein structure via two precisely controlled amino acid positions in the functional three-dimensional protein structure. The capability for site-directed force application can provide valuable information about both protein structure and directional protein mechanics. This protocol should be applicable to almost any protein that can be point mutated. Given correct setup of all necessary reagents, this protocol can be accomplished in fewer than 10 d.

Influence of substrate binding on the mechanical stability of mouse dihydrofolate reductase

We investigated the effect of substrate binding on the mechanical stability of mouse dihydrofolate reductase using single-molecule force spectroscopy by atomic force microscopy. We find that under mechanical forces dihydrofolate reductase unfolds via a metastable intermediate with lifetimes on the millisecond timescale. Based on the measured length increase of approximately 22 nm we suggest a structure for this intermediate with intact substrate binding sites. In the presence of the substrate analog methotrexate and the cofactor NADPH lifetimes of this intermediate are increased by up to a factor of two. Comparing mechanical and thermodynamic stabilization effects of substrate binding suggests mechanical stability is dominated by local interactions within the protein structure. These experiments demonstrate that protein mechanics can be used to probe the substrate binding status of an enzyme.