MADDD-seq, a novel massively parallel sequencing tool for simultaneous detection of DNA damage and mutations.

Our genome is exposed to a wide variety of DNA-damaging agents. If left unrepaired, this damage can be converted into mutations that promote carcinogenesis or the development of genetically inherited diseases. As a result, researchers and clinicians require tools that can detect DNA damage and mutations with exceptional sensitivity. In this study, we describe a massively parallel sequencing tool termed Mutation And DNA Damage Detection-seq (MADDD-seq) that is capable of detecting O6-methyl guanine lesions and mutations simultaneously, with a single assay. To illustrate the dual capabilities of MADDD-seq, we treated WT and DNA repair deficient yeast cells with the DNA-damaging agent MNNG and tracked DNA lesions and mutations over a 24-h time period. This approach allowed us to identify thousands of DNA adducts and mutations in a single sequencing run and gain deep insight into the kinetics of DNA repair and mutagenesis.

Astrocytic laminin-211 drives disseminated breast tumor cell dormancy in brain.

Although dormancy is thought to play a key role in the metastasis of breast tumor cells to the brain, our knowledge of the molecular mechanisms regulating disseminated tumor cell (DTC) dormancy in this organ is limited. Here using serial intravital imaging of dormant and metastatic triple-negative breast cancer lines, we identify escape from the single-cell or micrometastatic state as the rate-limiting step towards brain metastasis. We show that every DTC occupies a vascular niche, with quiescent DTCs residing on astrocyte endfeet. At these sites, astrocyte-deposited laminin-211 drives DTC quiescence by inducing the dystroglycan receptor to associate with yes-associated protein, thereby sequestering it from the nucleus and preventing its prometastatic functions. These findings identify a brain-specific mechanism of DTC dormancy and highlight the need for a more thorough understanding of tumor dormancy to develop therapeutic approaches that prevent brain metastasis.

Identifying Abundant Immunotherapy and Other Targets in Solid Tumors: Integrating RNA-seq and Mass Spectrometry Proteomics Data Sets.

RNA-seq and mass-spectrometry proteomics combined with growing data repositories have greatly increased the capacity to identify candidate proteins or protein sequence variants that share properties of ideal therapy targets, which include being abundant in cancer cells, absent or rare in adult organs (especially vital organs), and shared by many patient tumors. RNA-seq and fixed content arrays can identify genes that are overexpressed or misexpressed in cancer. RNA-seq is uniquely suited to identifying cancer-specific sequence variants. We review factors relevant for determining whether products of genes that are abundant or differentially abundant in RNA-seq are concordant or discordant with proteins that are identified as abundant or differentially abundant in mass-spectrometry proteomics assays.

The androgen-regulated protease TMPRSS2 activates a proteolytic cascade involving components of the tumor microenvironment and promotes prostate cancer metastasis.

UNLABELLED: TMPRSS2 is an androgen-regulated cell-surface serine protease expressed predominantly in prostate epithelium. TMPRSS2 is expressed highly in localized high-grade prostate cancers and in the majority of human prostate cancer metastases. Through the generation of mouse models with a targeted deletion of Tmprss2, we demonstrate that the activity of this protease regulates cancer cell invasion and metastasis to distant organs. By screening combinatorial peptide libraries, we identified a spectrum of TMPRSS2 substrates that include pro-hepatocyte growth factor (HGF). HGF activated by TMPRSS2 promoted c-MET receptor tyrosine kinase signaling, and initiated a proinvasive epithelial-to-mesenchymal transition phenotype. Chemical library screens identified a potent bioavailable TMPRSS2 inhibitor that suppressed prostate cancer metastasis in vivo. Together, these findings provide a mechanistic link between androgen-regulated signaling programs and prostate cancer metastasis that operate via context-dependent interactions with extracellular constituents of the tumor microenvironment. SIGNIFICANCE: The vast majority of prostate cancer deaths are due to metastasis. Loss of TMPRSS2 activity dramatically attenuated the metastatic phenotype through mechanisms involving the HGF-c-MET axis. Therapeutic approaches directed toward inhibiting TMPRSS2 may reduce the incidence or progression of metastasis in patients with prostate cancer.

Exome sequencing identifies a spectrum of mutation frequencies in advanced and lethal prostate cancers.

