Natural Hydrogels Support Kidney Organoid Generation and Promote in vitro Angiogenesis.

To date strategies aiming to modulate cell to extracellular matrix (ECM) interactions during organoid derivation remain largely unexplored. Here renal decellularized extracellular matrix (dECM) hydrogels are fabricated from porcine and human renal cortex as biomaterials to enrich cell-to-ECM crosstalk during the onset of kidney organoid differentiation from human pluripotent stem cells (hPSCs). Renal dECM-derived hydrogels are used in combination with hPSC-derived renal progenitor cells to define new approaches for 2D and 3D kidney organoid differentiation, demonstrating that in the presence of these biomaterials the resulting kidney organoids exhibit renal differentiation features, and the formation of an endogenous vascular component. Based on these observations, a new method to produce kidney organoids with vascular-like structures is achieved through the assembly of hPSC-derived endothelial-like organoids with kidney organoids in 3D. Major readouts of kidney differentiation and renal cell morphology are assessed exploiting these culture platforms as new models of nephrogenesis. Overall, this work shows that exploiting cell-to-ECM interactions during the onset of kidney differentiation from hPSCs facilitates and optimizes current approaches for kidney organoid derivation thereby increasing the utility of these unique culture cell platforms for personalized medicine. This article is protected by copyright. All rights reserved.

An atlas of cells in the human tonsil.

Palatine tonsils are secondary lymphoid organs (SLOs) representing the first line of immunological defense against inhaled or ingested pathogens. We generated an atlas of the human tonsil composed of >556,000 cells profiled across five different data modalities, including single-cell transcriptome, epigenome, proteome, and immune repertoire sequencing, as well as spatial transcriptomics. This census identified 121 cell types and states, defined developmental trajectories, and enabled an understanding of the functional units of the tonsil. Exemplarily, we stratified myeloid slan-like subtypes, established a BCL6 enhancer as locally active in follicle-associated T and B cells, and identified SIX5 as putative transcriptional regulator of plasma cell maturation. Analyses of a validation cohort confirmed the presence, annotation, and markers of tonsillar cell types and provided evidence of age-related compositional shifts. We demonstrate the value of this resource by annotating cells from B cell-derived mantle cell lymphomas, linking transcriptional heterogeneity to normal B cell differentiation states of the human tonsil.

Gene therapy restores the transcriptional program of hematopoietic stem cells in Fanconi anemia.

Clinical trials have shown that lentiviral-mediated gene therapy can ameliorate bone marrow failure (BMF) in nonconditioned Fanconi anemia (FA) patients resulting from the proliferative advantage of corrected FA hematopoietic stem and progenitor cells (HSPC). However, it is not yet known if gene therapy can revert affected molecular pathways in diseased HSPC. Single-cell RNA sequencing was performed in chimeric populations of corrected and uncorrected HSPC co-existing in the BM of gene therapy-treated FA patients. Our study demonstrates that gene therapy reverts the transcriptional signature of FA HSPC, which then resemble the transcriptional program of healthy donor HSPC. This includes a down-regulated expression of TGF-ß and p21, typically up-regulated in FA HSPC, and upregulation of DNA damage response and telomere maintenance pathways. Our results show for the first time the potential of gene therapy to rescue defects in the HSPC transcriptional program from patients with inherited diseases; in this case, in FA characterized by BMF and cancer predisposition.

The transcriptional regulator Sin3A balances IL-17A and Foxp3 expression in primary CD4 T cells.

The Sin3 transcriptional regulator homolog A (Sin3A) is the core member of a multiprotein chromatin-modifying complex. Its inactivation at the CD4/CD8 double-negative stage halts further thymocyte development. Among various functions, Sin3A regulates STAT3 transcriptional activity, central to the differentiation of Th17 cells active in inflammatory disorders and opportunistic infections. To further investigate the consequences of conditional Sin3A inactivation in more mature precursors and post-thymic T cell, we have generated CD4-Cre and CD4-CreERT2 Sin3AF/F mice. Sin3A inactivation in vivo hinders both thymocyte development and peripheral T-cell survival. In vitro, in Th17 skewing conditions, Sin3A-deficient cells proliferate and acquire memory markers and yet fail to properly upregulate Il17a, Il23r, and Il22. Instead, IL-2+ and FOXP3+ are mostly enriched for, and their inhibition partially rescues IL-17A+ T cells. Notably, Sin3A deletion also causes an enrichment of genes implicated in the mTORC1 signaling pathway, overt STAT3 activation, and aberrant cytoplasmic RORγt accumulation. Thus, together our data unveil a previously unappreciated role for Sin3A in shaping critical signaling events central to the acquisition of immunoregulatory T-cell phenotypes.

