Helicase-like transcription factor-deletion from the tumor microenvironment in a cell line-derived xenograft model of colorectal cancer reprogrammed the human transcriptome-S-nitroso-proteome to promote inflammation and redirect metastasis

Abstract

Department of Pathology, CHRISTUS Mother Frances Hospital, Tyler, Texas, United States of America Affiliation: Department of Pathology, Texas Tech University Health Sciences Center, Lubbock, Texas, United States of America Jannette M. Dufour Roles Resources, Writing – review & editing Affiliation: Department of Cell Biology & Biochemistry, Texas Tech University Health Sciences Center, Lubbock, Texas, United States of America Beverly S. Chilton Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Resources, Software, Supervision, Validation, Visualization, Writing – original draft * E-mail: [email protected] Affiliations Department of Cell Biology & Biochemistry, Texas Tech University Health Sciences Center, Lubbock, Texas, United States of America, Cancer Center, School of Medicine, Texas Tech University Health Sciences Center, Lubbock, Texas, United States of America ORCID logo https://orcid.org/0000-0002-8066-8697 Introduction Metastasis is responsible for nearly 90% of deaths from cancer [1]. Of the three nitric oxide synthase (NOS) isoforms, neuronal (nNOS, NOS1) and endothelial (eNOS, NOS3) are calcium-dependent isoforms that produce nanomolar concentrations of NO for seconds or minutes when activated. Results Experimental timeline ensured mouse survival There is a negative correlation between HLTF expression in tumor cells and survival in mice [33]. [...]when HLTF+/+ human HCT116 Red-FLuc cells were used to establish an orthotopic xenograft model in Hltf+/+ (control mice) and Hltf-deleted mice a timeline of 35 days was established during which primary tumor size and metastasis was assessed weekly by BLI. [...]three Hltf-deleted mice died. [...]the hazard ratio (Mantel-Haenszel) indicated Hltf-deleted mice at any time during the treatment protocol were 5-times more likely to die than control mice.

About the Authors:

Rebecca A. Helmer

Roles Formal analysis, Investigation, Methodology, Validation, Writing – review & editing

Current address: Garrison Independent School District, Garrison, Texas, United States of America

Affiliation: Department of Cell Biology & Biochemistry, Texas Tech University Health Sciences Center, Lubbock, Texas, United States of America

Raul Martinez-Zaguilan

Roles Data curation, Formal analysis, Investigation, Methodology, Software, Visualization, Writing – review & editing

Affiliation: Department of Cell Physiology & Molecular Biophysics, Texas Tech University Health Sciences Center, Lubbock, Texas, United States of America

Gurvinder Kaur

Roles Data curation, Formal analysis, Methodology, Validation, Visualization

Affiliation: Department of Medical Education, Texas Tech University Health Sciences Center, Lubbock, Texas, United States of America

Lisa A. Smith

Roles Investigation, Methodology, Validation, Visualization

Current address: Department of Pathology, CHRISTUS Mother Frances Hospital, Tyler, Texas, United States of America

Affiliation: Department of Pathology, Texas Tech University Health Sciences Center, Lubbock, Texas, United States of America

Jannette M. Dufour

Roles Resources, Writing – review & editing

Affiliation: Department of Cell Biology & Biochemistry, Texas Tech University Health Sciences Center, Lubbock, Texas, United States of America

Beverly S. Chilton

Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Resources, Software, Supervision, Validation, Visualization, Writing – original draft

* E-mail: [email protected]

Affiliations Department of Cell Biology & Biochemistry, Texas Tech University Health Sciences Center, Lubbock, Texas, United States of America, Cancer Center, School of Medicine, Texas Tech University Health Sciences Center, Lubbock, Texas, United States of America

ORCID logo https://orcid.org/0000-0002-8066-8697

Introduction

Metastasis is responsible for nearly 90% of deaths from cancer [1]. This estimate has been constant for more than 50 years, and CRC—the third most common cancer diagnosed in men and women worldwide [2]—is a prototypical example. Early screening tools, treatment modalities and lifestyle changes have improved survival. However, metastasis remains the main cause of CRC-related mortality mainly because 20–25% of CRCs are metastatic at initial diagnosis [3]. In the United States, the 5-year survival statistic for individuals with metastatic disease is 14%, compared to those with regional (71%) and localized (90%) CRC [4], and mortality rates are higher in men than women [5]. Understanding the mechanism responsible for CRC progression and metastasis is essential to the design of treatments to improve patient prognosis.

HLTF protein—the multidomain mammalian ortholog of yeast Rad5—has a DNA-binding domain, a HIRAN domain, and a C3HC4-type RING-finger embedded between SNF2 helicase motifs [6]. HLTF cDNA was cloned based on the protein DNA-binding capabilities. The cloning was independently done by three research groups, Sheridan et al. [7], Hayward-Lester et al. [8] and Ding et al. [9] who initiated the use of the HLTF name. HLTF mediates replication fork reversal to rescue a stalled replication fork under stress via its HIRAN domain [10], and is a ubiquitin ligase targeting PCNA for polyubiquitination via its RING domain [11]. HLTF is a putative tumor suppressor because of its role in DNA damage tolerance mechanisms [12]. Epigenetic silencing of HLTF has been implicated in the malignant transformation of adenomas to carcinomas [13] because the HLTF gene is silenced by hypermethylation in 43% of colon cancers [14]. Methylated HLTF DNA in serum is significantly correlated with tumor size, more aggressive tumors, advanced stage (III or IV) metastatic disease, including micro-metastasis, and shorter survival [15–18]. Progressive loss of HLTF expression from tumor cells correlates with tumor staging [19]. There are higher levels of HLTF protein in tumors from patients with early-stage 1 CRC compared to patients with late-stage IV disease. In contrast, immunohistochemistry showed negligible HLTF protein in fibroblasts of the tumor microenvironment (TME) regardless of tumor stage. These findings were reinforced by reference component RNAseq analysis of single-cell transcriptomes from eleven primary tumors and matched normal mucosa [20]. Three fibroblast populations—normal, myofibroblasts, and cancer-associated fibroblasts—were identified based on gene expression profiles. HLTF was not detected in any of the fibroblast populations despite measured changes in mRNA levels between normal and tumor cells. Epigenetic silencing of HLTF in tumor cells has long been the experimental focus to the exclusion of HLTF in the TME [21].

CRC progression is accompanied by overexpression of nitric oxide (NO)—a diffusible/pleiotropic regulator of many normal cellular processes—implicated in the activation of oncogenic signaling pathways [22, 23]. NO is endogenously generated by three different isoforms of the enzyme NO synthase (NOS) using L-arginine and molecular oxygen as substrates. Of the three nitric oxide synthase (NOS) isoforms, neuronal (nNOS, NOS1) and endothelial (eNOS, NOS3) are calcium-dependent isoforms that produce nanomolar concentrations of NO for seconds or minutes when activated. In contrast, the inducible (iNOS, NOS2) isoform is calcium-independent. and produces micromolar concentrations of NO for hours or days [24]. NO induces post-translational coupling of a nitroso moiety to a reactive cysteine leading to protein S-nitrosylation that affects protein-protein interactions [25]. This pathway is independent of the NO/soluble guanylate cyclase (sGC)/cGMP-dependent protein kinase (PKG) signaling pathway. The NO signal can be stored and propagated at nitrosyl adducts at specific cysteine sites of proteins. S-nitrosylation is purported to be site-selective. INOS-S100A8/A9 targets a conserved I/L-X-C-X2-D/E motif [26]. However, S-nitrosylation is a reversible covalent chemical reaction and therefore one of the more difficult post-translational modifications to study. Progress in this regard is owed to the three fundamental steps of the biotin-switch assay [27]: free sulfhydryls are chemically blocked, nitrosylated cysteines are selectively reduced, and a biotin adduct is switched for the NO adduct. Recent modifications substitute iodoTMT reagents [28] in the biotin-switch procedure and use mass spectrometry to identify and quantify specific S-nitrosylated residues.

To investigate the functional importance of negligible HLTF expression in fibroblasts of the TME, we developed an HCT116 cell line-derived xenograft model of metastatic CRC. Primary tumor xenografts were established in Hltf-deleted and control male mice by direct OCMI of HLTF+/+HCT116 Red-FLuc cells [29]. The HCT116 human colon carcinoma cell line was selected because it expresses HLTF [30], displays leading edge invasion in xenograft tumor model [31], and recapitulates the multi-step dissemination process to the liver and lungs [32]. However, HLTF+/+HCT116 cells in Hltf-deleted mice shifted their metastatic direction. Species-specific RNAseq analysis of primary CDX tumors arising from the same passage of HLTF+/+HCT116 Red-FLuc cells in the TME of Hltf-deleted and control mice revealed striking coordination of an inflammation pathway. Transactivation of iNOS expression and the iNOS-S100A8/A9 signaling axis in Hltf-deleted TME drove the comparison of S-nitroso-proteome of primary tumor with metastatic tumors in Hltf-deleted TMEs to reveal remarkable spatiotemporal continuity in S-nitrosylation signaling.

