Genomic Characterization of Endothelial Enhancers Reveals a Multifunctional Role for NR2F2 in Regulation of Arteriovenous Gene Expression

Intramyocardial Sprouting Tip Cells Specify Coronary Arterialization

BACKGROUND

The elaborate patterning of coronary arteries critically supports the high metabolic activity of the beating heart. How coronary endothelial cells coordinate hierarchical vascular remodeling and achieve arteriovenous specification remains largely unknown. Understanding the molecular and cellular cues that pattern coronary arteries is crucial to develop innovative therapeutic strategies that restore functional perfusion within the ischemic heart.

METHODS

Single-cell transcriptomics and histological validation were used to delineate heterogeneous transcriptional states of the developing and mature coronary endothelium with a focus on sprouting endothelium and arterial cell specification. Genetic lineage tracing and high-resolution 3-dimensional imaging were used to characterize the origin and mechanisms of coronary angiogenic sprouting, as well as to fate-map selective endothelial lineages. Integration of single-cell transcriptomic data from ischemic adult mouse hearts and human embryonic data served to assess the conservation of transcriptional states across development, disease, and species.

RESULTS

We discover that coronary arteries originate from cells that have previously transitioned through a specific tip cell phenotype. We identify nonoverlapping intramyocardial and subepicardial tip cell populations with differential gene expression profiles and regulatory pathways. Esm1-lineage tracing confirmed that intramyocardial tip cells selectively contribute to coronary arteries and endocardial tunnels, but not veins. Notably, prearterial cells are detected from development stages to adulthood, increasingly in response to ischemic injury, and in human embryos, suggesting that tip cell-to-artery specification is a conserved mechanism.

CONCLUSIONS

A tip cell-to-artery specification mechanism drives arterialization of the intramyocardial plexus and endocardial tunnels throughout life and is reactivated upon ischemic injury. Differential sprouting programs govern the formation and specification of the venous and arterial coronary plexus.

Machine Learning Uncovers Vascular Endothelial Cell Identity Genes by Expression Regulation Features in Single Cells

Deciphering cell identity genes is pivotal to understanding cell differentiation, development, and many diseases involving cell identity dysregulation. Here, we introduce SCIG, a machine-learning method to uncover cell identity genes in single cells. In alignment with recent reports that cell identity genes are regulated with unique epigenetic signatures, we found cell identity genes exhibit distinctive genetic sequence signatures, e.g., unique enrichment patterns of cis-regulatory elements. Using these genetic sequence signatures, along with gene expression information from single-cell RNA-seq data, enables SCIG to uncover the identity genes of a cell without a need for comparison to other cells. Cell identity gene score defined by SCIG surpassed expression value in network analysis to uncover master transcription factors regulating cell identity. Applying SCIG to the human endothelial cell atlas revealed that the tissue microenvironment is a critical supplement to master transcription factors for cell identity refinement. SCIG is publicly available at https://github.com/kaifuchenlab/SCIG , offering a valuable tool for advancing cell differentiation, development, and regenerative medicine research.

An Endothelial Cell Is Not Simply an Endothelial Cell

Endothelial cells (ECs) are a multifaceted component of the vascular system with roles in immunity, maintaining tissue fluid balance, and vascular tone. Dysregulation or dysfunction of ECs can have far-reaching implications, leading pathologies ranging from cardiovascular diseases, such as hypertension and atherosclerosis, ischemia, chronic kidney disease, blood-brain barrier integrity, dementia, and tumor metastasis. Recent advancements in regenerative medicine have highlighted the potential of stem cell-derived ECs, particularly from induced pluripotent stem cells, to treat ischemic tissues, as well as models of vascular integrity. This review summarizes what is known in the generation of ECs with an emphasis on tissue-specific ECs and EC subphenotypes important in the development of targeted cell-based therapies for patient treatment.

Genome-Wide Methylation Analysis Reveals a KCNK3-Prominent Causal Cascade on Hypertension

BACKGROUND

Despite advances in understanding hypertension's genetic structure, how noncoding genetic variants influence it remains unclear. Studying their interaction with DNA methylation is crucial to deciphering this complex disease's genetic mechanisms.

METHODS

We investigated the genetic and epigenetic interplay in hypertension using whole-genome bisulfite sequencing. Methylation profiling in 918 males revealed allele-specific methylation and methylation quantitative trait loci. We engineered rs1275988T/C mutant mice using CRISPR (clustered regularly interspaced short palindromic repeats)/Cas9 (CRISPR-associated protein 9), bred them for homozygosity, and subjected them to a high-salt diet. Telemetry captured their cardiovascular metrics. Protein-DNA interactions were elucidated using DNA pull-downs, mass spectrometry, and Western blots. A wire myograph assessed vascular function, and analysis of the Kcnk3 gene methylation highlighted the mutation's role in hypertension.

RESULTS

We discovered that DNA methylation-associated genetic effects, especially in non-cytosine-phosphate-guanine (non-CpG) island and noncoding distal regulatory regions, significantly contribute to hypertension predisposition. We identified distinct methylation quantitative trait locus patterns in the hypertensive population and observed that the onset of hypertension is influenced by the transmission of genetic effects through the demethylation process. By evidence-driven prioritization and in vivo experiments, we unearthed rs1275988 in a cell type-specific enhancer as a notable hypertension causal variant, intensifying hypertension through the modulation of local DNA methylation and consequential alterations in Kcnk3 gene expression and vascular remodeling. When exposed to a high-salt diet, mice with the rs1275988C/C genotype exhibited exacerbated hypertension and significant vascular remodeling, underscored by increased aortic wall thickness. The C allele of rs1275988 was associated with elevated DNA methylation levels, driving down the expression of the Kcnk3 gene by attenuating Nr2f2 (nuclear receptor subfamily 2 group F member 2) binding at the enhancer locus.

