Coronary blood vessels from distinct origins converge to equivalent states during mouse and human development

Ragini Phansalkar, Josephine Krieger, Mingming Zhao, Sai Saroja Kolluru, Robert C. Jones, Stephen R Quake, Irving Weissman, Daniel Bernstein, Virginia D. Winn, Gaetano D’Amato, Kristy Red-Horse

Preprint posted on 26 April 2021 https://www.biorxiv.org/content/10.1101/2021.04.25.441380v1

Article now published in eLife at http://dx.doi.org/10.7554/eLife.70246

All endothelial roads lead to “Rome”: understanding the cell plasticity of cardiac blood vessels

Selected by Yen Tran, Osvaldo Contreras

Categories: bioinformatics, cell biology, developmental biology

Background

The coronary arteries branch along the coronary groove of the heart and supply oxygenated blood to the myocardium and heart cells. Owing to their essential function, coronary blood vessel formation and maintenance regulate the normal functions of the cardiac muscle. When this “goes wrong,” for example, in diseased coronary blood vessels, cardiac function is severely affected. Coronary diseases (also known as coronary artery disease, or CAD) are rising in the industrialized world, representing the most common cause of heart disease in the United States. Because CAD severely impacts the quality of life of patients and has a significant socioeconomic impact (WHO), understanding the underlying mechanisms of blood vessel development and regulation in homeostasis and disease could open roads for novel and efficient medical applications.

How do coronary vessels develop? The formation of coronary vessels requires two principal cellular sources: sinus venosus (SV)- and endocardium (Endo)-derived endothelial cells (Lupu et al.). The sinus venosus is the inflow tract of the developing heart, and the endocardium is the layer lining inside the cardiac ventricles. At embryonic day E11.5, endothelial cells from SV migrate caudally and contribute to the vessel network in the outer myocardial wall. Conversely, endocardial cells extrapolate inside out and contribute to the vessels network of inner myocardial walls and the septum. Later, the two networks unite and remodel into mature capillaries, arteries and veins before connecting with the aorta and carrying blood flow. Previously, Su et al. used single-cell RNA sequencing (scRNA-seq) to characterize the transcriptomes of SV-derived endothelial cells. However, a direct comparison between Endo- and SV- derived ECs during coronary vessel formation has not been explored.

What is the biological question addressed in this preprint, and why is it important? In an exciting new preprint, Phansalkar et al. take advantage of scRNA-seq combined with endothelial cell lineage-tracing to explore for the first time the dual origins of coronary blood vessels during mouse and human development.

Key findings

To address the origin(s) of endothelial cells during coronary vessel formation, Phansalkar and colleagues integrated single-cell RNA sequencing with lineage tracing of the endocardium (Endo) and sinus venosus (SV) derived endothelial cells (ECs) using BmxCreER and ApjCreER mice, respectively. The authors collected ECs from two embryonic (E12 and E17.5) and adult stages (8 weeks of age). The single-cell data captured most endothelial cell types, and as expected, the proportion of each cell type changed over the developmental course. The E12 coronary vessels contained two capillary plexus (Cap1, Cap2) and one pre-artery (pre-Art) cluster. Later in development, E17.5 coronary ECs resolved two capillaries (Cap1 and Cap2), two arteries (Art1 & Art2), and one vein cluster. The authors found one artery, one vein, and two capillary clusters in adult coronary ECs.

Further characterization of the cell lineages of coronary ECs suggests that the origins and heterogeneity of SV and Endo fade over time. At E12, the Cap1 subpopulation shared molecular signatures with Endo, whereas the Cap2 subpopulation was enriched with SV marker genes. However, these patterns are not preserved at E17.5 and in adults. Remarkably, although E17.5 capillary clusters did not retain their origin’s signature, they were segregated by their location within the heart. Cap2 cells (Kcne3- and Car4+) are associated with the left and right ventricle free walls and dorsal side. Cap1 cells (Kcne3+ and Car4) are located in the septum and ventral side.

