Clinical relevance of polymorphic DNA copy number variation (CNV) in African American women with stage I-II breast cancer

Objectives. The chief goal of this study was to analyze copy number variation (CNV) in breast cancer tumors from 25 African American women with early stage breast cancer (BC) using molecular inversion probes (MIP) in order to: (1) compare the degree of CNV in tumors compared to normal lymph nodes, and (2) determine whether gains and/or losses of genes in specific chromosomes differ between pathologic subtypes of breast cancer defined by known prognostic markers, (3) determine whether gains/losses in CN are associated with known oncogenes or tumor suppressor genes, and (4) determine whether increased gains/losses in CN for specific chromosomes were associated with differences in breast cancer recurrence. ^ Methods. Twenty to 37 nanograms of DNA extracted from 25 formalin-fixed paraffin embedded (FFPE) tumor samples and matched normal lymph nodes were added to individual tubes. Oligonucleotide probes with recognition sequences at each terminus were hybridized with a genomic target sequence to form a circular structure. Probes are released from genomic DNA obtained from FFPE samples, and those which have been correctly "circularized" in the proper allele/nucleotide reaction combination are amplified using polymerase chain reaction (PCR) primers. Amplicons were fluorescently labeled and the tag sequences released from the genome homology regions by treatment with uracil-N-glycosylase to cleave the probe at the site where uracils are present, and detected using a complementary tag array developed by Affymetrix. ^ Results. Analysis of CN gains and losses from tumors and normal tissues showed marked differences in tumors with numerous chromosomes affected. Similar changes were not observed in normal lymph nodes. When tumors were stratified into four groups based on expression or lack of expression of the estrogen receptor and HER2/neu, distinct patterns of CNV for different chromosomes were observed. Gains or losses in CN for specific chromosomes correlated with amplifications/deletions of particular oncogenes or tumor suppressor genes (i.e. such as found on chromosome 17) known to be associated with aggressive tumor phenotype and poor prognosis. There was a trend for increases in CN observed for chromosome 17 to correlate inversely with time to recurrence of BC (p=0.14 for trend). CNV was also observed for chromosomes 5, 8, 10, 11, and 16, which are known sites for several breast cancer susceptibility alleles. ^ Conclusions. This study is the first to validate the MIP technique, to correlate differences in gene expression with known prognostic tumor markers, and to correlate significant increases/decreases in CN with known tumor markers associated with prognosis. The results of this study may have far reaching public health implications towards identifying new high-risk groups based on genomic differences in CNP, both with respect to prognosis and response to therapy, and to eventually identify new therapeutic targets for prevention and treatment of this disease. ^

KEAP1 E3 ligase-mediated downregulation of NF-kappaB signaling by targeting IKKbeta.

IkappaB kinase beta (IKKbeta) is involved in tumor development and progression through activation of the nuclear factor (NF)-kappaB pathway. However, the molecular mechanism that regulates IKKbeta degradation remains largely unknown. Here, we show that a Cullin 3 (CUL3)-based ubiquitin ligase, Kelch-like ECH-associated protein 1 (KEAP1), is responsible for IKKbeta ubiquitination. Depletion of KEAP1 led to the accumulation and stabilization of IKKbeta and to upregulation of NF-kappaB-derived tumor angiogenic factors. A systematic analysis of the CUL3, KEAP1, and RBX1 genomic loci revealed a high percentage of genome loss and missense mutations in human cancers that failed to facilitate IKKbeta degradation. Our results suggest that the dysregulation of KEAP1-mediated IKKbeta ubiquitination may contribute to tumorigenesis.

The Role of Acidic Fibroblast Growth Factor and Fibroblast Growth Factor Receptor 4 in High Grade Serous Carcinoma