Results

Experimental timeline ensured mouse survival

There is a negative correlation between HLTF expression in tumor cells and survival in mice [33]. Therefore, when HLTF+/+ human HCT116 Red-FLuc cells were used to establish an orthotopic xenograft model in Hltf+/+ (control mice) and Hltf-deleted mice a timeline of 35 days was established during which primary tumor size and metastasis was assessed weekly by BLI. At necropsy, primary tumors were ≤ 10 mm, i.e. below the maximum allowable tumor size of 20 mm tumor development/metastasis did not interfere with daily activities (eating, drinking, nest-building, locomotor). Additionally, none of the mice experienced bowel obstruction or changes in stooling or stool consistency. There was no evidence of blood associated with stooling. There was no evidence of anemia, i.e. loss of pink condition of mouse footpads. There was no evidence of compromised behavior, i.e. mice did not become lethargic and there was no incidence of piloerection, poor grooming or inability to thermoregulate. The mice experienced no sustained weight loss. Comparison survival curves for Hltf-deleted mice (n = 20) and control mice (n = 12) mice with the logrank (Mantel-Cox) test (Chi square 1.898, p = 0.1683), and the Gehan-Breslow-Wilcoxon test (Chi square 1.895, p = 0.1687), indicated Hltf-deletion had no significant effect on the mortality of the mice during the 35-day post-surgery timeline (Fig 1). However, three Hltf-deleted mice died. Thus, the hazard ratio (Mantel-Haenszel) indicated Hltf-deleted mice at any time during the treatment protocol were 5-times more likely to die than control mice. This approach determined an earlier humane endpoint of 28 days for Hltf-deleted mice in future experiments. The mice were treated alike and maintained under identical conditions.

[Figure omitted. See PDF.]

Fig 1. Kaplan-Meier survival plot.

Comparison of the cumulative survival curve of male Hltf-deleted (n = 20) and control (n = 12) mice shows the total probability of the mice surviving to the end of the study, i.e. 35 days post-surgery, was unaffected by genotype.

https://doi.org/10.1371/journal.pone.0251132.g001

HLTF expression is negatively associated with survival in human CRC

There is a negative correlation between HLTF expression and the progression of CRC in humans [11–16]. Previous studies have shown the progressive loss of HLTF from tumor cells correlates with negligible HLTF expression in the TME [19, 20]. The CDX model is predicated on the paucity of HLTF expression in the fibroblasts of the TME as shown in a well-differentiated adenocarcinoma (Fig 2) from a male patient that resulted in hemicolectomy.

[Figure omitted. See PDF.]

Fig 2. HLTF protein in a pT3 adenocarcinoma with no regional lymph node metastasis (pN0).

HLTF-immunostaining in well-differentiated (G1) cancer cells of an adenocarcinoma that arose from a tubulovillous adenoma (polyp). Tumor cells invaded through the muscularis propria into the peri-colorectal tissue at the hepatic flexure. There was no evidence of lymphovascular invasion. The photomicrographs in panels A-C show the cellular architecture of the tumor tissue composed of HLTF-positive cancer cells and HLTF-negative fibroblasts. A (4X magnification), B (10X magnification) and C (20X magnification). Antibodies to the HLTF N-terminus (residues 164–300) that are common to all known human HLTF proteins [21] in nuclear and cytoplasmic locations show the positive immunostain (brown) is predominately although not exclusively cytoplasmic. HLTF-negative fibroblasts (blue arrows) are devoid of nuclear and/or cytoplasmic staining.

https://doi.org/10.1371/journal.pone.0251132.g002

Hltf-deletion from the TME promotes metastasis in the CDX model

The goal of the first experiment was to establish primary tumor xenografts with HLTF expressing HCT116 Red-FLuc cells (easily detected by bioluminescence imaging—BLI) in control male mice by direct OCMI between the mucosa and the muscularis layers of the cecal wall (Fig 3). Care was taken to avoid microinjection into the lymphoid nodules in the distal part of the mouse cecum. Histopathological confirmation of tumor foci in the most clinically relevant [32] metastatic sites (e.g., colon-draining lymphatics, liver, lung and peritoneum) coincided with BLI. The goal of the second experiment was to establish primary tumor xenografts via OCMI in Hltf-deleted mice. Tumor take rates were 100% for primary xenograft tumors for Hltf-deleted (n-20) and control (n = 12) mice. However, a shift in the pattern of metastases (Fig 4) to regional lymph nodes, peritoneal carcinomatosis (PC), inguinal canal in close proximity to seminal vesicles and into the scrotum occurred in Hltf-deleted TME (Fig 5). PC is a late stage manifestation of CRC and considered lethal in humans [34].

[Figure omitted. See PDF.]

Fig 3. Direct orthotopic cell microinjections (OCMI).

Control mice received direct OCMI of 2x106 Hltf+/+ HCT116 Red-FLuc cells/10 μl into the submucosa of the variable size (3–4 cm) J-shaped cecum (A). In vivo BLI (one-second exposure) followed IP injection of D-luciferin (B). A 50 μl Hamilton syringe with 30GA/0.5 inch/30-degree custom needles eliminated retrograde leakage. Cell signal authenticity (C) was confirmed in vitro (cm = culture medium, luc = luciferin). BLI intensity (D) ranged from red (highest), through yellow, green, blue, and purple (lowest). Histopathological confirmation with Hematoxylin and Eosin staining of tumor foci in the most clinically relevant metastatic sites coincides with BLI (E) in control mice. BLI and histopathology showing progression in the same mouse.

https://doi.org/10.1371/journal.pone.0251132.g003

[Figure omitted. See PDF.]

Fig 4. Representative images of immunodetection, BLI, routine histochemistry and HLTF-immunohistochemistry.

Full-size HLTF protein (~115-kDa) was detected by Western blot of immunoprecipitated proteins from HCT116 Red-FLuc cells (A). The HLTF antibodies (residues 600–700) recognize all known and putative truncated mouse proteins [35]. BLI showed the metastatic properties (local invasion and distant colony formation) of primary tumor xenografts dramatically increased in an Hltf-deleted microenvironment (A). At necropsy, primary xenograft CDX tumors (B) were removed for routine histopathology (Hematoxylin and Eosin staining; 4X magnification). Primary xenograft CDX tumors (C) showed positive immunostaining for HLTF. Immunostaining in the Hltf+/+ TME contrasts with the Hltf-deleted TME (25 X magnification). The HLTF antibodies (residues 164–300) recognize all known human and mouse proteins [35].

https://doi.org/10.1371/journal.pone.0251132.g004

[Figure omitted. See PDF.]

Fig 5. Imaging of CDX metastasis in Hltf-deleted mice at necropsy.

Mice were anesthetized with isoflurane and metastatic CRC was detected by BLI (A) prior to necropsy. Immediately following euthanasia, mice were opened via a mid-sagittal incision and the primary CDX tumor (B1) was removed followed by lymph nodes near the stomach and spleen (B2). Metastatic sternal lymph nodes below (B3) and above (B4) the diaphragm were removed as was a peritoneal carcinomatosis (B5). Mice were imaged a second time (C) to guide removal of inguinal lymph nodes (D6) and metastatic tissue extending into the scrotum (D7). All tissues were evaluated for S-nitrosylation of human proteins.

https://doi.org/10.1371/journal.pone.0251132.g005

The dynamic interaction between cells of the tumor and the TME was decoded by species-specific comparative transcriptomics

Alternative splicing increases the complexity of HLTF gene expression in human and mouse cells. Because of its putative role in tumor development, alternative splicing in cells of the tumor and TME were evaluated. Two transcript variants encode the same HLTF protein in Homo sapiens [21]. Cuff.diff alternative splicing analysis identified transcript variant 1 (NM_003071.4), the longer transcript variant (5320-bp mRNA), as the only protein (Fig 4A) encoding mRNA transcript in primary CDX tumor cells. Hltf is alternatively spliced in mouse tissues [33, 35–37], and Cuff.diff alternative splicing analysis quantified the usage of each exon and each possible splice junction for Hltf in RNA-seq samples from control mouse TME. A 4:1 ratio of full-length 4956 nucleotide (nt) message isoform (NM_09210) to a long non-coding (Lnc) 3643 nt transcript variant 4 (Lnc-Hltf-4; NR_105047) was identified. This implicates Lnc-Hltf-4 in CRC progression.