Uncovering perturbations in human hematopoiesis associated with healthy aging and myeloid malignancies at single-cell resolution.

Early hematopoiesis is a continuous process in which hematopoietic stem and progenitor cells (HSPCs) gradually differentiate toward specific lineages. Aging and myeloid malignant transformation are characterized by changes in the composition and regulation of HSPCs. In this study, we used single-cell RNA sequencing (scRNA-seq) to characterize an enriched population of human HSPCs obtained from young and elderly healthy individuals.Based on their transcriptional profile, we identified changes in the proportions of progenitor compartments during aging, and differences in their functionality, as evidenced by gene set enrichment analysis. Trajectory inference revealed that altered gene expression dynamics accompanied cell differentiation, which could explain aging-associated changes in hematopoiesis. Next, we focused on key regulators of transcription by constructing gene regulatory networks (GRNs) and detected regulons that were specifically active in elderly individuals. Using previous findings in healthy cells as a reference, we analyzed scRNA-seq data obtained from patients with myelodysplastic syndrome (MDS) and detected specific alterations of the expression dynamics of genes involved in erythroid differentiation in all patients with MDS such as TRIB2. In addition, the comparison between transcriptional programs and GRNs regulating normal HSPCs and MDS HSPCs allowed identification of regulons that were specifically active in MDS cases such as SMAD1, HOXA6, POU2F2, and RUNX1 suggesting a role of these transcription factors (TFs) in the pathogenesis of the disease.In summary, we demonstrate that the combination of single-cell technologies with computational analysis tools enable the study of a variety of cellular mechanisms involved in complex biological systems such as early hematopoiesis and can be used to dissect perturbed differentiation trajectories associated with perturbations such as aging and malignant transformation. Furthermore, the identification of abnormal regulatory mechanisms associated with myeloid malignancies could be exploited for personalized therapeutic approaches in individual patients.

Our blood contains many different types of cells; red blood cells carry oxygen through the body, platelets help to stop bleeding and a variety of white blood cells fight infections. All of these critical components come from a pool of immature cells in bone marrow, which can develop and specialise into any of these. However, as we get older, these immature cells can accumulate damage, including mutations in specific genes. This increases the risk of diseases such as myelodysplastic syndromes (MDS), a type of cancer in which the cells cannot develop and the patient does not have enough healthy mature blood cells. The changes in gene activity in the immature cells have previously been studied using samples from young and elderly people, as well as individuals with MDS. These studies examined large numbers of cells together, revealing differences between young and elderly people, and individuals with MDS. However, this does not describe how the different types alter their behaviour. To address this, Ainciburu, Ezponda et al. used a technique called single-cell RNA sequencing to study the gene activity in individual immature blood cells. This revealed changes associated with maturation that may account for the different combinations of cell populations in younger and older people. The results confirmed findings from previous studies and suggested new genes involved in ageing or MDS. Ainciburu, Ezponda et al. used these results to create an analytical system that highlights gene activity differences in individual MDS patients that are independent of age-related changes. These results provide new insights that could help further research into the development of MDS and the ageing process. In addition, scientists could study other diseases using this approach of analysing individual patients' gene activity. In future, this could help to personalise clinical decisions on diagnosis and treatment.

Revealing cell populations catching the early stages of human embryo development in naive pluripotent stem cell cultures.

Naive human pluripotent stem cells (hPSCs) are defined as the in vitro counterpart of the human preimplantation embryo's epiblast and are used as a model system to study developmental processes. In this study, we report the discovery and characterization of distinct cell populations coexisting with epiblast-like cells in 5iLAF naive human induced PSC (hiPSC) cultures. It is noteworthy that these populations closely resemble different cell types of the human embryo at early developmental stages. While epiblast-like cells represent the main cell population, interestingly we detect a cell population with gene and transposable element expression profile closely resembling the totipotent eight-cell (8C)-stage human embryo, and three cell populations analogous to trophectoderm cells at different stages of their maturation process: transition, early, and mature stages. Moreover, we reveal the presence of cells resembling primitive endoderm. Thus, 5iLAF naive hiPSC cultures provide an excellent opportunity to model the earliest events of human embryogenesis, from the 8C stage to the peri-implantation period.

