Targeting Wnt signaling for improved glioma immunotherapy

Introduction

Despite aggressive standard-of-care therapy, including surgery, radiation, and chemotherapy, glioblastoma recurrence is almost inevitable and uniformly lethal. Activation of glioma-intrinsic Wnt/β-catenin signaling is associated with a poor prognosis and the proliferation of glioma stem-like cells, leading to malignant transformation and tumor progression. Impressive results in a subset of cancers have been obtained using immunotherapies including anti-CTLA4, anti-PD-1, and anti-PD-L1 or chimeric antigen receptor (CAR) T cell therapies. However, the heterogeneity of tumors, low mutational burden, single antigen targeting, and associated antigen escape contribute to non-responsiveness and potential tumor recurrence despite these therapeutic efforts. In the current study, we determined the effects of the small molecule, highly specific Wnt/CBP (CREB Binding Protein)/β-catenin antagonist ICG-001, on glioma tumor cells and the tumor microenvironment (TME)-including its effect on immune cell infiltration, blood vessel decompression, and metabolic changes.

Methods

Using multiple glioma patient-derived xenografts cell lines and murine tumors (GL261, K-Luc), we demonstrated in vitro cytostatic effects and a switch from proliferation to differentiation after treatment with ICG-001.

Results

In these glioma cell lines, we further demonstrated that ICG-001 downregulated the CBP/β-catenin target gene Survivin/BIRC5-a hallmark of Wnt/CBP/β-catenin inhibition. We found that in a syngeneic mouse model of glioma (K-luc), ICG-001 treatment enhanced tumor infiltration by CD3+ and CD8+ cells with increased expression of the vascular endothelial marker CD31 (PECAM-1). We also observed differential gene expression and induced immune cell infiltration in tumors pretreated with ICG-001 and then treated with CAR T cells as compared with single treatment groups or when ICG-001 treatment was administered after CAR T cell therapy.

Discussion

We conclude that specific Wnt/CBP/β-catenin antagonism results in pleotropic changes in the glioma TME, including glioma stem cell differentiation, modulation of the stroma, and immune cell activation and recruitment, thereby suggesting a possible role for enhancing immunotherapy in glioma patients.

A Deeply Quiescent Subset of CML LSC depend on FAO yet Avoid Deleterious ROS by Suppressing Mitochondrial Complex I

BACKGROUND AND OBJECTIVE

Disease relapse and therapy resistance remain serious impediments to treating cancer. Leukemia stem cells (LSC) are therapy resistant and the cause of relapse. A state of deep quiescence appears to enable cancer stem cells (CSC) to acquire new somatic mutations essential for disease progression and therapy resistance. Both normal hematopoietic stem cells (HSC) and LSC share many common features, thereby complicating the safe elimination of LSC. A recent study demonstrated that long lived normal oocytes exist without mitochondrial complex I (MC-1), expressing it in a developmentally regulated fashion, thereby mitigating their vulnerability to ROS. Quiescent CSC rely on mitochondrial FAO, without complex I expression, thereby avoiding the generation of damaging ROS, similar to long lived normal human stem cells. A deeper understanding of the biology of therapy resistance is important for the development of optimal strategies to attain complete leukemia cures.

METHODS

Here, using scRNA-sequencing and ATAC-seq on primary chronic myelogenous leukemia (CML) patient samples, combined with bioinformatics analyses, we further examine the heterogeneity of a previously characterized in vitro imatinib-selected CD34-CD38- CML LSC population. We utilized a series of functional analyses, including single-cell metabolomic and Seahorse analyses, to validate the existence of the deepest quiescent leukemia initiators (LI) subset.

RESULTS

Current study revealed heterogeneity of therapy resistant LSC in CML patients and their existence of two functionally distinct states. The most deeply quiescent LI suppress the expression of MC-1, yet are highly dependent on fatty acid oxidation (FAO) for their metabolic requirements and ATAC-seq demonstrated increased chromatin accessibility in this population, all consistent with an extremely primitive, quiescent stemness transcriptional signature. Importantly, the specific CREB binding protein (CBP)/β-catenin antagonist ICG-001 initiates the differentiation of LSC, including LI, decreases chromatin accessibility with differentiation and increasing expression of MC-1, CD34, CD38 and BCR-ABL1, thereby resensitizing them to imatinib.

CONCLUSION

We investigated the biological aspects related to LSC heterogeneity in CML patients and demonstrated the ability of specific small molecule CBP/β-catenin antagonists to safely eliminate deeply quiescent therapy resistant CSC. These observations may represent an attractive generalizable therapeutic strategy that could help develop better protocols to eradicate the quiescent LSC population.

E7386 is not a Specific CBP/β-Catenin Antagonist

BACKGROUND AND OBJECTIVE

The first clinically evaluated CBP/β-catenin antagonist, PRI-724, displayed an excellent safety profile administered intravenously via continuous infusion. Eisai recently disclosed a third-generation, orally available, reportedly CBP/β-catenin antagonist, E7386. However, several structural features and the reported cytotoxicity of E7386 were unexpected for a specific CBP/β-catenin antagonist. Therefore, we undertook a comparison of E7386 versus the highly specific bona fide CBP/β-catenin antagonists, ICG-001 and C82, the active agents derived from the prodrug PRI-724.

INTRODUCTION

CBP/β-catenin antagonists rebalance the equilibrium between CBP/β-catenin and p300/β-catenin dependent transcription and may be able to treat or prevent many diseases of aging via maintenance of somatic stem cell pool and regulating mitochondrial function and metabolism involved in differentiation and immune cell function. The safety, efficacy, and therapeutic potential of the specific CBP/β-catenin antagonists, ICG-001, and the second-generation compound, C82, the active agent derived from the pro-drug PRI-724, have been studied extensively in a variety of preclinical disease models and in the clinic for oncology and hepatic fibrosis. However, the lack of oral bioavailability has hampered the further development of PRI-724. Thus, Eisai recently proposed a third-generation, orally available, reportedly CBP/β-catenin antagonist E7386. Here, we have performed a comparative analysis of E7386 with the highly specific bona fide CBP/β-catenin antagonists, ICG-001 and C82.