Interestingly, differential gene expression of Cap1 and Cap2 reflected the distinct features of cardiac regions. Cap1 cells had high expression of hypoxia-related genes and tip markers, implying that the septum region may be hypoxic and require new blood vessel formation. In addition, Cap2 cells highly expressed shear stress-induced genes, suggesting that this region experiences a higher blood flow rate. Furthermore, most of the Cap1 cells in the septum and ventral areas are tdTomato+, indicating these cells derive from Endo lineage. In the adult, both lineages- and location-based cell heterogeneity disappeared.

Owing to Endo- and SV-derived endothelial cells having similar transcriptomic profiles in adult hearts, the authors sought to understand whether the cellular properties of Endo- and SV-derived ECs change under injury. Although endothelial cells proliferate at high rates, no significant differences were observed in the S-phase progression among both endothelial cell lineages, as evaluated by the incorporation of a nucleotide analog EdU. These findings suggest that endothelial lineage does not influence cell proliferation. In this context, it would be interesting to see the scRNA-seq following cardiac injury using the Bmx lineage-tracing in adult mice.

Next, using scRNA-seq analyses, the authors explored whether mouse and human coronary vessels share similar transcriptomic features, cellular states, and trajectories at different developmental stages. Smart-seq2 of two prominent human blood vessel fractions (CD31+ and CD36+) were combined and compared to the mouse scRNA-seq data through Label Transfer analysis. These results found shared developmental cues and endothelial cell lineages and states between the two species. Similar findings were found for the capillary-to-artery transition, as mouse endothelial cells mapped to the different human artery subclusters described by the authors. Remarkably, trajectory analyses of murine arterial cellular states suggested the progression of a single immature state -not found in adult artery endothelial cells- into two continuous and mature states, which also exist in the adult human heart. Using trajectory analysis, the authors also suggested that capillary endothelial cells give rise to arteries, as previously described in mouse heart development. The capillary-to-artery transition could be happening simultaneously at two different heart regions, reinforcing their previous findings concerning the defining role of heart location or microenvironment on blood vessel formation. However, detailed spatiotemporal studies of these cell subpopulations should corroborate these findings.

Motivated by the fact that coronary artery disease is a significant cause of death in the industrialized world, the authors closely analyzed the transcriptomes of the three different arterial cells during development, aiming to find potential genes that could shed light on disease progression and biomarkers. They found striking differences in their maturation states that correlated with differential gene expression and progression of artery arborization. In addition, their findings also highlighted the potential role of newly discovered transcription factors in Art subpopulations, including PRDM16, GATA2, and IRF6. Human-enriched genes involved in neurotransmission (e.g., GABA and Glutamate receptors) were not expressed in the mouse. The latter findings could open the way into potential new genes to target coronary diseases.

What we liked about this preprint

Most bioinformatic analyses of scRNA-seq data focus on the divergence state of cells during commitment and differentiation. Interestingly, here the authors use single-cell analysis and lineage tracing to address an opposite process, called convergent differentiation. In this context, two or more lineage origins contribute to the same endpoint cell state. Here, the authors addressed the converging state of endocardium and sinus venous-derived ECs onto capillary plexus. They found that the descendant capillary cells ‘remember’ their origins at the beginning of convergence but lose their origin memories during the converge trajectory and adopt location-based molecular properties when residing in different cardiac regions. By the end of this process, all capillary cells become indistinguishable by origin and location.

In conclusion, this preprint provides complete mouse and human single-cell datasets from the initiation of coronary vessel formation at embryonic stages to adulthood. The similarities between mouse and human coronary ECs single-cell data validate the usability of the mouse to understand human blood vessel development and speculatively might provide therapeutic windows to the high prevalence of coronary artery disease.

Future directions and questions to the authors

Intriguingly, the author chose 10X Genomics for coronary ECs at embryonic days E12 and E17.5 and Smart-seq2 for the remaining data. Why these different single-cell technologies? Given that Smart-seq2 captures full-length mRNA while 10X Genomics uses a 5’ or 3’-tag sequencing method, could their parallel use introduce bias in the number of genes recovered and subsequent comparisons?
The authors suggested that coronary capillary clusters at E12 are distinguished based on either SV- or Endo-derived gene signatures. What criteria did you use for selecting these gene signatures for the classification of Cap1 and Cap2?
In figure 3g, the authors quantified the number of Car4+ cells associated with tdTomato+ lineage and distinct cardiac regions at E17.5. The results clearly show an enrichment of tdTomato+/Car4- in the septum and ventral cardiac areas, whereas this location-based heterogeneity is lost in the adult. Do you have similar quantitative results in the adult?
Are there any differences in the proliferative or cell cycle progression capabilities of endothelial cell lineages (Endo- vs SV-derived) during heart development?
What do the authors think about the existence of different subclusters of immature (one) and mature (two) arteries with the absence or presence of smooth muscle? At what extent endothelial cell developmental trajectories are determined by smooth muscle cells or vice versa?

