Single-cell transcriptomic analysis of vascular endothelial cells in zebrafish embryos

Endothelial Cell Flow-Mediated Quiescence Is Temporally Regulated and Utilizes the Cell Cycle Inhibitor p27

BACKGROUND

Endothelial cells regulate their cell cycle as blood vessels remodel and transition to quiescence downstream of blood flow-induced mechanotransduction. Laminar blood flow leads to quiescence, but how flow-mediated quiescence is established and maintained is poorly understood.

METHODS

Primary human endothelial cells were exposed to laminar flow regimens and gene expression manipulations, and quiescence depth was analyzed via time-to-cell cycle reentry after flow cessation. Mouse and zebrafish endothelial expression patterns were examined via scRNA-seq (single-cell RNA sequencing) analysis, and mutant or morphant fish lacking p27 were analyzed for endothelial cell cycle regulation and in vivo cellular behaviors.

RESULTS

Arterial flow-exposed endothelial cells had a distinct transcriptome, and they first entered a deep quiescence, then transitioned to shallow quiescence under homeostatic maintenance conditions. In contrast, venous flow-exposed endothelial cells entered deep quiescence early that did not change with homeostasis. The cell cycle inhibitor p27 (CDKN1B) was required to establish endothelial flow-mediated quiescence, and expression levels positively correlated with quiescence depth. p27 loss in vivo led to endothelial cell cycle upregulation and ectopic sprouting, consistent with loss of quiescence. HES1 and ID3, transcriptional repressors of p27 upregulated by arterial flow, were required for quiescence depth changes and the reduced p27 levels associated with shallow quiescence.

CONCLUSIONS

Endothelial cell flow-mediated quiescence has unique properties and temporal regulation of quiescence depth that depends on the flow stimulus. These findings are consistent with a model whereby flow-mediated endothelial cell quiescence depth is temporally regulated downstream of p27 transcriptional regulation by HES1 and ID3. The findings are important in understanding endothelial cell quiescence misregulation that leads to vascular dysfunction and disease.

Blood flow regulates acvrl1 transcription via ligand-dependent Alk1 activity

Hereditary hemorrhagic telangiectasia (HHT) is an autosomal dominant disease characterized by the development of arteriovenous malformations (AVMs) that can result in significant morbidity and mortality. HHT is caused primarily by mutations in bone morphogenetic protein receptors ACVRL1/ALK1, a signaling receptor, or endoglin (ENG), an accessory receptor. Because overexpression of Acvrl1 prevents AVM development in both Acvrl1 and Eng null mice, enhancing ACVRL1 expression may be a promising approach to development of targeted therapies for HHT. Therefore, we sought to understand the molecular mechanism of ACVRL1 regulation. We previously demonstrated in zebrafish embryos that acvrl1 is predominantly expressed in arterial endothelial cells and that expression requires blood flow. Here, we document that flow dependence exhibits regional heterogeneity and that acvrl1 expression is rapidly restored after reinitiation of flow. Furthermore, we find that acvrl1 expression is significantly decreased in mutants that lack the circulating Alk1 ligand, Bmp10, and that, in the absence of flow, intravascular injection of BMP10 or the related ligand, BMP9, restores acvrl1 expression in an Alk1-dependent manner. Using a transgenic acvrl1:egfp reporter line, we find that flow and Bmp10 regulate acvrl1 at the level of transcription. Finally, we observe similar ALK1 ligand-dependent increases in ACVRL1 in human endothelial cells subjected to shear stress. These data suggest that ligand-dependent Alk1 activity acts downstream of blood flow to maintain or enhance acvrl1 expression via a positive feedback mechanism, and that ALK1 activating therapeutics may have dual functionality by increasing both ALK1 signaling flux and ACVRL1 expression.

Identification of New Markers of Angiogenic Sprouting Using Transcriptomics: New Role for RND3

BACKGROUND

New blood vessel formation requires endothelial cells to transition from a quiescent to an invasive phenotype. Transcriptional changes are vital for this switch, but a comprehensive genome-wide approach focused exclusively on endothelial cell sprout initiation has not been reported.

METHODS

Using a model of human endothelial cell sprout initiation, we developed a protocol to physically separate cells that initiate the process of new blood vessel formation (invading cells) from noninvading cells. We used this model to perform multiple transcriptomics analyses from independent donors to monitor endothelial gene expression changes.

RESULTS