CONCLUSIONS

Our research reveals new insights into the complex interplay between genetic variations and DNA methylation in hypertension. We underscore hypomethylation's potential in hypertension onset and identify rs1275988 as a causal variant in vascular remodeling. This work advances our understanding of hypertension's molecular mechanisms and encourages personalized health care strategies.

Cell Fate Regulation During the Development of Infantile Hemangioma

As the most common benign vascular tumor in infants, infantile hemangioma (IH) is characterized by rapid growth and vasculogenesis early in infancy, followed by spontaneous involution into fibrofatty tissues over time. Extensive evidence suggests that IH originates from hemangioma stem cells (HemSCs), a group of stem cells with clonal expansion and multi-directional differentiation capacity. However, the intricate mechanisms governing the cell fate transition of HemSCs during IH development remain elusive. Here we comprehensively examine the cellular composition of IH, emphasizing the nuanced properties of various IH cell types and their correlation with the clinical features of the tumor. We also summarize the current understanding of the regulatory pathways directing HemSC differentiation into endothelial cells (ECs), pericytes, and adipocytes throughout the stages of IH progression and involution. Furthermore, we discuss recent advances in unraveling the transcriptional and epigenetic regulation of EC and adipocyte development under physiological conditions, which offer crucial perspectives for understanding IH pathogenesis.

A Retinoic Acid:YAP1 signaling axis controls atrial lineage commitment

Vitamin A/Retinoic Acid (Vit A/RA) signaling is essential for heart development. In cardiac progenitor cells (CPCs), RA signaling induces the expression of atrial lineage genes while repressing ventricular genes, thereby promoting the acquisition of an atrial cardiomyocyte cell fate. To achieve this, RA coordinates a complex regulatory network of downstream effectors that is not fully identified. To address this gap, we applied a functional genomics approach (i.e scRNAseq and snATACseq) to untreated and RA-treated human embryonic stem cells (hESCs)-derived CPCs. Unbiased analysis revealed that the Hippo effectors YAP1 and TEAD4 are integrated with the atrial transcription factor enhancer network, and that YAP1 is necessary for activation of RA-enhancers in CPCs. Furthermore, in vivo analysis of control and conditionally YAP1 KO mouse embryos (Sox2-cre) revealed that the expression of atrial lineage genes, such as NR2F2, is compromised by YAP1 deletion in the CPCs of the second heart field. Accordingly, we found that YAP1 is required for the formation of an atrial chamber but is dispensable for the formation of a ventricle, in hESC-derived patterned cardiac organoids. Overall, our findings revealed that YAP1 is a non-canonical effector of RA signaling essential for the acquisition of atrial lineages during cardiogenesis.

A KLF2-BMPER-Smad1/5 checkpoint regulates high fluid shear stress-mediated artery remodeling

Vascular remodeling to match arterial diameter to tissue requirements commonly fails in ischemic disease. Endothelial cells sense fluid shear stress (FSS) from blood flow to maintain FSS within a narrow range in healthy vessels. Thus, high FSS induces vessel outward remodeling, but mechanisms are poorly understood. We previously reported that Smad1/5 is maximally activated at physiological FSS. Smad1/5 limits Akt activation, suggesting that inhibiting Smad1/5 may facilitate outward remodeling. Here we report that high FSS suppresses Smad1/5 by elevating KLF2, which induces the bone morphogenetic protein (BMP) pathway inhibitor, BMP-binding endothelial regulator (BMPER), thereby de-inhibiting Akt. In mice, surgically induced high FSS elevated BMPER expression, inactivated Smad1/5 and induced vessel outward remodeling. Endothelial BMPER deletion impaired blood flow recovery and vascular remodeling. Blocking endothelial cell Smad1/5 activation with BMP9/10 blocking antibodies improved vascular remodeling in mouse models of type 1 and type 2 diabetes. Suppression of Smad1/5 is thus a potential therapeutic approach for ischemic disease.

Differential endothelial cell cycle status in postnatal retinal vessels revealed using a novel PIP-FUCCI reporter and zonation analysis

Cell cycle regulation is critical to blood vessel formation and function, but how the endothelial cell cycle integrates with vascular regulation is not well-understood, and available dynamic cell cycle reporters do not precisely distinguish all cell cycle stage transitions in vivo. Here we characterized a recently developed improved cell cycle reporter (PIP-FUCCI) that precisely delineates S phase and the S/G2 transition. Live image analysis of primary endothelial cells revealed predicted temporal changes and well-defined stage transitions. A new inducible mouse cell cycle reporter allele was selectively expressed in postnatal retinal endothelial cells upon Cre-mediated activation and predicted endothelial cell cycle status. We developed a semi-automated zonation program to define endothelial cell cycle status in spatially defined and developmentally distinct retinal areas and found predicted cell cycle stage differences in arteries, veins, and remodeled and angiogenic capillaries. Surprisingly, the predicted dearth of S-phase proliferative tip cells relative to stalk cells at the vascular front was accompanied by an unexpected enrichment for endothelial tip and stalk cells in G2, suggesting G2 stalling as a contribution to tip-cell arrest and dynamics at the front. Thus, this improved reporter precisely defines endothelial cell cycle status in vivo and reveals novel G2 regulation that may contribute to unique aspects of blood vessel network expansion.

Atlas of cardiac endothelial cell enhancer elements linking the mineralocorticoid receptor to pathological gene expression

