Increased High Mobility Group Protein A2/SMAD3 Relates to Ovarian Cancer Progression

Introduction: The high mortality associated with ovarian cancer is generally related to the development of drug-resistant disease. HMGA2 protein, a member of the high-mobility group AT-hook (HMGA) family of non-histone chromatin binding factors, is overexpressed in high-grade serous ovarian and tubal carcinomas, though little is known about its contribution to disease progression and drug resistance. We sought to assess whether compositional changes in HMGA2 production were associated with ovarian cancer progression.

The distributions of intrinsic disorder propensity within the amino acid sequences of human HMGA2 (109 amino acids, UniProt ID: P52926) and SMAD3 (425 residues, UniProt ID: P84022) were determined using PONDR®VLXT [25], PONDR®VL3 [26], PONDR®VSL2B [27] and PONDR-FIT [28] algorithms. Potential amyloidogenic regions, regions responsible for amyloid formation and protein aggregation, were identified in HMGA2 protein using the FoldAmyloid software program [29]. Intrinsic regions of protein disorder within the human HMGA2 protein were predicted using the ANCHOR algorithm and the IUPred disorder prediction tool [30,31]. Potential protein binding partners of HMGA2 and SMAD3 were predicted using the String 9.1 platform [32]. Potential post-translational glycosylation modification sites within the human SMAD3 sequence were identified by the HMMpTM algorithm [33].
With prior University of South Florida Institutional Review Board committee approval, anonymized urine samples from women with benign ovarian disorders (N=3), patients with epithelial OC (N=13) and healthy controls (N=5) were released for this pilot research project. Where possible, anonymized matched H & E sections from paraffin blocks were reviewed to confirm the histologic diagnosis according to International Federation of Gynecology and Obstetrics (FIGO) scores. Urine samples were kept on ice following collection, centrifuged at 3000 × g, aliquoted and stored at -20 °C. For initial comparisons of urinary levels of HMGA2 by western immunblotting (WB) between healthy controls (N=5) and two sets of OC (N=5 each), urine samples were pooled to ensure sufficient amounts of protein for analyses. For comparisons of urinary levels of HMGA2 and SMAD3 from individual patients with benign gynecologic disease (N=3) and OC (N=3), equivalent volumes of urine were concentrated using Amicon Centriplus centrifugal filters (Millipore, Billerica, MA) to ensure sufficient amounts of protein for WB analyses.
Cell growth was measured by the CellTiter 96® AQ ueous One Solution Cell Proliferation MTS (Promega, Madison, WI) colorimetric assay. The assay was performed in 96 well microtiter plates according to manufacturers' instructions. Two thousand cells were plated on 96 well microtiter plates and selected for carboplatin (CB) resistance at intervals following bi-weekly treatments of 4 x 5 µM CB, followed by 4 x 10 µM CB. Absorbance at 490 nm was measured at 0, 24, 48, and 72 hours following treatment using an EL x 800 microplate reader (Bio-Tek Instruments, Winooski, VT) and the results expressed as the mean absorbance ± standard error (S.E.).
High Mobility Group protein A2 (HMGA2) is a member of the non-histone chromosomal high mobility group protein family. HMGA2 contains 3 DNA-binding domains, or A/T hooks, and acts as a transcriptional regulator [4]. HMGA2 is overexpressed in a number of cancers [5][6][7][8][9], including OC [10] and especially in highly aggressive Type II high-grade papillary serous carcinomas of the fallopian tube and ovary [11][12][13]. These OCs contain genetic instability distinguished by a p53 mutation signature [14]. Overexpression of HMGA2 is associated with increased tumor growth [15], metastasis [16], poor clinical outcome [17,18] and emergence of drug-resistant disease [19]. HMGA2 appears to mediate tumor progression by promoting epithelial-mesenchymal transition (EMT) [20]. Additionally, HMGA2 may confer drug resistance through increased DNA base excision repair due to inherent DNA lyase activity [21]. Since silencing HMGA2 results in reduced tumor growth [22][23][24], HMGA2 may represent an important therapeutic target. In order to better understand the mechanisms that may regulate OC development, we sought to determine whether compositional changes in HMGA2 are related to OC progression and emergence of drug resistance.

Computational analyses
Tissue culture

