# IPSC derived cardiac fibroblasts of DMD patients show compromised actin microfilaments, metabolic shift and pro-fibrotic phenotype

Salwa Soussi, Lesia Savchenko, Davide Rovina, Jason Iacovoni, Andrea Gottinger, Maxime Vialettes, Josè-Manuel Pioner, Andrea Farini, Sara Mallia, Martina Rabino, et al.

► To cite this version:

Salwa Soussi, Lesia Savchenko, Davide Rovina, Jason Iacovoni, Andrea Gottinger, et al. IPSC derived cardiac fibroblasts of DMD patients show compromised actin microfilaments, metabolic shift and pro-fibrotic phenotype. Biology Direct, 2023, 18 (1), pp.41. 10.1186/s13062-023-00398-2. hal-04302399

HAL Id: hal-04302399

https://hal.science/hal-04302399v1

Submitted on 26 Nov 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons CC BY 4.0 - Attribution - International License
RESEARCH

Open Access

<figure><img src="image_2.png" /></figure>

Check for
updates

# IPSC derived cardiac fibroblasts of DMD patients show compromised actin microfilaments, metabolic shift and pro-fibrotic phenotype

Salwa Soussi¹,⁷, Lesia Savchenko¹,⁷, Davide Rovina², Jason S. Iacovoni¹,³, Andrea Gottinger¹, Maxime Vialettes⁷, Josè-Manuel Pioner⁵, Andrea Farini⁶, Sara Mallia², Martina Rabino², Giulio Pompilio²,⁴, Angelo Parini¹,⁷, Olivier Lairez¹,⁷, Aoife Gowran² and Nathalie Pizzinat¹,⁷*

**Abstract**

Duchenne muscular dystrophy (DMD) is a severe form of muscular dystrophy caused by mutations in the dystrophin gene. We characterized which isoforms of dystrophin were expressed by human induced pluripotent stem cell (hiPSC)-derived cardiac fibroblasts obtained from control and DMD patients. Distinct dystrophin isoforms were observed; however, highest molecular weight isoform was absent in DMD patients carrying exon deletions or mutations in the dystrophin gene. The loss of the full-length dystrophin isoform in hiPSC-derived cardiac fibroblasts from DMD patients resulted in deficient formation of actin microfilaments and a metabolic switch from mitochondrial oxidation to glycolysis. The DMD hiPSC-derived cardiac fibroblasts exhibited a dysregulated mitochondria network and reduced mitochondrial respiration, with enhanced compensatory glycolysis to sustain cellular ATP production. This metabolic remodeling was associated with an exacerbated myofibroblast phenotype and increased fibroblast activation in response to pro fibrotic challenges. As cardiac fibrosis is a critical pathological feature of the DMD heart, the myofibroblast phenotype induced by the absence of dystrophin may contribute to deterioration in cardiac function. Our study highlights the relationship between cytoskeletal dynamics, metabolism of the cell and myofibroblast differentiation and provides a new mechanism by which inactivation of dystrophin in non-cardiomyocyte cells may increase the severity of cardiopathy.

**Keywords** Duchenne, Human induced pluripotent stem cell, Fibroblasts, Mitochondrial oxidation, Actin microfilaments

<footnote>*Correspondence:
Nathalie Pizzinat
nathalie.pizzinat@inserm.fr
Full list of author information is available at the end of the article</footnote>
## Introduction

Duchenne muscular dystrophy (DMD) is the most common form of muscular dystrophy caused by mutations in the DMD gene on the X chromosome (Xp21.1-Xp22), [1]. The different mutations observed in Duchenne patients (including exon deletions and duplications, as well as point mutations) that disrupt the translational reading frame result in the lack of the full-length form of Dystrophin, a component of the dystrophin glycoprotein complex (DGC). Dystrophin links the cytoskeleton and the extracellular matrix to help maintain muscle membrane integrity, in part by dissipating the forces of muscle contraction into the extracellular matrix. Thus, deficiencies in dystrophin result in chronic muscle damage, dysregulation of repair and progressive skeletal muscle weakness [2]. Cardiac muscle is also affected by dystrophin loss and often leads to dilated cardiomyopathies [3, 4]. Typically, cardiomyopathy begins in the first decade of life and evolves with progressive replacement of the myocardium by fibrous and fatty connective tissue resulting in end-stage dilated cardiomyopathy, which is a major cause of death in DMD patients [5]. Several mechanical disturbances have been observed in the DMD heart among them; increased heart rate variability, Ca²⁺ homeostasis dysregulation, tissue remodeling with fibro-fatty accumulation and cardiomyocyte death [6]. The dying cardiac cells trigger an inflammatory response that elicits immune cell infiltration and promotes fibrotic scar tissue, which impairs the compliance of the heart and worsens cardiac dysfunction. There may also be additional mechanisms by which dystrophin mutation engenders ventricular fibrosis that have yet to be uncovered.

Skeletal and cardiac myocytes expressed mainly the full-length (427 kDa) dystrophin isoform of the *DMD* gene that comprises 79 exons. The full-length protein consists of at least four main regions; an N-terminal domain plays an essential role in binding actin (ABD1), a central rod domain that harbors a second actin binding domain (ABD2), a cysteine-rich domain, which mainly interacts with the transmembrane protein dystroglycan and therefore indirectly with other DGC membrane proteins, followed by the C-terminal region [7, 8]. In addition, five smaller dystrophin isoforms, named according to their molecular weight (Dp260, Dp140, Dp116, Dp71 and Dp40), can be generated by the differential use of internal promoters. Furthermore, alternative splicing of specific exons, exon junctions or inclusion of intronic sequences can give rise to splice variants with modified carboxy termini [2]. Depending on the domain composition of the various isoforms, they have partially retained the capacity to bind to different proteins of the dystroglycan complex [9]. Although the expression of Dp71 is ubiquitous

in cells, the existence of splicing variants has been particularly documented in retina and brain [10–12]. In the heart, dystrophin is abundantly expressed in myocytes, and a large body of work has been devoted to studying the consequences of dystrophin mutation, whereas non-muscle cells have been little studied until now. In this study, we revealed the expression of dystrophin isoforms in hiPSC-derived cardiac fibroblasts. The loss of dystrophin in DMD patients modified actin cytoskeleton, induced a metabolic remodeling, exacerbated the myofibroblast phenotype and worsened fibroblast activation in response to pro fibrotic challenges. Our study highlights the relationship between cytoskeletal dynamics, metabolism of the cell and myofibroblast differentiation. Our work contributes to the understanding of the pathogenic mechanisms that occur in dystrophic cardiopathy.

## Material and methods

### Cell culture and treatments

The isolation of the human cells and the subsequent reprogramming into iPSC lines was performed in conformation with the Declaration of Helsinki. iPS cells from 3 DMD patients and 3 healthy donors were obtained from Centro Cardiologico Monzino, Milan Italy, as previously described [13–15]. The three Duchenne patients carried different *DMD* mutations: two exons deletions and a point mutation. In particular, DMD1 had a deletion of exons 49 and 50 [13], while DMD4 carried a bigger deletion from exon 45 to exon 50 [15]. DMD5 carried a point mutation c.6913-1 G > A; this mutation causes the deletion of a nucleotide in the mature transcript with a shift in the reading frame that causes a premature stop codon.

### hiPSC-EPI differentiation

Cell reprogramming into iPSC lines was previously described by Rovina et al. [15]. Briefly, iPSCs were seeded in Matrigel (Corning, 356,230) coated 6 well plates in StemFlex medium (ThermoFisher, A3349401) supplemented with RevitaCell (ThermoFisher, A2644501). On Day 0, medium was replaced with BPEL medium and cells were treated with 6 mM CHIR99021 (Bertin bioreagent, 13,122), 100 µg/ml Activin A (MiltenyiBiotec, 130-115-008), and 100 µg/ml BMP4 (R&D Systems, 5020-BP). On Day 3, cells were treated with 10 mM XAV939 (R&D Systems, 3748), 100 µg/ml BMP4 and 2 mM Retinoic acid (Sigma Aldrich, R2625). Medium is changed on Day6 with BPEL + 30 ng/ml BMP4 + 1 µM retinoic acid. After 3 days, cells were replated on fibronectin coated plates with BPEL medium supplemented with 20 mM SB 431542 (R&D Systems, 161).
## hiPSC-fibs differentiation

