RNASeqBrowser: a genome browser for simultaneous visualization of raw strand specific RNAseq reads and UCSC genome browser custom tracks

BACKGROUND

Strand specific RNAseq data is now more common in RNAseq projects. Visualizing RNAseq data has become an important matter in Analysis of sequencing data. The most widely used visualization tool is the UCSC genome browser that introduced the custom track concept that enabled researchers to simultaneously visualize gene expression at a particular locus from multiple experiments. Our objective of the software tool is to provide friendly interface for visualization of RNAseq datasets.

RESULTS

This paper introduces a visualization tool (RNASeqBrowser) that incorporates and extends the functionality of the UCSC genome browser. For example, RNASeqBrowser simultaneously displays read coverage, SNPs, InDels and raw read tracks with other BED and wiggle tracks -- all being dynamically built from the BAM file. Paired reads are also connected in the browser to enable easier identification of novel exon/intron borders and chimaeric transcripts. Strand specific RNAseq data is also supported by RNASeqBrowser that displays reads above (positive strand transcript) or below (negative strand transcripts) a central line. Finally, RNASeqBrowser was designed for ease of use for users with few bioinformatic skills, and incorporates the features of many genome browsers into one platform.

CONCLUSIONS

The features of RNASeqBrowser: (1) RNASeqBrowser integrates UCSC genome browser and NGS visualization tools such as IGV. It extends the functionality of the UCSC genome browser by adding several new types of tracks to show NGS data such as individual raw reads, SNPs and InDels. (2) RNASeqBrowser can dynamically generate RNA secondary structure. It is useful for identifying non-coding RNA such as miRNA. (3) Overlaying NGS wiggle data is helpful in displaying differential expression and is simple to implement in RNASeqBrowser. (4) NGS data accumulates a lot of raw reads. Thus, RNASeqBrowser collapses exact duplicate reads to reduce visualization space. Normal PC's can show many windows of NGS individual raw reads without much delay. (5) Multiple popup windows of individual raw reads provide users with more viewing space. This avoids existing approaches (such as IGV) which squeeze all raw reads into one window. This will be helpful for visualizing multiple datasets simultaneously. RNASeqBrowser and its manual are freely available at http://www.australianprostatecentre.org/research/software/rnaseqbrowser or http://sourceforge.net/projects/rnaseqbrowser/.

Targeting ASCT2-mediated glutamine uptake blocks prostate cancer growth and tumour development

Glutamine is conditionally essential in cancer cells, being utilized as a carbon and nitrogen source for macromolecule production, as well as for anaplerotic reactions fuelling the tricarboxylic acid (TCA) cycle. In this study, we demonstrated that the glutamine transporter ASCT2 (SLC1A5) is highly expressed in prostate cancer patient samples. Using LNCaP and PC-3 prostate cancer cell lines, we showed that chemical or shRNA-mediated inhibition of ASCT2 function in vitro decreases glutamine uptake, cell cycle progression through E2F transcription factors, mTORC1 pathway activation and cell growth. Chemical inhibition also reduces basal oxygen consumption and fatty acid synthesis, showing that downstream metabolic function is reliant on ASCT2-mediated glutamine uptake. Furthermore, shRNA knockdown of ASCT2 in PC-3 cell xenografts significantly inhibits tumour growth and metastasis in vivo, associated with the down-regulation of E2F cell cycle pathway proteins. In conclusion, ASCT2-mediated glutamine uptake is essential for multiple pathways regulating the cell cycle and cell growth, and is therefore a putative therapeutic target in prostate cancer.

Denhaminols A-H, dihydro-β-agarofurans from the endemic Australian rainforest plant Denhamia celastroides

Eight new dihydro-β-agarofurans, denhaminols A-H (1-8), were isolated from the leaves of the Australian rainforest tree Denhamia celastroides. The chemical structures of 1-8 were elucidated following analysis of 1D/2D NMR and MS data. The absolute configuration of denhaminol A (1) was determined by single-crystal X-ray crystallography. All compounds were evaluated for cytotoxic activity against the human prostate cancer cell line LNCaP, using live-cell imaging and metabolic assays. Denhaminols A (1) and G (7) were also tested for their effects on the lipid content of LNCaP cells. This is the first report of secondary metabolites from D. celastroides.

