Cellular heterogeneity profiling by hyaluronan probes reveals an invasive but slow-growing breast tumor subset

Tumor heterogeneity confounds cancer diagnosis and the outcome of therapy, necessitating analysis of tumor cell subsets within the tumor mass. Elevated expression of hyaluronan (HA) and HA receptors, receptor for HA-mediated motility (RHAMM)/HA-mediated motility receptor and cluster designation 44 (CD44), in breast tumors correlates with poor outcome. We hypothesized that a probe for detecting HA–HA receptor interactions may reveal breast cancer (BCa) cell heterogeneity relevant to tumor progression. A fluorescent HA (F-HA) probe containing a mixture of polymer sizes typical of tumor microenvironments (10–480 kDa), multiplexed profiling, and flow cytometry were used to monitor HA binding to BCa cell lines of different molecular subtypes. Formulae were developed to quantify binding heterogeneity and to measure invasion in vivo. Two subsets exhibiting differential binding (HA^−/low vs. HA^high) were isolated and characterized for morphology, growth, and invasion in culture and as xenografts in vivo. F-HA–binding amounts and degree of heterogeneity varied with BCa subtype, were highest in the malignant basal-like cell lines, and decreased upon reversion to a nonmalignant phenotype. Binding amounts correlated with CD44 and RHAMM displayed but binding heterogeneity appeared to arise from a differential ability of HA receptor-positive subpopulations to interact with F-HA. HA^high subpopulations exhibited significantly higher local invasion and lung micrometastases but, unexpectedly, lower proliferation than either unsorted parental cells or the HA^−/low subpopulation. Querying F-HA binding to aggressive tumor cells reveals a previously undetected form of heterogeneity that predicts invasive/metastatic behavior and that may aid both early identification of cancer patients susceptible to metastasis, and detection/therapy of invasive BCa subpopulations.Breast tumors display substantial heterogeneity driven by genetic and epigenetic mechanisms (1–3). These processes select and support tumor cell subpopulations with distinct phenotypes in proliferation, metastatic/invasive proclivity, and treatment susceptibility that contribute to clinical outcomes. Currently, there is a paucity of biomarkers to identify these subpopulations (3–12). Although detection of genetic heterogeneity may itself be a breast cancer (BCa) prognostic marker (3, 13–15), the phenotypes manifested from this diversity are context-dependent. Therefore, phenotypic markers provide additional powerful tools for biological information required to design diagnostics and therapeutics. Glycomic approaches have enormous potential for revealing tumor cell phenotypic heterogeneity because glycans are themselves highly heterogeneous and their complexity reflects the nutritional, microenvironmental, and genetic dynamics of the tumors (16–18).We used hyaluronan (HA) as a model carbohydrate ligand for probing heterogeneity in glycosaminoglycan–BCa cell receptor interactions. We reasoned this approach would reveal previously undetected cellular and functional heterogeneity linked to malignant progression because the diversity of cell glycosylation patterns, which can occur as covalent and noncovalent modifications of proteins and lipids as well as different sizes of such polysaccharides as HA, is unrivaled (16, 17, 19). In particular, tumor and wound microenvironments contain different sizes of HA polymers that bind differentially to cell receptors to activate signaling pathways regulating cell migration, invasion, survival, and proliferation (19–22).More than other related glycosaminoglycans, HA accumulation within BCa tumor cells and peritumor stroma is a predictor of poor outcome (23) and of the conversion of the preinvasive form of BCa, ductal carcinoma in situ, to an early invasive form of BCa (24). HA is a nonantigenic and large, relatively simple, unbranched polymer, but the manner in which it is metabolized is highly complex (19, 25). There are literally thousands of different HA sizes in remodeling microenvironments, including tumors. HA polymers bind to cells via at least six known receptors (16, 19, 20, 26–32). Two of these, cluster designation 44 (CD44) and receptor for HA-mediated motility/HA-mediated motility receptor (RHAMM/HMMR), form multivalent complexes with different ranges of HA sizes (19, 29, 33), and both receptors are implicated in BCa progression (19–21, 23, 29, 30, 33–36). Elevated CD44 expression in the peritumor stroma is associated with increased relapse (37), and in primary BCa cell subsets may contribute to tumor initiation and progression (38–40). Elevated RHAMM expression in BCa tumor subsets is a prognostic indicator of poor outcome and increased metastasis (22, 33, 41). RHAMM polymorphisms may also be a factor in BCa susceptibility (42, 43).We postulated that multivalent interactions resulting from mixture of a polydisperse population of fluorescent HA (F-HA) sizes, typical of those found in remodeling microenvironments of wounds and tumors (19, 20, 29), with cellular HA receptors would uncover a heterogeneous binding pattern useful for sorting tumor cells into distinct subsets. We interrogated the binding of F-HA to BCa lines of different molecular subtypes, and related binding/uptake patterns to CD44 and RHAMM display, and to tumor cell growth, invasion, and metastasis.