To catalog protein-altering mutations that may drive the development of prostate cancers and their progression to metastatic disease systematically, we performed whole-exome sequencing of 23 prostate cancers derived from 16 different lethal metastatic tumors and three high-grade primary carcinomas. All tumors were propagated in mice as xenografts, designated the LuCaP series, to model phenotypic variation, such as responses to cancer-directed therapeutics. Although corresponding normal tissue was not available for most tumors, we were able to take advantage of increasingly deep catalogs of human genetic variation to remove most germline variants. On average, each tumor genome contained ~200 novel nonsynonymous variants, of which the vast majority was specific to individual carcinomas. A subset of genes was recurrently altered across tumors derived from different individuals, including TP53, DLK2, GPC6, and SDF4. Unexpectedly, three prostate cancer genomes exhibited substantially higher mutation frequencies, with 2,000-4,000 novel coding variants per exome. A comparison of castration-resistant and castration-sensitive pairs of tumor lines derived from the same prostate cancer highlights mutations in the Wnt pathway as potentially contributing to the development of castration resistance. Collectively, our results indicate that point mutations arising in coding regions of advanced prostate cancers are common but, with notable exceptions, very few genes are mutated in a substantial fraction of tumors. We also report a previously undescribed subtype of prostate cancers exhibiting "hypermutated" genomes, with potential implications for resistance to cancer therapeutics. Our results also suggest that increasingly deep catalogs of human germline variation may challenge the necessity of sequencing matched tumor-normal pairs.

The retinoblastoma tumor suppressor controls androgen signaling and human prostate cancer progression.

Retinoblastoma (RB; encoded by RB1) is a tumor suppressor that is frequently disrupted in tumorigenesis and acts in multiple cell types to suppress cell cycle progression. The role of RB in tumor progression, however, is poorly defined. Here, we have identified a critical role for RB in protecting against tumor progression through regulation of targets distinct from cell cycle control. In analyses of human prostate cancer samples, RB loss was infrequently observed in primary disease and was predominantly associated with transition to the incurable, castration-resistant state. Further analyses revealed that loss of the RB1 locus may be a major mechanism of RB disruption and that loss of RB function was associated with poor clinical outcome. Modeling of RB dysfunction in vitro and in vivo revealed that RB controlled nuclear receptor networks critical for tumor progression and that it did so via E2F transcription factor 1-mediated regulation of androgen receptor (AR) expression and output. Through this pathway, RB depletion induced unchecked AR activity that underpinned therapeutic bypass and tumor progression. In agreement with these findings, disruption of the RB/E2F/nuclear receptor axis was frequently observed in the transition to therapy resistance in human disease. Together, these data reveal what we believe to be a new paradigm for RB function in controlling prostate tumor progression and lethal tumor phenotypes.

Combining results of microarray experiments: a rank aggregation approach.

As technology for microarray analysis becomes widespread, it is becoming increasingly important to be able to compare and combine the results of experiments that explore the same scientific question. In this article, we present a rank-aggregation approach for combining results from several microarray studies. The motivation for this approach is twofold; first, the final results of microarray studies are typically expressed as lists of genes, rank-ordered by a measure of the strength of evidence that they are functionally involved in the disease process, and second, using the information on this rank-ordered metric means that we do not have to concern ourselves with data on the actual expression levels, which may not be comparable across experiments. Our approach draws on methods for combining top-k lists from the computer science literature on meta-search. The meta-search problem shares several important features with that of combining microarray experiments, including the fact that there are typically few lists with many elements and the elements may not be common to all lists. We implement two meta-search algorithms, which use a Markov chain framework to convert pairwise preferences between list elements into a stationary distribution that represents an aggregate ranking (Dwork et al, 2001). We explore the behavior of the algorithms in hypothetical examples and a simulated dataset and compare their performance with that of an algorithm based on the order-statistics model of Thurstone (Thurstone, 1927). We apply all three algorithms to aggregate the results of five microarray studies of prostate cancer.

Substrates of the prostate-specific serine protease prostase/KLK4 defined by positional-scanning peptide libraries.

