Single-Cell RNA Sequencing Resolves Spatiotemporal Development of Pre-thymic Lymphoid Progenitors and Thymus Organogenesis in Human Embryos.

Generation of the first T lymphocytes in the human embryo involves the emergence, migration, and thymus seeding of lymphoid progenitors together with concomitant thymus organogenesis, which is the initial step to establish the entire adaptive immune system. However, the cellular and molecular programs regulating this process remain unclear. We constructed a single-cell transcriptional landscape of human early T lymphopoiesis by using cells from multiple hemogenic and hematopoietic sites spanning embryonic and fetal stages. Among heterogenous early thymic progenitors, one subtype shared common features with a subset of lymphoid progenitors in fetal liver that are known as thymus-seeding progenitors. Unbiased bioinformatics analysis identified a distinct type of pre-thymic lymphoid progenitors in the aorta-gonad-mesonephros (AGM) region. In parallel, we investigated thymic epithelial cell development and potential cell-cell interactions during thymus organogenesis. Together, our data provide insights into human early T lymphopoiesis that prospectively direct T lymphocyte regeneration, which might lead to development of clinical applications.

Deciphering human macrophage development at single-cell resolution.

Macrophages are the first cells of the nascent immune system to emerge during embryonic development. In mice, embryonic macrophages infiltrate developing organs, where they differentiate symbiotically into tissue-resident macrophages (TRMs)1. However, our understanding of the origins and specialization of macrophages in human embryos is limited. Here we isolated CD45+ haematopoietic cells from human embryos at Carnegie stages 11 to 23 and subjected them to transcriptomic profiling by single-cell RNA sequencing, followed by functional characterization of a population of CD45+CD34+CD44+ yolk sac-derived myeloid-biased progenitors (YSMPs) by single-cell culture. We also mapped macrophage heterogeneity across multiple anatomical sites and identified diverse subsets, including various types of embryonic TRM (in the head, liver, lung and skin). We further traced the specification trajectories of TRMs from either yolk sac-derived primitive macrophages or YSMP-derived embryonic liver monocytes using both transcriptomic and developmental staging information, with a focus on microglia. Finally, we evaluated the molecular similarities between embryonic TRMs and their adult counterparts. Our data represent a comprehensive characterization of the spatiotemporal dynamics of early macrophage development during human embryogenesis, providing a reference for future studies of the development and function of human TRMs.

Transcriptomic landscape of circulating mononuclear phagocytes in Langerhans cell histiocytosis at the single-cell level.

Langerhans cell histiocytosis (LCH) is an inflammatory myeloid neoplasm caused by aberrant activation of the mitogen-activated protein kinase (MAPK) pathway. Circulating myeloid cells from patients often carry disease-associated mutations and can be differentiated into langerinhigh LCH-like cells in vitro, but their detailed immune-phenotypic and molecular profiles are lacking and could shed key insights into disease biology. Here we recruited 217 pediatric LCH patients and took blood and tissue samples for BRAFV600E analysis. Immune-phenotyping of the circulating Lin-HLA-DR+ immune population in 49 of these patients revealed that decreased frequency of plasmacytoid dendritic cells was significantly linked to disease severity. By single-cell RNA sequencing of samples from 14 patients, we identified key changes in expression of RAS-MAPK-extracellular signal-regulated kinase (ERK) signaling-related genes and transcription factors in distinct members of the mononuclear phagocyte system in the presence of BRAFV600E. Moreover, treatment of patients with the BRAF inhibitor dabrafenib resulted in MAPK cascade inhibition, inflammation prevention, and regulation of cellular metabolism within mononuclear phagocytes. Finally, we also observed elevated expression of RAS-MAPK-ERK signaling-related genes in a CD207+CD1a+ cell subcluster in skin. Taken together, our data extend the molecular understanding of LCH biology at single-cell resolution, which might contribute to improvement of clinical diagnostics and therapeutics, and aid in the development of personalized medicine approaches.

Tracing haematopoietic stem cell formation at single-cell resolution.