J-Circos: an interactive Circos plotter

SUMMARY

Circos plots are graphical outputs that display three dimensional chromosomal interactions and fusion transcripts. However, the Circos plot tool is not an interactive visualization tool, but rather a figure generator. For example, it does not enable data to be added dynamically nor does it provide information for specific data points interactively. Recently, an R-based Circos tool (RCircos) has been developed to integrate Circos to R, but similarly, Rcircos can only be used to generate plots. Thus, we have developed a Circos plot tool (J-Circos) that is an interactive visualization tool that can plot Circos figures, as well as being able to dynamically add data to the figure, and providing information for specific data points using mouse hover display and zoom in/out functions. J-Circos uses the Java computer language to enable, it to be used on most operating systems (Windows, MacOS, Linux). Users can input data into J-Circos using flat data formats, as well as from the Graphical user interface (GUI). J-Circos will enable biologists to better study more complex chromosomal interactions and fusion transcripts that are otherwise difficult to visualize from next-generation sequencing data.

AVAILABILITY AND IMPLEMENTATION

J-circos and its manual are freely available at http://www.australianprostatecentre.org/research/software/jcircos

CONTACT

j.an@qut.edu.au

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Androgen-targeted therapy-induced epithelial mesenchymal plasticity and neuroendocrine transdifferentiation in prostate cancer: an opportunity for intervention

Androgens regulate biological pathways to promote proliferation, differentiation, and survival of benign and malignant prostate tissue. Androgen receptor (AR) targeted therapies exploit this dependence and are used in advanced prostate cancer to control disease progression. Contemporary treatment regimens involve sequential use of inhibitors of androgen synthesis or AR function. Although targeting the androgen axis has clear therapeutic benefit, its effectiveness is temporary, as prostate tumor cells adapt to survive and grow. The removal of androgens (androgen deprivation) has been shown to activate both epithelial-to-mesenchymal transition (EMT) and neuroendocrine transdifferentiation (NEtD) programs. EMT has established roles in promoting biological phenotypes associated with tumor progression (migration/invasion, tumor cell survival, cancer stem cell-like properties, resistance to radiation and chemotherapy) in multiple human cancer types. NEtD in prostate cancer is associated with resistance to therapy, visceral metastasis, and aggressive disease. Thus, activation of these programs via inhibition of the androgen axis provides a mechanism by which tumor cells can adapt to promote disease recurrence and progression. Brachyury, Axl, MEK, and Aurora kinase A are molecular drivers of these programs, and inhibitors are currently in clinical trials to determine therapeutic applications. Understanding tumor cell plasticity will be important in further defining the rational use of androgen-targeted therapies clinically and provides an opportunity for intervention to prolong survival of men with metastatic prostate cancer.

Differential effects of tissue culture coating substrates on prostate cancer cell adherence, morphology and behavior

Weak cell-surface adhesion of cell lines to tissue culture surfaces is a common problem and presents technical limitations to the design of experiments. To overcome this problem, various surface coating protocols have been developed. However, a comparative and precise real-time measurement of their impact on cell behavior has not been conducted. The prostate cancer cell line LNCaP, derived from a patient lymph node metastasis, is a commonly used model system in prostate cancer research. However, the cells’ characteristically weak attachment to the surface of tissue culture vessels and cover slips has impeded their manipulation and analysis and use in high throughput screening. To improve the adherence of LNCaP cells to the culture surface, we compared different coating reagents (poly-l-lysine, poly-l-ornithine, collagen type IV, fibronectin, and laminin) and culturing conditions and analyzed their impact on cell proliferation, adhesion, morphology, mobility and gene expression using real-time technologies. The results showed that fibronectin, poly-l-lysine and poly-l-ornithine improved LNCaP cells adherence and provoked cell morphology alterations, such as increase of nuclear and cellular area. These coating reagents also induced a higher expression of F-actin and reduced cell mobility. In contrast, laminin and collagen type IV did not improve adherence but promoted cell aggregation and affected cell morphology. Cells cultured in the presence of laminin displayed higher mobility than control cells. All the coating conditions significantly affected cell viability; however, they did not affect the expression of androgen receptor-regulated genes. Our comparative findings provide important insight for the selection of the ideal coating reagent and culture conditions for the cancer cell lines with respect to their effect on proliferation rate, attachment, morphology, migration, transcriptional response and cellular cytoskeleton arrangement.