METHODS

We utilized a series of previously validated biochemical and transcriptional assays to investigate the selective targeting of the CBP/β-catenin interaction in conjunction with global transcriptional profiling to compare the three small molecules, ICG-001, C82, and E7386.

RESULT

Our data cast significant doubt that the mechanism of action of E7386 is via specific CBP/β-catenin antagonism.

CONCLUSION

It can thus be concluded that E7386 is not a specific CBP/β-catenin antagonist.

Differential Kat3 Coactivator Usage Regulates Brain Metabolism and Neuronal Differentiation

INTRODUCTION

Our previous work has demonstrated significant effects on the oxidative stress response, mitochondrial function, and oxidative phosphorylation in the livers and intestines of p300 S89A knockin (S89AKI) mice. We now show that this mutation is also associated with brain metabolic defects and neuronal differentiation.

METHODS

p300 S89A edited P19 cells, and S89AKI mice demonstrated metabolic and neuronal differentiation defects based on proteomic, cell biological and PET imaging studies.

RESULTS

The metabolic and differentiation defects associated with the p300 S89A knockin mutation could be corrected both in vitro and in vivo utilizing the small molecule CBP/beta-catenin antagonist ICG-001.

CONCLUSION

Rebalancing the equilibrium between CBP/β-catenin versus p300/β-catenin associated transcription, utilizing the small molecule CBP/beta-catenin antagonist ICG-001, enhances mitochondrial oxidative phosphorylation, metabolic function, and neuronal differentiation and may be able to ameliorate the cognitive decline seen in neurodegenerative disorders, including Alzheimer's Disease.

Differential Kat3 Usage Orchestrates the Integration of Cellular Metabolism with Differentiation

Simple Summary

The coupling of metabolism with cellular status is critically important and highly evolutionarily conserved. However, how cells coordinate metabolism with transcription as they change their status is not clear. Utilizing multiomic and functional studies, we now demonstrate the dichotomous roles of the Kat3 coactivators CBP and p300 and, in particular, their extreme N-termini, in coordinating cellular metabolism with cell differentiation. Using multiple in vitro and in vivo systems, our study sheds new light on metabolic regulation in homeostasis and disease, including cancer.

Abstract

The integration of cellular status with metabolism is critically important and the coupling of energy production and cellular function is highly evolutionarily conserved. This has been demonstrated in stem cell biology, organismal, cellular and tissue differentiation and in immune cell biology. However, a molecular mechanism delineating how cells coordinate and couple metabolism with transcription as they navigate quiescence, growth, proliferation, differentiation and migration remains in its infancy. The extreme N-termini of the Kat3 coactivator family members, CBP and p300, by far the least homologous regions with only 66% identity, interact with members of the nuclear receptor family, interferon activated Stat1 and transcriptionally competent β-catenin, a critical component of the Wnt signaling pathway. We now wish to report based on multiomic and functional investigations, utilizing p300 knockdown, N-terminal p300 edited and p300 S89A edited cell lines and p300 S89A knockin mice, that the N-termini of the Kat3 coactivators provide a highly evolutionarily conserved hub to integrate multiple signaling cascades to coordinate cellular metabolism with the regulation of cellular status and function.

p300 Serine 89: A Critical Signaling Integrator and Its Effects on Intestinal Homeostasis and Repair

Simple Summary

Given their high degree of identity and even greater similarity at the amino acid level, Kat3 coactivators, CBP (Kat3A) and p300 (Kat3B), have long been considered redundant. We describe the generation of novel p300 S89A knock-in mice carrying a single site directed amino acid mutation in p300, changing the highly evolutionarily conserved serine 89 to alanine, thus enhancing Wnt/CBP/catenin signaling (at the expense of Wnt/p300/catenin signaling). p300 S89A knock-in mice exhibit multiple organ system, immunologic and metabolic differences, compared with their wild type counterparts. In particular, these p300 S89A knock-in mice are highly sensitive to intestinal injury resulting in colitis which is known to significantly predispose to colorectal cancer. Our results highlight the critical role of this region in p300 as a signaling nexus and provide further evidence that p300 and CBP are non-redundant, playing definite and distinctive roles in development and disease.

Abstract

Differential usage of Kat3 coactivators, CBP and p300, by β-catenin is a fundamental regulatory mechanism in stem cell maintenance and initiation of differentiation and repair. Based upon our earlier pharmacologic studies, p300 serine 89 (S89) is critical for controlling differential coactivator usage by β-catenin via post-translational phosphorylation in stem/progenitor populations, and appears to be a target for a number of kinase cascades. To further investigate mechanisms of signal integration effected by this domain, we generated p300 S89A knock-in mice. We show that S89A mice are extremely sensitive to intestinal insult resulting in colitis, which is known to significantly increase the risk of developing colorectal cancer. We demonstrate cell intrinsic differences, and microbiome compositional differences and differential immune responses, in intestine of S89A versus wild type mice. Genomic and proteomic analyses reveal pathway differences, including lipid metabolism, oxidative stress response, mitochondrial function and oxidative phosphorylation. The diverse effects on fundamental processes including epithelial differentiation, metabolism, immune response and microbiome colonization, all brought about by a single amino acid modification S89A, highlights the critical role of this region in p300 as a signaling nexus and the rationale for conservation of this residue and surrounding region for hundreds of million years of vertebrate evolution.

CLDN18.1 attenuates malignancy and related signaling pathways of lung adenocarcinoma in vivo and in vitro

