You are viewing the site in preview mode

Skip to main content

Alternative mRNA splicing in anthracycline-induced cardiomyopathy – a COG-ALTE03N1 report

Abstract

Background

Anthracycline-induced cardiomyopathy is a well-established adverse consequence in childhood cancer survivors. Altered mRNA expression in the peripheral blood has been found at the level of genes and pathways among anthracycline-exposed childhood cancer survivors with and without cardiomyopathy. However, the role of aberrant alternative splicing in anthracycline-induced cardiomyopathy remains unexplored. The present study examined if transcript-specific events, due to alternative splicing occur in anthracycline-exposed childhood cancer survivors with and without cardiomyopathy.

Methods

Participants were anthracycline-exposed childhood cancer survivors with cardiomyopathy (cases) matched with anthracycline-exposed childhood cancer survivors without cardiomyopathy (controls; matched on primary cancer diagnosis, year of diagnosis, and race/ethnicity). mRNA sequencing was performed on total RNA from peripheral blood in 32 cases and 32 matched controls. Event-level splicing tool, rMATS (replicate Multivariate Analysis of Transcript Splicing) was used for quantitative profiling of alternative splicing events.

Results

A total of 45 alternative splicing events in 36 genes were identified. Using a prioritization strategy to filter the alternative splicing events, intron retention in RPS24 and skipped exon of PFND5 showed differential expression of altered transcripts.

Conclusions

We identified specific alternative splicing events in anthracycline-exposed childhood cancer survivors with and without cardiomyopathy. Our findings suggest that differential alternative splicing events can provide additional insight into the peripheral blood transcriptomic landscape of anthracycline-induced cardiomyopathy.

Graphical abstract

Central Illustration. Aberrant alternative splicing and anthracycline-induced cardiomyopathy. This study sought to identify alternative splice variants that are differentially abundant between anthracycline-exposed childhood cancer survivors that developed cardiomyopathy (cases) versus those who did not (controls). We observed dysregulated alternative splicing of PFDN5 and RPS24 is associated with the development of cardiomyopathy.Splicing defects in PFDN5 impair cytoskeletal protein folding, while RPS24 dysregulation affects their translation, disrupting actin and tubulin homeostasis. Together, these alterations destabilize cardiomyocyte structure, contributing to sarcomere disorganization and the development of cardiomyopathy. Created in BioRender. Singh, P. (2025) https://BioRender.com/59zmgls.

Background

Anthracycline-induced cardiotoxicity represents a continuum from subclinical myocyte injury to congestive heart failure [1]. Studies of differential gene expression (DGE) in anthracycline-exposed cancer survivors employing microarrays [2] and RNA sequencing (RNA-Seq) [3, 4] from peripheral blood collected following cardiomyopathy diagnosis have identified molecular signatures and pathways associated with anthracycline-induced cardiomyopathy. However, DGE studies do not account for mRNA isoform/transcript diversity generated by differential or alternative splicing [5].

mRNA splicing is one of the regulatory mechanisms for gene expression and is essential to the flow of information from the DNA to protein in all eukaryotes [6]. More than 90% of human protein-coding genes produce multiple transcripts through alternative splicing [7, 8], resulting in diversification and expansion of the transcriptome and proteome [9]. Alternatively spliced isoforms may have related, distinct, or even opposing functions or subcellular localizations, or no protein product. There are five basic types of alternative splicing: alternative 5′ splice site (A5'SS), skipped exon (SE), mutually exclusive exons (MXE), retained intron (RI), and alternative 3′ splice site (A3'SS) [8]. Alternative splicing plays a role in biological processes such as cell differentiation and proliferation, organ development and stress response [10]. Aberrant splicing underlies many pathological processes including premature aging, infection, and inflammation, as well as in immune and metabolic disorders and cardiovascular disease [11,12,13]. Aberrant splicing of sarcomeric and ion channel genes has been seen in patients with cardiomyopathy in non-oncology settings [14,15,16]. The recent identification of several cardiac splice factors, such as RNA-binding motif protein 20 and 24, has provided insight as to how these splicing factors cause cardiomyopathy [17].

Both cardiomyopathy and anthracycline-induced cardiomyopathy share complex molecular mechanisms, [18] including oxidative stress, mitochondrial dysfunction, calcium dysregulation, and cell death, ultimately leading to cardiac dysfunction [19]. Alternative splicing plays a crucial role in cardiomyopathy by regulating key structural and functional proteins in the heart. In dilated cardiomyopathy (DCM), splicing of Titin (TTN) generates different isoforms, with the shorter, stiffer N2B isoform being optimal for adult heart function, while the longer, more compliant N2BA isoform predominates in fetal hearts and during heart failure [20]. Mutations in RNA binding motif protein 20 (RBM20), a key splicing regulator, shifts the balance toward N2BA, reducing cardiac contractility and contributing to DCM [21,22,23,24]. Changes in the splicing of OPA1, a gene regulating mitochondrial dynamics, can shift cardiomyocytes towards a pro-apoptotic phenotype, exacerbating cell death [25]. Li, et al. showed that calcium channel CaV1.2 is alternatively spliced in diabetes-induced cardiomyopathy [26]. CELF4 (CUGBP Elav-Like Family Member 4) is a member of the CELF family of RNA-binding proteins, which are key regulators of alternative splicing. It plays a critical role in developmental and tissue-specific splicing, especially in the heart. In the heart, CELF4 influences the splicing of key structural and contractile proteins, like TNNT2 (cardiac troponin T), which are essential for normal cardiac function. In transgenic MHC-CELFΔ mice, which express a dominant-negative CELF protein in the heart, impaired CELF activity leads to early-onset alternative splicing defects, resulting in cardiac hypertrophy, dilated cardiomyopathy, fibrosis, severe heart dysfunction, and premature death [27]. The condition is more pronounced in females, indicating possible sex-specific regulation of splicing. Importantly, restoring CELF function by overexpressing CUG-BP1 (a functional CELF protein) reverses the splicing defects and heart issues, confirming that CELF activity is vital for proper splicing and cardiac function in vivo. A specific SNP, rs1786814, in the CELF4 gene is associated with a tenfold increased risk of developing anthracycline-induced cardiomyopathy in survivors of childhood cancers, particularly in individuals with the ‘GG’ genotype [28]. This study also revealed a significant association between the presence of the ‘GG’ genotype and the coexistence of both the embryonic and adult splicing variants of TNNT2 in myocardial tissue from patients, suggesting that the continued expression of multiple troponin T isoforms enhances cardiotoxicity in response to treatment with high-dose anthracycline. Ragab et al., corroborated these findings that the rs1786814 polymorphism in CELF4 may influence alternative splicing of TNNT2, contributing to anthracycline-related cardiotoxicity in childhood cancer survivors [29]. They also showed that ‘GG’ genotype for rs1786814 was significantly associated with decreased ejection fraction and increased end-systolic diameter, indicating impaired cardiac function. Similarly, ‘GG’ genotype of rs17736312 in ROBO2 has been reported to be significantly associated with anthracycline-induced cardiomyopathy in survivors of childhood cancers [30]. ROBO1 and ROBO2 encode transmembrane Robo receptors that bind Slit ligands (SLIT). The Slit-Robo signaling pathway promotes cardiac fibrosis by interfering with the transforming growth factor-β1 [TGF-β1]/Smad pathway, resulting in disordered remodeling of the extracellular matrix and potentiating heart failure. Specifically, the Slit2-Robo1 pathway is activated in fibrotic heart tissues and acts as a cardiac fibrosis-promoting component [31, 32]. Distinct isoforms for ROBO2 [33] are reported and alternative splicing of ROBO1 has been shown in response to hypoxia in endothelial cells [34], albeit not in context of cardiomyopathy.