Identification of eusynstyelamide B as a potent cell cycle inhibitor following the generation and screening of an ascidian-derived extract library using a real time cell analyzer

Ascidians are marine invertebrates that have been a source of numerous cytotoxic compounds. Of the first six marine-derived drugs that made anticancer clinical trials, three originated from ascidian specimens. In order to identify new anti-neoplastic compounds, an ascidian extract library (143 samples) was generated and screened in MDA-MB-231 breast cancer cells using a real-time cell analyzer (RTCA). This resulted in 143 time-dependent cell response profiles (TCRP), which are read-outs of changes to the growth rate, morphology, and adhesive characteristics of the cell culture. Twenty-one extracts affected the TCRP of MDA-MB-231 cells and were further investigated regarding toxicity and specificity, as well as their effects on cell morphology and cell cycle. The results of these studies were used to prioritize extracts for bioassay-guided fractionation, which led to the isolation of the previously identified marine natural product, eusynstyelamide B (1). This bis-indole alkaloid was shown to display an IC₅₀ of 5 µM in MDA-MB-231 cells. Moreover, 1 caused a strong cell cycle arrest in G2/M and induced apoptosis after 72 h treatment, making this molecule an attractive candidate for further mechanism of action studies.

miRPlant: an integrated tool for identification of plant miRNA from RNA sequencing data

Background

Small RNA sequencing is commonly used to identify novel miRNAs and to determine their expression levels in plants. There are several miRNA identification tools for animals such as miRDeep, miRDeep2 and miRDeep*. miRDeep-P was developed to identify plant miRNA using miRDeep’s probabilistic model of miRNA biogenesis, but it depends on several third party tools and lacks a user-friendly interface. The objective of our miRPlant program is to predict novel plant miRNA, while providing a user-friendly interface with improved accuracy of prediction.

Result

We have developed a user-friendly plant miRNA prediction tool called miRPlant. We show using 16 plant miRNA datasets from four different plant species that miRPlant has at least a 10% improvement in accuracy compared to miRDeep-P, which is the most popular plant miRNA prediction tool. Furthermore, miRPlant uses a Graphical User Interface for data input and output, and identified miRNA are shown with all RNAseq reads in a hairpin diagram.

Conclusions

We have developed miRPlant which extends miRDeep* to various plant species by adopting suitable strategies to identify hairpin excision regions and hairpin structure filtering for plants. miRPlant does not require any third party tools such as mapping or RNA secondary structure prediction tools. miRPlant is also the first plant miRNA prediction tool that dynamically plots miRNA hairpin structure with small reads for identified novel miRNAs. This feature will enable biologists to visualize novel pre-miRNA structure and the location of small RNA reads relative to the hairpin. Moreover, miRPlant can be easily used by biologists with limited bioinformatics skills.

miRPlant and its manual are freely available at http://www.australianprostatecentre.org/research/software/mirplant or http://sourceforge.net/projects/mirplant/.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-275) contains supplementary material, which is available to authorized users.

Paracrine interactions between LNCaP prostate cancer cells and bioengineered bone in 3D in vitro culture reflect molecular changes during bone metastasis