Haematopoietic stem cells (HSCs) are derived early from embryonic precursors, such as haemogenic endothelial cells and pre-haematopoietic stem cells (pre-HSCs), the molecular identity of which still remains elusive. Here we use potent surface markers to capture the nascent pre-HSCs at high purity, as rigorously validated by single-cell-initiated serial transplantation. Then we apply single-cell RNA sequencing to analyse endothelial cells, CD45(-) and CD45(+) pre-HSCs in the aorta-gonad-mesonephros region, and HSCs in fetal liver. Pre-HSCs show unique features in transcriptional machinery, arterial signature, metabolism state, signalling pathway, and transcription factor network. Functionally, activation of mechanistic targets of rapamycin (mTOR) is shown to be indispensable for the emergence of HSCs but not haematopoietic progenitors. Transcriptome data-based functional analysis reveals remarkable heterogeneity in cell-cycle status of pre-HSCs. Finally, the core molecular signature of pre-HSCs is identified. Collectively, our work paves the way for dissection of complex molecular mechanisms regulating stepwise generation of HSCs in vivo, informing future efforts to engineer HSCs for clinical applications.

Dissecting transcriptional heterogeneity in primary gastric adenocarcinoma by single cell RNA sequencing.

OBJECTIVE: Tumour heterogeneity represents a major obstacle to accurate diagnosis and treatment in gastric adenocarcinoma (GA). Here, we report a systematic transcriptional atlas to delineate molecular and cellular heterogeneity in GA using single-cell RNA sequencing (scRNA-seq). DESIGN: We performed unbiased transcriptome-wide scRNA-seq analysis on 27 677 cells from 9 tumour and 3 non-tumour samples. Analysis results were validated using large-scale histological assays and bulk transcriptomic datasets. RESULTS: Our integrative analysis of tumour cells identified five cell subgroups with distinct expression profiles. A panel of differentiation-related genes reveals a high diversity of differentiation degrees within and between tumours. Low differentiation degrees can predict poor prognosis in GA. Among them, three subgroups exhibited different differentiation grade which corresponded well to histopathological features of Lauren's subtypes. Interestingly, the other two subgroups displayed unique transcriptome features. One subgroup expressing chief-cell markers (eg, LIPF and PGC) and RNF43 with Wnt/ß-catenin signalling pathway activated is consistent with the previously described entity fundic gland-type GA (chief cell-predominant, GA-FG-CCP). We further confirmed the presence of GA-FG-CCP in two public bulk datasets using transcriptomic profiles and histological images. The other subgroup specifically expressed immune-related signature genes (eg, LY6K and major histocompatibility complex class II) with the infection of Epstein-Barr virus. In addition, we also analysed non-malignant epithelium and provided molecular evidences for potential transition from gastric chief cells into MUC6+TFF2+ spasmolytic polypeptide expressing metaplasia. CONCLUSION: Altogether, our study offers valuable resource for deciphering gastric tumour heterogeneity, which will provide assistance for precision diagnosis and prognosis.

Single-cell transcriptomic architecture and intercellular crosstalk of human intrahepatic cholangiocarcinoma.

BACKGROUND & AIMS: Intrahepatic cholangiocarcinoma (ICC) is the second most common liver malignancy. ICC typically features remarkable cellular heterogeneity and a dense stromal reaction. Therefore, a comprehensive understanding of cellular diversity and the interplay between malignant cells and niche cells is essential to elucidate the mechanisms driving ICC progression and to develop therapeutic approaches. METHODS: Herein, we performed single-cell RNA sequencing (scRNA-seq) analysis on unselected viable cells from 8 human ICCs and adjacent samples to elucidate the comprehensive transcriptomic landscape and intercellular communication network. Additionally, we applied a negative selection strategy to enrich fibroblast populations in 2 other ICC samples to investigate fibroblast diversity. The results of the analyses were validated using multiplex immunofluorescence staining, bulk transcriptomic datasets, and functional in vitro and in vivo experiments. RESULTS: We sequenced a total of 56,871 single cells derived from human ICC and adjacent tissues and identified diverse tumor, immune, and stromal cells. Malignant cells displayed a high degree of inter-tumor heterogeneity. Moreover, tumor-infiltrating CD4 regulatory T cells exhibited highly immunosuppressive characteristics. We identified 6 distinct fibroblast subsets, of which the majority were CD146-positive vascular cancer-associated fibroblasts (vCAFs), with highly expressed microvasculature signatures and high levels of interleukin (IL)-6. Functional assays indicated that IL-6 secreted by vCAFs induced significant epigenetic alterations in ICC cells, particularly upregulating enhancer of zeste homolog 2 (EZH2) and thereby enhancing malignancy. Furthermore, ICC cell-derived exosomal miR-9-5p elicited high expression of IL-6 in vCAFs to promote tumor progression. CONCLUSIONS: Our single-cell transcriptomic dataset delineates the inter-tumor heterogeneity of human ICCs, underlining the importance of intercellular crosstalk between ICC cells and vCAFs, and revealing potential therapeutic targets. LAY SUMMARY: Intrahepatic cholangiocarcinoma is an aggressive and chemoresistant malignancy. Better understanding the complex transcriptional architecture and intercellular crosstalk of these tumors will help in the development of more effective therapies. Herein, we have identified important interactions between cancer cells and cancer-associated fibroblasts in the tumor stroma, which could have therapeutic implications.