Background: High grade serous carcinoma whether ovarian, tubal or primary peritoneal, continues to be the most lethal gynecologic malignancy in the USA. Although combination chemotherapy and aggressive surgical resection has improved survival in the past decade the majority of patients still succumb to chemo-resistant disease recurrence. It has recently been reported that amplification of 5q31-5q35.3 is associated with poor prognosis in patients with high grade serous ovarian carcinoma. Although the amplicon contains over 50 genes, it is notable for the presence of several members of the fibroblast growth factor signaling axis. In particular acidic fibroblast growth factor (FGF1) has been demonstrated to be one of the driving genes in mediating the observed prognostic effect of the amplicon in ovarian cancer patients. This study seeks to further validate the prognostic value of fibroblast growth receptor 4 (FGFR4), another candidate gene of the FGF/FGFR axis located in the same amplicon. The emphasis will be delineating the role the FGF1/FGFR4 signaling axis plays in high grade serous ovarian carcinoma; and test the feasibility of targeting the FGF1/FGFR4 axis therapeutically. Materials and Methods: Spearman and Pearson correlation studies on data generated from array CGH and transcriptome profiling analyses on 51 microdissected tumor samples were used to identify genes located on chromosome 5q31-35.3 that showed significant correlation between DNA and mRNA copy numbers. Significant correlation between FGF1 and FGFR4 DNA copy numbers was further validated by qPCR analysis on DNA isolated from 51 microdissected tumor samples. Immunolocalization and quantification of FGFR4 expression were performed on paraffin embedded tissue samples from 183 cases of high-grade serous ovarian carcinoma. The expression was then correlated with clinical data to assess impact on survival. The expression of FGF1 and FGFR4 in vitro was quantified by real-time PCR and western blotting in six high-grade serous ovarian carcinoma cell lines and compared to those in human ovarian surface epithelial cells to identify overexpression. The effect of FGF1 on these cell lines after serum starvation was quantified for in vitro cellular proliferation, migration/invasion, chemoresistance and survival utilizing a combination of commercially available colorimetric, fluorometric and electrical impedance assays. FGFR4 expression was then transiently silenced via siRNA transfection and the effects on response to FGF1, cellular proliferation, and migration were quantified. To identify relevant cellular pathways involved, responsive cell lines were transduced with different transcription response elements using the Cignal-Lenti reporter system and treated with FGF1 with and without transient FGFR4 knock down. This was followed by western blot confirmation for the relevant phosphoproteins. Anti-FGF1 antibodies and FGFR trap proteins were used to attempt inhibition of FGF mediated phenotypic changes and relevant signaling in vitro. Orthotopic intraperitoneal tumors were established in nude mice using serous cell lines that have been previously transfected with luciferase expressing constructs. The mice were then treated with FGFR trap protein. Tumor progression was then followed via bioluminescent imaging. The FGFR4 gene from 52 clinical samples was sequenced to screen for mutations. Results: FGFR4 DNA and mRNA copy numbers were significantly correlated and FGFR4 DNA copy number was significantly correlated with that of FGF1. Survival of patients with high FGFR4 expressing tumors was significantly shorter that those with low expression(median survival 28 vs 55 month p< 0.001) In a multivariate cox regression model FGFR expression significantly increased risk of death (HR 2.1, p<0.001). FGFR4 expression was significantly higher in all cell lines tested compared to HOSE, OVCA432 cell line in particular had very high expression suggesting amplification. FGF1 was also particularly overexpressed in OVCA432. FGF1 significantly increased cell survival after serum deprivation in all cell lines. Transient knock down of FGFR4 caused significant reduction in cell migration and proliferation in vitro and significantly decreased the proliferative effects of FGF1 in vitro. FGFR1, FGFR4 traps and anti-FGF1 antibodies did not show activity in vitro. OVCA432 transfected with the cignal lenti reporter system revealed significant activation of MAPK, NFkB and WNT pathways, western blotting confirmed the results. Reverse phase protein array (RPPA) analysis also showed activation of MAPK, AKT, WNT pathways and down regulation of E Cadherin. FGFR trap protein significantly reduced tumor growth in vivo in an orthotopic mouse model. Conclusions: Overexpression and amplification of several members of the FGF signaling axis present on the amplicon 5q31-35.3 is a negative prognostic indicator in high grade serous ovarian carcinoma and may drive poor survival associated with that amplicon. Activation of The FGF signaling pathway leads to downstream activation of MAPK, AKT, WNT and NFkB pathways leading to a more aggressive cancer phenotype with increased tumor growth, evasion of apoptosis and increased migration and invasion. Inhibition of FGF pathway in vivo via FGFR trap protein leads to significantly decreased tumor growth in an orthotopic mouse model.

