Search

Prostate cancer is initially responsive to androgen ablation therapy and progresses to androgen-unresponsive states that are refractory to treatment. The mechanism of this transition is unknown. A systems approach to disease begins with the quantitative delineation of the informational elements (mRNAs and proteins) in various disease states. We employed two recently developed high-throughput technologies, massively parallel signature sequencing (MPSS) and isotope-coded affinity tag, to gain a comprehensive picture of the changes in mRNA levels and more restricted analysis of protein levels, respectively, during the transition from androgen-dependent LNCaP (model for early-stage prostate cancer) to androgen-independent CL1 cells (model for late-stage prostate cancer). We sequenced >5 million MPSS signatures, obtained >142,000 tandem mass spectra, and built comprehensive MPSS and proteomic databases. The integrated mRNA and protein expression data revealed underlying functional differences between androgen-dependent and androgen-independent prostate cancer cells. The high sensitivity of MPSS enabled us to identify virtually all of the expressed transcripts and to quantify the changes in gene expression between these two cell states, including functionally important low-abundance mRNAs, such as those encoding transcription factors and signal transduction molecules. These data enable us to map the differences onto extant physiologic networks, creating perturbation networks that reflect prostate cancer progression. We found 37 BioCarta and 14 Kyoto Encyclopedia of Genes and Genomes pathways that are up-regulated and 23 BioCarta and 22 Kyoto Encyclopedia of Genes and Genomes pathways that are down-regulated in LNCaP cells versus CL1 cells. Our efforts represent a significant step toward a systems approach to understanding prostate cancer progression.

Supplementary Figure 1 from Evidence for the Presence of Disease-Perturbed Networks in Prostate Cancer Cells by Genomic and Proteomic Analyses: A Systems Approach to Disease

Supplementary Table S3 from Evidence for the Presence of Disease-Perturbed Networks in Prostate Cancer Cells by Genomic and Proteomic Analyses: A Systems Approach to Disease

Supplementary Table S7 from Evidence for the Presence of Disease-Perturbed Networks in Prostate Cancer Cells by Genomic and Proteomic Analyses: A Systems Approach to Disease

Supplementary Table S4 from Evidence for the Presence of Disease-Perturbed Networks in Prostate Cancer Cells by Genomic and Proteomic Analyses: A Systems Approach to Disease

Supplementary Table S2 from Evidence for the Presence of Disease-Perturbed Networks in Prostate Cancer Cells by Genomic and Proteomic Analyses: A Systems Approach to Disease

Supplementary Table S5 from Evidence for the Presence of Disease-Perturbed Networks in Prostate Cancer Cells by Genomic and Proteomic Analyses: A Systems Approach to Disease

Supplementary Table S1 from Evidence for the Presence of Disease-Perturbed Networks in Prostate Cancer Cells by Genomic and Proteomic Analyses: A Systems Approach to Disease

Supplementary Table S6 from Evidence for the Presence of Disease-Perturbed Networks in Prostate Cancer Cells by Genomic and Proteomic Analyses: A Systems Approach to Disease

Cancer diagnosis can be a complex process, which takes consideration of histopathological, clinical, immunophenotypic, and genetic features. Since non-random chromosomal translocations are specifically involved in the development of various cancers, the detection of these gene aberrations becomes increasingly important. In recent years, break-apart (or split-signal) fluorescence in situ hybridization (FISH) has emerged as an advantageous technique to detect gene translocations on tissue sections. However, FISH assays are technically challenging and require specialized fluorescence microscopes. Furthermore, the FISH signal is not stable for long term archiving due to photo bleaching. Our objective was to demonstrate the feasibility of brightfield break-apart in situ hybridization (ba-ISH) for anaplastic lymphoma kinase (ALK) and mucosa-associated lymphoid tissue translocation protein 1 (MALT1) genes as models. ALK or MALT1 break-apart probes were labeled with digoxigenin (DIG) or 2,4-dinitrophenyl (DNP) on both sides of a known gene breakpoint region and the hybridization sites were visualized with the combination of alkaline phosphatase (AP)-based blue and red detection. Therefore, normal genes are detected as purple dots by mixing blue and red colors while translocated genes are detected as isolated blue or red dots. Formalin-fixed, paraffin-embedded tonsil was used as control for the co-localized 5' and 3' probes. Gene translocations of ALK or MALT1 were detected as separate blue and red dots on ALCL and MALT lymphoma cases. Thus, ISH analyses of gene translocations can be conducted with a regular light microscope and the long term archiving of break-apart ISH slides can be achieved.

