Robust Identification of Developmentally Active Endothelial Enhancers in Zebrafish Using FANS-Assisted ATAC-Seq

Emerging Frontiers in Zebrafish Embryonic and Adult-Derived Cell Lines

Zebrafish (Danio rerio) has become a pivotal vertebrate model in biomedical research, renowned for its genetic similarity to humans, optical transparency, rapid embryonic development, and amenability to experimental manipulation. In recent years, the derivation of cell lines from zebrafish embryos has unlocked new possibilities for in vitro studies across developmental biology, toxicology, disease modeling, and genetic engineering. These embryo-derived cultures offer scalable, reproducible, and ethically favorable alternatives to in vivo approaches, enabling high-throughput screening and mechanistic exploration under defined conditions. This review provides a comprehensive overview of protocols for establishing and maintaining zebrafish embryonic cell lines, emphasizing culture conditions, pluripotency features, transfection strategies, and recent innovations such as genotype-defined mutant lines generated via CRISPR/Cas9 and feeder-free systems. We also highlight emerging applications in oncology, regenerative medicine, and functional genomics, positioning zebrafish cell lines as versatile platforms bridging animal models and next-generation in vitro systems. Its continued optimization holds promise for improved reproducibility, reduced animal use, and expanded translational impact in biomedical research.

Integration of ATAC-seq and RNA-seq reveals the dynamics of chromatin accessibility and gene expression in zoysiagrass response to drought

KEY MESSAGE

The 'X4' accession of zoysiagrass demonstrated superior drought tolerance compared to other accessions. Integration analysis of transcriptomics and epigenomics revealed a positive correlation between ATAC-seq peak intensity and gene expression levels. Six motifs involved in regulating drought responses were identified, which are similar to the domains of the ERF and C2H2 transcription factor families. Heterologous expression of Zja11G000860 in yeast enhanced tolerance to drought stress, allowing robust growth even at high PEG6000 concentrations. Zoysiagrass is renowned for its drought tolerance and serves as an exceptional domestic turfgrass in China. However, the changes in chromatin accessibility during drought in zoysiagrass are not well understood. We conducted a preliminary evaluation of the phenotypic changes in drought tolerance for six zoysiagrass cultivars, taking into account their growth characteristics and physiological traits under drought conditions. Additionally, we utilized an integrated multi-omics strategy, encompassing RNA sequencing (RNA-seq), Assay for Transposase Accessible Chromatin using high-throughput sequencing (ATAC-seq), and reverse transcription quantitative PCR (RT-qPCR) verification experiments, to gain deeper understanding of the chromatin accessibility patterns linked to gene expression in response to drought stress in zoysiagrass. Preliminary analysis of the trends in relative water content and proline content suggested that the variety 'X4' exhibited superior drought tolerance compared to the other five accessions. The KEGG pathway enrichment analysis revealed that zoysiagrass responded to environmental stress by regulating stress response and antioxidant defense pathways. Notably, the expression levels of genes Zja03G031540 and Zja11G000860 were significantly increased in the 'X4' zoysiagrass genotype, which exhibited improved drought tolerance, compared to the 'X1' zoysiagrass genotype with reduced drought tolerance. This study suggested that 63 high-confidence genes are related to drought stress, including Zja03G031540 and Zja11G000860. Additionally, six motifs regulating drought responses were unearthed. Furthermore, the heterologous expression of Zja11G000860 in yeast enhanced tolerance to drought stress. The study discovered a positive correlation between ATAC-seq peak intensity and gene expression levels. The expression of high-confidence genes was linked to zoysiagrass resistance evaluation and phenotypic traits, implying that these genes are involved in responding to external drought stress. This study combined ATAC-seq and RNA-seq technologies for the first time to identify drought-related gene expression in zoysiagrass, elucidating the grass adaptation to environmental stress and the regulatory mechanisms underlying stress responses, and laying the groundwork for zoysiagrass improvement and breeding.

Mapping the chromatin accessibility landscape of zebrafish embryogenesis at single-cell resolution by SPATAC-seq

Currently, the dynamic accessible elements that determine regulatory programs responsible for the unique identity and function of each cell type during zebrafish embryogenesis lack detailed study. Here we present SPATAC-seq: a split-pool ligation-based assay for transposase-accessible chromatin using sequencing. Using SPATAC-seq, we profiled chromatin accessibility in more than 800,000 individual nuclei across 20 developmental stages spanning the sphere stage to the early larval protruding mouth stage. Using this chromatin accessibility map, we identified 604 cell states and inferred their developmental relationships. We also identified 959,040 candidate cis-regulatory elements (cCREs) and delineated development-specific cCREs, as well as transcription factors defining diverse cell identities. Importantly, enhancer reporter assays confirmed that the majority of tested cCREs exhibited robust enhanced green fluorescent protein expression in restricted cell types or tissues. Finally, we explored gene regulatory programs that drive pigment and notochord cell differentiation. Our work provides a valuable open resource for exploring driver regulators of cell fate decisions in zebrafish embryogenesis.

Multiple cis-regulatory elements control prox1a expression in distinct lymphatic vascular beds

During embryonic development, lymphatic endothelial cell (LEC) precursors are distinguished from blood endothelial cells by the expression of Prospero-related homeobox 1 (Prox1), which is essential for lymphatic vasculature formation in mouse and zebrafish. Prox1 expression initiation precedes LEC sprouting and migration, serving as the marker of specified LECs. Despite its crucial role in lymphatic development, Prox1 upstream regulation in LECs remains to be uncovered. SOX18 and COUP-TFII are thought to regulate Prox1 in mice by binding its promoter region. However, the specific regulation of Prox1 expression in LECs remains to be studied in detail. Here, we used evolutionary conservation and chromatin accessibility to identify enhancers located in the proximity of zebrafish prox1a active in developing LECs. We confirmed the functional role of the identified sequences through CRISPR/Cas9 mutagenesis of a lymphatic valve enhancer. The deletion of this region results in impaired valve morphology and function. Overall, our results reveal an intricate control of prox1a expression through a collection of enhancers. Ray-finned fish-specific distal enhancers drive pan-lymphatic expression, whereas vertebrate-conserved proximal enhancers refine expression in functionally distinct subsets of lymphatic endothelium.

Construction of single-cell cross-species chromatin accessibility landscapes with combinatorial-hybridization-based ATAC-seq

Despite recent advances in single-cell genomics, the lack of maps for single-cell candidate cis-regulatory elements (cCREs) in non-mammal species has limited our exploration of conserved regulatory programs across vertebrates and invertebrates. Here, we developed a combinatorial-hybridization-based method for single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) named CH-ATAC-seq, enabling the construction of single-cell accessible chromatin landscapes for zebrafish, Drosophila, and earthworms (Eisenia andrei). By integrating scATAC censuses of humans, monkeys, and mice, we systematically identified 152 distinct main cell types and around 0.8 million cell-type-specific cCREs. Our analysis provided insights into the conservation of neural, muscle, and immune lineages across species, while epithelial cells exhibited a higher organ-origin heterogeneity. Additionally, a large-scale gene regulatory network (GRN) was constructed in four vertebrates by integrating scRNA-seq censuses. Overall, our study provides a valuable resource for comparative epigenomics, identifying the evolutionary conservation and divergence of gene regulation across different species.

Computational identification of key transcription factors for embryonic and postnatal Sox2+ dental epithelial stem cell