Claudins are a family of transmembrane proteins integral to the structure and function of tight junctions (TJ). Disruption of TJ and alterations in claudin expression are important features of invasive and metastatic cancer cells. Expression of CLDN18.1, the lung-specific isoform of CLDN18, is markedly decreased in lung adenocarcinoma (LuAd). Furthermore, we recently observed that aged Cldn18 ^-/- mice have increased propensity to develop LuAd. We now demonstrate that CLDN18.1 expression correlates inversely with promoter methylation and with LuAd patient mortality. In addition, when restored in LuAd cells that have lost expression, CLDN18.1 markedly attenuates malignant properties including xenograft tumor growth in vivo as well as cell proliferation, migration, invasion and anchorage-independent colony formation in vitro. Based on high throughput analyses of Cldn18 ^-/- murine lung alveolar epithelial type II cells, as well as CLDN18.1-repleted human LuAd cells, we hypothesized and subsequently confirmed by Western analysis that CLDN18.1 inhibits insulin-like growth factor-1 receptor (IGF-1R) and AKT phosphorylation. Consistent with recent data in Cldn18 ^-/- knockout mice, expression of CLDN18.1 in human LuAd cells also decreased expression of transcriptional co-activator with PDZ-binding motif (TAZ) and Yes-associated protein (YAP) and their target genes, contributing to its tumor suppressor activity. Moreover, analysis of LuAd cells in which YAP and/or TAZ are silenced with siRNA suggests that inhibition of TAZ, and possibly YAP, is also involved in CLDN18.1-mediated AKT inactivation. Taken together, these data indicate a tumor suppressor role for CLDN18.1 in LuAd mediated by a regulatory network that encompasses YAP/TAZ, IGF-1R and AKT signaling.

Estrogens and selective estrogen receptor modulators differentially antagonize Runx2 in ST2 mesenchymal progenitor cells

Estrogens attenuate bone turnover by inhibiting both osteoclasts and osteoblasts, in part through antagonizing Runx2. Apparently conflicting, stimulatory effects in osteoblast lineage cells, however, sway the balance between bone resorption and bone formation in favor of the latter. Consistent with this dualism, 17ß-estradiol (E2) both stimulates and inhibits Runx2 in a locus-specific manner, and here we provide evidence for such locus-specific regulation of Runx2 by E2 in vivo. We also demonstrate dual, negative and positive, regulation of Runx2-driven alkaline phosphatase (ALP) activity by increasing E2 concentrations in ST2 osteoblast progenitor cells. We further compared the effects of E2 to those of the Selective Estrogen Receptor Modulators (SERMs) raloxifene (ral) and lasofoxifene (las) and the phytoestrogen puerarin. We found that E2 at the physiological concentrations of 0.1-1 nM, as well as ral and las, but not puerarin, antagonize Runx2-driven ALP activity. At ≥10 nM, E2 and puerarin, but not ral or las, stimulate ALP relative to the activity measured at 0.1-1 nM. Contrasting the difference between E2 and SERMs in ST2 cells, they all shared a similar dose-response profile when inhibiting pre-osteoclast proliferation. That ral and las poorly mimic the locus- and concentration-dependent effects of E2 in mesenchymal progenitor cells may help explain their limited clinical efficacy.

Claudin-18-mediated YAP activity regulates lung stem and progenitor cell homeostasis and tumorigenesis

Claudins, the integral tight junction (TJ) proteins that regulate paracellular permeability and cell polarity, are frequently dysregulated in cancer; however, their role in neoplastic progression is unclear. Here, we demonstrated that knockout of Cldn18, a claudin family member highly expressed in lung alveolar epithelium, leads to lung enlargement, parenchymal expansion, increased abundance and proliferation of known distal lung progenitors, the alveolar epithelial type II (AT2) cells, activation of Yes-associated protein (YAP), increased organ size, and tumorigenesis in mice. Inhibition of YAP decreased proliferation and colony-forming efficiency (CFE) of Cldn18-/- AT2 cells and prevented increased lung size, while CLDN18 overexpression decreased YAP nuclear localization, cell proliferation, CFE, and YAP transcriptional activity. CLDN18 and YAP interacted and colocalized at cell-cell contacts, while loss of CLDN18 decreased YAP interaction with Hippo kinases p-LATS1/2. Additionally, Cldn18-/- mice had increased propensity to develop lung adenocarcinomas (LuAd) with age, and human LuAd showed stage-dependent reduction of CLDN18.1. These results establish CLDN18 as a regulator of YAP activity that serves to restrict organ size, progenitor cell proliferation, and tumorigenesis, and suggest a mechanism whereby TJ disruption may promote progenitor proliferation to enhance repair following injury.

Relationship between RUNX1 and AXIN1 in ER-negative versus ER-positive Breast Cancer

RUNX1 plays opposing roles in breast cancer: a tumor suppressor in estrogen receptor-positive (ER⁺) disease and an oncogenic role in ER-negative (ER^-) tumors. Potentially mediating the former, we have recently reported that RUNX1 prevents estrogen-driven suppression of the mRNA encoding the tumor suppressor AXIN1. Accordingly, AXIN1 protein expression was diminished upon RUNX1 silencing in ER⁺ breast cancer cells and was positively correlated with AXIN1 protein expression across tumors with high levels of ER. Here we report the surprising observation that RUNX1 and AXIN1 proteins are strongly correlated in ER^- tumors as well. However, this correlation is not attributable to regulation of AXIN1 by RUNX1 or vice versa. The unexpected correlation between RUNX1, playing an oncogenic role in ER^- breast cancer, and AXIN1, a well-established tumor suppressor hub, may be related to a high ratio between the expression of variant 2 and variant 1 (v2/v1) of AXIN1 in ER^- compared with ER⁺ breast cancer. Although both isoforms are similarly regulated by RUNX1 in estrogen-stimulated ER⁺ breast cancer cells, the higher v2/v1 ratio in ER^- disease is expected to weaken the tumor suppressor activity of AXIN1 in these tumors.

Glucocorticoids Hijack Runx2 to Stimulate Wif1 for Suppression of Osteoblast Growth and Differentiation

