Prostate cancer susceptibility gene HIST1H1A is a modulator of androgen receptor signaling and epithelial to mesenchymal transition.

In 2018, approximately 165,000 new prostate cancer (PC) cases will be diagnosed, and over 29,000 men will succumb to PC in the U.S. alone. The means of assessing outcome in the clinic are inaccurate, and there is a pressing need to more precisely identify men at risk of aggressive PC. We previously identified HIST1H1A as a susceptibility gene for aggressive PC. HIST1H1A encodes H1.1, a member of the linker histone family that is involved in chromatin organization and compaction. To understand the molecular basis of aggressive PC, we have characterized how germline variation modulates susceptibility to neuroendocrine differentiation, which is a form of aggressive PC. Immunohistochemistry studies revealed that HIST1H1A is over-expressed in normal human prostate tissue compared to prostate adenocarcinoma. Functional characterization of HIST1H1A in prostate LNCaP cells indicated that HIST1HA over-expression increased cell growth, as well as the expression of neuroendocrine and epithelial-to-mesenchymal markers in vitro. Assay for Transposase-Accessible Chromatin (ATAC-seq), which is used to assess chromatin compaction and thus the transcriptional availability of individual genomic regions, demonstrated that H1.1 plays a prominent role in modulating Wnt signaling pathway genes, which are implicated in prostate tumorigenesis. These results demonstrate that HIST1H1A is a modulator of aggressive PC susceptibility.

Modifier locus mapping of a transgenic F2 mouse population identifies CCDC115 as a novel aggressive prostate cancer modifier gene in humans.

BACKGROUND: It is well known that development of prostate cancer (PC) can be attributed to somatic mutations of the genome, acquired within proto-oncogenes or tumor-suppressor genes. What is less well understood is how germline variation contributes to disease aggressiveness in PC patients. To map germline modifiers of aggressive neuroendocrine PC, we generated a genetically diverse F2 intercross population using the transgenic TRAMP mouse model and the wild-derived WSB/EiJ (WSB) strain. The relevance of germline modifiers of aggressive PC identified in these mice was extensively correlated in human PC datasets and functionally validated in cell lines. RESULTS: Aggressive PC traits were quantified in a population of 30 week old (TRAMP x WSB) F2 mice (n = 307). Correlation of germline genotype with aggressive disease phenotype revealed seven modifier loci that were significantly associated with aggressive disease. RNA-seq were analyzed using cis-eQTL and trait correlation analyses to identify candidate genes within each of these loci. Analysis of 92 (TRAMP x WSB) F2 prostates revealed 25 candidate genes that harbored both a significant cis-eQTL and mRNA expression correlations with an aggressive PC trait. We further delineated these candidate genes based on their clinical relevance, by interrogating human PC GWAS and PC tumor gene expression datasets. We identified four genes (CCDC115, DNAJC10, RNF149, and STYXL1), which encompassed all of the following characteristics: 1) one or more germline variants associated with aggressive PC traits; 2) differential mRNA levels associated with aggressive PC traits; and 3) differential mRNA expression between normal and tumor tissue. Functional validation studies of these four genes using the human LNCaP prostate adenocarcinoma cell line revealed ectopic overexpression of CCDC115 can significantly impede cell growth in vitro and tumor growth in vivo. Furthermore, CCDC115 human prostate tumor expression was associated with better survival outcomes. CONCLUSION: We have demonstrated how modifier locus mapping in mouse models of PC, coupled with in silico analyses of human PC datasets, can reveal novel germline modifier genes of aggressive PC. We have also characterized CCDC115 as being associated with less aggressive PC in humans, placing it as a potential prognostic marker of aggressive PC.

Mapping Complex Traits in a Diversity Outbred F1 Mouse Population Identifies Germline Modifiers of Metastasis in Human Prostate Cancer.