While many reptiles can replace their tooth throughout life, human loss the tooth replacement capability after formation of the permanent teeth. It was thought that the difference in tooth regeneration capability depends on the persistence of a specialized dental epithelial structure, the dental lamina that contains dental epithelial stem cells (DESC). Currently, we know very little about DESC such as what genes are expressed or its chromatin accessibility profile. Multiple markers of DESC have been proposed such as Sox2 and Lgr5 . Few single cell RNA-seq experiments have been performed previously, but no obvious DESC cluster was identified in these scRNA-seq datasets, possible due to that the expression level of DESC markers such as Sox2 and Lgr5 is too low or the percentage of DESC is too low in whole tooth. We utilize a mouse line Sox2-GFP to enrich Sox2+ DESC and use Smart-Seq2 protocol and ATAC-seq protocol to generate transcriptome profile and chromatin accessibility profile of P2 Sox2+ DESC. Additionally, we generate transcriptome profile and chromatin accessibility profile of E11.5 Sox2+ dental lamina cells. With transcriptome profile and chromatin accessibility profile, we systematically identify potential key transcription factors for E11.5 Sox2+ cells and P2 Sox2+ cells. We identified transcription factors including Pitx2, Id3, Pitx1, Tbx1, Trp63, Nkx2-3, Grhl3, Dlx2, Runx1, Nfix, Zfp536 , etc potentially formed the core transcriptional regulatory networks of Sox2+ DESC in both embryonic and postnatal stages.

An injury-responsive mmp14b enhancer is required for heart regeneration

Mammals have limited capacity for heart regeneration, whereas zebrafish have extraordinary regeneration abilities. During zebrafish heart regeneration, endothelial cells promote cardiomyocyte cell cycle reentry and myocardial repair, but the mechanisms responsible for promoting an injury microenvironment conducive to regeneration remain incompletely defined. Here, we identify the matrix metalloproteinase Mmp14b as an essential regulator of heart regeneration. We identify a TEAD-dependent mmp14b endothelial enhancer induced by heart injury in zebrafish and mice, and we show that the enhancer is required for regeneration, supporting a role for Hippo signaling upstream of mmp14b. Last, we show that MMP-14 function in mice is important for the accumulation of Agrin, an essential regulator of neonatal mouse heart regeneration. These findings reveal mechanisms for extracellular matrix remodeling that promote heart regeneration.

Hand2 represses non-cardiac cell fates through chromatin remodeling at cis- regulatory elements

Developmental studies have revealed the importance of the transcription factor Hand2 in cardiac development. Hand2 promotes cardiac progenitor differentiation and epithelial maturation, while repressing other tissue types. The mechanisms underlying the promotion of cardiac fates are far better understood than those underlying the repression of alternative fates. Here, we assess Hand2-dependent changes in gene expression and chromatin remodeling in cardiac progenitors of zebrafish embryos. Cell-type specific transcriptome analysis shows a dual function for Hand2 in activation of cardiac differentiation genes and repression of pronephric pathways. We identify functional cis- regulatory elements whose chromatin accessibility are increased in hand2 mutant cells. These regulatory elements associate with non-cardiac gene expression, and drive reporter gene expression in tissues associated with Hand2-repressed genes. We find that functional Hand2 is sufficient to reduce non-cardiac reporter expression in cardiac lineages. Taken together, our data support a model of Hand2-dependent coordination of transcriptional programs, not only through transcriptional activation of cardiac and epithelial maturation genes, but also through repressive chromatin remodeling at the DNA regulatory elements of non-cardiac genes.

Regeneration and developmental enhancers are differentially compatible with minimal promoters

Enhancers and promoters are cis-regulatory elements that control gene expression. Enhancers are activated in a cell type-, tissue-, and condition-specific manner to stimulate promoter function and transcription. Zebrafish have emerged as a powerful animal model for examining the activities of enhancers derived from various species through transgenic enhancer assays, in which an enhancer is coupled with a minimal promoter. However, the efficiency of minimal promoters and their compatibility with multiple developmental and regeneration enhancers have not been systematically tested in zebrafish. Thus, we assessed the efficiency of six minimal promoters and comprehensively interrogated the compatibility of the promoters with developmental and regeneration enhancers. We found that the fos minimal promoter and Drosophila synthetic core promoter (DSCP) yielded high rates of leaky expression that may complicate the interpretation of enhancer assays. Notably, the adenovirus E1b promoter, the zebrafish lepb 0.8-kb (P0.8) and lepb 2-kb (P2) promoters, and a new zebrafish synthetic promoter (ZSP) that combines elements of the E1b and P0.8 promoters drove little or no ectopic expression, making them suitable for transgenic assays. We also found significant differences in compatibility among specific combinations of promoters and enhancers, indicating the importance of promoters as key regulatory elements determining the specificity of gene expression. Our study provides guidelines for transgenic enhancer assays in zebrafish to aid in the discovery of functional enhancers regulating development and regeneration.

A regulatory network of Sox and Six transcription factors initiate a cell fate transformation during hearing regeneration in adult zebrafish

Using adult zebrafish inner ears as a model for sensorineural regeneration, we ablated the mechanosensory receptors and characterized the single-cell epigenome and transcriptome at consecutive time points during hair cell regeneration. We utilized deep learning on the regeneration-induced open chromatin sequences and identified cell-specific transcription factor (TF) motif patterns. Enhancer activity correlated with gene expression and identified potential gene regulatory networks. A pattern of overlapping Sox- and Six-family TF gene expression and binding motifs was detected, suggesting a combinatorial program of TFs driving regeneration and cell identity. Pseudotime analysis of single-cell transcriptomic data suggested that support cells within the sensory epithelium changed cell identity to a “progenitor” cell population that could differentiate into hair cells. We identified a 2.6 kb DNA enhancer upstream of the sox2 promoter that, when deleted, showed a dominant phenotype that resulted in a hair-cell-regeneration-specific deficit in both the lateral line and adult inner ear.

Jimenez et al. interrogate the epigenomic and transcriptomic landscape of regenerating adult zebrafish inner-ear sensory epithelia. They show that the support-cell population transitions to an intermediate “progenitor” cell state that becomes new hair cells, and they demonstrate that the cell fate decisions may be driven by the coordinate regulation and spatial co-binding of Sox and Six transcription factors. By functionally validating a predicted regeneration-responsive enhancer upstream of sox2, they show that precise timing of sox2 expression is critical for hearing regeneration in zebrafish.

Chromatin conformation of human oral epithelium can identify orofacial cleft missing functional variants

Genome-wide association studies (GWASs) are the most widely used method to identify genetic risk loci associated with orofacial clefts (OFC). However, despite the increasing size of cohort, GWASs are still insufficient to detect all the heritability, suggesting there are more associations under the current stringent statistical threshold. In this study, we obtained an integrated epigenomic dataset based on the chromatin conformation of a human oral epithelial cell line (HIOEC) using RNA-seq, ATAC-seq, H3K27ac ChIP-seq, and DLO Hi-C. Presumably, this epigenomic dataset could reveal the missing functional variants located in the oral epithelial cell active enhancers/promoters along with their risk target genes, despite relatively less-stringent statistical association with OFC. Taken a non-syndromic cleft palate only (NSCPO) GWAS data of the Chinese Han population as an example, 3664 SNPs that cannot reach the strict significance threshold were subjected to this functional identification pipeline. In total, 254 potential risk SNPs residing in active cis-regulatory elements interacting with 1 718 promoters of oral epithelium-expressed genes were screened. Gapped k-mer machine learning based on enhancers interacting with epithelium-expressed genes along with in vivo and in vitro reporter assays were employed as functional validation. Among all the potential SNPs, we chose and confirmed that the risk alleles of rs560789 and rs174570 reduced the epithelial-specific enhancer activity by preventing the binding of transcription factors related to epithelial development. In summary, we established chromatin conformation datasets of human oral epithelial cells and provided a framework for testing and understanding how regulatory variants impart risk for clefts.