EPI-iPSC cells were seeded in vitronectin (ThermoFisher, A14700) coated 6 well plates in BPEL medium supplemented with 20 µg/ml FGF-2 (Peprotech, AF-100-18B). Medium was changed every two days for 10 days. At Day10, medium is switched to FGM3 (Promocell, C-39350) and changed every other day. Cells were used between passages 0 and 3. Profibrotic stimulations were performed by incubating iPSC derived cardiac fibroblasts with 10 ng/ml TGF beta or 100 ng/ml angiotensin II for 24 h in FGM3 containing 0.5% SVF.

## RNA sequencing and bioinformatics

RNA was isolated using ReliaPrep RNA Cell Miniprep System (Promega). Samples quantity and quality were evaluated on a Fragment analyzer (Advanced analytical). Libraries were generated with TruSeq stranded mRNA Illumina and quantified on a Fragment analyzer by the GeT-PlaGe Genotoul Toulouse-France.

Sequencing was performed by Integragen (France) via an Illumina NovaSeq™6000 S2 high throughput sequencer resulting in 30 to 38 million paired end reads (2×100 bp) per sample. Reads were cleaned with HTStream (v1.3.2) [https://github.com/s4hts/HTStream] and aligned with Salmon (v1.9.0) ([16] using the gcBias option, against the Gencode v40 human transcript database [17]).

The rest of the analysis was performed in R (v4.2.3) [R Core Team (2021). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. [https://www.R-project.org/]]. Tximport (v1.26.1) [18] was used to load in the Salmon mapping results. Functions from edgeR (v3.40.2) [19] were used to filter genes by expression and calculate normalization factors. limma (v3.54.2) [20] was used for the voom With Quality Weights function and then

<table><caption>**Table 1** Oligonucleotides used for qPCR</caption><thead><tr><td>Target</td><td>Sequence</td></tr></thead><tbody><tr><td>36B4 F</td><td>CAGATCACGTCATCGCACAAC</td></tr><tr><td>36B4 R</td><td>AAAAGGAGGTCTTCTCGGGC</td></tr><tr><td>α-SMA F</td><td>GGAGCAGCCCAGCCAAGC</td></tr><tr><td>α-SMA R</td><td>AGAGCCCAGAGCCATTGTCAC</td></tr><tr><td>COL1A F</td><td>GAGGGCCAAGACGAAGACATC</td></tr><tr><td>COL1A R</td><td>CAGATCACGTCATCGCACAAC</td></tr><tr><td>GAPDH F</td><td>AAGGTCGGAGTCAACGGATTT</td></tr><tr><td>GAPDH R</td><td>ATGAAGGGTCATTGATGGCA</td></tr><tr><td>GLUT1 F</td><td>GAGGGCCAAGACGAAGACATC</td></tr><tr><td>GLUT1 R</td><td>CAGATCACGTCATCGCACAAC</td></tr><tr><td>TNC F</td><td>CAACCATCACTGCCAAGTTCACAA</td></tr><tr><td>TNC R</td><td>GGGGGTCGCCAGGTAAGGAG</td></tr></tbody></table>

to perform cyclic Loess normalization with the “affy” method. Then we employed sPLS-DA from the mixOmics (v6.22.0) to select genes that we able to select genes differential with respect to the DMD and CON group factor. Figures were created with mixOmics plotIndiv, a custom ggplot2 (v3.4.2) [H. Wickham. ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag New York, 2016] adaptation of mixOmics' plotVar and pheatmap (v1.0.12) [https://cran.r-project.org/package=pheatmap].

Pathway enrichment was performed with EnrichR [21] using the KEGG 2021 Human database [22] and Human-Cyc 2016 database [23].

RNA-Seq data analysis was done by bioinformatic platform in Toulouse, functional annotation clustering and for up-regulated and down-regulated genes with a fold change of ≥ 2 and enrichment score of ≥ 1.

## Reverse transcription and real time PCR

Total RNA was purified using ReliaPrep RNA Cell Miniprep System (Promega) according to the manufacturer's instructions. RNA was quantified with NanoDrop (ND 2000, Thermofisher). The same RNA was retro-transcribed with the Applied biosystems Synthesis Kit (ThermoFisher Scientific, 4,368,814). Real time qPCR reactions were performed using TB Green Premix Ex Taq II (Takara) on a QuantStudio5 Real-Time PCR System (ThermoFisher) instrument. Gene expression levels were normalized to GAPDH and 36B4 housekeeping genes. Results were analyzed using the ΔΔCT method.

For mitochondrial DNA quantification, genomic DNA was purified on QIAamp DNA micro isolation kit (QIAGEN) following the manufacturer's instructions. ND1 and COX1 genes expression was normalized to nuclear DNA B2M using the ΔCT method.

## Seahorse analysis

hiPSC-fibs were seeded into vitronectin coated Seahorse 24 assay plates at a density of 80 000cells/well. Oxygen consumption rate (OCR) and extracellular acidification rate (ECAR) were assessed using a standard mitochondrial and glycolysis stress on the Seahorse Bioscience XF-24 analyzer (Agilent technologies). On the day of metabolic flux analysis, cells were washed with Seahorse XF base medium pH7.4 (Agilent technologies). The culture medium was replaced with 500 µl of Seahorse medium supplemented with 10 mM glucose, 1 mM sodium pyruvate and 2 mM glutamine and incubated 1 h in a CO2-free incubator at 37 °C. For XF glycolysis stress test, cells were incubated in DMEM with 2 mM glutamine. For the XF mito stress, inhibitors of the mitochondrial electron transport chain (oligomycin 1 µM, carbonyl cyanide 4-trifluoromethoxy phenylhydrazone (FCCP) 4 µM, rotenone and antimycin 1 µM)
<table><caption>**Table 2** Antibodies used for western blot and immunofluorescence</caption><thead><tr><th>Target</th><th>Provider</th><th>Ref</th></tr></thead><tbody><tr><td>Collagen I</td><td>Cell Signaling Technology</td><td>#72026</td></tr><tr><td>Dystrophin (Mandra1)</td><td>Santa Cruz Biotechnology</td><td>sc-73592</td></tr><tr><td>Dystrophin (7A10)</td><td>Santa Cruz Biotechnology</td><td>sc-33697</td></tr><tr><td>Dystrophin</td><td>Abnova</td><td>630-170</td></tr><tr><td>Dystrophin</td><td>Thermofisher</td><td>PA5-16734</td></tr><tr><td>DRP1</td><td>Abcam</td><td>ab184247</td></tr><tr><td>HK1</td><td>Santa Cruz Biotechnology</td><td>sc-46695</td></tr><tr><td>HK2</td><td>Santa Cruz Biotechnology</td><td>sc-374091</td></tr><tr><td>HSP60</td><td>Cell Signaling Technology</td><td>#12165</td></tr><tr><td>LDHA</td><td>Santa Cruz Biotechnology</td><td>sc-137243</td></tr><tr><td>LDHB</td><td>Santa Cruz Biotechnology</td><td>sc-100775</td></tr><tr><td>MFN2</td><td>Abcam</td><td>ab56889</td></tr><tr><td>Oxphos</td><td>Invitrogen</td><td>45-8099</td></tr><tr><td>PDHE1a</td><td>Santa Cruz Biotechnology</td><td>sc-377092</td></tr><tr><td>Phospho PDHE1a</td><td>Cell Signaling Technology</td><td>#31866</td></tr><tr><td>PKLR</td><td>Santa Cruz Biotechnology</td><td>sc-133222</td></tr><tr><td>RhoGDi</td><td>Santa Cruz Biotechnology</td><td>sc-365190</td></tr><tr><td>VDAC1</td><td>Santa Cruz Biotechnology</td><td>sc-390996</td></tr><tr><td>α-SMA</td><td>Sigma</td><td>A5228</td></tr><tr><td>β-Tubulin</td><td>Cell Signaling Technology</td><td>#2146</td></tr><tr><td>γ-Actin</td><td>Cell Signaling Technology</td><td>sc-65638</td></tr><tr><td>β-Actin</td><td>Santa Cruz Biotechnology</td><td>sc-47778</td></tr></tbody></table>

were sequentially injected to assess the OCR and the respiratory parameters. For XF glycolysis stress, glucose 10 mM, oligomycin 1 µM and 2-deoxy-D-glucose (2-DG) were injected to evaluate the ECAR and the glycolysis parameters. OCR and ECAR were automatically calculated by the Seahorse XF-24 software.

## Cellular ATP levels

Intracellular ATP levels were measured using the Cell Titer-Glo 2.0 cell viability assay (Promega). Briefly, $4 \times 10^4$ cells were loaded into each well in a 96well plate. To evaluate ATP generated from the oxidative phosphorylation or glycolysis, cells were treated with sodium iodoacetate 100 µM, oligomycin 1.6 µM

<table><caption>**Table 3** Antibodies used for flow cytometry</caption><thead><tr><th>Fluorochrome</th><th>Target</th><th>Clone</th><th>Provider</th><th>Ref</th></tr></thead><tbody><tr><td>BV421</td><td>CD90</td><td>5E10</td><td>Biolegend</td><td>328121</td></tr><tr><td>PE</td><td>CD140a (PDGFRα)</td><td>16A1</td><td>Biolegend</td><td>323505</td></tr><tr><td>AF488</td><td>CD29</td><td>TS2/16</td><td>Biolegend</td><td>303015</td></tr><tr><td>APC</td><td>CD105</td><td>43A3</td><td>Biolegend</td><td>323208</td></tr><tr><td>APC Vio700</td><td>CD45</td><td>REA747</td><td>Miltenyi</td><td>130-110-635</td></tr><tr><td>PerCp cy5.5</td><td>CD31</td><td>WM59</td><td>Biolegend</td><td>303131</td></tr></tbody></table>

antimycin A 100 µM (OAA) or PBS and incubated 1 h at 37 °C. 100 µl of Cell Titer Glo was added to each well and luminescence was measured using Infinite F500 reader (Tecan).

## Intracellular lactate levels

Intracellular L-lactate was measured using Lactate colorimetric assay kit (Sigma-Aldrich) according to the manufacturer's instructions. $1.6 \times 10^6$ cells were loaded per well in a 96well plate and mixed with master reaction containing enzyme mix, probe and buffer. After 30 min incubation, the colorimetric intensity was measured at 570 nm on Infinite F200 Pro (Tecan) (Table 1).

## Western blot analysis

Whole cell lysates were prepared using RIPA lysis buffer (Cell Signaling, 9806) containing protease and phosphatase inhibitor cocktail (Thermofisher scientific,78,420). Concentration of protein lysates were calculated using the BCA protein assay (ThermoFisher, 23,225). 10-50 µg were loaded onto a 12% polyacrylamide gel and subjected to gel electrophoresis. The resolved proteins were transferred with the aid of a Trans-Blot Turbo Transfer System (Bio-Rad) according to the manufacturer's recommendation (1A, 25 V, 30 min). For dystrophin, resolved proteins were transferred into 0.2 µm nitrocellulose membrane using wet transfer overnight at 4 °C at a constant current of 90 mA. The PVDF membrane was blocked in Tris-Buffered Saline 0.1% Tween TBST

(See figure on next page.)

**Fig. 1** Characterization of iPSC-derived cardiac fibroblasts from DMD and control patients. **A:** Western blot analysis of dystrophin protein expression in human iPSC-derived cardiac fibroblasts. 80 µg of iPSC-derived cardiac fibroblasts and 10 µg of heart were loaded, transferred on membrane and revealed with monoclonal antibody against the C-terminal end of the protein. Three major bands at 427 kDa, 71 kDa and 40 kDa were observed in control patients whereas DMD patients expressed only the 40 kDa and faint amount of 71 kDa. **B:** Representative Flow cytometry graphs obtained with iPSC-derived cardiac fibroblasts generated from DMD and control patients. Cells were first identified on a forward scatter/ side scatter (FSC-A/SSC-A) dot plot, doublets were removed and live cells were selected. Cells expressed different levels of CD29, CD105, CD90 (thy1.1) and PDGFRα. Endothelial and leucocyte markers, CD31 and CD45 respectively, were negative (not shown). **C:** Flow cytometry quantification of expressed markers Data are presented as mean ± SEM
A

<figure><img src="image_3.png" /></figure>

B

<figure><img src="image_4.png" /></figure>

C

<figure><img src="image_5.png" /><figcaption>**Fig. 1** (See legend on previous page.)</figcaption></figure>
containing 5% Bovine Serum Albumin BSA (Sigma) and washed with TBST. The membranes were then incubated with the appropriate antibodies (Table 2). After an overnight incubation period at 4 °C, the membranes were washed in TBST and incubated with appropriate secondary antibody conjugated to horseradish peroxidase for 1h30 at room temperature. Proteins were detected using Clarity western ECL Substrate (Bio-Rad) with the aid of ChemiDoc MP Imaging System (Bio-Rad).

## Flow cytometry

Phenotypic analysis of hiPSC-fibs was performed using multicolor flow cytometry. Briefly, nonspecific binding of antibodies by Fc receptors was blocked by incubating the cells with mouse serum diluted in FACS buffer (PBS 4% FBS and 2 mM EDTA). Cell surface markers were then stained with antibodies conjugated with fluorochromes (Table 3). Viability was assessed by incubation with Live-Dead Yellow (Invitrogen). The following isotype controls mouse IgG1k BV421, mouse IgG1k PE, mouse IgG1k FITC and mouse IgG1k APC were used for compensation setups. Samples were acquired on LSR Fortessa (BD Bioscience). Analysis of flow cytometric data was performed using FlowJo (TreeStar).

## Immunofluorescence staining

Cultured iPSC-CF on vitronectin coated glasses were fixed in paraformaldehyde 4% and permeabilized with Triton 0.1% (Sigma). The cells were blocked for 1h30 and incubated during 1 h with primary antibodies (Table 3) at room temperature, then secondary antibody for 1h30. Subsequently, cells were incubated with Phalloidin (1:400, A12381; Thermo Fisher Scientific) to stain F-actin. After washing steps, cells were marked with DAPI and slides were mounted using fluoromount mounting medium (Invitrogen). Images were taken with a LSM900 Confocal laser scanning microscope using the ×63 oil immersion objective lens. Immunofluorescence quantifications were made on 5–6 images for each

condition. Phalloidin and collagen quantification was performed using the mean fluorescence intensity (MFI) on Fiji. Mitochondrial network analysis quantification was performed using MiNA (Mitochondrial Network Analysis) workflow on Fiji, analyzing the "networks" and "mitochondrial footprint".

## Statistics

Data were obtained from three independent DMD patients and healthy donors. Results are expressed as means ± standard errors of the mean (SEM). Statistical analysis was performed using GraphPad Prism version 9. Statistical difference was analyzed by the student's T test when comparing two groups and a one-way ANOVA followed by the Tukey posttest for multiple comparisons. Values of *P* < 0.05 were considered statistically significant.

## Results

### The 427 kDa of dystrophin isoform is expressed in healthy human iPSC-derived cardiac fibroblasts (control hiPSC-fibs) and absent in hiPSC-fibs isolated from DMD patients

We explored isoforms of dystrophin expressed by hiPSC-derived cardiac fibroblasts (hiPSC-fibs) obtained from control and DMD patients with out of frame deletions of exons (DMD1 and DMD4) and a patient who has a spontaneous point mutation in intron 47 (DMD5) that leads to a frameshift and the formation of a premature stop signal [13–15].

Western blot analysis revealed that a 427kDa form of dystrophin, recognized in control hiPSC-fibs was absent in all DMD hiPSC-fibs patients. The 71 kDa isoform was expressed control hiPSC-fibs whereas a faint signal was observed in DMD hiPSC-fibs. In contrast the 42 kDa isoform was present in both control and DMD hiPSC-fibs (Fig. 1A). Similar results were obtained with affinity purified polyclonal antibodies (Additional file 1: Fig. S1). To explore the consequences of dystrophin deletion on the hiPSC-fibroblasts gene expression and phenotype, we first analyzed fibroblast markers expressed by control

(See figure on next page.)

**Fig. 2** Transcriptome analysis of iPSC-derived cardiac fibroblasts from DMD and control patients. sPLS-DA was used to select genes that discriminate between DMD (black) and CON (grey) samples. iPSC derived cardiac fibroblasts were obtained from 3 DMD patients and 3 healthy donors. One control and one DMD were run in duplicate. **A:** Individual plot shows that component 1 (x-axis) captured the variance between groups and that component 2 (y-axis) consisted mainly of inter-group variation. **B:** Variable (gene) plot of selected genes (keepX set to 400 per component) either up in DMD (red) or up in CON (blue). The x and y axes correspond to the sample placement of Fig. 2A. Black points were discarded either because the y-position was more significant than the x-position or if that gene had a vip-score of < 1. **C:** Heatmap of differentially expressed genes selected from Fig. 2B. The row annotation colors are as in Fig. 2B (Up in DMD (red) and Up in CON (blue)). The column annotation is colored as in Fig. 2A. Normalized RNA expression is row-scaled so that red and blue indicate expression above or below the row mean. EnrichR was used to find significantly enriched pathways with the DMD up-regulated genes, as selected by sPLS-DA, using both **D:** Kyoto encyclopedia of genes and genomesdatabase (KEGG2021) as well as the **E:** human metabolic pathways database (HumanCyc 2016). Bar length is indicative of significance (– Log10(P.Value)). Principal enriched pathways include TCA cycle metabolism, gluconeogenesisand glycolysis
A

<figure><img src="image_6.png" /></figure>

B

<figure><img src="image_7.png" /></figure>

C

<figure><img src="image_8.png" /></figure>

D

KEGG 2021

<figure><img src="image_9.png" /></figure>

HumanCyc 2016

<figure><img src="image_10.png" /><figcaption>**Fig. 2** (See legend on previous page.)</figcaption></figure>
and DMD hiPSC-fibs by flow cytometry. All patients' fibroblasts exhibited a fibroblast-like morphology and were negative for CD31 and CD45, regardless of genotype (Data not shown). The control and DMD hiPSC-fibs were all positive for CD29 and CD105 whereas about 60 ± 10% expressed Thy1.1, and 4 ± 1% expressed PDGFRα with no significant differences between the genotypes (Fig. 1B–C).

### RNA expression analysis revealed altered metabolic pathways in DMD hiPSC-fibroblasts

To select gene expression altered by the loss of dystrophin expression in DMD and control hiPSC-fibs, we used sPLS-DA (keeping 400 variables on 2 components) from mixOmics in order to select genes with the most importance in separating samples based on DMD/CON (Component 1) as seen in the individual (sample) plot show in Fig. 2A (left panel), with DMD samples in black and CON samples in light grey. The results of the mixOmics vip program was used to further select the best genes by filtering with vip scores for component 1 above 1.0 and requiring that the x position was greater than the y position within the variable (gene) plot shown in Fig. 2A (right panel). The final set of genes selected as upregulated in DMD (red) or CON (blue) samples were then clustered with ward.D clustering of the Euclidean distance in a heatmap (Fig. 2B, same red/blue coloring as in Fig. 2A on rows and black/grey coloring on columns). The heatmap is showing a total of 323 genes, 159 up in DMD hiPSC-fibs patients and 164 up in CON hiPSC-fibs patients. The 159 up in DMD genes were used with EnrichR to select significantly enriched pathways upregulated in DMD fibroblasts. The KEGG (Fig. 2C) and HumanCyc (Fig. 2D) databases revealed many common DMD up-regulated pathways related to glycolysis, gluconeogenesis and the TCA cycle (Additional file 1: Fig. S2). In addition, these results suggest DMD fibroblasts may be suffering hypoxic stress as we found genes to be upregulated from the HIF-1 signaling pathway.

### Loss of dystrophin in DMD hiPSC-fibs enhanced glycolysis

To further explore metabolic modifications, we first analyzed glycolysis in control and DMD hiPSC-fibs using the glycolysis stress test on the Seahorse analyser. DMD hiPSC-fibs showed up-regulated extracellular acidification, glycolytic rate, and glycolytic capacity compared to n control hiPSC-fibs, indicating an important up-regulation of glycolysis in the absence of dystrophin isoforms (Fig. 3A).

Then, the expression of the principal enzymes involved in glycolysis was analyzed by immunoblot. Hexokinases are the initial enzyme of the glycolytic pathway, catalyzing the phosphorylation of glucose by ATP to glucose-6-P. Immunoblots revealed some heterogeneity within both DMD and control hiPSC-fibs, however no significant changes in hexokinase I, the main isoform expressed in hiPSC-fibs was observed in DMD as compared to control hiPSC-fibs. Similarly the pyruvate kinase isoforms (PKM and PKLR) that produce pyruvate from its substrate phosphoenolpyruvate were not differentially expressed in DMD or control hiPSC-fibs (Fig. 3B). However, a light up-regulation of lactate dehydrogenase A and a down regulation of lactate dehydrogenase B were observed in DMD hiPSC-fibs as compared to control hiPSC-fibs (Fig. 3B). While LDH isoform B (LDHB) catalyzes the transformation of lactate into pyruvate, LDH-A catalyzes the inverse reaction. Under poor oxidative phosphorylation situations, the fermentation of pyruvate to lactate by LDH-A, regenerates NAD⁺ from NADH in order to maintain glycolysis. The enhanced capacity of DMD hiPSC-fibs to produce lactate was further confirmed by the increase in cellular lactate concentration (Fig. 3C). The prevalence of the glycolytic pathway in DMD hiPSC-fibs was further supported by an increased expression of glucose transporter 1 (GLUT-1), emphasizing a preponderant use of glucose for cellular metabolism in the DMD hiPSC-fibs as compared to control hiPSC-fibs (Fig. 3D).

(See figure on next page.)

**Fig. 3** iPSC-derived cardiac fibroblasts from DMD patients exhibit high glycolytic capacity. **A:** The extracellular acidification rate (ECAR) profile of iPSC-derived cardiac fibroblasts subjected to a glycolysis stress test and the individual parameters for glycolysis, glycolytic capacity and glycolytic reserve. iPSC-derived cardiac fibroblasts were seeded to 80,000 cells/well and incubated for 24 h. ECAR was measured under basal conditions followed by the sequential addition of 10 mM glucose, 0.5 µM oligomycin and 100 mM 2-deoxyglucose (2-DG). Each data point represents an ECAR measurement obtained from the six biological samples of different batches of iPSC-derived cardiac fibroblasts. Data are expressed as mean ± SEM***p < 0.001 by non parametric t test. **B:** The protein expression levels of pyruvate kinase (muscular isoform-PKM), pyruvate kinase (L/R isoform), lactate dehydrogenase A and B, were quantified in iPSC-derived cardiac fibroblasts by western blotting and normalized to beta tubulin. Data are expressed as mean ± SEM of replicates of the six biological samples. **C:** Quantification of Intracellular L-lactate in iPSC-derived cardiac fibroblasts. Data are presented in mean ± SEM of replicates of the six biological samples, performed independently on different differentiation batches. ***p < 0.001 by non parametric t test. **D:** RT-PCR analysis of glucose transporter 1 mRNA expression levels. **p < 0.01 by non parametric t-test. **E:** Western blot analysis of phosphorylation on serine 293 of the PDH-E1α subunit in iPSC-derived cardiac fibroblasts. Data are presented in mean ± SEM of replicates of the six biological samples, performed independently on different differentiation batches.
A

<figure><img src="image_11.png" /></figure>

Glycolysis

<figure><img src="image_12.png" /></figure>

Glycolytic capacity

<figure><img src="image_13.png" /></figure>

Glycolytic reserve

<figure><img src="image_14.png" /></figure>

B

<figure><img src="image_15.png" /></figure>

PKM

<figure><img src="image_16.png" /></figure>

PKLR

<figure><img src="image_17.png" /></figure>

HK1

<figure><img src="image_18.png" /></figure>

LDH A

<figure><img src="image_19.png" /></figure>

LDH B

<figure><img src="image_20.png" /></figure>

C

Lactate

<figure><img src="image_21.png" /></figure>

D

GLUT1

<figure><img src="image_22.png" /></figure>

E

<figure><img src="image_23.png" /></figure>

Phospho-PDH E1a/ Total PDH E1a

<figure><img src="image_24.png" /><figcaption>**Fig. 3** (See legend on previous page.)</figcaption></figure>
**DMD hiPSC-fibs showed reduced mitochondrial respiration**

The up-regulated lactate production observed in DMD hiPSC-fibs may reflect altered mitochondria function and/or impaired glycolytic pyruvate to acetyl coA conversion depleting the mitochondrial supply of with energetic substrates. Thus, we assessed the expression of the pyruvate dehydrogenase complex and the phosphorylation of the PDHE1a subunit, which reduces the enzymatic activity of the complex. The DMD hiPSC-fibs showed decreased phosphorylation of the PDH E1a suggesting a more efficient ability to decarboxylate pyruvate into acetyl-CoA and, thereby feeding the TCA cycle and producing energy by oxidative phosphorylation (OXPHOS) (Fig. 3E). We next evaluated mitochondrial respiration of DMD hiPSC-fibs using the cell mito stress test on the Seahorse analyzer. A down-regulation of mitochondrial respiration was observed in DMDs, as evidenced by a decrease in the basal and maximal respiration, a lower spare respiratory reserve and mitochondrial ATP production (Fig. 4A). As mitochondrial respiration was reduced in DMD hiPSC-fibs, we determined the protein expression for components of the mitochondrial electron transport chain (ETC). Concordantly, the protein abundance of complex I subunit (NDUFB8), complex II Subunit (SDHB), complex III subunit (UQCRC2), complex IV (COXIV) and ATPase Synthase (ATP5A) were significantly diminished in DMD hiPSC-fibs (Fig. 4B). We next assessed the cellular ATP levels in DMD and control hiPSC-fibs in a luciferase-based luminescence reaction. Despite the drastically reduced capacity to generate ATP by oxidative phosphorylation, the total cellular ATP level trended to decrease in DMD hiPSC-fibs as compared to controls (Fig. 4C). However, ATP derived from glycolysis was significantly up-regulated in DMD hiPSC-fibs, further strengthening the increased use of glycolysis in DMD hiPSC-fibs (Fig. 4C).

All together, these results indicate that the loss of dystrophin in hiPSC-fibs mitigates mitochondrial

respiration and induced a metabolic switch toward glycolysis.

**DMD hiPSC-fibs exhibited lower expression of mitochondrial complex and mitochondrial number**

We next determined whether the decreased expression of enzymes in the electron transport chain resulted from a reduction in the number of mitochondria in DMD. We first analyzed the expression of Heat shock protein 60 (HSP60), a mitochondrial chaperone and Voltage-Dependent Anion Channel (VDAC1), also used as a mitochondrial protein marker. Although some heterogeneity was observed within the DMD hiPSC-fibs, a significant decreased in HSP60 and VDAC1 expression was observed in the DMD hiPSC-fibs as compared to control hiPSC-fibs (Fig. 5A). In order to more precisely determine the mitochondria content within hiPSC-fibs, mitochondrial DNA to nuclear DNA ratio (mtDNA/nDNA) was assessed. The median content of mitochondrial genes COX1 and ND1 was reduced by about two-fold in DMD hiPSC-fibs, as compared to hiPSC-fibs obtained from control subjects (for mitochondrial COX1/ nuclear B2 microglobulin; 56 ± 4 in DMD vs 112 ± 21 in control and for mitochondrial ND1/ nuclear B2 microglobulin; 90 ± 3 in DMD vs 148 ± 17 in control, *p* < 0.05 by Mann Whitney test) (Fig. 5B).

The reduced number of mitochondria observed in DMD could result from altered mitochondrial dynamics, which control the quality of mitochondria and maintain mitochondrial function. The molecular mechanisms underlying these processes rely on fusion (mitofusin 1 and 2) and fission proteins (dynamin-related protein), Western blot analysis revealed no significant modifications in the protein expression of DRP1 and mitofusin2 between DMD and control hiPSC-fibs (Fig. 5C–D).

(See figure on next page.)

**Fig. 4** Dystrophin loss impairs mitochondrial bioenergetics. **A:** Mitochondrial respiration of iPSC-derived cardiac fibroblasts was measured in three different DMD patients and three control patients. Mitochondrial stress test shows the OCR trace measured sequentially under basal conditions, following inhibition of ATP synthase (with oligomycin), uncoupling the electron transporter chain (ETC) with FCCP and blocking complex I and III with rotenone and antimycin A (Left graph). Basal respiration, maximal respiration, ATP production and spare respiratory capacity are reported in the left panels. Three independent experiments were performed. Non parametric t-test was applied to identify statistically significant differences between DMD and control. \* *p* < 0.05, \*\* *p* < 0.01. **B:** Western blot of respiratory complex subunits (OXPHOS) and quantification of experimental replicates of the six biological samples obtained from multiple differentiation batches. Non parametric t-test was applied to determine the statistically significant differences between DMD and control. \* *p* < 0.05, \*\* *p* < 0.01. **C:** Quantification of ATP levels in iPSC-derived cardiac fibroblasts from DMD and control patients (left panel). Cells were incubated for 60 min with sodium iodoacetate, a glycolysis inhibitor, to assess glycolysis-dependent ATP production (right panel). Two technical replicates from each of the three biological samples were performed independently. Non parametric t-test was applied to determine the statistically significant differences between DMD and control. \* *p* < 0.05
A

<figure><img src="image_25.png" /></figure>

Basal respiration

<figure><img src="image_26.png" /></figure>

Maximal respiration

<figure><img src="image_27.png" /></figure>

ATP production

<figure><img src="image_28.png" /></figure>

Spare respiratory capacity

<figure><img src="image_29.png" /></figure>

B

<figure><img src="image_30.png" /></figure>

ATP synthase sub α

<figure><img src="image_31.png" /></figure>

Com III Core 2

<figure><img src="image_32.png" /></figure>

C

<figure><img src="image_33.png" /></figure>

Com II sub 3

<figure><img src="image_34.png" /></figure>

Com I GRIM19

<figure><img src="image_35.png" /></figure>

<figure><img src="image_36.png" /></figure>

Glycolysis generated ATP

<figure><img src="image_37.png" /><figcaption>**Fig. 4** (See legend on previous page.)</figcaption></figure>
These results showed a reduced mitochondrial content in DMD hiPSC-fibs that did not appear to be associated with impaired mitochondrial fusion or fission.

### DMD hiPSC-fibs exhibited spatial mitochondrial rearrangements and modification in actin microfilament number

We performed immunofluorescent staining with the anti-mitochondrial antibody hsp60 to assess the mitochondrial network. In contrast to control hiPSC-fibs, which showed mitochondrial distribution throughout the cytoplasm of the cell, the mitochondria of DMD hiPSC-fibs were mainly clustered around the nucleus of the cell with poor dispersion throughout the cytosol (Fig. 6A). The altered intracellular distribution of mitochondria in DMD hiPSC-fibs was characterized by a significantly reduced branching and elongation of the mitochondrial network (Fig. 6A). As mitochondria utilize the cytoskeletal network to move throughout the cell, we stained and visualized filamentous actin with Phalloidin. DMD hiPSC-fibs revealed a modified organization of actin filaments, with a reduced number and randomly oriented filament bundles as compared to control hiPSC-fibs, which showed well organized bundles (Fig. 6B). This loss of filamentous actin assembly and coherency was independent of β-actin and γ-actin protein content that appeared similar in both DMD and control hiPSC-fibs (Fig. 6C). Similarly, β-tubulin protein expression, a component of microtubules, was comparable in DMD and control hiPSC-fibs (Fig. 6C). However, microtubule structure in DMD hiPSC-fibs was slightly thicker than in control hiPSC-fibs, suggesting that F-actin disorganization in DMD interferes with α- and β-tubulin assembly (Fig. 6C).

All together these results showed that dystrophin plays an indispensable role in bundling filamentous actin and mitochondria distribution within the cell.

### DMD hiPSC-fibs exhibited myofibroblast phenotype and a high sensitivity to pro-fibotic stress

We further examined the consequences on cell mitochondrial and microfilament rearrangements in DMD hiPSC-fibs by investigating the expression of alpha-smooth muscle actin (α-SMA) and collagen 1 proteins

used as markers of myofibroblasts. DMD hiPSC-fibs express higher levels of α-SMA protein compared to control hiPSC-fibs (Fig. 7A). In addition, more cells were positive for collagen 1, suggesting that loss of dystrophin was associated with a change in cellular phenotype (Fig. 7B). Then, we evaluated if the DMD hiPSC-fibs were more prone to differentiate into myofibroblasts under profibrotic stress, despite deficient actin cytoskeleton network. Treatment with the profibrotic factors TGF β or angiotensin II induced a higher expression of alpha-smooth muscle actin, COL1A1 and tenascin C in DMD hiPSC-fibs than in control hiPSC-fibs (Fig. 7C–D). These results showed that the dystrophin plays a cell-intrinsic function in restraining fibroblast activation.

## Discussion

In this study, we revealed the expression of dystrophin isoforms in hiPSC derived cardiac fibroblasts. A 427 kDa dystrophin protein was only found in control patients and not in the DMD patients who have out of frame deletions of exons or a spontaneous point mutation that induce a frameshift and a premature stop signal. The lack of dystrophin in human iPSC derived fibroblasts triggers a cascade of noxious cellular effects leading to cytoskeletal remodeling, modified energy supply, myofibroblast phenotype and acute sensitivity to pro-fibrotic stress.

A previous study performed with the mdx mouse model of Duchenne Disease showed that microtubule disorganization and mitochondrial dysfunction occur concurrently in dystrophin-deficient muscle. However, the link between cytoskeletal disorganization and mitochondrial bioenergetics was not uncovered [24]. The modified actin cytoskeleton in DMD most likely results from the loss of the actin binding domains (ABD). Specific domains of dystrophin have been shown to interact with microtubules [25, 26], but the absence of dystrophin had little effect on microtubule structure, suggesting that dystrophin plays only a minor role in microtubule arrangement. In contrast, the disorganized F-actin, caused by the loss of dystrophin in cardiac fibroblasts appeared to muddle distribution of mitochondria within the cell, along with a reduction in mitochondria number. The organelle utilizes the cytoskeletal network to move throughout the cell, while microtubules are used for long-range transport, the

(See figure on next page.)

**Fig. 5** iPSC-derived cardiac fibroblasts from DMD patients showed a reduced number of mitochondria. **A:** Western blot (left) and quantification(right) of HSP60 and VDAC1 protein expression in iPSC-derived cardiac fibroblasts. Quantification was performed on six biological samples obtained from multiple differentiation batches. Non parametric t-test was applied to determine the statistically significant differences between DMD and control. \*p < 0.05, \*\*p < 0.01. **B:** The relative abundance of mtDNA genes COX1 and ND1 using nuclear gene β globulin as a reference, was measured by PCR in hiPSC-fibroblasts derived from control and DMD patients. Non parametric t-test was applied to determine the statistically significant differences between DMD and control. \*p < 0.05. Western blot quantification of MFN2 (C) and DRP1 (D) protein expression in iPSC-derived cardiac fibroblasts. Non parametric t-test was applied to determine the statistically significant differences between DMD and control
<figure><img src="image_38.png" /></figure>

<figure><img src="image_39.png" /></figure>

<figure><img src="image_40.png" /></figure>

<figure><img src="image_41.png" /></figure>

<figure><img src="image_42.png" /></figure>

<figure><img src="image_43.png" /></figure>

<figure><img src="image_44.png" /></figure>

<figure><img src="image_45.png" /><figcaption>**Fig. 5** (See legend on previous page.)</figcaption></figure>
actin network is used for short-range movements. Additionally, the correct positioning and docking of mitochondria at locations of nutrient abundance increases the efficiency of ATP generation [27]. The reshaped distribution and the reduced number of mitochondria in DMD hiPSC-derived cardiac fibroblasts blunted mitochondrial respiration, which was further confirmed by the lower expression of ETC complexes. A recent study reported impaired energy metabolism as well as abnormal mitochondrial structure and function in iPSC-derived cardiomyocytes generated from adult DMD patients, but not from a young DMD patient [28]. Similarly, mitochondrial function and biogenesis were found compromised in cardiomyocytes from the Duchenne mouse model mdx with “humanized” telomeres [29] suggesting that mitochondrial dysfunction accompanies dystrophin loss and may account for the etiology of dystrophic heart failure.

The use of glycolysis in strong preference to oxidative phosphorylation in human iPSC derived fibroblasts from DMD patients likely reflects a compensatory mechanism to counterbalance the dysfunctional mitochondria as global cell ATP production was maintained. The metabolic remodeling induced by the loss of dystrophin in DMDs may be the trigger for enhanced myofibroblast differentiation. Indeed, glycolysis was found to be a critical metabolic enabler of cardiac fibroblast activation whereas glucose oxidation alleviates it [30]. For example, inhibition of glycolysis blunted the differentiation of lung fibroblasts into myofibroblasts and attenuated profibrotic phenotypes in myofibroblasts isolated from the lungs of patients with idiopathic pulmonary fibrosis [31]. The contribution of HIF-1α pathway in the metabolic switch observed in DMD hiPSC-derived cardiac fibroblasts was not examined in the present study and needs to be further explored as over-expression of hypoxia-inducible factor 1-α, was shown to potentiate myofibroblast differentiation [32].

In addition to be constitutively activated, DMD hiPSC-derived cardiac fibroblasts have an enhanced fibrotic response when stimulated with TGF-β or angiotensin II. Studies have shown that dysregulation of the normal wound healing response may also potentially lead to fibrosis. Indeed, depletion of activated cardiac fibroblasts by immunological targeting alleviates cardiac fibrosis and rescued cardiac contraction showing that cardiac fibroblasts are important regulator of heart physiopathology [33]. Therefore, the accumulation of fibrotic tissue in the heart of DMD patients may not only be secondary to cardiomyocyte death, but results from inappropriate activation of fibroblasts.

Our study highlights the relationship between cytoskeletal dynamics, metabolism of the cell and myofibroblast differentiation. It provides a new mechanism by which inactivation of dystrophin in non-cardiomyocyte cells worsens cardiac fibrosis and leads to increased severity of cardiopathy.

Limitations concerning the isoforms expressed in control and DMD hiPSC-derived cardiac fibroblasts should be acknowledged in our study. Whereas the highest isoform of dystrophin protein was persistently found in control hiPSC, smaller isoforms such as 71 kDa or 40 kDa could be also observed in control and DMD hiPSC-fibs patients who have large deletions encompassing multiple exons or a spontaneous point mutation. The transcription of mRNAs encoding the 71 kDa (Dp71) and the 40 kDa (Dp40) isoforms of dystrophin are initiated from the same promoter in intron 62 [34]. Dp40 results from alternative splicing in the C terminal portion of the protein. The positive bands observed at 71 kDa and 40 kDa suggest that the expression of Dp71 and Dp40 isoforms might be ubiquitous, however fragments generated by proteolysis cannot be excluded.

(See figure on next page.)

**Fig. 6** iPSC-derived cardiac fibroblasts from DMD patients displayed abnormal microfilament morphology and mitochondrial network. **A:** Representative 63 × confocal images of hiPSC-fibroblasts stained with anti-HSP 60 antibody to show the mitochondrial network in iPSC-derived cardiac fibroblasts generated from control and DMD patients. The dashed white lines show plasma membrane of the cells. Quantification of mitochondrial network using MiNA(Mitochondrial Network Analysis) workflow within Fiji, for the “mitochondrial footprint” or mitochondrial dispersion area (left) and the mitochondrial “networks” (right). A t-test was applied to determine the statistically significant differences between DMD and control. \*\*\*p < 0.0001. **B:** Representative 63 × confocal images of phalloidin-stained filamentous actin and subsequent quantification. **C:** Western blot and quantification of beta actin, gamma actin, β tubulin and Rho GDI protein expression in iPSC-derived cardiac fibroblasts. Data are presented as mean ± SEM from analyses obtained from different batches of cardiac fibroblasts differentiated from iPSC.Non parametric t-test was applied to determine the statistically significant differences between DMD and control. \*\*\*p < 0.001. **D:** Representative 63 × confocal images of hiPSC-fibroblasts stained with anti β-tubulin antibody showing the microtubule network in iPSC-derived cardiac fibroblasts generated from control and DMD patients
A

CON

<figure><img src="image_46.png" /></figure>

DMD

<figure><img src="image_47.png" /></figure>

Mitochondrial dispersion
area

<figure><img src="image_48.png" /></figure>

Mitochondrial networks

<figure><img src="image_49.png" /></figure>

B

CON

<figure><img src="image_50.png" /></figure>

DMD

<figure><img src="image_51.png" /></figure>

Phalloidin stain

<figure><img src="image_52.png" /></figure>

C

<figure><img src="image_53.png" /></figure>

β-actin / β-tubulin

<figure><img src="image_54.png" /></figure>

γ-Actin / β-Tubulin

<figure><img src="image_55.png" /></figure>

β-Tubulin/Rho GDi

<figure><img src="image_56.png" /></figure>

D

CON

<figure><img src="image_57.png" /></figure>

DMD

<figure><img src="image_58.png" /></figure>

Tubulin networks

<figure><img src="image_59.png" /><figcaption>**Fig. 6** (See legend on previous page.)</figcaption></figure>
<figure><img src="image_60.png" /></figure>

<figure><img src="image_61.png" /></figure>

<figure><img src="image_62.png" /></figure>

<figure><img src="image_63.png" /></figure>

<figure><img src="image_64.png" /></figure>

<figure><img src="image_65.png" /></figure>

<figure><img src="image_66.png" /></figure>

<figure><img src="image_67.png" /></figure>

<figure><img src="image_68.png" /></figure>

<figure><img src="image_69.png" /></figure>

<figure><img src="image_70.png" /><figcaption>**Fig. 7** iPSC-derived cardiac fibroblasts from DMD patients display a myofibroblast phenotype and an increased sensitivity to stress. A: Representative immunoblot of α-sma expression obtained in iPSC-derived cardiac fibroblasts generated from DMD and control patients and quantification. Data are presented as mean ± SEM of the six biological samples. Replicates were performed independently on different differentiation batches. B: Immunofluorescence staining of collagen I in iPSC-derived cardiac fibroblasts generated from DMD and control patients and quantification of α-SMA, COL1A and TNC expression in iPSC derived cardiac fibroblast treated with 10 ng/ml TGFβ (C) or 100 ng/ml angiotensin II (D) for 24 h. One-way ANOVA was applied to determine statistically significant differences between the basal condition and TGFβ or angiotensin II treatments</figcaption></figure>
**Supplementary Information**

The online version contains supplementary material available at https://doi.org/10.1186/s13062-023-00398-2.

**Additional file 1.** For data supplement figure 1: Western blot analysis of dystrophin protein expression in human iPSC-derived cardiac fibroblasts with affinity purified polyclonal antibodies (3 to 4 technical replicates were performed). Immunodetection of b-actin was used as loading control. A: Expression of the 427kDa protein was retrieved with polyclonal antibodies raised against aa 410-450, B: Expression of the 427kDa, the 71kDa and the 40kDa protein isoforms was retrieved with polyclonal antibodies raised against C-terminal end. For data supplement figure 2: EnrichR was used to identify significantly up-regulated pathways in DMD hiPSC-fibs. List of top genes up regulated in DMD hiPSC-fibs. and associated with metabolic pathways from the KEGG and HumanCyc databases.

**Acknowledgements**

The authors thank the GeT-Santé facility (I2MC, Inserm, Génome et Transcrip-tome, GenoToul, Toulouse, France) for the advice and technical contribution. We thank the flow cytometry facility (I2MC- Toulouse) and A. Zakaroff and E. Riant for their help and advice, the microscopy facility R. Flores- Flores for his excellent assistance and advice. The authors thank Cecilia Ferrantini (University of Florence), Alessandra Fornaro (Careggi University Hospital, Florence) and Prof. Yvan Torrente (Fondazione IRCCS Ca'Granda Ospedale Maggiore Poli-clinico; University of Milan) for patient recruitment.

**Author contributions**

Conceptualization and writing, original draft preparation AGOL and NP; formal analysis, Data collection and methodology, SS, LS, DRJI, MV, SM, AG, MR, and NP; review JI, AP, GP, patient recruitment: JMP and AF. All authors have read and agreed to the published version of the manuscript.

**Funding**

This work was supported by grant from ERA CVD. L.S. was supported by grants from ANR-PAUSE and the FRM-PAUSE (Fondation pour la Recherche Médicale).

**Availability of data and materials**

All data used for the current study are available from the corresponding author (nathalie.pizzinat@inserm.fr) upon reasonable request. Raw data of RNA-seq are uploaded to GEO database with accession number GSE237014.

**Declarations**

**Competing interests**

The authors declare no competing interests.

**Author details**

¹National Institute of Health and Medical Research (INSERM), I2MC, U1297, Toulouse, France. ²Unit of Vascular Biology and Regenerative Medicine, Centro Cardiologico Monzino IRCCS, Milan, Italy. ³National Institute of Health and Medical Research (INSERM) U1297 I2MC, Bioinformatic Core Facility, I2MC, Toulouse, France. ⁴Department of Biomedical, Surgical and Dental Sciences, Università Degli Studi di Milano, Milan, Italy. ⁵Department of Biology, Univer-sity of Florence, Florence, Italy. ⁶Neurology Unit, Fondazione IRCCS Ca'Granda Ospedale Maggiore Policlinico, Milan, Italy. ⁷University Toulouse III, 118 route de Narbonne, 31062 Toulouse, CEDEX 9, Toulouse, France.

Received: 29 June 2023 Accepted: 21 July 2023

Published: 27 July 2023

**References**

1. Yiu EM, Kornberg AJ. Duchenne muscular dystrophy. J Paediatr Child Health. 2015;51(8):759-64.

2. Duan D, Goemans N, Takeda S, Mercuri E, Aartsma-Rus A. Duchenne muscular dystrophy. Nat Rev Dis Primers. 2021;7(1):13.

3. Lin B, Li Y, Han L, Kaplan AD, Ao Y, Kalra S, et al. Modeling and study of the mechanism of dilated cardiomyopathy using induced pluripotent stem cells derived from individuals with Duchenne muscular dystrophy. Dis Model Mech. 2015;8(5):457-66.

4. Kamdar F, Garry DJ. Dystrophin-deficient cardiomyopathy. J Am Coll Cardiol. 2016;67(21):2533-46.

5. Palladino A, Papa AA, Morra S, Russo V, Ergoli M, Rago A, et al. Are there real benefits to implanting cardiac devices in patients with end-stage dilated dystrophinopathic cardiomyopathy? Review of literature and personal results. Acta Myol. 2019;38(1):1-7.

6. Florczyk-Soluch U, Polak K, Dulak J. The multifaceted view of heart problem in Duchenne muscular dystrophy. Cell Mol Life Sci. 2021;78(14):5447-68.

7. Delalande O, Molza AE, Dos Santos MR, Cheron A, Pollet E, Raguenes-Nicol C, et al. Dystrophin's central domain forms a complex filament that becomes disorganized by in-frame deletions. J Biol Chem. 2018;293(18):6637-46.

8. Gao QQ, McNally EM. The dystrophin complex: structure, function, and implications for therapy. Compr Physiol. 2015;5(3):1223-39.

9. O'Brien KF, Kunkel LM. Dystrophin and muscular dystrophy: past, present, and future. Mol Genet Metab. 2001;74(1-2):75-88.

10. Barboni MTS, Joachimsthaler A, Roux MJ, Nagy ZZ, Ventura DF, Rendon A, et al. Retinal dystrophins and the retinopathy of Duchenne muscular dystrophy. Prog Retin Eye Res. 2022:101137.

11. Iannotti FA, Pagano E, Guardiola O, Adinolfi S, Saccone V, Consalvi S, et al. Genetic and pharmacological regulation of the endocannabinoid CB1 receptor in Duchenne muscular dystrophy. Nat Commun. 2018;9(1):3950.

12. Tadayoni R, Rendon A, Soria-Jasso LE, Cisneros B. Dystrophin Dp71: the smallest but multifunctional product of the Duchenne muscular dystrophy gene. Mol Neurobiol. 2012;45(1):43-60.

13. Spaltro G, Vigorelli V, Casalnuovo F, Spinelli P, Castiglioni E, Rovina D, et al. Derivation of the Duchenne muscular dystrophy patient-derived induced pluripotent stem cell line lacking DMD exons 49 and 50 (CCMi001DMD-A-3, Δ49, Δ50). Stem Cell Res. 2017;25:128-31.

14. Rovina D, Castiglioni E, Farini A, Bellicchi M, Gervasini C, Paganini S, et al. Establishment of a Duchenne muscular dystrophy patient-derived induced pluripotent stem cell line carrying a deletion of exons 51-53 of the dystrophin gene (CCMi003-A). Stem Cell Res. 2019;40: 101544.

15. Rovina D, Castiglioni E, Mallia S, Rabino M, Farini A, Belicchi M, et al. Reprogramming of dermal fibroblasts from a Duchenne muscular dystrophy patient carrying a deletion of exons 45-50 into an induced pluripotent stem cell line (CCMi005-A). Stem Cell Res. 2022;64: 102889.

16. Patro R, Duggal G, Love MI, Irizarry RA, Kingsford C. Salmon provides fast and bias-aware quantification of transcript expression. Nat Methods. 2017;14(4):417-9.

17. Frankish A, Diekhans M, Jungreis I, Lagarde J, Loveland JE, Mudge JM, et al. Gencode 2021. Nucleic Acids Res. 2021;49(D1):D916-23. https://doi.org/10.1093/nar/gkaa1087.

18. Soneson C, Love MI, and Robinson MD. Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences. F1000Res. 2015:4:1521.

19. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139-40.

20. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7): e47.

21. Chen EY, Tan CM, Kou Y, Duan Q, Wang Z, Meirelles GV, et al. Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinformatics. 2013;14:128.

22. Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27-30.

23. Romero P, Wagg J, Green ML, Kaiser D, Krummenacker M, Karp PD. Computational prediction of human metabolic pathways from the complete human genome. Genome Biol. 2005;6(1):R2.

24. Ramos SV, Hughes MC, Delfinis LJ, Bellissimo CA, Perry CGR. Mitochondrial bioenergetic dysfunction in the D2mdx model of Duchenne muscular dystrophy is associated with microtubule disorganization in skeletal muscle. PLoS One. 2020;15(10):e0237138.
25. Belanto JJ, Mader TL, Eckhoff MD, Strandjord DM, Banks GB, Gardner MK, et al. Microtubule binding distinguishes dystrophin from utrophin. Proc Natl Acad Sci USA. 2014;111(15):5723-8.

26. Belanto JJ, Olthoff JT, Mader, T L. Chamberlain, C M. Nelson, D M. McCourt, P. Talsness, D M. Gundersen GG, Lowe, Dawn A, Ervasti JM. Independent variability of microtubule perturbations associated with dystrophinopathy. Hum Mol Genet. 2016; 25(22): 4951-4961. https://doi.org/10.1093/hmg/ddw318.

27. Pekkurnaz G, Trinidad JC, Wang X, Kong D, Schwarz TL. Glucose regulates mitochondrial motility via Milton modification by O-GlcNAc transferase. Cell. 2014;158(1):54-68.

28. Willi L, Abramovich I, Fernandez-Garcia J, Agranovich B, Shulman M, Milman H, et al. Bioenergetic and metabolic impairments in induced pluripotent stem cell-derived cardiomyocytes generated from Duchenne muscular dystrophy patients. Int J Mol Sci. 2022;23(17):9808.

29. Chang AC, Ong SG, LaGory EL, Kraft PE, Giaccia AJ, Wu JC, et al. Telomere shortening and metabolic compromise underlie dystrophic cardiomyopathy. Proc Natl Acad Sci USA. 2016;113(46):13120-5.

30. Chen ZT, Gao QY, Wu MX, Wang M, Sun RL, Jiang Y, et al. Glycolysis inhibition alleviates cardiac fibrosis after myocardial infarction by suppressing cardiac fibroblast activation. Front Cardiovasc Med. 2021;8: 701745.

31. Xie N, Tan Z, Banerjee S, Cui H, Ge J, Liu RM, et al. Glycolytic reprogramming in myofibroblast differentiation and lung fibrosis. Am J Respir Crit Care Med. 2015;192(12):1462-74.

32. Goodwin J, Choi H, Hsieh MH, Neugent ML, Ahn JM, Hayenga HN, et al. Targeting hypoxia-inducible factor-1alpha/pyruvate dehydrogenase kinase 1 axis by dichloroacetate suppresses bleomycin-induced pulmonary fibrosis. Am J Respir Cell Mol Biol. 2018;58(2):216-31.

33. Aghajanian H, Kimura T, Rurik JG, Hancock AS, Leibowitz MS, Li L, et al. Targeting cardiac fibrosis with engineered T cells. Nature. 2019;573(7774):430-3.

34. Muntoni F, Torelli S, Ferlini A. Dystrophin and mutations: one gene, several proteins, multiple phenotypes. Lancet Neurol. 2003;2(12):731-40.

**Publisher's Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Ready to submit your research? Choose BMC and benefit from:

* fast, convenient online submission

* thorough peer review by experienced researchers in your field

* rapid publication on acceptance

* support for research data, including large and complex data types

* gold Open Access which fosters wider collaboration and increased citations

* maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

<figure><img src="image_71.png" /></figure>