Endothelial cells play crucial roles in physiology and are increasingly recognized as therapeutic targets in cardiovascular disease. Here, we analyzed the regulatory landscape of cardiac endothelial cells by assessing chromatin accessibility, histone modifications, and 3D chromatin organization and confirmed the functional relevance of enhancer-promoter interactions by CRISPRi-mediated enhancer silencing. We used this dataset to explore mechanisms of transcriptional regulation in cardiovascular disease and compared six different experimental models of heart failure, hypertension, or diabetes. Enhancers that regulate gene expression in diseased endothelial cells were enriched with binding sites for a distinct set of transcription factors, including the mineralocorticoid receptor (MR), a known drug target in heart failure and hypertension. For proof of concept, we applied endothelial cell-specific MR deletion in mice to confirm MR-dependent gene expression and predicted direct MR target genes. Overall, we have compiled here a comprehensive atlas of cardiac endothelial cell enhancer elements that provides insight into the role of transcription factors in cardiovascular disease.

Differential endothelial cell cycle status in postnatal retinal vessels revealed using a novel PIP-FUCCI reporter and zonation analysis

Cell cycle regulation is critical to blood vessel formation and function, but how the endothelial cell cycle integrates with vascular regulation is not well-understood, and available dynamic cell cycle reporters do not precisely distinguish all cell cycle stage transitions in vivo. Here we characterized a recently developed improved cell cycle reporter (PIP-FUCCI) that precisely delineates S phase and the S/G2 transition. Live image analysis of primary endothelial cells revealed predicted temporal changes and well-defined stage transitions. A new inducible mouse cell cycle reporter allele was selectively expressed in postnatal retinal endothelial cells upon Cre-mediated activation and predicted endothelial cell cycle status. We developed a semi-automated zonation program to define endothelial cell cycle status in spatially defined and developmentally distinct retinal areas and found predicted cell cycle stage differences in arteries, veins, and remodeled and angiogenic capillaries. Surprisingly, the predicted dearth of proliferative tip cells at the vascular front was accompanied by an unexpected enrichment for endothelial tip cells in G2, suggesting G2 stalling as a contribution to tip-cell arrest. Thus, this improved reporter precisely defines endothelial cell cycle status in vivo and reveals novel G2 regulation that may contribute to unique aspects of blood vessel network expansion.

Transcription factors regulating vasculogenesis and angiogenesis

Transcription factors (TFs) play a crucial role in regulating the dynamic and precise patterns of gene expression required for the initial specification of endothelial cells (ECs), and during endothelial growth and differentiation. While sharing many core features, ECs can be highly heterogeneous. Differential gene expression between ECs is essential to pattern the hierarchical vascular network into arteries, veins and capillaries, to drive angiogenic growth of new vessels, and to direct specialization in response to local signals. Unlike many other cell types, ECs have no single master regulator, instead relying on differing combinations of a necessarily limited repertoire of TFs to achieve tight spatial and temporal activation and repression of gene expression. Here, we will discuss the cohort of TFs known to be involved in directing gene expression during different stages of mammalian vasculogenesis and angiogenesis, with a primary focus on development.

Transcriptional regulators of arterial and venous identity in the developing mammalian embryo

The complex and hierarchical vascular network of arteries, veins, and capillaries features considerable endothelial heterogeneity, yet the regulatory pathways directing arteriovenous specification, differentiation, and identity are still not fully understood. Recent advances in analysis of endothelial-specific gene-regulatory elements, single-cell RNA sequencing, and cell lineage tracing have both emphasized the importance of transcriptional regulation in this process and shed considerable light on the mechanism and regulation of specification within the endothelium. In this review, we discuss recent advances in our understanding of how endothelial cells acquire arterial and venous identity and the role different transcription factors play in this process.

Gene regulation for inflammation and inflammation resolution differs between umbilical arterial and venous endothelial cells

Systemic inflammation affects the whole vasculature, yet whether arterial and venous endothelial cells differ in their abilities to mediate inflammation and to return to homeostasis after an inflammatory stimulus has not been addressed thoroughly. We assessed gene-expression profiles in isolated endothelial cells from human umbilical arteries (HUAEC) or veins (HUVEC) under basal conditions, after TNF-α stimulation and various time points after TNF-α removal to allow reinstatement of homeostasis. TNF-α regulates the expression of different sets of transcripts that are significantly changed only in HUAEC, only in HUVEC or changed in both. We identified three types of gene regulation, i.e. genes that were significantly regulated after 24 h of TNF-α stimulation but no longer when TNF-α was removed (homeostatic regulation), genes that maintained significantly regulated after TNF-α removal (not homeostatic regulation) and genes that were only significantly regulated when TNF-α was removed (post-regulation). HUAEC and HUVEC quantitatively differed in these types of gene regulation, with relatively more genes being post-regulated in HUAEC. In conclusion our data demonstrate that HUAEC and HUVEC respond intrinsically different to an inflammatory insult. Whether this holds true for all endothelial cells and its relevance for inflammatory insults in different organs during systemic inflammation warrants further studies.

Inflammation and DKK1-induced AKT activation contribute to endothelial dysfunction following NR2F2 loss

NR2F2 is expressed in endothelial cells (ECs) and Nr2f2 knockout produces lethal cardiovascular defects. In humans, reduced NR2F2 expression is associated with cardiovascular diseases including congenital heart disease and atherosclerosis. Here, NR2F2 silencing in human primary ECs led to inflammation, endothelial-to-mesenchymal transition (EndMT), proliferation, hypermigration, apoptosis-resistance, and increased production of reactive oxygen species. These changes were associated with STAT and AKT activation along with increased production of DKK1. Co-silencing DKK1 and NR2F2 prevented NR2F2-loss-induced STAT and AKT activation and reversed EndMT. Serum DKK1 concentrations were elevated in patients with pulmonary arterial hypertension (PAH) and DKK1 was secreted by ECs in response to in vitro loss of either BMPR2 or CAV1, which are genetic defects associated with the development of PAH. In human primary ECs, NR2F2 suppressed DKK1, whereas its loss conversely induced DKK1 and disrupted endothelial homeostasis, promoting phenotypic abnormalities associated with pathologic vascular remodeling. Activating NR2F2 or blocking DKK1 may be useful therapeutic targets for treating chronic vascular diseases associated with EC dysfunction.NEW & NOTEWORTHY NR2F2 loss in the endothelial lining of blood vessels is associated with cardiovascular disease. Here, NR2F2-silenced human endothelial cells were inflammatory, proliferative, hypermigratory, and apoptosis-resistant with increased oxidant stress and endothelial-to-mesenchymal transition. DKK1 was induced in NR2F2-silenced endothelial cells, while co-silencing NR2F2 and DKK1 prevented NR2F2-loss-associated abnormalities in endothelial signaling and phenotype. Activating NR2F2 or blocking DKK1 may be useful therapeutic targets for treating vascular diseases associated with endothelial dysfunction.

