Unraveling the genetic blueprint of doxorubicin-induced cardiotoxicity through systems genetics approaches

Figures

• 

  Fig. 1. Experimental scheme of phenotype collection, QTL mapping and MR analysis results in the BXD family of mice pre- and post-DOX treatment. ECHO, echocardiography; LV, left ventricular; QTL, quantitative trait locus; Chr, chromosome; LOD, logarithm of the odds

• 

  Fig. 2. The survival rates of BXD mice undergoing DOX treatment and the genetic mapping of these survival rates to the mouse genome. A The Kaplan–Meier curve illustrates the varying survival rates to DOX (Y-axis, %) among BXD strains. The X-axis represents days after DOX administration. Parental mice, B6 (black line) and D2 (red line), are indicated. The blue dashed lines denote BXD strains with longer survival rates than control B6 mice, while the green dashed lines represent BXD strains with shorter survival than D2 parental mice. B Manhattan plots display a significant 86–94 megabases (Mb) QTL on Chr10 (arrowhead) associated with survival rates from Day 6 to Day 10 in BXDs treated with DOX. The X-axis represents the chromosomal position in Mb on the mouse genome, while the Y-axis shows the peak LRS (likelihood ratio statistics) score. C A list of genetic variations in the top candidate genes within the Chr10 QTL identified that have genetic variations among BXD strains

• 

  Fig. 3. Rates of body weight (BW) loss in BXDs to DOX treatment and genetic mapping of BW loss to mouse genome. A Varied rates of mean BW loss in percent (%, Y-axis) among BXDs during 10 days of DOX treatment (X-axis). The BXD strains with significantly higher BW loss (green dashed lines) compared to parental D2 mice (red line) and BXDs with lower BW loss (blue dashed lines) than B6 control mice (black line) are indicated. B Values of BW in grams (g, Y-axis) on Day 10 post-DOX treatment among BXD mice (X-axis). Arrows indicate parental B6 and D2 mice. Asterisks (* = P < 0.05, ** = P < 0.01) indicate significant differences in BW compared with control B6 strains. C-D Manhattan plots showing a significant QTL of 52.5–59 Mb (arrowhead) on Chr19 associated with the BW loss mean in grams (g, C) and of 52.5–54.2 Mb in percent (%, D) in BXDs on Day 10. Y-axis, LRS and -logP scores

• 

  Fig. 4. Results of heart weight and echocardiographic assessment in BXDs treated with DOX. A-B Results of gross heart weight (HW) analysis in BXD mice on Day 10 post-DOX. Y-axis: Values of HW in grams (g, A) and heart weight to body weight (HW/BW, mg/g) ratio in arbitrary units (B) in BXD mice (X-axis). Arrows indicate parental B6 (black) and D2 (red) strains. C Representative echocardiographic long axis view images of BXD mouse heart at baseline (Day −1) and on Day 5 post-DOX treatment. Left ventricular (LV) walls and chamber length are indicated. D-F Results of comparative echocardiography assessment at the pre- and post-DOX treatment. Y-axis: Values of LV mass (mg, D), LV volumes (µl) at systole (s, E) and diastole (d, F) among BXDs (X-axis) at the baseline (Day −1, blue columns) compared to Day 5 post-DOX treatment (red columns). Asterisks denote significant differences (* = P < 0.05, ** = P < 0.01, *** = P < 0.001, **** = P < 0.0001) between Day −1 and Day 5

• 

  Fig. 5. Results of echocardiographic assessment of heart function in BXDs treated with DOX. Y-axis: A, ejection fraction (EF%); B fractional shortening (FS%); C, left atrial (LA) volumes (mm2); D heart rate (HR, beats/min); E cardiac output (CO, ml/min) among BXDs (X-axis) at baseline (Day −1, blue columns) compared to Day 5 post-DOX treatment (red columns). Asterisks denote significant differences between Day −1 and Day 5

• 

  Fig. 6. Results of GEMMA genetic mapping of echocardiography traits in BXDs in response to DOX. Manhattan plots showing a significant 103–120 Mb QTL (arrowhead): (A) on mouse genome (X-axis); and B on Chr14 associated with ejection fraction (EF%) and LV volumes at systole (LV Vol;s) and diastole (LV Vol;d). Y-axis indicates an -logP score

• 

  Fig. 7. Cardiac cell-type specific expression of the top candidate genes. The cell type specific expression based on single cell RNA sequencing was obtained from the Human Protein Atlas (HPA) database. The overlapping MR-test causal genes with other plots are not shown. Normalized TPM (nTPM) data corresponding to each cluster can be found in Supplementary Material 4.