As microenvironmental factors such as three-dimensionality and cell-matrix interactions are increasingly being acknowledged by cancer biologists, more complex 3D in vitro models are being developed to study tumorigenesis and cancer progression. To better understand the pathophysiology of bone metastasis, we have established and validated a 3D indirect co-culture model to investigate the paracrine interactions between prostate cancer (PCa) cells and human osteoblasts. Co-culture of the human PCa, LNCaP cells embedded within polyethylene glycol hydrogels with human osteoblasts in the form of a tissue engineered bone construct (TEB), resulted in reduced proliferation of LNCaP cells. LNCaP cells in both monoculture and co-culture were responsive to the androgen analog, R1881, as indicated by an increase in the expression (mRNA and/or protein induction) of androgen-regulated genes including prostate specific antigen and fatty acid synthase. Microarray gene expression analysis further revealed an up-regulation of bone markers and other genes associated with skeletal and vasculature development and a significant activation of transforming growth factor β1 downstream genes in LNCaP cells after co-culture with TEB. LNCaP cells co-cultured with TEB also unexpectedly showed similar changes in classical androgen-responsive genes under androgen-deprived conditions not seen in LNCaP monocultures. The molecular changes of LNCaP cells after co-culturing with TEBs suggest that osteoblasts exert a paracrine effect that may promote osteomimicry and modulate the expression of androgen-responsive genes in LNCaP cells. Taken together, we have presented a novel 3D in vitro model that allows the study of cellular and molecular changes occurring in PCa cells and osteoblasts that are relevant to metastatic colonization of bone. This unique in vitro model could also facilitate cancer biologists to dissect specific biological hypotheses via extensive genomic or proteomic assessments to further our understanding of the PCa-bone crosstalk.

PITX2 and non-canonical Wnt pathway interaction in metastatic prostate cancer

The non-canonical Wnt pathway, a regulator of cellular motility and morphology, is increasingly implicated in cancer metastasis. In a quantitative PCR array analysis of 84 Wnt pathway associated genes, both non-canonical and canonical pathways were activated in primary and metastatic tumors relative to normal prostate. Expression of the Wnt target gene PITX2 in a prostate cancer (PCa) bone metastasis was strikingly elevated over normal prostate (over 2,000-fold) and primary prostate cancer (over 200-fold). The elevation of PITX2 protein was also evident on tissue microarrays, with strong PITX2 immunostaining in PCa skeletal and, to a lesser degree, soft tissue metastases. PITX2 is associated with cell migration during normal tissue morphogenesis. In our studies, overexpression of individual PITX2A/B/C isoforms stimulated PC-3 PCa cell motility, with the PITX2A isoform imparting a specific motility advantage in the presence of non-canonical Wnt5a stimulation. Furthermore, PITX2 specific shRNA inhibited PC-3 cell migration toward bone cell derived chemoattractant. These experimental results support a pivotal role of PITX2A and non-canonical Wnt signaling in enhancement of PCa cell motility, suggest PITX2 involvement in homing of PCa to the skeleton, and are consistent with a role for PITX2 in PCa metastasis to soft and bone tissues. Our findings, which significantly expand previous evidence that PITX2 is associated with risk of PCa biochemical recurrence, indicate that variation in PITX2 expression accompanies and may promote prostate tumor progression and metastasis.

Exosomes in prostate cancer: putting together the pieces of a puzzle

Exosomes have been shown to act as mediators for cell to cell communication and as a potential source of biomarkers for many diseases, including prostate cancer. Exosomes are nanosized vesicles secreted by cells and consist of proteins normally found in multivesicular bodies, RNA, DNA and lipids. As a potential source of biomarkers, exosomes have attracted considerable attention, as their protein content resembles that of their cells of origin, even though it is noted that the proteins, miRNAs and lipids found in the exosomes are not a reflective stoichiometric sampling of the contents from the parent cells. While the biogenesis of exosomes in dendritic cells and platelets has been extensively characterized, much less is known about the biogenesis of exosomes in cancer cells. An understanding of the processes involved in prostate cancer will help to further elucidate the role of exosomes and other extracellular vesicles in prostate cancer progression and metastasis. There are few methodologies available for general isolation of exosomes, however validation of those methodologies is necessary to study the role of exosomal-derived biomarkers in various diseases. In this review, we discuss “exosomes” as a member of the family of extracellular vesicles and their potential to provide candidate biomarkers for prostate cancer.

Targeting amino acid transport in metastatic castration-resistant prostate cancer: effects on cell cycle, cell growth, and tumor development

BACKGROUND

L-type amino acid transporters (LATs) uptake neutral amino acids including L-leucine into cells, stimulating mammalian target of rapamycin complex 1 signaling and protein synthesis. LAT1 and LAT3 are overexpressed at different stages of prostate cancer, and they are responsible for increasing nutrients and stimulating cell growth.