Single-cell population analyses, single-cell cluster analyses, and bulk RNA sequencing revealed common transcriptomic changes associated with invading cells. We also found that collagenase digestion used to isolate single cells upregulated the Fos proto-oncogene transcription factor. Exclusion of Fos proto-oncogene expressing cells revealed a gene signature consistent with activation of signal transduction, morphogenesis, and immune responses. Many of the genes were previously shown to regulate angiogenesis and included multiple tip cell markers. Upregulation of SNAI1 (snail family transcriptional repressor 1), PTGS2 (prostaglandin synthase 2), and JUNB (JunB proto-oncogene) protein expression was confirmed in invading cells, and silencing JunB and SNAI1 significantly reduced invasion responses. Separate studies investigated rounding 3, also known as RhoE, which has not yet been implicated in angiogenesis. Silencing rounding 3 reduced endothelial invasion distance as well as filopodia length, fitting with a pathfinding role for rounding 3 via regulation of filopodial extensions. Analysis of in vivo retinal angiogenesis in Rnd3 heterozygous mice confirmed a decrease in filopodial length compared with wild-type littermates.

CONCLUSIONS

Validation of multiple genes, including rounding 3, revealed a functional role for this gene signature early in the angiogenic process. This study expands the list of genes associated with the acquisition of a tip cell phenotype during endothelial cell sprout initiation.

The development of brain pericytes requires expression of the transcription factor nkx3.1 in intermediate precursors

Brain pericytes are one of the critical cell types that regulate endothelial barrier function and activity, thus ensuring adequate blood flow to the brain. The genetic pathways guiding undifferentiated cells into mature pericytes are not well understood. We show here that pericyte precursor populations from both neural crest and head mesoderm of zebrafish express the transcription factor nkx3.1 develop into brain pericytes. We identify the gene signature of these precursors and show that an nkx3.1-, foxf2a-, and cxcl12b-expressing pericyte precursor population is present around the basilar artery prior to artery formation and pericyte recruitment. The precursors later spread throughout the brain and differentiate to express canonical pericyte markers. Cxcl12b-Cxcr4 signaling is required for pericyte attachment and differentiation. Further, both nkx3.1 and cxcl12b are necessary and sufficient in regulating pericyte number as loss inhibits and gain increases pericyte number. Through genetic experiments, we have defined a precursor population for brain pericytes and identified genes critical for their differentiation.

Direct male development in chromosomally ZZ zebrafish

The genetics of sex determination varies across taxa, sometimes even within a species. Major domesticated strains of zebrafish (Danio rerio), including AB and TU, lack a strong genetic sex determining locus, but strains more recently derived from nature, like Nadia (NA), possess a ZZ male/ZW female chromosomal sex-determination system. AB fish pass through a juvenile ovary stage, forming oocytes that survive in fish that become females but die in fish that become males. To understand mechanisms of gonad development in NA zebrafish, we studied histology and single cell transcriptomics in developing ZZ and ZW fish. ZW fish developed oocytes by 22 days post-fertilization (dpf) but ZZ fish directly formed testes, avoiding a juvenile ovary phase. Gonads of some ZW and WW fish, however, developed oocytes that died as the gonad became a testis, mimicking AB fish, suggesting that the gynogenetically derived AB strain is chromosomally WW. Single-cell RNA-seq of 19dpf gonads showed similar cell types in ZZ and ZW fish, including germ cells, precursors of gonadal support cells, steroidogenic cells, interstitial/stromal cells, and immune cells, consistent with a bipotential juvenile gonad. In contrast, scRNA-seq of 30dpf gonads revealed that cells in ZZ gonads had transcriptomes characteristic of testicular Sertoli, Leydig, and germ cells while ZW gonads had granulosa cells, theca cells, and developing oocytes. Hematopoietic and vascular cells were similar in both sex genotypes. These results show that juvenile NA zebrafish initially develop a bipotential gonad; that a factor on the NA W chromosome, or fewer than two Z chromosomes, is essential to initiate oocyte development; and without the W factor, or with two Z doses, NA gonads develop directly into testes without passing through the juvenile ovary stage. Sex determination in AB and TU strains mimics NA ZW and WW zebrafish, suggesting loss of the Z chromosome during domestication. Genetic analysis of the NA strain will facilitate our understanding of the evolution of sex determination mechanisms.

Endocardium gives rise to blood cells in zebrafish embryos

Previous studies have suggested that the endocardium contributes to hematopoiesis in murine embryos, although definitive evidence to demonstrate the hematopoietic potential of the endocardium is still missing. Here, we use a zebrafish embryonic model to test the emergence of hematopoietic progenitors from the endocardium. By using a combination of expression analysis, time-lapse imaging, and lineage-tracing approaches, we demonstrate that myeloid cells emerge from the endocardium in zebrafish embryos. Inhibition of Etv2/Etsrp or Scl/Tal1, two known master regulators of hematopoiesis and vasculogenesis, does not affect the emergence of endocardial-derived myeloid cells, while inhibition of Hedgehog signaling results in their reduction. Single-cell RNA sequencing analysis followed by experimental validation suggests that the endocardium is the major source of neutrophilic granulocytes. These findings will promote our understanding of alternative mechanisms involved in hematopoiesis, which are likely to be conserved between zebrafish and mammalian embryos.

Single-cell transcriptomics analysis of zebrafish brain reveals adverse effects of manganese on neurogenesis