mafba and mafbb differentially regulate lymphatic endothelial cell migration in topographically distinct manners

Lymphangiogenesis, formation of lymphatic vessels from pre-existing vessels, is a dynamic process that requires cell migration. Regardless of location, migrating lymphatic endothelial cell (LEC) progenitors probe their surroundings to form the lymphatic network. Lymphatic-development regulation requires the transcription factor MAFB in different species. Zebrafish Mafba, expressed in LEC progenitors, is essential for their migration in the trunk. However, the transcriptional mechanism that orchestrates LEC migration in different lymphatic endothelial beds remains elusive. Here, we uncover topographically different requirements of the two paralogs, Mafba and Mafbb, for LEC migration. Both mafba and mafbb are necessary for facial lymphatic development, but mafbb is dispensable for trunk lymphatic development. On the molecular level, we demonstrate a regulatory network where Vegfc-Vegfd-SoxF-Mafba-Mafbb is essential in facial lymphangiogenesis. We identify that mafba and mafbb tune the directionality of LEC migration and vessel morphogenesis that is ultimately necessary for lymphatic function.

Harnessing changes in open chromatin determined by ATAC-seq to generate insulin-responsive reporter constructs

BACKGROUND

Gene regulation is critical for proper cellular function. Next-generation sequencing technology has revealed the presence of regulatory networks that regulate gene expression and essential cellular functions. Studies investigating the epigenome have begun to uncover the complex mechanisms regulating transcription. Assay for transposase-accessible chromatin by sequencing (ATAC-seq) is quickly becoming the assay of choice for many epigenomic investigations. However, whether intervention-mediated changes in accessible chromatin determined by ATAC-seq can be harnessed to generate intervention-inducible reporter constructs has not been systematically assayed.

RESULTS

We used the insulin signaling pathway as a model to investigate chromatin regions and gene expression changes using ATAC- and RNA-seq in insulin-treated Drosophila S2 cells. We found correlations between ATAC- and RNA-seq data, especially when stratifying differentially-accessible chromatin regions by annotated feature type. In particular, our data demonstrated a weak but significant correlation between chromatin regions annotated to enhancers (1-2 kb from the transcription start site) and downstream gene expression. We cloned candidate enhancer regions upstream of luciferase and demonstrate insulin-inducibility of several of these reporters.

CONCLUSIONS

Insulin-induced chromatin accessibility determined by ATAC-seq reveals enhancer regions that drive insulin-inducible reporter gene expression.

Investigating chromatin accessibility during development and differentiation by ATAC-sequencing to guide the identification of cis-regulatory elements

Mapping accessible chromatin across time scales can give insights into its dynamic nature, for example during cellular differentiation and tissue or organism development. Analysis of such data can be utilised to identify functional cis-regulatory elements (CRE) and transcription factor binding sites and, when combined with transcriptomics, can reveal gene regulatory networks (GRNs) of expressed genes. Chromatin accessibility mapping is a powerful approach and can be performed using ATAC-sequencing (ATAC-seq), whereby Tn5 transposase inserts sequencing adaptors into genomic DNA to identify differentially accessible regions of chromatin in different cell populations. It requires low sample input and can be performed and analysed relatively quickly compared with other methods. The data generated from ATAC-seq, along with other genomic approaches, can help uncover chromatin packaging and potential cis-regulatory elements that may be responsible for gene expression. Here, we describe the ATAC-seq approach and give examples from mainly vertebrate embryonic development, where such datasets have identified the highly dynamic nature of chromatin, with differing landscapes between cellular precursors for different lineages.

Epigenetic Regulation of Endothelial Cell Lineages During Zebrafish Development-New Insights From Technical Advances

Epigenetic regulation is integral in orchestrating the spatiotemporal regulation of gene expression which underlies tissue development. The emergence of new tools to assess genome-wide epigenetic modifications has enabled significant advances in the field of vascular biology in zebrafish. Zebrafish represents a powerful model to investigate the activity of cis-regulatory elements in vivo by combining technologies such as ATAC-seq, ChIP-seq and CUT&Tag with the generation of transgenic lines and live imaging to validate the activity of these regulatory elements. Recently, this approach led to the identification and characterization of key enhancers of important vascular genes, such as gata2a, notch1b and dll4. In this review we will discuss how the latest technologies in epigenetics are being used in the zebrafish to determine chromatin states and assess the function of the cis-regulatory sequences that shape the zebrafish vascular network.

Using Light-Sheet Microscopy to Study Spontaneous Activity in the Developing Lateral-Line System

Hair cells are the sensory receptors in the auditory and vestibular systems of all vertebrates, and in the lateral-line system of aquatic vertebrates. The purpose of this work is to explore the zebrafish lateral-line system as a model to study and understand spontaneous activity in vivo. Our work applies genetically encoded calcium indicators along with light-sheet fluorescence microscopy to visualize spontaneous calcium activity in the developing lateral-line system. Consistent with our previous work, we show that spontaneous calcium activity is present in developing lateral-line hair cells. We now show that supporting cells that surround hair cells, and cholinergic efferent terminals that directly contact hair cells are also spontaneously active. Using two-color functional imaging we demonstrate that spontaneous activity in hair cells does not correlate with activity in either supporting cells or cholinergic terminals. We find that during lateral-line development, hair cells autonomously generate spontaneous events. Using localized calcium indicators, we show that within hair cells, spontaneous calcium activity occurs in two distinct domains—the mechanosensory bundle and the presynapse. Further, spontaneous activity in the mechanosensory bundle ultimately drives spontaneous calcium influx at the presynapse. Comprehensively, our results indicate that in developing lateral-line hair cells, autonomously generated spontaneous activity originates with spontaneous mechanosensory events.

Optimized assay for transposase-accessible chromatin by sequencing (ATAC-seq) library preparation from adult Drosophila melanogaster neurons

Assay for transposase-accessible chromatin by sequencing (ATAC-seq) is rapidly becoming the assay of choice to investigate chromatin-mediated gene regulation, largely because of low input requirements, a fast workflow, and the ability to interrogate the entire genome in an untargeted manner. Many studies using ATAC-seq use mammalian or human-derived tissues, and established protocols work well in these systems. However, ATAC-seq is not yet widely used in Drosophila. Vinegar flies present several advantages over mammalian systems that make them an excellent model for ATAC-seq studies, including abundant genetic tools that allow straightforward targeting, transgene expression, and genetic manipulation that are not available in mammalian models. Because current ATAC-seq protocols are not optimized to use flies, we developed an optimized workflow that accounts for several complicating factors present in Drosophila. We examined parameters affecting nuclei isolation, including input size, freezing time, washing, and possible confounds from retinal pigments. Then, we optimized the enzymatic steps of library construction to account for the smaller Drosophila genome size. Finally, we used our optimized protocol to generate ATAC-seq libraries that meet ENCODE quality metrics. Our optimized protocol enables extensive ATAC-seq experiments in Drosophila, thereby leveraging the advantages of this powerful model system to understand chromatin-mediated gene regulation.