It is unclear how standing genetic variation affects the prognosis of prostate cancer patients. To provide one controlled answer to this problem, we crossed a dominant, penetrant mouse model of prostate cancer to Diversity Outbred mice, a collection of animals that carries over 40 million SNPs. Integration of disease phenotype and SNP variation data in 493 F1 males identified a metastasis modifier locus on Chromosome 8 (LOD = 8.42); further analysis identified the genes Rwdd4, Cenpu, and Casp3 as functional effectors of this locus. Accordingly, analysis of over 5,300 prostate cancer patient samples revealed correlations between the presence of genetic variants at these loci, their expression levels, cancer aggressiveness, and patient survival. We also observed that ectopic overexpression of RWDD4 and CENPU increased the aggressiveness of two human prostate cancer cell lines. In aggregate, our approach demonstrates how well-characterized genetic variation in mice can be harnessed in conjunction with systems genetics approaches to identify and characterize germline modifiers of human disease processes.

SIN3A and SIN3B differentially regulate breast cancer metastasis.

SIN3 corepressor complexes play important roles in both normal development and breast cancer. Mammalian cells have two paralogs of SIN3 (SIN3A and SIN3B) that are encoded by distinct genes and have unique functions in many developmental processes. However, specific roles for SIN3A and SIN3B in breast cancer progression have not been characterized. We generated stable knockdown cells of SIN3 paralogs individually and in combination using three non-overlapping shRNA. Stable knockdown of SIN3B caused a significant decrease in transwell invasion through Matrigel and decreased the number of invasive colonies when grown in a 3D extracellular matrix. Conversely, stable knockdown of SIN3A significantly increased transwell invasion and increased the number of invasive colonies. These results were corroborated in vivo in which SIN3B knockdown significantly decreased and SIN3A knockdown increased experimental lung metastases. RNA sequencing was used to identify unique targets and biological pathways that were altered upon knockdown of SIN3A compared to SIN3B. Additionally, we analyzed microarray data sets to identify correlations of SIN3A and SIN3B expression with survival in patients with breast cancer. These data sets indicated that high mRNA expression of SIN3A as well as low mRNA expression of SIN3B correlates with longer relapse free survival specifically in patients with triple negative breast cancer which corresponds with our in vitro and in vivo data. These results demonstrate key functional differences between SIN3 paralogs in regulating the process of breast cancer metastasis and suggest metastasis suppressive roles of SIN3A and metastasis promoting roles of SIN3B.

GNL3 and SKA3 are novel prostate cancer metastasis susceptibility genes.

Prostate cancer (PC) is very common in developed countries. However, the molecular determinants of PC metastasis are unclear. Previously, we reported that germline variation influences metastasis in the C57BL/6-Tg(TRAMP)8247Ng/J (TRAMP) mouse model of PC. These mice develop prostate tumors similar to a subset of poor outcome, treatment-associated human PC tumors. Here, we used TRAMP mice to nominate candidate genes and validate their role in aggressive human PC in PC datasets and cell lines. Candidate metastasis susceptibility genes were identified through quantitative trait locus (QTL) mapping in 201 (TRAMP × PWK/PhJ) F2 males. Two metastasis-associated QTLs were identified; one on chromosome 12 (LOD = 5.86), and one on chromosome 14 (LOD = 4.41). Correlation analysis using microarray data from (TRAMP × PWK/PhJ) F2 prostate tumors identified 35 metastasis-associated transcripts within the two loci. The role of these genes in susceptibility to aggressive human PC was determined through in silico analysis using multiple datasets. First, analysis of candidate gene expression in two human PC datasets demonstrated that five candidate genes were associated with an increased risk of aggressive disease and lower disease-free survival. Second, four of these genes (GNL3, MAT1A, SKA3, and ZMYM5) harbored SNPs associated with aggressive tumorigenesis in the PLCO/CGEMS GWAS of 1172 PC patients. Finally, over-expression of GNL3 and SKA3 in the PC-3 human PC cell line decreased in vitro cell migration and invasion. This novel approach demonstrates how mouse models can be used to identify metastasis susceptibility genes, and gives new insight into the molecular mechanisms of fatal PC.

Necdin is a breast cancer metastasis suppressor that regulates the transcription of c-Myc.

