Altered Intercellular Communication and Extracellular Matrix Signaling as a Potential Disease Mechanism in Human Hypertrophic Cardiomyopathy

Study Patients, Sample Collection and Processing

A total of 9 patients with clinically documented HCM scheduled for surgical myectomy were approached for written informed consent to allow their tissue to be used for research. Those who consented underwent myectomy and tissue was collected. Myectomy sample processing was done as previously described ¹⁵. 100 mg of collected myectomy tissue was minced into 1 mm³ pieces, placed in 0.5 mL of CryoStor CS10 Freeze Media (STEMCELL Technologies), and stored in a MrFrosty (ThermoFisher) at 4°C for 10 minutes and then transferred to −80°C overnight. Bulk RNA was isolated from a piece of tissue using the Qiagen RNeasy Plus Micro kit and then assessed on the Agilent Bioanalyzer 2100. Samples with an RNA Integrity Number greater than 8.5 were used in library preparation. Sample collection was approved by the Tufts University/Medical Center Health Sciences Institutional Review Board under IRB protocol # 9487. All subjects gave their informed consent for inclusion before they participated in the study. The study was conducted in accordance with the Declaration of Helsinki. Patient characteristics were obtained from the medical record and are shown in Supplemental Table ST1. Tissue from organ donor patients without underlying cardiac disease was obtained and processed as described previously ¹⁵.

Nuclei isolation, library preparation, and sequencing

Generation of single nuclei sequencing libraries was performed as previously described ¹⁵. Cryopreserved samples were thawed at 37°C and placed on ice. Nuclei were isolated via Dounce homogenization as previously described ¹⁶. Homogenates were filtered through a Pluristrainer 10 µM cell strainer (Fisher Scientific) into a pre-chilled tube. Nuclei were pelleted by centrifuging at 500 x g for 5 min at 4°C. Nuclei pellets were washed and pelleted according to manufacturer protocol (10x Genomics). Nuclei were stained with trypan blue and counted on a hemocytometer to determine concentration prior to loading of the 10x Chromium device and samples were diluted to capture ∼10,000 nuclei. Nuclei were separated into Gel Bead Emulsion droplets using the 10x Chromium device according to the manufacturer protocol (10x Genomics). Sequencing libraries were prepared using the Chromium Single Cell 3’ reagent V2 kit according to manufacturer’s protocol. Libraries were multiplexed and sequenced on a NextSeq550, NovaSeq S2, or NovaSeq S4 (Illumina) to produce ∼50,000 reads per nucleus. Single nuclei RNAseq data for normal IVS tissue is available in the Gene Expression Omnibus database under accession number GSE161921 ¹⁵. Data for HCM IVS tissue is available under accession number GSE174691.

Clustering of Cells by Gene Expression Pattern and Assignment of Cell Type Identity

Sequencing reads were processed using Cell Ranger 6.0.1 ¹⁷. The gene expression matrix was subset to only include reads from the nuclear genome. Quality control (QC) filtering, clustering, dimensionality reduction, visualization, and differential gene expression were performed using the R package Seurat 3.5.0. Each dataset was filtered so that genes that were expressed in three nuclei or more were included in the final dataset. The dataset was further sublet to exclude nuclei that had fewer than 200 genes expressed to remove droplets containing only ambient RNA, and to exclude nuclei with greater than 2000 genes to remove droplets that contained two nuclei. Datasets were individually normalized and integrated using Seurat’s SCTransform development workflow to reduce batch effects ¹⁸. Optimal clustering resolution was determined using Clustree ¹⁹ to identify the resolution where the number of clusters stays stable and was determined to be 0.9 for the integrated dataset. Expression of known cell-specific gene markers were used to identify major cell types. Differentially expressed genes and their functions were used to identify clusters without an assigned identity, or to further refine the cell type.

Trajectory Analysis and Identification of Differentially Expressed Genes