Shear stress switches the association of endothelial enhancers from ETV/ETS to KLF transcription factor binding sites

Endothelial cells (ECs) lining blood vessels are exposed to mechanical forces, such as shear stress. These forces control many aspects of EC biology, including vascular tone, cell migration and proliferation. Despite a good understanding of the genes responding to shear stress, our insight into the transcriptional regulation of these genes is much more limited. Here, we set out to study alterations in the chromatin landscape of human umbilical vein endothelial cells (HUVEC) exposed to laminar shear stress. To do so, we performed ChIP-Seq for H3K27 acetylation, indicative of active enhancer elements and ATAC-Seq to mark regions of open chromatin in addition to RNA-Seq on HUVEC exposed to 6 h of laminar shear stress. Our results show a correlation of gained and lost enhancers with up and downregulated genes, respectively. DNA motif analysis revealed an over-representation of KLF transcription factor (TF) binding sites in gained enhancers, while lost enhancers contained more ETV/ETS motifs. We validated a subset of flow responsive enhancers using luciferase-based reporter constructs and CRISPR-Cas9 mediated genome editing. Lastly, we characterized the shear stress response in ECs of zebrafish embryos using RNA-Seq. Our results lay the groundwork for the exploration of shear stress responsive elements in controlling EC biology.

Environmental adaptation in fish induced changes in the regulatory region of fatty acid elongase gene, elovl5, involved in long-chain polyunsaturated fatty acid biosynthesis

Fish are the main source of long-chain polyunsaturated fatty acids (LC-PUFA) for human consumption. In the process of evolution via natural selection, adaptation to distinct environments has likely driven changes in the endogenous capacity for LC-PUFA biosynthesis between marine and freshwater fishes. However, the molecular mechanisms underlying adaptive changes in this metabolic pathway are poorly understood. Here, we compared the transcriptional regulation of elongation of very long chain fatty acids protein 5 (Elovl5), which is one of the critical enzymes in LC-PUFA biosynthesis pathway, in marine large yellow croaker (Larimichthys crocea) and freshwater rainbow trout (Oncorhynchus mykiss). Comparative transcriptomic and absolute mRNA quantification analyses revealed that the expression of elovl5 in rainbow trout was markedly higher than that in large yellow croaker. Correspondingly, the number of chromatin accessible areas in the regulatory region of elovl5 in rainbow trout was higher than in large yellow croaker, which revealed that chromatin accessibility in the regulatory region of elovl5 in rainbow trout was higher. Furthermore, the differences in sequence and activity of the elovl5 promoter were observed between rainbow trout and large yellow croaker, and transcription factors including CCAAT/enhancer-binding protein β (CEBPβ), GATA binding protein 3 (GATA3) and upstream stimulatory factor 2 (USF2) displayed different regulatory roles on elovl5 expression between the two species. We propose that changes in the gene regulatory region driven by natural selection likely play a key role in differences in elovl5 expression and the activity of Elovl5, which may influence the LC-PUFA biosynthesis capacities of rainbow trout and large yellow croaker. These findings may also provide opportunities to improve the quality of aquatic products and, consequently, human health.

Finding and Verifying Enhancers for Endothelial-Expressed Genes

Identification and analysis of enhancers for endothelial-expressed genes can provide crucial information regarding their upstream transcriptional regulators. However, enhancer identification can be challenging, particularly for people with limited access or experience of bioinformatics, and transgenic analysis of enhancer activity patterns can be prohibitively expensive. Here we describe how to use publicly available datasets displayed on the UCSC Genome Browser to identify putative endothelial enhancers for mammalian genes. Furthermore, we detail how to utilize mosaic Tol2-mediated transgenesis in zebrafish to verify whether a putative enhancer is capable of directing endothelial-specific patterns of gene expression.

Integrated molecular analysis identifies a conserved pericyte gene signature in zebrafish

Pericytes reside in capillary beds where they share a basement membrane with endothelial cells and regulate their function. However, little is known about embryonic pericyte development, in part, due to lack of specific molecular markers and genetic tools. Here, we applied single cell RNA-sequencing (scRNA-seq) of platelet derived growth factor beta (pdgfrb)-positive cells to molecularly characterize pericytes in zebrafish larvae. scRNA-seq revealed zebrafish cells expressing mouse pericyte gene orthologs, and comparison with bulk RNA-seq from wild-type and pdgfrb mutant larvae further refined a pericyte gene set. Subsequent integration with mouse pericyte scRNA-seq profiles revealed a core set of conserved pericyte genes. Using transgenic reporter lines, we validated pericyte expression of two genes identified in our analysis: NDUFA4 mitochondrial complex associated like 2a (ndufa4l2a), and potassium voltage-gated channel, Isk-related family, member 4 (kcne4). Both reporter lines exhibited pericyte expression in multiple anatomical locations, and kcne4 was also detected in a subset of vascular smooth muscle cells. Thus, our integrated molecular analysis revealed a molecular profile for zebrafish pericytes and allowed us to develop new tools to observe these cells in vivo.

Identification of downstream effectors of retinoic acid specifying the zebrafish pancreas by integrative genomics

Retinoic acid (RA) is a key signal for the specification of the pancreas. Still, the gene regulatory cascade triggered by RA in the endoderm remains poorly characterized. In this study, we investigated this regulatory network in zebrafish by combining RNA-seq, RAR ChIP-seq and ATAC-seq assays. By analysing the effect of RA and of the RA receptor (RAR) inverse-agonist BMS493 on the transcriptome and on the chromatin accessibility of endodermal cells, we identified a large set of genes and regulatory regions regulated by RA signalling. RAR ChIP-seq further defined the direct RAR target genes in zebrafish, including hox genes as well as several pancreatic regulators like mnx1, insm1b, hnf1ba and gata6. Comparison of zebrafish and murine RAR ChIP-seq data highlighted the conserved direct target genes and revealed that some RAR sites are under strong evolutionary constraints. Among them, a novel highly conserved RAR-induced enhancer was identified downstream of the HoxB locus and driving expression in the nervous system and in the gut in a RA-dependent manner. Finally, ATAC-seq data unveiled the role of the RAR-direct targets Hnf1ba and Gata6 in opening chromatin at many regulatory loci upon RA treatment.

Prmt5 promotes vascular morphogenesis independently of its methyltransferase activity

During development, the vertebrate vasculature undergoes major growth and remodeling. While the transcriptional cascade underlying blood vessel formation starts to be better characterized, little is known concerning the role and mode of action of epigenetic enzymes during this process. Here, we explored the role of the Protein Arginine Methyl Transferase Prmt5 in blood vessel formation as well as hematopoiesis using zebrafish as a model system. Through the combination of different prmt5 loss-of-function approaches we highlighted a key role of Prmt5 in both processes. Notably, we showed that Prmt5 promotes vascular morphogenesis through the transcriptional control of ETS transcription factors and adhesion proteins in endothelial cells. Interestingly, using a catalytic dead mutant of Prmt5 and a specific drug inhibitor, we found that while Prmt5 methyltransferase activity was required for blood cell formation, it was dispensable for vessel formation. Analyses of chromatin architecture impact on reporter genes expression and chromatin immunoprecipitation experiments led us to propose that Prmt5 regulates transcription by acting as a scaffold protein that facilitates chromatin looping to promote vascular morphogenesis.