Western Immunoblotting
For immunoprecipitation (IP), cells were washed with PBS and incubated on ice with radioimmunoprecipitation assay (RIPA) lysis buffer containing protease inhibitors. After 15 minutes, the cells were removed by scraping, the lysates were incubated for 60 minutes on ice and then the lysates were centrifuged at 10,000 x g for 10 minutes at 4 °C. Two µg of primary antibody were incubated with 1 mg of whole cell lysate for 2 hours at 4 °C. Protein A/G agarose suspension was added to lysate/antibody mixture followed by incubation at 4 °C on a rocker overnight. The suspension was collected by centrifugation at 500 x g for 2 minutes at 4 °C. The pellet was washed with RIPA buffer containing protease inhibitors followed by WB analysis.
For MTS assays, student's t-tests were performed to assess statistical difference between means of triplicates ± S.E. from three separate experiments.
Computational analyses of HMGA2 using four disorder predictors of the PONDR family indicated that HMGA2 is a highly disordered protein throughout its whole 109 amino acid length ( Figure 1A). Therefore, the HMGA2 protein does not have a stable secondary and/or tertiary structure, existing instead as a highly dynamic ensemble of multiple interconverting forms. Since intrinsically disordered proteins are associated with low aggregability, the FoldAmyloid program was used to assess the aggregation propensity of HMGA2. As indicated in Figure 1B, failure to achieve threshold status confirmed that HMGA2 has a low aggregation propensity consistent with properties of a disordered protein.

Results and Discussion
Immunoprecipitation Statistics HMGA2 is a highly disordered protein capable of binding many proteins A B

Figure 1: HMGA2 is an intrinsically disordered protein. A)
The disorder propensity of human HMGA2 (UNiProt ID:P52926) was evaluated by the members of the PONDR family of disorder predictors. Scores above 0.5 (arrow) correspond to disordered residues/regions. B) Failure to surpass threshold status (dotted line) with FoldAmyloid program analysis indicated that HMGA2 is essentially non-amyloidogenic and has a low propensity to self-aggregate. C) The ANCHOR prediction tool was employed to identify potential protein-binding sites in HMGA2. Three regions encompassing amino acid residues 1-33, 52-73 and 78-93 were identified as potential protein binding sites within HMGA2. D) Analysis of the interactivity of human HMGA2 by String produces the network of predicted associations for a particular group of proteins. The network nodes are proteins, whereas the edges represent the predicted or known functional associations. An edge may be drawn with up to 7 differently colored lines that represent the existence of the seven types of evidence used in predicting the associations. A red line indicates the presence of fusion evidence; a green line -neighborhood evidence; a blue line -co-occurrence evidence; a purple line -experimental evidence; a yellow line -text mining evidence; a light blue line -database evidence; a black line -co-expression evidence [32].
The intrinsic plasticity of such a highly disordered protein also implies that HMGA2, as typical of intrinsically disordered proteins, can interact with a multitude of protein binding partners [35]. This assumption was supported by ANCHOR-based analysis which identified three discrete regions (residues 1-33, 52-73 and 78-93) as potential protein binding sites within HMGA2 ( Figure 1C). Interestingly, these sites correspond to the three HMGA2 AT hooks located at amino acids 25-33, 45-51 and 73-83 [36], respectively. Lastly, results of the String-based analysis show that HMGA2 can interact with numerous protein partners, including SMAD 3 and SMAD 4 proteins ( Figure 1D). SMAD3 belongs to a protein family involved in receptor-mediated transforming growth factor-β (TGF-β) signal transduction [37]. Activation by TGF-β results in SMAD2 and SMAD3 phosphorylation with subsequent trimerization between SMAD2, SMAD3 and SMAD4. Additionally, HMGA2 has been shown to physically interact with SMAD3, and to a lesser extent with SMAD2 and SMAD4, in human embryonic kidney cells [38]. The SMAD complex then translocates into the nucleus and binds to DNA to drive transcription for genes involved in EMT [39]. SMAD3 is also a moderately disordered protein with predicted disordered regions C D IUPred ANCHOR encompassing approximately 20% of its sequence (Figure 2A). Interestingly, the disordered regions of SMAD3 encompass portions of its amino MH1 domain (amino acid residues 31-131), which, with its carboxy MH2 (amino acid residues 226-403) domains [40] are responsible for binding HMGA2 [38]. Analyses of potential post-translational modifications in SMAD3 by HMMpTM software program predicted glycosylation of SMAD3 at amino acid residues 339, 344, 345, 357, 388 and 418 associated with a SMAD3 transmembrane domain at residues 296 -314 ( Figure 2B). Therefore, HMGA2 could be recruited to cellular membranes via its interaction with SMAD3. Finally, Figure 2C shows that similar to the HMGA2 protein, SMAD3 is predicted to have a large number of potential binding partners, including HMGA2, and, therefore, also serves as an important hub protein. Consequently, while TGF-β-mediated signaling activates the formation and nuclear translocation of an HMGA2/SMAD complex to drive EMT, computational analyses of an HMGA2/SMAD3 complex suggests that this complex could be recruited to the cell membrane for possible extracellular transport.
Given reports of HMGA2 overexpression in OC, an HMGA2/SMAD3 complex could be associated with OC progression. To expand our understanding of the role of HMGA2 in OC, we surveyed HMGA2 protein levels by WB among HIOSE cells and cancer cell lines ( Figure 3). Normal cells (HIOSE 1816-575, 117, 118 and MCC3) failed to show any detectable levels of HMGA2 at the expected 18 kDa HMGA2 size consistent with published reports that HMGA2 is not expressed in normal somatic cells [41]. The colon cancer cell line, SW626, also demonstrated undetectable levels of HMGA2 which may reflect degree of tumor progression. A recent report by Rizzi et al. [42] showed that HMGA2 negative colon cancers surrounded by HMGA2 positive stromal fibroblasts were associated with improved clinical outcome compared to HMGA2 positive colon cancers in an environment of HMGA2 negative stromal fibroblasts. A possible role of genetic and epigenetic events in cancer associated stroma has been similarly described for prostate cancer development [43] and underscores the role of stroma for tumor progression. High expression of HMGA2 was noted in 5 out of 7 OC cell lines, consistent with the literature reporting elevated levels of HMGA2 in OC [10][11][12][13]. Interestingly, among pairs of drug-sensitive parental OC cell lines (A2870s, OV2008) and their respective drug-resistant daughter cell lines (A2780cp, C13), we found that levels of HMGA2 were increased up to 850x in drug-resistant daughter cell lines ( Figure  3A), thereby further supporting overexpression of HMGA2 with OC progression.