References

Lupu, IE., De Val, S. & Smart, N. Coronary vessel formation in development and disease: mechanisms and insights for therapy. Nat Rev Cardiol 17, 790–806 (2020). https://doi.org/10.1038/s41569-020-0400-1

Su, T., Stanley, G., Sinha, R. et al. Single-cell analysis of early progenitor cells that build coronary arteries. Nature 559, 356–362 (2018). https://doi.org/10.1038/s41586-018-0288-7

World Health Organization (WHO): https://www.who.int/en/news-room/fact-sheets/detail/cardiovascular-diseases-(cvds). Accessed June the 3rd, 2021.

Tags: cell differentiation, coronary artery disease, coronary vessels, development, endothelial cells, single-cell rnaseq

Posted on: 4 June 2021 , updated on: 7 June 2021

doi: https://doi.org/10.1242/prelights.29382

Read preprint

(No Ratings Yet)

Author's response

Kristy Red-Horse shared

Question 1: Intriguingly, the author chose 10X Genomics for coronary ECs at embryonic days E12 and E17.5 and Smart-seq2 for the remaining data. Why these different single-cell technologies? Given that Smart-seq2 captures full-length mRNA while 10X Genomics uses a 5’ or 3’-tag sequencing method, could their parallel use introduce bias in the number of genes recovered and subsequent comparisons?

While using 10X Genomics is a cost-effective way to perform large-scale single-cell RNA sequencing studies, we chose to use Smart-seq2 for our adult mouse and human datasets for two reasons:

First, we knew from preliminary analysis of datasets including Tabula Muris that coronary endothelial cells, particularly capillary endothelial cells (ECs), would be more homogenous in adults than in embryos. Thus, we chose to use Smart-seq2, which is a more sensitive technology (Ding et al, Nature Biotechnology, 2020) in order to detect even smaller differences between sinus venosus (SV) and endocardium (Endo) derived cells in adults. Nevertheless, we still did not find any significant differences between the lineages. However, we did see about 10 genes differentially expressed between the two adult capillary clusters, and these genes were also differentially expressed between the two capillary clusters at e17.5 (one of our 10X datasets). This indicates that some small remnant of e17.5 heterogeneity can be detected in adults, and also that the two different technologies used were able to detect a similar set of differentially expressed genes. Finally, since we did not make direct comparisons of the stages, but instead compared lineages and clusters within each individual dataset, we do not believe that bias in the methods used significantly influenced our overall conclusions.

Second, using Smart-seq2 to produce a dataset of developing human coronary ECs allows us to maximize the new information gained from these rare samples. We were able to produce as a resource for the field a very high-quality dataset with coverage of up to >1 million reads/cell and up to 6000 genes/cell. Additionally, because Smart-seq2 captures full-length cDNA, this will enable future studies of alternative splicing in these cells.

Question 2: The authors suggested that coronary capillary clusters at E12 are distinguished based on either SV- or Endo-derived gene signatures. What criteria did you use for selecting these gene signatures for the classification of Cap1 and Cap2?

The SV and Endo gene signatures were determined by making a direct comparison of the SV and Endo clusters in the e12 dataset (see Fig S2 and S6), and identifying differentially expressed genes that were enriched in either the SV or Endo and passed some thresholds as specified in the methods.

Question 3: In figure 3g, the authors quantified the number of Car4+ cells associated with tdTomato+ lineage and distinct cardiac regions at E17.5. The results clearly show an enrichment of tdTomato+/Car4- in the septum and ventral cardiac areas, whereas this location-based heterogeneity is lost in the adult. Do you have similar quantitative results in the adult?

We found an average of 87% Car4+ ECs in the LV, and an average of 84% Car4+ ECs in the septum (n=2).

