Determining the genetic basis of anthracycline-cardiotoxicity by response QTL mapping in induced cardiomyocytes

  1. David A Knowles  Is a corresponding author
  2. Courtney K Burrows
  3. John D Blischak
  4. Kristen M Patterson
  5. Daniel J Serie
  6. Nadine Norton
  7. Carole Ober
  8. Jonathan K Pritchard
  9. Yoav Gilad
  1. Stanford University, United States
  2. University of Chicago, United States
  3. Mayo Clinic, United States

Abstract

Anthracycline-induced cardiotoxicity (ACT) is a key limiting factor in setting optimal chemotherapy regimes, with almost half of patients expected to develop congestive heart failure given high doses. However, the genetic basis of sensitivity to anthracyclines remains unclear. We created a panel of iPSC-derived cardiomyocytes from 45 individuals and performed RNA-seq after 24h exposure to varying doxorubicin dosages. The transcriptomic response is substantial: the majority of genes are differentially expressed and over 6000 genes show evidence of differential splicing, the later driven by reduced splicing fidelity in the presence of doxorubicin. We show that inter-individual variation in transcriptional response is predictive of in vitro cell damage, which in turn is associated with in vivo ACT risk. We detect 447 response-expression QTLs and 42 response-splicing QTLs, which are enriched in lower ACT GWAS p-values, supporting the in vivo relevance of our map of genetic regulation of cellular response to anthracyclines.

Data availability

All the custom analysis scripts used for this project are available at https://github.com/davidaknowles/dox (Knowles and Blischak, 2017). The suez response eQTL mapping R package is available at https://github.com/davidaknowles/suez (Knowles, 2017). The following data are available as Supplementary Data: 1) differential expression cluster assignments, 2) significant (5% FDR) eQTLs and sQTLs, 3) differential splicing results, 4) levels of cardiac troponin and the predicted transcriptomic response. In addition to the Supplementary Data included with this paper, further results are hosted at Dryad (doi:10.5061/dryad.r5t8d04) including 1) gene-by-sample matrix of RNA-seq quantification (log counts per million), 2) LeafCutter intron excision quantification 3) p-values for all tested eQTLs, reQTLs, sQTLs, and rsQTLs, 4) RARG variant response and marginal trans-eQTLs, 5) RIN, RNA concentration and other technical covariates, 6) embryoid body imaging for all iPSC lines. The RNA-seq FASTQ files will be added to the dbGaP database (Tryka et al., 2014) under dbGaP accession phs000185 (https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000185). The genotype data files cannot be shared because releasing genotype data from a subset of individuals in the pedigree would enable the reconstruction of genotypes of other members of the pedigree, which would violate the original protocol approved by the research ethics board (Livne et al., 2015). The summary statistics for the ACT GWAS were given to us by the authors of the study (Schneider et al., 2016; Serie et al. 2017).

The following data sets were generated
    1. Knowles DA
    (2018) Expression and splicing quantification, eQTLs and sQTLs
    Available at Dryad Digital Repository under a CC0 Public Domain Dedication.

Article and author information

Author details

  1. David A Knowles

    Department of Genetics, Stanford University, Stanford, United States
    For correspondence
    knowles84@gmail.com
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-7408-146X
  2. Courtney K Burrows

    Department of Human Genetics, University of Chicago, Chicago, United States
    Competing interests
    The authors declare that no competing interests exist.
  3. John D Blischak

    Department of Human Genetics, University of Chicago, Chicago, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-2634-9879
  4. Kristen M Patterson

    Department of Human Genetics, University of Chicago, Chicago, United States
    Competing interests
    The authors declare that no competing interests exist.
  5. Daniel J Serie

    Department of Health Sciences Research, Mayo Clinic, Jacksonville, United States
    Competing interests
    The authors declare that no competing interests exist.
  6. Nadine Norton

    Department of Cancer Biology, Mayo Clinic, Jacksonville, United States
    Competing interests
    The authors declare that no competing interests exist.
  7. Carole Ober

    Department of Human Genetics, University of Chicago, Chicago, United States
    Competing interests
    The authors declare that no competing interests exist.
  8. Jonathan K Pritchard

    Department of Genetics, Stanford University, Stanford, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-8828-5236
  9. Yoav Gilad

    Department of Human Genetics, University of Chicago, Chicago, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-8284-8926