Inhibition of Runx2 is one of many mechanisms that suppress bone formation in glucocorticoid (GC)-induced osteoporosis (GIO). We profiled mRNA expression in ST2/Rx2(dox) cells after treatment with doxycycline (dox; to induce Runx2) and/or the synthetic GC dexamethasone (dex). As expected, dex typically antagonized Runx2-driven transcription. Select genes, however, were synergistic stimulated and this was confirmed by RT-qPCR. Among the genes synergistically stimulated by GCs and Runx2 was Wnt inhibitory Factor 1 (Wif1), and Wif1 protein was readily detectable in medium conditioned by cultures co-treated with dox and dex, but neither alone. Cooperation between Runx2 and GCs in stimulating Wif1 was also observed in primary preosteoblast cultures. GCs strongly inhibited dox-driven alkaline phosphatase (ALP) activity in control ST2/Rx2(dox) cells, but not in cells in which Wif1 was silenced. Unlike its anti-mitogenic activity in committed osteoblasts, induction of Runx2 transiently increased the percentage of cells in S-phase and accelerated proliferation in the ST2 mesenchymal pluripotent cell culture model. Furthermore, like the inhibition of Runx2-driven ALP activity, dex antagonized the transient mitogenic effect of Runx2 in ST2/Rx2(dox) cultures, and this inhibition eased upon Wif1 silencing. Plausibly, homeostatic feedback loops that rely on Runx2 activation to compensate for bone loss in GIO are thwarted, exacerbating disease progression through stimulation of Wif1. J. Cell. Physiol. 232: 145-153, 2017. © 2016 Wiley Periodicals, Inc.

RUNX1 prevents oestrogen-mediated AXIN1 suppression and β-catenin activation in ER-positive breast cancer

Recent high-throughput studies revealed recurrent RUNX1 mutations in breast cancer, specifically in oestrogen receptor-positive (ER⁺) tumours. However, mechanisms underlying the implied RUNX1-mediated tumour suppression remain elusive. Here, by depleting mammary epithelial cells of RUNX1 in vivo and in vitro, we demonstrate combinatorial regulation of AXIN1 by RUNX1 and oestrogen. RUNX1 and ER occupy adjacent elements in AXIN1's second intron, and RUNX1 antagonizes oestrogen-mediated AXIN1 suppression. Accordingly, RNA-seq and immunohistochemical analyses demonstrate an ER-dependent correlation between RUNX1 and AXIN1 in tumour biopsies. RUNX1 loss in ER⁺ mammary epithelial cells increases β-catenin, deregulates mitosis and stimulates cell proliferation and expression of stem cell markers. However, it does not stimulate LEF/TCF, c-Myc or CCND1, and it does not accelerate G1/S cell cycle phase transition. Finally, RUNX1 loss-mediated deregulation of β-catenin and mitosis is ameliorated by AXIN1 stabilization in vitro, highlighting AXIN1 as a potential target for the management of ER⁺ breast cancer.

The tumour suppressor RUNX1 is often lost or mutated in oestrogen receptor-positive breast cancers. In this study, the authors demonstrate that the loss of RUNX1 unleashes oestrogen-mediated inhibition of AXIN1, a negative regulator of β-catenin, resulting in β-catenin signalling-mediated cancer cell proliferation and mitosis deregulation.

Estrogens antagonize RUNX2-mediated osteoblast-driven osteoclastogenesis through regulating RANKL membrane association

In addition to its thoroughly investigated role in bone formation, the osteoblast master transcription factor RUNX2 also promotes osteoclastogenesis and bone resorption. Here we demonstrate that 17β-estradiol (E2), strongly inhibits RUNX2-mediated osteoblast-driven osteoclastogenesis in co-cultures. Towards deciphering the underlying mechanism, we induced premature expression of RUNX2 in primary murine pre-osteoblasts, which resulted in robust differentiation of co-cultured splenocytes into mature osteoclasts. This was attributable to RUNX2-mediated increase in RANKL secretion, determined by ELISA, as well as to RUNX2-mediated increase in RANKL association with the osteoblast membrane, demonstrated using confocal fluorescence microscopy. The increased association with the osteoblast membrane was recapitulated by transiently expressed GFP-RANKL. E2 abolished the RUNX2-mediated increase in membrane-associated RANKL and GFP-RANKL, as well as the concomitant osteoclastogenesis. RUNX2-mediated RANKL cellular redistribution was attributable in part to a decrease in Opg expression, but E2 did not influence Opg expression either in the presence or absence of RUNX2. Diminution of RUNX2-mediated osteoclastogenesis by E2 occurred regardless of whether the pre-osteoclasts were derived from wild type or estrogen receptor alpha (ERα)-knockout mice, suggesting that activated ERα inhibited osteoblast-driven osteoclastogenesis by acting in osteoblasts, possibly targeting RUNX2. Indeed, microarray analysis demonstrated global attenuation of the RUNX2 response by E2, including abrogation of Pstpip2 expression, which likely plays a critical role in membrane trafficking. Finally, the selective ER modulators (SERMs) tamoxifen and raloxifene mimicked E2 in abrogating the stimulatory effect of osteoblastic RUNX2 on osteoclast differentiation in the co-culture assay. Thus, E2 antagonizes RUNX2-mediated RANKL trafficking and subsequent osteoclastogenesis. Targeting RUNX2 and/or downstream mechanisms that regulate RANKL trafficking may lead to the development of improved SERMs and possibly non-hormonal therapeutic approaches to high turnover bone disease.

Differential effects of RUNX2 on the androgen receptor in prostate cancer: synergistic stimulation of a gene set exemplified by SNAI2 and subsequent invasiveness