Question 4: Are there any differences in the proliferative or cell cycle progression capabilities of endothelial cell lineages (Endo- vs SV-derived) during heart development?

We focused on non-proliferating cells for our direct comparison of the SV and Endo lineages in order to avoid any unwanted sources of variation. However, we did note that at e12, the majority of SV-derived ECs are proliferating while the Endo-derived ECs are largely non-proliferating. This could be due to the Endo-derived cells being relatively “immature” as vascular ECs compared to the SV-derived cells, which already started off as part of a vessel. It could also be that Endo-derived ECs need to complete their transition away from the Endo before they can begin proliferating. By e17.5, there is no difference in the percentage of proliferating cells either between the SV and Endo lineages (43% and 42%, respectively), or between Cap1 and Cap2 and including or excluding proliferating cells in the analysis yields the same results (see Fig S4).

Question 5: What do the authors think about the existence of different subclusters of immature (one) and mature (two) arteries with the absence or presence of smooth muscle? At what extent endothelial cell developmental trajectories are determined by smooth muscle cells or vice versa?

In both developing mouse and human, we see multiple clusters of artery ECs that seem to correspond to the maturity of the artery cells, in that one cluster (the most mature) is expressing GJA5 (Cx40) and the others are not. In human, we also found that the GJA5-expressing cells (which also uniquely express JAG1) are generally in the larger, more proximal arteries, while the GJA5-negative artery ECs are generally in smaller, more distal arteries. Both of these can be covered by smooth muscle, although some of the GJA5-negative vessels are not. Whether all smooth-muscle covered arteries eventually become GJA5-positive in adult humans is unknown, as is the extent of any difference in remodeling or pruning between arteries that are or are not covered in smooth muscle at these mid-gestation stages. Previous work by Drs. Amber Stratman and Brant Weinstein in zebrafish has shown that vascular smooth muscle cells (vSMCs) can stabilize vessels (Stratman et al, Development, 2017) and that expression of Cxcr4 in artery ECs is part of a pathways that recruits vSMCs (Stratman et al, Communications Biology, 2020). In our human dataset, CXCR4 is most highly expressed in Art3, our immature artery cluster, suggesting that this could represent a set of newly-formed arteries in the process of vSMC recruitment.

Have your say Cancel reply

This site uses Akismet to reduce spam. Learn how your comment data is processed.

Sign up to customise the site to your preferences and to receive alerts

Also in the bioinformatics category:

Expressive modeling and fast simulation for dynamic compartments

Till Köster, Philipp Henning, Tom Warnke, et al.

Selected by Benjamin Dominik Maier

Transcriptional profiling of human brain cortex identifies novel lncRNA-mediated networks dysregulated in amyotrophic lateral sclerosis

Alessandro Palma, Monica Ballarino

Selected by Julio Molina Pineda

Spatial transcriptomics elucidates medulla niche supporting germinal center response in myasthenia gravis thymoma

Yoshiaki Yasumizu, Makoto Kinoshita, Martin Jinye Zhang, et al.

Selected by Jessica Chevallier

Also in the cell biology category:

Fetal brain response to maternal inflammation requires microglia

Bridget Elaine LaMonica Ostrem, Nuria Dominguez Iturza, Jeffrey Stogsdill, et al.

Selected by Manuel Lessi

Alteration of long and short-term hematopoietic stem cell ratio causes myeloid-biased hematopoiesis

Katsuyuki Nishi, Taro Sakamaki, Akiomi Nagasaka, et al.

Selected by Manuel Vicente, Anaisa Ferreira

Clusters of lineage-specific genes are anchored by ZNF274 in repressive perinucleolar compartments

Martina Begnis, Julien Duc, Sandra Offner, et al.

Selected by Silvia Carvalho

Also in the developmental biology category:

Alteration of long and short-term hematopoietic stem cell ratio causes myeloid-biased hematopoiesis

Katsuyuki Nishi, Taro Sakamaki, Akiomi Nagasaka, et al.

Selected by Manuel Vicente, Anaisa Ferreira

Temporal constraints on enhancer usage shape the regulation of limb gene transcription

Raquel Rouco, Antonella Rauseo, Guillaume Sapin, et al.