Trajectory analysis was performed using Monocle3 ²⁰ to determine the relationship between subtypes of cells identified in our clustering analysis. Since our data do not represent a developmental time course, we determined the root nodes for each subtype by hierarchical clustering prior to generating trajectories and assigning pseudotime to each cell. Each cell type was analyzed in three cohorts: 1. Normal and HCM cells together; 2. Normal cells alone; 3. HCM cells alone. Differentially expressed genes over trajectory paths in uniform manifold approximation and projection (UMAP) space (i.e., spatial autocorrelation) was performed in Monocle3 using Moran’s I statistic. Moran’s I statistic is a value that varies from −1 to 1, where - 1 indicates perfect dispersion, 0 indicates no spatial autocorrelation, and 1 indicates perfect positive autocorrelation (i.e., nearby cells in have similar gene expression values in focal region of UMAP space). For each Normal and HCM cell type, a gene was determined to be differentially expressed over space if the associated Moran’s I statistic value was positive, paired with a significant adjusted p-value ≤ 0.05, and expressed in ≥ 1% of associated cells. Since many genes showed differential expression over space, further conservative filtering was performed in which genes with Moran’s I statistic available in a single class (i.e., Normal or HCM) were filtered by Moran’s I statistic values >0.1. For genes with Moran’s I statistics available in both classes (i.e., Normal and HCM), genes were filtered by an absolute difference >0.1. Gene Ontology (GO) analysis of molecular function and biological process associated with differentially expressed genes was done using the online tools at uniprot.org/ uniprotkb ²¹.

Analysis of Ligand-Receptor Pair Gene Expression to Discover Intercellular Communication Pathways