Effects of blood sample handling procedures on measurable interleukin 6 in plasma and serum.

INTRODUCTION: Interleukin-6(IL-6) measurement is used as a biomarker in medical diagnosis, therapy, and prognosis in various diseases. However, several pre-analytical factors may yield a false IL-6 result. In this study, we set out to investigate the effects of corrected blood sample handling procedures on measurable IL-6. METHOD: EDTA plasma and serum samples were collected from 45 healthy individuals. The participants were divided into three groups to perform different handling procedures. Different centrifugal timing, storage temperature, and time were executed on the samples. The changed trends of IL-6 levels were analyzed. RESULTS: At baseline, while the paired plasma and serum IL-6 values had a good correlation, the plasma levels were higher than serum. In general, the unseparated EDTA plasma kept steady with time. With the increase in storage temperature and time, a more pronounced rise in unseparated serum IL-6 was observed. Nevertheless, the samples in Group 3 which centrifuged and separated immediately kept stable after a different temperature and longtime storage. CONCLUSION: Sample types, centrifugal timing, storage temperature, and time may affect the IL-6 levels. A standard blood sample handling procedure should be performed to ensure the accuracy and stability of IL-6 values.

ATF4 plays a pivotal role in the development of functional hematopoietic stem cells in mouse fetal liver.

The fetal liver (FL) serves as a predominant site for expansion of functional hematopoietic stem cells (HSCs) during mouse embryogenesis. However, the mechanisms for HSC development in FL remain poorly understood. In this study, we demonstrate that deletion of activating transcription factor 4 (ATF4) significantly impaired hematopoietic development and reduced HSC self-renewal in FL. In contrast, generation of the first HSC population in the aorta-gonad-mesonephros region was not affected. The migration activity of ATF4(-/-) HSCs was moderately reduced. Interestingly, the HSC-supporting ability of both endothelial and stromal cells in FL was significantly compromised in the absence of ATF4. Gene profiling using RNA-seq revealed downregulated expression of a panel of cytokines in ATF4(-/-) stromal cells, including angiopoietin-like protein 3 (Angptl3) and vascular endothelial growth factor A (VEGFA). Addition of Angptl3, but not VEGFA, partially rescued the repopulating defect of ATF4(-/-) HSCs in the culture. Furthermore, chromatin immunoprecipitation assay in conjunction with silencing RNA-mediated silencing and complementary DNA overexpression showed transcriptional control of Angptl3 by ATF4. To summarize, ATF4 plays a pivotal role in functional expansion and repopulating efficiency of HSCs in developing FL, and it acts through upregulating transcription of cytokines such as Angptl3 in the microenvironment.

Preparation of poly(vinylphenylboronic acid) chain grafted poly(glycidylmethacrylate-co-ethylenedimethacrylate) beads for the selective enrichment of glycoprotein.

Surface-initiated atom transfer radical polymerization was successfully used to prepare 4-vinylphenylboronic acid functionalized poly(glycidylmethacrylate-co-ethylenedimethacrylate) beads for the selective enrichment of glycoprotein from complex biological samples in this study. The modified bead surfaces were characterized using Fourier transform infrared spectroscopy and X-ray photoelectron spectroscopy. The sorption behaviors, including adsorption isotherms, incubation time, and pH effect, were investigated. The results demonstrated that the boronated beads have a high affinity for glycoprotein, which is due to the well-defined boronic acid brushes on the beads surfaces. Furthermore, the polyvinylphenylboronic acid grafted poly(glycidylmethacrylate-co-ethylenedimethacrylate) beads were used to efficiently enrich and purify glycoprotein from real egg white samples and α-fetoprotein from human serum samples. The mass spectrometry results demonstrated that the polyvinylphenylboronic acid grafted poly(glycidylmethacrylate-co-ethylenedimethacrylate) beads are a suitable material for the enrichment of glycosylated protein from complex biological samples.