METHODS

We examined LAT3 protein expression in human prostate cancer tissue microarrays. LAT function was inhibited using a leucine analog (BCH) in androgen-dependent and -independent environments, with gene expression analyzed by microarray. A PC-3 xenograft mouse model was used to study the effects of inhibiting LAT1 and LAT3 expression. Results were analyzed with the Mann-Whitney U or Fisher exact tests. All statistical tests were two-sided.

RESULTS

LAT3 protein was expressed at all stages of prostate cancer, with a statistically significant decrease in expression after 4-7 months of neoadjuvant hormone therapy (4-7 month mean = 1.571; 95% confidence interval = 1.155 to 1.987 vs 0 month = 2.098; 95% confidence interval = 1.962 to 2.235; P = .0187). Inhibition of LAT function led to activating transcription factor 4-mediated upregulation of amino acid transporters including ASCT1, ASCT2, and 4F2hc, all of which were also regulated via the androgen receptor. LAT inhibition suppressed M-phase cell cycle genes regulated by E2F family transcription factors including critical castration-resistant prostate cancer regulatory genes UBE2C, CDC20, and CDK1. In silico analysis of BCH-downregulated genes showed that 90.9% are statistically significantly upregulated in metastatic castration-resistant prostate cancer. Finally, LAT1 or LAT3 knockdown in xenografts inhibited tumor growth, cell cycle progression, and spontaneous metastasis in vivo.

CONCLUSION

Inhibition of LAT transporters may provide a novel therapeutic target in metastatic castration-resistant prostate cancer, via suppression of mammalian target of rapamycin complex 1 activity and M-phase cell cycle genes.

Ribosomal frameshifting used in influenza A virus expression occurs within the sequence UCC_UUU_CGU and is in the +1 direction

Programmed ribosomal frameshifting is used in the expression of many virus genes and some cellular genes. In eukaryotic systems, the most well-characterized mechanism involves -1 tandem tRNA slippage on an X_XXY_YYZ motif. By contrast, the mechanisms involved in programmed +1 (or -2) slippage are more varied and often poorly characterized. Recently, a novel gene, PA-X, was discovered in influenza A virus and found to be expressed via a shift to the +1 reading frame. Here, we identify, by mass spectrometric analysis, both the site (UCC_UUU_CGU) and direction (+1) of the frameshifting that is involved in PA-X expression. Related sites are identified in other virus genes that have previously been proposed to be expressed via +1 frameshifting. As these viruses infect insects (chronic bee paralysis virus), plants (fijiviruses and amalgamaviruses) and vertebrates (influenza A virus), such motifs may form a new class of +1 frameshift-inducing sequences that are active in diverse eukaryotes.

IGF2 increases de novo steroidogenesis in prostate cancer cells

IGF2 is a mitogenic foetal growth factor commonly over-expressed in cancers, including prostate cancer (PC). We recently demonstrated that insulin can activate de novo steroidogenesis in PC cells, a major pathway for reactivation of androgen pathways and PC progression. IGF2 can activate the IGF1 receptor (IGF1R) or insulin receptor (INSR) or hybrids of these two receptors. We therefore hypothesized that IGF2 may contribute to PC progression via de novo steroidogenesis. IGF2 mRNA but not IGF2 receptor mRNA expression was increased in patient samples during progression to castrate-resistant PC as was immunoreactivity to INSR and IGF1R antibodies. Treatment of androgen receptor (AR)-positive PC cell lines LNCaP and 22RV1 with IGF2 for 48 h resulted in increased expression of steroidogenic enzyme mRNA and protein, including steroid acute regulatory protein (StAR), cytochrome p450 family member (CYP)17A1, aldo-keto reductase family member (AKR)1C3 and hydroxysteroid dehydrogenase (HSD)17B3. IGF2 treatment resulted in increased steady state steroid levels and increased de novo steroidogenesis resulting in AR activation as demonstrated by PSA mRNA induction. Inhibition of the IGF1R/INSR signalling axis attenuated the effects of IGF2 on steroid hormone synthesis. We present a potential mechanism for prostatic IGF2 contributing to PC progression by inducing steroidogenesis and that IGF2 signalling and related pathways present attractive targets for PC therapy.