Manganese (Mn) is considered as an important environmental risk factor for Parkinson's disease. Excessive exposure to Mn can damage various neural cells and affect the neurogenesis, resulting in neurological dysfunction. However, the specific mechanisms of Mn exposure affecting neurogenesis have not been well understood, including compositional changes and heterogeneity of various neural cells. Zebrafish have been successfully used as a neurotoxicity model due to its homology with mammals in several key regions of the brain, as well as its advantages such as small size. We performed single-cell RNA sequencing of zebrafish brains from normal and Mn-exposed groups. Our results suggested that low levels of Mn exposure activated neurogenesis in the zebrafish brain, including promoting the proliferation of neural progenitor cells and differentiation to newborn neurons and oligodendrocytes, while high levels of Mn exposure inhibited neurogenesis and neural function. Mn could affect neurogenesis through specific molecular pathways. In addition, Mn regulated intercellular communication and affected cellular communication in neural cells through specific signaling pathways. Taken together, our study elucidates the cellular composition of the zebrafish brain and adds to the understanding of the mechanisms involved in Mn-induced neurogenesis damage.

Direct Male Development in Chromosomally ZZ Zebrafish

The genetics of sex determination varies across taxa, sometimes even within a species. Major domesticated strains of zebrafish ( Danio rerio ), including AB and TU, lack a strong genetic sex determining locus, but strains more recently derived from nature, like Nadia (NA), possess a ZZ male/ZW female chromosomal sex-determination system. AB strain fish pass through a juvenile ovary stage, forming oocytes that survive in fish that become females but die in fish that become males. To understand mechanisms of gonad development in NA zebrafish, we studied histology and single cell transcriptomics in developing ZZ and ZW fish. ZW fish developed oocytes by 22 days post-fertilization (dpf) but ZZ fish directly formed testes, avoiding a juvenile ovary phase. Gonads of some ZW and WW fish, however, developed oocytes that died as the gonad became a testis, mimicking AB fish, suggesting that the gynogenetically derived AB strain is chromosomally WW. Single-cell RNA-seq of 19dpf gonads showed similar cell types in ZZ and ZW fish, including germ cells, precursors of gonadal support cells, steroidogenic cells, interstitial/stromal cells, and immune cells, consistent with a bipotential juvenile gonad. In contrast, scRNA-seq of 30dpf gonads revealed that cells in ZZ gonads had transcriptomes characteristic of testicular Sertoli, Leydig, and germ cells while ZW gonads had granulosa cells, theca cells, and developing oocytes. Hematopoietic and vascular cells were similar in both sex genotypes. These results show that juvenile NA zebrafish initially develop a bipotential gonad; that a factor on the NA W chromosome or fewer than two Z chromosomes is essential to initiate oocyte development; and without the W factor or with two Z doses, NA gonads develop directly into testes without passing through the juvenile ovary stage. Sex determination in AB and TU strains mimics NA ZW and WW zebrafish, suggesting loss of the Z chromosome during domestication. Genetic analysis of the NA strain will facilitate our understanding of the evolution of sex determination mechanisms.

Sinus venosus adaptation models prolonged cardiovascular disease and reveals insights into evolutionary transitions of the vertebrate heart

How two-chambered hearts in basal vertebrates have evolved from single-chamber hearts found in ancestral chordates remains unclear. Here, we show that the teleost sinus venosus (SV) is a chamber-like vessel comprised of an outer layer of smooth muscle cells. We find that in adult zebrafish nr2f1a mutants, which lack atria, the SV comes to physically resemble the thicker bulbus arteriosus (BA) at the arterial pole of the heart through an adaptive, hypertensive response involving smooth muscle proliferation due to aberrant hemodynamic flow. Single cell transcriptomics show that smooth muscle and endothelial cell populations within the adapting SV also take on arterial signatures. Bulk transcriptomics of the blood sinuses flanking the tunicate heart reinforce a model of greater equivalency in ancestral chordate BA and SV precursors. Our data simultaneously reveal that secondary complications from congenital heart defects can develop in adult zebrafish similar to those in humans and that the foundation of equivalency between flanking auxiliary vessels may remain latent within basal vertebrate hearts.

Disentangling brain vasculature in neurogenesis and neurodegeneration using single-cell transcriptomics

The vasculature is increasingly recognized to impact brain function in health and disease across the life span. During embryonic brain development, angiogenesis and neurogenesis are tightly coupled, coordinating the proliferation, differentiation, and migration of neural and glial progenitors. In the adult brain, neurovascular interactions continue to play essential roles in maintaining brain function and homeostasis. This review focuses on recent advances that leverage single-cell transcriptomics of vascular cells to uncover their subtypes, their organization and zonation in the embryonic and adult brain, and how dysfunction in neurovascular and gliovascular interactions contributes to the pathogenesis of neurodegenerative diseases. Finally, we highlight key challenges for future research in neurovascular biology.