Funding

NIH Office of the Director (HL092206)

  • Yoav Gilad

Howard Hughes Medical Institute

  • Jonathan K Pritchard

NIH Office of the Director (HG008140)

  • Jonathan K Pritchard

NIH Office of the Director (HG009431)

  • Jonathan K Pritchard

NIH Office of the Director (TL1 TR 432-7)

  • Courtney K Burrows

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Ethics

Human subjects: Human Subjects work was approved by the University of Chicago IRB (protocol 10-416-B). Written informed consent was obtained from all participants.

Copyright

© 2018, Knowles et al.

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 4,334
    views
  • 521
    downloads
  • 94
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. David A Knowles
  2. Courtney K Burrows
  3. John D Blischak
  4. Kristen M Patterson
  5. Daniel J Serie
  6. Nadine Norton
  7. Carole Ober
  8. Jonathan K Pritchard
  9. Yoav Gilad
(2018)
Determining the genetic basis of anthracycline-cardiotoxicity by response QTL mapping in induced cardiomyocytes
eLife 7:e33480.
https://doi.org/10.7554/eLife.33480

Share this article

https://doi.org/10.7554/eLife.33480

Further reading

    1. Genetics and Genomics
    Wenjing Liu, Shujin Li ... Xianjun Zhu
    Research Article

    Familial exudative vitreoretinopathy (FEVR) is a severe genetic disorder characterized by incomplete vascularization of the peripheral retina and associated symptoms that can lead to vision loss. However, the underlying genetic causes of approximately 50% of FEVR cases remain unknown. Here, we report two heterozygous variants in calcyphosine-like gene (CAPSL) that is associated with FEVR. Both variants exhibited compromised CAPSL protein expression. Vascular endothelial cell (EC)-specific inactivation of Capsl resulted in delayed radial/vertical vascular progression, compromised endothelial proliferation/migration, recapitulating the human FEVR phenotypes. CAPSL-depleted human retinal microvascular endothelial cells (HRECs) exhibited impaired tube formation, decreased cell proliferation, disrupted cell polarity establishment, and filopodia/lamellipodia formation, as well as disrupted collective cell migration. Transcriptomic and proteomic profiling revealed that CAPSL abolition inhibited the MYC signaling axis, in which the expression of core MYC targeted genes were profoundly decreased. Furthermore, a combined analysis of CAPSL-depleted HRECs and c-MYC-depleted human umbilical vein endothelial cells uncovered similar transcription patterns. Collectively, this study reports a novel FEVR-associated candidate gene, CAPSL, which provides valuable information for genetic counseling of FEVR. This study also reveals that compromised CAPSL function may cause FEVR through MYC axis, shedding light on the potential involvement of MYC signaling in the pathogenesis of FEVR.

    1. Developmental Biology
    2. Genetics and Genomics
    Mitchell Bestry, Alexander N Larcombe ... David Martino
    Research Article

    Alcohol consumption in pregnancy can affect genome regulation in the developing offspring but results have been contradictory. We employed a physiologically relevant murine model of short-term moderate prenatal alcohol exposure (PAE) resembling common patterns of alcohol consumption in pregnancy in humans. Early moderate PAE was sufficient to affect site-specific DNA methylation in newborn pups without altering behavioural outcomes in adult littermates. Whole-genome bisulfite sequencing of neonatal brain and liver revealed stochastic influence on DNA methylation that was mostly tissue-specific, with some perturbations likely originating as early as gastrulation. DNA methylation differences were enriched in non-coding genomic regions with regulatory potential indicative of broad effects of alcohol on genome regulation. Replication studies in human cohorts with fetal alcohol spectrum disorder suggested some effects were metastable at genes linked to disease-relevant traits including facial morphology, intelligence, educational attainment, autism, and schizophrenia. In our murine model, a maternal diet high in folate and choline protected against some of the damaging effects of early moderate PAE on DNA methylation. Our studies demonstrate that early moderate exposure is sufficient to affect fetal genome regulation even in the absence of overt phenotypic changes and highlight a role for preventative maternal dietary interventions.