Cloning of a novel insulin-regulated ghrelin transcript in prostate cancer

Ghrelin is a multifunctional hormone, with roles in stimulating appetite and regulating energy balance, insulin secretion and glucose homoeostasis. The ghrelin gene locus (GHRL) is highly complex and gives rise to a range of novel transcripts derived from alternative first exons and internally spliced exons. The wild-type transcript encodes a 117 amino acid preprohormone that is processed to yield the 28 amino acid peptide ghrelin. Here, we identified insulin-responsive transcription corresponding to cryptic exons in intron 2 of the human ghrelin gene. A transcript, termed in2c-ghrelin (intron 2-cryptic), was cloned from the testis and the LNCaP prostate cancer cell line. This transcript may encode an 83 amino acid preproghrelin isoform that codes for ghrelin, but not obestatin. It is expressed in a limited number of normal tissues and in tumours of the prostate, testis, breast and ovary. Finally, we confirmed that in2c-ghrelin transcript expression, as well as the recently described in1-ghrelin transcript, is significantly upregulated by insulin in cultured prostate cancer cells. Metabolic syndrome and hyperinsulinaemia have been associated with prostate cancer risk and progression. This may be particularly significant after androgen deprivation therapy for prostate cancer, which induces hyperinsulinaemia, and this could contribute to castrate-resistant prostate cancer growth. We have previously demonstrated that ghrelin stimulates prostate cancer cell line proliferation in vitro. This study is the first description of insulin regulation of a ghrelin transcript in cancer and should provide further impetus for studies into the expression, regulation and function of ghrelin gene products.

mCOPA: analysis of heterogeneous features in cancer expression data

Background

Cancer outlier profile analysis (COPA) has proven to be an effective approach to analyzing cancer expression data, leading to the discovery of the TMPRSS2 and ETS family gene fusion events in prostate cancer. However, the original COPA algorithm did not identify down-regulated outliers, and the currently available R package implementing the method is similarly restricted to the analysis of over-expressed outliers. Here we present a modified outlier detection method, mCOPA, which contains refinements to the outlier-detection algorithm, identifies both over- and under-expressed outliers, is freely available, and can be applied to any expression dataset.

Results

We compare our method to other feature-selection approaches, and demonstrate that mCOPA frequently selects more-informative features than do differential expression or variance-based feature selection approaches, and is able to recover observed clinical subtypes more consistently. We demonstrate the application of mCOPA to prostate cancer expression data, and explore the use of outliers in clustering, pathway analysis, and the identification of tumour suppressors. We analyse the under-expressed outliers to identify known and novel prostate cancer tumour suppressor genes, validating these against data in Oncomine and the Cancer Gene Index. We also demonstrate how a combination of outlier analysis and pathway analysis can identify molecular mechanisms disrupted in individual tumours.

Conclusions

We demonstrate that mCOPA offers advantages, compared to differential expression or variance, in selecting outlier features, and that the features so selected are better able to assign samples to clinically annotated subtypes. Further, we show that the biology explored by outlier analysis differs from that uncovered in differential expression or variance analysis. mCOPA is an important new tool for the exploration of cancer datasets and the discovery of new cancer subtypes, and can be combined with pathway and functional analysis approaches to discover mechanisms underpinning heterogeneity in cancers.

miRDeep*: an integrated application tool for miRNA identification from RNA sequencing data

miRDeep and its varieties are widely used to quantify known and novel micro RNA (miRNA) from small RNA sequencing (RNAseq). This article describes miRDeep*, our integrated miRNA identification tool, which is modeled off miRDeep, but the precision of detecting novel miRNAs is improved by introducing new strategies to identify precursor miRNAs. miRDeep* has a user-friendly graphic interface and accepts raw data in FastQ and Sequence Alignment Map (SAM) or the binary equivalent (BAM) format. Known and novel miRNA expression levels, as measured by the number of reads, are displayed in an interface, which shows each RNAseq read relative to the pre-miRNA hairpin. The secondary pre-miRNA structure and read locations for each predicted miRNA are shown and kept in a separate figure file. Moreover, the target genes of known and novel miRNAs are predicted using the TargetScan algorithm, and the targets are ranked according to the confidence score. miRDeep* is an integrated standalone application where sequence alignment, pre-miRNA secondary structure calculation and graphical display are purely Java coded. This application tool can be executed using a normal personal computer with 1.5 GB of memory. Further, we show that miRDeep* outperformed existing miRNA prediction tools using our LNCaP and other small RNAseq datasets. miRDeep* is freely available online at http://www.australianprostatecentre.org/research/software/mirdeep-star.