Selective solid-phase extraction of artificial chemicals from milk samples using multiple-template surface molecularly imprinted polymers.

A novel multiple-template surface molecularly imprinted polymer (MTMIP) was synthesized using ofloxacin and 17ß-estradiol as templates and modified monodispersed poly(glycidylmethacrylate-co-ethylene dimethacrylate) (PGMA/EDMA ) beads as the support material. Static adsorption, solid-phase extraction and high-performance liquid chromatography were performed to investigate the adsorption properties and selective recognition characteristics of the polymer templates and their structural analogs. The maximum binding capacities of ofloxacin and 17ß-estradiol on the MTMIP were 9.0 and 6.6 mg/g, respectively. Compared with the corresponding nonimprinted polymer, the MTMIP exhibited a much higher adsorption performance and selectivity toward three quinolones and three estrogens, which are common drug residues in food. The MTMIP served as a simple and effective pretreatment method and could be successfully applied to the simultaneous analysis of multiple target components in complex samples. Furthermore, the MTMIP may find useful applications as a solid-phase absorbent in the simultaneous determination of trace quinolones and estrogens in milk samples, as the recoveries were in the range 77.6-98.0%.

Hematopoietic stem cell heterogeneity in non-human primates revealed by five-lineage output bias analysis.

Understanding hematopoietic stem cell (HSC) heterogeneity is crucial for treating malignant blood disorders. Compared with mice, we have limited knowledge of the heterogeneity of human HSCs. Fortunately, non-human primates (NHPs) have become the best animal models for studying human HSCs. Here, we employed a public dataset derived from NHP autologous bone marrow transplantation, and focused on a total of 820 HSC clones with reconstitution capacity of all available five lineages (granulocyte, monocyte, B cell, T cell, and natural killer cell) at two time points (11/12 and/or 42/43 months). Intriguingly, unsupervised clustering on these clones revealed six HSC subtypes, including a lymphoid/myeloid balanced (LM-balanced) subtype and five single-lineage-biased subtypes. We also observed that the subtypes of these HSC clones might change over time, and a given subtype could transition into any one of the other five subtypes, albeit with a certain degree of selectivity. Particularly, each of the six subtypes was more likely to turn into lymphoid-biased rather than myeloid-biased ones. Additionally, our five-lineage classification method exhibited strong correlation with traditional lymphoid/myeloid bias classification method. Specifically, our granulocyte- and monocyte-biased subtypes were predominantly attributed to α-HSCs, while LM-balanced, B cell-biased, and T cell-biased subtypes were primarily associated with ß-HSCs. The γ-HSCs were composed of a small subset of B cell-biased and T cell-biased subtypes. In summary, our five-lineage classification identifies more finely tuned HSC subtypes based on lineage output bias. These findings enrich our understanding of HSC heterogeneity in NHPs and provide important insights for human research.

Human yolk sac-derived innate lymphoid-biased multipotent progenitors emerge prior to hematopoietic stem cell formation.

Hematopoietic stem cell (HSC)-independent lymphopoiesis has been elucidated in murine embryos. However, our understanding regarding human embryonic counterparts remains limited. Here, we demonstrated the presence of human yolk sac-derived lymphoid-biased progenitors (YSLPs) expressing CD34, IL7R, LTB, and IRF8 at Carnegie stage 10, much earlier than the first HSC emergence. The number and lymphopoietic potential of these progenitors were both significantly higher in the yolk sac than the embryo proper at this early stage. Importantly, single-cell/bulk culture and CITE-seq have elucidated the tendency of YSLP to differentiate into innate lymphoid cells and dendritic cells. Notably, lymphoid progenitors in fetal liver before and after HSC seeding displayed distinct transcriptional features, with the former closely resembling those of YSLPs. Overall, our data identified the origin, potential, and migratory dynamics of innate lymphoid-biased multipotent progenitors in human yolk sac before HSC emergence, providing insights for understanding the stepwise establishment of innate immune system in humans.

Nupr1 Negatively Regulates Endothelial to Hematopoietic Transition in the Aorta-Gonad-Mesonephros Region.