The tissue-specific transcriptional landscape underlines the involvement of endothelial cells in health and disease

Endothelial cells (ECs) that line vascular and lymphatic vessels are being increasingly recognized as important to organ function in health and disease. ECs participate not only in the trafficking of gases, metabolites, and cells between the bloodstream and tissues but also in the angiocrine-based induction of heterogeneous parenchymal cells, which are unique to their specific tissue functions. The molecular mechanisms regulating EC heterogeneity between and within different tissues are modeled during embryogenesis and become fully established in adults. Any changes in adult tissue homeostasis induced by aging, stress conditions, and various noxae may reshape EC heterogeneity and induce specific transcriptional features that condition a functional phenotype. Heterogeneity is sustained via specific genetic programs organized through the combinatory effects of a discrete number of transcription factors (TFs) that, at the single tissue-level, constitute dynamic networks that are post-transcriptionally and epigenetically regulated. This review is focused on outlining the TF-based networks involved in EC specialization and physiological and pathological stressors thought to modify their architecture.

Single-cell transcriptomic heterogeneity between conduit and resistance mesenteric arteries in rats

The endothelium contains morphologically similar cells throughout the vasculature, but individual cells along the length of a single vascular tree or in different regional circulations function dissimilarly. When observations made in large arteries are extrapolated to explain the function of endothelial cells (ECs) in the resistance vasculature, only a fraction of these observations are consistent between artery sizes. To what extent endothelial (EC) and vascular smooth muscle cells (VSMCs) from different arteriolar segments of the same tissue differ phenotypically at the single-cell level remains unknown. Therefore, single-cell RNA-seq (10x Genomics) was performed using a 10X Genomics Chromium system. Cells were enzymatically digested from large (>300 µm) and small (<150 µm) mesenteric arteries from nine adult male Sprague-Dawley rats, pooled to create six samples (3 rats/sample, 3 samples/group). After normalized integration, the dataset was scaled before unsupervised cell clustering and cluster visualization using UMAP plots. Differential gene expression analysis allowed us to infer the biological identity of different clusters. Our analysis revealed 630 and 641 differentially expressed genes (DEGs) between conduit and resistance arteries for ECs and VSMCs, respectively. Gene ontology analysis (GO-Biological Processes, GOBP) of scRNA-seq data discovered 562 and 270 pathways for ECs and VSMCs, respectively, that differed between large and small arteries. We identified eight and seven unique ECs and VSMCs subpopulations, respectively, with DEGs and pathways identified for each cluster. These results and this dataset allow the discovery and support of novel hypotheses needed to identify mechanisms that determine the phenotypic heterogeneity between conduit and resistance arteries.

Cell lineage-specific mitochondrial resilience during mammalian organogenesis

Mitochondrial activity differs markedly between organs, but it is not known how and when this arises. Here we show that cell lineage-specific expression profiles involving essential mitochondrial genes emerge at an early stage in mouse development, including tissue-specific isoforms present before organ formation. However, the nuclear transcriptional signatures were not independent of organelle function. Genetically disrupting intra-mitochondrial protein synthesis with two different mtDNA mutations induced cell lineage-specific compensatory responses, including molecular pathways not previously implicated in organellar maintenance. We saw downregulation of genes whose expression is known to exacerbate the effects of exogenous mitochondrial toxins, indicating a transcriptional adaptation to mitochondrial dysfunction during embryonic development. The compensatory pathways were both tissue and mutation specific and under the control of transcription factors which promote organelle resilience. These are likely to contribute to the tissue specificity which characterizes human mitochondrial diseases and are potential targets for organ-directed treatments.

Transcriptomic Heterogeneity of Skin Across Different Anatomic Sites