Metastasis is the primary cause of death in breast cancer. Earlier studies using a mammary tumorigenesis mouse model identified Necdin (Ndn)as a germline modifier of metastasis. Differential expression of Ndn induces a gene-expression signature that predicts prognosis in human breast cancer. Additionally, a non-synonymous germline single nucleotide polymorphism (T50C; V17A) in Ndn distinguishes mouse strains with differing metastatic capacities. To better understand how hereditary factors influence metastasis in breast cancer, we characterized NDN-mediated transcription. Haplotype analysis in a well-characterized breast cancer cohort revealed that NDN germline variation is associated with both NDN expression levels and patient outcome. To examine the role of NDN in mammary tumor metastasis and transcriptional regulation, mouse mammary tumor cell lines stably over-expressing either the wildtype 50T or variant 50C Ndn allele were generated. Cells over-expressing Ndn 50T, but not Ndn 50C, exhibited significant decrease in cell invasiveness and pulmonary metastases compared to control cells. Transcriptome analyses identified a 71-gene expression signature that distinguishes cells over-expressing the two Ndn allelic variants. Furthermore, ChIP assays revealed c-Myc, a target gene of NDN, to be differentially regulated by the allelic variants. These data demonstrate that NDN and the T50C allele regulate gene expression and metastasis efficiency.

Phenotype-driven chemical screening in zebrafish for compounds that inhibit collective cell migration identifies multiple pathways potentially involved in metastatic invasion.

In the last decade, high-throughput chemical screening has become the dominant approach for discovering novel compounds with therapeutic properties. Automated screening using in vitro or cultured cell assays have yielded thousands of candidate drugs for a variety of biological targets, but these approaches have not resulted in an increase in drug discovery despite major increases in expenditures. In contrast, phenotype-driven screens have shown a much stronger success rate, which is why we developed an in vivo assay using transgenic zebrafish with a GFP-marked migrating posterior lateral line primordium (PLLp) to identify compounds that influence collective cell migration. We then conducted a high-throughput screen using a compound library of 2160 annotated bioactive synthetic compounds and 800 natural products to identify molecules that block normal PLLp migration. We identified 165 compounds that interfere with primordium migration without overt toxicity in vivo. Selected compounds were confirmed in their migration-blocking activity by using additional assays for cell migration. We then proved the screen to be successful in identifying anti-metastatic compounds active in vivo by performing orthotopic tumor implantation assays in mice. We demonstrated that the Src inhibitor SU6656, identified in our screen, can be used to suppress the metastatic capacity of a highly aggressive mammary tumor cell line. Finally, we used CRISPR/Cas9-targeted mutagenesis in zebrafish to genetically validate predicted targets of compounds. This approach demonstrates that the migrating PLLp in zebrafish can be used for large-scale, high-throughput screening for compounds that inhibit collective cell migration and, potentially, anti-metastatic compounds.

A systems genetics approach identifies CXCL14, ITGAX, and LPCAT2 as novel aggressive prostate cancer susceptibility genes.

Although prostate cancer typically runs an indolent course, a subset of men develop aggressive, fatal forms of this disease. We hypothesize that germline variation modulates susceptibility to aggressive prostate cancer. The goal of this work is to identify susceptibility genes using the C57BL/6-Tg(TRAMP)8247Ng/J (TRAMP) mouse model of neuroendocrine prostate cancer. Quantitative trait locus (QTL) mapping was performed in transgene-positive (TRAMPxNOD/ShiLtJ) F2 intercross males (n = 228), which facilitated identification of 11 loci associated with aggressive disease development. Microarray data derived from 126 (TRAMPxNOD/ShiLtJ) F2 primary tumors were used to prioritize candidate genes within QTLs, with candidate genes deemed as being high priority when possessing both high levels of expression-trait correlation and a proximal expression QTL. This process enabled the identification of 35 aggressive prostate tumorigenesis candidate genes. The role of these genes in aggressive forms of human prostate cancer was investigated using two concurrent approaches. First, logistic regression analysis in two human prostate gene expression datasets revealed that expression levels of five genes (CXCL14, ITGAX, LPCAT2, RNASEH2A, and ZNF322) were positively correlated with aggressive prostate cancer and two genes (CCL19 and HIST1H1A) were protective for aggressive prostate cancer. Higher than average levels of expression of the five genes that were positively correlated with aggressive disease were consistently associated with patient outcome in both human prostate cancer tumor gene expression datasets. Second, three of these five genes (CXCL14, ITGAX, and LPCAT2) harbored polymorphisms associated with aggressive disease development in a human GWAS cohort consisting of 1,172 prostate cancer patients. This study is the first example of using a systems genetics approach to successfully identify novel susceptibility genes for aggressive prostate cancer. Such approaches will facilitate the identification of novel germline factors driving aggressive disease susceptibility and allow for new insights into these deadly forms of prostate cancer.