In the aorta of mid-gestational mouse embryos, a specialized endothelial subpopulation termed hemogenic endothelial cells (HECs) develops into hematopoietic stem and progenitor cells (HSPCs), through a conserved process of endothelial-to-hematopoietic transition (EHT). EHT is tightly controlled by multiple intrinsic and extrinsic mechanisms. Nevertheless, the molecular regulators restraining this process remain poorly understood. Here, it is uncovered that, one of the previously identified HEC signature genes, Nupr1, negatively regulates the EHT process. Nupr1 deletion in endothelial cells results in increased HSPC generation in the aorta-gonad-mesonephros region. Furthermore, single-cell transcriptomics combined with serial functional assays reveals that loss of Nupr1 promotes the EHT process by promoting the specification of hematopoiesis-primed functional HECs and strengthening their subsequent hematopoietic differentiation potential toward HSPCs. This study further finds that the proinflammatory cytokine, tumor necrosis factor α (TNF-α), is significantly upregulated in Nupr1-deficient HECs, and the use of a specific TNF-α neutralizing antibody partially reduces excessive HSPC generation in the explant cultures from Nupr1-deficient embryos. This study identifies a novel negative regulator of EHT and the findings indicate that Nupr1 is a new potential target for future hematopoietic stem cell regeneration research.

Divergent expression of Neurl3 from hemogenic endothelial cells to hematopoietic stem progenitor cells during development.

Prior to the generation of hematopoietic stem cells (HSCs) from the hemogenic endothelial cells (HECs) mainly in the dorsal aorta in midgestational mouse embryos, multiple hematopoietic progenitors including erythro-myeloid progenitors and lymphoid progenitors are generated from yolk sac HECs. These HSC-independent hematopoietic progenitors have recently been identified as major contributors to functional blood cell production until birth. However, little is known about yolk sac HECs. Here, combining integrative analyses of multiple single-cell RNA-sequencing datasets and functional assays, we reveal that Neurl3-EGFP, in addition to marking the continuum throughout the ontogeny of HSCs from HECs, can also serve as a single enrichment marker for yolk sac HECs. Moreover, while yolk sac HECs have much weaker arterial characteristics than either arterial endothelial cells in the yolk sac or HECs within the embryo proper, the lymphoid potential of yolk sac HECs is largely confined to the arterial-biased subpopulation featured by the Unc5b expression. Interestingly, the B lymphoid potential of hematopoietic progenitors, but not for myeloid potentials, is exclusively detected in Neurl3-negative subpopulations in midgestational embryos. Taken together, these findings enhance our understanding of blood birth from yolk sac HECs and provide theoretical basis and candidate reporters for monitoring step-wise hematopoietic differentiation.

Single-cell dissection of cellular and molecular features underlying human cervical squamous cell carcinoma initiation and progression.

Cervical squamous cell carcinoma (CESC) is a prototypical human cancer with well-characterized pathological stages of initiation and progression. However, high-resolution knowledge of the transcriptional programs underlying each stage of CESC is lacking, and important questions remain. We performed single-cell RNA sequencing of 76,911 individual cells from 13 samples of human cervical tissues at various stages of malignancy, illuminating the transcriptional tumorigenic trajectory of cervical epithelial cells and revealing key factors involved in CESC initiation and progression. In addition, we found significant correlations between the abundance of specific myeloid, lymphoid, and endothelial cell populations and the progression of CESC, which were also associated with patients' prognosis. Last, we demonstrated the tumor-promoting function of matrix cancer-associated fibroblasts via the NRG1-ERBB3 pathway in CESC. This study provides a valuable resource and deeper insights into CESC initiation and progression, which is helpful in refining CESC diagnosis and for the design of optimal treatment strategies.

Cloning and characterization of SCIRR69: a novel transcriptional factor belonging to the CREB/ATF family.

The complete cDNA sequence of a novel gene, SCIRR69 (spinal cord injury and regeneration related no. 69 gene), was obtained by RACE technique. It codes for a protein of 521 amino acid residues homologous to human CREB3l2 (also known as BBF2H7) and mouse CREB3l2. The protein contains a basic DNA binding and leucine zipper dimerization (B-ZIP) motif and a hydrophobic region representing a putative transmembrane domain, similar to the structure of other CREB/ATF transcription factors. Monoclonal antibody against SCIRR69 was developed and could recognize the SCIRR69 protein in both native and denatured forms. Constructing of SCIRR69 fusion proteins with the GAL4 DNA-binding domain disclosed that SCIRR69 functioned as a transcriptional activator and its N-terminal 60 amino acids accounted for the activation ability. SCIRR69 resides in the cytoplasm of primary neurons, whereas neuron damage by incision led to the cleavage and translocation from the cytoplasm to the nucleus. These results suggest that SCIRR69 is activated by proteolytic cleavage at the transmembrane domain in response to neuron damage and its amino-terminal cytoplasmic domain translocates into the nucleus to activate the transcription of target genes.

