Empowering the biomedical research community: Innovative SAS deployment on the All of Us Researcher Workbench.

OBJECTIVES: The All of Us Research Program is a precision medicine initiative aimed at establishing a vast, diverse biomedical database accessible through a cloud-based data analysis platform, the Researcher Workbench (RW). Our goal was to empower the research community by co-designing the implementation of SAS in the RW alongside researchers to enable broader use of All of Us data. MATERIALS AND METHODS: Researchers from various fields and with different SAS experience levels participated in co-designing the SAS implementation through user experience interviews. RESULTS: Feedback and lessons learned from user testing informed the final design of the SAS application. DISCUSSION: The co-design approach is critical for reducing technical barriers, broadening All of Us data use, and enhancing the user experience for data analysis on the RW. CONCLUSION: Our co-design approach successfully tailored the implementation of the SAS application to researchers' needs. This approach may inform future software implementations on the RW.

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.