Cuff.diff FPKM tracking files were analyzed with iPathwayGuide (Advaita Bioinformatics)–using the q value of 0.05 for statistical significance and a log fold change of expression with an absolute value of at least 0.06–30 differentially expressed (DE) genes were identified out of a total of 15,857 genes with measured expression in the mouse TME (S1 Table). Similarly, for human tumor analysis–using the q value of 0.05 for statistical significance and a log fold change of expression with an absolute value of at least 0.06–151 DE genes were identified out of a total of 12,358 genes with measured expression in the human tumor (S2 Table).

Hypoxia is an important regulator in CRC, and Hltf is a transcriptional regulator of murine hypoxia inducible factor-1α (Hif-1α) in heart [36]. However, transcriptional availability of HIF1α in the CDX tumors was unaffected by the presence or absence of Hltf in the TME. In addition, of the 151 DE genes, only 14 of ~2450 TRANSFAC predicted HIF1α targets were either up (n = 3) or down (n = 11) regulated, and no functional links between them were identified by iPathwayGuide. These findings argue Hltf-deletion from the TME had little or no effect on the oxygenated/metabolic state of CDX tumors. In fact, cytokine-receptor interaction (KEGG: 04060) was the top biological pathway in CDX tumors (p = 0.002) and mouse TME (p = 2.546e-5) with Bonferroni corrected p-values. Increased mRNA (Fig 6A) and protein (Fig 6B and 6C) expression of the pro-inflammatory cytokines, interleukin 8 (IL8) and intercellular adhesion molecule-1 (ICAM-1) was identified in CDX tumors with an Hltf-deleted TME. Transactivation of interleukin 6 (IL6) and all components of iNOS-S100A8/A9 signaling in the Hltf-deleted TME (Fig 6D) coincided with dramatic increases in the mRNA abundance of haptoglobin (HP) and SERPINA3, and decreased hepatocyte growth factor (HGF) in CDX tumors (Fig 6E). Compartmentalized expression of all components of the iNOS-S100A8/A9 signaling pathway in the Hltf-deleted TME are shown in Fig 7.

[Figure omitted. See PDF.]

Fig 6. Gene measured expression bar plots and confirmatory immunohistochemistry.

All DE genes are ranked based on their absolute log fold change. Upregulated genes are shown in red and down regulated genes are shown in blue. Box and whisker plots on the left of each histogram summarize the distribution of all DE genes in a specific pathway or annotated to a GO term out of all of the target genes. The box represents the first quartile, the median and the third quartile. Increased mRNA (A) and protein expression of pro-inflammatory cytokines, IL8 (B) and ICAM-1 (C), occurs in human tumors with an Hltf-deleted TME (40X magnification). There is no IL8 staining in the TME because IL8 is unique to the human genome, i.e. there is no IL8 homologue in the mouse. The log fold increase in message for IL6 and all members of the iNOS-S100A8/A9 signaling axis was documented in the mouse Hltf-deleted TME (D). Log fold changes in message for HP (increase, p = 0.01), SERPINA3 (increase, p = 0.01) and HGF (decrease, p = 0.043) in CDX tumors from Hltf-deleted TME are highlighted with yellow dots in a volcano plot (E). The horizontal axis is the log fold change, and the vertical axis is the negative base-10 logarithm of the p-value. The red-dotted lines represent the threshold. The upper regulated genes (positive long fold change) are shown in red, while the down-regulated genes are blue. Unaffected genes are black.

https://doi.org/10.1371/journal.pone.0251132.g006

[Figure omitted. See PDF.]

Fig 7. Immunohistochemistry for the iNOS-S100A8/A9 signaling axis.

Increased S100A8/A9 was exclusive to the Hltf-deleted mouse TME. S100A8 and S100A9 typically form heterodimers (S100A8/A9, calprotectin) and homodimers are not generally detectable. Increased S100A8/A9 expression coincides with increased iNOS protein levels in the Hltf-deleted mouse TME. Increased endogenous iNOS expression, a major mediator of inflammation, correlates with poor patient survival in CRC due to increased metastasis.

https://doi.org/10.1371/journal.pone.0251132.g007

iNOS-S100A8/A9 targets in CDX tumors

We tested the hypothesis that iNOS derived nitric oxide (NO) production in the Hltf-deleted TME promoted metastasis in the CDX model via S-nitrosylation. A snapshot of the entire CDX tumor S-nitroso-proteome in the presence/absence of Hltf in the TME was obtained when iodoTMT enriched S-nitrosylated human proteins were interrogated by nanoLC-MS/MS. NanoLC-MS/MS-based protein identification provided a comprehensive collection of 136 SNO-proteins in CDX tumors with an Hltf-deleted TME compared to 178 SNO-proteins in CDX tumors with control TME (S1 File). Protein-protein interaction analysis with PANTHER designated cytoskeletal regulation by the Rho GTPase (P00016) the major pathway for CDX tumors from both Hltf-deleted and control TME. However, a secondary pathway of inflammation-mediated by chemokine and cytokine signaling (P00031) increased from 11 protein candidates to 18 protein candidates in CDX tumors with the Hltf-deleted TME. The emPAI value (0.28) for POTEE, the primate-specific POTE (prostate, ovary, testis and embryo expressed) ankyrin domain family member E that promotes CRC growth [38], indicates this protein is a general target of S-nitrosylation in CDX tumors from Hltf-deleted and control TME.

Differential expression of TMT-labeled proteins was studied by 2DE with Cy5 label followed by 2D-DIGE (S2 File) and MALDI-TOF/TOF (S3 File). IodoTMT-labeled proteins from Hltf-deleted and control mice were covalently tagged with either Cy3 (Hltf-deleted) or Cy2 (control) for 2D-DIGE (Fig 8). A gel with three fluorescent images (Cy5/Cy2/Cy3) was processed for spot-codetection and protein identification by MALDI-TOF/TOF with a standard database search for human proteins with MASCOT. Proteins-peptides with a total ion C.I.% = 100 are presented in the protein ID summary.

[Figure omitted. See PDF.]

Fig 8. 2D-DIGE gel electrophoresis.

For purposes of data presentation, only select spots are indicated in the gel. The pH range 4 to 9 is indicated on the x-axis (red horizontal line). Molecular mass (20–180 KD) is indicated on the y-axis. Differentially expressed proteins were subjected to MALDI-TOF/TOF analysis.

https://doi.org/10.1371/journal.pone.0251132.g008

2D-DIGE and protein identification by MALDI-TOF/TOF led to the authentication of iNOS-S100A8/A9 consensus site-specific S-nitrosylated proteins in CDX tumors from Hltf-deleted TME. As shown in Fig 8, spot 21 that is unique to CDX from Hltf-deleted TME has authentic S-nitrosylation of Cys638 in the consensus sequence motif (I/L-X-C-X-X-D/E) of three members (E/F/I) of the primate-specific POTE gene family comprised of 14 members and subdivided into four groups. Group 3 POTE-actin genes (POTEs E/F/I/J/KP) encode proteins that are actin chimeras with a full-sized long inverted repeat, seven ankyrin repeats and a C-terminal coiled-coil domain. The consensus iNOS-S100A8/A9 target sequence (aa 636–641) immediately precedes the start of the coiled-coil domain (aa 642–698). POTEE-SNO in CDX tumors is both a general S-nitrosylation target and an iNOS-S100A8/A9 site-specific target in the Hltf-deleted TME. When additional very weak spots (30–35) were added to the analysis, authentic S-nitrosylation of Cys238 in the consensus sequence motif (I/L-X-C-X-X-D/E) was identified in the tripartite motif (TRIM) family member TRIM52. Dysregulation of TRIM family members characterized by a tripartite motif–RING domain, one or two B-box domains, and a coiled-coil domain–has been implicated in CRC cell proliferation [39]. This newly identified iNOS-S100A8/A9 target sequence (aa 236–241) is located in the B-box domain. In contrast, all members of the ERM family—ezrin, radixin, and moesin—were found in spot 31 as general targets of S-nitrosylation (Cys284). There was no evidence of iNOS-S100A8/A9 site-specific S-nitrosylation (Cys117). 2D-DIGE and protein identification by MALDI-TOF/TOF led to the authentication of S-nitrosylation of 53 individual cysteines in half-site motifs of either I/L-X-C or C-X-X-D/E in proteins from CDX tumors from both Hltf-deleted and control TME (Table 1).

[Figure omitted. See PDF.]

Table 1. S-nitrosylated cysteine residues (yellow) in half-site motifs of either I/L-X-C or C-X-X-D/E from CDX tumors from Hltf-deleted (white) and control (green) TME.

https://doi.org/10.1371/journal.pone.0251132.t001

Continuity between primary CDX tumors and metastatic sites

