Single cell RNA sequencing to dissect transcriptional signature of cardiac fibroblasts activation in experimental autoimmune myocarditis
Authors: Stefanska M , Tkacz K , Bielawska M , Weglarczyk K , Szatanek R , Siedlar M , Kania G , Blyszczuk P 
Published: November 2023

Abstract

Introduction
Dilated cardiomyopathy (DCM) is a condition that evolves from inflammation of the myocardium, and it is often associated with a poor prognosis and fatal outcome. Experimental autoimmune myocarditis (EAM) represents a CD4+ T cell-dependent animal model of acute myocarditis followed by the development of post-inflammatory DCM.

Purpose
The aim of the study was to investigate the transcriptional changes in cardiac fibroblasts during disease progression in a mouse model of EAM at the single cell level.

Methods
6–8-week-old BALB/c mice (or reporter mice expressing EGFP under the collagen type I promoter) were immunized with αMyHC peptide together with Complete Freunds’ Adjuvant. Live cardiac cells were sorted from healthy hearts (d0), at acute myocarditis (d19), during resolution of inflammation (d25), and at the DCM stage (d35) using the BD FACSAria II Cell Sorter. Single-cell cDNA libraries were prepared using 10x Genomics technology, and the sequencing was performed on an Illumina NovaSeq 6000 platform, targeting 50000 reads per cell. Single cell RNA sequencing data was analyzed using RStudio software and Seurat packages. At d0 and d19 of EAM, bulk RNAseq was performed on live cardiac EGFP+CD45-CD34- fibroblasts.

Results
Following a quality control check, we included in the analysis 25058 cells. We identified 13 distinct main cell populations, including fibroblasts, lymphoid and myeloid immune cells, endothelial cells, smooth muscle cells, and cardiomyocytes (Fig. 1A). Differential expression analysis revealed that fibroblasts in acute myocarditis (d19) showed activation of genes related to immune processes such as cytokine-cytokine receptor interaction, IL-17, and TNF-alpha signaling pathways. Subset analysis on 4901 cardiac fibroblasts resulted in 5 distinct fibroblast subsets, FB1–FB5 (Fig. 1B). A major fibroblast subset found in the heart (FB1) expanded during the course of EAM (Fig. 1C). Of note, cells expressing the myofibroblastic signature (FB4) were identified at d19 only (Fig. 1C). Fibroblast subset FB1 activated expression of proinflammatory genes, such as Ly6a, Ly6c1, Cxcl14, Ccl2, and Ccl7, particularly during the inflammatory phase of EAM (d19, Fig. 1D). Surprisingly, the expression of profibrotic genes, including Col1a1, Col3a1, and Postn in FB1, was downregulated during the progression of EAM (Fig. 1E). Bulk RNA sequencing analysis confirmed that cardiac fibroblasts in acute myocarditis (d19) activated immune processes (chemokine production), but also showed dysregulated expression of genes involved in cytoskeletal re-organization and extracellular matrix turnover.

Conclusions
Our data indicated that acute inflammation in the heart activated a pro-inflammatory transcriptional response in cardiac fibroblasts and induced a fibroblast-to-myofibroblast transition. Understanding changes in transcriptomic signatures on a single cell level can help to develop new therapies in inflammatory heart diseases.