Changes to androgen signaling during prostate carcinogenesis are associated with both inhibition of cellular differentiation and promotion of malignant phenotypes. The androgen receptor (AR)-binding transcription factor RUNX2 has been linked to prostate cancer progression but the underlying mechanisms have not been fully defined. In this study, we investigated the genome-wide influence of RUNX2 on androgen-induced gene expression and AR DNA binding in prostate cancer cells. RUNX2 inhibited the androgen response partly by promoting the dissociation of AR from its target genes such as the tumor suppressor NKX3-1. However, AR activity persists in the presence of RUNX2 at other AR target genes, some of which are cooperatively stimulated by androgen and RUNX2 signaling. These genes are associated with putative enhancers co-occupied by AR and RUNX2. One such gene, the invasion-promoting Snail family transcription factor SNAI2, was co-activated by AR and RUNX2. Indeed, these two transcription factors together, but neither alone stimulated prostate cancer cell invasiveness, which could be abolished by SNAI2 silencing. Furthermore, an immunohistochemical analysis of SNAI2 in archived primary prostate cancer specimens revealed a correlation with the RUNX2 histoscore, and simultaneous strong staining for SNAI2, RUNX2, and AR (but not any pair alone) was associated with disease recurrence. Overall, our findings suggest cooperation between AR and RUNX in the stimulation of oncogenes such as SNAI2, which might be targeted for individualized prostate cancer therapy.

Prolactin-induced protein (PIP) regulates proliferation of luminal A type breast cancer cells in an estrogen-independent manner

Prolactin-induced Protein (PIP), an aspartyl protease unessential for normal mammalian cell function, is required for the proliferation and invasion of some breast cancer (BCa) cell types. Because PIP expression is particularly high in the Luminal A BCa subtype, we investigated the roles of PIP in the related T47D BCa cell line. Nucleic acid and antibody arrays were employed to screen effects of PIP silencing on global gene expression and activation of receptor tyrosine kinases (RTKs), respectively. Expression of PIP-stimulated genes, as defined in the T47D cell culture model, was well correlated with the expression of PIP itself across a cohort of 557 mRNA profiles of diverse BCa tumors, and bioinformatics analysis revealed cJUN and cMYC as major nodes in the PIP-dependent gene network. Among 71 RTKs tested, PIP silencing resulted in decreased phosphorylation of focal adhesion kinase (FAK), ephrin B3 (EphB3), FYN, and hemopoietic cell kinase (HCK). Ablation of PIP also abrogated serum-induced activation of the downstream serine/threonine kinases AKT, ERK1/2, and JNK1. Consistent with these results, PIP-depleted cells exhibited defects in adhesion to fibronectin, cytoskeletal stress fiber assembly and protein secretion. In addition, PIP silencing abrogated the mitogenic response of T47D BCa cells to estradiol (E2). The dependence of BCa cell proliferation was unrelated, however, to estrogen signaling because: 1) PIP silencing did not affect the transcriptional response of estrogen target genes to hormone treatment, and 2) PIP was required for the proliferation of tamoxifen-resistant BCa cells. Pharmacological inhibition of PIP may therefore serve the bases for both augmentation of existing therapies for hormone-dependent tumors and the development of novel therapeutic approaches for hormone-resistant BCa.

Runx2 controls a feed-forward loop between androgen and prolactin-induced protein (PIP) in stimulating T47D cell proliferation

Prolactin-Induced Protein (PIP) is a small polypeptide expressed by breast and prostate cancer (BCa, PCa) cells. However, both the regulation of PIP expression and its function in cancer cells are poorly understood. Using BCa and PCa cells, we found that Runx2, a pro-metastatic transcription factor, functionally interacts with the Androgen Receptor (AR) to regulate PIP expression. Runx2 expression in C4-2B PCa cells synergized with AR to promote PIP expression, whereas its knockdown in T47D BCa cells abrogated basal as well as hormone stimulated PIP expression. Chromatin immunoprecipitation (ChIP) assays showed that Runx2 and AR co-occupied an enhancer element located ∼11 kb upstream of the PIP open reading frame, and that Runx2 facilitated AR recruitment to the enhancer. PIP knockdown in T47D cells compromised DHT-stimulated expression of multiple AR target genes including PSA, FKBP5, FASN, and SGK1. The inhibition of AR activity due to loss of PIP was attributable at least in part to abrogation of its nuclear translocation. PIP knockdown also suppressed T47D cell proliferation driven by either serum growth factors or dihydrotestosterone (DHT). Our data suggest that Runx2 controls a positive feedback loop between androgen signaling and PIP, and pharmacological inhibition of PIP may be useful to treat PIP positive tumors.

Regulation of breast cancer metastasis by Runx2 and estrogen signaling: the role of SNAI2

Introduction

In contrast to its role in breast cancer (BCa) initiation, estrogen signaling has a protective effect in later stages, where estrogen receptor (ER)α loss associates with aggressive metastatic disease. We asked whether the beneficial effect of estrogen signaling in late-stage BCa is attributable to the recently reported estrogen-mediated antagonism of the pro-metastatic transcription factor Runx2.

Methods

MCF7/Rx2^dox breast cancer cells were engineered with a lentivirus expressing Runx2 in response to doxycycline (dox). Cells treated with dox and/or estradiol (E2) were subjected to genome-wide expression profiling, RT-qPCR analysis of specific genes, and Matrigel™ invasion assays. Knockdown of genes of interest was performed using lentiviruses expressing appropriate shRNAs, either constitutively or in response to dox. Gene expression in BCa tumors was investigated using a cohort of 557 patients compiled from publicly available datasets. Association of gene expression with clinical metastasis was assessed by dichotomizing patients into those expressing genes of interest at either high or low levels, and comparing the respective Kaplan-Meier curves of metastasis-free survival.

Results