BACKGROUND: Prostase/KLK4 is a member of the human kallikrein (KLK) gene family that is expressed in prostate epithelial cells under the regulation of androgenic hormones. In this study, we sought to characterize the substrate specificity of KLK4 in order to gain insight into potential physiological roles of the enzyme. METHODS: A chimeric form of KLK4 was constructed in which the pro-region of KLK4 was replaced with the signal and propeptide sequence of trypsinogen (proT-KLK4) to create an activation site susceptible to enterokinase cleavage. proT-KLK4 was expressed in Drosophila S2 cells, purified, and activated with enterokinase to generate mature KLK4. The extended substrate specificity of KLK4 was defined by screening tetrapeptide positional scanning synthetic combinatorial libraries (PS-SCL). RESULTS: The preferred P1-P4 positions as determined by PS-SCL were: P1-Arg; P2-Gln/Leu/Val; P3-Gln/Ser/Val; P4-Ile/Val. The trypsin-like specificity of KLK4 was further confirmed using synthetic chromogenic peptides. Based upon the optimal cleavage site residues, a database search for potential KLK4 substrates identified several proteins with potential roles mediating normal prostate physiology or neoplastic growth including KLK3/PSA, parathyroid hormone-related peptide (PTHrP), and members of the bone morphogenetic protein (BMP) family. Recombinant KLK4 was able to activate pro-PSA/KLK3 and degrade members of the insulin-like growth factor (IGF) binding protein (IGFBP) family. CONCLUSIONS: These results identify potential KLK4 substrates that may serve to define the role of this protease in normal prostate physiology, and facilitate studies of the consequences of KLK4 expression in pathological conditions.

Characterization and comparative analyses of transcriptomes from the normal and neoplastic human prostate.

BACKGROUND: The prostate gland is a highly specialized organ with functional attributes that serve to enhance the fertility of mammalian species. Pathological processes affecting the prostate include benign prostate hypertrophy and prostate carcinoma; diseases that account for major morbidity and mortality in middle-aged and elderly men. To facilitate studies of biological processes uniquely represented in the prostate and assess molecular alterations associated with prostate carcinoma, we sought to establish the diversity of gene expression in the normal and neoplastic prostate through the compilation and analysis of a prostate transcriptome. METHODS: We assembled and annotated ESTs derived from prostate cDNA libraries that were either produced in our laboratory or available from public sequence repositories such as CGAP, dbEST, and Unigene. Determinations of differential gene expression between the normal prostate, other normal tissues, and neoplastic prostate tissues was performed using statistical algorithms. Confirmation of differential expression was performed by quantitative PCR and Northern analysis. RESULTS: A total of 99,448 high-quality ESTs were assembled and annotated to produce a prostate transcriptome comprised of 24,580 distinct TUs. Comparative analyses of gene expression levels identified 61 TUs with exclusive expression in the prostate and 45 TUs with high levels of expression in the prostate relative to at least 25 other normal tissues (P > 0.99). Comparative analyses of ESTs derived from neoplastic prostate tissues identified 75 genes with dysregulated expression in cancer (P > 0.99). CONCLUSIONS: The human prostate expresses a diverse repertoire of genes that reflect a functionally complex organ. The identification of genes with prostate-restricted or enhanced expression may provide additional insights into the biochemical processes that interact to form the developmental, signaling, and metabolic pathways of the normal and neoplastic gland.

Expressed sequence tag profiling identifies developmental and anatomic partitioning of gene expression in the mouse prostate.

BACKGROUND: The prostate gland is an organ with highly specialized functional attributes that serves to enhance the fertility of mammalian species. Much of the information pertaining to normal and pathological conditions affecting the prostate has been obtained through extensive developmental, biochemical and genetic analyses of rodent species. Although important insights can be obtained through detailed anatomical and histological assessments of mouse and rat models, further mechanistic explanations are greatly aided through studies of gene and protein expression. RESULTS: In this article we characterize the repertoire of genes expressed in the normal developing mouse prostate through the analysis of 50,562 expressed sequence tags derived from 14 mouse prostate cDNA libraries. Sequence assemblies and annotations identified 15,009 unique transcriptional units of which more than 600 represent high quality assemblies without corresponding annotations in public gene expression databases. Quantitative analyses demonstrate distinct anatomical and developmental partitioning of prostate gene expression. This finding may assist in the interpretation of comparative studies between human and mouse and guide the development of new transgenic murine disease models. The identification of several novel genes is reported, including a new member of the beta-defensin gene family with prostate-restricted expression. CONCLUSIONS: These findings suggest a potential role for the prostate as a defensive barrier for entry of pathogens into the genitourinary tract and, further, serve to emphasize the utility of the continued evaluation of transcriptomes from a diverse repertoire of tissues and cell types.