BACKGROUND Fifteen human tissue kallikrein (KLK) genes have been identified as a cluster on chromosome 19. KLK expression is associated with various human diseases including cancers. Noncoding RNAs such as PCA3/DD3 and PCGEM1 have been identified in prostate cancer cells. METHODS Using massively parallel signature sequencing (MPSS) technology, RT-PCR, and 5′ rapid amplification of cDNA ends (RACE), we identified and cloned a novel gene that maps to the KLK locus. RESULTS We have characterized this gene, named as KLK31P by the HUGO Gene Nomenclature Committee, as an unprocessed KLK pseudogene. It contains five exons, two of which are KLK-derived while the rest are “exonized” interspersed repeats. KLK31P is expressed abundantly in prostate tissues and is androgen regulated. KLK31P is expressed at lower levels in localized and metastatic prostate cancer cells than in normal prostate cells. CONCLUSIONS KLK31P is a novel androgen regulated and transcribed pseudogene of kallikreins that may play a role in prostate carcinogenesis or maintenance. Prostate 66: 936–944, 2006. © 2006 Wiley-Liss, Inc.

A major goal in genomics is to understand how genes are regulated in different tissues, stages of development, diseases, and species. Mapping DNase I hypersensitive (HS) sites within nuclear chromatin is a powerful and well-established method of identifying many different types of regulatory elements, but in the past it has been limited to analysis of single loci. We have recently described a protocol to generate a genome-wide library of DNase HS sites. Here, we report high-throughput analysis, using massively parallel signature sequencing (MPSS), of 230,000 tags from a DNase library generated from quiescent human CD4+ T cells. Of the tags that uniquely map to the genome, we identified 14,190 clusters of sequences that group within close proximity to each other. By using a real-time PCR strategy, we determined that the majority of these clusters represent valid DNase HS sites. Approximately 80% of these DNase HS sites uniquely map within one or more annotated regions of the genome believed to contain regulatory elements, including regions 2 kb upstream of genes, CpG islands, and highly conserved sequences. Most DNase HS sites identified in CD4+ T cells are also HS in CD8+ T cells, B cells, hepatocytes, human umbilical vein endothelial cells (HUVECs), and HeLa cells. However, ∼10% of the DNase HS sites are lymphocyte specific, indicating that this procedure can identify gene regulatory elements that control cell type specificity. This strategy, which can be applied to any cell line or tissue, will enable a better understanding of how chromatin structure dictates cell function and fate.

We have used massively parallel signature sequencing (MPSS) to sample the transcriptomes of 32 normal human tissues to an unprecedented depth, thus documenting the patterns of expression of almost 20,000 genes with high sensitivity and specificity. The data confirm the widely held belief that differences in gene expression between cell and tissue types are largely determined by transcripts derived from a limited number of tissue-specific genes, rather than by combinations of more promiscuously expressed genes. Expression of a little more than half of all known human genes seems to account for both the common requirements and the specific functions of the tissues sampled. A classification of tissues based on patterns of gene expression largely reproduces classifications based on anatomical and biochemical properties. The unbiased sampling of the human transcriptome achieved by MPSS supports the idea that most human genes have been mapped, if not functionally characterized. This data set should prove useful for the identification of tissue-specific genes, for the study of global changes induced by pathological conditions, and for the definition of a minimal set of genes necessary for basic cell maintenance. The data are available on the Web at http://mpss.licr.org and http://sgb.lynxgen.com.

The mouse Fused locus encodes a protein that has been implicated in the regulation of embryonic axis formation. The protein, which has been named Axin to distinguish it from the product of the unrelated Drosophila melanogaster gene fused , contains regions of similarity to the RGS (regulators of G-protein signaling) family of proteins as well as to dishevelled, a protein that acts downstream of Wingless in D. melanogaster . Loss-of-function mutations at Fused lead to lethality between days 8 and 10 of gestation. Three dominant mutations result in a kinked tail in heterozygotes. Two of the dominant mutations, Fused and Knobbly , result from insertions of intracisternal A particle retrotransposons into the gene. The insertion in Fused , within the sixth intron, creates a gene that produces wild-type transcripts as well as mutant transcripts that initiate at both the authentic promoter and the 3′-most long terminal repeat of the insertion. Knobbly , an insertion of the retrotransposon into exon 7, precludes the production of wild-type protein. Thus the Fused homozygote is viable whereas Knobbly is a recessive embryonic lethal. In both mutants the dominant kink-tailed phenotype is likely to result from the synthesis of similar amino-terminal fragments of Axin protein that would contain the RGS domain, but lack the dishevelled domain.

