Clinical PerspectiveWhat Is New?
What Are the Clinical Implications?
IntroductionThe mammalian heart is composed of a complex and heterogeneous network of cells. Although cardiomyocytes (CMs), the contractile cells of the heart, constitute the majority of cardiac cell mass, nonmyocytes far outnumber CMs.1 Little is known about how the individual cellular components of the heart operate together as an integrated unit, how this cellular ecosystem is altered in the context of physiological stress, and how disparate cell populations may contribute to the functional decline that accompanies cardiac injury, aging, hypertension, or obesity. Moreover, although biological sex is recognized as an important variable in cardiovascular homeostasis and disease,2 little is known about how biological sex affects transcriptional responses of distinct cardiac cell populations to chronic tissue stress. To study the cardiac cellulome, the network of cells that forms the heart, in the context of chronic stress, we used an angiotensin II (AngII) cardiac hypertrophy model and applied a novel method for simultaneous transcriptional profiling of single CM nuclei and nonmyocyte cells. Our study revealed shared and distinct cellular pathways driving AngII-induced cardiac fibrosis. Early fibrosis developed independently of smooth muscle actin–expressing (ACTA2+) myofibroblasts, a cell type conventionally recognized as a key mediator of cardiac fibrosis. This study also documents extensive sexual dimorphism in the gene profiles of cardiac cells, underscoring how biological sex influences responses to tissue stress. This work provides a valuable resource for cardiac cell biology research and offers important insights into the orchestrated cellular and molecular mechanisms that drive cardiac fibrosis and heart failure. MethodsFor a detailed description of methods see Expanded Methods in the Data Supplement. All animal procedures were in accordance with institutional guidelines and approved by The Jackson Laboratory Institutional Animal Care and Use Committee or La Trobe University Animal Ethics Committee (AEC16-93). All sequencing data have been made publicly available at ArrayExpress and can be accessed by using accession number E-MTAB-8810. Statistical AnalysisThe methods used are described in detail in Expanded Methods in the Data Supplement and in figure legends (where appropriate). All P values considered statistically significant were ≤0.05. ResultsPrevious large-scale single-cell transcriptomic studies exploring cardiac cellular diversity have failed to profile key cardiac cell populations including CMs 3–6 because of a reliance on droplet-based sequencing approaches that are unable to process large cells. Moreover, extraction of CMs from rodent hearts typically requires retroaortic perfusion with proteases to liberate cells from the extracellular matrix (ECM); this is a relatively time-consuming process that has limited applicability for cell preparation in transcriptomic studies.7 To overcome these limitations, we developed a novel experimental framework to isolate and prepare all cardiac cells for single-cell RNA sequencing (scRNA-seq; Figure 1A). The goal of this approach was to profile all cardiac cell types with maximum possible throughput and accuracy, subject to limitations based on biological characteristics and current scRNA-seq instrumentation. Building on recent innovations in CM isolation,8 we developed a perfusion-based tissue dissociation protocol that enables simultaneous isolation of cardiac cells from multiple hearts (see the Methods section; Figure 1A). To address the challenge of CM size, CMs were denuded of their cell bodies to liberate nuclei for droplet-based transcriptional profiling. As proof of concept, both CM nuclei and nonmyocyte nuclei were isolated independently and pooled for sequencing (Figure I in the Data Supplement). Isolation of cells using this methodology yielded similar nonmyocyte cell proportions to protocols previously used,4 although total cell yield and the number of extracted endothelial cells was lower (Figure II in the Data Supplement). Although isolation of nuclei overcame critical limitations related to CM size, for all further experimentation the nonmyocyte fraction was preserved as whole cells to provide access to total cellular RNA and to allow control of input cell proportions for scRNA-seq by fluorescence-activated cell sorting using cell surface markers (Figure 1A). After isolation of both CM nuclei and nonmyocyte cells, scRNA-seq was performed on a mixture of myocytes and nonmyocytes from untreated female and male hearts. We limited the relative number of endothelial cells and CM nuclei to minimize oversampling of these abundant cell populations at the expense of resolving less frequent cell types (Figure 1A). After the clustering of 7474 cells derived from 2 mice of each sex (Figure 1B), examination of established marker genes and top enriched genes in each cell population (Figure 1C and 1D; Table I in the Data Supplement) revealed a wide array of cell types that were present in all samples (Figure III in the Data Supplement). These include fibroblasts (Pdgfra, Col1a1), pericytes (Pdgfrb, Vtn), smooth muscle cells (Acta2, Myh11), Schwann cells (Plp1, Kcna1), endothelial cells (Pecam, Ly6c1), macrophages (Fcgr1, Csf1r), and other immune cell populations (granulocytes, B cells, T cells, and natural killer cells). A subset of fibroblasts (Fibroblast-Wif1) showed distinct gene patterns including active Wnt signaling and specific of Wif1, in line with our previous report4 and the findings of Farbehi et al.5 Additional less well-characterized cell populations included CMs (Ttn, Myh6), epicardial cells (Msln, Upk3b), lymphatic endothelial cells (Lyve1, Cldn5), and endocardial cells (Npr3, Cytl1; Figure 1B through 1D). In this study we examined pathological remodeling of the heart after 2 weeks of chronic AngII administration, which corresponds to a relatively early stage of disease development (Figure 2A). An increase in heart mass and fibrotic area was noted in both female and male mice treated with AngII; however, hemodynamic changes and functional impairment were principally restricted to male mice that had a reduction of ejection fraction (Figure IV in the Data Supplement), in line with trends in human populations.9 Histological analyses of hypertrophic mouse hearts showed a greater cross-sectional area and higher levels of fibrosis in both female and male hearts (Figure 2B and 2C; Figure V in the Data Supplement). Together these observations describe the pathological remodeling of the heart with AngII-induced pathological hypertrophy in both female and male contexts. To examine cellular and molecular changes that accompany induction of fibrosis, we performed scRNA-seq on isolated myocytes and nonmyocytes in a total of 12 hearts from sham- and AngII-treated mice of both sexes. Combining these data with cells from the 4 untreated female and male hearts described earlier, we obtained a total of 29 558 cells and performed clustering to identify a diverse array of cardiac cell types and subtypes (Figure 2D; Table II in the Data Supplement). Examination of cell abundances within this scRNA-seq data set revealed several cell populations that changed in relative prevalence in response to AngII (Figure 2E and 2F). Most notably, 2 fibroblast subpopulations that were not present at an appreciable number in unstressed hearts (Figure 2E and 2F; Figure VI in the Data Supplement) increased dramatically in hearts of mice treated with AngII and expressed high levels of periostin (Figure 2E). We refer to these cell populations as Fibroblast-Cilp and Fibroblast-Thbs4, reflecting markers that are highly expressed within these populations, with Fibroblast-Cilp being the largest cell population in this fibrosis data set. Subclustering fibroblasts did not significantly fragment fibroblast subpopulations further, whereas subclustering endothelial cells and dendritic-like cells yielded subpopulations that shifted in abundance with AngII treatment (Figure VII in the Data Supplement). The limited input of specific cell populations for scRNA-seq analysis might restrict the detection of shifts in abundances of certain cell populations. To examine cell composition using an orthogonal method, we performed flow cytometric analysis of control and AngII-treated mouse hearts (Figure VIII in the Data Supplement). Macrophages significantly increased in abundance in males, compared with controls. However, there was no significant increase in fibroblasts detected by flow cytometry (Figure VIII in the Data Supplement), suggesting that Fibroblast-Cilp and Fibroblast-Thbs4 represent a cell state arising primarily from resident fibroblasts rather than infiltration or proliferation of cells. Single-cell transcriptional profiling of cell cycle states revealed no significant increase in proliferation of fibroblast populations or in other cell types in response to AngII (Figure IX in the Data Supplement). AngII treatment induced gene changes in all cardiac cell populations (Figure 3A; Figure X in the Data Supplement; Table III in the Data Supplement). Most notably, fibroblast populations up- or downregulated the greatest number of genes, with macrophage and endothelial cell subsets also exhibiting high responsiveness. Examination of the top 10 genes upregulated in response to AngII within each cell population revealed several patterns relevant to hypertrophy (Figure 3B; Table III in the Data Supplement). Among the top genes upregulated on AngII treatment were transcripts for Cilp, Comp, Crlf1, Postn, and Thbs4, some of which have been previously implicated in human heart failure.10–12 These genes tend to be most highly expressed in only a few cell types but exhibit a significant increase in across many cell types under AngII treatment. Genes that were upregulated in clusters of cells or across a limited number of cell populations included Ms4a7 (macrophages) and Gm13889 (pericytes and smooth muscle cells). Upregulated genes highly restricted to 1 cell type included Ccl21a (in lymphatic endothelial cells), and Nppa and Nppb (in CMs), which are well-established biomarkers of heart failure.13 Gene ontology (GO) analysis of AngII-upregulated genes revealed that the top 3 enriched GO terms in most cardiac cell types contribute to remodeling of the cardiac ECM (Figure 3C; Table IV in the Data Supplement). Among the top genes downregulated in response to AngII, we observed a number of cell type–specific downregulated genes, in particular, in epicardial cells and multiple genes related to ATP synthesis and mitochondrial activity (Figure XIA in the Data Supplement; Table III in the Data Supplement). Downregulated genes also formed coherent gene programs, most notably in fibroblasts, related to protein translation and inflammation (Figure XIB in the Data Supplement; Table IV in the Data Supplement). Cellular abundances of Fibroblast-Cilp and Fibroblast-Thbs4 changed most dramatically on treatment with AngII (Figure 2F). To determine the effect of these cell populations on cardiac remodeling, we first compared differentially expressed genes between Fibroblast-Cilp and Fibroblast-Thbs4 collectively with all other control cells within this data set, cardiac fibroblasts from control animals, and fibroblasts from AngII-treated animals (Figure 4A through 4C, respectively; Tables V and VI in the Data Supplement). Features common to basic fibroblast biology defined Fibroblast-Cilp and Fibroblast-Thbs4 in comparison with all cardiac cells (eg, collagen and ECM organization). When contrasted with all fibroblasts from control animals or AngII-treated cohorts (Figure 4B and 4C, respectively), collagen-remodeling genes were the primary distinguishing features of Fibroblast-Cilp and Fibroblast-Thbs4. Comparison of Fibroblast-Cilp and Fibroblast-Thbs4 also showed higher levels of ECM-remodeling genes in Fibroblast-Cilp than in Fibroblast-Thbs4 (Figure 4D; Tables V and VI in the Data Supplement). Examination of total transcripts corresponding to ECM-remodeling genes suggests that Fibroblast-Cilp and Fibroblast-Thbs4 are key contributors to ECM remodeling, with the Fibroblast-Cilp population contributing greater summed transcript abundance than all other fibroblast populations in AngII-treated mice (Figure 4E; Figure XII in the Data Supplement). Fibroblast-Cilp and Fibroblast-Thbs4 were also distinguished by the upregulation of genes such as Thbs4, Fmod, Cthrc1, and Cilp2 (Figure 4F; Table V in the Data Supplement). No genes were uniquely and highly expressed in Fibroblast-Cilp that were not also expressed in Fibroblast-Thbs4. THBS4+ cells in AngII-treated hearts were aggregated within foci and could also be identified as discrete cells (Figure 4Gi and 4Gii, respectively; Movie I in the Data Supplement; Figure XIII in the Data Supplement). An absence of THBS4+ cells surrounding perivascular fibrotic lesions (Figure 4Giii) suggested that Fibroblast-Thbs4 are principally involved in interstitial fibrosis. Last, in accordance with transcriptomic analyses, we confirmed that neither Fibroblast-Cilp nor Fibroblast-Thbs4 corresponds to ACTA2+ myofibroblasts (Figure 4H and 4I). These findings suggest that Fibroblast-Cilp and Fibroblast-Thbs4 may be transcriptionally similar to the matrifibrocyte, a specialized fibroblast state found in the mature scar after myocardial infarction.14 Some of the matrifibrocyte markers Comp, Sfrp2, and Wisp26,14 was observed in Fibroblast-Cilp/Thbs4 cells, although these genes were expressed at varying levels within the cluster, and cells with high of each marker did not colocalize. These genes are less specific for Fibroblast-Cilp and Fibroblast-Thbs4 than Cilp, Ddah1, or Thbs4 (Figure XIV in the Data Supplement). To explore the origins of Fibroblast-Cilp and Fibroblast-Thbs4, we examined dynamic patterns of transcriptional activation and repression in detail in all fibroblasts using RNA velocity.15,16 This analysis suggested that the Fibroblast-Thbs4 population predominantly arises from Fibroblast-Cilp, whereas Fibroblast-Cilp has broader origins from Fibro5 and Fibro6 (Figure 5A). Use of 2 independent pseudo–time-based methods yielded results that were also concordant with these conclusions (Figure XV in the Data Supplement). Some transcripts expressed in both Fibroblast-Cilp and Fibroblast-Thbs4, such as Cilp and Postn, showed gradients in velocity across Fibroblast-Cilp into Fibroblast-Thbs4 (Figure 5B and 5C), suggesting that Fibroblast-Thbs4 cells represent a cell state characterized by transcriptional activation of a unique set of genes. To determine the effect of AngII treatment on cardiac intercellular signaling, we mapped ligands and cognate receptors expressed by various cell populations (Figure 6; Table VII in the Data Supplement). As previously reported,4 a similar analysis revealed extensive intercellular communication in unstressed hearts. Signaling by fibroblasts, along with macrophages, represented key features of the interstitial cardiac niche at baseline (Figure 6A; Figure XVI in the Data Supplement). Dendritic-like cells, epicardial cells, and endocardial cells were also significant populations forming extensive potential signaling networks within a diverse array of cell types (Figure 6A; Figure XVI in the Data Supplement). In contrast, consistent with previous analyses,4 patrolling cell populations, granulocytes, natural killer cells, and lymphocytes, were among the least interactive cell populations. Few interactions observed for CMs may be attributable to a lower number of genes and reads in this cell population attributable to transcriptional profiling of nuclei rather than whole cells. Examination of factors essential for supporting specific cell populations points to important trophic intercellular interactions (Figure 6B). In addition to previously reported patterns for Igf1, Il34, Ngf, and Ntf3 (Figure 6B),4 epicardial cells were identified as the most highly enriched for Csf1 transcripts, required for macrophage vitality, whereas CMs were the primary producers of Vegfa that supports endothelial cell growth (Figure 6B). Although the role of the epicardium in supporting macrophage trafficking in fetal development has been previously established,17 this result suggests that epicardial cells are the most macrophage-trophic cell population by Csf1/Csf1r signaling (Figure 6B).The spatial distribution of macrophage subsets relative to the epicardium is distinct (Figure 6C), and further accentuated in older animals,18 underscoring the epicardium and subepicardial space as a macrophage-trophic environment. Cell–cell signaling in response to AngII-induced hypertrophy was increased in all cell types studied, with fibroblasts increasing the greatest number of connections (Figure 6D; Table VIII in the Data Supplement). GO enrichment analysis of upregulated connections between cell types revealed that the most frequently increased ligands involved cellular response to amino acid stimulus and collagen fibril organization (Figure 6E; Table IX in the Data Supplement); upregulated receptors included those interacting with ECM (Figure 6F; Table IX in the Data Supplement). Ligands involved in angiogenesis were upregulated in multiple cell types, whereas endothelial cells upregulated cognate receptors (Figure 6F; Table IX in the Data Supplement). Together these analyses suggest a concerted effort by multiple cell types in the cardiac cellulome to alter the cardiac niche and promote ECM remodeling. Conversely, GO terms related to downregulated ligands suggests dampening of inflammation in the cardiac microenvironment (Figure XVII in the Data Supplement; Table IX in the Data Supplement), with almost all enriched GO terms related to inflammation or leukocyte/monocyte trafficking. It is intriguing that other GO terms enriched within downregulated ligands and receptors were associated with angiogenesis, suggesting that the up- and downregulation of receptors within endothelial cells may reflect changes in cellular phenotype and shifts in sensitivity to angiotrophic cues. Analysis of hemodynamic and cardiac functional parameters, and quantification of tissue fibrosis, as well, suggested extensive sexual dimorphism in the development of pathological remodeling in response to AngII. Many upregulated genes are modulated by estradiol (Figure XVIII in the Data Supplement), and sexually dimorphic differences in cardiac cellularity after AngII treatment were detected by flow cytometry (Figure VIII in the Data Supplement). Moreover, fibrosis developed at specific anatomic loci to varying degrees in a sex-dependent manner (Figure 7A). To explore sex-specific gene patterns within specific populations of the cardiac cellulome, we isolated female and male cells from control and AngII-treated groups based on the of genes expected to differ between sexes (chromosome X and Y genes outside the pseudoautosomal region; Figure 7B; Figure XIX in the Data Supplement). Comparison of proportions of female and male cell populations in control and AngII-treated hearts revealed no differences in cell abundance, with the exception of the Fibroblast-Thbs4 population, where we observed ≈2-fold higher abundance of these cells in females than in males (Figure 7C). Comparison of sex-specific gene differences within cell populations from both control and AngII-treated hearts revealed highly sexually dimorphic gene in almost all cell populations, with the greatest level of sexually dimorphic gene occurring within fibroblast subpopulations 2 and 4 to 6 in control hearts, and fibroblast subpopulations 2 to 6 and Fibroblast-Cilp of AngII-treated mice (Figure 7D and 7E; Tables X and XI in the Data Supplement). This is in line with the high enrichment of sex hormone receptors expressed by these cells (Figure XX in the Data Supplement). Sexual dimorphism in the of genes upregulated after AngII treatment (Figure 7E and 7F; Figures XXI and XXII in the Data Supplement; Table XI in the Data Supplement) was predominantly detected in fibroblast populations (with the exception of Fibroblast-Wif1 and Fibroblast-Thbs4), with very few genes differentially regulated in other cell types (Figure 7E; Table XI in the Data Supplement). Dimorphisms in were less evident in genes downregulated after cardiac fibrosis (Figure XXIII in the Data Supplement). Examination of genes with ≥2-fold difference between sexes in cells of control hearts revealed a heterogeneous collection of genes (Figure XXIVA in the Data Supplement; Table X in the Data Supplement). These included genes encoding proteins such as GNGT2, CAMLG, TDP1, and NEK7 whose orthologs have been previously identified as sexually dimorphic in human heart failure.19 GO enrichment analysis of sexually dimorphic genes in control hearts showed an enrichment of genes corresponding to antigen processing in multiple female cell populations and protein folding, stress, and angiogenesis in multiple male cell types (Figure XXIVB in the Data Supplement). Consistent with sexually dimorphic regulation of the endothelium, flow cytometric analysis of endothelial cells showed greater numbers of endothelial cells in male mouse hearts (Figure XXV in the Data Supplement). Genes such as Vegfa were upregulated in male (but not female) CMs and Fibroblast-Thbs4 populations from AngII-treated mice, further suggesting sexually dimorphic angiogenic stimulation (Figure 7E). Among genes upregulated or downregulated in either sex (Figure XXII in the Data Supplement), female hearts had an enrichment of ECM-remodeling genes among the upregulated genes in Fibro6 and Fibroblast-Cilp, whereas males showed upregulation of genes involved in protein folding or responses to endoplasmic reticulum stress (Figure 7E). We also observed upregulation of ECM-organizing genes in male Fibro6 cells that were distinct from those upregulated in female Fibro6 or Fibroblast-Cilp populations, underscoring a potential difference in the characteristics of ECM generated in the 2 sexes on tissue stress. GO enrichment analysis revealed no coherent biological themes shared among downregulated genes, principally in fibroblasts, whose was sexually dimorphic. To determine whether Fibroblast-Cilp and Fibroblast-Thbs4 populations are also present in the context of human pathological remodeling, we analyzed bulk cardiac RNA-seq data from 85 women and 172 men between the ages of 40 and 69 profiled by the Genotype-Tissue Project. of the NPPB transcript was examined to segregate humans without and with putative hypertrophy; this gene codes for the cardiac hypertrophy biomarker B-type natriuretic peptide. High NPPB was used as a rough proxy for potential hypertrophy, and samples with the highest (top 20%) NPPB (without regard to age or sex) were classified as putatively hypertrophic (Figure XXVIA in the Data Supplement). of POSTN, TGFB1, and collagens was significantly elevated in all hypertrophy groups in comparison with corresponding no-hypertrophy controls (uncorrected P<1e-07 for all genes; Figure XXVIB in the Data Supplement) confirming that human hypertrophy cohorts also exhibit other markers of cardiac stress. of THBS4, DDAH1, FMOD, CTHRC1, and CILP, 5 genes that are enriched in Fibroblast-Cilp and Fibroblast-Thbs4, was significantly increased in hypertrophic human hearts (1.4e-09 < P< 2e-03 for all genes, uncorrected; Figure XXVID in the Data Supplement). The overall concordance of these observations is striking because of the limitations of this underpowered statistical analysis, including imperfect status of NPPB transcript as a biomarker for heart failure, the inconsistency of human heart sample isolation, and age-dependent changes in the systemic milieu such as the onset of menopause. Nevertheless, our observations suggest that cells corresponding to Fibroblast-Cilp and Fibroblast-Thbs4 may be present in the context of human cardiac stress. DiscussionAlthough a great deal of data exist concerning individual cellular components that contribute to cardiac dysfunction, less is known about how these components interact during concerted pathological responses to tissue stress. To clarify the roles of distinct heart cell populations in the context of chronic maladaptive cardiac stress responses leading to fibrosis, we developed a strategy enabling unbiased analysis of myocyte and nonmyocyte transcriptomes in parallel using scRNA-seq. This generated a high-resolution atlas of the cardiac cellulome with and without fibrosis. Virtually all cell types in the heart have detectable transcriptional responses to tissue stress, with most cells participating in remodeling of the ECM. Specifically, 2 activated fibroblast populations emerge in response to stress, presumably to remodel the extracellular environment by direct ECM deposition and paracrine regulation of the heart cell niche. The extensive sexual dimorphism of cell-specific gene in resting and maladaptive hearts implicates sex-specific mechanisms driving cardiac physiology and stress responses. Overall, this study provides insight into the cellular and molecular drivers of cardiac pathological remodeling and fibrosis. As a resource for researchers, the extensive data set presented here comprises a deep and comprehensive map of the cardiac cellulome. These data were derived from 8 individual scRNA-seq libraries, with each library composed of female and male cells from a single animal of each sex. Cells from the 16 individual biological specimens can be isolated in silico to moderate effects of outliers and to examine the effect of biological sex on cardiac homeostasis and disease, which is an underexamined area in cardiac cell biology and physiology. 2,20 Our analysis attained a detailed portrait of the intercellular circuitry supporting the heart cell network at homeostasis and in response to the induction of pathological hypertrophy. For example, this data set implicates epicardial cells as critically important for supporting macrophage vitality and trafficking in the adult heart. Fibroblast-derived Csf1 and mural cell–derived Il34 are both growth factors that support macrophages that we previously detected in the heart4; however, the current analysis found that epicardial cells are the most macrophage-trophic cardiac cell population as a result of their high production of Csf1. Moreover, we identified distinct spatial distributions of macrophage subsets at the epicardium, suggesting additional unidentified regulation of macrophage positioning and trafficking. These results further support the important role of epicardial cells in macrophage homeostasis, in accordance with previous studies demonstrating the role of the epicardium in macrophage colonization of the heart.17 Also notable was the identification of CMs as the principal producers of Vegfa, which is a dominant endothelial and angiogenesis growth factor. Although we have not examined the precise isoforms of Vegfa expressed by cardiac cells, our analysis indicates that the cell type perhaps most dependent on tissue perfusion is the key driver of vascular innervation of the heart. Examination of ligand genes upregulated after induction of hypertrophy suggests the cardiac cellular niche is altered to support ECM remodeling and angiogenesis. This corresponded with upregulation of receptor genes mediating cell adhesion in a wide variety of cells, whereas angiogenesis-related receptors were upregulated in endothelial cells. GO enrichment analysis of all genes upregulated after induction of fibrosis indicated that a primary response of most cell types is the activation of genes involved in ECM remodeling. However, Fibroblast-Cilp and Fibroblast-Thbs4 appear to be the most fibrogenic cell populations based on the of transcripts related to ECM remodeling GO terms. Their fibrogenic phenotype is further accentuated in comparison with other fibroblast populations from AngII-treated mice that extensively activate ECM remodeling and collagen genes. Fibroblast-Cilp and Fibroblast-Thbs4 do not constitute classically defined myofibroblast populations, widely accepted as the main protagonists driving fibrosis.21 Recent scRNA-seq studies examining myocardial infarction have readily detected ACTA2+ myofibroblasts 5,6,22,23. The experimental conditions of these reports represent examples of tissue trauma where myofibroblast activation may mediate wound healing. In contrast, Acta2 transcript was not detected in Fibroblast-Cilp or in Fibroblast-Thbs4 above background transcriptomic levels, or at the protein level by histology. A lack of smooth muscle α-actin would be consistent with the absence of a need for structural support in a necrotic area. Fibrosis without myofibroblast conversion has been reported in a model of diabetic cardiomyopathy,24 and this fibrosis data set underscores that accumulation of extensive fibrosis in the heart is not dependent on ACTA2+ myofibroblasts. Fibroblast-Cilp and Fibroblast-Thbs4 exhibit a closer resemblance to matrifibrocytes that reside in a mature scar after myocardial infarction.14 Although of key matrifibrocyte marker genes, Comp, Sfrp2, or Wisp2, was not highly enriched in cells defined as Fibroblast-Cilp and Fibroblast-Thbs4 in this study, the functional inferences made of these cell populations are comparable. Therefore, cells analogous to Fibroblast-Cilp and Fibroblast-Thbs4 may promote fibrosis in multiple cardiac stress contexts. Extending on these findings, we examined of transcriptional signatures resembling Fibroblast-Cilp and Fibroblast-Thbs4 in the human context. Specifically, using human heart bulk RNA-seq data from normal and putatively hypertrophic cardiac samples, we found that genes with restricted to Fibroblast-Cilp and Fibroblast-Thbs4 appear to be upregulated in hypertrophy and correlate with the cardiac stress marker NPPB. This finding suggests that future studies targeting Fibroblast-Cilp and Fibroblast-Thbs4, or genes specifically activated within these cell subsets, may yield therapeutically relevant outcomes for addressing human disease. Analysis of downregulated genes indicates that a key feature of fibrosis induction is dampening of inflammation and leukocyte infiltration. Downregulated genes included CC chemokine ligands, which are key monocyte chemoattractants. Although AngII infusion caused an increase in macrophage abundance, observed by flow cytometry, this was not accompanied by significant increases in proinflammatory mediators. This suggests that the cardiac cellulome functions to dampen inflammation to regulate inflammatory responses. Furthermore, we found that mitochondrial bioenergetics are subdued in many cell types, in particular, in endothelial cells and macrophages, in accordance with the association between mitochondrial activity and heart failure. The current study also identified numerous cell type–specific genes that were robustly upregulated in hypertrophy (such as Rapsn in smooth muscle cells or Fastkd1 in epicardial cells) with unknown implications for cardiac pathology. Similar relationships were observed among downregulated genes (eg, Thap8 in macrophages or Has1 in the epicardium). These genes constitute appealing subjects for further research to explore their in further detail and their validity as potential therapeutic targets. Analysis of the abundance of Fibroblast-Cilp and Fibroblast-Thbs4 cell populations suggests sex differences in cellular abundance (Figure 7C), which is of particular interest in light of known sexual dimorphisms in cardiac stiffness.25 Indeed, analysis of cardiac fibrosis in AngII-treated mice revealed differences in fibrotic area in addition to the distribution of fibrotic loci. Females were characterized by less fibrosis in the left ventricle wall and the absence of fibrosis in the right ventricle walls, and less perivascular fibrosis, as well. In contrast, greater fibrosis accumulated in the female interventricular septum than in males, although a functional decline was only observed in males. In both the unstressed and stressed contexts, sexual dimorphisms were observed in cardiac gene across almost all cell types, with particularly robust differences in fibroblasts that are the major sex hormone receptor–expressing cells in the heart. In the context of cardiac fibrosis, female and male cardiac fibroblasts regulated distinct sets of genes that are associated with ECM remodeling, suggesting that cardiac ECM quality is biologically distinct in the 2 sexes. A key sexually dimorphic characteristic is the upregulation of genes related to unfolded protein responses (UPR) and endoplasmic reticulum stress in multiple male cell populations. UPRs may explain, at least in part, 2 sexually dimorphic features observed in this fibrosis data set. First, an increase in UPR may contribute to the increased fibrosis observed in male hearts (Figure 7A). Endoplasmic reticulum stress is observed in multiple models of tissue fibrosis including models of maladaptive cardiac remodeling,26 and attenuation of UPR inhibits fibrosis.27–29 Second, UPR may contribute to downregulation of translation-related genes in males. UPR may elicit this by multiple means including IRE1-mediated translational control of the transcription factor Xbp1, and regulation of multiple genes including eukaryotic initiation factor (Eif) genes and others involved in protein translation.26 Incidentally, Xbp1 is upregulated in multiple male cardiac cell types and is associated with tissue fibrosis in the heart and elsewhere.30,31 In addition to the link between protein translation and UPR, translation is also supported by estrogens in several species and in both neoplastic and nonneoplastic tissues.32–34 In tumor cells, inhibition of estrogen receptor α induces UPR.32 Moreover, estrogen has been recently shown to regulate Eif genes to promote protein synthesis.35 Many Eif genes are upregulated after AngII infusion across multiple female cell populations in addition to many ribosomal protein genes. Together, these observations suggest that female cells are better primed for protein synthesis, and avoidance of tissue fibrosis and maladaptive remodeling, as well. Conversely, the enrichment of genes implicated in angiogenesis was observed in multiple male cardiac cell populations including Fibro 2 and Fibro 6 in unstressed mice. Accordingly, increased numbers of endothelial cells were found in male hearts than in female hearts. In addition, significantly higher Vegfa was detected in male CMs and Fibroblast-Thbs4 cells after fibrosis induction. Although estrogen has been reported to enhance angiogenesis in a range of contexts,36 more recent research suggests that loss of estrogen receptor promotes angiogenesis.37 Moreover, testosterone induces angiogenesis in a sex-dependent manner, further underscoring the importance of the hormonal milieu.38 Last, UPR is also strongly linked to upregulation of proangiogenic genes,39 resulting in a more proangiogenic environment in male hearts. Although this study provides a framework for understanding cardiac cellular responses to a physiological stress, further work is required to demonstrate a causal relationship between specific cell types and genes and fibrosis. Moreover, limitations common to droplet-based single-cell sequencing studies, including detection of a limited number of genes and cells for each cell population, cell stress resulting from tissue dissociation, and the inability to analyze rare and difficult to isolate cell types, may have impeded the detection of important pathways for disease development. These limitations are particularly pertinent for our analysis of cardiomyocytes, where only nuclei, and, therefore, only a fraction of cardiomyocyte RNA, were analyzed. We also noted a significant correlation between cell population size and number of differentially regulated genes identified by comparing groups corresponding to AngII-treated and nontreated, and female and male. This suggests that our capacity to detect AngII-responsive and sexually dimorphic genes for smaller cell populations is limited; nevertheless, abundant cell populations are of particular interest in this initial study given their capacity to contribute to stress-induced cardiac remodeling. Last, although female and male mice were housed separately, the estrus cycles of female mice were not analyzed. This may affect the sexually dimorphic genes identified in this study. Nevertheless, the findings of this study offer unique insights into the multifaceted mechanisms driving maladaptive remodeling of the heart, including important components of this process that are sexually dimorphic. They also highlight the utility of the single-cell transcriptomics data and the overall approach presented in this study for exploring the cardiac cellulome and its changes in the context of chronic tissue stress. Further analysis of cardiac cell networks, with integration of data from different developmental and disease states, time points, and genetic backgrounds, will advance our understanding of cardiac cellularity in health and disease. AcknowledgmentsWe acknowledge the use of JAX Flow Cytometry Core, Microscopy Core, Surgical Services, Histology Service, and Single-Cell Biology Laboratory, and the Monash Micro Imaging facility at the Alfred Research Alliance for provision of microscopy instrumentation and training. The data used for the analyses described in this article were obtained from the Genotype-Tissue Portal on February 21, 2019. Sources of FundingThis work is supported by National Health and Medical Research Council (Australia) Ideas Grant (GNT1188503) to Dr Pinto. JAX Cores are supported by the Jackson Laboratory Cancer Center Core grant and the Leducq Foundation Transatlantic Network of Excellence in Cardiac Research to Dr Rosenthal. The Genotype-Tissue (GTEx) Project was supported by the Common Fund of the Office of the Director of the National Institutes of Health, and by National Cancer Institute, National Human Genome Research Institute, National Heart, Lung, and Blood Institute, National Institute on Drug Abuse, National Institute of Mental Health, and National Institute of Neurological Disorders and Stroke. DisclosuresNone. Supplemental MaterialsExpanded Methods Data Supplement Figures I–XXVII Data Supplement Tables I–XII Data Supplement Movie I Data Supplement Excel Files I–XIII References 40–50 Footnotes*M.A. McLellan and Dr Skelly contributed equally. Sources of Funding, see page 1462 The Data Supplement is available with this article at https://www./doi/suppl/10.1161/CIRCULATIONAHA.119.045115. Alexander R. Pinto, BS, PhD, Baker Heart and Diabetes Institute, 75 Commercial Rd, Prahran, Victoria, Australia 3004. Email alex. References
|
|
来自: 昵称37064826 > 《待分类》