Comparative assessment of S-nitrosylation in primary CDX tumors and five metastatic sites was undertaken with iodoTMT-switch labeling and nanoLC-MS/MS. As shown in Fig 5, BLI was used to carefully isolate tumors and lymph nodes, and proteins were selectively labeled with iodoTMT isobaric label reagents as follows: iodoTMT-126, lymph nodes near stomach and spleen; iodoTMT-127, primary tumors; iodoTMT-128, peritoneal carcinomatosis; iodoTMT-129, tumors at the sternum (below and above the diaphragm); iodo-TMT130 lymph nodes in the inguinal region; and iodoTMT-131, lymph nodes extending into the scrotum.

As shown in S4 File, 331 proteins were identified in primary CDX tumors and distant metastatic sites. Of that total, 60% were identified in all distant metastatic sites with different ratios of expression compared to primary CDX tumors. In contrast, 40% of the proteins from primary CDX tumors were found in some but not all sites. Protein-protein interaction analysis amongst the shared SNO-proteins with PANTHER designated cytoskeletal regulation by the Rho GTPase (P00016) the major pathway for continuity between CDX tumors and metastatic tumors/lymph nodes in Hltf-deleted TME. The previously identified secondary pathway of inflammation-mediated by chemokine and cytokine signaling (P00031) now includes SNO-myosins and provides continuity between the primary CDX tumor and all distant locations. Functional enrichment assessment with STRING included actin filament-based process, movement, and binding in conjunction with muscle filament sliding and contraction (Fig 9). These findings strongly implicate S-nitrosylation of actins and myosins under physiological conditions in the metastatic dissemination of CDX tumors in an Hltf-deleted TME. Additionally, S-nitrosylation of Cys813 in UNC45 myosin chaperone B protein was identified in an iNOS-S1008A/A9-specific motif. UNC45B, a co-chaperone for HSP90 required for folding and accumulation of type II myosins [40] has three tetratricopeptide repeats followed by three armadillo (ARM) repeats. ARM 3 at position 751–790 ends prior to the solitary SNO site at Cys813. This SNO-protein is present in all distant metastatic sites except the inguinal region. The NF-kB subunit REL is an example of convergence of transcriptomic-S-nitroso-proteomic signaling. REL is transcriptionally activated (log fold increase 1.768, p = 0.01) in primary CDX tumors from an Hltf-deleted TME, and the encoded protein is an S-nitrosylation target (Cys143) in all metastatic sites. However, expression is increased (14-fold) in tissue from the inguinal region compared with the primary CDX tumor. Vinculin—a general adhesion protein—dramatically distinguished the metastatic peritoneal carcinomatosis (39-fold) and the lymph nodes extending into the scrotum (168-fold) from the other metastatic sites. However, vinculin appears to be a co-immunoprecipitate because we were unable to authenticate an S-nitrosylated cysteine residue.

[Figure omitted. See PDF.]

Fig 9. Protein-protein interaction (PPI) networks.

Cluster analysis of the proteins in the primary tumors and all metastatic sites had a PPI enrichment p-value of <1.0e-16 indicating the shared SNO-proteins are biologically connected as a group. The reliability of the putative interactions was increased by setting the confidence score at the highest level of 0.900 with a medium (5%) setting for FDR stringency. For secondary analysis, only databased sourced interactions were selected for the network, and all disconnected nodes in the network were eliminated. GO term affiliated with this analysis are as follows. GO:0030049 (orange) muscle filament sliding (1.29e-09), GO:0007018 (blue) microtubule-band movement (5.15e-08); GO:00030048 (yellow) actin filament-based movement (7.48e-10); GO:0006928 (green) movement of cell or subcellular component (7.73e-05).

https://doi.org/10.1371/journal.pone.0251132.g009

Discussion

The CDX model—male HCT-116 cells in male mice—is predicated on the fact that male gender is a risk factor in metastatic CRC [41]. General risk factors include age, family history of CRC or inflammatory bowel disease, alcoholism, smoking, obesity, and diabetes. Regardless of lifestyle-specific risks, at any age, the incidence of CRC is higher in men [42]. Hormones in the form of oral contraception and hormone replacement therapy protect women from CRC [43]. There is gender disparity in the anatomical location for CRC. Men have a higher risk of developing left-sided (distal to splenic flexure) CRC compared to women who have a higher risk of developing right-sided (proximal to splenic flexure) CRC [44]. Epigenetic silencing of HLTF occurs more often in men [18]. Experimental models such as this CDX model are invaluable as a CRC metastasis model due to the small sample size and sample amount of relevant human tissue, and the limited availability of matched primary and metastatic samples from individual patients. Not to mention the paucity of normal control tissue.

RNAseq and transcriptomic analyses have identified LncRNAs [45, 46] that are gaining importance in cancer genomics studies [47]. LncRNAs—nonprotein coding transcripts >200 nucleotides—are evolutionarily conserved and have been implicated in regulatory specificity through interactions with proteins and/or RNA. LNCipedia (v 5.2) lists 10 Lnc-HLTF transcripts of varying lengths. Lnc-HLTF-1, Lnc-HLTF-2 and Lnc-HLTF-3 are cancer-related although no function is currently assigned to them. In contrast, Lnc-HLTF-5 is elevated in thoracic aorta tissue of patients with hypertension, and positively correlated with expanding ascending aortic diameter [48]. This study is the first to report murine Lnc-Hltf-4 in the Hltf+/+ TME of a CDX tumor model. It remains to be determined if silencing Lnc-Hltf-4 expression in the Hltf-deleted TME actively promoted metastasis.

In this study, the proinflammatory mediator S100A8/A9 in combination with species-specific expression of proinflammatory cytokines IL6 (mouse) and IL8 (human) established a pro-metastatic primary tumor niche in the Hltf-deleted TME. Changes in CDX gene expression further promoted metastasis as ICAM-1 increases cancer cell invasion/intravasation into the microvasculature [49] and mediates peritoneal carcinomatosis [50] HP promotes colorectal cancer cell motility [51] and SERPINA3 promotes tumor cell migration and invasion [52]. Decreased HGF may be partly responsible for reduced liver metastasis in this model [53]. Combinatorial induction of IL6 and S100A8/A9 in the TME and ICAM-1 and IL8 in the CDX tumor comprise a positive feedback loop that drives inflammation [54]. These findings support the development of nano-therapeutic tumor-specific and TME-targeted anti-inflammatory therapy for CRC [55, 56].

Increased HP known to preserve vascular NO signaling [57] occurred in conjunction with increased iNOS in the Hltf-deleted TME. INOS contributed to the metastatic phenotype via S-nitrosylation of cytoskeletal target proteins (actin, myosin). One of the newest components of the TME under investigation is the mechanical microenvironment [58] with an emphasis on actin [59] especially in invadopodia formation [60] where S-nitrosylation likely plays a regulatory role. The metastatic phenotype was promoted via iNOS-S100A8/A9 site-specific S-nitrosylation of UNC45 myosin chaperone B protein that engages in non-muscle myosin (myosin II) assembly [40]. Most intriguing is the enrichment of vinculin in the two sites of peritoneal metastasis [34]—metastatic peritoneal carcinomatosis [61, 62] and testicular/scrotal metastasis [63]—considered rare and life-threatening. Vinculin, a coimmunoprecipitate—an intracellular F-actin-binding protein and a mechanotransducer in cell-cell and cell-matrix adhesions [64, 65]—has been implicated in tumor progression [66].

Functional protein posttranslational modifications include phosphorylation, ubiquitination and S-nitrosylation. Despite the ongoing search for S-nitrosylation motifs [67], iNOS-S100A8/A9 signaling is the only known process of S-nitrosylation that is site-selective [26]. There are more than 100 putative targets of iNOS-S100A8/A9 nitrosylation. Of those targets, the five with experimentally validated consensus motifs are Annexin A5, Ezrin, GAPDH, Moesin and Vimentin [26]. As a result of this study, POTEE, TRIM52 and UN45B can be added to the list of proteins with authentic SNO sites in the conserved I/L-X-C-X2-D/E motif indicative of iNOS-S100A8/A9-mediated S-nitrosylation. Moreover, it appears there is iNOS-S100A8/A9 half-site selective S-nitrosylation. Jai et al. [26] authenticated the iNOS-S100A8/A9 consensus site with in vitro mutagenesis experiments in which the target cysteine was changed, or knockdown experiments (siRNA) in which individual components of the heteroduplex (S100A8/A9) were eliminated. However, in this in vivo situation, in the presence of both S100A8 and S100A9 proteins, half-site S-nitrosylation was authenticated. Global S-nitrosylation has a role in the regulation of gene transcription [68]. Zinc-fingers motifs are labile to NO where the thio-ligands are vulnerable to S-nitrosylation that can result in zinc release and collapse of the ring-motif resulting in loss of function [69]. HLTF is a zinc finger protein with a highly conserved C3HC4 ring motif (aa760-801); however, we found no evidence of human HLTF-SNO in primary tumors from Hltf+/+TME. Collectively, these findings support ongoing efforts to find selective iNOS inhibitors as chemo-preventive agents against CRC [23].