However, the role of alternative splicing in anthracycline-induced cardiomyopathy among childhood cancer survivors remains unexplored. We addressed this gap by conducting a differential splicing analysis of RNA-Seq data in anthracycline-exposed childhood cancer survivors with and without cardiomyopathy.

Methods

Study design

Participants were enrolled to a Children’s Oncology Group (COG) study ALTE03N1 (Key Adverse Events after Childhood Cancer, NCT00082745). The multicenter trial (through the Children’s Oncology Group) was approved by the National Cancer Institute Clinical Therapy Evaluation Program on 2/26/2004. The study was initiated at City of Hope National Medical Centre as coordinating center in October 2003 and moved to University of Alabama at Birmingham in January 2015. The participants included in this study were enrolled from 11/9/2005–11/17/2017. Each participating institution obtained local ethics committee approval prior to enrolling participants. Cases consisted of childhood cancer survivors who developed cardiomyopathy after exposure to anthracyclines. COG member institutions contributed participants to the study after obtaining approval from local institutional review boards. Written informed consent/assent was obtained from patients and/or parents/legal guardians for participants < 18y. City of Hope (IRB-03066) and the University of Alabama at Birmingham Institutional Review Board (IRB-150115006) approved all experimental protocols and methods. All methods were performed in accordance with the ethical standards of City of Hope and University of Alabama at Birmingham Institutional Review Board and with the 1964 Helsinki Declaration. For each case, one anthracycline-exposed survivor with no signs or symptoms of cardiomyopathy was randomly selected as a control from the same COG cohort, matched on primary cancer diagnosis, year of diagnosis (± 5 years), and race/ethnicity. The matching allowed control for confounding factors, by ensuring that cases and controls are similar on key characteristics allowing for more precise estimates of association between exposure (anthracycline) and outcome (cardiomyopathy). The selected controls also needed a longer duration of cardiomyopathy-free follow-up compared with the time from cancer diagnosis to cardiomyopathy for the corresponding case. Participants provided peripheral blood samples in PAXgene blood RNA tubes for germline RNA.

As described in previous publications [2,3,4, 28, 30, 35,36,37,38], case definition was based on echocardiographic parameters: left ventricular ejection fraction (LVEF) ≤ 40% and/or fractional shortening (SF) ≤ 28%. Presence (or absence) of signs or symptoms suggestive of congestive heart failure were documented and patients were classified as “symptomatic or asymptomatic”. Lifetime anthracycline exposure was calculated by multiplying the cumulative dose (mg/m2) of individual anthracyclines by a factor that reflects the drug’s cardiotoxic potential and then summing the results [39]. Radiation to the chest with the heart in the field was captured as a yes/no variable. A composite binary variable for cardiovascular risk factors (CVRFs) (yes [presence of any of the following: diabetes, hypertension, dyslipidemia]; no [absence of all CVRFs]) was ascertained through self-report.

RNA isolation, library construction, sequencing, and differential gene and transcript expression

Bioinformatic processing was performed as previously described [3]. Briefly, STAR [RRID:SCR_004463] was used to align the raw RNA-Seq fastq reads to the human reference genome from Gencode (GRCh38 p7 Release M25) [RRID:SCR_014966] [40]. For individual case–control comparisons, Cufflinks [RRID:SCR_014597)] was used on the aligned reads to assemble transcripts, estimate their abundance and test for differential expression and regulation [41, 42]. Cuffmerge [RRID:SCR_015688] was utilized to generate a transcriptome assembly and merge transcript data from all the samples. Finally, Cuffdiff [RRID:SCR_001647] was utilized to test for differences in transcript expression between cases and controls. We considered adjusted p-values that preserved the false discovery rate (FDR) at < 0.05 to be evidence of significant differential transcript expression.

The data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Series accession number GSE218276 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE218276).

Detection of differential splicing by rMATS