Runx2 induced epithelial-mesenchymal transition (EMT) evidenced by acquisition of a fibroblastic morphology, decreased expression of E-cadherin, increased expression of vimentin and invasiveness. Runx2 stimulated SNAI2 expression in a WNT- and transforming growth factor (TGF)β-dependent manner, and knockdown of SNAI2 abrogated the pro-metastatic activities of Runx2. E2 antagonized the pro-metastatic activities of Runx2, including SNAI2 upregulation. In primary BCa tumors, Runx2 activity, SNAI2 expression, and metastasis were positively correlated, and SNAI2 expression was negatively correlated with ERα. However, the negative correlation between SNAI2 and ERα in bone-seeking BCa cells was weaker than the respective negative correlation in tumors seeking lung. Furthermore, the absence of ERα in primary tumors was associated with lung- and brain- but not with bone metastasis, and tumor biopsies from bone metastatic sites displayed the unusual combination of high Runx2/SNAI2 and high ERα expression.

Conclusions

E2 antagonizes Runx2-induced EMT and invasiveness of BCa cells, partly through attenuating expression of SNAI2, a Runx2 target required for mediating its pro-metastatic property. That ERα loss promotes non-osseous metastasis by unleashing Runx2/SNAI2 is supported by the negative correlation observed in corresponding tumors. Unknown mechanisms in bone-seeking BCa allow high Runx2/SNAI2 expression despite high ERα level

Opposing effects of Runx2 and estradiol on breast cancer cell proliferation: in vitro identification of reciprocally regulated gene signature related to clinical letrozole responsiveness

PURPOSE

To assess the clinical significance of the interaction between estrogen and Runx2 signaling, previously shown in vitro.

EXPERIMENTAL DESIGN

MCF7/Rx2(dox) breast cancer cells were treated with estradiol and/or doxycycline to induce Runx2, and global gene expression was profiled to define genes regulated by estradiol, Runx2, or both. Anchorage-independent growth was assessed by soft-agar colony formation assays. Expression of gene sets defined using the MCF7/Rx2(dox) system was analyzed in pre- and on-treatment biopsies from hormone receptor-positive patients undergoing neoadjuvant letrozole treatment in two independent studies, and short-term changes in gene expression were correlated with tumor size reduction or Ki67 index at surgery.

RESULTS

Reflecting its oncogenic property, estradiol strongly promoted soft-agar colony formation, whereas Runx2 blocked this process suggesting tumor suppressor property. Transcriptome analysis of MCF7/Rx2(dox) cells treated with estradiol and/or doxycycline showed reciprocal attenuation of Runx2 and estrogen signaling. Correspondingly in breast cancer tumors, expression of estradiol- and Runx2-regulated genes was inversely correlated, and letrozole increased expression of Runx2-stimulated genes, as defined in the MCF7/Rx2(dox) model. Of particular interest was a gene set upregulated by estradiol and downregulated by Runx2 in vitro; its short-term response to letrozole treatment associated with tumor size reduction and Ki67 index at surgery better than other estradiol-regulated gene sets.

CONCLUSION

This work provides clinical evidence for the importance of antagonism between Runx2 and E2 signaling in breast cancer. Likely sensing the tension between them, letrozole responsiveness of a genomic node, positively regulated by estradiol and negatively regulated by Runx2 in vitro, best correlated with the clinical efficacy of letrozole treatment.

Segmental duplication as one of the driving forces underlying the diversity of the human immunoglobulin heavy chain variable gene region

Background

Segmental duplication and deletion were implicated for a region containing the human immunoglobulin heavy chain variable (IGHV) gene segments, 1.9III/hv3005 (possible allelic variants of IGHV3-30) and hv3019b9 (a possible allelic variant of IGHV3-33). However, very little is known about the ranges of the duplication and the polymorphic region. This is mainly because of the difficulty associated with distinguishing between allelic and paralogous sequences in the IGHV region containing extensive repetitive sequences. Inability to separate the two parental haploid genomes in the subjects is another serious barrier. To address these issues, unique DNA sequence tags evenly distributed within and flanking the duplicated region implicated by the previous studies were selected. The selected tags in single sperm from six unrelated healthy donors were amplified by multiplex PCR followed by microarray detection. In this way, individual haplotypes of different parental origins in the sperm donors could be analyzed separately and precisely. The identified polymorphic region was further analyzed at the nucleotide sequence level using sequences from the three human genomic sequence assemblies in the database.

Results

A large polymorphic region was identified using the selected sequence tags. Four of the 12 haplotypes were shown to contain consecutively undetectable tags spanning in a variable range. Detailed analysis of sequences from the genomic sequence assemblies revealed two large duplicate sequence blocks of 24,696 bp and 24,387 bp, respectively, and an incomplete copy of 961 bp in this region. It contains up to 13 IGHV gene segments depending on haplotypes. A polymorphic region was found to be located within the duplicated blocks. The variants of this polymorphism unusually diverged at the nucleotide sequence level and in IGHV gene segment number, composition and organization, indicating a limited selection pressure in general. However, the divergence level within the gene segments is significantly different from that in the intergenic regions indicating that these regions may have been subject to different selection pressures and that the IGHV gene segments in this region are functionally important.

Conclusions

Non-reciprocal genetic rearrangements associated with large duplicate sequence blocks could substantially contribute to the IGHV region diversity. Since the resulting polymorphisms may affect the number, composition and organization of the gene segments in this region, it may have significant impact on the function of the IGHV gene segment repertoire, antibody diversity, and therefore, the immune system. Because one of the gene segments, 3-30 (1.9III), is associated with autoimmune diseases, it could be of diagnostic significance to learn about the variants in the haplotypes by using the multiplex haplotype analysis system used in the present study with DNA sequence tags specific for the variants of all gene segments in this region.

Generation of neural crest progenitors from human embryonic stem cells

The neural crest (NC) is a transient population of multipotent progenitors arising at the lateral edge of the neural plate in vertebrate embryos, which then migrate throughout the body to generate diverse array of tissues such as the peripheral nervous system, skin melanocytes, and craniofacial cartilage, bone and teeth. The transient nature of neural crest stem cells make extremely challenging to study the biology of these important cells. In humans induction and differentiation of embryonic NC occurs very early, within a few weeks of fertilization giving rise to technical and ethical issues surrounding isolation of early embryonic tissues and therefore severely limiting the study of human NC cells. For that reason our current knowledge of the biology of NC mostly derives from the studies of lower organisms. Recent progress in human embryonic stem cell research provides a unique opportunity for generation of a useful source of cells for basic developmental studies. The development of cost-effective, time and labor efficient improved differentiation protocols for the production of human NC cells is a critical step toward a better understanding of NC biology.