Blood vessel formation is an essential developmental process required for the survival of all vertebrates. The vascular anatomy and the mechanisms involved in vessel formation are highly conserved among vertebrates. Hence, we used zebrafish as a model, to decipher the role and the mode of action of Prmt5, an enzyme known to regulate gene expression, in vascular morphogenesis and in blood cell formation in vivo. Using different approaches, we highlighted a key role of Prmt5 during both processes. However, we found that while blood cell formation required Prmt5 enzymatic activity, vascular morphogenesis was independent on its activity. Prmt5 has been proposed as a therapeutic target in many diseases, including cancer. Yet, we show here that Prmt5 acts at least in part independently of its methyltransferase activity to regulate vascular morphogenesis. By shedding light on a mechanism of action of Prmt5 that will be insensitive to enzymatic inhibitors, our data calls forth the design of alternative drugs. In addition, this non-canonical function of Prmt5 may have a more pervasive role than previously thought in physiological conditions, i.e. during development, but also in pathological situations such as in tumor angiogenesis and certainly deserves more attention in the future.

Functional in vivo characterization of sox10 enhancers in neural crest and melanoma development

The role of a neural crest developmental transcriptional program, which critically involves Sox10 upregulation, is a key conserved aspect of melanoma initiation in both humans and zebrafish, yet transcriptional regulation of sox10 expression is incompletely understood. Here we used ATAC-Seq analysis of multiple zebrafish melanoma tumors to identify recurrently open chromatin domains as putative melanoma-specific sox10 enhancers. Screening in vivo with EGFP reporter constructs revealed 9 of 11 putative sox10 enhancers with embryonic activity in zebrafish. Focusing on the most active enhancer region in melanoma, we identified a region 23 kilobases upstream of sox10, termed peak5, that drives EGFP reporter expression in a subset of neural crest cells, Kolmer-Agduhr neurons, and early melanoma patches and tumors with high specificity. A ~200 base pair region, conserved in Cyprinidae, within peak5 is required for transgenic reporter activity in neural crest and melanoma. This region contains dimeric SoxE/Sox10 dimeric binding sites essential for peak5 neural crest and melanoma activity. We show that deletion of the endogenous peak5 conserved genomic locus decreases embryonic sox10 expression and disrupts adult stripe patterning in our melanoma model background. Our work demonstrates the power of linking developmental and cancer models to better understand neural crest identity in melanoma.

Zebrafish mafbb Mutants Display Osteoclast Over-Activation and Bone Deformity Resembling Osteolysis in MCTO Patients

Multicentric carpotarsal osteolysis (MCTO) is a rare skeletal dysplasia with osteolysis at the carpal and tarsal bones. Heterozygous missense mutations in the transcription factor MAFB are found in patients with MCTO. MAFB is reported to negatively regulate osteoclastogenesis in vitro. However, the in vivo function of MAFB and its relation to MCTO remains unknown. In this study, we generated zebrafish MAFB homolog mafbb mutant utilizing CRISPR/Cas9 technology. Mafbb deficient zebrafish demonstrated enhanced osteoclast cell differentiation and abnormal cartilage and bone development resembling MCTO patients. It is known that osteoclasts are hematopoietic cells derived from macrophages. Loss of mafbb caused selective expansion of definitive macrophages and myeloid cells, supporting that mafbb restricts myeloid differentiation in vivo. We also demonstrate that MAFB MCTO mutations failed to rescue the defective osteoclastogenesis in mafbb^-/- embryos, but did not affect osteoclast cells in wild type embryos. The mechanism of MCTO mutations is likely haploinsufficiency. Zebrafish mafbb mutant provides a useful model to study the function of MAFB in osteoclastogenesis and the related MCTO disease.

New tools for 'ZEBRA-FISHING'

Zebrafish has been established as a classical model for developmental studies, yet in the past years, with the explosion of novel technological methods, the use of zebrafish as a model has expanded. One of the prominent fields that took advantage of zebrafish as a model organism early on is hematopoiesis, the process of blood cell generation from hematopoietic stem and progenitor cells (HSPCs). In zebrafish, HSPCs are born early during development in the aorta-gonad-mesonephros region and then translocate to the caudal hematopoietic tissue, where they expand and finally take residence in the kidney marrow. This journey is tightly regulated at multiple levels from extracellular signals to chromatin. In order to delineate the mechanistic underpinnings of this process, next-generation sequencing techniques could be an important ally. Here, we describe genome-wide approaches that have been undertaken to delineate zebrafish hematopoiesis.

Retrograde Mitochondrial Transport Is Essential for Organelle Distribution and Health in Zebrafish Neurons

In neurons, mitochondria are transported by molecular motors throughout the cell to form and maintain functional neural connections. These organelles have many critical functions in neurons and are of high interest as their dysfunction is associated with disease. While the mechanics and impact of anterograde mitochondrial movement toward axon terminals are beginning to be understood, the frequency and function of retrograde (cell body directed) mitochondrial transport in neurons are still largely unexplored. While existing evidence indicates that some mitochondria are retrogradely transported for degradation in the cell body, the precise impact of disrupting retrograde transport on the organelles and the axon was unknown. Using long-term, in vivo imaging, we examined mitochondrial motility in zebrafish sensory and motor axons. We show that retrograde transport of mitochondria from axon terminals allows replacement of the axon terminal population within a day. By tracking these organelles, we show that not all mitochondria that leave the axon terminal are degraded; rather, they persist over several days. Disrupting retrograde mitochondrial flux in neurons leads to accumulation of aged organelles in axon terminals and loss of cell body mitochondria. Assays of neural circuit activity demonstrated that disrupting mitochondrial transport and function has no effect on sensory axon terminal activity but does negatively impact motor neuron axons. Taken together, our work supports a previously unappreciated role for retrograde mitochondrial transport in the maintenance of a homeostatic distribution of mitochondria in neurons and illustrates the downstream effects of disrupting this process on sensory and motor circuits.

SIGNIFICANCE STATEMENT Disrupted mitochondrial transport has been linked to neurodegenerative disease. Retrograde transport of this organelle has been implicated in turnover of aged organelles through lysosomal degradation in the cell body. Consistent with this, we provide evidence that retrograde mitochondrial transport is important for removing aged organelles from axons; however, we show that these organelles are not solely degraded, rather they persist in neurons for days. Disrupting retrograde mitochondrial transport impacts the homeostatic distribution of mitochondria throughout the neuron and the function of motor, but not sensory, axon synapses. Together, our work shows the conserved reliance on retrograde mitochondrial transport for maintaining a healthy mitochondrial pool in neurons and illustrates the disparate effects of disrupting this process on sensory versus motor circuits.

A map of cis-regulatory elements and 3D genome structures in zebrafish