The transcription factor DDIT3 is a potential driver of dyserythropoiesis in myelodysplastic syndromes.

Myelodysplastic syndromes (MDS) are hematopoietic stem cell (HSC) malignancies characterized by ineffective hematopoiesis, with increased incidence in older individuals. Here we analyze the transcriptome of human HSCs purified from young and older healthy adults, as well as MDS patients, identifying transcriptional alterations following different patterns of expression. While aging-associated lesions seem to predispose HSCs to myeloid transformation, disease-specific alterations may trigger MDS development. Among MDS-specific lesions, we detect the upregulation of the transcription factor DNA Damage Inducible Transcript 3 (DDIT3). Overexpression of DDIT3 in human healthy HSCs induces an MDS-like transcriptional state, and dyserythropoiesis, an effect associated with a failure in the activation of transcriptional programs required for normal erythroid differentiation. Moreover, DDIT3 knockdown in CD34+ cells from MDS patients with anemia is able to restore erythropoiesis. These results identify DDIT3 as a driver of dyserythropoiesis, and a potential therapeutic target to restore the inefficient erythroid differentiation characterizing MDS patients.

CAR density influences antitumoral efficacy of BCMA CAR T cells and correlates with clinical outcome.

Identification of new markers associated with long-term efficacy in patients treated with CAR T cells is a current medical need, particularly in diseases such as multiple myeloma. In this study, we address the impact of CAR density on the functionality of BCMA CAR T cells. Functional and transcriptional studies demonstrate that CAR T cells with high expression of the CAR construct show an increased tonic signaling with up-regulation of exhaustion markers and increased in vitro cytotoxicity but a decrease in in vivo BM infiltration. Characterization of gene regulatory networks using scRNA-seq identified regulons associated to activation and exhaustion up-regulated in CARHigh T cells, providing mechanistic insights behind differential functionality of these cells. Last, we demonstrate that patients treated with CAR T cell products enriched in CARHigh T cells show a significantly worse clinical response in several hematological malignancies. In summary, our work demonstrates that CAR density plays an important role in CAR T activity with notable impact on clinical response.

One-Step In Vitro Generation of ETV2-Null Pig Embryos.

Each year, tens of thousands of people worldwide die of end-stage organ failure due to the limited availability of organs for use in transplantation. To meet this clinical demand, one of the last frontiers of regenerative medicine is the generation of humanized organs in pigs from pluripotent stem cells (PSCs) via blastocyst complementation. For this, organ-disabled pig models are needed. As endothelial cells (ECs) play a critical role in xenotransplantation rejection in every organ, we aimed to produce hematoendothelial-disabled pig embryos targeting the master transcription factor ETV2 via CRISPR-Cas9-mediated genome modification. In this study, we designed five different guide RNAs (gRNAs) against the DNA-binding domain of the porcine ETV2 gene, which were tested on porcine fibroblasts in vitro. Four out of five guides showed cleavage capacity and, subsequently, these four guides were microinjected individually as ribonucleoprotein complexes (RNPs) into one-cell-stage porcine embryos. Next, we combined the two gRNAs that showed the highest targeting efficiency and microinjected them at higher concentrations. Under these conditions, we significantly improved the rate of biallelic mutation. Hence, here, we describe an efficient one-step method for the generation of hematoendothelial-disabled pig embryos via CRISPR-Cas9 microinjection in zygotes. This model could be used in experimentation related to the in vivo generation of humanized organs.

A diabetic milieu increases ACE2 expression and cellular susceptibility to SARS-CoV-2 infections in human kidney organoids and patient cells.

It is not well understood why diabetic individuals are more prone to develop severe COVID-19. To this, we here established a human kidney organoid model promoting early hallmarks of diabetic kidney disease development. Upon SARS-CoV-2 infection, diabetic-like kidney organoids exhibited higher viral loads compared with their control counterparts. Genetic deletion of the angiotensin-converting enzyme 2 (ACE2) in kidney organoids under control or diabetic-like conditions prevented viral detection. Moreover, cells isolated from kidney biopsies from diabetic patients exhibited altered mitochondrial respiration and enhanced glycolysis, resulting in higher SARS-CoV-2 infections compared with non-diabetic cells. Conversely, the exposure of patient cells to dichloroacetate (DCA), an inhibitor of aerobic glycolysis, resulted in reduced SARS-CoV-2 infections. Our results provide insights into the identification of diabetic-induced metabolic programming in the kidney as a critical event increasing SARS-CoV-2 infection susceptibility, opening the door to the identification of new interventions in COVID-19 pathogenesis targeting energy metabolism.