Metastasis-associated protein ribosomal RNA processing 1 homolog B (RRP1B) modulates metastasis through regulation of histone methylation.

UNLABELLED: Overexpression of ribosomal RNA processing 1 homolog B (RRP1B) induces a transcriptional profile that accurately predicts patient outcome in breast cancer. However, the mechanism by which RRP1B modulates transcription is unclear. Here, the chromatin-binding properties of RRP1B were examined to define how it regulates metastasis-associated transcription. To identify genome-wide RRP1B-binding sites, high-throughput ChIP-seq was performed in the human breast cancer cell line MDA-MB-231 and HeLa cells using antibodies against endogenous RRP1B. Global changes in repressive marks such as histone H3 lysine 9 trimethylation (H3K9me3) were also examined by ChIP-seq. Analysis of these samples identified 339 binding regions in MDA-MB-231 cells and 689 RRP1B-binding regions in HeLa cells. Among these, 136 regions were common to both cell lines. Gene expression analyses of these RRP1B-binding regions revealed that transcriptional repression is the primary result of RRP1B binding to chromatin. ChIP-reChIP assays demonstrated that RRP1B co-occupies loci with decreased gene expression with the heterochromatin-associated proteins, tripartite motif-containing protein 28 (TRIM28/KAP1), and heterochromatin protein 1-α (CBX5/HP1α). RRP1B occupancy at these loci was also associated with higher H3K9me3 levels, indicative of heterochromatinization mediated by the TRIM28/HP1α complex. In addition, RRP1B upregulation, which is associated with metastasis suppression, induced global changes in histone methylation. IMPLICATIONS: RRP1B, a breast cancer metastasis suppressor, regulates gene expression through heterochromatinization and transcriptional repression, which helps our understanding of mechanisms that drive prognostic gene expression in human breast cancer.

Parathyroid hormone activation of matrix metalloproteinase-13 transcription requires the histone acetyltransferase activity of p300 and PCAF and p300-dependent acetylation of PCAF.

Parathyroid hormone (PTH) regulates the transcription of many genes involved in bone remodeling in osteoblasts. One of these genes is matrix metalloproteinase-13 (MMP-13), which is involved in bone remodeling and early stages of endochondral bone formation. We have previously shown that Mmp-13 gene expression is highly induced by PTH treatment in osteoblastic UMR 106-01 cells, as well as primary osteoblasts. Here, we show that p300/CBP-associated factor (PCAF), in addition to p300 and Runx2, is required for PTH activation of Mmp-13 transcription. PCAF was increasingly recruited to the MMP-13 proximal promoter region after PTH treatment, and this was associated with an increase in RNA polymerase II recruitment and histone acetylation. In addition, PTH treatment increased the acetylation of PCAF, a process that required p300. Knockdown of PCAF, p300, or Runx2 by siRNA decreased Mmp-13 mRNA expression after PTH treatment in both UMR 106-01 cells and primary osteoblasts. We found that there is a mutual dependence between p300 and PCAF to be recruited to the Mmp-13 promoter after PTH treatment. In promoter-reporter assays, p300 and PCAF had an additive effect on PTH stimulation of MMP-13 promoter activity, and this required their histone acetyltransferase activity. Our findings demonstrate that PCAF acts downstream of PTH signaling as a transcriptional coactivator that is required for PTH stimulation of MMP-13 transcription. PCAF cooperates with p300 and Runx2 to mediate PTH activation of MMP-13 transcription.