To quantify potential cardiac cell-cell communication in Normal and HCM hearts, cell communication networks were plotted in igraph ²² and compared on the basis of ligand-receptor pair gene expression. Our cell-cell communication networks were derived as described previously ¹¹, using a list of 2557 human ligand-receptor pairs ²³. A ligand or receptor for each cell type or cluster was considered expressed if the corresponding gene showed an above zero gene count in ≥ 20% of cells in our snRNAseq data. We initially analyzed the potential signaling interactions between the 10 cell types identified in our snRNAseq data. Lines in our cell networks connect two cell types and represent expressed human ligand-receptor pairs (i.e., potential cell-cell communication between a broadcasting (ligand) and recipient (receptor) cell types. Line color in our networks represents the broadcasting ligand source. Line thickness is proportional to the number of uniquely expressed ligand-receptor pairs. Cell-cell communication networks were also analyzed by fibroblast cluster along with other cell types and also by fibroblast clusters and cardiomyocyte clusters. GO analysis of differentially expressed ligand receptor pairs was performed using the R package clusterProfiler ²⁴.

Statistics

Patient characteristics

Surgical myectomy tissue was obtained from nine HCM patients with severely symptomatic LVOT obstruction. Patient characteristics are summarized in Supplemental Table ST1. Patients gave informed consent for their myectomy tissue to be used in research. The patients varied in age from 37 to 76. Three of nine were female. One of the patients carried a pathogenic Mybpc3 mutation, 1 patient had a likely pathogenic KRAS mutation, and the remaining seven patients had no known mutations pathogenic for HCM. The patients and datasets for the normal human IVS have been previously described ¹⁵.

Single nuclei RNA-sequencing of human HCM myectomy tissue reveals extensive cardiomyocyte and fibroblast diversity in the interventricular septum but does not reveal a disease-specific population of cells

After sequencing and initial data processing with Cell Ranger software ¹⁷, each sample dataset was processed further to remove called nuclei that were likely droplets with only ambient RNA, or droplets that contained two nuclei. The nine HCM datasets and four donor heart datasets were combined into one dataset using the Seurat Integration function ¹⁸. The final combined dataset included 109,823 nuclei from HCM hearts and 24,858 nuclei from donor hearts. Clustering of the integrated dataset revealed 34 cell populations within the IVS, and this was visualized using the dimensionality reduction algorithm UMAP (Fig. 1A), where each point represents a single nucleus colored by cluster identity. Visualizing the integrated dataset by Normal and HCM datasets reveals that all 34 clusters are present in each case, with no cell populations (i.e., clusters) specific to either condition (Fig. 1B).

• Download figure
• Open in new tab

Figure 1. Single nuclei RNA-seq of Normal and HCM IVS Tissue Reveals Cellular Diversity but No HCM-associated Cell Types.

A. UMAP plot of 34 distinct cell clusters identified in the combined Normal and HCM dataset. B. Separating the cell clusters by disease label does not identify HCM-specific cell types. C. Listing of cell markers used to assign cell types. D. UMAP plot of combined dataset showing cell identity labels instead of cluster numbers.

Cell identities were assigned to each cluster using known gene markers of expected cell types, differentially expressed genes, gene ontology using GOstats ²⁵, and Ingenuity Pathway Analysis ²⁶. Figure 1C lists the gene markers used to identify each cell type and illustrates each cluster’s unique expression of gene markers. Upon assigning cell types, similar cell types were expectedly positioned close to each other in UMAP space (Fig. 1D). Interestingly, we see thirteen separate cardiomyocyte populations and eight different fibroblast populations, revealing significant cardiomyocyte and fibroblast diversity in the IVS (Fig. 1C, 1D). Other cell types identified included endothelial, lymphatic, smooth muscle, pericyte, neuronal, stromal, and immune cell populations (macrophages, lymphocytes). Among normal donor heart samples, assigned cardiomyocytes made up approximately 41.0% (mean, range 30.1 to 61.6%) of the total cell population. Similarly, among HCM myectomy samples cardiomyocytes made up approximately 51.6% (mean, range 35.5-63.0%, p-value not significant) of the total cell population. Noncardiomyocytes make up the remaining cell population. HCM does not appear to be associated with a shift in the relative proportion of cells in the heart.

Trajectory and differential gene expression analysis reveals HCM-associated perturbations

Clustering analysis revealed great diversity for both cardiomyocytes and fibroblasts. To gain insight into potential relationships among the different cell populations for each cell type, clusters were projected onto pseudotime trajectory analysis using Monocle3 ²⁰. Trajectory analysis was performed on Normal cells only, HCM cells only, and both conditions combined for each of the 10 cell types identified in the above analysis. In single-cell trajectory analysis, a trajectory is a computed path that describes a cell type’s biological progression through a dynamic process. Prior to building trajectory paths, the root node, or beginning, of each trajectory path was determined through hierarchical clustering of our single nuclei data. Then, trajectory paths for each assigned cell type (assigned in Seurat) were constructed in UMAP space using Monocle3 with Normal and HCM data together and individually. Once trajectory paths were established, pseudotime could be assigned to each cell. For all cell types, root nodes showed consistent placement among Normal only groups, HCM only groups, and Normal and HCM groups. Trajectory paths also showed similar basic structures among Normal only groups, HCM only groups, and Normal and HCM groups. Small differences in trajectory paths for a single cell type among groups are likely due to different sample sizes in Normal and HCM groups (e.g., Normal cardiomyocyte nuclei n=10131, HCM cardiomyocyte nuclei n=57847). Therefore, we were not able to distinguish any meaningful changes among trajectory paths. Representative trajectories and root nodes are shown for the two most common cell types, cardiomyocytes and fibroblasts, in Supplemental Figure S1. These findings suggest that the relationships between the various subtypes of each cell type do not vary significantly in HCM heart tissue in comparison to Normal heart tissue.

Establishment of trajectories facilitates the analysis of differential gene expression along the trajectory and between conditions along the trajectory. We analyzed differential gene expression for each cell type along their trajectories for the Normal and HCM populations through spatial autocorrelation. Moran’s I statistic and adjusted p-values were calculated for each gene in Normal only and HCM only groups within each cell type. When paired with a significant p-value (≤0.05), a Moran’s I statistic value near zero indicated no spatial autocorrelation and a value near 1 indicated perfect positive autocorrelation in which a gene is expressed in a focal region of the UMAP space. Results showed between 91 and 3589 genes are differentially expressed along trajectory paths among cell types and between 4 and 1823 differentially expressed genes overlapped between Normal and HCM groups among cell types (Supplemental Table ST2).

Since many genes showed differential expression over space, further conservative filtering was performed to identify genes of interest. Genes with Moran’s I statistic available for a single condition, Normal or HCM, were filtered by a Moran’s I statistic value above 0.1, while genes with Moran’s I statistics available for both classes were filtered by an absolute difference >0.1. For each cell type, 0 to 41 genes passed the filter with 176 unique genes in total passing the filter among all cell types (Supplemental Table ST3). The expression of these filtered genes was then plotted in UMAP space for their associated cell type in Normal and HCM conditions as shown for a subset of identified genes in cardiomyocytes (Supplemental Figure S2). Visual analysis of the UMAP plots revealed 28 genes with pronounced differences in their spatial expression between Normal and HCM conditions (summary in Table 1). Of these 28 genes, 6 (ACTA1, ACTC1, MYHL7, MYL2, TNNI3, TNNI2) are already known to be associated with human HCM. Extracellular matrix genes (COL1A1, COL6A1, COL6A2, CYR61, POSTN) are also differentially expressed in Normal vs HCM fibroblasts, suggesting a change in ECM composition. GO Enrichment analysis of these 28 differentially expressed genes revealed significant enrichment for molecular functions involving the ECM, biological processes involving muscle development and cellular components involving the ECM (Supplemental Fig. S3A). Analysis of genes from this list showing increased expression in obstructive HCM revealed significant enrichment in molecular functions involving actin binding, biological processes involving muscle filament sliding, contraction and movement and cellular components involving the sarcomere, myofibril and contractile fiber (Supplemental Fig. S3B). Analysis of genes from this list showing decreased expression in HCM revealed enrichment for molecular functions involving the ECM, biological processes involving the response to transforming growth factor-β (TGF-β), steroid hormones, the immune response and ECM organization, and cellular components involving the ECM (Supplemental Fig. S3C).

Analysis of ligand-receptor pair gene expression indicates profound alteration of intercellular communication in HCM

HCM is associated with a general decrease in cell-cell communication but an increase in fibroblast to lymphocyte communication

To quantify potential cardiac cell-cell communication in Normal and HCM IVS tissue, we quantified the number of possible expressed ligand-receptor pairs among cell types as previously described ¹¹. We examined the expression of a curated list of 2557 human ligand-receptor pairs ²³ in which ligands and receptors were considered expressed if their associated gene was detectable in ≥ 20% of cells in a cell type. Notably, quantification of cell-cell communication in this study represents potential communication as we account for only expressed ligand-receptor pairs and not the position or boundaries of cell types. Results showed that Normal cardiac IVS tissue demonstrated an overall dense intercellular communication network among our 10 identified cell types (Fig 2). HCM cardiac tissue, however, unexpectedly demonstrated an overall sparse intercellular communication network in which potential communication was drastically reduced compared to the Normal tissue (Fig 2A). Quantitatively, the total number of expressed ligand-receptor pairs in the HCM tissue was more than two times less than in Normal tissue (546 and 1138 respectively, Fig. 2B). Analysis of the broadcast ligands by individual cell types demonstrates that 1) cardiomyocytes show a broad decrease in both paracrine and autocrine ligand broadcasting in HCM (Fig. 2C, D) and 2) fibroblasts broadcast the greatest number of ligands to the other cell types in both Normal and HCM tissue. Fibroblast communication is particularly high with fibroblast, endothelial, pericyte, smooth muscle, neuronal, and cardiomyocyte cells (Fig. 2). While intercellular communication was generally reduced in HCM tissue compared to Normal tissue, there was one case of increased HCM communication from fibroblasts to lymphocytes (1 expressed ligand-receptor pairs in Normal tissue to 13 in HCM tissue; Fig 2C, D).