Evidence for natural selection on leukocyte immunoglobulin-like receptors for HLA class I in Northeast Asians

Human leukocyte antigen (HLA) plays a critical role in innate and adaptive immunity and is a well-known example of genes under natural selection. However, the genetic aspect of receptors recognizing HLA molecules has not yet been fully elucidated. Leukocyte immunoglobulin (Ig)-like receptors (LILRs) are a family of HLA class I-recognizing receptors comprising activating and inhibitory forms. We previously reported that the allele frequency of the 6.7 kb LILRA3 deletion is extremely high (71%) in the Japanese population, and we identified premature termination codon (PTC)-containing alleles. In this study, we observed a wide distribution of the high deletion frequency in Northeast Asians (84% in Korean Chinese, 79% in Man Chinese, 56% in Mongolian, and 76% in Buryat populations). Genotyping of the four HapMap populations revealed that LILRA3 alleles were in strong linkage disequilibrium with LILRB2 alleles in Northeast Asians. In addition, PTC-containing LILRA3 alleles were detected in Northeast Asians but not in non-Northeast Asians. Furthermore, flow-cytometric analysis revealed that the LILRB2 allele frequent in Northeast Asians was significantly associated with low levels of expression. F(ST) and extended-haplotype-homozygosity analysis for the HapMap populations provided evidence of positive selection acting on the LILRA3 and LILRB2 loci. Taken together, our results suggest that both the nonfunctional LILRA3 alleles and the low-expressing LILRB2 alleles identified in our study have increased in Northeast Asians because of natural selection. Our findings, therefore, lead us to speculate that not only HLA class I ligands but also their receptors might be sensitive to the local environment.

Laser-assisted microdissection (LAM) in developmental biology

The analysis of gene expression in developing organs is a valuable tool for the assessment of genetic fingerprints during the various stages of differentiation. Complex processes in developing tissues are particularly difficult to understand in terms of biochemical phenomena. Laser-assisted microdissection (LAM) allows the efficient and precise capture of cells or groups of cells from developing tissues in sufficient quantities and within the context of time and space to permit the subsequent molecular characterization of the targeted tissue. The technique development has dramatically increased the ease of isolating specific cells which, together with progress in tissue preparation and microextraction protocols, allows for broad-range down-stream applications in the fields of genomics, transcriptomics and proteomics. This review gives an overview of the LAM technology and its application in developmental biology.

Expression profiling of BEN regulated genes in mouse embryonic fibroblasts

BEN is a member of the TFII-I family of helix-loop-helix transcription factors. Both TFII-I and BEN are involved in gene regulation through interactions with tissue-specific transcription factors and chromatin remodeling complexes. Identification of the downstream target genes of TFII-I proteins is critical in delineating the regulatory effects of these proteins. In this study, we conducted a microarray analysis to determine gene expression alterations following the overexpression of BEN in primary mouse embryonic fibroblasts (MEFs). We found the BEN-dependent modulation in the expression of large groups of genes representing a wide variety of functional categories including genes important in the immune response, cell cycle, transcriptional regulation and cell signaling. A set of genes identified by the microarray analysis was validated by independent real-time PCR analysis. Among upregulated genes were Shrm, Tgfb2, Ube2l6, G1p2, Ccl7 while downregulated genes were Folr1, Tgfbr2, Csrp2, and Dlk1. These results support a versatile function of TFII-I proteins in vertebrate physiology and lead to an increased understanding of the BEN-dependent molecular events.

Gene expression analysis of TFII-I modulated genes in mouse embryonic fibroblasts

TFII-I is a founding member of a family of helix-loop-helix transcription factors involved in modulation of genes through interaction with various nuclear factors and chromatin remodeling complexes. Recent studies indicate that TFII-I performs important function in cell physiology and mouse embryogenesis. In order to understand its molecular role, TFII-I was overexpressed in primary mouse embryonic fibroblasts (MEFs) and alterations in gene expression were monitored with a mouse 16 K oligonucleotide microarray. These studies allowed us to identify genes that lie downstream of TFII-I-dependent pathways. Among the modulated candidates were genes involved in the immunity response, catalytic activity, signaling pathways and transcriptional regulation. Expression of several candidates including those for the interferon-stimulated protein (G1p2), small inducible cytokine A7 (Ccl7), ubiquitin-conjugating enzyme 8 (Ube2l6), cysteine-rich protein (Csrp2) and Drosophila delta-like 1 homolog (Dlk1) were confirmed by real-time PCR. The obtained results suggest that TFII-I participates in multiple signaling and regulatory pathways in MEFs.

Strong correlation between meiotic crossovers and haplotype structure in a 2.5-Mb region on the long arm of chromosome 21

Although the haplotype structure of the human genome has been studied in great detail, very little is known about the mechanisms underlying its formation. To investigate the role of meiotic recombination on haplotype block formation, single nucleotide polymorphisms were selected at a high density from a 2.5-Mb region of human chromosome 21. Direct analysis of meiotic recombination by high-throughput multiplex genotyping of 662 single sperm identifies 41 recombinants. The crossovers were nonrandomly distributed within 16 small areas. All, except one, of these crossovers fall in areas where the haplotype structure exhibits breakdown, displaying a strong statistically positive association between crossovers and haplotype block breaks. The data also indicate a particular clustered distribution of recombination hotspots within the region. This finding supports the hypothesis that meiotic recombination makes a primary contribution to haplotype block formation in the human genome.

Genetic diversity of the human immunoglobulin heavy chain VH region