The bone marrow niche regulates redox and energy balance in MLL::AF9 leukemia stem cells.

Eradicating leukemia requires a deep understanding of the interaction between leukemic cells and their protective microenvironment. The CXCL12/CXCR4 axis has been postulated as a critical pathway dictating leukemia stem cell (LSC) chemoresistance in AML due to its role in controlling cellular egress from the marrow. Nevertheless, the cellular source of CXCL12 in the acute myeloid leukemia (AML) microenvironment and the mechanism by which CXCL12 exerts its protective role in vivo remain unresolved. Here, we show that CXCL12 produced by Prx1+ mesenchymal cells but not by mature osteolineage cells provide the necessary cues for the maintenance of LSCs in the marrow of an MLL::AF9-induced AML model. Prx1+ cells promote survival of LSCs by modulating energy metabolism and the REDOX balance in LSCs. Deletion of Cxcl12 leads to the accumulation of reactive oxygen species and DNA damage in LSCs, impairing their ability to perpetuate leukemia in transplantation experiments, a defect that can be attenuated by antioxidant therapy. Importantly, our data suggest that this phenomenon appears to be conserved in human patients. Hence, we have identified Prx1+ mesenchymal cells as an integral part of the complex niche-AML metabolic intertwining, pointing towards CXCL12/CXCR4 as a target to eradicate parenchymal LSCs in AML.

In vivo CRISPR-Cas9 inhibition of hepatic LDH as treatment of primary hyperoxaluria.

Genome-editing strategies, especially CRISPR-Cas9 systems, have substantially increased the efficiency of innovative therapeutic approaches for monogenic diseases such as primary hyperoxalurias (PHs). We have previously demonstrated that inhibition of glycolate oxidase using CRISPR-Cas9 systems represents a promising therapeutic option for PH type I (PH1). Here, we extended our work evaluating the efficacy of liver-specific inhibition of lactate dehydrogenase (LDH), a key enzyme responsible for converting glyoxylate to oxalate; this strategy would not be limited to PH1, being applicable to other PH subtypes. In this work, we demonstrate a liver-specific inhibition of LDH that resulted in a drastic reduction of LDH levels in the liver of PH1 and PH3 mice after a single-dose delivery of AAV8 vectors expressing the CRISPR-Cas9 system, resulting in reduced urine oxalate levels and kidney damage without signs of toxicity. Deep sequencing analysis revealed that this approach was safe and specific, with no off-targets detected in the liver of treated animals and no on-target/off-tissue events. Altogether, our data provide evidence that in vivo genome editing using CRISPR-Cas9 systems would represent a valuable tool for improved therapeutic approaches for PH.

Preneoplastic somatic mutations including MYD88L265P in lymphoplasmacytic lymphoma.

Normal cell counterparts of solid and myeloid tumors accumulate mutations years before disease onset; whether this occurs in B lymphocytes before lymphoma remains uncertain. We sequenced multiple stages of the B lineage in elderly individuals and patients with lymphoplasmacytic lymphoma, a singular disease for studying lymphomagenesis because of the high prevalence of mutated MYD88. We observed similar accumulation of random mutations in B lineages from both cohorts and unexpectedly found MYD88L265P in normal precursor and mature B lymphocytes from patients with lymphoma. We uncovered genetic and transcriptional pathways driving malignant transformation and leveraged these to model lymphoplasmacytic lymphoma in mice, based on mutated MYD88 in B cell precursors and BCL2 overexpression. Thus, MYD88L265P is a preneoplastic event, which challenges the current understanding of lymphomagenesis and may have implications for early detection of B cell lymphomas.

Tumor cells in light-chain amyloidosis and myeloma show distinct transcriptional rewiring of normal plasma cell development.