• Download figure
• Open in new tab

Figure 2. Intercellular communication networks are reduced in HCM.

A. Cell-cell communication networks between cardiac cell types in normal control (left) and HCM (right) conditions. Line color indicates ligand broadcast by the cell population with the same color. Lines connect to cell types which expressed cognate receptors. Line thickness is proportional to the number of uniquely expressed ligand-receptor pairs. Loops indicates communication within a cell type. B. Quantity of ligands and receptors in expressed ligand-receptor pairs described by cell type and condition (Normal or HCM). C., D. Cell-cell communication networks broken down by cell type in normal control (C) and HCM (D) conditions. Figure formatting follows panel A. Numbers indicate the quantity of uniquely expressed ligand-receptor pairs between the broadcasting cell type (expressing ligand) and receiving cell type (expressing receptor).

To assess the molecular functions potentially affected by changes in ligand-receptor signaling among our 10 identified cell types, GO enrichment analysis was performed on the ligands and receptors in expressed ligand-receptor pairs from both Normal and HCM tissue (Fig. 3). We observed a global decrease in nearly all identified ligand functions in the HCM heart tissue compared to the Normal heart tissue, particularly in functions such as extracellular matrix structural constituent (55% total ligand count decrease), extracellular matrix structural constituent conferring tensile strength (54% total ligand decrease), growth factor binding (63% total ligand count decrease), protease binding (46% total ligand count decrease), and platelet-derived growth factor binding (63% total ligand count decrease; Fig. 3A). As expected, the molecular functions associated with changes in receptor expression mirrored those observed for changes in ligand expression (Supplemental Fig. S4A). HCM cell types including cardiomyocytes, endothelial, macrophages and lymphocytes showed disproportionately large decreases in the above listed ligand processes ranging from 67-100% decreases in total ligand counts associated with a process compared to Normal ligand counts (Figure 3B-E). Two of these ligand processes (growth factor binding and platelet-derived growth factor binding) were completely lost in the HCM condition in cardiomyocytes, endothelial, macrophages and lymphocytes (Figure 3B-E). These trends found in ligand GO enrichment analysis were also identified in receptor GO enrichment analysis, except in lymphocytes (Supplemental Fig. S4B-F). In Normal tissue, lymphocyte receptors do not overrepresent any GO molecular functions, but in HCM tissue they show overrepresentation of receptors associated with coreceptor activity, cytokine binding, exogenous protein binding, integrin binding, etc., suggesting important changes in lymphocyte function. Together, ligand and receptor findings suggest that growth, the extracellular matrix, protein breakdown and lymphocyte function may be highly affected in HCM cardiac cells. Since fibroblast to lymphocyte signaling was the only pathway associated with increased ligand-receptor interactions in HCM heart tissue, we examined the specific ligand-receptor interactions in fibroblast to lymphocyte signaling that were increased (Supplemental Table ST4). In HCM, lymphocytes gained ITGB1 receptor expression, allowing communication with several cognate ligands expressed in fibroblasts (COL1A1, COL1A2, COL3A1, COL6A1, COL6A2, COL6A3, FBLN1, FN1, HSPG2, VCAN). HCM fibroblasts also gained expression of COL1A1 and COL1A2 ligands which enabled communication with the CD36 lymphocyte receptor. These findings further underscore the importance of the extracellular matrix in HCM.