Molecular characterization of prostatic small-cell neuroendocrine carcinoma.

OBJECTIVES: A subset of prostate carcinomas is composed predominantly, even exclusively, of neuroendocrine (NE) cells. In this report, we sought to characterize the gene expression profile of a prostate small cell NE carcinoma by assessing the diversity and abundance of transcripts in the LuCaP 49 prostate small cell carcinoma xenograft. METHODS: We constructed a cDNA library (PRCA3) from the LuCap 49 prostate small cell xenograft. Single pass DNA sequencing of randomly selected cDNA clones followed by sequence assembly and annotation produced a library of Expressed Sequence Tags (ESTs) representing the LuCaP 49 transcriptome. Comparative sequence analysis with ESTs derived from prostate adenocarcinoma libraries was performed using statistical algorithms designed to identify differentially expressed sequences. Putative NE cell-specific genes were further examined by Northern analysis. RESULTS: Sequence assembly and analysis identified 1,447 distinct genes expressed in the LuCaP 49 cDNA library. These include cDNAs encoding the NE markers secretogranin (SCG2), CD24, and ENO2. Northern analysis revealed that three additional genes, ASCL1, INA, and SV2B are expressed in LuCaP 49 but not in various prostate cancer cell lines or xenografts. Fifteen genes were identified with a statistical probability (P > 0.9) of being up-regulated in LuCaP 49 small cell carcinoma relative to prostate adenocarcinoma (two primary prostate adenocarcinomas and the LNCaP prostate adenocarcinoma cell line). CONCLUSIONS: Prostate small cell carcinoma expresses a diverse repertoire of genes that reflect characteristics of their NE cell of origin. ASCL1, INA, and SV2B are potential molecular markers for small cell NE tumors and NE cells of the prostate. This small cell NE carcinoma gene expression profile may yield insights into the development, progression, and treatment of subtypes of prostate cancer.

Maelstrom, a Drosophila spindle-class gene, encodes a protein that colocalizes with Vasa and RDE1/AGO1 homolog, Aubergine, in nuage.

A hallmark of germline cells across the animal kingdom is the presence of perinuclear, electron-dense granules called nuage. In many species examined, Vasa, a DEAD-box RNA helicase, is found in these morphologically distinct particles. Despite its evolutionary conservation, the function of nuage remains obscure. We have characterized a null allele of maelstrom (mael) and shown that Maelstrom protein is localized to nuage in a Vasa-dependent manner. By phenotypic characterization, we have defined maelstrom as a spindle-class gene that affects Vasa modification. In a nuclear transport assay, we have determined that Maelstrom shuttles between the nucleus and cytoplasm, which may indicate a nuclear origin for nuage components. Interestingly, Maelstrom, but not Vasa, depends on two genes involved in RNAi phenomena, aubergine and spindle-E (spn-E), for its nuage localization. Furthermore, maelstrom mutant ovaries show mislocalization of two proteins involved in the microRNA and/or RNAi pathways, Dicer and Argonaute2, suggesting a potential connection between nuage and the microRNA-pathway.

A neuroendocrine/small cell prostate carcinoma xenograft-LuCaP 49.

The late stages of progression of prostate carcinoma are typically characterized by an androgen-insensitive, rapidly proliferative state. Some late-stage tumors are composed predominantly of neuroendocrine cells. Virtually no animal models of a neuroendocrine/small cell variant of prostate carcinoma are available for experimental studies. We report a human neuroendocrine/small cell prostate carcinoma xenograft that was developed from a nodal metastasis of a human prostate carcinoma and that has been propagated as serial subcutaneous implants in severe combined immunodeficient mice for >4 years. Designated LuCaP 49, all tumor passages exhibit a neuroendocrine/small cell carcinoma phenotype-insensitivity to androgen deprivation, expression of neuroendocrine proteins, lack of expression of prostate-specific antigen or androgen receptor, and an unusually rapid growth (a doubling time of 6.5 days) for prostate cancer xenografts. Genetically this tumor exhibits loss of heterozygosity for the short arm of chromosome 8 and has a complex karyotype. This xenograft should prove to be useful in the investigation of mechanisms underlying the androgen-insensitive state of progressive prostate carcinoma.