The zebrafish (Danio rerio) has been widely used in the study of human disease and development, and about 70% of the protein-coding genes are conserved between the two species¹. However, studies in zebrafish remain constrained by the sparse annotation of functional control elements in the zebrafish genome. Here we performed RNA sequencing, assay for transposase-accessible chromatin using sequencing (ATAC-seq), chromatin immunoprecipitation with sequencing, whole-genome bisulfite sequencing, and chromosome conformation capture (Hi-C) experiments in up to eleven adult and two embryonic tissues to generate a comprehensive map of transcriptomes, cis-regulatory elements, heterochromatin, methylomes and 3D genome organization in the zebrafish Tübingen reference strain. A comparison of zebrafish, human and mouse regulatory elements enabled the identification of both evolutionarily conserved and species-specific regulatory sequences and networks. We observed enrichment of evolutionary breakpoints at topologically associating domain boundaries, which were correlated with strong histone H3 lysine 4 trimethylation (H3K4me3) and CCCTC-binding factor (CTCF) signals. We performed single-cell ATAC-seq in zebrafish brain, which delineated 25 different clusters of cell types. By combining long-read DNA sequencing and Hi-C, we assembled the sex-determining chromosome 4 de novo. Overall, our work provides an additional epigenomic anchor for the functional annotation of vertebrate genomes and the study of evolutionarily conserved elements of 3D genome organization.

Distinctive regulatory architectures of germline-active and somatic genes in C. elegans

RNA profiling has provided increasingly detailed knowledge of gene expression patterns, yet the different regulatory architectures that drive them are not well understood. To address this, we profiled and compared transcriptional and regulatory element activities across five tissues of Caenorhabditis elegans, covering ∼90% of cells. We find that the majority of promoters and enhancers have tissue-specific accessibility, and we discover regulatory grammars associated with ubiquitous, germline, and somatic tissue–specific gene expression patterns. In addition, we find that germline-active and soma-specific promoters have distinct features. Germline-active promoters have well-positioned +1 and −1 nucleosomes associated with a periodic 10-bp WW signal (W = A/T). Somatic tissue–specific promoters lack positioned nucleosomes and this signal, have wide nucleosome-depleted regions, and are more enriched for core promoter elements, which largely differ between tissues. We observe the 10-bp periodic WW signal at ubiquitous promoters in other animals, suggesting it is an ancient conserved signal. Our results show fundamental differences in regulatory architectures of germline and somatic tissue–specific genes, uncover regulatory rules for generating diverse gene expression patterns, and provide a tissue-specific resource for future studies.

An improved zebrafish transcriptome annotation for sensitive and comprehensive detection of cell type-specific genes

The zebrafish is ideal for studying embryogenesis and is increasingly applied to model human disease. In these contexts, RNA-sequencing (RNA-seq) provides mechanistic insights by identifying transcriptome changes between experimental conditions. Application of RNA-seq relies on accurate transcript annotation for a genome of interest. Here, we find discrepancies in analysis from RNA-seq datasets quantified using Ensembl and RefSeq zebrafish annotations. These issues were due, in part, to variably annotated 3' untranslated regions and thousands of gene models missing from each annotation. Since these discrepancies could compromise downstream analyses and biological reproducibility, we built a more comprehensive zebrafish transcriptome annotation that addresses these deficiencies. Our annotation improves detection of cell type-specific genes in both bulk and single cell RNA-seq datasets, where it also improves resolution of cell clustering. Thus, we demonstrate that our new transcriptome annotation can outperform existing annotations, providing an important resource for zebrafish researchers.

Zebrafish models of skeletal dysplasia induced by cholesterol biosynthesis deficiency

Human disorders of the post-squalene cholesterol biosynthesis pathway frequently result in skeletal abnormalities, yet our understanding of the mechanisms involved is limited. In a forward-genetic approach, we have found that a late-onset skeletal mutant, named kolibernu7 , is the result of a cis-acting regulatory mutation leading to loss of methylsterol monooxygenase 1 (msmo1) expression within pre-hypertrophic chondrocytes. Generated msmo1nu81 knockdown mutation resulted in lethality at larval stage. We demonstrated that this is a result of both cholesterol deprivation and sterol intermediate accumulation by creating a mutation eliminating activity of Lanosterol synthase (Lss). Our results indicate that double lssnu60;msmo1nu81 and single lssnu60 mutants survive significantly longer than msmo1nu81 homozygotes. Liver-specific restoration of either Msmo1 or Lss in corresponding mutant backgrounds suppresses larval lethality. Rescued mutants develop dramatic skeletal abnormalities, with a loss of Msmo1 activity resulting in a more-severe patterning defect of a near-complete loss of hypertrophic chondrocytes marked by col10a1a expression. Our analysis suggests that hypertrophic chondrocytes depend on endogenous cholesterol synthesis, and blocking C4 demethylation exacerbates the cholesterol deficiency phenotype. Our findings offer new insight into the genetic control of bone development and provide new zebrafish models for human disorders of the cholesterol biosynthesis pathway.

Analysis of zebrafish periderm enhancers facilitates identification of a regulatory variant near human KRT8/18

Genome-wide association studies for non-syndromic orofacial clefting (OFC) have identified single nucleotide polymorphisms (SNPs) at loci where the presumed risk-relevant gene is expressed in oral periderm. The functional subsets of such SNPs are difficult to predict because the sequence underpinnings of periderm enhancers are unknown. We applied ATAC-seq to models of human palate periderm, including zebrafish periderm, mouse embryonic palate epithelia, and a human oral epithelium cell line, and to complementary mesenchymal cell types. We identified sets of enhancers specific to the epithelial cells and trained gapped-kmer support-vector-machine classifiers on these sets. We used the classifiers to predict the effects of 14 OFC-associated SNPs at 12q13 near KRT18. All the classifiers picked the same SNP as having the strongest effect, but the significance was highest with the classifier trained on zebrafish periderm. Reporter and deletion analyses support this SNP as lying within a periderm enhancer regulating KRT18/KRT8 expression.

Comprehensive epigenome characterization reveals diverse transcriptional regulation across human vascular endothelial cells

Background

Endothelial cells (ECs) make up the innermost layer throughout the entire vasculature. Their phenotypes and physiological functions are initially regulated by developmental signals and extracellular stimuli. The underlying molecular mechanisms responsible for the diverse phenotypes of ECs from different organs are not well understood.

Results

To characterize the transcriptomic and epigenomic landscape in the vascular system, we cataloged gene expression and active histone marks in nine types of human ECs (generating 148 genome-wide datasets) and carried out a comprehensive analysis with chromatin interaction data. We developed a robust procedure for comparative epigenome analysis that circumvents variations at the level of the individual and technical noise derived from sample preparation under various conditions. Through this approach, we identified 3765 EC-specific enhancers, some of which were associated with disease-associated genetic variations. We also identified various candidate marker genes for each EC type. We found that the nine EC types can be divided into two subgroups, corresponding to those with upper-body origins and lower-body origins, based on their epigenomic landscape. Epigenomic variations were highly correlated with gene expression patterns, but also provided unique information. Most of the deferentially expressed genes and enhancers were cooperatively enriched in more than one EC type, suggesting that the distinct combinations of multiple genes play key roles in the diverse phenotypes across EC types. Notably, many homeobox genes were differentially expressed across EC types, and their expression was correlated with the relative position of each organ in the body. This reflects the developmental origins of ECs and their roles in angiogenesis, vasculogenesis and wound healing.

Conclusions

This comprehensive analysis of epigenome characterization of EC types reveals diverse transcriptional regulation across human vascular systems. These datasets provide a valuable resource for understanding the vascular system and associated diseases.

Valves Are a Conserved Feature of the Zebrafish Lymphatic System