Multiomic studies, including RNA sequencing, single-cell RNA sequencing, and epigenomics, can provide insight into the connection between anatomically heterogeneous gene expression profile of the skin and dermatoses-predisposed sites, in which RNA sequencing is essential. Therefore, in this study, 159 skin samples collected mainly from discarded normal skin tissue during surgical treatment for benign skin tumors were used for RNA sequencing. On the basis of cluster analysis, the skin was divided into four regions, with each region showing specific physiological characteristics through differentially expressed gene analysis. The results showed that the head and neck region, perineum, and palmoplantar area were closely associated with lipid metabolism, hormone metabolism, blood circulation, and related neural regulation, respectively. Transcription factor enrichment indicated that different regions were associated with the development of adjacent tissues. Specifically, the head and neck region, trunk and extremities, perineum, and palmoplantar area were associated with the central nervous, axial, urogenital, and vascular systems, respectively. The results were imported into an open website (https://dermvis.github.io/) for retrieval. Our transcriptomic data elucidated that human skin exhibits transcriptomic heterogeneity reflecting physiological and developmental variation at different anatomic sites and provided guidance for further studies on skin development and dermatoses predisposed sites.

Vascular endothelial cell development and diversity

Vascular endothelial cells form the inner layer of blood vessels where they have a key role in the development and maintenance of the functional circulatory system and provide paracrine support to surrounding non-vascular cells. Technical advances in the past 5 years in single-cell genomics and in in vivo genetic labelling have facilitated greater insights into endothelial cell development, plasticity and heterogeneity. These advances have also contributed to a new understanding of the timing of endothelial cell subtype differentiation and its relationship to the cell cycle. Identification of novel tissue-specific gene expression patterns in endothelial cells has led to the discovery of crucial signalling pathways and new interactions with other cell types that have key roles in both tissue maintenance and disease pathology. In this Review, we describe the latest findings in vascular endothelial cell development and diversity, which are often supported by large-scale, single-cell studies, and discuss the implications of these findings for vascular medicine. In addition, we highlight how techniques such as single-cell multimodal omics, which have become increasingly sophisticated over the past 2 years, are being utilized to study normal vascular physiology as well as functional perturbations in disease.

ETS transcription factors: Multifaceted players from cancer progression to tumor immunity

The E26 transformation specific (ETS) family comprises 28 transcription factors, the majority of which are involved in tumor initiation and development. Serving as a group of functionally heterogeneous gene regulators, ETS factors possess a structurally conserved DNA-binding domain. As one of the most prominent families of transcription factors that control diverse cellular functions, ETS activation is modulated by multiple intracellular signaling pathways and post-translational modifications. Disturbances in ETS activity often lead to abnormal changes in oncogenicity, including cancer cell survival, growth, proliferation, metastasis, genetic instability, cell metabolism, and tumor immunity. This review systematically addresses the basics and advances in studying ETS factors, from their tumor relevance to clinical translational utility, with a particular focus on elucidating the role of ETS family in tumor immunity, aiming to decipher the vital role and clinical potential of regulation of ETS factors in the cancer field.

Impaired Human Cardiac Cell Development due to NOTCH1 Deficiency

BACKGROUND

NOTCH1 pathogenic variants are implicated in multiple types of congenital heart defects including hypoplastic left heart syndrome, where the left ventricle is underdeveloped. It is unknown how NOTCH1 regulates human cardiac cell lineage determination and cardiomyocyte proliferation. In addition, mechanisms by which NOTCH1 pathogenic variants lead to ventricular hypoplasia in hypoplastic left heart syndrome remain elusive.

METHODS

CRISPR (Clustered Regularly Interspaced Short Palindromic Repeats)/Cas9 genome editing was utilized to delete NOTCH1 in human induced pluripotent stem cells. Cardiac differentiation was carried out by sequential modulation of WNT signaling, and NOTCH1 knockout and wild-type differentiating cells were collected at day 0, 2, 5, 10, 14, and 30 for single-cell RNA-seq.

RESULTS

Human NOTCH1 knockout induced pluripotent stem cells are able to generate functional cardiomyocytes and endothelial cells, suggesting that NOTCH1 is not required for mesoderm differentiation and cardiovascular development in vitro. However, disruption of NOTCH1 blocks human ventricular-like cardiomyocyte differentiation but promotes atrial-like cardiomyocyte generation through shortening the action potential duration. NOTCH1 deficiency leads to defective proliferation of early human cardiomyocytes, and transcriptomic analysis indicates that pathways involved in cell cycle progression and mitosis are downregulated in NOTCH1 knockout cardiomyocytes. Single-cell transcriptomic analysis reveals abnormal cell lineage determination of cardiac mesoderm, which is manifested by the biased differentiation toward epicardial and second heart field progenitors at the expense of first heart field progenitors in NOTCH1 knockout cell populations.

CONCLUSIONS

NOTCH1 is essential for human ventricular-like cardiomyocyte differentiation and proliferation through balancing cell fate determination of cardiac mesoderm and modulating cell cycle progression. Because first heart field progenitors primarily contribute to the left ventricle, we speculate that pathogenic NOTCH1 variants lead to biased differentiation of first heart field progenitors, blocked ventricular-like cardiomyocyte differentiation, and defective cardiomyocyte proliferation, which collaboratively contribute to left ventricular hypoplasia in hypoplastic left heart syndrome.

Single-cell RNA sequencing reveals resident progenitor and vascularization-associated cell subpopulations in rat annulus fibrosus

Background

One of the main causes of low back pain is intervertebral disc degeneration (IDD). Annulus fibrosus (AF) is important for the integrity and functions of the intervertebral disc (IVD). However, the resident functional cell components such as progenitors and vascularization-associated cells in AF are yet to be fully identified.

Purpose

Identification of functional AF cell subpopulations including resident progenitors and vascularization-associated cells.

Methods

In this study, the single-cell RNA sequencing data of rat IVDs from a public database were analyzed using Seurat for cell clustering, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) for functional analysis, StemID for stem cell identification, Monocle and RNA velocity for pseudotime differentiation trajectory validation, single-cell regulatory network inference and clustering (SCENIC) for gene regulatory network (GRN) analysis, and CellChat for cell-cell interaction analysis. Immunostaining on normal and degenerated rat IVDs, as well as human AF, was used for validations.

Results

From the data analysis, seven AF cell clusters were identified, including two newly discovered functional clusters, the Grem1 ⁺ subpopulation and the Lum ⁺ ​subpopulation. The Grem1 ⁺ subpopulation had progenitor characteristics, while the Lum ⁺ ​subpopulation was associated with vascularization during IDD. The GRN analysis showed that Sox9 and Id1 were among the key regulators in the Grem1 ⁺ subpopulation, and Nr2f2 and Creb5 could be responsible for the vascularization function in the Lum ⁺ ​subpopulation. Cell-cell interaction analysis revealed highly regulated cellular communications between these cells, and multiple signaling networks including PDGF and MIF signaling pathways were involved in the interactions.

Conclusions

Our results revealed two new functional AF cell subpopulations, with stemness and vascularization induction potential, respectively.

The Translational potential of this article

These findings complement our knowledge about IVDs, especially the AF, and in return provide potential cell source and regulation targets for IDD treatment and tissue repair. The existence of the cell subpopulations was also validated in human AF, which strengthen the clinical relevance of the findings.