The program of androgen-responsive genes in neoplastic prostate epithelium.

The human prostate gland is an important target organ of androgenic hormones. Testosterone and dihydrotestosterone interact with the androgen receptor to regulate vital aspects of prostate growth and function including cellular proliferation, differentiation, apoptosis, metabolism, and secretory activity. Our objective in this study was to characterize the temporal program of transcription that reflects the cellular response to androgens and to identify specific androgen-regulated genes (ARGs) or gene networks that participate in these responses. We used cDNA microarrays representing about 20,000 distinct human genes to profile androgen-responsive transcripts in the LNCaP adenocarcinoma cell line and identified 146 genes with transcript alterations more than 3-fold. Of these, 103 encode proteins with described functional roles, and 43 represent transcripts that have yet to be characterized. Temporal gene expression profiles grouped the ARGs into four distinct cohorts. Five uncharacterized ARGs demonstrated exclusive or high expression levels in the prostate relative to other tissues studied. A search of available DNA sequence upstream of 28 ARGs identified 25 with homology to the androgen response-element consensus-binding motif. These results identify previously uncharacterized and unsuspected genes whose expression levels are directly or indirectly regulated by androgens; further, they provide a comprehensive temporal view of the transcriptional program of human androgen-responsive cells.

Digital expression profiles of the prostate androgen-response program.

The androgen receptor (AR) and cognate ligands regulate vital aspects of prostate cellular growth and function including proliferation, differentiation, apoptosis, lipid metabolism, and secretory action. In addition, the AR pathway also influences pathological processes of the prostate such as benign prostatic hypertrophy and prostate carcinogenesis. The pivotal role of androgens and the AR in prostate biology prompted this study with the objective of identifying molecular mediators of androgen action. Our approach was designed to compare transcriptomes of the LNCaP prostate cancer cell line under conditions of androgen depletion and androgen stimulation by generating and comparing collections of expressed sequence tags (ESTs). A total of 4400 ESTs were produced from LNCaP cDNA libraries and these ESTs assembled into 2486 distinct transcripts. Rigorous statistical analysis of the expression profiles indicated that 17 genes exhibited a high probability (P>0.9) of androgen-regulated expression. Northern analysis confirmed that the expression of KLK3/PSA, FKBP5, KRT18, DKFZP564K247, DDX15, and HSP90 is regulated by androgen exposure. Of these, only KLK3/PSA is known to be androgen-regulated while the other genes represent new members of the androgen-response program in prostate epithelium. LNCaP gene expression profiles defined by two independent experiments using the serial analysis of gene expression (SAGE) method were compared with the EST profiles. Distinctly different expression patterns were produced from each dataset. These results are indicative of the sensitivity of the methods to experimental conditions and demonstrate the power and the statistical limitations of digital expression analyses.

The human (PEDB) and mouse (mPEDB) Prostate Expression Databases.

The Prostate Expression Databases (PEDB and mPEDB) are online resources designed to allow researchers to access and analyze gene expression information derived from the human and murine prostate, respectively. Human PEDB archives more than 84 000 Expressed Sequence Tags (ESTs) from 38 prostate cDNA libraries in a curated relational database that provides detailed library information including tissue source, library construction methods, sequence diversity and sequence abundance. The differential expression of each EST species can be viewed across all libraries using a Virtual Expression Analysis Tool (VEAT), a graphical user interface written in Java for intra- and inter-library sequence comparisons. Recent enhancements to PEDB include (i) the development of a murine prostate expression database, mPEDB, that complements the human gene expression information in PEDB, (ii) the assembly of a non-redundant sequence set or 'prostate unigene' that represents the diversity of gene expression in the prostate, and (iii) an expanded search tool that supports both text-based and BLAST queries. PEDB and mPEDB are accessible via the World Wide Web at http://www.pedb.org and http://www.mpedb.org.