Abstract
Atrial cardiomyopathy is characterized by altered atrial structures and the genetic basis underlying the disorders remains inadequately explored. TIE1 variants or loss of function mutations were reported in a subset of lymphedema patients, and it is unknown whether the patients have also cardiac defects in addition to the lymphatic abnormality. We show in this study that endothelial Tie1 and Tek are highly expressed in endocardial cells of atria by the single cell RNA-seq analysis. TIE1 deficiency led to the disruption of atrial morphogenesis with minor defects in the ventricles. The bulk RNA-seq analysis of hearts at the four-chambered stage revealed that gene transcripts related to endothelial cell development and cardiac trabeculation were reduced in the Tie1 mutant mice compared with littermate controls. This was further confirmed by the RNA-seq analysis of atria and ventricles separately, showing more trabecular genes downregulated in the atria including Tek, upon the loss of Tie1. Consistent with the scRNA-seq data, we found that Tie1 and Tek transcripts were higher in atria than in ventricles. Furthermore, the endothelial deletion of Tek resulted in the defective formation of cardiac trabeculae, particularly in atria. Consistently, the loss of Tie1 combined with one null allele of Tek disrupted both atrial and ventricular trabeculation. Surprisingly, defects with the atrial chamber morphogenesis were already detectable 48 hours later upon the induced endothelial loss of Tie1 plus the Tek heterozygous deletion, implying a critical role of endocardium in the organization of the atrial internal muscular network. The synergy of TIE1 and TIE2 in the remodeling of atrial trabeculae was further confirmed at the postnatal stage, while TIE1 insufficiency alone had no obvious effect. Together, findings from this study imply that TIE1 is differentially required for the atrial and ventricular development and acts in synergy with TIE2 to regulate the endocardial cell-coordinated atrial internal muscular network assembly.
Introduction
The endocardium, as a specialized type of endothelial cells (EC) lining the cardiac lumen, plays an important role in vertebrate heart morphogenesis 1,2. Cardiac trabeculation relies on the coordinated interaction of endocardium and myocardium 3,4. The endocardium senses mechanical forces and provides signals for the myocardial morphogenesis 5. When endocardial cells fail to form buds, the assembly and extension process of trabeculae are disrupted 2,6. Defective trabecular formation could lead to embryonic lethality or congenital heart disease 7–9. Cell lineage studies have also shown that the endocardium is an important source of the coronary vascular endothelial cells, and gives rise to most mesenchymal cells that constitute the cardiac cushions and subsequently contribute to the heart valves 10. Atrial internal muscular networks, mainly pectinate muscles running in an anterolateral direction, are found on the anterior surface of both right and left atrial walls and corresponding auricles. In contrast to the better understanding of ventricular trabecula morphogenesis, the dynamic process and regulatory networks of atrial trabeculation are inadequately elucidated. Likewise, the atrial cardiomyopathy is characterized by altered atrial structure as well as function, and pathological mechanisms underlying the disorders are largely unknown.
Cardiac morphogenesis is a complex process involving multiple pathways to orchestrate the reciprocal paracrine signaling between the endocardial and other cardiac cells. Specifically, myocardial VEGFA activates VEGFR2 signaling for coronary vascular growth and also participate in the regulation of cardiac trabeculation 11. Endocardial cell-derived NRG1 regulates myocardium development via its receptors ERBB2/4 in cardiomyocytes. Mice null for Nrg1 or its receptors display hypoplastic ventricular wall lacking normal trabeculation 12. The aberrant NOTCH1 signaling was found to disrupt the heart morphogenesis 13,14. NRG1 has been show to regulate VEGFA expression in the myocardium, and VEGFA-VEGFR2 signaling is required for the NOTCH1 pathway restriction in endocardium during trabeculation. NOTCH1 signaling was shown to promote extracellular matrix (ECM) degradation during the formation of endocardial projections while NRG1 promoted myocardial ECM synthesis for trabecular rearrangement during the process of cardiac trabeculation 15. It has recently been shown that endocardium and myocardium interacted directly through tunneling nanotube–like structures extending from cardiomyocytes to endocardial cells, contributing to NOTCH1 activation 16. Factors involved in the synthesis or degradation of cardiac jelly, a layer of extracellular matrix between endocardium and myocardium, also play essential roles in the process of chamber morphogenesis, including glycosaminoglycan proteins such as versican and ADAMTS family proteins 17,18. The ANGPT1-TIE2 pathway plays a crucial role in cardiovascular development 19–22. ANGPT1 is expressed by cardiomyocytes, particularly in atrial myocardium and ventricular trabeculae, and participates in the atrial morphogenesis by regulating the spatiotemporal degradation of cardiac jelly 19. The endocardial specific deletion of Tek resulted in the fewer but thicker ventricular trabeculae as well as impaired endocardial sprouting, potentially via the paracrine suppression of retinoic acid signaling and proliferation in trabecular cardiomyocytes 23.
Endothelial TIE1 shares high homology with TIE2 and could form TIE1/TIE2 heterodimers 24. TIE1 has been shown to participate in the regulation of blood vascular and lymphatic network formation and also aortic valve development 25–29. It has recently shown that TIE1 and TIE2, via the PI3K-AKT-mediated signaling pathway for the protein stability of COUP-TFII, regulate venous specification 28,30. Loss of the transcription factor COUP-TFII was reported to disrupt both the vein and atrial development 31–33. Tie1 mutations or variants have also been linked to lymphatic disorders 34,35. It remains to be characterized whether Tie1 mutations could affect heart structure and function, and how TIE1 participates in the process of cardiac wall morphogenesis. To investigate the role of TIE1 and TIE2 in the process of endocardial cell-mediated chamber formation, we employed genetic mouse models targeting Tie1, Tek (encoding TIE2) or both at different developmental stages. Based on the analysis of the published single cell RNA sequencing data (scRNA-seq) by Feng et al. 36 and the RNA-seq analysis of atria and ventricles from the Tie1 mutant and littermate control mice, we found in this study that Tie1 and Tek expression were higher in atrial than ventricular endocardial cells. TIE1 is differentially required for the atrial and ventricular trabeculation, acting in synergy with TIE2.
Materials and Methods
Mouse models
All animal procedures performed conform the guidelines from Directive 2010/63/EU of the European Parliament on the protection of animals used for scientific purposes; and the experiment protocols were approved by the Animal Care Committee of Soochow and Nanjing University Animal Center (MARC-AP#YH2/SUDA20250507A03). All the mice used in this study were housed in a SPF (specific pathogen free) animal facility with a 12/12 hours dark/light cycle and were free to food and water access. Normal mouse diet (Suzhou Shuangshi Experimental Animal Feed Technology, Co, Ltd) and cage bedding (Suzhou Baitai Laboratory Equipment, Co, Ltd) were used. Two genetically modified mouse models targeting Tie1 gene were employed in this study. One mouse line targets TIE1 intracellular kinase domain (ICD), with exon 15 and exon 16 floxed, Tie1ICDFlox/Flox37. The other line is a knockout first mouse model established from EUCOMM embryonic stem cells (EPD0735-3B07) targeting Tie1 gene (Tie1tm1a/tm1a), in which targeting cassette is recombined downstream of exon 7 (with exon 8 floxed) 28. Tek knockout mouse model was generated as previously reported 38 and was crossbred with Tie1ΔICD/ΔICD mouse model to obtain Tie1ICD and Tek double knockout mice (Tie1ΔICD/ΔICD;Tek+/-). To generate mice with endothelial cell–specific gene deletion mouse models, we employed the Cdh5-CreERT2 mouse line 39. In all the phenotype analysis, wildtype or heterozygous littermates were used as controls. For the genotyping of Tie1 and Tek mouse lines, the primers used were as previously described 28,37. The genetic background of Tie1tm1a/tm1a is on C57BL/6 N, and the other lines are on C57BL/6J or C57BL/6J/SV129. Mice were sacrificed by asphyxiation with rising concentration of carbon dioxide gas, followed by cervical dislocation and tissue collection.
Induced gene deletion
Induction of gene deletion was performed as previously described by the tamoxifen treatment 28,38,40. Briefly, pregnant mice were treated with tamoxifen (Sigma-Aldrich, T5648-5G) at E10.5-E11.5, E12.5-14.5 or E12.5-15.5 (1 mg/per mouse for 2-4 consecutive days by intraperitoneal injection) and embryos were collected for analysis at the specified stages. For the postnatal studies, new-born pups were treated by four daily intragastric injections of tamoxifen starting from postnatal day 1 (P1, 60 μg/per mice daily). Tissues were collected for analysis at postnatal day 21. The genotypes of the Tie1/Tek knockout and control mice are as follows: Tie1ICDFlox/-;Cdh5-CreERT2 named as Tie1ICDiECKO, TekFlox/-;Cdh5-CreERT2 named as TekiECKO, Tie1ICDFlox/-;Cdh5-CreERT2;Tek+/- named as Tie1ICDiECKO;Tek+/-.
Embryo preparation
To investigate the developmental process of cardiac chamber formation, female mice were mated in the late afternoon and vaginal plugs were checked in the morning of the following day. The embryonic stages are estimated considering midday of the day on which the vaginal plug is present as embryonic day 0.5 (E0.5). Tail tips from mutant and control embryos were used for genotyping, and heart tissues from the embryos were collected for the subsequent histology and RNA analysis. For the RNA preparation, heart tissues were harvested and placed in cold PBS. The residual blood from the hearts was cleaned and heart tissues were snap-frozen in liquid nitrogen. For the separate analysis of atria and ventricles, heart tissues were collected and dissected before snap-freezing.
Immunostaining
For the frozen tissue section staining, tissues were collected and fixed in 4% paraformaldehyde for 2 hours at 4°C, followed by the incubation in 20% sucrose overnight before being embedded in optimal cutting temperature (OCT) compound. Frozen sections (10 µm) were used for the immunostaining analysis of cardiac chamber structures and coronary vessels. For the visualization of atrial internal muscular network-trabeculae structures, the whole-mount immunostaining was performed. Heart tissues were fixed in 4% paraformaldehyde overnight at 4°C, followed by blocking with 3%(w/v) skim milk in PBS-TX (0.3% Triton X-100) and incubating with the primary antibodies. The antibodies used were: hamster anti-mouse PECAM1(MAB1398Z; Millipore); rat anti-mouse CD31 (553370; BD Pharmigen), rat anti-mouse Endomucin (14-5851; eBioscience), goat-anti-mouse DLL4 (AF1389; R&D) and goat anti-human TIE1 (AF619; R&D). Appropriate Alexa 488(Invitrogen), Cy3(Jackson ImmunoReaearch) conjugated secondary antibodies were used. All fluorescently labeled samples were mounted and analyzed with a confocal microscope (Olympus Flueview 3000) or Leica MZ16F fluorescent dissection microscope. For comparison, the parameters were kept consistent for all the confocal microscopic imaging in this study.
Bulk RNA sequencing analysis of atria and ventricles
For the bulk RNA-seq analysis of heart tissue, the whole hearts from mutant and control embryos were collected. Atria and ventricles were separated and tissues from two embryos with the same genotype (Tie1tm1a/tm1a and control) were pooled for the RNA preparation and further analysis. According to the manufacturer’s instructions, total RNA was extracted from the tissues using Trizol (Invitrogen). The procedures for the RNA-seq analysis were performed as previously reported 28. Briefly, sequencing was performed on DNBSEQ platform with PE150 (read length, BGI-Shenzhen, China), and the sequencing data were processed according to the predefined analysis pipeline. The differential expression analysis was conducted using the R package DESeq2 (v 1.42.1). Functional enrichment analysis was performed using clusterProfiler (v4.10.1) and org.Mm.eg.db (v 3.18.0). The hallmark gene sets used in the GSEA (Gene Set Enrichment Analysis), including the hallmarks for trabeculation and EC migration, are mainly defined according to the transcriptome atlas of murine endothelial cells, the molecular signatures database (mh.all.v2024.1.Mm.symbols) and the related GO terms 28,41,42. The raw and processed data of RNA-seq analysis will be deposited in GEO after the publication.
Single cell RNA-seq data analysis
The single-cell RNA sequencing data by Feng et al. 36 were processed and analyzed with the Seurat package (v5.2.1) following the same parameters and criteria as described to ensure consistency and comparability. Major cell types were annotated based on established marker genes: pan-cardiomyocyte (Ttn), atrial cardiomyocytes (Sln), ventricular cardiomyocyte (Myl2), pan-endothelial cells (Pecam1), endocardial cells (Npr3), coronary endothelial cells (Fabp4), epicardial cells (Wt1, Tbx18, and Aldh1a2), fibroblast-like cells (Postn), mural cells (Pdgfrb), red blood cells (Hba-a1), macrophages (Adgre1). Following the identification of endocardial cells, subsequent analysis was performed using the Seurat package and normalized by the “LogNormalize” method, highly variable genes were identified with the ‘FindVariableFeatures” function (selection.method = “vst”, nfeatures = 2000). The optimal number of principal components (dims = 1:22) was determined by the ElbowPlot method. Cells were then clustered using the “FindNeighbors” and “FindClusters” functions (dims= 1:16, resolution = 0.5). Finally, Unsupervised clustering of cells will be performed by “RunUMAP” and gene expression was visualized with “FeaturePlot”. Ligand–receptor interactions between cardiomyocytes and endocardial cells during the embryonic stage (E10.5-E18.5) were identified using the CellChat package (v1.6.1). Atrial and ventricular cell subsets were analyzed separately using the mouse ligand-receptor database (CellChatDB.mouse) and communication probabilities were computed (min.cells = 10) to identify significant interactions.
Quantitative analysis of cardiac wall and trabeculae
For the quantification of trabecular structures in atria and ventricles, the EMCN positive areas was measured. The PECAM1 positive area was quantified to indicate the coronary vessels within the ventricular wall. The cardiac wall thickness was also measured as described in the supplemental figure 2. For the quantitative analysis of dorsal coronary veins, images were taken under a fluorescence microscope and kept constant near the venous sinus for all of the samples. Horizontal lines were evenly laid on the images, and the diameters of veins crossed with these lines were measured and analyzed using Image Pro Plus (Media Cybernetics, Inc., Bethesda, MD). The ratio of heart to embryo length was used to represent heart size. The starting point for the heart measurement was the junction of aorta and pulmonary artery, and the junction of left and right ventricles at the apex as the end point. All the quantification was measured and analyzed by using Image Pro Plus. Furthermore, we also performed the measurement of cardiac trabecular complexity by the Fractal analysis. Fractal dimensions were derived from confocal images of hearts from Tie1 mutants and control mice. Quantification was carried out using ImageJ with the Fraclac plugin (including 6 grid orientations).
Statistical analysis
Data analyses were performed using GraphPad Prism (version 8.0). For the 2-group comparison, the unpaired t test was performed with Welch correction if data passed the D’Agostino-Pearson normality test, or the unpaired nonparametric Mann-Whitney U test was applied. Data are expressed as mean ± SD. All statistical tests were 2-sided.
Results
Higher TIE1 and TIE2 in atrial than ventricular endocardial cells
To dissect the transcriptional dynamics of key endothelial regulators during the murine heart development, we analyzed the published single-cell transcriptomic dataset spanning the murine embryonic development (E10.5-E18.5) 36. We found that Tie1, and Tek were highly expressed in endocardial cells (Pecam1+, Npr3+) while the coronary EC gene such as the coronary endothelial marker Fabp4 was lowly expressed in endocardial cells (Fig. 1A-B, Supplemental Fig. 1A-C). Interestingly, we found that the transcript levels of Tie1, and Tek were higher in the atrial endocardial cells than those of ventricles, especially at embryonic stages up to E14.5 when the structures of four-chambered hearts were established (Table 1). In contrast, it appears that the Angpt1 transcript level was higher in atrial than ventricular cardiomyocytes (Fig. 1A-B, Table 1).

Analysis of TIE1 expression in endocardial cells by the single cell RNA-seq data and its role in cardiac trabeculation during embryogenesis.
A-B. Analysis of the single-cell RNA-seq data of embryonic hearts from E10.5 to E18.5 by Feng et al. Nat Commun 2022 (DOI: 10.1038/s41467-022-35691-7). UMAP visualization of 10 distinct cell clusters in the murine hearts during embryogenesis (E10.5 to E18.5, A). Tie1 and Tek are highly expressed in the endocardial cell population (Npr3⁺) as shown in feature plots. Angpt1 are highly expressed in the atrial cardiomyocytes (B). Black arrows point to the coronary EC population (Fabp4+). Red arrows point to the endocardial cell population (Npr3⁺). Blue arrows point to the atrial cardiomyocytes. The quantitative results are shown in Table 1. C-J. Analysis of heart trabeculae of Tie1tm1a/tm1a and littermate controls by immunostaining for EMCN (endomucin; green), PECAM1 (platelet endothelial cell adhesion molecule 1; red) and DAPI (blue) at E12.5 (C-D) and E14.5 (G-H). Note that Tie1tm1a/tm1a mice displayed defective trabeculation in the atrium (white arrows) and sparse trabecular structures in the ventricles (white arrows) at E12.5, with more severe defects by E14.5, lacking atrial trabeculae (white arrows) and sparser ventricular trabeculae (white arrows). Asterisks in C and G point to the delayed formation of interventricular septum in Tie1tm1a/tm1a mice. Arrows point to the trabecula-associated endocardium. E-F. Quantification of the EMCN-positive area in the atrial and ventricular trabeculae at E12.5 (E, LA: Tie1tm1a/tm1a: 0.63±0.33, n=6; Control: 1.00±0.22, n=8; P=0.054. RA: Tie1tm1a/tm1a: 0.77±0.26, n=6; Control: 1.00±0.13, n=8; P=0.075. LV: Tie1tm1a/tm1a: 0.86±0.16, n=6; Control: 1.00±0.046, n=8; P=0.054. RV: Tie1tm1a/tm1a: 0.91±0.11, n=6; Control: 1.00±0.12, n=8; P=0.13.). Quantification of the ventricular wall thickness (F, LV: Tie1tm1a/tm1a: 24.87±3.70 µm, n=6; Control: 27.67±4.52 µm, n=8; P=0.41. RV: Tie1tm1a/tm1a: 26.01±8.42 µm, n=6; Control: 32.84±11.21 µm, n=8; P=0.28.). I-J. Quantification of the EMCN-positive area in the atrial and ventricular trabeculae at E14.5 (I, LA: Tie1tm1a/tm1a: 0.40±0.16, n=7; Control: 1.00±0.079, n=7; P=0.00060. RA: Tie1tm1a/tm1a: 0.36±0.17, n=7; Control: 1.00±0.092, n=7; P=0.00060. LV: Tie1tm1a/tm1a: 0.81±0.22, n=7; Control: 1.00±0.033, n=7; P=0.036. RV: Tie1tm1a/tm1a: 0.77±0.23, n=7; Control: 1.00±0.040, n=7; P=0.026.). Quantification of the PECAM1-positive area in ventricular walls (J, LV: Tie1tm1a/tm1a: 0.91±0.15, n=4; Control: 1.00±0.013, n=6; P=0.56. RV: Tie1tm1a/tm1a: 0.64±0.096, n=4; Control: 1.00±0.024, n=6; P=0.0095.). The quantification data of mutant mice was normalized against that of littermate control mice. In E-J, statistical analysis was performed using the unpaired nonparametric Mann-Whitney U test. Values are represented as means ± SD of at least three independent technical replicates. LA, Left atria, RA, Right atria, LV, Left ventricle, RV, Right ventricle, Scale bar: Scale bar: 400 μm in C and G, 50 μm in D and H.

The transcript levels of Tie1, Tek in atrial and ventricular endocardial cells and Angpt1 in atrial and ventricular cardiomyocytes during the embryonic stage from E10.5-E18.5 (based on the published scRNA-seq data, Feng et al., Nat Commun 2022; DOI: 10.1038/s41467-022-35691-7).
Values are represented as means ± SD and the number referred to cell counts.
To investigate mechanisms underlying the murine atrial trabeculation, we first analyzed the dynamic processes of cardiac chamber morphogenesis by the immunostaining analysis of endomucin (EMCN) for trabeculae and PECAM1 for coronary vessels (embryonic day 10.5-E18.5). The quantification methods for atrial and ventricular trabecular areas, wall thickness and vessel density were shown in Supplemental Fig. 2A. As shown in Supplemental Fig. 2B, the heart tube became looped to form the preliminary chamber structures at E10.5, including the initiation of interventricular septum. While there appeared to have some wavy edges or the budding of endocardium in atria (E10.5), trabecular structures were readily detectable in ventricles at this stage. Therefore, atrial wall morphogenesis occurs later than that of ventricles, which initiates the trabeculation process around E8.5 when endocardial cells migrate through cardiac jelly to form anchor points with compact myocardium 15. We found that pectinate muscles, defined as the atrial trabeculae, became detectable in the atria at E12.5 (Supplemental Fig. 2B). The formation of interventricular septum was completed by E14.5, establishing the four-chambered heart. By E16.5 or later stages (Supplemental Fig. 2C), the atrial and ventricular trabeculae gradually formed a complex network. As shown in Supplemental Fig. 2D-G (Supplemental Table 1), the trabecular areas in both the atria and ventricles gradually increased, being greater in the right atria (RA) than that in the left atria (LA, Supplemental Fig. 2D). There was no obvious difference in the trabecular areas between the left and right ventricles (Supplemental Fig. 2E). The ventricular wall thickness increased rapidly after E12.5, with the left ventricular wall being thicker than that of the right ventricular wall by E18.5 (Supplemental Fig. 2F), accompanied by the increase of coronary vessel density (Supplemental Fig. 2G).
Differential requirement of TIE1 in atrial and ventricular trabeculation
To study the role and mechanism underlying TIE1 function in cardiac chamber morphogenesis, we employed two genetically modified mouse models targeting Tie1, including Tie1tm1a/tm1a and Tie1ΔICD/ΔICD. The effects of Tie1 gene deletion on atrial and ventricular trabeculae were analyzed at different developmental stages (including E12.5, E14.5 and E17.5). TIE1 deficiency led to an impaired trabeculation in the atria, characterized by a decrease of trabecular structures at E12.5, while the formation of ventricular trabeculae was not obviously affected at the same stage (Fig. 1C-D). Shown in Fig. 1E was the quantification of trabecular areas in atria and ventricles. The thickness of left and right ventricular walls at this stage was shown in Fig. 1F. Consistently, the atrial trabeculae were almost absent in the Tie1 mutant mice at E14.5, while trabeculae became sparse in both left and right ventricles (Fig. 1G-I). There was also a trend of decrease in the coronary vessel density upon TIE1 deficiency (Fig. 1J). By the later stages of embryonic development (E17.5), TIE1 deficiency resulted in more severe cardiac defects, including a significant reduction in heart size, absence of atrial trabeculae, and disorganized trabecular structures in ventricles, together with a decrease in the ventricular wall thickness as well as the coronary vessel density (Fig. 2A-F). The defective atrial trabecula formation was also validated in another Tie1 mutant line targeting its intracellular kinase domain (Tie1ΔICD/ΔICD, Supplemental Fig. 3A-E). Consistent with the above results, TIE1 deficiency resulted in the lack of atrial trabeculation while trabeculae in ventricles became sparse in Tie1 mutants at E14.5. There was also a trend of decrease in ventricular wall thickness as well as the coronary vessel density in the Tie1ΔICD/ΔICD mice. Consistently, the atrial trabecular defects could also be reflected by the fractal analysis of cardiac trabecular complexity in atria from Tie1 mutant and control mice at E14.5 and E17.5 using the Fraclac Plugin in ImageJ, while the relatively minor defects in the ventricles were not detected by the method (Tie1tm1a/tm1a, Supplemental Fig. 4; Supplemental Table 2).

Suppression of cardiac trabeculation-related gene expression by the loss of TIE1.
A-B. Analysis of heart trabeculae of Tie1tm1a/tm1a and littermate controls by immunostaining for EMCN (green), PECAM1 (red) and DAPI (blue) at E17.5. Note that the ventricular trabeculae of Tie1tm1a/tm1a mice become sparse compared with those of control mice (white arrows). Arrows point to the trabecula-associated endocardium. C-F. Quantification of the EMCN-positive area in the atrial and ventricular trabeculae (C, LA: Tie1tm1a/tm1a: 0.56±0.16, n=6; Control: 1.00±0.22, n=8; P=0.0043. RA: Tie1tm1a/tm1a: 0.57±0.33, n=6; Control: 1.00±0.24, n=8; P=0.026. D, LV: Tie1tm1a/tm1a: 0.86±0.18, n=6; Control: 1.00±0.082, n=8; P=0.054. RV: Tie1tm1a/tm1a: 0.61±0.18, n=6; Control: 1.00±0.10, n=8; P=0.00070). Quantification of the PECAM1-positive area in ventricular walls (E, LV: Tie1tm1a/tm1a: 0.66±0.14, n=6; Control: 1.00±0.11, n=8; P=0.00070. RV: Tie1tm1a/tm1a: 0.55±0.10, n=6; Control: 1.00±0.088, n=8; P=0.00070.). Quantification of the ventricular wall thickness (F, LV: Tie1tm1a/tm1a: 101.97 ± 22.44 µm, n=6; Control: 174.09 ± 22.84 µm, n=8; P=0.0013. RV: Tie1tm1a/tm1a: 87.14±26.20 µm, n=6; Control: 137.35±24.06 µm, n=8; P=0.0080.) at E17.5. Except the ventricular wall thickness, the quantification data of mutant mice was normalized against that of littermate control mice. G-L. RNA-seq analysis of whole hearts of Tie1tm1a/tm1a and littermate controls (E14.5, n=4 per group). G. PCA analysis showed that Tie1tm1a/tm1a and control groups were divided into two clusters. H. Volcano plot visualization of the differentially expressed genes (DEGs, P-value ≤ 0.05, |Log2Foldchange| > 0). I. Heatmap of significantly downregulated cardiac trabeculation genes (P-value ≤ 0.05). J. GSEA analysis showed the decrease of enrichment in trabeculation, venous development, and endothelial cell migration pathways. K-L. GO analysis revealed that gene related to cardiac chamber morphogenesis and trabeculae development, extracellular matrix, cell adhesion and endothelial migration/differentiation were significantly downregulated (K), while cardiac contraction and apoptosis related genes were significantly upregulated (L). In C-F, statistical analysis was performed using the unpaired nonparametric Mann-Whitney U test. Values are represented as means ± SD of at least three independent technical replicates. LA: Left atria, RA: Right atria, LV: Left ventricle, RV: Right ventricle; Scale bar: 400 μm in A, 50 μm in B; c-c: cell-cell, p-m: plasma membrane in K.
Consistent with the previous observation 28, the lack of TIE1 also disrupted the coronary vein formation at the later embryonic stages (E14.5 and E17.5), including the decrease of vein diameter and the increase of vein-associated angiogenic sprouting (Supplemental Fig. 5A-D).
Alteration of atrial versus ventricular gene regulatory network upon TIE1 deficiency
To explore the mechanism underlying TIE1 in heart development, we performed the RNA sequencing (RNA-seq) analysis of hearts to examine the alteration of cardiac transcriptome upon TIE1 deficiency. Lack of TIE1 led to a dramatic change in the overall transcriptional profile (Fig. 2G-L), including the suppression of gene expression related to the endocardial / endothelial development and cardiac chamber morphogenesis, such as Tek (normalized by Pecam1, Tie1tm1a/tm1a: 0.89 ± 0.031, n=4; WT: 1.00 ± 0.12, n=4; P = 0.057.). Specifically, pathways related to extracellular matrix, cell adhesion and endothelial cell migration were downregulated (Fig. 2K), while the cardiac contraction and cell apoptosis pathways were upregulated (Fig. 2L). This was further confirmed by the GSEA analysis using the custom gene sets (Fig. 2J). The gene sets included the following entries based on the GO database (HALLMARK_TRABECULATION, HALLMARK_EC_MIGRATION) and the fifty-three entries from the previous publication 28. Consistent with previous reports 28, TIE1 deficiency led to the downregulation of venous genes. In addition, key regulatory factors for cardiac trabecular development were downregulated, including Adgrg6 (adhesion G protein-coupled receptor G6), Nrg1 (neuregulin 1), and Notch1 (Notch receptor 1) and Dll1 (Delta Like Canonical Notch Ligand 1). Notably, we found that the expression of the TIE receptor family member Tek was also significantly reduced (Fig. 2I), suggesting that TIE1 acts, at least partially, via TIE2 in the regulation of cardiac trabecular morphogenesis.
To better visualize the trabecular structures, particularly the atrial internal muscular network, we performed the wholemount immunostaining of the hearts at the four-chambered stage. As shown in Fig. 3A, TIE1 deficiency led to the absence of the atrial trabecular network, with minor defects in the ventricles including a trend of sparse trabeculae (the 3D views shown in Supplemental Fig. 6). To further investigate the differential requirement of TIE1 in atrial and ventricular trabeculation, we collected atria and ventricle tissues separately from Tie1 mutant and control mice (E14.5) for the RNA-seq analysis. We found that the number of differentially expressed genes was greater in atria than in the ventricles after the loss of TIE1 (Fig. 3B-G). Similar to the results from the whole heart RNA-seq analysis, the transcript levels of genes related to atrial development was decreased. Although a similar trend of decrease in genes involved in the cardiac morphogenesis were also observed in the ventricles, the enrichment levels and the number of genes involved were much lower in ventricles than in atria. Specifically, 17 trabecula-related genes were downregulated in the atria, while there was only 8 in the ventricles. Notably, Tek and Dll1 were downregulated in both the atria and ventricles while Notch1 only in atria of Tie1 mutant mice (Fig. 3F-G; Supplemental Fig. 3F-G; Supplemental Table 3 and Supplemental Table 4). Interestingly, the expression of endothelial Tie1, Tek and Notch1 was higher in atria than in ventricles (Fig. 3H, Supplemental Fig. 1D, Supplemental Table 5). The top GO terms enriched for the up- or down-regulated genes were related to the biological processes of heart development, including those in trabecular formation, extracellular matrix organization and binding, intercellular adhesion and endothelial cell migration in both atria and ventricles (Fig. 3I-J). In addition, the analysis of scRNA-seq data for the ligand-receptor interactions showed that the atrial_EC (endocardial cells) and atrial_CM (cardiomyocytes) exhibited stronger ANGPT1-TIE2 signaling compared with those in ventricles (Supplementary Fig. 1E). This may account for the differential requirement of TIE1 and also TIE2 in the process of atrial versus ventricular trabeculation.

Differential effect of Tie1 deficiency on atrial versus ventricular trabeculation related gene network.
A. Visualization of heart trabeculae structure of Tie1tm1a/tm1a and littermate controls by whole-mount immunostaining for EMCN (green), PECAM1 (red) and DAPI (blue) at E14.5. Note that Tie1 deletion resulted in the loss of atrial trabecular network with minor defects with ventricular trabecular structures (white arrows). Arrows point to the trabecula-associated endocardium. B-J. RNA-seq analysis of atrium and ventricles of Tie1tm1a/tm1a and littermate controls (E14.5, n=3 per group). B. PCA analysis results showed clear separation between the atria and ventricles of Tie1tm1a/tm1a and control mice. C. Differential expression analysis (P-value ≤ 0.05, |log2(Foldchange)| ≥ 0.15) identified over 693 down-regulated and 948 up-regulated atrial genes, while 144 down-regulated and 190 up-regulated genes in ventricles. D-E. Volcano plots of atrial (D) and ventricular (E) differential expressed genes (DEGs). F-G. Heatmap showed significantly downregulated (P-value ≤ 0.05) heart trabeculation genes in atrium (F) and ventricles (G), including Tek, Dll1 and Notch1. H. Quantitation of Tie1 and Tek transcript levels (normalized by Pecam1) in atrium and ventricles of Tie1tm1a/tm1a and littermate controls (Tie1: Ventricles: 0.72 ± 0.021, n = 5; Atria: 1.00 ± 0.026; P = 0.0079. Tek: Ventricle: 0.55 ± 0.024; Atria: 1.00 ± 0.048; P = 0.0079). I-J. GO enrichment analysis of the DEGs in atrium and ventricles including the heart trabecula morphogenesis, extracellular matrix, cell adhesion and endothelial migration related process down-regulated while genes related to the cardiac contraction, apoptosis process up-regulated. There was a more obvious enrichment detected in the atria. In H, statistical analysis was performed using the unpaired nonparametric Mann-Whitney U test. Values are represented as means ± SD of at least three technical replicates. LA: left atria, RA: right atria, LV: left ventricle, RV: right ventricle; Scale bar: 50 μm in A; c-c: cell-cell, p-m: plasma membrane in I-J.
Stronger ANGPT1-TIE2 interactions among atrial endocardium and myocardium
Using the published scRNA-seq data 36, we also analyzed the ligand-receptor interactions among endocardial cells and cardiomyocytes in atria and ventricles respectively during the embryonic stage (E10.5-E18.5). We found that the ANGPT1-TIE2 signaling between atrial cardiomyocytes and endocardial cells was stronger than that in ventricles (Supplemental Fig. 1E). Consistently, the crucial role of endothelial TIE2 in the atrial trabeculation was validated in this study by employing a conditional knockout mouse model targeting Tek (also known as Tie2 gene; TekFlox/-; Cdh5-CreERT2; named TekiECKO). As the complete knockout of Tek gene leads to embryonic lethality before E10.5, we first performed the induced endothelial Tek deletion starting from E12.5 (Fig. 4A-B). We found that TIE2 insufficiency led to an obvious decrease of trabeculae in the atrium 48 h later (E14.5), while there were no obvious defects observed in the ventricular trabeculae at this stage (Fig. 4C-E). The coronary vessel density in the ventricular wall as well as the wall thickness, including left and right ventricles, were not obviously altered when examined 48 h later (Fig. 4F). In separate experiments with the induced gene deletion from E12.5 to E14.5, embryos were collected for the analysis at 72 hours (E15.5) or later (E17.5). The formation of atrial trabeculae was completely missing as shown in Fig. 4A, G-H. The defective trabecula formation was also detected in the ventricles, especially at later embryonic stages (E17.5, Fig. 4A, G-H). TIE2 insufficiency also led to the defective formation of the coronary veins as previously shown in skin and mesentery 30 (Supplemental Fig. 5E, G).

Impaired atrial trabecular formation after TIE2 insufficiency.
A. Tamoxifen intraperitoneal (i.p.) administration and analysis scheme. B-D. Analysis of heart trabeculae of TekiECKO and littermate controls by immunostaining for EMCN (green), PECAM1 (red) and DAPI (blue) 48 hours later (tamoxifen treatment starting at E12.5). E-F. Quantification of the EMCN-positive area in the atrial and ventricular trabeculae (E, LA: TekiECKO: 0.60±0.20, n=9; Control: 1.00±0.097, n=8; P=0.00010. RA: TekiECKO: 0.64±0.23, n=9; Control: 1.00±0.22, n=8; P=0.0047. LV: TekiECKO: 0.85±0.15, n=9; Control: 1.00±0.057, n=8; P=0.020. RV: TekiECKO: 0.90±0.16, n=9; Control: 1.00±0.089, n=8; P=0.13.), quantification of the PECAM1-positive area in ventricular walls (F, LV: TekiECKO: 0.92±0.25, n=7; Control: 1.00±0.048, n=6; P=0.81. RV: TekiECKO: 1.22±0.53, n=7; Control: 1.00±0.13, n=6; P=0.28.) and quantification of the ventricular wall thickness (F, LV: TekiECKO: 62.35±11.70 µm, n=9; Control: 62.59±10.26 µm, n=8; P=0.96. RV: TekiECKO: 63.87± 9.19 µm, n=9; Control: 58.87±11.02 µm, n=8; P=0.32.) at E14.5. Except the ventricular wall thickness, the quantification data of mutant mice was normalized against that of littermate control mice. G-H. Analysis of heart trabeculae of TekiECKO and littermate controls by immunostaining for EMCN (green), PECAM1 (red) and DAPI (blue) 72 hours after the induced gene deletion or later (tamoxifen treatment starting at E12.5). Arrows point to the trabecula-associated endocardium. In E-F, statistical analysis was performed using the unpaired t test. Values are represented as means ± SD of at least three independent technical replicates. LA: left atria, RA: right atria, LV: left ventricle, RV: right ventricle; Scale bar: 500 μm (embryo) and 400 μm (heart) in B, 400 μm in C, 50 μm in D, G, H.
Synergy of TIE1 and TIE2 in atrial trabeculation
To investigate the synergistic effects of TIE1 and TIE2 in the regulation of cardiac morphogenesis, we generated a doubly knockout model targeting Tie1 and Tek (Tie1ΔICD/ΔICD; Tek+/-). There were no obvious defects observed with the atrial or ventricular structures upon Tie1 deficiency alone at E10.5 (Fig. 5A-B). Intriguingly, loss of Tie1 combined with the heterozygous deletion of Tek suppressed both atria and ventricle development (E10.5), including reduced heart size, abnormal atrial chamber formation, and fewer trabecular structures in the ventricles (Fig. 5C-D), leading to embryonic lethality before E12.5.

Developmental requirement of TIE1 and TIE2 in atrial and ventricular chamber morphogenesis.
A-D. Analysis of heart trabeculae of Tie1ΔICD/ΔICD, Tie1ΔICD/ΔICD;Tek+/- and littermate controls by immunostaining for EMCN (green), PECAM1 (red) and DAPI (blue) at E10.5. Note that there were no obvious defects observed with the atrial or ventricular structures upon Tie1 deficiency alone at E10.5, while loss of Tie1 combined with the heterozygous deletion of Tek suppressed both atria and ventricle development. White arrows point to the trabecular structures. Quantification of the EMCN-positive area in the ventricular trabeculae at E10.5, and the quantification data of mutant mice was normalized against that of littermate control mice (B, LV: Tie1ΔICD/ΔICD: 0.86±0.11, n=4; Control: 1.00±0.14, n=5; P=0.29. RV: Tie1ΔICD/ΔICD: 0.79±0.19, n=4; Control: 1.00±0.15, n=5; P=0.19. D, Tie1ΔICD/ΔICD;Tek+/-: 0.64±0.22, n=6; Control: 1.00±0.16, n=7; P=0.0047). E. For the induced gene deletion, tamoxifen intraperitoneal (i.p.) administration (starting from E10.5) and analysis scheme (at E14.5-E16.5). F. The deletion efficiency of TIE1 was examined by the immunostaining of skins from Tie1ICDiECKO;Tek+/- and littermate control mice at E16.5 for PECAM1 (green) and TIE1 (red). G-H. Analysis of heart trabeculae of Tie1ICDiECKO;Tek+/- mice and littermate controls by immunostaining for EMCN (green), PECAM1 (red) and DAPI (blue) at E14.5. Note that the atrial trabeculae structure was missing, while the ventricular trabeculae became sparse (white arrows). Arrows point to the trabecula-associated endocardium. In B and D, statistical analysis was performed using the unpaired nonparametric Mann-Whitney U test. Values are represented as means ± SD of at least three independent technical replicates. LA: left atria, RA: right atria, LV: left ventricle, RV: right ventricle. Scale bar: 200 μm in A(left) and C(left), 30 μm in A(right) and C (right), 400 μm in G, 200 μm in F and 50 μm in H.
To further explore the synergy of TIE1/TIE2 in the heart development at later stages of embryonic development, we generated an inducible mouse model with the endothelial deletion of Tie1 combined with Tek heterozygous knockout (Tie1ICDFlox/- ;Cdh5-CreERT2;Tek+/-, named Tie1ICDiECKO;Tek+/-). The induced gene deletion was performed by the tamoxifen treatment starting from E10.5, and mice were sacrificed for analysis after 96 hours (E14.5). The gene deletion efficiency was confirmed by the immunostaining for PECAM1 plus TIE1. TIE1 insufficiency combined with the heterozygous loss of TIE2 led to the suppression of atrial trabeculation, while the trabeculae at both right and left ventricle became sparse and irregular (Fig. 5E-H). Furthermore, when the induced gene deletion was performed at E12.5, the defective atrial trabeculation were already detectable in the mutant mice 48 hours later (E14.5 Tie1ICDiECKO;Tek+/-; Fig. 6 A-D). Consistent with the observation with the endothelial Tek deletion for 48 hours, there was no obvious defects with the ventricular trabeculation, coronary vessel density and ventricular wall thickness (Tie1ICDiECKO;Tek+/-; Fig. 6D-E, Table 2). The disruption of cardiac trabeculation became more obvious when analyzed 72-96 hours later, particularly in atria but less severe with ventricles (Tie1ICDiECKO; Tek+/-, E15.5-E16.5; Fig. 6 F-G, Table 3). There was also a trend of decrease in coronary vessel density and ventricular wall thickness as shown at later stages (E16.5, Fig. 6 F, H, Table 3). In addition, the TIE1/TIE2 insufficiency also resulted in the disruption of cardiac vein development as shown in (Supplemental Fig. 5F, H), as previously reported in skin and retina 28.

Synergistic role of TIE1 and TIE2 in the atrial trabeculation.
A. Tamoxifen intraperitoneal (i.p.) administration and analysis scheme. B-H. The induced gene deletion by tamoxifen started at E12.5. Analysis of heart trabeculae of Tie1ICDiECKO;Tek+/-and littermate controls by immunostaining for EMCN (green), PECAM1 (red) and DAPI (blue) at E14.5 (B-C) and E16.5 (F). Note that the atrial trabeculae were absent and the ventricular trabeculae were slightly sparse in the doubly mutant mice at E14.5 (B-C), and the phenotype became more severe at E16.5 (F). Quantification of the EMCN-positive area in the atrial and ventricular trabeculae (D, E14.5; G, E16.5), the PECAM1-positive area in ventricular walls and the ventricular wall thickness (E, E14.5; H, E16.5). Arrows point to the trabecula-associated endocardium. The quantification data of mutant mice was normalized against that of littermate control mice as shown in Table 2 and Table 3. Arrows in C and F point to the trabecula-associated endocardium. In D, E, G and H, statistical analysis was performed using the unpaired nonparametric Mann-Whitney U test. Values are represented as means ± SD of at least three independent technical replicates. LA: left atria, RA: right atria, LV: left ventricle, RV: right ventricle. Scale bar: 400 μm in B, 50 μm in C and F.


Quantitation of the trabecular area (EMCN+), blood vessel density in ventricular wall (PECAM1+) and ventricular wall thickness in Tie1ICDiECKO;Tek+/- and littermate controls (at E14.5, 48 hours after the induced gene deletion).
Statistical analysis was performed using the unpaired nonparametric Mann-Whitney U test. Values are represented as means ± SD of at least three independent technical replicates.

Quantitation of the trabecular area (EMCN+), coronary vessel density in ventricular wall (PECAM1+) and ventricular wall thickness in Tie1ICDiECKO;Tek+/- and littermate controls 96 hours after the induced gene deletion (E16.5; tamoxifen treatment from E12.5-E13.5).
Statistical analysis was performed using the unpaired nonparametric Mann-Whitney U test. Values are represented as means ± SD of at least three independent technical replicates.
Requirement of TIE receptors for the postnatal atrial trabecular remodeling
The role of TIE1 and its synergy with TIE2 in cardiac trabecular development was further confirmed at the neonatal stage (Tie1ICDiECKO or Tie1ICDiECKO; Tek+/-) with the induced gene deletion from postnatal day 1 (P1-4) and the analysis performed at P21. The Tie1 deletion was confirmed by the vascular defects in retinas of Tie1ICDiECKO or Tie1ICDiECKO; Tek+/-mice as previously reported (Fig. 7B, E) 28. TIE1 insufficiency alone had no obvious impact on the heart weight and also the atrial or ventricular trabeculation analyzed within three weeks (Fig. 7 A-D). Consistent with the studies at embryonic stages, the induced postnatal gene deletion (Tie1ICDiECKO;Tek+/-) showed the retarded heart development by the quantification of heart-body weight (Fig. 7 F, G), and showed abnormal atrial trabeculation while the trabeculae in the ventricles had no obvious defects (Fig. 7F, H). There was a significant decrease of trabecular structures (white arrows in Fig. 7H) in atria of Tie1ICDiECKO; Tek+/- mice compare with that of the control mice (LA, Tie1ICDiECKO; Tek+/-: 0.69±0.16, n=11; Control: 1.00±0.28, n=16; P=0.0022. RA, Tie1ICDiECKO; Tek+/-: 0.77±0.15, n=11; Control: 1.00±0.11, n=16; P=0.0002). The effects of TIE1 and TIE2 deficiency or insufficiency on heart trabecular morphogenesis were summarized in Fig. 7I. Based on the previous findings about the TIE1/TIE2 in venous EC specification 28,30, we propose that TIE1 acts in synergy with TIE2 in the specification of atrial endocardial cells (Fig. 7I). TIE1 deficiency or TIE1/TIE2 insufficiency results in a more severe defects with the atrial trabeculation, suggesting the differential requirement of TIE1 and TIE2 in the endocardium-mediated organization of the atrial internal muscular network.

Abnormal atrial trabecular remodeling after the postnatal deletion of Tie1 combined with one null allele of Tek.
A. Tamoxifen intragastric (i.g.) administration and the analysis scheme. B-H. Analysis of blood vessels in the retinas of Tie1ICDiECKO (B), Tie1ICDiECKO;Tek+/−(E) and control mice by whole-mount immunostaining for PECAM1 (green), DLL4 (red) at P21. Note that hemangioma-like vascular tufts observed in Tie1ICDiECKO and Tie1ICDiECKO;Tek+/−mice. The retinal vascular defects were used to confirm the efficiency of induced gene deletion. Analysis of heart trabeculae in Tie1ICDiECKO (C), Tie1ICDiECKO;Tek+/−(F) and control mice by the immunostaining for EMCN (green), PECAM1 (red) and DAPI (blue) at P21. TIE1 insufficiency alone had no obvious effect on heart size (D, Tie1ICDiECKO: 0.54±0.021, n=7; Control: 0.54±0.014, n=5; P=0.91.) and atrial trabecular structures (C). The Tie1ICDiECKO;Tek+/- mice displayed smaller hearts (G, Tie1ICDiECKO;Tek+/-: 0.49±0.086, n=8; Control: 0.59±0.10, n=11; P=0.035.) and abnormal atrial trabecular morphology (F,H). White arrows in H point to the trabecula-associated endocardium in the atria. I. Schematic illustration of TIE1, TIE2 and their synergy in cardiac chamber morphogenesis. TIE1 deficiency, TIE2 or TIE1/TIE2 insufficiency results in a more severe defect with the atrial trabeculation than that of ventricles. Based on the previous findings about the TIE1/TIE2 in venous EC specification 28,30 and findings from this study, we propose that TIE1 acts in synergy with TIE2 in the specification and remodeling of atrial endocardium, and that there is a differential requirement of TIE receptors in atrial versus ventricular wall morphogenesis. In D and G, statistical analysis was performed using the unpaired nonparametric Mann-Whitney U test (D) and the unpaired t test (G). Values are represented as means ± SD of at least three independent technical replicates. Scale bar: 50 μm in B and E, 1 mm in C and F, 200 μm in H.
Discussion
Developmental defects or pathological remodeling of atria may result in the atrial dysfunction leading to the heart failure 43. In contrast to the better understanding of ventricular chamber morphogenesis, the dynamic process and regulatory networks of the atrial internal muscular network formation (known as atrial trabeculation) are inadequately elucidated. TIE1 variants or loss of function mutations were reported in a subset of lymphedema patients 34,35, and it is unknown whether the patients have also cardiac defects in addition to the lymphatic abnormality. We show in this study that the atrial trabeculation occurred later than in ventricles. Based on the analysis of the published scRNA-seq data by Feng et al. 36 and the bulk RNA-seq analysis of atria and ventricles separately in this study, we showed that Tie1 and its homologue Tek were highly expressed in atrial endocardial cells than those of ventricles. By employing the genetically modified mouse models targeting Tie1, Tek and both, we found that TIE1 was differentially required in the process of atrial versus ventricular trabeculation and acted in synergy with TIE2 in the process of cardiac trabeculation, especially in atria.
Endocardial versus vascular effects of TIE1
Cardiac chamber morphogenesis involves the coordinated endocardial-myocardial interaction governed by a plethora of cardiovascular regulators 2. It is worth noting that most of the reported endothelial regulators have important roles in both heart and vascular development. ANGPT1 has been shown to participate in the regulation of atria and cardiac vein development 19,44. We have demonstrated that TIE1 and TIE2 are required for the vein specification via regulating the protein stability of COUPTFII 30. Disruption of COUPTFII led to the defective atrial and venous formation 31,33. Defective cardiac venous formation was also observed in this study in Tie1 or Tek mutant mice at later embryonic stages. It is therefore important to dissect whether cardiac defects in the mutants are a direct effect or secondary to the coronary vascular defects. By performing a serial analysis of the developing hearts during embryogenesis (including E10.5, E12,5, E14.5, E16.5, E17.5), we found in this study that TIE1 deficiency resulted in an impaired atrial trabeculation readily detectable at embryonic day 12.5 (E12.5), while there were only minor defects observed with the ventricular trabeculae at the same stage. In addition, Tek deficiency leads to the embryonic lethality before E10.5 30,45. Interestingly, the defective atrial trabeculae formation was detected within 48 hours after the induced deletion of endothelial Tek (starting from E12.5), while there were no obvious defects with coronary vessels at the same stage. This was further confirmed in the Tie1/Tek doubly mutant mice with defective atrial trabeculation detected within 48 hours after the induced deletion of Tie1 combined with one null allele of Tek. Lack of TIE1 did not produce obvious developmental defects in mutant embryos by E10.5 as previously reported 26,28. Therefore, the differential requirement of TIE1 may be related to the specification of atrial versus ventricular endocardium, which involves a synergistic role of TIE2, as previously reported in vein development 28,30. Based on the above findings, it is likely that TIE1 and TIE2 are directly required for the endocardial cell-mediated organization of atrial trabecular structures. In addition, defects with the ventricular trabeculation and a trend of thinner ventricle walls were also observed at later stages after the Tie1 deficiency or induced Tie1/Tek insufficiency. The cardiac trabeculation contributes to the thickening and maturation of the ventricular wall by being incorporated into the compact myocardium 4,8. It has been reported that endocardial cells contribute to the formation of coronary vascular network formation 10. The defective ventricular wall morphogenesis may also be related to the decreased coronary vessel growth in Tie1 or Tie1/Tek mutant mice as observed in this study.
Synergy of TIE1 and TIE2 in cardiac trabeculation
Loss of TIE1 alone did not produce an obvious defect with atrial chamber formation at initial stages of heart development (e.g. E10.5). However, the loss of Tie1 combined with the heterozygous deletion of Tek disrupted the atrial and ventricle development at the same stage (E10.5). The synergy of TIE1 and TIE2 in heart development was further verified at later embryogenesis as well as the postnatal studies. The defects with the atrial trabeculation were already detectable 48 hours later in the induced knockout of Tie1 plus the heterzygous Tek deletion (starting at E12.5, Tie1ICDiECKO; Tek+/-). This was further confirmed at the neonatal with the induced gene deletion from postnatal day 1 and the analysis 3-weeks later stage (Tie1ICDiECKO; Tek+/-), while TIE1 insufficiency alone (Tie1ICDiECKO) had no obvious impact on the atrial or ventricular trabecula development. TIE1 and TIE2 could form heterodimers, participate in the formation of endothelial cell-cell adherens junction as well as the EC migration and survival 21,22,46,47. Furthermore, TIE1 and TIE2 act in a synergistic manner to specify the vein formation via regulating the protein stability of the venous factor COUPTFII 28,30. As COUPTFII is crucial for the atrial development 31,33, this could provide some insights into the synergy of TIE1/TIE2 in the regulation of cardiac trabeculation especially in atria, likely involving COUPTFII-mediated regulation of the atrial endocardial cell specification.
Differential requirement of TIE receptors in atria and ventricles
So far, little is known about mechanisms underlying the differential regulation of atrial versus ventricular endocardial development. We found in this study that Tie1 and Tek were expressed at higher levels in atrial endocardial cells than those of ventricles by the RNA-seq analysis at the single cell level or the separate analysis of atria and ventricles. By the analysis of scRNA-seq data from embryonic hearts for the ligand-receptor interactions, we showed that the atrial endocardial cells and cardiomyocytes exhibited stronger ANGPT1-TIE2 signaling compared with those in ventricles. This suggests that the differential requirement of TIE1 and TIE2 in cardiac chamber morphogenesis may, at least partly, be related to the differential expression levels of TIE receptors in atrial versus ventricular endocardial cells. In addition, we also found that Dll1 was mainly expressed by endocardial cells from the analysis of single cell RNA-seq data 36, and that Dll1 was among the most downregulated genes in the hearts from Tie1 deficient mice. A similar trend of decreased expression with Notch1 was also detected in atria but not in ventricles of Tie1 deficient mice. As NOTCH1 signaling promotes the extracellular matrix degradation during the formation of endocardial projections that are critical for individualization of trabecular units 13,15, this may also contribute to the differential requirement of TIE1 as well as TIE2 in atrial versus ventricular trabeculation.
In summary, findings from this study show that endothelial TIE1 is crucial for the atrial trabecular formation and play synergistic roles with TIE2 in the process of cardiac chamber morphogenesis. Tie1 deficiency resulted in the significant decrease of Tek, Dll1 and Notch1. This suggests that the differential requirement of TIE receptors in the process of atrial and ventricular trabeculation may involve the DLL1-NOTCH1 pathway. Mechanisms underlying the differential regulation of atrial versus ventricular formation during the heart development warrants further investigation. Knowledge in this direction will facilitate the development of novel interventions for the atria-related cardiomyopathy.
Data availability
RNA sequence data will be deposited in GEO and information provided later.
Acknowledgements
The authors thank the staff in Animal facility of Soochow University for technical assistance.
Additional information
Funding
This work was supported by grants from the National Key R&D Program of China (2021YFA0805000), the National Natural Science Foundation of China (82470518, 82401544), China Postdoctoral Science Foundation (CPSF, 2024M762300), the Natural Science Foundation of Jiangsu Province (BK20240786), the Project of State Key Laboratory of Radiation Medicine and Protection (No. GZN120 20 02), and the Priority Academic Program Development of Jiangsu Higher Education Institutions.
Funding
MOST | National Key Research and Development Program of China (NKPs) (2021YFA0805000)
Yulong He
MOST | National Natural Science Foundation of China (NSFC) (82470518)
Yulong He
MOST | National Natural Science Foundation of China (NSFC) (82401544)
Xudong Cao
China Postdoctoral Science Foundation (2024M762300)
Taotao Li
JST | Natural Science Foundation of Jiangsu Province (Jiangsu Natural Science Foundation) (BK20240786)
Xudong Cao
Project of State Key Laboratory of Radiation Medicine and Protection (GZN120 20 02)
Yulong He
Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD)
Yulong He
Additional files
Supplemental Figures and Tables
References
- 1.The deployment of cell lineages that form the mammalian heartNat Rev Cardiol 15:705–724https://doi.org/10.1038/s41569-018-0086-9PubMedGoogle Scholar
- 2.Importance of myocyte-nonmyocyte interactions in cardiac development and diseaseCirc Res 110:1023–1034https://doi.org/10.1161/CIRCRESAHA.111.243899PubMedGoogle Scholar
- 3.Cardiac chamber formation: development, genes, and evolutionPhysiol Rev 83:1223–1267https://doi.org/10.1152/physrev.00006.2003PubMedGoogle Scholar
- 4.Cardiac trabeculation in vertebrates: Convergent evolution or evolutionary adaptations associated with heart complexity?Seminars in cell & developmental biology 172https://doi.org/10.1016/j.semcdb.2025.103622PubMedGoogle Scholar
- 5.The force within: endocardial development, mechanotransduction and signalling during cardiac morphogenesisDevelopment 143:373–386https://doi.org/10.1242/dev.131425PubMedGoogle Scholar
- 6.Endocardial/endothelial angiocrines regulate cardiomyocyte development and maturation and induce features of ventricular non-compactionEur Heart J 42:4264–4276https://doi.org/10.1093/eurheartj/ehab298PubMedGoogle Scholar
- 7.Formation and Malformation of Cardiac Trabeculae: Biological Basis, Clinical Significance, and Special Yield of Magnetic Resonance Imaging in AssessmentThe Canadian journal of cardiology 31:1325–1337https://doi.org/10.1016/j.cjca.2015.07.003PubMedGoogle Scholar
- 8.Sculpting the heart: Cellular mechanisms shaping valves and trabeculaeCurrent opinion in cell biology 73:26–34https://doi.org/10.1016/j.ceb.2021.04.009PubMedGoogle Scholar
- 9.Embryonic cardiac chamber maturation: Trabeculation, conduction, and cardiomyocyte proliferationAm J Med Genet C Semin Med Genet :157–168https://doi.org/10.1002/ajmg.c.31366PubMedGoogle Scholar
- 10.The deployment of cell lineages that form the mammalian heartNature reviews Cardiology 15:705–724https://doi.org/10.1038/s41569-018-0086-9PubMedGoogle Scholar
- 11.Endocardial cells form the coronary arteries by angiogenesis through myocardial-endocardial VEGF signalingCell 151:1083–1096https://doi.org/10.1016/j.cell.2012.10.023PubMedGoogle Scholar
- 12.Neuregulin 1 sustains the gene regulatory network in both trabecular and nontrabecular myocardiumCirc Res 107:715–727https://doi.org/10.1161/CIRCRESAHA.110.218693PubMedGoogle Scholar
- 13.Notch signaling is essential for ventricular chamber developmentDev Cell 12:415–429https://doi.org/10.1016/j.devcel.2006.12.011PubMedGoogle Scholar
- 14.Activation of Notch1 signaling in cardiogenic mesoderm induces abnormal heart morphogenesis in mouseDevelopment 133:1625–1634https://doi.org/10.1242/dev.02344PubMedGoogle Scholar
- 15.Control of cardiac jelly dynamics by NOTCH1 and NRG1 defines the building plan for trabeculationNature 557:439–445https://doi.org/10.1038/s41586-018-0110-6PubMedGoogle Scholar
- 16.Tunneling nanotube-like structures regulate distant cellular interactions during heart formationScience 387:eadd3417https://doi.org/10.1126/science.add3417PubMedGoogle Scholar
- 17.Extracellular matrix and heart developmentBirth Defects Res A Clin Mol Teratol 91:535–550https://doi.org/10.1002/bdra.20810PubMedGoogle Scholar
- 18.Reduced versican cleavage due to Adamts9 haploinsufficiency is associated with cardiac and aortic anomaliesMatrix biology : journal of the International Society for Matrix Biology 29:304–316https://doi.org/10.1016/j.matbio.2010.01.005PubMedGoogle Scholar
- 19.Myocardial Angiopoietin-1 Controls Atrial Chamber Morphogenesis by Spatiotemporal Degradation of Cardiac JellyCell Rep 23:2455–2466https://doi.org/10.1016/j.celrep.2018.04.080PubMedGoogle Scholar
- 20.Tie2 regulates endocardial sprouting and myocardial trabeculationJCI Insight 5https://doi.org/10.1172/jci.insight.96002PubMedGoogle Scholar
- 21.Control of vascular morphogenesis and homeostasis through the angiopoietin-Tie systemNature reviews Molecular cell biology 10:165–177https://doi.org/10.1038/nrm2639PubMedGoogle Scholar
- 22.Therapeutic targeting of the angiopoietin-TIE pathwayNature reviews Drug discovery 16:635–661https://doi.org/10.1038/nrd.2016.278PubMedGoogle Scholar
- 23.Tie2 regulates endocardial sprouting and myocardial trabeculationJCI Insight 5https://doi.org/10.1172/jci.insight.96002PubMedGoogle Scholar
- 24.Structural basis of Tie2 activation and Tie2/Tie1 heterodimerizationProc Natl Acad Sci U S A 114:4376–4381https://doi.org/10.1073/pnas.1616166114PubMedGoogle Scholar
- 25.The receptor tyrosine kinase TIE is required for integrity and survival of vascular endothelial cellsThe EMBO journal 14:5884–5891https://doi.org/10.1002/j.1460-2075.1995.tb00276.xPubMedGoogle Scholar
- 26.Genetic dissection of tie pathway in mouse lymphatic maturation and valve developmentArterioscler Thromb Vasc Biol 34:1221–1230https://doi.org/10.1161/ATVBAHA.113.302923PubMedGoogle Scholar
- 27.Abnormal embryonic lymphatic vessel development in Tie1 hypomorphic miceDevelopment 137:1285–1295https://doi.org/10.1242/dev.043380PubMedGoogle Scholar
- 28.Endothelial TIE1 Restricts Angiogenic Sprouting to Coordinate Vein Assembly in Synergy With Its Homologue TIE2Arterioscler Thromb Vasc Biol 43:e323–e338https://doi.org/10.1161/ATVBAHA.122.318860PubMedGoogle Scholar
- 29.Loss of flow responsive Tie1 results in ImpairedAortic valve remodelingDev Biol https://doi.org/10.1016/j.ydbio.2019.07.011PubMedGoogle Scholar
- 30.Angiopoietin receptor Tie2 is required for vein specification and maintenance via regulating COUP-TFIIeLife 5:e21032https://doi.org/10.7554/eLife.21032PubMedGoogle Scholar
- 31.The orphan nuclear receptor COUP-TFII is required for angiogenesis and heart developmentGenes Dev 13:1037–1049https://doi.org/10.1101/gad.13.8.1037PubMedGoogle Scholar
- 32.Suppression of Notch signalling by the COUP-TFII transcription factor regulates vein identityNature 435:98–104https://doi.org/10.1038/nature03511PubMedGoogle Scholar
- 33.Atrial identity is determined by a COUP-TFII regulatory networkDev Cell 25:417–426https://doi.org/10.1016/j.devcel.2013.04.017PubMedGoogle Scholar
- 34.TIE1 as a Candidate Gene for Lymphatic Malformations with or without LymphedemaInternational journal of molecular sciences 21https://doi.org/10.3390/ijms21186780PubMedGoogle Scholar
- 35.Loss-of-function mutations of the TIE1 receptor tyrosine kinase cause late-onset primary lymphedemaJ Clin Invest 134https://doi.org/10.1172/JCI173586PubMedGoogle Scholar
- 36.Single-cell transcriptomic analysis identifies murine heart molecular features at embryonic and neonatal stagesNature communications 13:7960https://doi.org/10.1038/s41467-022-35691-7PubMedGoogle Scholar
- 37.Genetic Dissection of Tie Pathway in Mouse Lymphatic Maturation and Valve DevelopmentArterioscl Throm Vas 34:1221–1230https://doi.org/10.1161/Atvbaha.113.302923PubMedGoogle Scholar
- 38.Angiopoietin receptor Tie2 is required for vein specification and maintenance via regulating COUP-TFIIeLife 5https://doi.org/10.7554/eLife.21032PubMedGoogle Scholar
- 39.Neurons limit angiogenesis by titrating VEGF in retinaCell 159:584–596https://doi.org/10.1016/j.cell.2014.09.025PubMedGoogle Scholar
- 40.A Genetically Engineered Mouse Model of Venous Anomaly and Retinal Angioma-like Vascular MalformationBio Protoc 11:e4117https://doi.org/10.21769/BioProtoc.4117PubMedGoogle Scholar
- 41.Single-Cell Transcriptome Atlas of Murine Endothelial CellsCell 180:764–779https://doi.org/10.1016/j.cell.2020.01.015PubMedGoogle Scholar
- 42.The Molecular Signatures Database (MSigDB) hallmark gene set collectionCell Syst 1:417–425https://doi.org/10.1016/j.cels.2015.12.004PubMedGoogle Scholar
- 43.Genetics of atrial cardiomyopathyCurrent opinion in cardiology 34:275–281https://doi.org/10.1097/HCO.0000000000000610PubMedGoogle Scholar
- 44.Myocardium-derived angiopoietin-1 is essential for coronary vein formation in the developing heartNature communications 5:4552https://doi.org/10.1038/ncomms5552PubMedGoogle Scholar
- 45.Distinct roles of the receptor tyrosine kinases Tie-1 and Tie-2 in blood vessel formationNature 376:70–74https://doi.org/10.1038/376070a0PubMedGoogle Scholar
- 46.Angiopoietins assemble distinct Tie2 signalling complexes in endothelial cell-cell and cell-matrix contactsNat Cell Biol 10:527–537https://doi.org/10.1038/ncb1715PubMedGoogle Scholar
- 47.Differential function of Tie2 at cell-cell contacts and cell-substratum contacts regulated by angiopoietin-1Nat Cell Biol 10:513–526https://doi.org/10.1038/ncb1714PubMedGoogle Scholar
Article and author information
Author information
Version history
- Preprint posted:
- Sent for peer review:
- Reviewed Preprint version 1:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.111156. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2026, Ding et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
- views
- 0
- downloads
- 0
- citations
- 0
Views, downloads and citations are aggregated across all versions of this paper published by eLife.