aiMeRA: A generic modular response analysis R package and its application to estrogen and retinoic acid receptors crosstalk

Modular response analysis (MRA) is a widely used modeling technique to uncover coupling strengths in molecular networks under a steady-state condition by means of perturbation experiments. We propose an extension of this methodology to search genomic data for new associations with a network modeled by MRA and to improve the predictive accuracy of MRA models. These extensions are illustrated by exploring the cross talk between estrogen and retinoic acid receptors, two nuclear receptors implicated in several hormone-driven cancers such as breast. We also present a novel, rigorous and elegant mathematical derivation of MRA equations, which is the foundation of this work and of an R package that is freely available at https://github.com/bioinfo-ircm/aiMeRA/. This mathematical analysis should facilitate MRA understanding by newcomers. Author summary Estrogen and retinoic acid receptors play an important role in several hormone-driven cancers and share co-regulators and co-repressors that modulate their transcription factor activity. The literature shows evidence for crosstalk between these two receptors and suggests that spatial competition on the promoters could be a mechanism. We used MRA to explore the possibility that key co-repressors, i.e., NRIP1 (RIP140) and LCoR could also mediate crosstalk by exploiting new quantitative (qPCR) and RNA sequencing data. The transcription factor role of the receptors and the availability of genome-wide data enabled us to explore extensions of the MRA methodology to explore genome-wide data sets a posteriori, searching for genes associated with a molecular network that was sampled by perturbation experiments. Despite nearly two decades of use, we felt that MRA lacked a systematic mathematical derivation. We present here an elegant and rather simple analysis that should greatly facilitate newcomers’ understanding of MRA details. Moreover, an easy-to-use R package is released that should make MRA accessible to biology labs without mathematical expertise. Quantitative data are embedded in the R package and RNA sequencing data are available from GEO.

[1]  Tapesh Santra,et al.  Reconstructing static and dynamic models of signaling pathways using Modular Response Analysis , 2018, Current Opinion in Systems Biology.

[2]  H. Leon Harter,et al.  Tables of Range and Studentized Range , 1960 .

[3]  Sarah L Vowler,et al.  Cooperative interaction between retinoic acid receptor-alpha and estrogen receptor in breast cancer. , 2010, Genes & development.

[4]  P. Kushner,et al.  Nuclear factor RIP140 modulates transcriptional activation by the estrogen receptor. , 1995, The EMBO journal.

[5]  M. Spinella,et al.  Negative Feedback at the Level of Nuclear Receptor Coregulation , 2003, Journal of Biological Chemistry.

[6]  Bertram Klinger,et al.  Modelling Signalling Networks from Perturbation Data , 2018, bioRxiv.

[7]  Trey Ideker,et al.  Cytoscape 2.8: new features for data integration and network visualization , 2010, Bioinform..

[8]  N. Harbeck,et al.  Importance of RIP140 and LCoR Sub-Cellular Localization for Their Association With Breast Cancer Aggressiveness and Patient Survival1 , 2018, Translational oncology.

[9]  C. Rodríguez,et al.  Complex regulation of LCoR signaling in breast cancer cells , 2017, Oncogene.

[10]  B. Katzenellenbogen,et al.  Integrative genomics of gene and metabolic regulation by estrogen receptors a and b , and their coregulators , 2013 .

[11]  Anthony C. Davison,et al.  Bootstrap Methods and Their Application , 1998 .

[12]  V. Giguère,et al.  Functional genomics identifies a mechanism for estrogen activation of the retinoic acid receptor alpha1 gene in breast cancer cells. , 2005, Molecular endocrinology.

[13]  K. White,et al.  Genomic Antagonism between Retinoic Acid and Estrogen Signaling in Breast Cancer , 2009, Cell.

[14]  Paul Theodor Pyl,et al.  HTSeq—a Python framework to work with high-throughput sequencing data , 2014, bioRxiv.

[15]  Cole Trapnell,et al.  TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions , 2013, Genome Biology.

[16]  E. Margeat,et al.  Negative regulation of estrogen signaling by ERβ and RIP140 in ovarian cancer cells. , 2013, Molecular endocrinology.

[17]  Margaret Warner,et al.  Estrogen receptors: how do they signal and what are their targets. , 2007, Physiological reviews.

[18]  Eduardo Sontag,et al.  Untangling the wires: A strategy to trace functional interactions in signaling and gene networks , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[19]  M. Spinella,et al.  Limiting Effects of RIP140 in Estrogen Signaling , 2005, Journal of Biological Chemistry.

[20]  Donald J. Wheeler,et al.  Understanding Statistical Process Control , 1986 .

[21]  John H. White,et al.  Ligand-dependent nuclear receptor corepressor LCoR functions by histone deacetylase-dependent and -independent mechanisms. , 2003, Molecular cell.

[22]  Davis J. McCarthy,et al.  Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation , 2012, Nucleic acids research.

[23]  S. Mader,et al.  The N-terminal of the estrogen receptor (ERα) mediates transcriptional cross-talk with the retinoic acid receptor in human breast cancer cells , 2003, The Journal of Steroid Biochemistry and Molecular Biology.

[24]  Bertram Klinger,et al.  Supplementary Materials for : Network quantification of EGFR signaling unveils potential for targeted combination therapy , 2022 .

[25]  J. Nicolas,et al.  Reporter cell lines are useful tools for monitoring biological activity of nuclear receptor ligands. , 2001, Luminescence : the journal of biological and chemical luminescence.