Phenotypic characterization of prostate cancer LNCaP cells cultured within a bioengineered microenvironment

Biophysical and biochemical properties of the microenvironment regulate cellular responses such as growth, differentiation, morphogenesis and migration in normal and cancer cells. Since two-dimensional (2D) cultures lack the essential characteristics of the native cellular microenvironment, three-dimensional (3D) cultures have been developed to better mimic the natural extracellular matrix. To date, 3D culture systems have relied mostly on collagen and Matrigel™ hydrogels, allowing only limited control over matrix stiffness, proteolytic degradability, and ligand density. In contrast, bioengineered hydrogels allow us to independently tune and systematically investigate the influence of these parameters on cell growth and differentiation. In this study, polyethylene glycol (PEG) hydrogels, functionalized with the Arginine-glycine-aspartic acid (RGD) motifs, common cell-binding motifs in extracellular matrix proteins, and matrix metalloproteinase (MMP) cleavage sites, were characterized regarding their stiffness, diffusive properties, and ability to support growth of androgen-dependent LNCaP prostate cancer cells. We found that the mechanical properties modulated the growth kinetics of LNCaP cells in the PEG hydrogel. At culture periods of 28 days, LNCaP cells underwent morphogenic changes, forming tumor-like structures in 3D culture, with hypoxic and apoptotic cores. We further compared protein and gene expression levels between 3D and 2D cultures upon stimulation with the synthetic androgen R1881. Interestingly, the kinetics of R1881 stimulated androgen receptor (AR) nuclear translocation differed between 2D and 3D cultures when observed by immunofluorescent staining. Furthermore, microarray studies revealed that changes in expression levels of androgen responsive genes upon R1881 treatment differed greatly between 2D and 3D cultures. Taken together, culturing LNCaP cells in the tunable PEG hydrogels reveals differences in the cellular responses to androgen stimulation between the 2D and 3D environments. Therefore, we suggest that the presented 3D culture system represents a powerful tool for high throughput prostate cancer drug testing that recapitulates tumor microenvironment.

Seed-based biclustering of gene expression data

BACKGROUND

Accumulated biological research outcomes show that biological functions do not depend on individual genes, but on complex gene networks. Microarray data are widely used to cluster genes according to their expression levels across experimental conditions. However, functionally related genes generally do not show coherent expression across all conditions since any given cellular process is active only under a subset of conditions. Biclustering finds gene clusters that have similar expression levels across a subset of conditions. This paper proposes a seed-based algorithm that identifies coherent genes in an exhaustive, but efficient manner.

METHODS

In order to find the biclusters in a gene expression dataset, we exhaustively select combinations of genes and conditions as seeds to create candidate bicluster tables. The tables have two columns (a) a gene set, and (b) the conditions on which the gene set have dissimilar expression levels to the seed. First, the genes with less than the maximum number of dissimilar conditions are identified and a table of these genes is created. Second, the rows that have the same dissimilar conditions are grouped together. Third, the table is sorted in ascending order based on the number of dissimilar conditions. Finally, beginning with the first row of the table, a test is run repeatedly to determine whether the cardinality of the gene set in the row is greater than the minimum threshold number of genes in a bicluster. If so, a bicluster is outputted and the corresponding row is removed from the table. Repeating this process, all biclusters in the table are systematically identified until the table becomes empty.

CONCLUSIONS

This paper presents a novel biclustering algorithm for the identification of additive biclusters. Since it involves exhaustively testing combinations of genes and conditions, the additive biclusters can be found more readily.

Understanding requirements of novel healthcare information systems for management of advanced prostate cancer