NR2F2 Regulates Cell Proliferation and Immunomodulation in Whartons’ Jelly Stem Cells

(1) Background: Wharton’s Jelly stem cells (WJ-MSCs) are multipotent mesenchymal stem cells that can proliferate rapidly and have low immunogenicity. Therefore, WJ-MSCs have gained considerable attention in the fields of immunomodulation and disease treatment and have entered clinical trials for the treatment of various diseases. Therefore, it is crucial to study the underlying mechanisms of WJ-MSCs proliferation, immune regulation, and disease treatment. Nuclear Receptor Subfamily 2 Group F Member 2 (NR2F2) is a transcription factor that is involved in the regulation of many different genes. However, it remains unknown how NR2F2 regulates stem cell identity in WJ-MSCs. (2) Methods: We used RNAi technology to knock down NR2F2 in WJ-MSCs, and studied the regulatory role of NR2F2 in WJ-MSCs by MTT, flow cytometry, RNA-seq, and other methods. We also utilized a co-culture system in which NR2F2-depleted WJ-MSCs with MH7A and HCT116/HepG2 were used to investigate the role of NR2F2 in immunomodulation and the inhibition of cancer cell growth. (3) Results: NR2F2 knockdown resulted in decreased expressions of Cyclin D1 and CDK4, slower cell proliferation, and increased expressions of IL6 and IL8. Furthermore, Cyclin D1, CDK4, and inflammatory factors were increased in human rheumatoid fibroblast-like synoviocyte line MH7A if co-cultured with NR2F2 depleted WJ-MSCs. In addition, we observed increased p53, decreased BCL-2, and increased cell apoptosis in liver cancer cell line HepG2 if co-cultured with NR2F2-depleted WJ-MSCs. (4) Conclusions: NR2F2 not only plays an important role in the cell cycle and immune regulation of WJ-MSCs but also has potential effects on the WJ-MSCs treatment of related diseases.

TREM2 macrophages induced by human lipids drive inflammation in acne lesions

Acne affects 1 in 10 people globally, often resulting in disfigurement. The disease involves excess production of lipids, particularly squalene, increased growth of Cutibacterium acnes, and a host inflammatory response with foamy macrophages. By combining single-cell and spatial RNA sequencing as well as ultrahigh-resolution Seq-Scope analyses of early acne lesions on back skin, we identified TREM2 macrophages expressing lipid metabolism and proinflammatory gene programs in proximity to hair follicle epithelium expressing squalene epoxidase. We established that the addition of squalene induced differentiation of TREM2 macrophages in vitro, which were unable to kill C. acnes. The addition of squalene to macrophages inhibited induction of oxidative enzymes and scavenged oxygen free radicals, providing an explanation for the efficacy of topical benzoyl peroxide in the clinical treatment of acne. The present work has elucidated the mechanisms by which TREM2 macrophages and unsaturated lipids, similar to their involvement in atherosclerosis, may contribute to the pathogenesis of acne.

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.

"Enhancing" mechanosensing: Enhancers and enhancer-derived long non-coding RNAs in endothelial response to flow

Endothelial cells (ECs), uniquely localized and strategically forming the inner lining of vascular wall, constitute the largest cell surface by area in the human body. The dynamic sensing and response of ECs to mechanical cues, especially shear stress, is crucial for maintenance of vascular homeostasis. It is well recognized that different flow patterns associated with atheroprotective vs atheroprone regions in the arterial tree, result in distinct EC functional phenotypes with differential transcriptome profiles. Mounting evidence has demonstrated an integrative and essential regulatory role of non-coding genome in EC biology. In particular, recent studies have begun to reveal the importance of enhancers and enhancer-derived transcripts in flow-regulated EC gene expression and function. In this minireview, we summarize studies in this area and discuss examples in support of the emerging importance of enhancers and enhancer(-derived) long non-coding RNAs (elncRNAs) in EC mechanosensing, with a focus on flow-responsive EC transcription. Finally, we will provide perspective and discuss standing questions to elucidate the role of these novel regulators in EC mechanobiology.

Hyperglycaemia up-regulates placental growth factor (PlGF) expression and secretion in endothelial cells via suppression of PI3 kinase-Akt signalling and activation of FOXO1

Placenta growth factor (PlGF) is a pro-inflammatory angiogenic mediator that promotes many pathologies including diabetic complications and atherosclerosis. Widespread endothelial dysfunction precedes the onset of these conditions. As very little is known of the mechanism(s) controlling PlGF expression in pathology we investigated the role of hyperglycaemia in the regulation of PlGF production in endothelial cells. Hyperglycaemia stimulated PlGF secretion in cultured primary endothelial cells, which was suppressed by IGF-1-mediated PI3K/Akt activation. Inhibition of PI3K activity resulted in significant PlGF mRNA up-regulation and protein secretion. Similarly, loss or inhibition of Akt activity significantly increased basal PlGF expression and prevented any further PlGF secretion in hyperglycaemia. Conversely, constitutive Akt activation blocked PlGF secretion irrespective of upstream PI3K activity demonstrating that Akt is a central regulator of PlGF expression. Knock-down of the Forkhead box O-1 (FOXO1) transcription factor, which is negatively regulated by Akt, suppressed both basal and hyperglycaemia-induced PlGF secretion, whilst FOXO1 gain-of-function up-regulated PlGF in vitro and in vivo. FOXO1 association to a FOXO binding sequence identified in the PlGF promoter also increased in hyperglycaemia. This study identifies the PI3K/Akt/FOXO1 signalling axis as a key regulator of PlGF expression and unifying pathway by which PlGF may contribute to common disorders characterised by endothelial dysfunction, providing a target for therapy.

ETS-dependent enhancers for endothelial-specific expression of serum/glucocorticoid-regulated kinase 1 during mouse embryo development