Conclusions

In this study, we show for the first time that Hltf-deletion from the TME promotes inflammation in the shared TME-primary tumor niche. Induction of iNOS in the TME produced general and iNOS-S100A8/A9 site-specific S-nitrosylation of previously unidentified human tumor proteins. We establish continuity between global S-nitrosylation of proteins in pathways related to cytoskeleton and motility under physiological conditions in the metastatic dissemination of CDX tumors in an Hltf-deleted TME. We provide the first evidence of cross-talk between increased gene transcription and S-nitrosylation of the encoded protein. This new role for Hltf-deletion in NO-mediated protein S-nitrosylation to promote metastasis extends our understanding of Hltf as a tumor suppressor.

Materials and methods

Reagents and kits

Abcam (Cambridge, MA) was the source of the following antibodies: mouse monoclonal (ab22506) to S100A9 + Calprotectin (S100A8/A9 complex), mouse monoclonal (ab18672) to IL8, rabbit polyclonal anti-iNOS (ab15323), rabbit monoclonal anti-ICAM-1 (ab109361), goat anti-rabbit IgG H&L (HRP; ab97051), goat anti-mouse (HRP; ab205719); and mouse on mouse polymer IHC kit (ab 127055). For immunoprecipitation and Western blotting, Abcam was the source of rabbit polyclonal anti-HLTF (ab17984) to human HLTF aa 600–700 that reacts with mouse and human. For immunohistochemistry, Sigma (St Louis, MO) was the source of rabbit polyclonal anti-HLTF (HPA015284) to human HLTF aa 164–300 that reacts with mouse and human. Biotinylated goat-anti-mouse (BA-9200) and goat anti-rabbit (BA-1000) IgG antibodies and the ABC-enzyme complex were purchased from Vector Laboratories, Inc. (Burlingame, CA). Harris Modified Hematoxylin (HHS16) was purchased from MilliporeSigma (St. Louis, MO). XenoLight D-luciferin potassium salt (122799) was purchased from PerkinElmer (Waltham, MA). ThermoScientific (Waltham, MA) was the source of the following Pierce™ reagents: S-nitrosylation Western blot kit (90105), HENS buffer (90106), mouse anti-TMTantibody (90075), immobilized anti-TMT resin (90076), TMT elution buffer (90104) and, iodoTMTsixplex™ label reagent set (90101). BioRad (Hercules, CA) was the source of 7.5% Mini-PROTEAN TGX precast protein gels (4561024), Clarity Western ECL-substrate (170–5060), Precision Protein StrepTactin-HRP conjugate (1610381) and Kaleidoscope SDS-PAGE Standards (1610324). PerkinElmer was the source of Bioware® Brite Cell Line HCT116 Red-FLuc (BW124318). All protocols are accessible in protocols.io (dx.doi.org/10.17504/protocols.io.bs5xng7n).

Cell culture

Bioware® Brite Cell Line HCT116 Red-FLuc is a cell line derived from the parental cell line (ATCC, CCL-247) from adult male colorectal carcinoma by stable transduction with red-shifted lentivirus containing firefly luciferase from Luciola Italica (Red-FLuc) under the control of human ubiquitin C promoter. HCT116 Red-FLuc cells were confirmed to be pathogen free by the IMPACT Profile I (PCR) at the University of Missouri Research Animal Diagnostic and Investigative Laboratory. HCT116 Red-FLuc cells have a mutation in codon 13 of the ras proto-oncogene, and express transforming growth factor beta 1 (TGFβ1) and tumor protein 53 (TP53), HCT116 Red-FLuc cells grown in McCoy’s 5a Modified Medium (ATCC, 30–2007) supplemented with 10% fetal bovine serum and puromycin (2μg/mL) have an average doubling time of 16 hours. For each experiment, cell stocks in liquid nitrogen at passage 2 were thawed using T25 flasks, expanded in T150 flasks for 2 days, passaged overnight and harvested at 70–75% confluency. This protocol provided a unified framework for comparative gene expression analysis.

Hltf-deleted and control mice