Effective management of chronic diseases is a global health priority. A healthcare information system offers opportunities to address challenges of chronic disease management. However, the requirements of health information systems are often not well understood. The accuracy of requirements has a direct impact on the successful design and implementation of a health information system. Our research describes methods used to understand the requirements of health information systems for advanced prostate cancer management. The research conducted a survey to identify heterogeneous sources of clinical records. Our research showed that the General Practitioner was the common source of patient's clinical records (41%) followed by the Urologist (14%) and other clinicians (14%). Our research describes a method to identify diverse data sources and proposes a novel patient journey browser prototype that integrates disparate data sources.

Androgen receptor and nutrient signaling pathways coordinate the demand for increased amino acid transport during prostate cancer progression

L-Type amino acid transporters such as LAT1 and LAT3 mediate the uptake of essential amino acids. Here, we report that prostate cancer cells coordinate the expression of LAT1 and LAT3 to maintain sufficient levels of leucine needed for mTORC1 signaling and cell growth. Inhibiting LAT function was sufficient to decrease cell growth and mTORC1 signaling in prostate cancer cells. These cells maintained levels of amino acid influx through androgen receptor-mediated regulation of LAT3 expression and ATF4 regulation of LAT1 expression after amino acid deprivation. These responses remained intact in primary prostate cancer, as indicated by high levels of LAT3 in primary disease, and by increased levels of LAT1 after hormone ablation and in metastatic lesions. Taken together, our results show how prostate cancer cells respond to demands for increased essential amino acids by coordinately activating amino acid transporter pathways vital for tumor outgrowth.

Insulin increases de novo steroidogenesis in prostate cancer cells

Androgen-dependent pathways regulate maintenance and growth of normal and malignant prostate tissues. Androgen deprivation therapy (ADT) exploits this dependence and is used to treat metastatic prostate cancer; however, regression initially seen with ADT gives way to development of incurable castration-resistant prostate cancer (CRPC). Although ADT generates a therapeutic response, it is also associated with a pattern of metabolic alterations consistent with metabolic syndrome including elevated circulating insulin. Because CRPC cells are capable of synthesizing androgens de novo, we hypothesized that insulin may also influence steroidogenesis in CRPC. In this study, we examined this hypothesis by evaluating the effect of insulin on steroid synthesis in prostate cancer cell lines. Treatment with 10 nmol/L insulin increased mRNA and protein expression of steroidogenesis enzymes and upregulated the insulin receptor substrate insulin receptor substrate 2 (IRS-2). Similarly, insulin treatment upregulated intracellular testosterone levels and secreted androgens, with the concentrations of steroids observed similar to the levels reported in prostate cancer patients. With similar potency to dihydrotestosterone, insulin treatment resulted in increased mRNA expression of prostate-specific antigen. CRPC progression also correlated with increased expression of IRS-2 and insulin receptor in vivo. Taken together, our findings support the hypothesis that the elevated insulin levels associated with therapeutic castration may exacerbate progression of prostate cancer to incurable CRPC in part by enhancing steroidogenesis.

Interactions between human osteoblasts and prostate cancer cells in a novel 3D in vitro model

Cell-cell and cell-matrix interactions play a major role in tumor morphogenesis and cancer metastasis. Therefore, it is crucial to create a model with a biomimetic microenvironment that allows such interactions to fully represent the pathophysiology of a disease for an in vitro study. This is achievable by using three-dimensional (3D) models instead of conventional two-dimensional (2D) cultures with the aid of tissue engineering technology. We are now able to better address the complex intercellular interactions underlying prostate cancer (CaP) bone metastasis through such models. In this study, we assessed the interaction of CaP cells and human osteoblasts (hOBs) within a tissue engineered bone (TEB) construct. Consistent with other in vivo studies, our findings show that intercellular and CaP cell-bone matrix interactions lead to elevated levels of matrix metalloproteinases, steroidogenic enzymes and the CaP biomarker, prostate specific antigen (PSA); all associated with CaP metastasis. Hence, it highlights the physiological relevance of this model. We believe that this model will provide new insights for understanding of the previously poorly understood molecular mechanisms of bone metastasis, which will foster further translational studies, and ultimately offer a potential tool for drug screening.