Spontaneous mutations at the mouse Fused (Fu) locus cause dominant skeletal and neurological defects and recessive lethal embryonic defects including neuroectodermal abnormalities and axial duplications. Here, we describe a new allele at the Fu locus caused by a transgenic insertional mutation, H epsilon 46. Embryos homozygous for the H epsilon 46 insertion die at day 9-10 post coitum and display phenotypic defects similar to those associated with Fu alleles. The H epsilon 46 locus was cloned and shown to contain a 20-kb deletion at the site of transgene insertion with no other detectable rearrangements. Genomic probes from the H epsilon 46 locus were mapped to a genetic locus closely linked to Fu on chromosome 17 and were hybridized to a YAC contig covering the FuKi critical region. Compound heterozygotes between H epsilon 46 and FuKi were inviable and displayed abnormalities at the same stage of embryogenesis as do homozygotes for either of the two mutations, demonstrating that these two recessive lethal mutations belong to the same complementation group. A genomic probe from the wild-type H epsilon 46 locus detected a transcript that is disrupted by the transgenic insertion, representing a candidate for the wild-type allele of Fused.

We have created several transgenic mouse strains that bear the human lambda light chain gene driven by its own promoter and a mouse immunoglobulin heavy chain enhancer. The transgene is expressed in many tissues, with particularly high levels of expression in the bone marrow, thymus, spleen, and lymph nodes. One of these transgenic lines, B-less, displays a dramatic phenotype characterized by an acute susceptibility to bacterial and viral infections. Analysis of this strain shows it to be profoundly deficient in both immature (pre-B) and mature B cells, as well as in circulating immunoglobulin. The pre-B and B cell defects are cell autonomous, as judged by cell culture and bone marrow graft chimeras. Despite this B cell deficiency, the T cell lineage appears grossly normal as assessed by flow cytometric analysis and by its response to mitogen stimulation. Since an independently derived transgenic strain bearing the same human lambda construct displays a partial B-less phenotype, it is likely that the B lineage deficiency is due to a dominant effect of transgene expression rather than to the insertional perturbation of an endogenous mouse gene. It is interesting that the deficiency phenotype is fully expressed in the FVB/N genetic background, but is suppressed in F1 hybrids formed between the FVB/N and C57BL/6 inbred strains. Evidently, there are one or more dominant genetic suppressors of B-less in the C57BL/6 genome.

Transcription factors play a key role in integrating and modulating biological information. In this study, we comprehensively measured the changing abundances of mRNAs over a time course of activation of human peripheral-blood-derived mononuclear cells (“macrophages”) with lipopolysaccharide. Global and dynamic analysis of transcription factors in response to a physiological stimulus has yet to be achieved in a human system, and our efforts significantly advanced this goal. We used multiple global high-throughput technologies for measuring mRNA levels, including massively parallel signature sequencing and GeneChip microarrays. We identified 92 of 1,288 known human transcription factors as having significantly measurable changes during our 24-h time course. At least 42 of these changes were previously unidentified in this system. Our data demonstrate that some transcription factors operate in a functional range below 10 transcripts per cell, whereas others operate in a range three orders of magnitude greater. The highly reproducible response of many mRNAs indicates feedback control. A broad range of activation kinetics was observed; thus, combinatorial regulation by small subsets of transcription factors would permit almost any timing input to cis -regulatory elements controlling gene transcription.

Massively parallel signature sequencing is an ultra-high throughput sequencing technology. It can simultaneously sequence millions of sequence tags, and, therefore, is ideal for whole genome analysis. When applied to expression profiling, it reveals almost every transcript in the sample and provides its accurate expression level. This chapter describes the technology and its application in establishing stem cell transcriptome databases.

Massively parallel signature sequencing is an ultra-high throughput sequencing technology. It can simultaneously sequence millions of sequence tags, and, therefore, is ideal for whole genome analysis. When applied to expression profiling, it reveals almost every transcript in the sample and provides its accurate expression level. This chapter describes the technology and its application in establishing stem cell transcriptome databases.

HTLV-I is the etiologic agent of adult T-cell leukemia (ATL), a fatal T-cell malignancy that is associated with profound immunosuppression. In this study, comprehensive gene expression profiling was performed using massively parallel signature sequencing (MPSS) to investigate virus-host interactions in acutely HTLV-1 transformed cells. The analysis revealed the modulation of numerous genes across different functional classes, many of which have not been previously implicated in HTLV-1 transformation or ATL. Differences in the transcriptomes of transformed cell lines were observed that have provided clues on how different clonal populations of cells respond to virus transformation. Quantitation of HTLV-1 transcription was possible, thus making MPSS a useful tool to study emerging pathogens and unknown microbial causes of human diseases. Importantly, overexpression of GITR, an activation marker that has not been previously reported to be upregulated by HTLV-1-infection or in transformed/leukemic cells and that is associated with the suppressor phenotype of CD4+CD25+ regulatory T-cells (Tregs), was also observed. The deep and quantitative gene expression profile generated by MPSS should provide additional leads for discovery research that can be applied to better understand the pathobiology of HTLV-1 transformation and ATL as well as to developing new therapies.