Global Hltf-deleted mice were developed in collaboration with genOway (Lyon, France) as previously described [35], and bred to be fully congenic (N11) on the C57BL/6J genomic background [36]. Global Hltf-deleted mice present a neonatal lethal phenotype. However, when the Hltf-deletion line was bred (IACUC# 02007) into the recombinase activating gene 2 (Rag2)/common gamma (IL2rg) double knockout background [37], i.e. mice lacking lymphocytes (NK-, T- B-; alymphoid), the perinatal lethal phenotype was eliminated—clearly showing the perinatal lethal phenotype requires an immune component.

Immune-deficient mice were housed with a 12:12 light/dark cycle with access to food and water ad libitum and bedding was changed 2–3 times/week. Routine testing of sentinel mice ensured the colony was disease free. All studies and the anticipated mortality were conducted in accord with the NIH Guidelines for the Care and Use of Laboratory Animals, as reviewed and approved by the Animal Care and Use Committee at Texas Tech University Health Sciences Center (NIH Assurance of Compliance A3056-01; USDA Certification 74-R-0050, Customer 1481, S1 Checklist). TTUHSC’s IACUC (# 02009) specifically approved this study.

The orthotopic HLTF+/+HCT116 xenograft model was established as follows: randomly selected six- to eight-week old Hltf-deleted (n = 20) and HLTF+/+ (n = 12) male Rag2-/-IL2rg-/- mice received direct orthotopic cell microinjections (OCMI) of HLTF+/+HCT116 Red-Fluc cells (2x106 cells/10 μl) between the mucosa and the muscularis layers of the cecal wall. Hereafter the mice were designated Hltf-deleted and control. All surgery was performed with isoflurane (Isothesia) and the SomnoSuite® Low-Flow anesthesia system (Kent Scientific) with far infrared warming pads during surgery and recovery. Additional efforts to minimize suffering included an IP injection of Buprenorphine (Buprenex, 0.1 mg/kg) prior to surgery to manage incisional pain followed by a second dose 4–8 hours later. The cecum was exteriorized via a small midline laparotomy on the vertical linea alba to eliminate bleeding. Non-invasive bioluminescence imaging (BLI) with an IVIS Spectrum In Vivo Imaging System was used to validate the quality and accuracy of the injection, and to track and quantify tumor growth and metastasis. Histopathology at necropsy confirmed placement of the inoculum. Mouse behavior and well-being were monitored daily. Tumor growth/metastasis was monitored weekly with BLI.

Techniques

Postmortem analysis.

Necropsy was performed at 35 days post CDX establishment. Mice under continuous isothesia were imaged and killed immediately (< 15 seconds) by cervical dislocation. Primary tumor xenografts and metastatic tumors were quickly removed, rinsed in physiological saline and either flash frozen for biochemical evaluation (RNA-seq, Western blotting, MALDI-TOF/TOF MS, nanoLC-MS/MS) or fixed in formalin and processed for either routine histopathology or immunohistochemistry.

Immunohistochemistry.

Tissue blocks were serially sectioned (3–4 μm). Two sections were placed on each slide and deparaffinized prior to staining. Beginning with the first slide, sections on every fifth slide were stained with hematoxylin and eosin (H&E) for evaluation by light microscopy. Sections on alternate slides were processed for immunohistochemistry with heat-induced epitope retrieval. Two tissue sections per slide facilitated the use of one section for positive immunostaining, and the serial section for negative (minus primary antibody) control staining. All primary antibodies (1:50) were paired with an appropriate HRP-conjugated secondary antibody (1:200) depending upon the species in which the primary antibody was generated. Nuclei were counterstained (blue) with hematoxylin.

Immunoprecipitation and Western blotting.

Immunoprecipitation and Western analysis were performed as previously reported [33, 35]. Briefly, whole cell lysates from two T150 flasks of HCT116 Red-FLuc cells (80% confluent) were immunoprecipitated with rabbit polyclonal anti-HLTF (ab17984) to human HLTF aa 600–700 at a concentration of 5 ug/ml. Western blotting was achieved with the same primary antibody (1:5000) followed by HRP-conjugated mouse anti-rabbit (1:5000). Signal was detected by chemiluminescence with the Clarity Western ECL Substrate Kit.

Tumor transcriptome analysis (RNA-seq)

Primary tumor xenografts (1 per individual mouse x 3 biological replicates for Hltf-deleted and control male mice = 6 total samples) were flash frozen and sent to Otogenetics Corp. (Norcross, GA) for RNA-seq assays as previously described [33, 35–37]. Briefly, total RNA was isolated, and evaluated for its integrity and purity with an Agilent Bioanalyzer (Table 2). RNA samples were rRNA-depleted prior to random-primed cDNA preparation/QC, Illumina library preparation/QC; PE100-125 and HiSeq2500 sequencing at a minimum of 40 million reads. Paired-end 106 nucleotide reads were subjected to species-specific mapping against reference genomes for mouse (mm10) and human (hg38). Differential expression (DE) analysis and alternative splicing analysis for each species was done with cufflinks.cuffdiff (2.2.1). Data were imported into iPathwayGuide (Advaita Corporation, Plymouth, MI) and evaluated using the q-value of 0.05 for statistical significance and a log-fold change (logFC) of expression with an absolute value ≤ 0.6. Cuffdiff generated q-values are adjusted p-values that consider the false discovery rate (FDR). The q-value is an essential statistic when measuring thousands of gene expression levels from a relatively small sample set because it has a greater ability (power) to identify significant changes in gene expression. All RNA-seq data in this publication are accessible through NCBI’s Gene Expression Omnibus (GEO) Series accession number GSE161961 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE161961).

[Figure omitted. See PDF.]

Table 2. Sample quality control and RNA-seq outcome.

https://doi.org/10.1371/journal.pone.0251132.t002

With iPathwayGuide, data were analyzed in the context of pathways obtained from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (Release 90.0+/05-29, May 19), gene ontologies from the Gene Ontology Consortium database (2019-Apr26), miRNAs from the miRBase (Release 22.1, October 2018) and TARGETSCAN (Targetscan version: Mouse:7.2, Human:7.2) databases, network of regulatory relations from BioGRID: Biological General Repository for Interaction Datasets v3.5.171. March 25th, 2019, and diseases from the KEGG database (Release 90.0+/05-29, May 19).

Tumor S-nitroso-proteome analysis

Step-wise tumor S-nitroso-proteomic analysis leading to iNOS-S100A8/A9 site-specific analysis was performed with iodoTMT-switch labeling [28], affinity enrichment and high-resolution LC-MS/MS analysis by Applied Biomics, Inc (Hayward, CA). Briefly, primary tumor xenografts from Hltf-deleted (n = 2 mice) and control (n = 2 mice) male mice were homogenized/sonicated (Polytron) in 4 volumes HENS buffer (100mM HEPES, pH 7.8, 1 mM EDTA, 0.1 mM Neocuproine, 1% SDS). Protein concentrations were determined (OD280) with NanoDrop One (ThermoFisher), and adjusted to a final concentration of 2 μg/μl with HENS buffer. Two hundred microgram samples were incubated for 30 minutes at room temperature with 20 mM methyl methanethiosulfonate (MMTS) in dimethylformamide (DMF) to block free cysteine thiols. Proteins were precipitated with 6 volumes pre-chilled (-20°C) acetone for a minimum of 60 minutes to remove MMTS, pelleted by centrifugation (10,000 x g) for 10 minutes at 4°C, and dried for 10 minutes. Following sample resuspension in HENS buffer (100 μl), S-nitrosylated cysteines were selectively reduced with ascorbate (protected from light) and irreversibly labeled with iodoTMTzero reagent for 2 hours at room temperature. Proteins were precipitated with 6 volumes pre-chilled (-20°C) acetone for a minimum of 60 minutes, pelleted by centrifugation (10,000 x g) for 10 minutes at 4°C, and dried for 10 minutes. Samples were resuspended in HENS buffer (100 μl), added to anti-TMT resin and incubated overnight with end-over-end mixing at 4°C. Resin was washed 3X with 1XTBS, 3X with water, and eluted with 4 volumes TMT elution buffer. Eluates were frozen and lyophilized to dryness. In experiment 1, global S-nitrosylation of CDX tumors from TMT-tagged Hltf-deleted and control samples was analyzed by NanoLC-MS/MS with a data base search for human genes in SwissProt using MASCOT. Exponentially Modified Protein Abundance Index (emPAI defined as 10PAI-1) values, where PAI (Protein Abundance Index) is the ratio of observed to observable peptides, is approximately proportional to the logarithm of protein concentration and indicates absolute protein abundance.

In experiment 2, Hltf-deleted and control TMT-tagged samples were resolved by two-dimensional gel electrophoresis (2DE) with mouse anti-TMT as the primary antibody (1 μg/ml final concentration) and Cy5-labeled donkey anti-mouse IgG (1:2000) as the second antibody. In experiment 3, TMT-labeled proteins were co-labeled with different color fluorescent dyes, i.e. Cy3 for Hltf-deleted and Cy2 for control for two-dimensional difference gel electrophoresis (2D-DIGE), which detected as little as 1.0 fmol of protein in each sample. Spots were excised and proteins identified by matrix-assisted laser-desorption ionization time-of-flight (MALDI-TOF/TOF). Database search was for human proteins in SWISSProt using MASCOT.

In experiment 4, unmodified cysteines in proteins (1 mg/ml) from primary tumor xenografts (n = 4 mice) and metastatic tumors (n = 2 mice) from Hltf-deleted mice were blocked with MMTS, selectively reduced with ascorbate (protected from light), and individually tagged with iodoTMTsixplex isobaric reagent such that a unique mass reporter (126–131 Da) in the low-mass region of the MS/MS spectrum was generated for samples from six locations. Equal amounts of six different samples (1 mg each) were combined into a single sample for enrichment with anti-TMT resin and NanoLC-MS/MS analysis. Relative expression for each protein fragment, i.e. the average ratio(s) for the protein, together with the number of peptide ratios that contributed (N) and the geometric standard deviation (SDgeo) were calculated.

In all experiments, peptides labeled with iodoTMT were quantified. Off-target (non-cysteine) labeling was <5%, and authentic cysteine nitrosylation as well as authentic iNOS/S100A8/A9 consensus sequences [26] were confirmed by peptide sequence. Gene symbols for differentially expressed proteins were input to the Protein ANalysis THrough Evolutionary Relationships (PANTHER) database [70] for functional classification and pathway identification. Protein-protein networks (cluster analysis) constructed with STRING [71] with EMPAI values identified direct (physical) and indirect (functional) interactions.

Supporting information

S1 Checklist. ARRIVE guidelines 2.0 (fillable).

https://doi.org/10.1371/journal.pone.0251132.s001

(PDF)

S1 Table. Met-analysis gene summary identified 30 DE genes in mouse TME.

Cuffdiff, part of Cufflinks, uses RPKM values to calculate changes in gene expression. Cuffdiff data were input into iPathwayGuide. Met-analysis calculated long fold change (logFC) and an adjusted p-value (adjpv) for each comparison. The minus sign in the column for logFC shows 12 genes were downregulated and 18 genes were upregulated.

https://doi.org/10.1371/journal.pone.0251132.s002

(XLSX)

S2 Table. Met-analysis gene summary identified 151 DE genes in CDX tumor cells.

Cuffdiff, part of Cufflinks, uses RPKM values to calculate changes in gene expression. Cuffdiff data were input into iPathwayGuide. Met-analysis calculated long fold change (logFC) and an adjusted p-value (adjpv) for each comparison. The minus sign in the column for logFC shows 136 genes were downregulated and 15 genes were upregulated.

https://doi.org/10.1371/journal.pone.0251132.s003

(CSV)

S1 File. Excel files for control, test (Hltf-deleted TME) and summary of the first experiment in tumor S-nitroso-proteome analysis.

Control and test files each contain three sections: section 1 is the legend (how to read the data sheets); section 2 is the protein peptide summary with peptide sequence, and section 3 is the list of significant hits (95% confidence). The summary file contains a comparative list of proteins with emPAI values for control or test or both.

https://doi.org/10.1371/journal.pone.0251132.s004

(ZIP)

S2 File. Data from the second and third experiment in tumor S-nitroso-proteome analysis.

Western blot analysis of test and control samples, 2D-DIGE analysis, and spot identification. Spots were excised and proteins identified by matrix-assisted laser-desorption ionization time-of-flight (MALDI-TOF/TOF).

https://doi.org/10.1371/journal.pone.0251132.s005

(ZIP)

S3 File. Excel files for spot analysis by MALDI-TOF/TOF.

Analysis was conducted in two phases in which spots with a unique or strong signal were analyzed first followed by analysis of spots with a very weak signal Two files each have a section on how to read MS report, a protein peptide summary and a protein ID summary. Only high confidence matches are reported.

https://doi.org/10.1371/journal.pone.0251132.s006

(ZIP)

S4 File. Excel files for proteins shared by primary and metastatic tumors.

There is a legend showing how to read the file, a peptide list and a protein list. The protein list contains ratio data show the relationship of protein expression in the secondary metastatic sites to the primary CDX tumor in the Hltf-deleted TME.

https://doi.org/10.1371/journal.pone.0251132.s007

(ZIP)

Citation: Helmer RA, Martinez-Zaguilan R, Kaur G, Smith LA, Dufour JM, Chilton BS (2021) Helicase-like transcription factor-deletion from the tumor microenvironment in a cell line-derived xenograft model of colorectal cancer reprogrammed the human transcriptome-S-nitroso-proteome to promote inflammation and redirect metastasis. PLoS ONE 16(5): e0251132. https://doi.org/10.1371/journal.pone.0251132

References

1. Seyfriend TN, Huysentruyt LC. On the origin of cancer metastasis. Crit Rev Oncog 2013;18(1–2):43–73. pmid:23237552

2. Keum N, Giovannucci E. Global burden of colorectal cancer: emerging trends, risk factors and prevention strategies. Nat Rev Gastroenterol Hepatol. 2019;16(12):713–732. pmid:31455888

3. Xu W, He Y, Wang Y, Li X, Young J, Loannidis JPA, et al. Risk factors and risk prediction models for colorectal cancer metastasis and recurrence: an umbrella review of systematic reviews and meta-analysis of observational studies. BCM Medicine 2020;18(172): pmid:32586325

4. Lichtenstern CR, Ngu RK, Salapour S, Karin M. Immunotherapy, inflammation and colorectal cancer. Cells. 2020;9(3):618–635

5. Keum N, Giovannucci E. Global burden of colorectal cancer: emerging trends, risk factors and prevention strategies. Nat Rev Gastroenterol Hepatol. 2019;16:713–732, pmid:31455888

6. Hishiki A, Hara K, Ikegaya Y, Yokoyama H, Shimizu , Sato M, et al. Structure of a novel DNA-binding domain of helicase-like transcription factor (HLTF) and its functional implication in DNA damage tolerance. J Biol Chem. 2015;290(21):13215–13223. pmid:25858588

7. Sheridan PL, Schorpp M, Voz ML, Jones KA. Cloning of an SNF2/SWI2-related protein that binds specifically to the SPH motifs of the SV40 enhancer and to the HIV-1 promoter. J Biol Chem. 1995;270(9):4575–4587. pmid:7876228

8. Hayward-Lester A, Hewetson A, Beale EG, Oefner PJ, Doris PA, Chilton BS. Cloning, characterization and steroid-dependent posttranslational processing of RUSH-1α and β, two uteroglobin promoter binding proteins. Mol Endocrinol 1996;10:1335–1349. pmid:8923460

9. Ding H, Descheemaeker K, Marynen P, Nelles L, Carvalho T, Carmo-Fonesca M, et al. Characterization of helicase-like transcription factor involved in the suppression of the human plasminogen activator inhibitor-1 gene. DNA Cell Biol. 1996;15:429–442. pmid:8672239

10. Kile AC, Chavez DA, Bacal J, Eldirany S, Korzhnev DM, Bezsonova I, et al. HLTF’s ancient HIRAN domain binds 3’ DNA ends to drive replication fork reversal. Mol Cell 2015;58(6):1090–1100. pmid:26051180

11. Unk I, Hajdu I, Fatyol K, Hurwitz J, Yoon JH, Prakash L, et al. Human HLTF functions as a ubiquitin ligase for proliferating cell nuclear antigen polyubiquitination. Proc Natl Acad Sci USA 2008;105(10):3768–3773. pmid:18316726

12. Seelinger M, Otterlei M. Helicase-like transcription factor HLTF and E3´ ubiquitin ligase SHPRH confer DNA damage tolerance through direct interactions with proliferating cell nuclear antigen (PCNA). Int J Mol Sci 2020;21(3):693–706.

13. Moinova HR, Chen WD, Shen L, Smiraglia D, Olechnowicz J, Ravi L, et al. HLTF gene silencing in human colon cancer. Proc Natl Acad Sci USA. 2002;99:4562–4567. pmid:11904375

14. Leung WK, To KF, Man EP, Chan MW, Hui AJ, Ng SS, et al. Detection of hypermethylated DNA or cyclooxygenase-2 messenger RNA in fecal samples of patients with colorectal cancer or polyps. 2007; Am J Gastroenterol. 102(5):1070–1076.

15. Wallner , Herbst A, Behrens A, Crispin A, Stieber P, Goke B, et al. Methylation of serum DNA is an independent prognostic marker in colorectal cancer. Clin Cancer Res. 2006;12(24):7347–7352. pmid:17189406

16. Herbst A, Wallner M, Rahmig K, Stieber P, Crispin A, Lamerz R, et al. Methylation of helicase-like transcription factor in serum of patients with colorectal cancer is an independent predictor of disease recurrence. Eur J Gastroenterol Hepatol. 2009;21(5):565–569. pmid:19282772

17. Philipp AB, Nagel D, Stieber P, Lamerz R, Thalhammer I, Herbst A, et al. Circulating cell-free methylated DNA and lactate dehydrogenase release in colorectal cancer. BMC Cancer. 2014;14:245–253 pmid:24708595

18. Philipp AB, Stieber P, Nagel D, Neumann J, et al. Prognostic role of methylated free circulating DNA in colorectal cancer. Int J Cancer 2012;131:2308–2319. pmid:22362391

19. Liu L, Liu H, Zhou Y, He J, Liu Q, Wang J, et al. HLTF suppresses the migration and invasion of colorectal cancer cells via TGF–β/SMAD signaling in vitro. Int J Oncol. 2018;53:2780–2788. pmid:30320371

20. Li H, Courtois ET, Sengupta D, Tan Y, Chen KH, Goh JJL, et al. Reference component analysis of single-cell transcriptomics elucidates cellular heterogeneity in human colorectal tumors. Nat Genet. 2017;49(5):708–718. pmid:28319088

21. Dhont L, Mascaux C, Belayew A. The helicase-like transcription factor (HLTF) in cancer: loss of function or oncomorphic conversion of a tumor suppressor? Cell Mol Life Sci. 2016;73(1):129–147. pmid:26472339

22. Somasundaram V, Basudhar D, Bharadwaj G, No JH, Ridnour LA, Cheng RYS, et al. Molecular mechanisms of nitric oxide in cancer progression, signal transduction, and metabolism. Antioxid Redox Signal. 2019;30(8):1124–1143. pmid:29634348

23. Wang H, Wang L, Xie Z, Zhou S, Li Y, Zhou Y, et al. Nitric oxide (NO) and NO synthases (NOS)-based targeted therapy for colon cancer. Cancers (Basel) 2020;12(7):1881–1905.

24. Zheng FV, Zheng X, Walia Y, Sharma V, Letson J, Furuta S. S-nitrosylation: an emerging paradigm of redox signaling. Antioxidants (Basel) 2019;8(9):404–435. pmid:31533268

25. Stomberski CT, Hess DT, Stamler JS. Protein S-nitrosylation: Determinants of specificity and enzymatic regulation pf S-nitrosothiol-based signaling. Antioxid Redox Signal 2019;30(10):1331–1351. pmid:29130312

26. Jai J, Arif A, Terenzi F, Willard B, Plow EF, Hazen SL, et al. Target-selective protein S-nitrosylation by sequence motif recognition. Cell 2014;159(3):623–634. pmid:25417112

27. Jaffrey SR, Snyder SH. The biotin switch method for the detection of S-nitrosylated proteins. Sci STKE 2001;(86):pl1 pmid:11752655

28. Qu Z, Meng F, Bomgarden RD, Viner RI, Li J, Rogers JC, et al. Proteomic quantification and site-mapping of S-nitrosylated proteins using isobaric iodoTMT reagents. J Proteome Res 2014;13:3200–3211. pmid:24926564

29. Cespedes MV, Espina C, Carcia-Cabezas MA, Trias M, Boluda A, Gomez del Pulgar MT, et al. Orthotopic microinjection of human colon cancer cells in nude mice induces tumor foci in all clinically relevant metastatic sites. Am J Pathol 2007;170(3):1077–1085 pmid:17322390

30. Sandu S, Wu X, Nabi Z, Rastegar M, Kung S, Mai S, et al. Loss of HLTF function promotes intestinal carcinogenesis. Mol Cancer. 2012;11:18 pmid:22452792

31. Rajput A, San Martin DS, Rose R, Beko A, LeVea C, Sharratt E, et al. Characterization of HCT116 human colon cancer cells in an orthotopic model. J Surg Res. 2008;147(2):276–281. pmid:17961596

32. Fernandez Y, Foradada L, Garcia-Aranda N, Mancilla S, Suarez-Lopez L, Cespedes MV, et al. Bioluminescent imaging of animal models for human colorectal cancer tumor growth and metastatic dissemination to clinically significant sites. J Mol Bio & Mol Imaging 2015;2(2):1019–1023.

33. Helmer RA, Kaur G, Smith LA, Chilton BS. Helicase-like transcription factor (Hltf) gene-deletion promotes oxidative phosphorylation (OXPHOS) in colorectal tumors of AOM/DSS-treated mice. PlosOne 2019;14(8):e0221751. pmid:31461471

34. Coccolini F, Gheza F, Lotti M, Virzi S, Lusco D, Ghermandi C, et al. Peritoneal carcinomatosis. World J Gastroenterol 2013;19(41):6979–6994. pmid:24222942

35. Helmer RA, Foreman O, Dertien JS, Panchoo M, Bhakta SM, Chilton BS. Role of helicase-like transcription factor (Hltf) in the G2/M transition and apoptosis in brain. PLoS One. 2013;8(6):e66799. pmid:23826137

36. Helmer RA, Martinez-Zaguilan R, Dertien JS, Fulford C, Foreman O, Peiris V, et al. Helicase-like transcription factor (Hltf) regulates G2/M transition, Wt1/Gata4/Hif-1α cardiac transcription networks, and collagen biosynthesis. PLoS One. 2013 Nov 20;8(11):e80461. pmid:24278285 eCollection 2013.

37. Kaur G, Helmer RA, Smith LA, Martinez-Zaguilan R, Dufour JM, Chilton BS. Alternative splicing of helicase-like transcription factor (Hltf): Intron retention dependent activation of immune tolerance at the feto-maternal interface. PLos One 2018;13(7):e0200211 pmid:29975766

38. Shen Z, Feng X, Fang Y, Li Y, Li Z, Zhan Y, et al. POTEE drives colorectal cancer development via regulating SPHJ1/p65 signaling. Cell Death Dis 2019;10:863–878. pmid:31723122

39. Pan S, Deng Y, Fu J, Zhang Y, Zhang Z, Ru X, et al. TRIM52 promotes colorectal cancer cell proliferation through the STAT3 signaling. Cancer Cell Int. 2019;19(57). pmid:30918473

40. Hansen L, Comyn S, Mang Y, Lind-Thomsen A, Myhre L, Jean F, et al. The myosin chaperone UNC45B is involved in lens development and autosomal juvenile cataract. Eur J Hum Genet. 2014:22(11):1290–1297 pmid:24549050

41. Kim SE, Park HY, Yoon H. Lee JE, Kim N, Sung MK. Sex- and gender-specific disparities in colorectal cancer risk. World J Gastroenterol. 2015;21(17):5167–5175. pmid:25954090

42. Hendifar A, Yang D, Lenz F, Lurje G, Pohl A, Lenz C, et al. Gender disparities in metastatic colorectal cancer survival. Clin Cancer Res. 2009;15(20):6391–6397. pmid:19789331

43. Amitay EL, Carr PR, Jansen L, Alwers E, Roth W, Helpel E, et al. Postmenopausal hormone replacement therapy and colorectal cancer risk by molecular subtypes and pathways. Int J Cancer. 2020;147(4):1–9. pmid:31943160

44. Lee MS, Menter DG, Kopetz S. Right versus left colon cancer biology: Integrating the consensus molecular subtypes. J Natl Compr Canc Netw. 2017;15(3)”411–419. pmid:28275039

45. Batista PJ, Cha ng HY. Long noncoding RNAs: cellular address codes in development and disease. Cell 2013;152(6):1298–1307. pmid:23498938

46. Gutschner T, Diederichs S. The hallmark of cancer: a long non-coding RNA point of view. RNA Biology. 2012;9(6):703–719. pmid:22664915

47. Carlevaro-Fita J, Lanzos A, Feuerbach L, Hong C, Mas-Ponte D, Pedersen JS, et al. Cancer LncRNA census reveals evidence of deep functional conservation of long noncoding RNAs is tumorigenesis. Commun Biol 2020;3(56): pmid:32024996

48. Li Y, Liu Y, Liu S, Wu F, Li S, Yang F, et al. Differential expression profile of long non-coding RNAs in human thoracic aortic aneurysm. J Cell Biochem. 2018;119(10):7991–7997. pmid:29323743

49. Sokelaand G, Schumacher U. The functional role of integrins during intra- and extravasation within the metastatic cascade. Mol Cancer 2019;18:12. pmid:30657059

50. Jung M, Yanping Y, McCloskey JE, Zaman M, Vedvyas Y, Zhang X, et al., Chimeric antigen receptor T cell therapy targeting ICAM-1 in gastric cancer. Mol Ther Oncolytics 2020;18(25):587–601. pmid:32995483

51. Marino-Crespo O, Cervas-Alvarez E, Harding AL, Murdoch C, Fernandez-Briera A, Gil-Martin E. Haptoglobin expression in human colorectal cancer. Histol Histopathol 2019;34(8):953–963. pmid:30864745

52. Cao LL, Pei XF, Qiao X, Yu J, Ye H, Xi CL, et al. SERPINA3 silencing inhibits the migration, invasion, and liver metastasis of colon cancer cells. Dig Dis Sci. 2018;63(9):2309–2319. pmid:29855767

53. Yao JF, Li XJ, Yan LK, He S, Zheng JB, Wang XR, et al. Role of HGF/c-Met in the treatment of colorectal cancer with liver metastasis. J Biochem Mol Toxicol. 2019;33(6):e22316. pmid:30897285

54. Kany S, Vollrath JT, Relja B. Cytokines in inflammatory disease. Int J Mol Sci 2019;20(23):6008–6038. pmid:31795299

55. Tiwara A, Saraf S, Jain A, Panda PK, Verma A, Jain SK. Basics to advances in nanotherapy of colorectal cancer. Drug Deliiv and Transl Res. 2020;10:319–338.

56. Vinchhi P, Patel MM. Triumph against cancer: invading colorectal cancer with nanotechology. Expert Opin Drug Deliv. 2021; pmid:33567909

57. Schaer CA, Deuer JW, Schildknecht D, Mahmoudi L, Garcia-Rubio I, Owczarek C, et al. Haptoglobin preserves vascular nitric oxide signaling during hemolysis. Am J Respir Crit Car Med. 2016;193(10):1111–1122. pmid:26694989

58. Li C, Wang . Mechanical tumor microenvironment and transduction: cytoskeleton mediates cancer cell invasion and metastasis. Int J Biol Sci. 2020;16(12):2014–2018. pmid:32549750

59. Azadu S, Shadpour MT. Chapter Five–The microenvironmental and cytoskeletal remodeling in tumour cell invasion. Int Rev Cell Mol Biol. 2020;365:257–289.

60. Masi I, Caprara V, Bagnato A, Rosano L. Tumor cellular and microenvironmental clues controlling invadopodia formation. Front Cell Dev Biol. 2020 Oct 15;8:584181. pmid:33178698 eCollection 2020.

61. Bern RN, Ford JM. 2019 Management of peritoneal malignancies. Textbook of Gastrointestinal Oncology In: Yalcin S., Philip P. (eds) Textbook of Gastrointestinal Oncology. Springer, Cham. https://doi.org/10.1007/978-3-030-18890-0_23

62. Pretzsch E, Bosch F, Neumann J, Ganschow P, Bazhin A, Guba M, et al. Mechanisms of metastasis in colorectal cancer and metastatic organotrophism: hematogenous versus peritoneal spread. J Oncol 2019;2019:7407190. pmid:31641356

63. Smit JA, Baas JM, Verheijen PM. Metastasis of colorectal carcinoma in the testis: first sign of peritoneal disease. J Surg Case Rep 2019;2019(10):rjz270 pmid:31687127

64. Bays JL, DeMali KA. Vinculin in cell-cell and cell-matrix adhesions. Cell Mol Life Sci 2017;74(16):2999–309 pmid:28401269

65. Goldmann WH, Auernheimer V, Thievessen I, Fabry B. Vinculin, cell mechanics and tumour cell invasion. Cell Biol Int. 2013;37:397–405. pmid:23494917

66. Rubashkin MG, Cassereau L, Bainer R, DuFort CC, Yui Y, Ou G, et al. Force engages vinculin and promoter tumor progression by enhancing PI3K activation of phosphatidylinositol (3,4,5)-triphosphate. Cancer Res 2014;74(17):4697–4611

67. Mnatsakanyan R, Markoutsa S, Walbrunn K, Roos A, Verhelst SHL, Zahedi RP. Proteome-wide detection of S-nitrosylation targets and motifs using bioorthogonal cleavable-linker-based enrichment and switch technique. Nat Commun 2019;10:2195–2206. pmid:31097712

68. Sha Y, Marshall HE. S-nitrosylation in the regulation of gene transcription. Biochim Biophys Acta 2012;1820(6):701–711. pmid:21640163

69. Kluska K, Adamczyk J, Krezel A. Metal binding properties, stability and reactivity of zinc fingers. Coord Chem Rev. 2018;367(15):18–64.

70. Mi H, Muruganujan A, Ebert D, Huang X, Thomas PD. PANTHER version 14: more genomes, a new PANTHER GO-slim and improvements in enrichment analysis tools. Nuc Acids Res 2019;47:D419–D426. pmid:30407594

71. Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, et al. STRING v11:protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nuc Acids Res 2019;47:D607–D613. pmid:30476243