Serum/glucocorticoid-regulated kinase 1 (SGK1) is predominantly expressed in endothelial cells of mouse embryos, and Sgk1 null mice show embryonic lethality due to impaired vascular formation. However, how the SGK1 expression is controlled in developing vasculature remains unknown. In this study, we first identified a proximal endothelial enhancer through lacZ reporter mouse analyses. The mouse Sgk1 proximal enhancer was narrowed down to the 5' region of the major transcription initiation site, while a human corresponding region possessed relatively weak activity. We then searched for distal enhancer candidates using in silico analyses of publicly available databases for DNase accessibility, RNA polymerase association and chromatin modification. A region approximately 500 kb distant from the human SGK1 gene was conserved in the mouse, and the mouse and human genomic fragments drove transcription restricted to embryonic endothelial cells. Minimal fragments of both proximal and distal enhancers had consensus binding elements for the ETS transcription factors, which were essential for the responsiveness to ERG, FLI1 and ETS1 proteins in luciferase assays and the endothelial lacZ reporter expression in mouse embryos. These results suggest that endothelial SGK1 expression in embryonic vasculature is maintained through at least two ETS-regulated enhancers located in the proximal and distal regions.

Novel insights for the role of nitric oxide in placental vascular function during and beyond pregnancy

More than 30 years have passed since endothelial nitric oxide synthesis was described using the umbilical artery and vein endothelium. That seminal report set the cornerstone for unveiling the molecular aspects of endothelial function. In parallel, the understanding of placental physiology has gained growing interest, due to its crucial role in intrauterine development, with considerable long-term health consequences. This review discusses the evidence for nitric oxide (NO) as a critical player of placental development and function, with a special focus on endothelial nitric oxide synthase (eNOS) vascular effects. Also, the regulation of eNOS-dependent vascular responses in normal pregnancy and pregnancy-related diseases and their impact on prenatal and postnatal vascular health are discussed. Recent and compelling evidence has reinforced that eNOS regulation results from a complex network of processes, with novel data concerning mechanisms such as mechano-sensing, epigenetic, posttranslational modifications, and the expression of NO- and l-arginine-related pathways. In this regard, most of these mechanisms are expressed in an arterial-venous-specific manner and reflect traits of the fetal systemic circulation. Several studies using umbilical endothelial cells are not aimed to understand placental function but general endothelial function, reinforcing the influence of the placenta on general knowledge in physiology.

Specific arterio-venous transcriptomic and ncRNA-RNA interactions in human umbilical endothelial cells: A meta-analysis

Whether arterial-venous differences of primary endothelial cells commonly used for vascular research are preserved in vitro remains under debate. To address this issue, a meta-analysis of Affymetrix transcriptomic data sets from human umbilical artery (HUAECs) and vein (HUVEC) endothelial cells was performed. The meta-analysis showed 2,742 transcripts differentially expressed (false discovery rate <0.05), of which 78% were downregulated in HUVECs. Comparisons with RNA-seq data sets showed high levels of agreement and correlation (p < 0.0001), identifying 84 arterial-venous identity markers. Functional analysis revealed enrichment of key vascular processes in HUAECs/HUVECs, including nitric oxide- (NO) and hypoxia-related genes, as well as differences in miRNA- and ncRNA-mRNA interaction profiles. A proof of concept of these findings in primary cells exposed to hypoxia in vitro and in vivo confirmed the arterial-venous differences in NO-related genes and miRNAs. Altogether, these data defined a cross-platform arterial-venous transcript profile for cultured HUAEC-HUVEC and support a preserved identity involving key vascular pathways post-transcriptionally regulated in vitro.

•

Transcriptional differences among HUAEC and HUVEC are preserved in culture

•

These differences occur even after correcting for experimental conditions

•

The heterogenous regulation affects NO- and hypoxia-related genes

•

Cell-specific ncRNA/mRNA interactions are found

ETS factors are required but not sufficient for specific patterns of enhancer activity in different endothelial subtypes

Correct vascular differentiation requires distinct patterns of gene expression in different subtypes of endothelial cells. Members of the ETS transcription factor family are essential for the transcriptional activation of arterial and angiogenesis-specific gene regulatory elements, leading to the hypothesis that they play lineage-defining roles in arterial and angiogenic differentiation directly downstream of VEGFA signalling. However, an alternative explanation is that ETS binding at enhancers and promoters is a general requirement for activation of many endothelial genes regardless of expression pattern, with subtype-specificity provided by additional factors. Here we use analysis of Ephb4 and Coup-TFII (Nr2f2) vein-specific enhancers to demonstrate that ETS factors are equally essential for vein, arterial and angiogenic-specific enhancer activity patterns. Further, we show that ETS factor binding at these vein-specific enhancers is enriched by VEGFA signalling, similar to that seen at arterial and angiogenic enhancers. However, while arterial and angiogenic enhancers can be activated by VEGFA in vivo, the Ephb4 and Coup-TFII venous enhancers are not, suggesting that the specificity of VEGFA-induced arterial and angiogenic enhancer activity occurs via non-ETS transcription factors. These results support a model in which ETS factors are not the primary regulators of specific patterns of gene expression in different endothelial subtypes.

•

Vein-specific enhancers can contain essential ETS motifs.

•

VEGFA induced an increase in ETS binding at vein, arterial and angiogenic enhancers.

•

VEGFA stimulation cannot induce vein-specific enhancer activity.

Single-cell RNA sequencing reveals regulation of fetal ovary development in the monkey (Macaca fascicularis)

