Transcriptome analysis of serous ovarian cancers identifies differentially expressed chromosome 3 genes

Bayesian Gene Selection Based on Pathway Information and Network-Constrained Regularization

High-throughput data make it possible to study expression levels of thousands of genes simultaneously under a particular condition. However, only few of the genes are discriminatively expressed. How to identify these biomarkers precisely is significant for disease diagnosis, prognosis, and therapy. Many studies utilized pathway information to identify the biomarkers. However, most of these studies only incorporate the group information while the pathway structural information is ignored. In this paper, we proposed a Bayesian gene selection with a network-constrained regularization method, which can incorporate the pathway structural information as priors to perform gene selection. All the priors are conjugated; thus, the parameters can be estimated effectively through Gibbs sampling. We present the application of our method on 6 microarray datasets, comparing with Bayesian Lasso, Bayesian Elastic Net, and Bayesian Fused Lasso. The results show that our method performs better than other Bayesian methods and pathway structural information can improve the result.

TMEM158 promotes pancreatic cancer aggressiveness by activation of TGFβ1 and PI3K/AKT signaling pathway

Pancreatic cancer (PC) is one of the most deadly digestive cancers world-wide, with a dismal five-year survival rate of <8%. Upregulation of transmembrane protein 158 (TMEM158) is known to facilitate the progression of several carcinomas. However, little is known concerning the potential roles of TMEM158 in PC. Herein, we first found that TMEM158 was significantly upregulated in PC samples as well as PC cell lines. The overexpression of TMEM158 was significantly correlated with advanced clinicopathologic features (including tumor size, TNM stage, and blood vessel invasion) and poorer prognosis of patients with PC in clinic. Evidenced based on a series of loss- and gain-of-function assays uncovered that TMEM158 enhanced PC cell proliferation, migration, and invasion by stimulating the progression of cell cycle, epithelial-mesenchymal transition, and MMP-2/9 production. Furthermore, mechanism-related investigations disclosed that activation of TGFβ1 and PI3K/AKT signal might be responsible for TMEM158-triggered PC aggressiveness. Collectively, TMEM158 was upregulated in PC and promoted PC cell proliferation, migration, and invasion through the activation of TGFβ1 and PI3K/AKT signaling pathways, highlighting its potential as a tumor promoter and a therapeutic target for PC.

Proteome characterization of human pancreatic cyst fluid from intraductal papillary mucinous neoplasm by liquid chromatography/tandem mass spectrometry

RATIONALE

In recent years, the molecular components of pancreatic cyst fluid have been used for diagnosis and prognosis. Because the protein markers that are currently used in clinical tests are unreliable, proteomic studies to find new protein markers are being conducted. However, such researches have been limited due to the complexity of pancreatic cyst fluid and the immaturity of proteomic techniques.

METHODS

To overcome these limitations and provide a pancreatic cyst proteome dataset, we examined cyst fluid proteome with tandem mass spectrometry. The proteomic analysis was performed using a Orbitrap-based mass spectrometer (Q-Exactive) coupled with a 50-cm-long nano-liquid chromatography column. Protein mutations were identified using mutation sequence database search.

RESULTS

A total of 5850 protein groups were identified from microliters of cyst fluid. Among those, 3934 protein groups were reported for the first time in pancreatic cyst fluid. Although high-abundance proteins were not depleted in the experiment, our dataset detected almost all pancreatic tumor markers such as mucin family members, S100 proteins, and CEA-related proteins. In addition, 590 protein mutation marker candidates were discovered.

CONCLUSIONS

