Global transcriptomic profiling of microcystin-LR or -RR treated hepatocytes (HepaRG).

The canonical mode of action (MOA) of microcystins (MC) is the inhibition of protein phosphatases, but complete characterization of toxicity pathways is lacking. The existence of over 200 MC congeners complicates risk estimates worldwide. This work employed RNA-seq to provide an unbiased and comprehensive characterization of cellular targets and impacted cellular processes of hepatocytes exposed to either MC-LR or MC-RR congeners. The human hepatocyte cell line, HepaRG, was treated with three concentrations of MC-LR or -RR for 2 h. Significant reduction in cell survival was observed in LR1000 and LR100 treatments whereas no acute toxicity was observed in any MR-RR treatment. RNA-seq was performed on all treatments of MC-LR and -RR. Differentially expressed genes and pathways associated with oxidative and endoplasmic reticulum (ER) stress, and the unfolded protein response (UPR) were highly enriched by both congeners as were inflammatory pathways. Genes associated with both apoptotic and inflammatory pathways were enriched in LR1000. We present a model of MC toxicity that immediately causes oxidative stress and leads to ER stress and the activation of the UPR. Differential activation of the three arms of the UPR and the kinetics of JNK activation ultimately determine whether cell survival or apoptosis is favored. Extracellular exosomes were enrichment of by both congeners, suggesting a previously unidentified mechanism for MC-dependent extracellular signaling. The complement system was enriched only in MC-RR treatments, suggesting congener-specific differences in cellular effects. This study provided an unbiased snapshot of the early systemic hepatocyte response to MC-LR and MC-RR congeners and may explain differences in toxicity among MC congeners.