Runx2 recruits p300 to mediate parathyroid hormone's effects on histone acetylation and transcriptional activation of the matrix metalloproteinase-13 gene.

PTH regulates transcription of a number of genes involved in bone remodeling and calcium homeostasis. We have previously shown that the matrix metalloproteinase-13 (MMP-13) gene is induced by PTH in osteoblastic cells as a secondary response through the protein kinase A pathway requiring the runt domain and activator protein 1 binding sites of the proximal promoter. Here, we investigated the changes PTH causes in histone acetylation in this region (which contains the only deoxyribonuclease-hypersensitive sites in the promoter) leading to MMP-13 gene activation in these cells. Chromatin immunoprecipitation experiments revealed that PTH rapidly increased histone H4 acetylation followed by histone H3 acetylation associated with the different regions of the MMP-13 proximal promoter. The hormone also stimulated p300 histone acetyl transferase activity and increased p300 bound to the MMP-13 proximal promoter, and this required protein synthesis. Upon PTH treatment, Runx2, already bound to the runt domain site of the MMP-13 promoter, interacted with p300, which then acetylated histones H4 and H3. The knockdown of either Runx2 or p300 by RNA interference reduced PTH-induced acetylation of histones H3 and H4, association of p300 with the MMP-13 promoter, and resultant MMP-13 gene transcription. Overall, our studies suggest that without altering the gross chromatin structure, PTH stimulates acetylation of histones H3 and H4 via recruitment of p300 to Runx2 bound to the MMP-13 promoter, resulting in gene activation. This work establishes the molecular basis of transcriptional regulation in osteoblasts by PTH, a hormone acting through a G-protein coupled receptor.

Parathyroid hormone signaling in bone and kidney.

PURPOSE OF REVIEW: Parathyroid hormone (PTH) maintains a physiological balance of calcium and phosphate concentrations by binding to its receptor on the plasma membrane of cells in bone and kidney. It signals through multiple pathways, including protein kinase A and protein kinase C, although a preference for certain pathways is apparent in each organ and function. Here, we will review the recent advancements regarding PTH signaling in bone and kidney. RECENT FINDINGS: Wnt proteins have been reported as important regulators of bone metabolism in both PTH-dependent and independent pathways. Recent studies emphasize its role as a mediator of PTH signaling, as PTH treatment increased the expression of wnt4 and sfrp4 and decreased the expression of Wnt inhibitors such as Sost and sclerostin, leading to an increase in Wnt signaling. In kidney, sodium-hydrogen exchanger regulatory factor 1, originally known for its role in the retention of NaPi-IIa at the apical membrane, was shown to have multiple roles in PTH signaling, both as a mediator and regulator. SUMMARY: PTH activates a number of different signaling pathways by binding to a single receptor in bone and kidney. Recent studies demonstrate the involvement of novel factors as well as additional roles for previously identified downstream factors of PTH.

Identification and characterization of Runx2 phosphorylation sites involved in matrix metalloproteinase-13 promoter activation.

Matrix metalloproteinase-13 (MMP-13) plays a critical role in parathyroid hormone (PTH)-induced bone resorption. PTH acts via protein kinase A (PKA) to phosphorylate and stimulate the transactivation of Runx2 for MMP-13 promoter activation. We show here that PTH stimulated Runx2 phosphorylation in rat osteoblastic cells. Runx2 was phosphorylated on serine 28 and threonine 340 after 8-bromo cyclic adenosine mono phosphate (8-Br-cAMP) treatment. We further demonstrate that in the presence of 8-Br-cAMP, the wild-type Runx2 construct stimulated MMP-13 promoter activity, while the Runx2 construct having mutations at three phosphorylation sites (S28, S347 and T340) was unable to stimulate MMP-13 promoter activity. Thus, we have identified the Runx2 phosphorylation sites necessary for PKA stimulated MMP-13 promoter activation and this event may be critical for bone remodeling.