Functional Validation of Doxorubicin-Induced Cardiotoxicity-Related Genes

Figures

• 

  Figure 1. Prioritization of Doxorubicin-Induced Cardiotoxicity–Associated Loci Table showing the 38 doxorubicin-induced cardiotoxicity–associated loci ranked on the basis of the single-nucleotide polymorphism (SNP) with highest P value from their respective publications.23, 24, 25, 26, 27, 28, 29, 30 The heat map of cardiac tissue expression shows the expression of anthracycline-induced cardiotoxicity–associated genes in adult human heart (n = 2), in fetal human heart (n = 2), and in human induced pluripotent stem cell–derived cardiomyocytes (hiPSC-CMs) (n = 7) by RNA sequencing. n = number of distinct patient-specific samples.

• 

  Figure 2. Validation of KOs of 36 Genes Significantly Associated With Doxorubicin-Induced Cardiotoxicity (A) Validation of successful clustered regularly interspaced short palindromic repeats/Cas9–based gene knockout (KO) in human induced pluripotent stem cell–derived cardiomyocytes (hiPSC-CMs) using quantitative reverse transcriptase polymerase chain reaction. The data presented here are the relative expression of each studied gene in its own KO lines relative to isogenic control (ISO). (B) Effect of gene KOs on hiPSC-CM viability after doxorubicin treatments (72 hours). Each data point represents 1 median lethal dose (LD50) calculation on the basis of an individual experimental replicate derived from a 5-log doxorubicin dosing. Central bar represents mean. Error bars represent ±SEM. P values (∗P < 0.05, ∗∗P < 0.01, ∗∗∗P < 0.001, and ∗∗∗∗P < 0.0001) are derived from differences between gene KO and ISO (Mann-Whitney U test). Raw dose-response curves from which these LD50 values are derived are provided in Supplemental Figures 4 to 6.

• 

  Central Illustration. Mechanistic Implications of Relevant Candidate Genes in DOX-Induced Cardiotoxicity The different mechanisms of cardiotoxicity after exposure to doxorubicin (DOX) are highlighted in colored boxes. AA = arachidonic acid; ABC = adenosine triphosphate–binding cassette; EET = epoxyeicosatrienoic acid; mRNA = messenger RNA; ROS = reactive oxygen species; SLC = solute carrier; SR = sarcoplasmic reticulum.

• 

  Figure 3. Reactive Oxygen Species Production, DNA Damage, and Iron Uptake (A-C) Hydrogen peroxide levels measured by ROS-Glo assay (luminescence) in hiPSC-CMs after doxorubicin treatments (24 hours). (D) Representative images for γH2AX immunofluorescent staining in hiPSC-CMs after treatments with doxorubicin (24 hours, 1 and 3 μM). (E) Quantification of DNA damage on the basis of γH2AX staining in hiPSC-CMs using flow cytometry (ISO, n = 5; PRDM2-KO, n = 5; and MLH1-KO, n = 5). (F) Effect of HFE knockout on hiPSC-CM iron uptake (ISO, n = 6; and HFE-KO, n = 6) measured using calcein staining. Error bars represent ±SEM. n = full independent experimental replicates. ∗P < 0.05, ∗∗P ≤ 0.01, ∗∗∗P < 0.001, and ∗∗∗∗P < 0.0001 by Mann-Whitney U test (A-C) and 2-way analysis of variance (E and F). DAPI = 4′,6-diamidino-2-phenylindole; IC50 = half maximal inhibitory concentration; other abbreviations as in Figure 2.

• 

  Figure 4. Doxorubicin Uptake, Calcium Handling, and Contractility (A) Doxorubicin uptake in hiPSC-CMs with knockouts of SLC and ABC transporters (isotype, n = 11; SLC28A3-KO, n = 4; SLC22A17-KO, n = 6; SLC28A1-KO, n = 4; ABCC2-KO, n = 4; ABCB4-KO, n = 7; ABCC5-KO, n = 4; ABCC9-KO, n = 4; and ABCC10-KO, n = 4). n = full independent experimental replicates. (B) Effect of knocking out MYH7 and CELF4 on calcium transients. Left: representative calcium transients. Middle left: full width at half maximum (FWHM). Middle right: calcium transient duration 75% (CTD75). Right: decay time. (C) Contractility analysis using impedance measurement in RIN3-KO, ZFN521-KO, and CYP2J2-KO hiPSC-CMs (ISO, n = 41; RIN3-KO, n = 48; ZFN521-KO, n = 48; and CYP2J2-KO, n = 89). n = number of assessed wells of 96-well plate derived from at least 3 independent rounds of differentiation (B and C). Error bars represent ±SEM. ∗P < 0.05, ∗∗P ≤ 0.01, ∗∗∗P < 0.001, and ∗∗∗∗P < 0.0001 by 2-way analysis of variance (A) and 1-way analysis of variance (B and C). Abbreviations as in Figure 2.