Levels of HMGA2 increase with ovarian transformation and cancer progression in culture
The reported molecular weights for HMGA2 and SMAD3 are 18 and 52 kDa, respectively. While disordered proteins typically display aberrant electrophoretic migration as larger molecular weight proteins due to their high ratio of charged to hydrophobic residues which subsequently interferes with SDS binding [44], we noted additional higher molecular weight HMGA2 protein bands at 100 and 200 kDa on WB, generally limited to OC cell lines ( Figure 3B). Additional WB analyses from OV2008 and C13 cell lysates not only confirmed increased levels of HMGA2 protein in drug-resistant cells ( Figure 4A, B), but also confirmed elevated levels of SMAD3 present in western immunblots at 100 and 200 kDa ( Figure 4A ,B) suggesting that HMGA2 and SMAD3 might be present in a complex in OC cell cultures. Levels of high molecular weight HMGA2 and SMAD3 were 98% and 30% greater, respectively, in the drug-resistant C13 cells compared to their parental drug-sensitive OV2008 OC cells. To differentiate between non-specific binding and the presence of larger molecular weight protein species associated with HGMGA2/SMAD3 protein complexes, immunoprecipitation (IP) was performed followed by WB on OV2008 and C13 cell lysates. We found that IP for HMGA2 co-localized with SMAD3 and IP for SMAD3 co-localized with HMGA2 at molecular weights of 100 and 200 kDa ( Figure 4C) confirming for the first time the presence of HMGA2/SMAD3 as components of a complex in cultured OC cells. Consequently, our data are in agreement with others [38] and underscore the possible contribution of an HMGA2/SMAD3 complex for tumor progression by enhanced EMT. About 20% of OC tumors exhibit primary resistance to platinum compounds and the majority of women who respond to initial chemotherapy relapse with emergence of drug resistant disease [2,3]. To determine whether the HMGA2/SMAD3 complex has a clinical impact related to drug-resistance, we measured levels of HMGA2, SMAD3 and their downstream target, SNAIL1 [38,39], in a temporal manner during selection of drug resistance in cell culture. The drug sensitive OC cell line, OV90, was sequentially Levels of HMGA2/SMAD3 and SNAIL1 increase with emergence of drug-resistance A B C treated over 4 months with biweekly treatments of 4x5 µM CB followed by biweekly treatments of 4x10 µM CB. Selection for drug-resistance was accompanied with a 70% increase in cell growth ( Figure 5A) compared with the parental cell line. When challenged with a single dose of 25 µM CB, the parental cells remained responsive to CB-induced cytotoxicity and demonstrated a 35% reduction in cell viability. In contrast, when challenged with a single dose of 25 µM CB, cells selected for drug resistance demonstrated only 23% reduction in cell viability compared to their respective untreated controls. Further, cells selected for drug resistance and treated with a single dose of 25 µM CB demonstrated overall growth more than twice that of parental cells similarly treated with a single dose of 25 µM CB. Enhanced genotoxic survival and cell proliferation appeared related to HMGA2/SMAD3 complex formation ( Figure 5B, C) in which increased levels of 200 kDa SMAD3 preceded overexpression of 200 kDa HMGA2. This is in keeping with previous reports that SMAD3 and SMDA4 are crucial for induction of HMGA2 transcription [38]. In addition, increased HMGA2 and SMAD3 protein expression was related to enhanced SNAIL1 expression ( Figure 5D). Consequently, increased production of an HMGA2/SMAD3 complex may be related to drug resistance and cancer progression since enhanced expression of SNAIL1 has been reported to promote emergence of drug-resistant disease and increased cancer cell proliferation by maintaining an EMT phenotype [45]. Journal of Gynecology Research   [46], suggesting the possibility of extracellular HMGA2 localization. Given that our computational data suggested that HMGA2 could be recruited to the cell membrane by SMAD3, we sought to determine whether HMGA2 and HMGA2 complexes could be found extracellularly in OC. Conditioned media (CM) from OV2008 and C13 cells were analyzed by WB for expression of large molecular weight (200 kDa) HMGA2/SMA3 bands. We not only found 200 kDa HMGA2/SMAD3 in the CM suggesting that HMGA2/SMAD3 complex can be secreted by OC cells, but noted that there were 88% and 10% greater levels of 200 kDa HMGA2 and SMAD3, respectively, in the CM of C13 cells compared to OV2008 cells ( Figure 6A).