Although light-chain amyloidosis (AL) and multiple myeloma (MM) are characterized by tumor plasma cell (PC) expansion in bone marrow (BM), their clinical presentation differs. Previous attempts to identify unique pathogenic mechanisms behind such differences were unsuccessful, and no studies have investigated the differentiation stage of tumor PCs in patients with AL and MM. We sought to define a transcriptional atlas of normal PC development in secondary lymphoid organs (SLOs), peripheral blood (PB), and BM for comparison with the transcriptional programs (TPs) of tumor PCs in AL, MM, and monoclonal gammopathy of undetermined significance (MGUS). Based on bulk and single-cell RNA sequencing, we observed 13 TPs during transition of normal PCs throughout SLOs, PB, and BM. We further noted the following: CD39 outperforms CD19 to discriminate newborn from long-lived BM-PCs; tumor PCs expressed the most advantageous TPs of normal PC differentiation; AL shares greater similarity to SLO-PCs whereas MM is transcriptionally closer to PB-PCs and newborn BM-PCs; patients with AL and MM enriched in immature TPs had inferior survival; and protein N-linked glycosylation-related TPs are upregulated in AL. Collectively, we provide a novel resource to understand normal PC development and the transcriptional reorganization of AL and other monoclonal gammopathies.

Immunological Biomarkers of Fatal COVID-19: A Study of 868 Patients.

Information on the immunopathobiology of coronavirus disease 2019 (COVID-19) is rapidly increasing; however, there remains a need to identify immune features predictive of fatal outcome. This large-scale study characterized immune responses to severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) infection using multidimensional flow cytometry, with the aim of identifying high-risk immune biomarkers. Holistic and unbiased analyses of 17 immune cell-types were conducted on 1,075 peripheral blood samples obtained from 868 COVID-19 patients and on samples from 24 patients presenting with non-SARS-CoV-2 infections and 36 healthy donors. Immune profiles of COVID-19 patients were significantly different from those of age-matched healthy donors but generally similar to those of patients with non-SARS-CoV-2 infections. Unsupervised clustering analysis revealed three immunotypes during SARS-CoV-2 infection; immunotype 1 (14% of patients) was characterized by significantly lower percentages of all immune cell-types except neutrophils and circulating plasma cells, and was significantly associated with severe disease. Reduced B-cell percentage was most strongly associated with risk of death. On multivariate analysis incorporating age and comorbidities, B-cell and non-classical monocyte percentages were independent prognostic factors for survival in training (n=513) and validation (n=355) cohorts. Therefore, reduced percentages of B-cells and non-classical monocytes are high-risk immune biomarkers for risk-stratification of COVID-19 patients.

Design and Synthesis of Novel Epigenetic Inhibitors Targeting Histone Deacetylases, DNA Methyltransferase 1, and Lysine Methyltransferase G9a with In Vivo Efficacy in Multiple Myeloma.

Concomitant inhibition of key epigenetic pathways involved in silencing tumor suppressor genes has been recognized as a promising strategy for cancer therapy. Herein, we report a first-in-class series of quinoline-based analogues that simultaneously inhibit histone deacetylases (from a low nanomolar range) and DNA methyltransferase-1 (from a mid-nanomolar range, IC50 < 200 nM). Additionally, lysine methyltransferase G9a inhibitory activity is achieved (from a low nanomolar range) by introduction of a key lysine mimic group at the 7-position of the quinoline ring. The corresponding epigenetic functional cellular responses are observed: histone-3 acetylation, DNA hypomethylation, and decreased histone-3 methylation at lysine-9. These chemical probes, multitarget epigenetic inhibitors, were validated against the multiple myeloma cell line MM1.S, demonstrating promising in vitro activity of 12a (CM-444) with GI50 of 32 nM, an adequate therapeutic window (>1 log unit), and a suitable pharmacokinetic profile. In vivo, 12a achieved significant antitumor efficacy in a xenograft mouse model of human multiple myeloma.

Characterization of complete lncRNAs transcriptome reveals the functional and clinical impact of lncRNAs in multiple myeloma.