Selected by María Mariner-Faulí

OGT prevents DNA demethylation and suppresses the expression of transposable elements in heterochromatin by restraining TET activity genome-wide

Hugo Sepulveda, Xiang Li, Xiaojing Yue, et al.

Selected by Mansi

preLists in the bioinformatics category:

‘In preprints’ from Development 2022-2023

A list of the preprints featured in Development's 'In preprints' articles between 2022-2023

Coronary blood vessels from distinct origins converge to equivalent states during mouse and human development

Background

Key findings

What we liked about this preprint

Future directions and questions to the authors

References

Share this:

Have your say Cancel reply

Sign up to customise the site to your preferences and to receive alerts

Also in the bioinformatics category:

Expressive modeling and fast simulation for dynamic compartments

Transcriptional profiling of human brain cortex identifies novel lncRNA-mediated networks dysregulated in amyotrophic lateral sclerosis

Spatial transcriptomics elucidates medulla niche supporting germinal center response in myasthenia gravis thymoma

Also in the cell biology category:

Fetal brain response to maternal inflammation requires microglia

Alteration of long and short-term hematopoietic stem cell ratio causes myeloid-biased hematopoiesis

Clusters of lineage-specific genes are anchored by ZNF274 in repressive perinucleolar compartments

Also in the developmental biology category:

Alteration of long and short-term hematopoietic stem cell ratio causes myeloid-biased hematopoiesis

Temporal constraints on enhancer usage shape the regulation of limb gene transcription

OGT prevents DNA demethylation and suppresses the expression of transposable elements in heterochromatin by restraining TET activity genome-wide

preLists in the bioinformatics category:

‘In preprints’ from Development 2022-2023

9th International Symposium on the Biology of Vertebrate Sex Determination

Alumni picks – preLights 5th Birthday

Fibroblasts

Single Cell Biology 2020

Antimicrobials: Discovery, clinical use, and development of resistance

Also in the cell biology category:

BSCB-Biochemical Society 2024 Cell Migration meeting

‘In preprints’ from Development 2022-2023

preLights peer support – preprints of interest

The Society for Developmental Biology 82nd Annual Meeting

CSHL 87th Symposium: Stem Cells

Journal of Cell Science meeting ‘Imaging Cell Dynamics’

9th International Symposium on the Biology of Vertebrate Sex Determination

Alumni picks – preLights 5th Birthday

CellBio 2022 – An ASCB/EMBO Meeting

Fibroblasts

EMBL Synthetic Morphogenesis: From Gene Circuits to Tissue Architecture (2021)

FENS 2020

Planar Cell Polarity – PCP

BioMalPar XVI: Biology and Pathology of the Malaria Parasite

Cell Polarity

TAGC 2020

3D Gastruloids

ECFG15 – Fungal biology

ASCB EMBO Annual Meeting 2019

EMBL Seeing is Believing – Imaging the Molecular Processes of Life

Autophagy

Lung Disease and Regeneration

Cellular metabolism

BSCB/BSDB Annual Meeting 2019

MitoList

Biophysical Society Annual Meeting 2019

ASCB/EMBO Annual Meeting 2018

Also in the developmental biology category:

GfE/ DSDB meeting 2024

‘In preprints’ from Development 2022-2023

preLights peer support – preprints of interest

The Society for Developmental Biology 82nd Annual Meeting

CSHL 87th Symposium: Stem Cells

Journal of Cell Science meeting ‘Imaging Cell Dynamics’

9th International Symposium on the Biology of Vertebrate Sex Determination

Alumni picks – preLights 5th Birthday

CellBio 2022 – An ASCB/EMBO Meeting

2nd Conference of the Visegrád Group Society for Developmental Biology

Fibroblasts

EMBL Synthetic Morphogenesis: From Gene Circuits to Tissue Architecture (2021)

EMBL Conference: From functional genomics to systems biology

Single Cell Biology 2020

Society for Developmental Biology 79th Annual Meeting

FENS 2020

Planar Cell Polarity – PCP

Cell Polarity

TAGC 2020

3D Gastruloids

ASCB EMBO Annual Meeting 2019

EDBC Alicante 2019

EMBL Seeing is Believing – Imaging the Molecular Processes of Life