The lymphatic system comprises blind-ended tubes that collect interstitial fluid and return it to the circulatory system. In mammals, unidirectional lymphatic flow is driven by muscle contraction working in conjunction with valves. Accordingly, defective lymphatic valve morphogenesis results in backflow leading to edema. In fish species, studies dating to the 18^th century failed to identify lymphatic valves, a precedent that currently persists, raising the question of whether the zebrafish could be used to study the development of these structures. Here, we provide functional and morphological evidence of valves in the zebrafish lymphatic system. Electron microscopy revealed valve ultrastructure similar to mammals, while live imaging using transgenic lines identified the developmental origins of lymphatic valve progenitors. Zebrafish embryos bearing mutations in genes required for mammalian valve morphogenesis show defective lymphatic valve formation and edema. Together, our observations provide a foundation from which to further investigate lymphatic valve formation in zebrafish.

Stem cell safe harbor: the hematopoietic stem cell niche in zebrafish

Each stem cell resides in a highly specialized anatomic location known as the niche that protects and regulates stem cell function. The importance of the niche in hematopoiesis has long been appreciated in transplantation, but without methods to observe activity in vivo, the components and mechanisms of the hematopoietic niche have remained incompletely understood. Zebrafish have emerged over the past few decades as an answer to this. Use of zebrafish to study the hematopoietic niche has enabled discovery of novel cell-cell interactions, as well as chemical and genetic regulators of hematopoietic stem cells. Mastery of niche components may improve therapeutic efforts to direct differentiation of hematopoietic stem cells from pluripotent cells, sustain stem cells in culture, or improve stem cell transplant.

Broad Heterochromatic Domains Open in Gonocyte Development Prior to De Novo DNA Methylation

Facultative heterochromatin forms and reorganizes in response to external stimuli. However, how the initial establishment of such a chromatin state is regulated in cell-cycle-arrested cells remains unexplored. Mouse gonocytes are arrested male germ cells, at which stage the genome-wide DNA methylome forms. Here, we discovered transiently accessible heterochromatin domains of several megabases in size in gonocytes and named them differentially accessible domains (DADs). Open DADs formed in gene desert and gene cluster regions, primarily at transposons, with the reprogramming of histone marks, suggesting DADs as facultative heterochromatin. De novo DNA methylation took place with two waves in gonocytes: the first region specific and the second genome-wide. DADs were resistant to the first wave and their opening preceded the second wave. In addition, the higher-order chromosome architecture was reorganized with less defined chromosome compartments in gonocytes. These findings suggest that multiple layers of chromatin reprogramming facilitate de novo DNA methylation.

Detect accessible chromatin using ATAC-sequencing, from principle to applications

Background

Chromatin accessibility is crucial for gene expression regulation in specific cells and in multiple biological processes. Assay for Transposase Accessible Chromatin with high-throughput sequencing (ATAC-seq) is an effective way to reveal chromatin accessibility at a genome-wide level. Through ATAC-seq, produced reads from a small number of cells reflect accessible regions that correspond to nucleosome positioning and transcription factor binding sites, due to probing hyperactive Tn5 transposase to DNA sequence.

Conclusion

In this review, we summarize both principle and features of ATAC-seq, highlight its applications in basic and clinical research. ATAC-seq has generated comprehensive chromatin accessible maps, and is becoming a powerful tool to understand dynamic gene expression regulation in stem cells, early embryos and tumors.

foxc1 is required for embryonic head vascular smooth muscle differentiation in zebrafish

Vascular smooth muscle of the head derives from neural crest, but developmental mechanisms and early transcriptional drivers of the vSMC lineage are not well characterized. We find that in early development, the transcription factor foxc1b is expressed in mesenchymal cells that associate with the vascular endothelium. Using timelapse imaging, we observe that foxc1b expressing mesenchymal cells differentiate into acta2 expressing vascular mural cells. We show that in zebrafish, while foxc1b is co-expressed in acta2 positive smooth muscle cells that associate with large diameter vessels, it is not co-expressed in capillaries where pdgfrβ positive pericytes are located. In addition to being an early marker of the lineage, foxc1 is essential for vSMC differentiation; we find that foxc1 loss of function mutants have defective vSMC differentiation and that early genetic ablation of foxc1b or acta2 expressing populations blocks vSMC differentiation. Furthermore, foxc1 is expressed upstream of acta2 and is required for acta2 expression in vSMCs. Using RNA-Seq we determine an enriched intersectional gene expression profile using dual expression of foxc1b and acta2 to identify novel vSMC markers. Taken together, our data suggests that foxc1 is a marker of vSMCs and plays a critical functional role in promoting their differentiation.

Packaging development: how chromatin controls transcription in zebrafish embryogenesis

How developmental gene expression is activated, co-ordinated and maintained is one of the biggest questions in developmental biology. While transcription factors lead the way in directing developmental gene expression, their accessibility to the correct repertoire of genes can depend on other factors such as DNA methylation, the presence of particular histone variants and post-translational modifications of histones. Collectively, factors that modify DNA or affect its packaging and accessibility contribute to a chromatin landscape that helps to control the timely expression of developmental genes. Zebrafish, perhaps better known for their strength as a model of embryology and organogenesis during development, are coming to the fore as a powerful model for interpreting the role played by chromatin in gene expression. Several recent advances have shown that zebrafish exhibit both similarities and differences to other models (and humans) in the way that they employ chromatin mechanisms of gene regulation. Here, I review how chromatin influences developmental transcriptional programmes during early zebrafish development, patterning and organogenesis. Lastly, I briefly highlight the importance of zebrafish chromatin research towards the understanding of human disease and transgenerational inheritance.

PAN-INTACT enables direct isolation of lineage-specific nuclei from fibrous tissues

Recent studies have highlighted the extraordinary cell type diversity that exists within mammalian organs, yet the molecular drivers of such heterogeneity remain elusive. To address this issue, much attention has been focused on profiling the transcriptome and epigenome of individual cell types. However, standard cell type isolation methods based on surface or fluorescent markers remain problematic for cells residing within organs with significant connective tissue. Since the nucleus contains both genomic and transcriptomic information, the isolation of nuclei tagged in specific cell types (INTACT) method provides an attractive solution. Although INTACT has been successfully applied to plants, flies, zebrafish, frogs, and mouse brain and adipose tissue, broad use across mammalian organs remains challenging. Here we describe the PAN-INTACT method, which can be used to isolate cell type specific nuclei from fibrous mouse organs, which are particularly problematic. As a proof-of-concept, we demonstrate successful isolation of cell type-specific nuclei from the mouse heart, which contains substantial connective tissue and harbors multiple cell types, including cardiomyocytes, fibroblasts, endothelial cells, and epicardial cells. Compared to established techniques, PAN-INTACT allows more rapid isolation of cardiac nuclei to facilitate downstream applications. We show cell type-specific isolation of nuclei from the hearts of Nkx2-5^Cre/+; R26^{Sun1-2xsf-GFP-6xmyc/+} mice, which we confirm by expression of lineage markers. Furthermore, we perform Assay for Transposase Accessible Chromatin (ATAC)-Seq to provide high-fidelity chromatin accessibility maps of Nkx2-5⁺ nuclei. To extend the applicability of PAN-INTACT, we also demonstrate successful isolation of Wt1⁺ podocytes from adult kidney. Taken together, our data suggest that PAN-INTACT is broadly applicable for profiling the transcriptional and epigenetic landscape of specific cell types. Thus, we envision that our method can be used to systematically probe mechanistic details of cell type-specific functions within individual organs of intact mice.

Recent insights into vascular development from studies in zebrafish

PURPOSE OF REVIEW