OC cells secrete HMGA2/SMAD3
For clinical validation, we measured HMGA2 and SMAD3 protein levels by WB in urinary samples obtained from healthy controls and women with gynecological disorders. In order to obtain sufficient amounts of protein for WD, initial studies comparing urinary levels of 18 kDa HMGA2 among healthy controls and women with OC necessitated pooling urine samples. In comparison to a positive control consisting of a C13 OC cell lysate, negligible levels of HMGA2 were found in the urine of healthy donors ( Figure 6B). However, among 2 pooled urinary samples (representing a total of 10 patients), HMGA2 protein levels were 10 fold greater than healthy controls. Likewise, greater levels of high MW (200 kDa) HMGA2/SMAD3 proteins were found in the urine of OC patients (>1000 times) compared to urine from women with benign disease ( Figure 6C).
Given the relationship between levels of the HMGA2/SMAD3 complex with OC progression, further studies of urinary levels of HMGA2/SMAD3 alone or in combination with other urinary biomarkers [47][48][49] are warranted. These may offer a non-invasive, convenient and cost-effective prognostic tool for this lethal gynecologic malignancy. This may be especially relevant when assessing co-morbidity factors associated with OC. For example, a recent study by Califano et al. demonstrated a link between high body mass index (BMI) and overexpression of HMGA2 in OC tissues with poor clinical outcome and poor response to chemotherapy [50]. Importantly, since the HMGA2/SMAD3 complex is secreted by OC cells and is present in urine, changes in the levels or components of the circulating HMGA2/SMAD3 complex may have clinical impact. Lastly, by beginning to understand the components and mechanisms of an HMGA2/SMAD3 complex that drives SNAIL1-mediated drug resistance in OC, it may be possible to abrogate complex formation for therapeutic intervention. Taken together, the identification of an HMGA2/SMAD3 complex in OC may provide a new target for clinical management as urinary levels of such a complex may be of diagnostic and prognostic relevance in this gynecologic malignancy. Interference with the HMGA2/SMAD3 complex itself, then, may be therapeutically advantageous.

Conclusion
This work was supported, in part, by a Program Project Developmental Grant (PDD/USF/01.2011) from the Ovarian Cancer Research Fund to SVN and PAK.