Longitudinal dynamics of the tumor hypoxia response: From enzyme activity to biological phenotype.

Poor oxygenation (hypoxia) is a common spatially heterogeneous feature of human tumors. Biological responses to tumor hypoxia are orchestrated by the decreased activity of oxygen-dependent enzymes. The affinity of these enzymes for oxygen positions them along a continuum of oxygen sensing that defines their roles in launching reactive and adaptive cellular responses. These responses encompass regulation of all steps in the central dogma, with rapid perturbation of the metabolome and proteome followed by more persistent reprogramming of the transcriptome and epigenome. Core hypoxia response genes and pathways are commonly regulated at multiple inflection points, fine-tuning the dependencies on oxygen concentration and hypoxia duration. Ultimately, shifts in the activity of oxygen-sensing enzymes directly or indirectly endow cells with intrinsic hypoxia tolerance and drive processes that are associated with aggressive phenotypes in cancer including angiogenesis, migration, invasion, immune evasion, epithelial mesenchymal transition, and stemness.

Oxygen-independent disulfide bond formation in VEGF-A and CA9.

Low levels of oxygen (hypoxia) occurs in many (patho)physiological situations. Adaptation to hypoxia is in part mediated by proteins expressed in the extracellular space that mature in the endoplasmic reticulum (ER) prior to traversing the secretory pathway. The majority of such ER cargo proteins require disulfide bonds for structural stability. Disulfide bonds are formed co- and posttranslationally in a redox relay that requires a terminal electron acceptor such as oxygen. We have previously demonstrated that some ER cargo proteins such as low-density lipoprotein receptor (LDLR) and influenza hemagglutinin (Flu-HA) are unable to complete disulfide bond formation in the absence of oxygen, limiting their ability to pass ER quality control and their ultimate expression. Here, using radioactive pulse-chase immunoprecipitation analysis, we demonstrate that hypoxia-induced ER cargo proteins such as carbonic anhydrase 9 (CA9) and vascular endothelial growth factor A (VEGF-A) complete disulfide bond formation and mature with similar kinetics under hypoxia and normoxia. A global in silico analysis of ER cargo revealed that hypoxia-induced proteins on average contain fewer free cysteines and shorter-range disulfide bonds in comparison to other ER cargo proteins. These data demonstrate the existence of alternative electron acceptors to oxygen for disulfide bond formation in cellulo. However, the ability of different proteins to utilize an oxygen-independent pathway for disulfide bond formation varies widely, contributing to differential gene expression in hypoxia. The superior ability of hypoxia-induced proteins such as VEGF-A and CA9 to mature in hypoxia may be conferred by a simpler disulfide architecture.

Comparative toxicoproteogenomics of mouse and rat liver identifies TCDD-resistance genes.

The aryl hydrocarbon receptor (AHR) mediates many toxic effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). However, the AHR alone does not explain the widely different outcomes among organisms. To identify the other factors involved, we evaluated three transgenic mouse lines, each expressing a different rat AHR isoform (rWT, DEL, and INS) providing widely different resistance to TCDD toxicity, as well as C57BL/6 and DBA/2 mice which exhibit a ~ tenfold divergence in TCDD sensitivity (exposures of 5-1000 µg/kg TCDD). We supplement these with whole-genome sequencing, together with transcriptomic and proteomic analyses of the corresponding rat models, Long-Evans (L-E) and Han/Wistar (H/W) rats (having a ~ 1000-fold difference in their TCDD sensitivities; 100 µg/kg TCDD), to identify genes associated with TCDD-response phenotypes. Overall, we identified up to 50% of genes with altered mRNA abundance following TCDD exposure are associated with a single AHR isoform (33.8%, 11.7%, 5.2% and 0.3% of 3076 genes altered unique to rWT, DEL, C57BL/6 and INS respectively following 1000 µg/kg TCDD). Hepatic Pxdc1 was significantly repressed in all three TCDD-sensitive animal models (C57BL/6 and rWT mice, and L-E rat) after TCDD exposure. Three genes, including Cxxc5, Sugp1 and Hgfac, demonstrated different AHRE-1 (full) motif occurrences within their promoter regions between rat strains, as well as different patterns of mRNA abundance. Several hepatic proteins showed parallel up- or downward alterations with their RNAs, with three genes (SNRK, IGTP and IMPA2) showing consistent, strain-dependent changes. These data show the value of integrating genomic, transcriptomic and proteomic evidence across multi-species models in toxicologic studies.