Massively Parallel Signature Sequencing (MPSS), a recently developed high-throughput transcription profiling technology, has the ability to profile almost every transcript in a sample without requiring prior knowledge of the sequence of the transcribed genes. As is the case with DNA microarrays, effective data analysis depends crucially on understanding how noise affects measurements. We analyze the sources of noise in MPSS and present a quantitative model describing the variability between replicate MPSS assays. We use this model to construct statistical hypotheses that test whether an observed change in gene expression in a pair-wise comparison is significant. This analysis is then extended to the determination of the significance of changes in expression levels measured over the course of a time series of measurements. We apply these analytic techniques to the study of a time series of MPSS gene expression measurements on LPS-stimulated macrophages. To evaluate our statistical significance metrics, we compare our results with published data on macrophage activation measured by using Affymetrix GeneChips.

Background Pooled human embryonic stem cells (hESC) cell lines were profiled to obtain a comprehensive list of genes common to undifferentiated human embryonic stem cells. Results Pooled hESC lines were profiled to obtain a comprehensive list of genes common to human ES cells. Massively parallel signature sequencing (MPSS) of approximately three million signature tags (signatures) identified close to eleven thousand unique transcripts, of which approximately 25% were uncharacterised or novel genes. Expression of previously identified ES cell markers was confirmed and multiple genes not known to be expressed by ES cells were identified by comparing with public SAGE databases, EST libraries and parallel analysis by microarray and RT-PCR. Chromosomal mapping of expressed genes failed to identify major hotspots and confirmed expression of genes that map to the X and Y chromosome. Comparison with published data sets confirmed the validity of the analysis and the depth and power of MPSS. Conclusions Overall, our analysis provides a molecular signature of genes expressed by undifferentiated ES cells that can be used to monitor the state of ES cells isolated by different laboratories using independent methods and maintained under differing culture conditions

Mutations at the mouse Fused locus have pleiotropic developmental effects, including the formation of axial duplications in homozygous embryos. The product of the Fused locus, Axin, displays similarities to RGS (Regulators of G-Protein Signaling) and Dishevelled proteins. Mutant Fused alleles that cause axial duplications disrupt the major mRNA, suggesting that Axin negatively regulates the response to an axis-inducing signal. Injection of Axin mRNA into Xenopus embryos inhibits dorsal axis formation by interfering with signaling through the Wnt pathway. Furthermore, ventral injection of an Axin mRNA lacking the RGS domain induces an ectopic axis, apparently through a dominant-negative mechanism. Thus, Axin is a novel inhibitor of Wnt signaling and regulates an early step in embryonic axis formation in mammals and amphibians.

We have cloned and characterized a c-myc (now designated MYC) oncogene that had been translocated into the mu switch region of the immunoglobulin heavy chain locus in a Burkitt lymphoma cell line. The breakpoint of the translocation occurs within the first intron of the c-myc gene, thereby separating the untranslocated first exon from the two coding exons. Transcription from the translocated gene arises from a cryptic promoter within the first intron, which produces a 438-nucleotide untranslated 5' region. The amino acid sequence of the protein encoded by the c-myc gene has been substantially altered. In particular, a compensating set of frameshift mutations alters a string of 24 amino acids in a region of the protein tightly conserved in human, mouse, and chicken c-myc genes as well as in the human N-myc and L-myc oncogenes. Despite this, the mutated gene retains a reduced transforming ability in a rat embryo fibroblast focus-formation assay.

From a lambda phage gene library we have isolated phage containing the gene encoding human preproparathyroid hormone. The phage were isolated by using both the plaque-hybridization technique and the in vivo recombination-selection technique. The human preproparathyroid hormone gene contains two intervening sequences that separate the gene into a 5' noncoding domain, a "prepro" sequence domain, and a domain containing the parathyroid hormone sequence and the 3' noncoding region. The gene is approximately 4,200 base pairs long. Restriction endonuclease analysis of human leukocyte DNA shows that the haploid human genome contains one copy of the preproparathyroid hormone gene. A 14-base-pair sequence of alternating purines and pyrimidines that has the potential of adopting the Z-DNA conformation lies 134 base pairs upstream from the presumed site of initiation of transcription.

In contrast to other human tumors in which the c-myc gene and its transcript are greatly amplified, careful analysis of t(8;14) Burkitt cell lines indicates that the c-myc transcript is marginally, and in some cases not at all, increased by comparison to control lymphoblastoid cell lines. Instead, there is a more subtle alteration in the expression of the translocated c-myc gene characterized by a shift in promoter utilization and an apparent insensitivity to the regulation that inactivates the normal c-myc allele within these same cells. In some Burkitt cell lines, such deregulation might be because of the loss of a putative control region through removal of the large dual promoter/leader segment of the c-myc gene. In other cell lines, however, this deregulation may be explained by somatic mutations that occur within the putative control region even though it is located many hundreds of bases from the translocation breakpoint.