• Download figure
• Open in new tab

Figure 3. Bar plot representing the total count of ligands (in expressed ligand-receptor pairs) associated with different cellular processes in Normal and HCM IVS Cells.

Bar color distinguishes ligand count in normal or HCM conditions. A. Comparison of molecular processes across all cell types. B. Comparison in cardiomyocytes. C. Endothelial cells. D. Macrophages. E. Lymphocytes.

Fibroblast subtypes show subtype-specific changes in cell-cell communication

To better understand the interaction between the different fibroblast populations and other cell types in the IVS, we examined the potential ligand-receptor communication through cell network plots for all fibroblast subtypes alongside the other previously identified 9 cell types (Fig. 4A, C, D). In Normal tissue, all fibroblast subtypes broadcast ligands extensively. In HCM tissue, fibroblast subtypes show a drastically reduced number of broadcasting ligands compared to Normal tissue, (1694 versus 4573 respectively; Fig. 4). There was a disproportionate reduction in ITGB1 and LRP1 receptor expression in fibroblast clusters 7 and 8 (Fig. 4B), which disabled communication with several ligands (ITGβ1 to COL1A1, COL1A2, COL3A1, COL6A1, COL6A3, FBLN1, FN1, HSPG2, and VCAN ligands; LRP1 to APP, C3, CALR, CTGF, HSP90B1, LRPAP1, PSAP, SERPING1, SERPINE2, TFPI).

• Download figure
• Open in new tab

Figure 4. Cell-cell communication networks between fibroblast subtypes and other heart cells in normal control and HCM conditions.

A. Comparison of Normal (left) and HCM (right) communication networks. Line color indicates ligand broadcast by the cell population with the same color. Lines connect to cell types which expressed cognate receptors. Line thickness is proportional to the number of uniquely expressed ligand-receptor pairs. Loops indicates communication within a cell type. B. Percent expression of receptors ITGB1 and LRP1 by fibroblast clusters (7 and 8) in normal and HCM conditions. C, D. Cell-cell communication networks broken down by cell type and fibroblast cluster in normal control (C) and HCM (D) conditions. Figure formatting follows panel A. Numbers indicate the quantity of uniquely expressed ligand-receptor pairs between the broadcasting cell type (expressing ligand) and receiving cell type (expressing receptor).

Cardiomyocyte subtypes and fibroblast subtypes together show cell subtype-specific alterations in cell-cell communication