Germ cells are vital for reproduction and heredity. However, the mechanisms underlying female germ cell development in primates, especially in late embryonic stages, remain elusive. Here, we performed single-cell RNA sequencing of 12,471 cells from whole fetal ovaries, and explored the communications between germ cells and niche cells. We depicted the two waves of oogenesis at single-cell resolution and demonstrated that progenitor theca cells exhibit similar characteristics to Leydig cells in fetal monkey ovaries. Notably, we found that ZGLP1 displays differentially expressed patterns between mouse and monkey, which is not overlapped with NANOG in monkey germ cells, suggesting its role in meiosis entry but not in activating oogenic program in primates. Furthermore, the majority of germ cell clusters that sharply express PRDM9 and SPO11 might undergo apoptosis after cyst breakdown, leading to germ cell attrition. Overall, our work provides new insights into the molecular and cellular basis of primate fetal ovary development at single-cell resolution.

Arterialization requires the timely suppression of cell growth

The formation of arteries is thought to occur by the induction of a highly conserved arterial genetic programme in a subset of vessels that will later experience an increase in oxygenated blood flow^1,2. The initial steps of arterial specification require both the VEGF and Notch signalling pathways^3–5. Here, we combine inducible genetic mosaics and transcriptomics to modulate and define the function of these signalling pathways in cell proliferation, arteriovenous differentiation and mobilization. We show that endothelial cells with high levels of VEGF or Notch signalling are intrinsically biased to mobilize and form arteries; however, they are not genetically pre-determined, and can also form veins. Mechanistically, we found that increased levels of VEGF and Notch signalling in pre-arterial capillaries suppresses MYC-dependent metabolic and cell-cycle activities, and promotes the incorporation of endothelial cells into arteries. Mosaic lineage-tracing studies showed that endothelial cells that lack the Notch–RBPJ transcriptional activator complex rarely form arteries; however, these cells regained the ability to form arteries when the function of MYC was suppressed. Thus, the development of arteries does not require the direct induction of a Notch-dependent arterial differentiation programme, but instead depends on the timely suppression of endothelial cell-cycle progression and metabolism, a process that precedes arterial mobilization and complete differentiation.

MiR-21-5p directly contributes to regulating eNOS expression in human artery endothelial cells under normoxia and hypoxia

Clinical conditions associated with hypoxia and oxidative stress, such as fetal growth restriction (FGR), results in endothelial dysfunction. Previous reports show that changes in eNOS expression under these conditions are tightly controlled by DNA methylation and histone posttranslational modifications. However, the contribution of an orchestrating epigenetic mechanism, such as miRNAs, on the NO-related genes expression has not been addressed. We aimed to determine the levels of miRNAs highly expressed in normal endothelial cells (EC), miR-21 and miR-126, in FGR human umbilical artery EC (HUAEC), and their effects on hypoxia-dependent regulation of both, NO-related and oxidative stress-related genes. Results were validated by transcriptome analysis of HUAEC cultured under chronic low oxygen conditions. Cultured FGR-HUAEC showed decreased hsa-miR-21, DDAH1, SOD1, and NRF2, but increased miR-126, NOX4, and eNOS levels, compared with controls. MiR-21-5p levels in FGR were associated with increased hg-miR-21 gene promoter methylation, with no changes in hg-miR-126 gene promoter methylation. HUAEC exposed to hypoxia showed a transient increase in eNOS and DDAH11, paralleled by decrease miR-21-5p levels, but no changes in miR-126-3p and the other genes under study. Transcriptome profiling showed an inverse relationship among miR-21 and several transcripts targeted by miR-21 in HUAEC exposed to hypoxia, meanwhile miR-21-5p-mimic decreased eNOS and DDAH1 transcripts stability, blocking their induction by hypoxia. Consequently, FGR programs a hypoxia-related miRNA that contributes to the regulation of the NO pathway, involving a direct effect of miR-21-5p on eNOS transcript stability, not previously reported. Moreover, hypoxia downregulates miR-21-5p, contributing to increasing the expression of NO-related genes in arterial endothelial cells.

DolphinNext: a distributed data processing platform for high throughput genomics

BACKGROUND

The emergence of high throughput technologies that produce vast amounts of genomic data, such as next-generation sequencing (NGS) is transforming biological research. The dramatic increase in the volume of data, the variety and continuous change of data processing tools, algorithms and databases make analysis the main bottleneck for scientific discovery. The processing of high throughput datasets typically involves many different computational programs, each of which performs a specific step in a pipeline. Given the wide range of applications and organizational infrastructures, there is a great need for highly parallel, flexible, portable, and reproducible data processing frameworks. Several platforms currently exist for the design and execution of complex pipelines. Unfortunately, current platforms lack the necessary combination of parallelism, portability, flexibility and/or reproducibility that are required by the current research environment. To address these shortcomings, workflow frameworks that provide a platform to develop and share portable pipelines have recently arisen. We complement these new platforms by providing a graphical user interface to create, maintain, and execute complex pipelines. Such a platform will simplify robust and reproducible workflow creation for non-technical users as well as provide a robust platform to maintain pipelines for large organizations.

RESULTS

To simplify development, maintenance, and execution of complex pipelines we created DolphinNext. DolphinNext facilitates building and deployment of complex pipelines using a modular approach implemented in a graphical interface that relies on the powerful Nextflow workflow framework by providing 1. A drag and drop user interface that visualizes pipelines and allows users to create pipelines without familiarity in underlying programming languages. 2. Modules to execute and monitor pipelines in distributed computing environments such as high-performance clusters and/or cloud 3. Reproducible pipelines with version tracking and stand-alone versions that can be run independently. 4. Modular process design with process revisioning support to increase reusability and pipeline development efficiency. 5. Pipeline sharing with GitHub and automated testing 6. Extensive reports with R-markdown and shiny support for interactive data visualization and analysis.

CONCLUSION

DolphinNext is a flexible, intuitive, web-based data processing and analysis platform that enables creating, deploying, sharing, and executing complex Nextflow pipelines with extensive revisioning and interactive reporting to enhance reproducible results.