Genome-wide algorithm for detecting CNV associations with diseases

SNP genotyping arrays have been developed to characterize single-nucleotide polymorphisms (SNPs) and DNA copy number variations (CNVs). The quality of the inferences about copy number can be affected by many factors including batch effects, DNA sample preparation, signal processing, and analytical approach. Nonparametric and model-based statistical algorithms have been developed to detect CNVs from SNP genotyping data. However, these algorithms lack specificity to detect small CNVs due to the high false positive rate when calling CNVs based on the intensity values. Association tests based on detected CNVs therefore lack power even if the CNVs affecting disease risk are common. In this research, by combining an existing Hidden Markov Model (HMM) and the logistic regression model, a new genome-wide logistic regression algorithm was developed to detect CNV associations with diseases. We showed that the new algorithm is more sensitive and can be more powerful in detecting CNV associations with diseases than an existing popular algorithm, especially when the CNV association signal is weak and a limited number of SNPs are located in the CNV.^

THEORETICAL STUDIES ON THE METHODS OF RECONSTRUCTING PHYLOGENETIC TREES FROM DNA SEQUENCE DATA (MOLECULAR EVOLUTION, HOMINOID EVOLUTION, COMPUTER SIMULATION)

(1) A mathematical theory for computing the probabilities of various nucleotide configurations is developed, and the probability of obtaining the correct phylogenetic tree (model tree) from sequence data is evaluated for six phylogenetic tree-making methods (UPGMA, distance Wagner method, transformed distance method, Fitch-Margoliash's method, maximum parsimony method, and compatibility method). The number of nucleotides (m*) necessary to obtain the correct tree with a probability of 95% is estimated with special reference to the human, chimpanzee, and gorilla divergence. m* is at least 4,200, but the availability of outgroup species greatly reduces m* for all methods except UPGMA. m* increases if transitions occur more frequently than transversions as in the case of mitochondrial DNA. (2) A new tree-making method called the neighbor-joining method is proposed. This method is applicable either for distance data or character state data. Computer simulation has shown that the neighbor-joining method is generally better than UPGMA, Farris' method, Li's method, and modified Farris method on recovering the true topology when distance data are used. A related method, the simultaneous partitioning method, is also discussed. (3) The maximum likelihood (ML) method for phylogeny reconstruction under the assumption of both constant and varying evolutionary rates is studied, and a new algorithm for obtaining the ML tree is presented. This method gives a tree similar to that obtained by UPGMA when constant evolutionary rate is assumed, whereas it gives a tree similar to that obtained by the maximum parsimony tree and the neighbor-joining method when varying evolutionary rate is assumed. ^

EXPRESSION OF THE HUMAN PLACENTAL LACTOGEN GENES

The placenta is the site of synthesis of various peptide and steroid hormones related to pregnancy. Human placental lactogen (hPL) is the predominant peptide hormone secreted by term placenta and its synthesis is tissue-specific and coupled to placenta development. The objective of this work was to study the structure and expression of the hPL.^ Poly(A('+))RNA from human term placenta was translated in a mouse-derived cell-free system. A major band corresponding to pre-hPL and a minor band comigrating with mature hPL, represent (TURN)15% of the total radioactively labeled proteins. Analysis of the poly(A('+))RNA showed a prominent band at approximately 860 nucleotides. A corresponding band was observed in Northern blots of total RNA, hybridized with {('32)P}-labeled recombinant plasmid containing a portion of hPL cDNA. Similar analyses of nuclear RNA showed at least four additional bands at 990, 1200, 1460 and 1760 nucleotides, respectively, which are likely precursors of hPL mRNA. Poly(A('+))RNA was used to construct a cDNA library, of which approximately 5% of the clones were found to hybridize to hPL DNA sequences. Heteroduplexes constructed between a clone containing a 815 bp hPL cDNA insert and a hPL genomic DNA clone revealed four small intervening sequences which can account for the lengths observed in hnRNA molecules.^ Recombinant plasmid HCS-pBR322 containing a 550 bp insert of a cDNA transcript of human placental lactogen (hPL) mRNA was ('3)H-labeled an hybridized in situ to human chromosome preparations. These experiments allowed assignment of the hPL and growth hormone (hGH) genes, which have over 90% nucleotide homology in their coding sequences, to band q22-24 of chromosome 17. A gene copy number experiment showed that both genes are present in (TURN)3 copies per haploid genome.^ Experiments were designed to determine if all members of the hPL gene cluster, consisting of four non-allelic genes, are transcribed in term placenta. Advantage was taken of differences in restriction endonuclease sites in the coding portions of the different hPL genes, to distinguish the putative cDNAs of the transcriptionally active genes. Two genes were found to be represented in the cDNA library and their cDNA transcripts were isolated and characterized. Three independent methods showed that their corresponding mRNAs are about equally represented in the hPL mRNA population. The two cDNAs code for prehPL proteins which differ at a single amino acid position. However the secreted hPLs have identical amino acid sequences. A tetramer insertion duplication was found in a palindrome area of the 3' untranslated region of one of the hPL mRNAs. ^

Molecular alterations in nucleotide excision repair genes in malignant glioma and their relevance to malignant progression and drug resistance

Recently, it has become apparent that DNA repair mechanisms are involved in the malignant progression and resistance to therapy of gliomas. Many investigators have shown that increased levels of O6-methyl guanine DNA alkyltransferase, a DNA monoalkyl adduct repair enzyme, are correlated with resistance of malignant glioma cell lines to nitrosourea-based chemotherapy. Three important DNA excision repair genes ERCC1 (excision repair cross complementation group 1), ERCC2 (excision repair cross complementation group 2), and ERCC6 (excision repair cross complementation group 6) have been studied in human tumors. Gene copy number variation of ERCC1 and ERCC2 has been observed in primary glioma tissues. A number of reports describing a relationship between ERCC1 gene alterations and resistance to anti-cancer drugs have been also described. The levels of ERCC1 gene expression, however, have not been correlated with drug resistance in gliomas. The expression of ERCC6 gene transcribes has been shown to vary with tissue types and to be highest in the brain. There have been no comprehensive studies so far, however, of ERCC6 gene expression and molecular alterations in malignant glioma. This project examined the ERCC1 expression levels and correlated them with cisplatin resistance in malignant glioma cell lines. We also examined the molecular alterations of ERCC6 gene in primary glioma tissues and cells and analyzed whether these alterations are related to tumor progression and chemotherapy resistance. Our results indicate the presence of mutations and/or deletions in exons II and V of the ERCC6 gene, and these alterations are more frequent in exon II. Furthermore, the mutations and/or deletions in exon II were shown to be associated with increased malignant grade of gliomas. The results on the Levels of ERCC1 gene transcripts showed that expression levels correlate with cisplatin resistance. The increase in ERCC1 mRNA induced by cisplatin could be down-regulated by cyclosporin A and herbimycin A. The results of this study are likely to provide useful information for clinical treatment of human gliomas. ^

Mitochondrial defects in primary leukemia cells: Biological consequences and therapeutic implications

Mitochondria are actively engaged in the production of cellular energy sources, generation of reactive oxygen species (ROS), and regulation of apoptosis. Mitochondrial DNA (mtDNA) mutations/deletions and other mitochondrial abnormalities have been implicated in many diseases, especially cancer. Despite this, the roles that these defects play in cancer development, drug sensitivity, and disease progression still remain to be elucidated. The major objective of this investigation was to evaluate the mechanistic relationship between mitochondrial defects and alterations in free radical generation and chemosensitivity in primary chronic lymphocytic leukemia (CLL) cells. This study revealed that the mtDNA mutation frequency and basal superoxide generation are both significantly higher in primary cells from CLL patients with a history of chemotherapy as compared to cells from their untreated counterparts. CLL cells from refractory patients tended to have high mutation frequencies. The data suggest that chemotherapy with DNA-damaging agents may cause mtDNA mutations, which are associated with increased ROS generation and reduced drug sensitivity. Subsequent analyses demonstrated that CLL cells contain significantly more mitochondria than normal lymphocytes. This abnormal accumulation of mitochondria was linked to increased expression of nuclear respiratory factor-1 and mitochondrial transcription factor A, two key free radical-regulated mitochondrial biogenesis factors. Further analysis showed that mitochondrial content may have therapeutic implications since patient cells with high mitochondrial mass display significantly reduced in vitro sensitivity to fludarabine, a frontline agent in CLL therapy. The reduced in vitro and in vivo sensitivity to fludarabine observed in CLL cells with mitochondrial defects highlights the need for novel therapeutic strategies for the treatment of refractory disease. Brefeldin A, an inhibitor of endoplasmic reticulum (ER) to Golgi protein transport that is being developed as an anticancer agent, effectively induces apoptosis in fludarabine-refractory CLL cells through a secretory stress-mediated mechanism involving intracellular sequestration of pro-survival secretory factors. Taken together, these data indicate that mitochondrial defects in CLL cells are associated with alterations in free radical generation, mitochondrial biogenesis activity, and chemosensitivity. Abrogation of survival signaling by blocking ER to Golgi protein transport may be a promising therapeutic strategy for the treatment of CLL patients that respond poorly to conventional chemotherapy. ^

High-resolution microarray analysis of chromosome 20q in human colon cancer metastasis model systems

Amplification of human chromosome 20q DNA is the most frequently occurring chromosomal abnormality detected in sporadic colorectal carcinomas and shows significant correlation with liver metastases. Through comprehensive high-resolution microarray comparative genomic hybridization and microarray gene expression profiling, we have characterized chromosome 20q amplicon genes associated with human colorectal cancer metastasis in two in vitro metastasis model systems. The results revealed increasing complexity of the 20q genomic profile from the primary tumor-derived cell lines to the lymph node and liver metastasis derived cell lines. Expression analysis of chromosome 20q revealed a subset of over expressed genes residing within the regions of genomic copy number gain in all the tumor cell lines, suggesting these are Chromosome 20q copy number responsive genes. Bases on their preferential expression levels in the model system cell lines and known biological function, four of the over expressed genes mapping to the common intervals of genomic copy gain were considered the most promising candidate colorectal metastasis-associated genes. Validation of genomic copy number and expression array data was carried out on these genes, with one gene, DNMT3B, standing out as expressed at a relatively higher levels in the metastasis-derived cell lines compared with their primary-derived counterparts in both the models systems analyzed. The data provide evidence for the role of chromosome 20q genes with low copy gain and elevated expression in the clonal evolution of metastatic cells and suggests that such genes may serve as early biomarkers of metastatic potential. The data also support the utility of the combined microarray comparative genomic hybridization and expression array analysis for identifying copy number responsive genes in areas of low DNA copy gain in cancer cells. ^

Mitochondrial dysfunction leads to NOX activation: A novel mechanism to maintain high glycolysis in cancer cells

Increased glycolysis and oxidative stress are common features of cancer cells. These metabolic alterations are associated with mitochondrial dysfunction and can be caused by mitochondrial DNA (mtDNA) mutations, oncogenic signals, loss of tumor suppressor, and tumor tissue hypoxia. It is well established that mitochondria play central roles in energy metabolism, maintenance of redox balance, and regulation of apoptosis. However, the biochemical and molecular mechanisms that maintain high glycolysis in cancer cells (the Warburg effect) with mitochondrial dysfunction and oxidative stress remain to be determined. The major goals of this study were to establish a unique experimental system in which the mitochondrial respiratory function can be regulated as desired, and to use this system to investigate the mechanistic link between mitochondrial dysfunction and the Warburg effect along with oxidative stress in cancer cells. To achieve these goals, I have established a tetracycline-inducible system in which a dominant negative form of mitochondrial DNA polymerase y (POLGdn) expression could be regulated by tetracycline; thus controlling mitochondrial respiratory function. Using this cell system, I demonstrated that POLGdn expression resulted in mitochondrial dysfunction through decreasing mtDNA content, depletion of mtDNA encoded mRNA and protein expression. This process was mediated by TFAM proteasome degradation. Mitochondrial dysfunction mediated by POLGdn expression led to a significant increase in cellular glycolysis and oxidative stress. Surprisingly, mitochondrial dysfunction also resulted in increased NAD(P)H oxidase (NOX) enzyme activity, which was shown to be essential for maintaining high glycolysis. Chemical Inhibition of NOX activity by diphenyliodonium (DPI) preferentially impacted the survival of mitochondrial defective cells. The colon cancer HCT116-/- cells that have lost transcriptional regulation of the mitochondrial assembling enzyme SCO2, leading to compromised mitochondrial respiratory function, were found to have increased NOX activity and were highly sensitive to DPI treatment. Ovarian epithelial cells with Ras transformation also exhibited an increase in NOX gene expression and NOX enzyme activity, rendering the cells sensitive to DPI inhibition especially under hypoxic condition. These data together suggest that NOX plays a novel role in maintaining high glycolysis in cancer cells with mitochondrial defects, and that NOX may be a potential target for cancer therapy. ^

Mechanisms of HER2/{\it neu\/} proto-oncogene overexpression in breast cancer cells

Overexpression and amplification of HER2/neu have been documented in many primary tumors, most notably in breast. Not only do approximately 30% of breast cancer patients carry tumors that overexpress the gene, but those that do generally have shorter overall and disease-free survival times than patients with tumors expressing low levels of HER2/neu. Thus, overexpression of HER2/neu plays an important role in the pathogenesis of breast cancer. We have examined the mechanisms that result in HER2/neu overexpression in breast cancer by using, as a model system, established breast cancer cell lines that express much higher levels of HER2/neu mRNA than normal breast tissue while maintaining a near normal HER2/neu gene copy number. Nuclear run-on experiments indicate that the breast cancer cell lines MDA-MB453, BT483, and BT474 have an increased HER2/neu gene transcription rate. By using HER2/neu promoter-CAT constructs, we have found that the enhanced HER2/neu transcription rate in MDA-MB453 cells is due to activation of the gene in trans, while the enhanced transcription rate in BT483 cells is due to activation of the gene in either trans or cis. In BT474 cells, transcriptional upregulation is primarily due to gene amplification. Since the levels of increased transcription are not as high as the levels of HER2/neu mRNA in any of these three lines, post-transcriptional deregulation that increases HER2/neu expression must also be functioning in these cells. The half-life of HER2/neu mRNA was measured and found to be equivalent in these lines as in a control. Thus, the post-transcriptional deregulation is not increased stability of the HER2/neu transcript.^ Much work has been performed in characterizing the altered trans-acting factor involved in increased HER2/neu transcription in MDA-MB453 cells. Using promoter deletion constructs linked to a reporter gene, the region responsive to this factor was localized in the rat neu promoter. When human HER2/neu promoter constructs were used, the homologous sequence in the human promoter was identified. Furthermore, a number of protein/DNA complexes are detected when these promoter regions are used in gel mobility shift assays. UV-crosslinking experiments indicate DNA-binding proteins of roughly 110 kDa, 70 kDa, and 35 kDa are capable of interacting with the human promoter element. ^