Given the importance of sarcomere dysfunction in HCM and the observed reduction in ligand broadcasting of HCM cardiomyocytes (Fig. 2D), we examined the cell communication networks between fibroblast subtypes (i.e., clusters, n=8) and cardiomyocyte subtypes (n=13), in Normal and HCM tissue (Fig. 5). The total number of expressed ligand-receptor pairs in the Normal tissue was 2.6 times greater than in the HCM tissue (6975 and 2683 expressed ligand-receptor pairs respectively). Compared to the Normal tissue, HCM communication expectedly decreased drastically in most pathways (415 of 441 total pathways, Figure 5A-C). More so, 10.8% of pathways decreased to zero communication in HCM tissue. All zeroed communication pathways involved losses normally broadcasted from cardiomyocytes (clusters 2-5, 7, 11; Figure 5C, Supplemental Table ST5). Zeroed expression ligands (ligands that fell below the 20% expression threshold) in HCM cardiomyocytes included TIMP1, COL1A2, COL3A1, FN1, LAMA2, COL1A2, CALM2, MFGE8, LPL, PSAP, SERPING1, COL6A1, COL6A2, VEGFA, SORBS1 and ADAM17. GO analysis of molecular functions associated with these ligands included platelet-derived growth factor binding, extracellular matrix structural constituent conferring tensile strength, collagen binding, protease binding and integrin binding (data not shown).

• Download figure
• Open in new tab

Figure 5. Cell-cell communication networks between cardiac fibroblast and cardiomyocyte subtypes in normal control and HCM conditions.

A. Overall communication network between cardiomyocytes and fibroblasts. Line color indicates ligand broadcast by the cell population with the same color. Lines connect to cell types which expressed cognate receptors. Line thickness is proportional to the number of uniquely expressed ligand-receptor pairs. Loops indicates communication within a cell type. B, C. Cell-cell communication networks broken down by fibroblast cluster and cardiomyocyte cluster in normal control (B) and HCM (C) conditions. Numbers indicate the quantity of uniquely expressed ligand-receptor pairs between the broadcasting cell type (expressing ligand) and receiving cell type (expressing receptor).

The largest decreases in ligand-receptor pairs between Normal and HCM communication pathways (decreases >30 ligand-receptor pairs; Figure 4) involved fibroblast cluster 4 as fibroblast cluster 4 ligands were lost in HCM tissue compared to normal tissue (C3 CALR, COL18A1, COL4A1, COL5A1, COL5A2, CXCL12, CYR61, FBLN1, FBN1, HSP90B1, IGF1, LAMA2, LAMB1, LAMC1, LGALS3BP, LRPAP1, NID1, PDGFD, PTN, TFPI, THBS2, TIMP2, TNC; Supplemental Table ST6). The reduction in L-R pairs in fibroblast cluster 4 affected communication with other fibroblast clusters (Supplemental Figure ST7). The lost ligands are involved in vital molecular functions including integrin binding, extracellular matrix structural constituent, and extracellular matrix structural constituent conferring tensile strength among others, thereby suggesting fundamental alterations in these functions in HCM (Supplemental Table ST6).

While intercellular communication generally decreased in HCM tissue compared to Normal tissue, some specific cluster communication paths showed an opposite trend with increased HCM communication. HCM communication increased 3.2 times from fibroblasts clusters 2-5 to cardiomyocyte clusters 1, 2 and 13 compared to Normal tissue (Fig. 5c). This increase in potential communication can be attributed to 1) expression of the receptor ITGβ1 in cardiomyocyte clusters 1, 2 and 13 in HCM cells only, and 2) expression of several of ITGβ1’s cognate ligands in HCM fibroblast clusters (COL3A1, COL6A1, COL6A2, COL6A3, FBLN1, FBN1, FN1, HSPG2, LAMA2, LAMB1, LAMC1, LGALS3BP, TIMP2, VCAN). These cognate ligands were not expressed in Normal fibroblast clusters 2-5 (Supplemental Tables ST8). Another example of increased HCM communication (2.3 fold increase) relative to Normal communication involves communication from cardiomyocyte cluster 8 to cardiomyocyte clusters 1, 2 and 13 (Supplemental Table ST9).This increased communication can again be attributed to 1) expression of the of the receptor ITGβ1 in receiving clusters and 2) expression of several of ITGβ1’s cognate ligands in the broadcasting cardiomyocytes in cluster 8 (COL3A1, COL6A1, COL6A2, FN1, HSPG2, LAMA2, LGALS3BP). These cognate ligands were not expressed in Normal fibroblasts in cluster 8.