Decoding lymphomyeloid divergence and immune hyporesponsiveness in G-CSF-primed human bone marrow by single-cell RNA-seq.

Granulocyte colony-stimulating factor (G-CSF) has been widely used to mobilize bone marrow hematopoietic stem/progenitor cells for transplantation in the treatment of hematological malignancies for decades. Additionally, G-CSF is also accepted as an essential mediator in immune regulation, leading to reduced graft-versus-host disease following transplantation. Despite the important clinical roles of G-CSF, a comprehensive, unbiased, and high-resolution survey into the cellular and molecular ecosystem of the human G-CSF-primed bone marrow (G-BM) is lacking so far. Here, we employed single-cell RNA sequencing to profile hematopoietic cells in human bone marrow from two healthy donors before and after 5-day G-CSF administration. Through unbiased bioinformatics analysis, our data systematically showed the alterations in the transcriptional landscape of hematopoietic cells in G-BM, and revealed that G-CSF-induced myeloid-biased differentiation initiated from the stage of lymphoid-primed multipotent progenitors. We also illustrated the cellular and molecular basis of hyporesponsiveness of T cells and natural killer (NK) cells caused by G-CSF stimulation, including the potential direct mechanisms and indirect regulations mediated by ligand-receptor interactions. Taken together, our data extend the understanding of lymphomyeloid divergence and potential mechanisms involved in hyporesponsiveness of T and NK cells in human G-BM, which might provide basis for optimization of stem cell transplantation in hematological malignancy treatment.

Pre-configuring chromatin architecture with histone modifications guides hematopoietic stem cell formation in mouse embryos.

The gene activity underlying cell differentiation is regulated by a diverse set of transcription factors (TFs), histone modifications, chromatin structures and more. Although definitive hematopoietic stem cells (HSCs) are known to emerge via endothelial-to-hematopoietic transition (EHT), how the multi-layered epigenome is sequentially unfolded in a small portion of endothelial cells (ECs) transitioning into the hematopoietic fate remains elusive. With optimized low-input itChIP-seq and Hi-C assays, we performed multi-omics dissection of the HSC ontogeny trajectory across early arterial ECs (eAECs), hemogenic endothelial cells (HECs), pre-HSCs and long-term HSCs (LT-HSCs) in mouse embryos. Interestingly, HSC regulatory regions are already pre-configurated with active histone modifications as early as eAECs, preceding chromatin looping dynamics within topologically associating domains. Chromatin looping structures between enhancers and promoters only become gradually strengthened over time. Notably, RUNX1, a master TF for hematopoiesis, enriched at half of these loops is observed early from eAECs through pre-HSCs but its enrichment further increases in HSCs. RUNX1 and co-TFs together constitute a central, progressively intensified enhancer-promoter interactions. Thus, our study provides a framework to decipher how temporal epigenomic configurations fulfill cell lineage specification during development.

Heterogeneity in endothelial cells and widespread venous arterialization during early vascular development in mammals.

Arteriogenesis rather than unspecialized capillary expansion is critical for restoring effective circulation to compromised tissues in patients. Deciphering the origin and specification of arterial endothelial cells during embryonic development will shed light on the understanding of adult arteriogenesis. However, during early embryonic angiogenesis, the process of endothelial diversification and molecular events underlying arteriovenous fate settling remain largely unresolved in mammals. Here, we constructed the single-cell transcriptomic landscape of vascular endothelial cells (VECs) during the time window for the occurrence of key vasculogenic and angiogenic events in both mouse and human embryos. We uncovered two distinct arterial VEC types, the major artery VECs and arterial plexus VECs, and unexpectedly divergent arteriovenous characteristics among VECs that are located in morphologically undistinguishable vascular plexus intra-embryonically. Using computational prediction and further lineage tracing of venous-featured VECs with a newly developed Nr2f2CrexER mouse model and a dual recombinase-mediated intersectional genetic approach, we revealed early and widespread arterialization from the capillaries with considerable venous characteristics. Altogether, our findings provide unprecedented and comprehensive details of endothelial heterogeneity and lineage relationships at early angiogenesis stages, and establish a new model regarding the arteriogenesis behaviors of early intra-embryonic vasculatures.