Zebrafish has provided a powerful platform to study vascular biology over the past 25 years, owing to their distinct advantages for imaging and genetic manipulation. In this review, we summarize recent progress in vascular biology with particular emphasis on vascular development in zebrafish.

RECENT FINDINGS

The advent of transcription activator-like effector nuclease and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 genome-editing technologies has dramatically facilitated reverse genetic approaches in zebrafish, as in other models. Here, we highlight recent studies on vascular development in zebrafish which mainly employed forward or reverse genetics combined with high-resolution imaging. These studies have advanced our understanding of diverse areas in vascular biology, including transcriptional regulation of endothelial cell differentiation, endothelial cell signaling during angiogenesis and lymphangiogenesis, vascular bed-specific developmental mechanisms, and perivascular cell recruitment.

SUMMARY

The unique attributes of the zebrafish model have allowed critical cellular and molecular insights into fundamental mechanisms of vascular development. Knowledge acquired through recent zebrafish work further advances our understanding of basic mechanisms underlying vascular morphogenesis, maintenance, and homeostasis. Ultimately, insights provided by the zebrafish model will help to understand the genetic, cellular, and molecular underpinnings of human vascular malformations and diseases.

Genome-wide analysis of chromatin accessibility using ATAC-seq

Programs of gene transcription are controlled by cis-acting DNA elements, including enhancers, silencers, and promoters. Local accessibility of chromatin has proven to be a highly informative structural feature for identifying such regulatory elements, which tend to be relatively open due to their interactions with proteins. Recently, ATAC-seq (assay for transposase-accessible chromatin using sequencing) has emerged as one of the most powerful approaches for genome-wide chromatin accessibility profiling. This method assesses DNA accessibility using hyperactive Tn5 transposase, which simultaneously cuts DNA and inserts sequencing adaptors, preferentially in regions of open chromatin. ATAC-seq is a relatively simple procedure which can be applied to only a few thousand cells. It is well-suited to developing embryos of sea urchins and other echinoderms, which are a prominent experimental model for understanding the genomic control of animal development. In this chapter, we present a protocol for applying ATAC-seq to embryonic cells of sea urchins.

Bloody Zebrafish: Novel Methods in Normal and Malignant Hematopoiesis

Hematopoiesis is an optimal system for studying stem cell maintenance and lineage differentiation under physiological and pathological conditions. In vertebrate organisms, billions of differentiated hematopoietic cells need to be continuously produced to replenish the blood cell pool. Disruptions in this process have immediate consequences for oxygen transport, responses against pathogens, maintenance of hemostasis and vascular integrity. Zebrafish is a widely used and well-established model for studying the hematopoietic system. Several new hematopoietic regulators were identified in genetic and chemical screens using the zebrafish model. Moreover, zebrafish enables in vivo imaging of hematopoietic stem cell generation and differentiation during embryogenesis, and adulthood. Finally, zebrafish has been used to model hematopoietic diseases. Recent technological advances in single-cell transcriptome analysis, epigenetic regulation, proteomics, metabolomics, and processing of large data sets promise to transform the current understanding of normal, abnormal, and malignant hematopoiesis. In this perspective, we discuss how the zebrafish model has proven beneficial for studying physiological and pathological hematopoiesis and how these novel technologies are transforming the field.

Transcriptional and epigenomic landscapes of CNS and non-CNS vascular endothelial cells

Vascular endothelial cell (EC) function depends on appropriate organ-specific molecular and cellular specializations. To explore genomic mechanisms that control this specialization, we have analyzed and compared the transcriptome, accessible chromatin, and DNA methylome landscapes from mouse brain, liver, lung, and kidney ECs. Analysis of transcription factor (TF) gene expression and TF motifs at candidate cis-regulatory elements reveals both shared and organ-specific EC regulatory networks. In the embryo, only those ECs that are adjacent to or within the central nervous system (CNS) exhibit canonical Wnt signaling, which correlates precisely with blood-brain barrier (BBB) differentiation and Zic3 expression. In the early postnatal brain, single-cell RNA-seq of purified ECs reveals (1) close relationships between veins and mitotic cells and between arteries and tip cells, (2) a division of capillary ECs into vein-like and artery-like classes, and (3) new endothelial subtype markers, including new validated tip cell markers.

ATACseqQC: a Bioconductor package for post-alignment quality assessment of ATAC-seq data

BACKGROUND

ATAC-seq (Assays for Transposase-Accessible Chromatin using sequencing) is a recently developed technique for genome-wide analysis of chromatin accessibility. Compared to earlier methods for assaying chromatin accessibility, ATAC-seq is faster and easier to perform, does not require cross-linking, has higher signal to noise ratio, and can be performed on small cell numbers. However, to ensure a successful ATAC-seq experiment, step-by-step quality assurance processes, including both wet lab quality control and in silico quality assessment, are essential. While several tools have been developed or adopted for assessing read quality, identifying nucleosome occupancy and accessible regions from ATAC-seq data, none of the tools provide a comprehensive set of functionalities for preprocessing and quality assessment of aligned ATAC-seq datasets.

RESULTS

We have developed a Bioconductor package, ATACseqQC, for easily generating various diagnostic plots to help researchers quickly assess the quality of their ATAC-seq data. In addition, this package contains functions to preprocess aligned ATAC-seq data for subsequent peak calling. Here we demonstrate the utilities of our package using 25 publicly available ATAC-seq datasets from four studies. We also provide guidelines on what the diagnostic plots should look like for an ideal ATAC-seq dataset.

CONCLUSIONS

This software package has been used successfully for preprocessing and assessing several in-house and public ATAC-seq datasets. Diagnostic plots generated by this package will facilitate the quality assessment of ATAC-seq data, and help researchers to evaluate their own ATAC-seq experiments as well as select high-quality ATAC-seq datasets from public repositories such as GEO to avoid generating hypotheses or drawing conclusions from low-quality ATAC-seq experiments. The software, source code, and documentation are freely available as a Bioconductor package at https://bioconductor.org/packages/release/bioc/html/ATACseqQC.html .

Epigenetics in teleost fish: From molecular mechanisms to physiological phenotypes

While the field of epigenetics is increasingly recognized to contribute to the emergence of phenotypes in mammalian research models across different developmental and generational timescales, the comparative biology of epigenetics in the large and physiologically diverse vertebrate infraclass of teleost fish remains comparatively understudied. The cypriniform zebrafish and the salmoniform rainbow trout and Atlantic salmon represent two especially important teleost orders, because they offer the unique possibility to comparatively investigate the role of epigenetic regulation in 3R and 4R duplicated genomes. In addition to their sequenced genomes, these teleost species are well-characterized model species for development and physiology, and therefore allow for an investigation of the role of epigenetic modifications in the emergence of physiological phenotypes during an organism's lifespan and in subsequent generations. This review aims firstly to describe the evolution of the repertoire of genes involved in key molecular epigenetic pathways including histone modifications, DNA methylation and microRNAs in zebrafish, rainbow trout, and Atlantic salmon, and secondly, to discuss recent advances in research highlighting a role for molecular epigenetics in shaping physiological phenotypes in these and other teleost models. Finally, by discussing themes and current limitations of the emerging field of teleost epigenetics from both theoretical and technical points of view, we will highlight future research needs and discuss how epigenetics will not only help address basic research questions in comparative teleost physiology, but also inform translational research including aquaculture, aquatic toxicology, and human disease.