The human immunoglobulin heavy chain VH region is one of the most complex regions in the human genome. The high level of diversity of this region has been shown by a number of studies. However, because of the limitations of the conventional experimental methods, it has been difficult to learn the extent of the diversity and the underlying mechanisms. This review describes a number of new genetic approaches developed in the authors' laboratory. By using these approaches, significant progress has been made in assigning different VH sequences to their respective loci, in learning the diversity of gene segment number and composition among the VH haplotypes, and in learning VH gene segment organization in individual haplotypes. Information obtained toward this direction could help in understanding the mechanisms underlying VH region diversity and the biological impact of the VH region diversity.

Wide distribution of the MICA-MICB null haplotype in East Asians

Stress-inducible MICA (MHC class I chain-related A) is known to bind to NKG2D, which is one of the natural killer (NK) cell receptors, and plays a role in immune surveillance. We have reported that a MICA-MICB null haplotype is in linkage disequilibrium with HLA-B*4801 in the Japanese population. In the haplotype, an approximately 100-kb deletion, including the entire MICA gene, was observed and MICB possessed a premature stop codon. In this study, a multiplex polymerase chain reaction (PCR) method was developed for detecting the MICA deletion. MICB alleles were typed by PCR-single-strand conformation polymorphism (SSCP) method and direct sequencing. The frequency of the MICA-MICB null haplotype was 3.7% on the average, and was strongly associated with HLA-B48 in seven East Asian populations. It was presumed that the stop codon of MICB gene generated after the large-scale deletion. The wide distribution of the null haplotype at polymorphic frequencies suggests that the haplotype has been conservatively maintained because of some selective advantage.

Diversity of HLA-B17 alleles and haplotypes in East Asians and a novel Cw6 allele (Cw*0604) associated with B*5701

The distribution of HLA-B17 alleles and their association with HLA-A, -C and -DRB1 alleles were investigated in seven East Asian populations Japanese, South Korean, Chinese-Korean, Man, Northern Han, Mongolian and Buryat populations). The B17 alleles were identified from genomic DNA using group-specific polymerase chain reaction (PCR) followed by hybridization with sequence-specific oligonucleotide probes (SSOP). In all of these East Asian populations, except Japanese and Chinese-Koreans, B*5701 was detected and strongly associated with A*0101, Cw*0602 and DRB1*0701. In contrast, B*5801 was detected in all the seven populations and strongly associated with A*3303, Cw*0302, DRB1*0301 and DRB1*1302. The A*3303-Cw*0302-B*5801-DRB1*1302 haplotype was observed in South Korean, Chinese-Korean, Buryat and Japanese populations, while A*3303-Cw*0302-B*5801-DRB1*0301 was predominantly observed in the Mongolian population. A similar haplotype, A*0101-Cw*0302-B*5801-DRB1*1302, was observed in the Buryat population. A novel Cw6 allele, Cw*0604, was identified in the Man population. This Cw allele was observed on the haplotype A*0101-B*5701-DRB1*0701. Thus, we confirmed, at the sequence level, that the common haplotypes carrying B*5701 and B*5801 have been conserved and shared in East Asian populations.

Diversity of HLA-B61 alleles and haplotypes in East Asians and Spanish Gypsies

The distribution of HLA-B61 alleles and their association with HLA-C and DRB1 alleles were investigated in six East Asian populations (South Korean, Chinese Korean, Man (Manchu), Northern Han, Mongolian and Buryat) and Spanish Gypsies and compared to our previous report on the Japanese population. The alleles were identified using a group-specific polymerase chain reaction (PCR) and genomic DNA followed by hybridization with sequence-specific oligonucleotide probes (SSOP). Both HLA-B*4002 and B*4006 were commonly detected in the South Korean, Chinese Korean, Man, Northern Han and Japanese populations, while HLA-B*4002 was predominant in the Mongolian and Buryat populations. Strong associations of B*4002 with Cw*0304 and of B*4006 with Cw*0801 were commonly observed in these East Asian populations. In contrast, in Spanish Gypsies, only HLA-B*4006 was found and the allele exhibited a strong association with Cw*1502. HLA-B*4003 was also identified in the South Korean, Chinese Korean, Northern Han, Mongolian and Japanese populations at relatively low frequencies, and exhibited an association with Cw*0304. Moreover, the association of these B61 alleles with the DRB1 alleles revealed considerable diversity among the different populations. HLA-B*4004 and B*4009 were not observed in these populations. Consequently, the frequencies of the B61 alleles varied among the different East Asian populations, but the individual B61 alleles were carried by specific haplotypes often regardless of the ethnic differences.

The HLA system in the population of Mongolia

The frequencies of HLA-A, B, C and DRB1 alleles and haplotype frequencies for HLA-A, B, C and A, B, DR loci were studied in 201 healthy unrelated Khalkha-Mongolians. The most common class I antigens were A24 (25.8%), A2 (23.4%), B61 (12.6%), B51 (8.5%), Cw10 (14.2%), Cw9 (14.1%), Cw7 (13.3%), Cw1 (11.8%) and Cw6 (10.2%). A total of 35 DRB1 alleles were identified in this group of samples. The most frequent DRB1 allele was DRB1*0301 (11.1%), followed by DRB1*0701 (9.7%) and DRB1*1101 (8.5%). One novel DRB1 allele (14MV) and three rare types, DRB1*1111, DRB1*1504 and DRB1*1412, recently described in Jews, the Dai minority of China and Japanese, respectively, were identified in Mongolians. The phylogenetic tree constructed by UPG method put Mongolians in the Northeast Asian cluster. A comparison of three locus haplotype distributions with world populations, revealed that Mongolians share several characteristics in common with other Mongoloids as well as with Caucasoids and Amerindians. The most common A, B and DR haplotype in Mongolians, A33-B58-DR3, was shared with Thai, Thai Chinese and Singapore Chinese. These data support that unique genetic background of Mongolians played a major role in ethnic formation and differentiation of Mongoloid populations.