Multiple myeloma (MM) is an incurable disease, whose clinical heterogeneity makes its management challenging, highlighting the need for biological features to guide improved therapies. Deregulation of specific long non-coding RNAs (lncRNAs) has been shown in MM, nevertheless, the complete lncRNA transcriptome has not yet been elucidated. In this work, we identified 40,511 novel lncRNAs in MM samples. lncRNAs accounted for 82% of the MM transcriptome and were more heterogeneously expressed than coding genes. A total of 10,351 overexpressed and 9,535 downregulated lncRNAs were identified in MM patients when compared with normal bone-marrow plasma cells. Transcriptional dynamics study of lncRNAs in the context of normal B-cell maturation revealed 989 lncRNAs with exclusive expression in MM, among which 89 showed de novo epigenomic activation. Knockdown studies on one of these lncRNAs, SMILO (specific myeloma intergenic long non-coding RNA), resulted in reduced proliferation and induction of apoptosis of MM cells, and activation of the interferon pathway. We also showed that the expression of lncRNAs, together with clinical and genetic risk alterations, stratified MM patients into several progression-free survival and overall survival groups. In summary, our global analysis of the lncRNAs transcriptome reveals the presence of specific lncRNAs associated with the biological and clinical behavior of the disease.

Deep MRD profiling defines outcome and unveils different modes of treatment resistance in standard- and high-risk myeloma.

Patients with multiple myeloma (MM) carrying standard- or high-risk cytogenetic abnormalities (CAs) achieve similar complete response (CR) rates, but the later have inferior progression-free survival (PFS). This questions the legitimacy of CR as a treatment endpoint and represents a biological conundrum regarding the nature of tumor reservoirs that persist after therapy in high-risk MM. We used next-generation flow (NGF) cytometry to evaluate measurable residual disease (MRD) in MM patients with standard- vs high-risk CAs (n = 300 and 90, respectively) enrolled in the PETHEMA/GEM2012MENOS65 trial, and to identify mechanisms that determine MRD resistance in both patient subgroups (n = 40). The 36-month PFS rates were higher than 90% in patients with standard- or high-risk CAs achieving undetectable MRD. Persistent MRD resulted in a median PFS of âˆ¼3 and 2 years in patients with standard- and high-risk CAs, respectively. Further use of NGF to isolate MRD, followed by whole-exome sequencing of paired diagnostic and MRD tumor cells, revealed greater clonal selection in patients with standard-risk CAs, higher genomic instability with acquisition of new mutations in high-risk MM, and no unifying genetic event driving MRD resistance. Conversely, RNA sequencing of diagnostic and MRD tumor cells uncovered the selection of MRD clones with singular transcriptional programs and reactive oxygen species-mediated MRD resistance in high-risk MM. Our study supports undetectable MRD as a treatment endpoint for patients with MM who have high-risk CAs and proposes characterizing MRD clones to understand and overcome MRD resistance. This trial is registered at www.clinicaltrials.gov as #NCT01916252.

Immunologic characterization of COVID-19 patients with hematological cancer.

Not available.

Single-Cell RNA Sequencing Analysis Reveals a Crucial Role for CTHRC1 (Collagen Triple Helix Repeat Containing 1) Cardiac Fibroblasts After Myocardial Infarction.

BACKGROUND: Cardiac fibroblasts (CFs) have a central role in the ventricular remodeling process associated with different types of fibrosis. Recent studies have shown that fibroblasts do not respond homogeneously to heart injury. Because of the limited set of bona fide fibroblast markers, a proper characterization of fibroblast population heterogeneity in response to cardiac damage is lacking. The purpose of this study was to define CF heterogeneity during ventricular remodeling and the underlying mechanisms that regulate CF function. METHODS: Collagen1α1-GFP (green fluorescent protein)-positive CFs were characterized after myocardial infarction (MI) by single-cell and bulk RNA sequencing, assay for transposase-accessible chromatin sequencing, and functional assays. Swine and patient samples were studied using bulk RNA sequencing. RESULTS: We identified and characterized a unique CF subpopulation that emerges after MI in mice. These activated fibroblasts exhibit a clear profibrotic signature, express high levels of Cthrc1 (collagen triple helix repeat containing 1), and localize into the scar. Noncanonical transforming growth factor-ß signaling and different transcription factors including SOX9 are important regulators mediating their response to cardiac injury. Absence of CTHRC1 results in pronounced lethality attributable to ventricular rupture. A population of CFs with a similar transcriptome was identified in a swine model of MI and in heart tissue from patients with MI and dilated cardiomyopathy. CONCLUSIONS: We report CF heterogeneity and their dynamics during the course of MI and redefine the CFs that respond to cardiac injury and participate in myocardial remodeling. Our study identifies CTHRC1 as a novel regulator of the healing scar process and a target for future translational studies.