We provide a comprehensive cyst proteome dataset that includes cystic cellular proteins and mutated proteins. Our findings would serve as a rich resource for further IPMN studies and clinical applications. The MS data have been deposited in the ProteomeXchange with identifier PXD005671 (http://proteomecentral.proteomexchange.org/dataset/PXD005671).

Low levels of IGFBP7 expression in high-grade serous ovarian carcinoma is associated with patient outcome

BACKGROUND

Insulin-like growth factor binding protein 7 (IGFBP7) has been suggested to act as a tumour suppressor gene in various human cancers, yet its role in epithelial ovarian cancer (EOC) has not yet been investigated. We previously observed that IGFBP7 was one of several genes found significantly upregulated in an EOC cell line model rendered non-tumourigenic as consequence of genetic manipulation. The aim of the present study was to investigate the role of IGFBP7 in high-grade serous ovarian carcinomas (HGSC), the most common type of EOC.

METHODS

We analysed IGFBP7 gene expression in 11 normal ovarian surface epithelial cells (NOSE), 79 high-grade serous ovarian carcinomas (HGSC), and seven EOC cell lines using a custom gene expression array platform. IGFBP7 mRNA expression profiles were also extracted from publicly available databases. Protein expression was assessed by immunohistochemistry of 175 HGSC and 10 normal fallopian tube samples using tissue microarray and related to disease outcome. We used EOC cells to investigate possible mechanisms of gene inactivation and describe various in vitro growth effects of exposing EOC cell lines to human recombinant IGFBP7 protein and conditioned media.

RESULTS

All HGSCs exhibited IGFBP7 expression levels that were significantly (p = 0.001) lower than the mean of the expression value of NOSE samples and that of a whole ovary sample. IGFBP7 gene and protein expression were lower in tumourigenic EOC cell lines relative to a non-tumourigenic EOC cell line. None of the EOC cell lines harboured a somatic mutation in IGFBP7, although loss of heterozygosity (LOH) of the IGFBP7 locus and epigenetic methylation silencing of the IGFBP7 promoter was observed in two of the cell lines exhibiting loss of gene/protein expression. In vitro functional assays revealed an alteration of the EOC cell migration capacity. Protein expression analysis of HGSC samples revealed that the large majority of tumour cores (72.6%) showed low or absence of IGFBP7 staining and revealed a significant correlation between IGFBP7 protein expression and a prolonged overall survival (p = 0.044).

CONCLUSION

The low levels of IGFPB7 in HGSC relative to normal tissues, and association with survival are consistent with a purported role in tumour suppressor pathways.

Screening the Molecular Targets of Ovarian Cancer Based on Bioinformatics Analysis

Aims and background

Ovarian cancer (OC) is the most lethal gynecologic malignancy. This study aims to explore the molecular mechanisms of OC and identify potential molecular targets for OC treatment.

Methods and Study Design

Microarray gene expression data (GSE14407) including 12 normal ovarian surface epithelia samples and 12 OC epithelia samples were downloaded from Gene Expression Omnibus database. Differentially expressed genes (DEGs) between 2 kinds of ovarian tissue were identified by using limma package in R language (|log2 fold change| >1 and false discovery rate [FDR] <0.05). Protein-protein interactions (PPIs) and known OC-related genes were screened from COXPRESdb and GenBank database, respectively. Furthermore, PPI network of top 10 upregulated DEGs and top 10 downregulated DEGs was constructed and visualized through Cytoscape software. Finally, for the genes involved in PPI network, functional enrichment analysis was performed by using DAVID (FDR <0.05).

Results

In total, 1136 DEGs were identified, including 544 downregulated and 592 upregulated DEGs. Then, PPI network was constructed, and DEGs CDKN2A, MUC1, OGN, ZIC1, SOX17, and TFAP2A interacted with known OC-related genes CDK4, EGFR/JUN, SRC, CLI1, CTNNB1, and TP53, respectively. Moreover, functions about oxygen transport and embryonic development were enriched by the genes involved in the network of downregulated DEGs.

Conclusions

We propose that 4 DEGs (OGN, ZIC1, SOX17, and TFAP2A) and 2 functions (oxygen transport and embryonic development) might play a role in the development of OC. These 4 DEGs and known OC-related genes might serve as therapeutic targets for OC. Further studies are required to validate these predictions.

Transcriptional regulation induced by cAMP elevation in mouse Schwann cells

In peripheral nerves, Schwann cell development is regulated by a variety of signals. Some of the aspects of Schwann cell differentiation can be reproduced in vitro in response to forskolin, an adenylyl cyclase activator elevating intracellular cAMP levels. Herein, the effect of forskolin treatment was investigated by a comprehensive genome-wide expression study on primary mouse Schwann cell cultures. Additional to myelin-related genes, many so far unconsidered genes were ascertained to be modulated by forskolin. One of the strongest differentially regulated gene transcripts was the transcription factor Olig1 (oligodendrocyte transcription factor 1), whose mRNA expression levels were reduced in treated Schwann cells. Olig1 protein was localized in myelinating and nonmyelinating Schwann cells within the sciatic nerve as well as in primary Schwann cells, proposing it as a novel transcription factor of the Schwann cell lineage. Data analysis further revealed that a number of differentially expressed genes in forskolin-treated Schwann cells were associated with the ECM (extracellular matrix), underlining its importance during Schwann cell differentiation in vitro. Comparison of samples derived from postnatal sciatic nerves and from both treated and untreated Schwann cell cultures showed considerable differences in gene expression between in vivo and in vitro, allowing us to separate Schwann cell autonomous from tissue-related changes. The whole data set of the cell culture microarray study is provided to offer an interactive search tool for genes of interest.

Identification of eight candidate target genes of the recurrent 3p12-p14 loss in cervical cancer by integrative genomic profiling

The pathogenetic role, including its target genes, of the recurrent 3p12-p14 loss in cervical cancer has remained unclear. To determine the onset of the event during carcinogenesis, we used microarray techniques and found that the loss was the most frequent 3p event, occurring in 61% of 92 invasive carcinomas, in only 2% of 43 high-grade intraepithelial lesions (CIN2/3), and in 33% of 6 CIN3 lesions adjacent to invasive carcinomas, suggesting a role in acquisition of invasiveness or early during the invasive phase. We performed an integrative DNA copy number and expression analysis of 77 invasive carcinomas, where all genes within the recurrent region were included. We selected eight genes, THOC7, PSMD6, SLC25A26, TMF1, RYBP, SHQ1, EBLN2, and GBE1, which were highly down-regulated in cases with loss, as confirmed at the protein level for RYBP and TMF1 by immunohistochemistry. The eight genes were subjected to network analysis based on the expression profiles, revealing interaction partners of proteins encoded by the genes that were coordinately regulated in tumours with loss. Several partners were shared among the eight genes, indicating crosstalk in their signalling. Gene ontology analysis showed enrichment of biological processes such as apoptosis, proliferation, and stress response in the network and suggested a relationship between down-regulation of the eight genes and activation of tumourigenic pathways. Survival analysis showed prognostic impact of the eight-gene signature that was confirmed in a validation cohort of 74 patients and was independent of clinical parameters. These results support the role of the eight candidate genes as targets of the 3p12-p14 loss in cervical cancer and suggest that the strong selection advantage of the loss during carcinogenesis might be caused by a synergetic effect of several tumourigenic processes controlled by these targets.

TMEM158 and FBLP1 as novel marker genes of cisplatin sensitivity in non-small cell lung cancer cells

Even after development of molecular targeting therapies, platinum-based chemotherapy is still a standard care for treatment of locally advanced non-small cell lung cancer (NSCLC). So far, critical molecular markers capable to predict the therapeutic response in NSCLC patients remain undetermined. We here attempted to identify novel biomarker genes for cisplatin (CDDP) for a tailored therapy. Initial screening to explorer association of IC(50) values of CDDP obtained by MTT assay and gene expression levels measured with oligonucleotide microarray and real-time RT-PCR provided 6 candidate genes, namely, NUBPL, C9orf30, ZNF12, TMEM158, GSK3B, and FBLP1 using 9 lung cancer cells consisting of 3 small and 6 NSCLC cells. These 6 genes together with 5 reported biomarkers, i.e., GSTP1, ERCC1, BRCA1, FRAP1, and RRM1, were subjected to a linear regression analysis using 12 NSCLC cell lines including 6 additional NSCLC cells: only FBLP1 and TMEM158 genes showed positive associations with statistical significances (P = .016 and .026, respectively). The biological significance of these genes was explored by in vitro experiments: Knockdown experiments in PC-9/CDDP cells revealed that the reduced expression of TMEM158 significantly decreased the chemo-resistance against CDDP (P <.0001), while 2 transformants of PC-6 cells stably over-expressing FBLP1 resulted in an enhanced resistance to CDDP (P = .004 and P = .001). Furthermore, a stepwise multiple regression analysis demonstrated the best prediction formula could be fixed when we used expression data of TMEM158 and FBLP1 (R(2) = 0.755, P = .0018). TMEM158 and FBLP1 may be powerful predictive biomarkers for CDDP therapy in NSCLC.

FKBP10/FKBP65 expression in high-grade ovarian serous carcinoma and its association with patient outcome

The frequent loss of chromosome 17 in epithelial ovarian carcinomas (EOC), particularly high-grade serous carcinomas (HGSC), has been attributed to the disruption of TP53 (at 17p13.1) and other chromosome 17 genes suspected to play a role in tumour suppressor pathways. In a transcriptome analysis of HGSC, we showed underexpression of a number of chromosome 17 genes, which included FKBP10 (at 17q21.1) and collagen I α 1 (COL1A1; at 17q21.33). FKBP10 codes for the immunophilin FKBP65 and is suspected to act as a chaperone for COL1A1. We have investigated FKBP10 (gene) and FKBP65 (protein) expression in HGSC samples and EOC cell lines that differ in their tumourigenic potential. COL1A1 expression was also investigated given the purported function of FKBP65. RT-PCR analysis verified underexpression of FKBP10 and COL1A1 in HGSCs (n=14) and six tumourigenic EOC cell lines, relative to normal ovarian surface epithelial cells and a non-tumourigenic EOC cell line. Immunohistochemistry analyses of 196 HGSC samples using tissue microarrays revealed variable staining intensities in the epithelial tumour component where only 7.8% and 1.0% of samples stained intensely for FKBP65 and COL1A1, respectively. Variable staining intensities were also observed for the stromal component where 23.6% and 24.1% stained intensely for FKBP65 and COL1A1, respectively. There was no significant correlation of staining intensity of either protein with disease stage. Staining of FKBP65 was clearly visible in normal epithelial cells of the ovarian surface and fallopian tube. There was a significant correlation between absence of FKBP65 staining in the epithelial cell component of the tumour and prolonged overall survival (p<0.001). Our results suggest that underexpression of FKBP65 protein is characteristic of HGSCs and that this expression profile may be linked to molecular pathways associated with an unfavourable outcome in cancer patients.

VGLL3 expression is associated with a tumor suppressor phenotype in epithelial ovarian cancer

Previous studies have implicated vestigial like 3 (VGLL3), a chromosome 3p12.3 gene that encodes a putative transcription co-factor, as a candidate tumor suppressor gene (TSG) in high-grade serous ovarian carcinomas (HGSC), the most common type of epithelial ovarian cancer. A complementation analysis based on microcell-mediated chromosome transfer (MMCT) using a centric fragment of chromosome 3 (der3p12-q12.1) into the OV-90 ovarian cancer cell line haploinsufficient for 3p and lacking VGLL3 expression was performed to assess the effect on tumorigenic potential and growth characteristics. Genetic characterization of the derived MMCT hybrids revealed that only the hybrid that contained an intact VGLL3 locus exhibited alterations of tumorigenic potential in a nude mouse xenograft model and various in vitro growth characteristics. Only stable OV-90 transfectant clones expressing low levels of VGLL3 were derived. These clones exhibited an altered cytoplasmic morphology characterized by numerous single membrane bound multivesicular-bodies (MVB) that were not attributed to autophagy. Overexpression of VGLL3 in OV-90 was achieved using a lentivirus-based tetracycline inducible gene expression system, which also resulted in MVB formation in the infected cell population. Though there was no significant differences in various in vitro and in vivo growth characteristics in a comparison of VGLL3-expressing clones with empty vector transfectant controls, loss of VGLL3 expression was observed in tumors derived from mouse xenograft models. VGLL3 gene and protein expression was significantly reduced in HGSC samples (>98%, p < 0.05) relative to either normal ovarian surface epithelial cells or epithelial cells of the fallopian tube, possible tissues of origin of HGSC. Also, there appeared to be to be more cases with higher staining levels in stromal tissue component from HGSC cases that had a prolonged disease-free survival. The results taken together suggest that VGLL3 is involved in tumor suppressor pathways, a feature that is characterized by the absence of VGLL3 expression in HGSC samples.

Derivation and characterization of matched cell lines from primary and recurrent serous ovarian cancer

Background

Cell line models have proven to be effective tools to investigate a variety of ovarian cancer features. Due to the limited number of cell lines, particularly of the serous subtype, the heterogeneity of the disease, and the lack of cell lines that model disease progression, there is a need to further develop cell line resources available for research. This study describes nine cell lines derived from three ovarian cancer cases that were established at initial diagnosis and at subsequent relapse after chemotherapy.

Methods

The cell lines from three women diagnosed with high-grade serous ovarian cancer (1369, 2295 and 3133) were derived from solid tumor (TOV) and ascites (OV), at specific time points at diagnosis and relapse (R). Primary treatment was a combination of paclitaxel/carboplatin (1369, 3133), or cisplatin/topotecan (2295). Second line treatment included doxorubicin, gemcitabine and topotecan. In addition to molecular characterization (p53, HER2), the cell lines were characterized based on cell growth characteristics including spheroid growth, migration potential, and anchorage independence. The in vivo tumorigenicity potential of the cell lines was measured. Response to paclitaxel and carboplatin was assessed using a clonogenic assay.

Results

All cell lines had either a nonsense or missense TP53 mutations. The ability to form compact spheroids or aggregates was observed in six of nine cell lines. Limited ability for migration and anchorage independence was observed. The OV3133(R) cell line, formed tumors at subcutaneous sites in SCID mice. Based on IC₅₀ values and dose response curves, there was clear evidence of acquired resistance to carboplatin for TOV2295(R) and OV2295(R2) cell lines.

Conclusion

The study identified nine new high-grade serous ovarian cancer cell lines, derived before and after chemotherapy that provides a unique resource for investigating the evolution of this common histopathological subtype of ovarian cancer.

Chromosome 3 anomalies investigated by genome wide SNP analysis of benign, low malignant potential and low grade ovarian serous tumours

Ovarian carcinomas exhibit extensive heterogeneity, and their etiology remains unknown. Histological and genetic evidence has led to the proposal that low grade ovarian serous carcinomas (LGOSC) have a different etiology than high grade carcinomas (HGOSC), arising from serous tumours of low malignant potential (LMP). Common regions of chromosome (chr) 3 loss have been observed in all types of serous ovarian tumours, including benign, suggesting that these regions contain genes important in the development of all ovarian serous carcinomas. A high-density genome-wide genotyping bead array technology, which assayed >600,000 markers, was applied to a panel of serous benign and LMP tumours and a small set of LGOSC, to characterize somatic events associated with the most indolent forms of ovarian disease. The genomic patterns inferred were related to TP53, KRAS and BRAF mutations. An increasing frequency of genomic anomalies was observed with pathology of disease: 3/22 (13.6%) benign cases, 40/53 (75.5%) LMP cases and 10/11 (90.9%) LGOSC cases. Low frequencies of chr3 anomalies occurred in all tumour types. Runs of homozygosity were most commonly observed on chr3, with the 3p12-p11 candidate tumour suppressor region the most frequently homozygous region in the genome. An LMP harboured a homozygous deletion on chr6 which created a GOPC-ROS1 fusion gene, previously reported as oncogenic in other cancer types. Somatic TP53, KRAS and BRAF mutations were not observed in benign tumours. KRAS-mutation positive LMP cases displayed significantly more chromosomal aberrations than BRAF-mutation positive or KRAS and BRAF mutation negative cases. Gain of 12p, which harbours the KRAS gene, was particularly evident. A pathology review reclassified all TP53-mutation positive LGOSC cases, some of which acquired a HGOSC status. Taken together, our results support the view that LGOSC could arise from serous benign and LMP tumours, but does not exclude the possibility that HGOSC may derive from LMP tumours.

Subtype specific elevated expression of hyaluronidase-1 (HYAL-1) in epithelial ovarian cancer

Background

Epithelial ovarian cancer (EOC) is morphologically heterogeneous being classified as serous, endometrioid, clear cell, or mucinous. Molecular genetic analysis has suggested a role for tumor suppressor genes located at chromosome 3p in serous EOC pathogenesis. Our objective was to evaluate the expression of HYAL1, located at chromosome 3p21.3, in these EOC subtypes, and to investigate its correlation with the expression of steroid hormone receptors.

Methodology/Principal Findings

We determined the mRNA expression of HYAL1, estrogen receptor (ER)-α, ERβ and progesterone receptor (PR) in EOC tumor samples and cell lines using quantitative RT-PCR. We also examined the expression of these genes in a publicly available microarray dataset. HYAL-1 enzyme activity was measured in EOC cell lines and in plasma samples from patients. We found that HYAL1 mRNA expression was elevated in clear cell and mucinous EOC tissue samples, but not in serous and endometrioid samples, normal ovaries or benign tumors. Similar results were obtained by two different techniques and with tissue sample cohorts from two independent institutions. Concordantly, HYAL1 mRNA levels and enzymatic activity were elevated only in EOC cell lines derived from clear cell and mucinous subtypes. We also showed that HYAL1 mRNA was inversely correlated to that of ERα specifically in clear cell and mucinous EOCs. Additionally, ectopic expression of ERα in a clear cell EOC cell line (ER- and PR-negative) induced 50% reduction of HYAL1 mRNA expression, supporting a role of ERα in HYAL1 gene regulation. Significantly, HYAL-1 activity was also high in the plasma of patients with these EOC subtypes.

Conclusions/Significance

This is the first report showing high HYAL-1 levels in EOC and demonstrating HYAL1 gene repression by ERα. Our results identify Hyaluronidase-1 as a potential target/biomarker for clear cell and mucinous EOCs and especially in tumors with low ERα levels.

Identification of ChIP-seq mapped targets of HP1β due to bombesin/GRP receptor activation

Epithelial cells lining the adult colon do not normally express gastrin-releasing peptide (GRP) or its receptor (GRPR). In contrast, GRP/GRPR can be aberrantly expressed in human colorectal cancer (CRC) including Caco-2 cells. We have previously shown that GRPR activation results in the up-regulation of HP1β, an epigenetic modifier of gene transcription. The aim of this study was to identify the genes whose expression is altered by HP1β subsequent to GRPR activation. We determined HP1β binding positions throughout the genome using chromatin immunoprecipitation followed by massively parallel DNA sequencing (ChIP-seq). After exposure to GRP, we identified 9,625 genomic positions occupied by HP1β. We performed gene microarray analysis on Caco-2 cells in the absence and presence of a GRPR specific antagonist as well as siRNA to HP1β. The expression of 97 genes was altered subsequent to GRPR antagonism, while the expression of 473 genes was altered by HP1β siRNA exposure. When these data were evaluated in concert with our ChIP-seq findings, 9 genes showed evidence of possible altered expression as a function of GRPR signaling via HP1β. Of these, genomic PCR of immunoprecipitated chromatin demonstrated that GRPR signaling affected the expression of IL1RAPL2, FAM13A, GBE1, PLK3, and SLCO1B3. These findings provide the first evidence by which GRPR aberrantly expressed in CRC might affect tumor progression.

Characterization of the molecular differences between ovarian endometrioid carcinoma and ovarian serous carcinoma

The histopathological diagnosis of high-grade endometrioid and serous carcinoma of the ovary is poorly reproducible under the current morphology based classification system, especially for anaplastic, high-grade tumours. The transcription factor Wilms' tumour-1 (WT1) is differentially expressed among the gynaecological epithelia from which epithelial ovarian cancers (EOCs) are believed to originate. In EOCs, WT1 protein is observed in the majority of serous carcinomas and in up to 30% of endometrioid carcinomas. It is unclear whether the latter is a reflection of the actual incidence of WT1 protein expression in endometrioid carcinomas, or whether a significant number of high-grade serous carcinomas have been misclassified as endometrioid carcinoma. Several genetic aberrations are reported to occur in EOCs. These include mutation of the TP53 gene, aberrant activation of beta-catenin signalling and loss of PTEN protein expression, among others. It is unclear whether these aberrations are histotype-specific. The aim of this study was to better define the molecular characteristics of serous and endometrioid carcinomas in an attempt to address the problems with the current histopathological classification methods. Gene expression profiles were analysed to identify reproducible gene expression phenotypes for endometrioid and serous carcinomas. Tissue microarrays (TMA) were used to assess the incidence of TP53, beta-catenin and PTEN aberrations in order to correlate their occurrence with WT1 as an immunohistochemistry based biomarker of serous histotype. It was found that nuclear WT1 protein expression can identify misclassified high-grade endometrioid carcinomas and these tumours should be reassigned to serous histotype. Although low-grade endometrioid carcinomas rarely progress to high-grade carcinomas, a combined WT1-negative, TP53-positive immunophenotype may identify an uncommon high-grade subtype of ovarian endometrioid carcinoma. GEO database: array data accession number GSE6008.

Characterization of the 3p12.3-pcen region associated with tumor suppression in a novel ovarian cancer cell line model genetically modified by chromosome 3 fragment transfer

The genetic analysis of nontumorigenic radiation hybrids generated by transfer of chromosome 3 fragments into the tumorigenic OV-90 ovarian cancer cell line identified the 3p12.3-pcen region as a candidate tumor suppressor gene (TSG) locus. In the present study, polymorphic microsatellite repeat analysis of the hybrids further defined the 3p12.3-pcen interval to a 16.1 Mb common region containing 12 known or hypothetical genes: 3ptel-ROBO2-ROBO1-GBE1-CADM2-VGLL3-CHMP2B-POU1F1-HTR1F-CGGBP1-ZNF654-C3orf38-EPHA3-3pcen. Seven of these genes, ROBO1, GBE1, VGLL3, CHMP2B, CGGBP1, ZNF654, and C3orf38, exhibited gene expression in the hybrids, placing them as top TSG candidates for further analysis. The expression of all but one (VGLL3) of these genes was also detected in the parental OV-90 cell line. Mutations were not identified in a comparative sequence analysis of the predicted protein coding regions of these candidates in OV-90 and donor normal chromosome 3 contig. However, the nondeleterious sequence variants identified in the transcribed regions distinguished parent of origin alleles for ROBO1, VGLL3, CHMP2B, and CGGBP1 and cDNA sequencing of the hybrids revealed biallelic expression of these genes. Interestingly, underexpression of VGLL3 and ZNF654 were observed in malignant ovarian tumor samples as compared with primary cultures of normal ovarian surface epithelial cells or benign ovarian tumors, and this occurred regardless of allelic content of 3p12.3-pcen. The results taken together suggest that dysregulation of VGLL3 and/or ZNF654 expression may have affected pathways important in ovarian tumorigenesis which was offset by the transfer of chromosome 3 fragments in OV-90, a cell line hemizygous for 3p.

Gene dosage, expression, and ontology analysis identifies driver genes in the carcinogenesis and chemoradioresistance of cervical cancer

Integrative analysis of gene dosage, expression, and ontology (GO) data was performed to discover driver genes in the carcinogenesis and chemoradioresistance of cervical cancers. Gene dosage and expression profiles of 102 locally advanced cervical cancers were generated by microarray techniques. Fifty-two of these patients were also analyzed with the Illumina expression method to confirm the gene expression results. An independent cohort of 41 patients was used for validation of gene expressions associated with clinical outcome. Statistical analysis identified 29 recurrent gains and losses and 3 losses (on 3p, 13q, 21q) associated with poor outcome after chemoradiotherapy. The intratumor heterogeneity, assessed from the gene dosage profiles, was low for these alterations, showing that they had emerged prior to many other alterations and probably were early events in carcinogenesis. Integration of the alterations with gene expression and GO data identified genes that were regulated by the alterations and revealed five biological processes that were significantly overrepresented among the affected genes: apoptosis, metabolism, macromolecule localization, translation, and transcription. Four genes on 3p (RYBP, GBE1) and 13q (FAM48A, MED4) correlated with outcome at both the gene dosage and expression level and were satisfactorily validated in the independent cohort. These integrated analyses yielded 57 candidate drivers of 24 genetic events, including novel loci responsible for chemoradioresistance. Further mapping of the connections among genetic events, drivers, and biological processes suggested that each individual event stimulates specific processes in carcinogenesis through the coordinated control of multiple genes. The present results may provide novel therapeutic opportunities of both early and advanced stage cervical cancers.

The chemiluminescence based Ziplex automated workstation focus array reproduces ovarian cancer Affymetrix GeneChip expression profiles

Background

As gene expression signatures may serve as biomarkers, there is a need to develop technologies based on mRNA expression patterns that are adaptable for translational research. Xceed Molecular has recently developed a Ziplex^® technology, that can assay for gene expression of a discrete number of genes as a focused array. The present study has evaluated the reproducibility of the Ziplex system as applied to ovarian cancer research of genes shown to exhibit distinct expression profiles initially assessed by Affymetrix GeneChip^® analyses.

Methods

The new chemiluminescence-based Ziplex^® gene expression array technology was evaluated for the expression of 93 genes selected based on their Affymetrix GeneChip^® profiles as applied to ovarian cancer research. Probe design was based on the Affymetrix target sequence that favors the 3' UTR of transcripts in order to maximize reproducibility across platforms. Gene expression analysis was performed using the Ziplex Automated Workstation. Statistical analyses were performed to evaluate reproducibility of both the magnitude of expression and differences between normal and tumor samples by correlation analyses, fold change differences and statistical significance testing.

Results

Expressions of 82 of 93 (88.2%) genes were highly correlated (p < 0.01) in a comparison of the two platforms. Overall, 75 of 93 (80.6%) genes exhibited consistent results in normal versus tumor tissue comparisons for both platforms (p < 0.001). The fold change differences were concordant for 87 of 93 (94%) genes, where there was agreement between the platforms regarding statistical significance for 71 (76%) of 87 genes. There was a strong agreement between the two platforms as shown by comparisons of log₂ fold differences of gene expression between tumor versus normal samples (R = 0.93) and by Bland-Altman analysis, where greater than 90% of expression values fell within the 95% limits of agreement.

Conclusion

Overall concordance of gene expression patterns based on correlations, statistical significance between tumor and normal ovary data, and fold changes was consistent between the Ziplex and Affymetrix platforms. The reproducibility and ease-of-use of the technology suggests that the Ziplex array is a suitable platform for translational research.

Reprogramming of the transcriptome in a novel chromosome 3 transfer tumor suppressor ovarian cancer cell line model affected molecular networks that are characteristic of ovarian cancer

Tumor suppression as a consequence of the transfer of chromosome 3p fragments was previously observed in a novel epithelial ovarian cancer (EOC) OV-90 cell line model harboring loss of 3p. Microarray analysis revealed that tumor suppression was associated with a modified transcriptome. To investigate the relevance of the altered transcriptome, the differentially expressed genes identified by Affymetrix analysis in the 3p transfer studies, were integrated with a comparative microarray analysis of normal ovarian surface epithelial (NOSE) cells and malignant ovarian (TOV) cancers. Data from 219 significantly differentially expressed genes exhibited patterns in the direction predicted by the analysis of 3p transfer study. The 30 genes with the highest statistically significant differences (P < 1 x 10(-8)) in expression were found consistently differentially expressed between NOSE and TOV samples. The investigation of these genes in benign serous ovarian tumors and EOC cell lines also exhibited predictable expression patterns. Within the group of differentially expressed genes were SPARC, DAB2, CP, EVI1, ELF3, and EHD2, known to play a role in ovarian cancer, genes implicated in other cancers, such as GREM1 and GLIPR1, as well as genes not previously reported in a cancer context such as AKAP2 and ATAD4. A number of the differentially expressed genes are implicated in the TGF-beta signaling pathway. These findings suggest that the reprogramming of the transcriptome that occurred as a consequence of the chromosome 3 transfer and tumor suppression affected molecular networks that are characteristic of ovarian carcinogenesis thus validating our novel ovarian cancer cell line model.

Influence of monolayer, spheroid, and tumor growth conditions on chromosome 3 gene expression in tumorigenic epithelial ovarian cancer cell lines

Background

Expression microarray analyses of epithelial ovarian cancer (EOC) cell lines may be exploited to elucidate genetic and epigenetic events important in this disease. A possible variable is the influence of growth conditions on discerning candidates. The present study examined the influence of growth conditions on the expression of chromosome 3 genes in the tumorigenic EOC cell lines, OV-90, TOV-21G and TOV-112D using Affymetrix GeneChip^® HG-U133A expression microarray analysis.

Methods

Chromosome 3 gene expression profiles (n = 1147 probe sets, representing 735 genes) were extracted from U133A expression microarray analyses of the EOC cell lines OV-90, TOV-21G and TOV-112D that were grown as monolayers, spheroids or nude mouse xenografts and monolayers derived from these tumors. Hierarchical cluster analysis was performed to compare chromosome 3 transcriptome patterns of each growth condition. Differentially expressed genes were identified and characterized by two-way comparative analyses of fold-differences in gene expression between monolayer cultures and each of the other growth conditions, and between the maximum and minimum values of expression of all growth conditions for each EOC cell line.

Results

An overall high degree of similarity (> 90%) in gene expression was observed when expression values of alternative growth conditions were compared within each EOC cell line group. Two-way comparative analysis of each EOC cell line grown in an alternative condition relative to the monolayer culture showed that overall less than 15% of probe sets exhibited at least a 3-fold difference in expression profile. Less than 23% of probe sets exhibited greater than 3-fold differences in gene expression in comparisons of the maximum and minimum value of expression of all growth conditions within each EOC cell line group. The majority of these differences were less than 5-fold. There were 17 genes in common which were differentially expressed in all EOC cell lines. However, the patterns of expression of these genes were not necessarily the same for each growth condition when one cell line was compared with another.

Conclusion

The various alternative in vivo and in vitro growth conditions of tumorigenic EOC cell lines appeared to modestly influence the global chromosome 3 transcriptome supporting the notion that the in vitro cell line models are a viable option for testing gene candidates.