rMATS (v4.1.2, Update 12/17/2021; https://rnaseq-mats.sourceforge.io/rmats4.1.2/)) was used for alternative splicing analysis of the sorted, aligned sequence files [43,44,45]. rMATS was utilized in a fixed mode to detect five known and annotated splice events: A3'SS, A5'SS, MXE, RI and SE. Two versions of the rMATS results were generated. One version evaluated splicing variants with only those reads that spanned splice junctions (JC) and the second version included reads that spanned splice junctions and additionally those reads placed fully on the adjacent, alternatively spliced exon region (reads on target) (JCEC). rMATS quantifies splicing events as PSI (Percentage Spliced In), which is a ratio of reads specific to exon inclusion isoforms divided by the sum of reads specific to exon inclusion and exclusion isoforms. rMATS calculates the difference between PSI values (ΔPSI) between two groups under study, which serves as an effect size measure. rMATS uses a likelihood-ratio test to assess the statistical significance of ΔPSI between two groups, providing p-values and false discovery rate (FDR). The raw rMATS results were filtered using an FDR threshold of < 0.05 and inclusion level difference or delta percent spliced in (ΔPSI) > 0.1 or ΔPSI < −0.1 as the cutoff. For SE, A5'SS, A3'SS and RI, the results were imported based on JCEC. For MXE, the results were imported based on JC.

Visualization of alternative splicing in integrative genome browser

Aligned reads from participants were visualized by Sashimi plots generated using IGV genome browser [RRID:SCR_011793] [46, 47]. Read densities across exons and junction reads were plotted as ‘arcs’ that were annotated with the raw number of junctions reads present in each sample.

Gene and transcript prioritization

Prioritization was based on genes and respective transcript expression levels obtained for whole blood from GTEx portal [48]. We followed a conservative cutoff with stringent thresholds as our samples contained no replicates. Highly expressed transcripts with transcript per million (TPM) values ≥ 20 were examined further to see if alternative splicing events correlated with differential transcript expression. To evaluate the biological plausibility of the transcripts deemed differentially expressed between cases and controls, strength of evidence from literature was utilized. Tissue- and cell-specific isoform expression for each gene in human tissues was obtained from the GTEx Portal (http://gtexportal.org) and GTEx Analysis Release V8 (dbGaP Accession phs000424.v8.p2).

Gene set enrichment analyses

We performed the Gene-Disease Association dataset (GAD) disease enrichment analyses via the Database for Annotation, Visualization, and Integrated Discovery (DAVID) tool (RRID:SCR_001881; https://david.ncifcrf.gov/home.jsp) using default settings [49, 50].

Results

Patient characteristics

The median age at primary cancer diagnosis for the 32 cases and 32 matched controls was 9.5 years and 10.7 years, respectively (Table 1). Cases received a higher dose of chest radiation (1335.6 cGy vs. 524.1 cGy; P = 0.04) and were more likely to have had a CVRF (37.5% vs. 3.1%; P < 0.001). The median time between cancer diagnosis and cardiomyopathy for cases was 2.4 years (IQR: 0.6–7.7); controls were followed for a significantly longer period (median, 9.4 years; P < 0.001).

Table 1 Characteristics of anthracycline-exposed childhood cancer survivors

Summary of alternative splicing events

A total of 45 alternative splicing events were identified in 36 genes. Difference between PSI values (ΔPSI) between Cases and Controls, serves as an effect size measure which is shown as ‘Inclusion level difference’ (Table 2). SE events were the most frequent alternative splicing events, followed by RI, A5SS and A3SS events; only one MXE event was detected (Table 2). Six alternative splicing events (BISPR, PPP3 CB-AS1, CD27-AS1, SNHG8 and LINC00892) were in long-ncRNAs and one in a pseudogene (CASP4LP), while the remaining 38 events were in the protein-coding genes (84.4%).

Table 2 Thirty-six genes exhibiting evidence of alternative splicing

Functional annotation

Table 2 shows the results from the GAD analysis to identify human diseases associated with the 36 genes. Twenty-one of the 36 genes showed a disease association. The disease classes showing relevant associations included ‘CARDIOVASCULAR’, ‘METABOLIC’ and ‘PHARMACOGENOMIC’. ‘Myocardial Infarction’, ‘Cholesterol’, ‘Diabetes’, ‘Coronary Artery Disease’ and ‘Cardiomyopathy’ were the most relevant diseases associated with the alternatively spliced genes (Table 3).

Table 3 Disease annotation of 21 (of 36) genes using the Gene-Disease Associations dataset (GAD) from Database for Annotation, Visualization, and Integrated Discovery (DAVID)

Prioritization strategy for alternative splicing events

We excluded genes with a TPM value < 20 of the most abundant transcript/isoform in whole blood. We checked the remaining genes (PFDN5, RBM38, CD300A, DDX3X and RPS24) manually using IGV to identify genes with a change in expressed alternative splicing patterns between the cases and controls. To this end, we tested if there was a visible difference in used splice sites based on junction information, or a change in exon and intron expression. We checked whether the observed alternative splicing event was visually dominant in cases vs. controls and if differences in inclusion levels or inclusion junction counts (IJC) and skipped junction counts (SJC) for each case/control pair correlated with differential transcript expression. Alternative splicing events in two genes: RPS24 and PFND5 met the above criteria.

Sashimi plots depicting retained intron in RPS24 (S1 Fig) for representative samples are shown in Fig. 1. We found that the intron retention level was inversely correlated with the differential expression of RPS24 transcript ENST00000372360. Overall controls showed a higher level of intron retention compared to cases (Fig. 2). Most mRNA derived from retained intron event is degraded via nonsense-mediated decay (NMD) resulting in the observed decreased expression of the transcript in controls compared to cases. However, we observed that the alternative splicing event mainly caused dysregulation of RPS24 at transcript level rather than gene expression level (S1 Table and S3 Fig).

Fig. 1
figure 1

Sashimi plots in IGV genome browser depicting retained intron in RPS24 for representative samples. Cases are shown in red and matched controls are shown in blue. Genomic coordinates are plotted on x-axis and read density (whose value is configurable via IGV) on y-axis, and mRNA isoforms quantified are shown on bottom (exons in blue, introns as lines with arrow heads). Data for the plot were taken from three sets of paired samples. The plot highlights the differential splicing of the intron, which is present in the controls, but mostly absent in the case samples. Exon coverage max set to the same level for all samples. Junction coverage minimum >10. It is interesting to note that there is an abrupt decrease in the density of last exon in controls (marked by arrows)

Fig. 2
figure 2

Intron retention level is inversely correlated with the differential expression of RPS24 transcript ENST00000372360. The expression level of the transcript measured as FPKM (Fragments per Kilobase of transcript per Million mapped reads) is plotted on the Y-axis. The matched case–control sets are plotted side by side on the X-axis. Cases are indicated in red, and control subjects are indicated in green. The samples are ordered by pairwise inclusion level difference from rMATS to show correlation with transcript expression. Inclusion level difference of < −0.1 indicates that intron inclusion level in cases is lower compared to controls and inclusion level difference of > + 0.1 indicates higher level of inclusion in cases compared to controls. Transcript ENST00000613865 showed a similar pattern, albeit at lower FPKM values

Exon skipping was observed in PFDN5 (Fig. 3 and S2 Fig). ENST00000551018, the transcript with skipped exon showed higher expression in cases compared to controls (Fig. 4). The alternative splicing event caused dysregulation of PFDN5 both at transcript level and at gene expression level (Fig. 4 and S1 Table).

Fig. 3
figure 3

Sashimi plots in IGV genome browser depicting exon skipping in PFDN5. The plot highlights the differential splicing of the exon, which is largely present in the controls, but mostly absent in the case samples. Exon coverage max set to the same level for all samples

Fig. 4
figure 4

Correlation between difference in IJC in paired cases and controls and transcript expression for PFDN5 transcript ENST00000551018. The expression level of the transcript measured as FPKM is plotted on the Y-axis. The matched case–control sets are plotted side by side on the X-axis. Cases are indicated in red, and control subjects are indicated in green. The samples are ordered by pairwise inclusion junction count (IJC) difference from rMATS to show correlation with transcript expression

Discussion

Aberrant alternative splicing has been widely reported in a variety of disease states in the general population, including heart failure [13,14,15]. We performed an exploration of alternative splicing in whole blood samples from anthracycline-exposed childhood cancer survivors with and without cardiomyopathy, which enabled a quantitative assessment of alternative splicing events and a qualitative identification of genes affected by alternative splicing. We identified intron retention in RPS24 and exon skipping in PFDN5 with expression-level differences for transcripts harboring retained introns and skipped exons. RPS24 transcript ENST00000372360 and PFDN5 transcript ENST00000551018 had higher expression in patients with cardiomyopathy with direct correlation of intron retention and exon skipping to lower transcript level expression in controls. While RPS24 and PFDN5 are not traditionally recognized as major cardiac genes, their roles in protein synthesis, folding, and cellular stress responses suggest they could influence cardiac physiology and pathology.

Intron retention is an alternative splicing method whereby introns, rather than being spliced out, are retained in mature mRNAs. Intron retention interrupts the main open reading frame and may lead to inclusion of premature termination codons, whereby intron-retaining isoforms are often rapidly degraded by the NMD pathway, resulting in downregulation or no gene expression [51]. Exon skipping or ES is the most common alternative splicing event whereby a particular exon is omitted from an alternatively spliced mRNA [52].

Ribosomes are a complex of ribosomal RNA and ribosomal proteins that function as a machinery for mRNA translation and protein synthesis. Elevated protein synthesis rates are characteristic of proliferating cells that need new cellular constituents [53]. The eukaryotic ribosome enzyme system is composed of four ribosomal RNAs (rRNAs) and 80 ribosomal proteins (RPs). The 40S ribosomal protein S24 (RPS24) is one of the ribosomal proteins. Studies have shown that RPs have additional ribosomal functions including roles in DNA repair, replication, proliferation, apoptosis and chemoresistance [54]. Ribosomal dysfunction has been implicated in a variety of developmental disorders including congenital heart disease [55]. Several RPs have been implicated in the development and progression of cardiovascular disease in the general population [56, 57]. The Minute syndrome in Drosophila, is associated with RP haploinsufficiency, and is characterized by developmental delay, impaired growth, poor fertility, and cardiac dysfunction [58]. RNAi-mediated knockdown of RpS24 in Drosophila cardiac tissue reveals that the hearts of these larvae cease to contract within 48 h with accompanying cardiac atrophy and breakage of cardiac collagen [58]. These phenotypes suggest that RpS24 is essential for cardiac integrity. Interestingly, mutations in RPS24 cause Diamond Blackfan anemia; with a high prevalence of congenital heart disease (~ 30%) [59, 60]. In a previous study, characterization of mCpG in heart failure showed that RPS24 is associated with heart failure. [61] Kerry et al., have recently shown that alternative splicing of RPS24 results in long vs. short variants, where the long variant produces a more stable protein isoform that aids in hypoxic cell survival. [62] A splice altering variant in RPS24 (chr10:79800375G > A) [63] has been identified in patients with atrial septal defect, another congenital heart defect.

The Prefoldin Subunit 5 (PFDN5) encodes one of the six subunits of prefoldin, a molecular chaperone complex that binds and stabilizes newly synthesized polypeptides and regulates the folding of nascent actin and tubulin monomers, essential for cardiomyocyte integrity [64]. Splicing alterations in PFDN5 may impair the proper folding of actin and tubulin, leading to cytoskeletal abnormalities that contribute to heart failure and arrhythmias. Furthermore, since PFDNs contribute to cellular adaptive response to stress, its mis-splicing could disrupt cellular stress responses, making cardiomyocytes more susceptible to ischemic damage [65, 66]. Zhang et al. [67], showed that PFDN5 was upregulated in patients with chronic heart failure and is a promising biomarker for the prediction of heart failure. Similar to our study, Li et al. [68], showed that both RPS24 and PFDN5 were upregulated in individuals from the general population with heart failure. Chen et. al., identified RPS24 and PFDN5 as key hub genes that are dysregulated in hypertrophic cardiomyopathy samples compared to healthy controls [69]. Additionally, whole exome sequencing data showed that RPS24 mutations were associated with heart failure.

Hypoxia is a key regulator of cardiac hypertrophy and hypoxia also induces hypoxia-inducible factor 1-alpha (HIF1A) that in turn induces alternative splicing [34]. Alternative splicing events in RPS24 transcripts have been reported to be altered by hypoxia and favors hypoxic cell survival [70]. In mammals, the PFDN complex including PFDN5 binds to nascent actin and tubulin cytoskeletal proteins to deliver them to the chaperonin CCT to promote their folding [65]. In the heart, both actin and tubulin are crucial components of the cytoskeleton, playing vital roles in cardiac function and structure, with actin forming the sarcomeric units for contraction and tubulin forming microtubules that support cell shape and transport. In cardiomyopathy, the actin and tubulin cytoskeleton undergo remodeling, including changes in the density, stability, and post-translational modifications of microtubules [71]. Splicing alterations resulting in differential expression of splicing isoforms in RPS24 and PFDN5 could compromise protein homeostasis, cardiac structure, and stress response, ultimately contributing to cardiomyopathy (Central Illustration).

Many studies have examined alternative splicing in cardiovascular diseases [27, 28, 72, 73] and various genome editing and small molecules have been demonstrated to be able to correct alternative splicing patterns [74,75,76]. Manipulating dysregulated RNAs to adjust expression could improve the course of anthracycline-induced cardiomyopathy [76]. Identifying specific pathological alternative splicing targets is key to offering new insights into their therapeutic potential.

Limitations

Ideally, alternative splicing should be measured in the affected tissue (i.e., cardiac tissue). However, obtaining heart biopsies from cancer survivors is logistically challenging and not without risk. Recent evidence has shown that peripheral blood reflects the transcriptomic signature of other tissues including heart [77,78,79,80,81]. Prevalent case–control studies by the very nature of their design exclude fatal endpoints from the case set. Presence of survival bias risks under-ascertainment of genes associated with high lethality, with consequent underestimation of disease risk effect size for their alternative splicing events associated with both increased disease risk and disease-associated lethality. Further research should focus on replicating current findings in a larger sample size.

Conclusions

Our findings provide an insight into the altered splicing landscape of anthracycline-induced cardiomyopathy. We show that alternative splicing is a prominent feature in the blood obtained from patients with anthracycline-induced cardiomyopathy. We identified alternative transcripts of RPS24 and PFDN5 genes among childhood cancer survivors with anthracycline-induced cardiomyopathy. These genes have been previously implicated in the pathogenesis of cardiovascular diseases in the general population. Anthracycline-induced cardiomyopathy is a complex disease, where the relationship between gene expression and phenotypes is subject to various genetic and epigenetic influences prior to its clinical manifestation; examining alternative splicing events adds a piece to the whole puzzle. We therefore propose that examining alternative splicing should be included as part of the gene expression analysis as it provides additional insight into the transcriptomic landscape and could potentially allow a more accurate prediction of the functional consequences of detected changes in gene expression.

Data availability

The data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Series accession number GSE218276 https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc = GSE218276.

Abbreviations

COG:

Children’s Oncology Group

mRNA:

Messenger RiboNucleic Acid

RNA-Seq:

RNA Sequencing

AS:

Alternative Splicing

DGE:

Differential Gene Expression

A5SS:

Alternative 5′ Splice Site

SE:

Skipped Exon

MXE:

Mutually Exclusive Exon

RI:

Retained Intron

A3SS:

Alternative 3′ Splice Site

rMATS:

Replicate Multivariate Analysis of Transcript Splicing

PSI:

Percent Spliced In

IGV:

Integrative Genomics Viewer

TPM:

Transcript Per Million

GTEx:

Genotype-Tissue Expression

NMD:

Nonsense-mediated mRNA decay

rRNA:

Ribosomal RNA

FPKM:

Fragments per Kilobase of transcript per Million mapped reads

References

  1. Groarke JD, Nohria A. Anthracycline cardiotoxicity: a new paradigm for an old classic. Circulation. 2015;131(22):1946–9.

    Article  PubMed  Google Scholar 

  2. Singh P, Wang X, Hageman L, Chen Y, Magdy T, Landier W, et al. Association of GSTM1 null variant with anthracycline-related cardiomyopathy after childhood cancer-A Children’s Oncology Group ALTE03N1 report. Cancer. 2020;126(17):4051–8.

    Article  CAS  PubMed  Google Scholar 

  3. Singh P, Crossman DK, Zhou L, Wang X, Sharafeldin N, Hageman L, et al. Haptoglobin Gene Expression and Anthracycline-Related Cardiomyopathy in Childhood Cancer Survivors: A COG-ALTE03N1 Report. JACC CardioOncol. 2023;5(3):392–401.

    Article  PubMed Central  PubMed  Google Scholar 

  4. Singh P, Shah DA, Jouni M, Cejas RB, Crossman DK, Magdy T, et al. Altered Peripheral Blood Gene Expression in Childhood Cancer Survivors With Anthracycline-Induced Cardiomyopathy - A COG-ALTE03N1 Report. J Am Heart Assoc. 2023;12(19): e029954.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  5. Kelemen O, Convertini P, Zhang Z, Wen Y, Shen M, Falaleeva M, et al. Function of alternative splicing. Gene. 2013;514(1):1–30.

    Article  CAS  PubMed  Google Scholar 

  6. Shi Y. Mechanistic insights into precursor messenger RNA splicing by the spliceosome. Nat Rev Mol Cell Biol. 2017;18(11):655–70.

    Article  CAS  PubMed  Google Scholar 

  7. Pan Q, Shai O, Lee LJ, Frey BJ, Blencowe BJ. Deep surveying of alternative splicing complexity in the human transcriptome by high-throughput sequencing. Nat Genet. 2008;40(12):1413–5.

    Article  CAS  PubMed  Google Scholar 

  8. Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C, et al. Alternative isoform regulation in human tissue transcriptomes. Nature. 2008;456(7221):470–6.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  9. Birzele F, Csaba G, Zimmer R. Alternative splicing and protein structure evolution. Nucleic Acids Res. 2008;36(2):550–8.

    Article  CAS  PubMed  Google Scholar 

  10. Baralle FE, Giudice J. Alternative splicing as a regulator of development and tissue identity. Nat Rev Mol Cell Biol. 2017;18(7):437–51.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  11. Scotti MM, Swanson MS. RNA mis-splicing in disease. Nat Rev Genet. 2016;17(1):19–32.

    Article  CAS  PubMed  Google Scholar 

  12. Ren P, Lu L, Cai S, Chen J, Lin W, Han F. Alternative Splicing: A New Cause and Potential Therapeutic Target in Autoimmune Disease. Front Immunol. 2021;12: 713540.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  13. Hasimbegovic E, Schweiger V, Kastner N, Spannbauer A, Traxler D, Lukovic D, et al. Alternative Splicing in Cardiovascular Disease-A Survey of Recent Findings. Genes (Basel). 2021;12(9):1457.

  14. Lara-Pezzi E, Gomez-Salinero J, Gatto A, Garcia-Pavia P. The alternative heart: impact of alternative splicing in heart disease. J Cardiovasc Transl Res. 2013;6(6):945–55.

    Article  PubMed  Google Scholar 

  15. van den Hoogenhof MM, Pinto YM, Creemers EE. RNA Splicing: Regulation and Dysregulation in the Heart. Circ Res. 2016;118(3):454–68.

    Article  PubMed  Google Scholar 

  16. Zhu C, Chen Z, Guo W. Pre-mRNA mis-splicing of sarcomeric genes in heart failure. Biochim Biophys Acta Mol Basis Dis. 2017;1863(8):2056–63.

    Article  CAS  PubMed  Google Scholar 

  17. Beqqali A. Alternative splicing in cardiomyopathy. Biophys Rev. 2018;10(4):1061–71.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  18. Tanwar SS, Dwivedi S, Khan S, Sharma S. Cardiomyopathies and a brief insight into DOX-induced cardiomyopathy. Egypt Heart J. 2025;77(1):29.

    Article  PubMed Central  PubMed  Google Scholar 

  19. Xie S, Sun Y, Zhao X, Xiao Y, Zhou F, Lin L, et al. An update of the molecular mechanisms underlying anthracycline induced cardiotoxicity. Front Pharmacol. 2024;15:1406247.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  20. Radke MH, Peng J, Wu Y, McNabb M, Nelson OL, Granzier H, et al. Targeted deletion of titin N2B region leads to diastolic dysfunction and cardiac atrophy. Proc Natl Acad Sci U S A. 2007;104(9):3444–9.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  21. Vad OB, Angeli E, Liss M, Ahlberg G, Andreasen L, Christophersen IE, et al. Loss of Cardiac Splicing Regulator RBM20 Is Associated With Early-Onset Atrial Fibrillation. JACC Basic Transl Sci. 2024;9(2):163–80.

    Article  PubMed  Google Scholar 

  22. Pitt GS, Long Y. Mutations of Splicing Regulator RBM20 in Atrial Fibrillation. JACC Basic Transl Sci. 2024;9(2):181–4.

    Article  PubMed Central  PubMed  Google Scholar 

  23. Ma J, Lu L, Guo W, Ren J, Yang J. Emerging Role for RBM20 and its Splicing Substrates in Cardiac Function and Heart Failure. Curr Pharm Des. 2016;22(31):4744–51.

    Article  CAS  PubMed  Google Scholar 

  24. Briganti F, Wang Z. Alternative Splicing in the Heart: The Therapeutic Potential of Regulating the Regulators. Int J Mol Sci. 2024;25(23):13023.

  25. Olichon A, Elachouri G, Baricault L, Delettre C, Belenguer P, Lenaers G. OPA1 alternate splicing uncouples an evolutionary conserved function in mitochondrial fusion from a vertebrate restricted function in apoptosis. Cell Death Differ. 2007;14(4):682–92.

    Article  CAS  PubMed  Google Scholar 

  26. Li P, Qin D, Chen T, Hou W, Song X, Yin S, et al. Dysregulated Rbfox2 produces aberrant splicing of Ca(V)1.2 calcium channel in diabetes-induced cardiac hypertrophy. Cardiovasc Diabetol. 2023;22(1):168.

  27. Ladd AN, Taffet G, Hartley C, Kearney DL, Cooper TA. Cardiac tissue-specific repression of CELF activity disrupts alternative splicing and causes cardiomyopathy. Mol Cell Biol. 2005;25(14):6267–78.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  28. Wang X, Sun CL, Quinones-Lombrana A, Singh P, Landier W, Hageman L, et al. CELF4 Variant and Anthracycline-Related Cardiomyopathy: A Children’s Oncology Group Genome-Wide Association Study. J Clin Oncol. 2016;34(8):863–70.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  29. Ragab SM, El-Hawy MA, El-Hefnawy SM, El-Deeb HMA, Elfalah AS, Mahmoud AA. CELF4 (rs1786814) gene polymorphism and speckle-tracking Echocardiography for cardiovascular complications in childhood cancer survivors. Pediatr Res. 2024. https://doiorg.publicaciones.saludcastillayleon.es/10.1038/s41390-024-03400-3. Epub ahead of print.

  30. Wang X, Singh P, Zhou L, Sharafeldin N, Landier W, Hageman L, et al. Genome-Wide Association Study Identifies ROBO2 as a Novel Susceptibility Gene for Anthracycline-Related Cardiomyopathy in Childhood Cancer Survivors. J Clin Oncol. 2023;41(9):1758–69.

    Article  CAS  PubMed  Google Scholar 

  31. Liu Y, Yin Z, Xu X, Liu C, Duan X, Song Q, et al. Crosstalk between the activated Slit2-Robo1 pathway and TGF-beta1 signalling promotes cardiac fibrosis. ESC Heart Fail. 2021;8(1):447–60.

    Article  PubMed  Google Scholar 

  32. Zhao J, Mommersteeg MTM. Slit-Robo signalling in heart development. Cardiovasc Res. 2018;114(6):794–804.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  33. Yue Y, Grossmann B, Galetzka D, Zechner U, Haaf T. Isolation and differential expression of two isoforms of the ROBO2/Robo2 axon guidance receptor gene in humans and mice. Genomics. 2006;88(6):772–8.

    Article  CAS  PubMed  Google Scholar 

  34. Weigand JE, Boeckel JN, Gellert P, Dimmeler S. Hypoxia-induced alternative splicing in endothelial cells. PLoS ONE. 2012;7(8): e42697.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  35. Blanco JG, Sun CL, Landier W, Chen L, Esparza-Duran D, Leisenring W, et al. Anthracycline-related cardiomyopathy after childhood cancer: role of polymorphisms in carbonyl reductase genes–a report from the Children’s Oncology Group. J Clin Oncol. 2012;30(13):1415–21.

    Article  CAS  PubMed  Google Scholar 

  36. Sharafeldin N, Zhou L, Singh P, Crossman DK, Wang X, Hageman L, et al. Gene-Level Analysis of Anthracycline-Induced Cardiomyopathy in Cancer Survivors: A Report From COG-ALTE03N1, BMTSS, and CCSS. JACC CardioOncol. 2023;5(6):807–18.

    Article  PubMed Central  PubMed  Google Scholar 

  37. Singh P, Zhou L, Shah DA, Cejas RB, Crossman DK, Jouni M, et al. Identification of novel hypermethylated or hypomethylated CpG sites and genes associated with anthracycline-induced cardiomyopathy. Sci Rep. 2023;13(1):12683.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  38. Wang X, Liu W, Sun CL, Armenian SH, Hakonarson H, Hageman L, et al. Hyaluronan synthase 3 variant and anthracycline-related cardiomyopathy: a report from the children’s oncology group. J Clin Oncol. 2014;32(7):647–53.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  39. Feijen EAM, Leisenring WM, Stratton KL, Ness KK, van der Pal HJH, van Dalen EC, et al. Derivation of Anthracycline and Anthraquinone Equivalence Ratios to Doxorubicin for Late-Onset Cardiotoxicity. JAMA Oncol. 2019;5(6):864–71.

    Article  PubMed Central  PubMed  Google Scholar 

  40. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2013;29(1):15–21.

    Article  CAS  PubMed  Google Scholar 

  41. Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7(3):562–78.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  42. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511–5.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  43. Park JW, Tokheim C, Shen S, Xing Y. Identifying differential alternative splicing events from RNA sequencing data using RNASeq-MATS. Methods Mol Biol. 2013;1038:171–9.

    Article  CAS  PubMed  Google Scholar 

  44. Shen S, Park JW, Huang J, Dittmar KA, Lu ZX, Zhou Q, et al. MATS: a Bayesian framework for flexible detection of differential alternative splicing from RNA-Seq data. Nucleic Acids Res. 2012;40(8): e61.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  45. Shen S, Park JW, Lu ZX, Lin L, Henry MD, Wu YN, et al. rMATS: robust and flexible detection of differential alternative splicing from replicate RNA-Seq data. Proc Natl Acad Sci U S A. 2014;111(51):E5593–601.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  46. Robinson JT, Thorvaldsdottir H, Winckler W, Guttman M, Lander ES, Getz G, et al. Integrative genomics viewer. Nat Biotechnol. 2011;29(1):24–6.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  47. Katz Y, Wang ET, Silterra J, Schwartz S, Wong B, Thorvaldsdottir H, et al. Quantitative visualization of alternative exon expression from RNA-seq data. Bioinformatics. 2015;31(14):2400–2.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  48. Stanfill AG, Cao X. Enhancing Research Through the Use of the Genotype-Tissue Expression (GTEx) Database. Biol Res Nurs. 2021;23(3):533–40.

    Article  PubMed  Google Scholar 

  49. Sherman BT, Hao M, Qiu J, Jiao X, Baseler MW, Lane HC, et al. DAVID: a web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic Acids Res. 2022;50(W1):W216–21.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  50. da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.

    Article  CAS  PubMed  Google Scholar 

  51. Kumari A, Sedehizadeh S, Brook JD, Kozlowski P, Wojciechowska M. Differential fates of introns in gene expression due to global alternative splicing. Hum Genet. 2022;141(1):31–47.

    Article  CAS  PubMed  Google Scholar 

  52. Kim P, Yang M, Yiya K, Zhao W, Zhou X. ExonSkipDB: functional annotation of exon skipping event in human. Nucleic Acids Res. 2020;48(D1):D896–907.

    CAS  PubMed  Google Scholar 

  53. Ruggero D, Pandolfi PP. Does the ribosome translate cancer? Nat Rev Cancer. 2003;3(3):179–92.

  54. Zhou X, Liao WJ, Liao JM, Liao P, Lu H. Ribosomal proteins: functions beyond the ribosome. J Mol Cell Biol. 2015;7(2):92–104.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  55. Narla A, Ebert BL. Ribosomopathies: human disorders of ribosome dysfunction. Blood. 2010;115(16):3196–205.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  56. Alexander SJ, Woodling NS, Yedvobnick B. Insertional inactivation of the L13a ribosomal protein gene of Drosophila melanogaster identifies a new Minute locus. Gene. 2006;368:46–52.

    Article  CAS  PubMed  Google Scholar 

  57. Casad ME, Abraham D, Kim IM, Frangakis S, Dong B, Lin N, et al. Cardiomyopathy is associated with ribosomal protein gene haplo-insufficiency in Drosophila melanogaster. Genetics. 2011;189(3):861–70.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  58. Nim HT, Dang L, Thiyagarajah H, Bakopoulos D, See M, Charitakis N, et al. A cis-regulatory-directed pipeline for the identification of genes involved in cardiac development and disease. Genome Biol. 2021;22(1):335.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  59. Vlachos A, Blanc L, Lipton JM. Diamond Blackfan anemia: a model for the translational approach to understanding human disease. Expert Rev Hematol. 2014;7(3):359–72.

    Article  CAS  PubMed  Google Scholar 

  60. Vlachos A, Osorio DS, Atsidaftos E, Kang J, Lababidi ML, Seiden HS, et al. Increased Prevalence of Congenital Heart Disease in Children With Diamond Blackfan Anemia Suggests Unrecognized Diamond Blackfan Anemia as a Cause of Congenital Heart Disease in the General Population: A Report of the Diamond Blackfan Anemia Registry. Circ Genom Precis Med. 2018;11(5): e002044.

    Article  PubMed Central  PubMed  Google Scholar 

  61. Gilsbach R, Schwaderer M, Preissl S, Gruning BA, Kranzhofer D, Schneider P, et al. Distinct epigenetic programs regulate cardiac myocyte development and disease in the human heart in vivo. Nat Commun. 2018;9(1):391.

    Article  PubMed Central  PubMed  Google Scholar 

  62. Kerry J, Specker EJ, Mizzoni M, Brumwell A, Fell L, Goodbrand J, et al. Autophagy-dependent alternative splicing of ribosomal protein S24 produces a more stable isoform that aids in hypoxic cell survival. FEBS Lett. 2024;598(5):503–20.

  63. Jang MY, Patel PN, Pereira AC, Willcox JAL, Haghighi A, Tai AC, et al. Contribution of Previously Unrecognized RNA Splice-Altering Variants to Congenital Heart Disease. Circ Genom Precis Med. 2023;16(3):224–31.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  64. Liang J, Xia L, Oyang L, Lin J, Tan S, Yi P, et al. The functions and mechanisms of prefoldin complex and prefoldin-subunits. Cell Biosci. 2020;10:87.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  65. Herranz-Montoya I, Park S, Djouder N. A comprehensive analysis of prefoldins and their implication in cancer. iScience. 2021;24(11):103273.

  66. Tahmaz I, Shahmoradi Ghahe S, Topf U. Prefoldin Function in Cellular Protein Homeostasis and Human Diseases. Front Cell Dev Biol. 2021;9: 816214.

    Article  PubMed  Google Scholar 

  67. Shikuan Zhang XL, Xiang Xu, Qingxia Wei, Hanlu Li, Tong Liu and Chunlei L. Liu. Genomic and RNA-Seq Profiling of Patients with Heart Failure Identified Alterations in Phosphorylation and Immune Signaling Pathways. PREPRINT (Version 1) available at Research Square 2021. https://doiorg.publicaciones.saludcastillayleon.es/10.21203/rs.3.rs-625496/v1.

  68. Li X, Shen Y, Xu X, Guo G, Chen Y, Wei Q, et al. Genomic and RNA-Seq profiling of patients with HFrEF unraveled OAS1 mutation and aggressive expression. Int J Cardiol. 2023;375:44–54.

    Article  PubMed  Google Scholar 

  69. Chen P, Yawar W, Farooqui AR, Ali S, Lathiya N, Ghous Z, et al. Transcriptomics data integration and analysis to uncover hallmark genes in hypertrophic cardiomyopathy. Am J Transl Res. 2024;16(2):637–53.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  70. Kerry J, Specker EJ, Mizzoni M, Brumwell A, Fell L, Goodbrand J, et al. Autophagy-dependent alternative splicing of ribosomal protein S24 produces a more stable isoform that aids in hypoxic cell survival. FEBS Lett. 2024;598(5):503–20.

    Article  CAS  PubMed  Google Scholar 

  71. Caporizzo MA, Prosser BL. The microtubule cytoskeleton in cardiac mechanics and heart failure. Nat Rev Cardiol. 2022;19(6):364–78.

    Article  PubMed Central  PubMed  Google Scholar 

  72. Yang J, Hung LH, Licht T, Kostin S, Looso M, Khrameeva E, et al. RBM24 is a major regulator of muscle-specific alternative splicing. Dev Cell. 2014;31(1):87–99.

    Article  CAS  PubMed  Google Scholar 

  73. Lu SH, Lee KZ, Hsu PW, Su LY, Yeh YC, Pan CY, et al. Alternative Splicing Mediated by RNA-Binding Protein RBM24 Facilitates Cardiac Myofibrillogenesis in a Differentiation Stage-Specific Manner. Circ Res. 2022;130(1):112–29.

    Article  CAS  PubMed  Google Scholar 

  74. Nishiyama T, Zhang Y, Cui M, Li H, Sanchez-Ortiz E, McAnally JR, et al. Precise genomic editing of pathogenic mutations in RBM20 rescues dilated cardiomyopathy. Sci Transl Med. 2022;14(672):eade1633.

  75. Liu N, Olson EN. CRISPR Modeling and Correction of Cardiovascular Disease. Circ Res. 2022;130(12):1827–50.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  76. Cao J, Wei Z, Nie Y, Chen HZ. Therapeutic potential of alternative splicing in cardiovascular diseases. EBioMedicine. 2024;101: 104995.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  77. Wingrove JA, Daniels SE, Sehnert AJ, Tingley W, Elashoff MR, Rosenberg S, et al. Correlation of peripheral-blood gene expression with the extent of coronary artery stenosis. Circ Cardiovasc Genet. 2008;1(1):31–8.

    Article  CAS  PubMed  Google Scholar 

  78. Sinnaeve PR, Donahue MP, Grass P, Seo D, Vonderscher J, Chibout SD, et al. Gene expression patterns in peripheral blood correlate with the extent of coronary artery disease. PLoS ONE. 2009;4(9): e7037.

    Article  PubMed Central  PubMed  Google Scholar 

  79. Devaux Y. Transcriptome of blood cells as a reservoir of cardiovascular biomarkers. Biochim Biophys Acta Mol Cell Res. 2017;1864(1):209–16.

    Article  CAS  PubMed  Google Scholar 

  80. Bondar G, Cadeiras M, Wisniewski N, Maque J, Chittoor J, Chang E, et al. Comparison of whole blood and peripheral blood mononuclear cell gene expression for evaluation of the perioperative inflammatory response in patients with advanced heart failure. PLoS ONE. 2014;9(12): e115097.

    Article  PubMed Central  PubMed  Google Scholar 

  81. Basu M, Wang K, Ruppin E, Hannenhalli S. Predicting tissue-specific gene expression from whole blood transcriptome. Sci Adv. 2021;7(14).

Download references

Acknowledgements

The Genotype-Tissue Expression (GTEx) Project was supported by the Common Fund of the Office of the Director of the National Institutes of Health, and by NCI, NHGRI, NHLBI, NIDA, NIMH, and NINDS. The data used for the analyses described in this manuscript were obtained from the GTEx Portal on 08/12/2023. Juw Won Park for help with rMATS analysis.

Disclaimer

The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Funding

This research is supported by the National Cancer Institute (R35 CA220502; Principal Investigator [PI], S. Bhatia), Leukemia and Lymphoma Society (6563–19; PI, S. Bhatia), and the V Foundation (DT2019-010; PI, S. Bhatia). The Children’s Oncology Group study (COG-ALTE03 N1 [NCT00082745]; PI, S. Bhatia) reported here is supported by the National Clinical Trials Network Operations Center Grant (U10 CA180886; PI, D.S. Hawkins), the National Clinical Trials Network Statistics & Data Center Grant (U10 CA180899; PI, Alonzo), the National Cancer Institute Community Oncology Research Program Grant (UG1 CA189955; PI, Roth, and the Community Clinical Oncology Program Grant (U10 CA095861; PI, Pollock), and the St. Baldrick’s Foundation through an unrestricted grant.

Author information

Authors and Affiliations

Authors

Contributions

P.S., D.C. and S.B. conceived and designed the study; P.S. and D.C. acquired and performed data analysis; P.S., D.C., C.C., P.T., N.S., X.W., L.Z., L.H., S.A., F.B., D.H., F.K., M.H., J.N., J.G., W.L., and S.B. interpreted the results, drafted, and revised the final manuscript.

Corresponding author

Correspondence to Purnima Singh.

Ethics declarations

Ethics approval and consent to participate

COG member institutions contributed participants to the study after obtaining approval from local institutional review boards. Written informed consent/assent was obtained from patients and/or parents/legal guardians. City of Hope (IRB-03066) and The University of Alabama at Birmingham Institutional Review Board (IRB-150115006) approved all experimental protocols and methods. All methods were performed in accordance with the ethical standards of City of Hope and University of Alabama at Birmingham Institutional Review Board and with the 1964 Helsinki Declaration.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Singh, P., Crossman, D.K., Cheng, C. et al. Alternative mRNA splicing in anthracycline-induced cardiomyopathy – a COG-ALTE03N1 report. Cardio-Oncology 11, 47 (2025). https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s40959-025-00345-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doiorg.publicaciones.saludcastillayleon.es/10.1186/s40959-025-00345-2

Keywords