A transcriptomic computational analysis of mastic oil-treated Lewis lung carcinomas reveals molecular mechanisms targeting tumor cell growth and survival

Background

DNA copy number alterations are frequently observed in ovarian cancer, but it remains a challenge to identify the most relevant alterations and the specific causal genes in those regions.

Methods

We obtained high-resolution 500K SNP array data for 52 ovarian tumors and identified the most statistically significant minimal genomic regions with the most prevalent and highest-level copy number alterations (recurrent CNAs). Within a region of recurrent CNA, comparison of expression levels in tumors with a given CNA to tumors lacking that CNA and to whole normal ovary samples was used to select genes with CNA-specific expression patterns. A public expression array data set of laser capture micro-dissected (LCM) non-malignant fallopian tube epithelia and LCM ovarian serous adenocarcinoma was used to evaluate the effect of cell-type mixture biases.

Results

Fourteen recurrent deletions were detected on chromosomes 4, 6, 9, 12, 13, 15, 16, 17, 18, 22 and most prevalently on X and 8. Copy number and expression data suggest several apoptosis mediators as candidate drivers of the 8p deletions. Sixteen recurrent gains were identified on chromosomes 1, 2, 3, 5, 8, 10, 12, 15, 17, 19, and 20, with the most prevalent gains localized to 8q and 3q. Within the 8q amplicon, PVT1, but not MYC, was strongly over-expressed relative to tumors lacking this CNA and showed over-expression relative to normal ovary. Likewise, the cell polarity regulators PRKCI and ECT2 were identified as putative drivers of two distinct amplicons on 3q. Co-occurrence analyses suggested potential synergistic or antagonistic relationships between recurrent CNAs. Genes within regions of recurrent CNA showed an enrichment of Cancer Census genes, particularly when filtered for CNA-specific expression.

Conclusion

These analyses provide detailed views of ovarian cancer genomic changes and highlight the benefits of using multiple reference sample types for the evaluation of CNA-specific expression changes.     

Spatial and temporal evolution of distal 10q deletion,a prognostically unfavorable event in diffuse low-grade gliomas

Background

The disease course of patients with diffuse low-grade glioma is notoriously unpredictable. Temporal and spatially distinct samples may provide insight into the evolution of clinically relevant copy number aberrations (CNAs). The purpose of this study is to identify CNAs that are indicative of aggressive tumor behavior and can thereby complement the prognostically favorable 1p/19q co-deletion.

Results

Genome-wide, 50 base pair single-end sequencing was performed to detect CNAs in a clinically well-characterized cohort of 98 formalin-fixed paraffin-embedded low-grade gliomas. CNAs are correlated with overall survival as an endpoint. Seventy-five additional samples from spatially distinct regions and paired recurrent tumors of the discovery cohort were analyzed to interrogate the intratumoral heterogeneity and spatial evolution. Loss of 10q25.2-qter is a frequent subclonal event and significantly correlates with an unfavorable prognosis. A significant correlation is furthermore observed in a validation set of 126 and confirmation set of 184 patients. Loss of 10q25.2-qter arises in a longitudinal manner in paired recurrent tumor specimens, whereas the prognostically favorable 1p/19q co-deletion is the only CNA that is stable across spatial regions and recurrent tumors.

Conclusions

CNAs in low-grade gliomas display extensive intratumoral heterogeneity. Distal loss of 10q is a late onset event and a marker for reduced overall survival in low-grade glioma patients. Intratumoral heterogeneity and higher frequencies of distal 10q loss in recurrences suggest this event is involved in outgrowth to the recurrent tumor.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0471-6) contains supplementary material, which is available to authorized users.     

Latent Factor Analysis to Discover Pathway-Associated Putative Segmental Aneuploidies in Human Cancers

Tumor microenvironmental stresses, such as hypoxia and lactic acidosis, play important roles in tumor progression. Although gene signatures reflecting the influence of these stresses are powerful approaches to link expression with phenotypes, they do not fully reflect the complexity of human cancers. Here, we describe the use of latent factor models to further dissect the stress gene signatures in a breast cancer expression dataset. The genes in these latent factors are coordinately expressed in tumors and depict distinct, interacting components of the biological processes. The genes in several latent factors are highly enriched in chromosomal locations. When these factors are analyzed in independent datasets with gene expression and array CGH data, the expression values of these factors are highly correlated with copy number alterations (CNAs) of the corresponding BAC clones in both the cell lines and tumors. Therefore, variation in the expression of these pathway-associated factors is at least partially caused by variation in gene dosage and CNAs among breast cancers. We have also found the expression of two latent factors without any chromosomal enrichment is highly associated with 12q CNA, likely an instance of “trans”-variations in which CNA leads to the variations in gene expression outside of the CNA region. In addition, we have found that factor 26 (1q CNA) is negatively correlated with HIF-1α protein and hypoxia pathways in breast tumors and cell lines. This agrees with, and for the first time links, known good prognosis associated with both a low hypoxia signature and the presence of CNA in this region. Taken together, these results suggest the possibility that tumor segmental aneuploidy makes significant contributions to variation in the lactic acidosis/hypoxia gene signatures in human cancers and demonstrate that latent factor analysis is a powerful means to uncover such a linkage.     

Chromosomal imbalances in primary and metastatic melanomas revealed by comparative genomic hybridization

Characteristic genetic changes underlying the metastatic progression of malignant melanoma is incompletely understood. The goal of our study was to explore specific chromosomal alterations associated with the aggressive behavior of this neoplasm. Comparative genomic hybridization was performed to screen and compare genomic imbalances present in primary and metastatic melanomas. Sixteen primary and 12 metastatic specimens were analyzed. We found that the pattern of chromosomal aberrations is similar in the two subgroups; however, alterations present only in primary and/or metastatic tumors were also discovered. The mean number of genetic changes was 6.3 (range 1-14) in primary and 7.8 (range 1-16) in metastatic lesions. Frequent losses involved 9p and 10q, whereas gains most often occurred at 1q, 6p, 7q, and 8q. Distinct, high-level amplifications were mapped to 1p12-p21 and 1p22-p31 in both tumor types. Amplification of 4q12-q13.1, 7q21.3-qter and 8q23-qter were detected only in primary tumors. The 20q13-qter amplicon was present in a metastatic tumor. The number of genetic alterations were significantly higher in primary tumors which developed metastases within one year after the surgery compared to tumors without metastasis during this time period. Fluorescence in situ hybridization with centromeric and locus-specific probes was applied to validate CGH results on a subset of tumors. Comparison of FISH and CGH data gave good correlation. The aggressive behavior of melanoma is associated with accumulation of multiple genetic alterations. Chromosome regions, which differ in the primary and metastatic lesions, may represent potential targets to identify metastases-related chromosomal alterations.     

Methylation microarray-based detection of clinical copy-number aberrations in CLL benchmarked to standard FISH analysis

Copy-number aberrations (CNAs) are assessed using FISH analysis in diagnostics of chronic lymphocytic leukemia (CLL), but CNAs can also be extrapolated from Illumina BeadChips developed for genome-wide methylation microarray screening. Increasing numbers of microarray data-sets are available from diagnostic samples, making it useful to assess the potential in CNA diagnostics.We benchmarked the limitations of CNA testing from two Illumina BeadChips (EPIC and 450k) and using two common packages for analysis (conumee and ChAMP) to FISH-based assessment of 11q, 13q, and 17p deletions in 202 CLL samples.Overall, the two packages predicted CNAs with similar accuracy regardless of the microarray type, but lower than FISH-based assessment. We showed that the bioinformatics analysis needs to be adjusted to the specific CNA, as no general settings were identified. Altogether, we were able to predict CNAs using methylation microarray data, however, with limited accuracy, making FISH-based assessment of deletions the superior diagnostic choice.     

High Resolution Genome-Wide Analysis of Chromosomal Alterations in Burkitt's Lymphoma

Additional chromosomal abnormalities are currently detected in Burkitt''s lymphoma. They play major roles in the progression of BL and in prognosis. The genes involved remain elusive. A whole-genome oligonucleotide array CGH analysis correlated with karyotype and FISH was performed in a set of 27 Burkitt''s lymphoma-derived cell lines and primary tumors. More than half of the 145 CNAs<2 Mb were mapped to Mendelian CNVs, including GSTT1, glutathione s-transferase and BIRC6, an anti-apoptotic protein, possibly predisposing to some cancers. Somatic cell line-specific CNVs localized to the IG locus were consistently observed with the 244 K aCGH platform. Among 136 CNAs >2 Mb, gains were found in 1q (12/27), 13q (7/27), 7q (6/27), 8q(4/27), 2p (3/27), 11q (2/27) and 15q (2/27). Losses were found in 3p (5/27), 4p (4/27), 4q (4/27), 9p (4/27), 13q (4/27), 6p (3/27), 17p (3/27), 6q (2/27),11pterp13 (2/27) and 14q12q21.3 (2/27). Twenty one minimal critical regions (MCR), (range 0.04–71.36 Mb), were delineated in tumors and cell lines. Three MCRs were localized to 1q. The proximal one was mapped to 1q21.1q25.2 with a 6.3 Mb amplicon (1q21.1q21.3) harboring BCA2 and PIAS3. In the other 2 MCRs, 1q32.1 and 1q44, MDM4 and AKT3 appeared as possible drivers of these gains respectively. The 13q31.3q32.1 <89.58–96.81> MCR contained an amplicon and ABCC4 might be the driver of this amplicon. The 40 Kb 2p16.1 <60.96–61> MCR was the smallest gained MCR and specifically encompassed the REL oncogene which is already implicated in B cell lymphomas. The most frequently deleted MCR was 3p14.1 <60.43–60.53> that removed the fifth exon of FHIT. Further investigations which combined gene expression and functional studies are essential to understand the lymphomagenesis mechanism and for the development of more effective, targeted therapeutic strategies.     

Genomic profiling of submucosal-invasive gastric cancer by array-based comparative genomic hybridization

Genomic copy number aberrations (CNAs) in gastric cancer have already been extensively characterized by array comparative genomic hybridization (array CGH) analysis. However, involvement of genomic CNAs in the process of submucosal invasion and lymph node metastasis in early gastric cancer is still poorly understood. In this study, to address this issue, we collected a total of 59 tumor samples from 27 patients with submucosal-invasive gastric cancers (SMGC), analyzed their genomic profiles by array CGH, and compared them between paired samples of mucosal (MU) and submucosal (SM) invasion (23 pairs), and SM invasion and lymph node (LN) metastasis (9 pairs). Initially, we hypothesized that acquisition of specific CNA(s) is important for these processes. However, we observed no significant difference in the number of genomic CNAs between paired MU and SM, and between paired SM and LN. Furthermore, we were unable to find any CNAs specifically associated with SM invasion or LN metastasis. Among the 23 cases analyzed, 15 had some similar pattern of genomic profiling between SM and MU. Interestingly, 13 of the 15 cases also showed some differences in genomic profiles. These results suggest that the majority of SMGCs are composed of heterogeneous subpopulations derived from the same clonal origin. Comparison of genomic CNAs between SMGCs with and without LN metastasis revealed that gain of 11q13, 11q14, 11q22, 14q32 and amplification of 17q21 were more frequent in metastatic SMGCs, suggesting that these CNAs are related to LN metastasis of early gastric cancer. In conclusion, our data suggest that generation of genetically distinct subclones, rather than acquisition of specific CNA at MU, is integral to the process of submucosal invasion, and that subclones that acquire gain of 11q13, 11q14, 11q22, 14q32 or amplification of 17q21 are likely to become metastatic.     

A meta-analysis of array-CGH studies implicates antiviral immunity pathways in the development of hepatocellular carcinoma

Background

The development and progression of hepatocellular carcinoma (HCC) is significantly correlated to the accumulation of genomic alterations. Array-based comparative genomic hybridization (array CGH) has been applied to a wide range of tumors including HCCs for the genome-wide high resolution screening of DNA copy number changes. However, the relevant chromosomal variations that play a central role in the development of HCC still are not fully elucidated.

Methods

In present study, in order to further characterize the copy number alterations (CNAs) important to HCC development, we conducted a meta-analysis of four published independent array-CGH datasets including total 159 samples.

Results

Eighty five significant gains (frequency ≥25%) were mostly mapped to five broad chromosomal regions including 1q, 6p, 8q, 17q and 20p, as well as two narrow regions 5p15.33 and 9q34.2-34.3. Eighty eight significant losses (frequency ≥25%) were most frequently present in 4q, 6q, 8p, 9p, 13q, 14q, 16q, and 17p. Significant correlations existed between chromosomal aberrations either located on the same chromosome or the different chromosomes. HCCs with different etiologies largely exhibited surprisingly similar profiles of chromosomal aberrations with only a few exceptions. Furthermore, the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis indicated that the genes affected by these chromosomal aberrations were significantly enriched in 31 canonical pathways with the highest enrichment observed for antiviral immunity pathways.

Conclusions

Taken together, our findings provide novel and important clues for the implications of antiviral immunity-related gene pathways in the pathogenesis and progression of HCC.     

Akt-Induced Phosphorylation of N-CoR at Serine 1450 Contributes to Its Misfolded Conformational Dependent Loss (MCDL) in Acute Myeloid Leukemia of the M5 Subtype

The nuclear receptor co-repressor (N-CoR) is a key component of the generic co-repressor complex that plays an important role in the control of cellular growth and differentiation. As shown by us recently, the growth suppressive function of N-CoR largely relies on its capacity to repress Flt3, a key regulator of cellular gorwth during normal and malignant hematopoesis. We further demonstrated how de-repression of Flt3 due to the misfolded conformation dependent loss (MCDL) of N-CoR contributed to malignant growth in acute myeloid leukemia (AML). However, the molecular mechanism underlying the MCDL of N-CoR and its implication in AML pathogenesis is not fully understood. Here, we report that Akt-induced phosphorylation of N-CoR at the consensus Akt motif is crucial for its misfolding and subsequent loss in AML (AML-M5). N-CoR displayed significantly higher level of serine specific phosphorylation in almost all AML-M5 derived cells and was subjected to processing by AML-M5 specific aberrant protease activity. To identify the kinase linked to N-CoR phosphorylation, a library of activated kinases was screened with the extracts of AML cells; leading to the identification of Akt as the putative kinase linked to N-CoR phosphorylation. Consistent with this finding, a constitutively active Akt consistently phosphorylated N-CoR leading to its misfolding; while the therapeutic and genetic ablation of Akt largely abrogated the MCDL of N-CoR in AML-M5 cells. Site directed mutagenic analysis of N-CoR identified serine 1450 as the crucial residue whose phosphorylation by Akt was essential for the misfolding and loss of N-CoR protein. Moreover, Akt-induced phosphorylation of N-CoR contributed to the de-repression of Flt3, suggesting a cross talk between Akt signaling and N-CoR misfolding pathway in the pathogenesis of AML-M5. The N-CoR misfolding pathway could be the common downstream thread of pleiotropic Akt signaling activated by various oncogenic insults in some subtypes of leukemia and solid tumors.     

Copy Number Analysis Identifies Novel Interactions Between Genomic Loci in Ovarian Cancer

Ovarian cancer is a heterogeneous disease displaying complex genomic alterations, and consequently, it has been difficult to determine the most relevant copy number alterations with the scale of studies to date. We obtained genome-wide copy number alteration (CNA) data from four different SNP array platforms, with a final data set of 398 ovarian tumours, mostly of the serous histological subtype. Frequent CNA aberrations targeted many thousands of genes. However, high-level amplicons and homozygous deletions enabled filtering of this list to the most relevant. The large data set enabled refinement of minimal regions and identification of rare amplicons such as at 1p34 and 20q11. We performed a novel co-occurrence analysis to assess cooperation and exclusivity of CNAs and analysed their relationship to patient outcome. Positive associations were identified between gains on 19 and 20q, gain of 20q and loss of X, and between several regions of loss, particularly 17q. We found weak correlations of CNA at genomic loci such as 19q12 with clinical outcome. We also assessed genomic instability measures and found a correlation of the number of higher amplitude gains with poorer overall survival. By assembling the largest collection of ovarian copy number data to date, we have been able to identify the most frequent aberrations and their interactions.     

A classification model for distinguishing copy number variants from cancer-related alterations

Background  

Both somatic copy number alterations (CNAs) and germline copy number variants (CNVs) that are prevalent in healthy individuals can appear as recurrent changes in comparative genomic hybridization (CGH) analyses of tumors. In order to identify important cancer genes CNAs and CNVs must be distinguished. Although the Database of Genomic Variants (DGV) contains a list of all known CNVs, there is no standard methodology to use the database effectively.     

NPM1 Deletion Is Associated with Gross Chromosomal Rearrangements in Leukemia

Background

NPM1 gene at chromosome 5q35 is involved in recurrent translocations in leukemia and lymphoma. It also undergoes mutations in 60% of adult acute myeloid leukemia (AML) cases with normal karyotype. The incidence and significance of NPM1 deletion in human leukemia have not been elucidated.

Methodology and Principal Findings

Bone marrow samples from 145 patients with myelodysplastic syndromes (MDS) and AML were included in this study. Cytogenetically 43 cases had isolated 5q-, 84 cases had 5q- plus other changes and 18 cases had complex karyotype without 5q deletion. FISH and direct sequencing investigated the NPM1 gene. NPM1 deletion was an uncommon event in the “5q- syndrome” but occurred in over 40% of cases with high risk MDS/AML with complex karyotypes and 5q loss. It originated from large 5q chromosome deletions. Simultaneous exon 12 mutations were never found. NPM1 gene status was related to the pattern of complex cytogenetic aberrations. NPM1 haploinsufficiency was significantly associated with monosomies (p<0.001) and gross chromosomal rearrangements, i.e., markers, rings, and double minutes (p<0.001), while NPM1 disomy was associated with structural changes (p = 0.013). Interestingly, in complex karyotypes with 5q- TP53 deletion and/or mutations are not specifically associated with NPM1 deletion.

Conclusions and Significance

NPM1/5q35 deletion is a consistent event in MDS/AML with a 5q-/-5 in complex karyotypes. NPM1 deletion and NPM1 exon 12 mutations appear to be mutually exclusive and are associated with two distinct cytogenetic subsets of MDS and AML.     

Loss of 9p21 is embedded in a complex but consistent pattern of genomic imbalances in oral squamous cell carcinomas

35 oral squamous cell carcinomas examined previously by comparative genomic hybridization (CGH) exhibited 5 up to 47 copy number alterations (CNAs). 13 of those cases showed a loss of parts of the short arm of chromosome 9, band p21 being affected in all of these cases. A highly complex but strikingly consistent pattern of genomic imbalances with an average 31.5 CNAs per tumor was associated with this deletion, and gains clearly dominated over losses of genomic material. Comparable patterns, however, could also be found in tumors with a high number of CNAs (24 CNAs) but without the deletion. Low numbers of imbalances were accompanied by low consistency of the CNA patterns. None of these latter cases showed the deletion 9p21. 66.7% of the dim(9p21)-positive tumors were of class pT4 (vs. 22% in dim(9p21)-negative cases), 77% of stage III or IV (vs. 47% in the group without the deletion), but only 8% of the dim(9p21)-positive tumors were classified as grade 3 (vs. 41% in the negative group). Other clinicopathologic features like prevalence of relapse, or survival time could not be as clearly associated with the deletion. For instance, short relapse-free survival was clearly associated with a high number of CNAs, rather independent of presence or absence of dim(9p21) in the affected tumor. From these findings it is concluded that previously found associations of 9p21 deletion with clinical parameters can reasonably be estimated only in the context of the pattern and complexity of the genomic imbalances accompanying this chromosomal loss in the examined tumors.     

DNA profiling by array comparative genomic hybridization (CGH) of peripheral blood mononuclear cells (PBMC) and tumor tissue cell in non-small cell lung cancer (NSCLC)

Lung tumor cell DNA copy number alteration (CNA) was expected to display specific patterns such as a large-scale amplification or deletion of chromosomal arms, as previously published data have reported. Peripheral blood mononuclear cell (PBMC) CNA however, was expected to show normal variations in cancer patients as well as healthy individuals, and has thus been used as normal control DNA samples in various published studies. We performed array CGH to measure and compare genetic changes in terms of the CNA of PBMC samples as well as DNA isolated from tumor tissue samples, obtained from 24 non-small cell lung cancer patients. Contradictory to expectations, our studies showed that the PBMC CNA also showed chromosomal variant regions. The list included well-known tumor-associated NTRK1, FGF8, TP53, and TGFβ1 genes and potentially novel oncogenes such as THPO (3q27.1), JMJD1B, and EGR1 (5q31.2), which was investigated in this study. The results of this study highlighted the connection between PBMC and tumor cell genomic DNA in lung cancer patients. However, the application of these studies to cancer prognosis may pose a challenge due to the large amount of information contained in genetic predisposition and family history that has to be processed for useful downstream clinical applications.     

DNA copy number changes in human malignant fibrous histiocytomas by array comparative genomic hybridisation

Background

Malignant fibrous histiocytomas (MFHs), or undifferentiated pleomorphic sarcomas, are in general high-grade tumours with extensive chromosomal aberrations. In order to identify recurrent chromosomal regions of gain and loss, as well as novel gene targets of potential importance for MFH development and/or progression, we have analysed DNA copy number changes in 33 MFHs using microarray-based comparative genomic hybridisation (array CGH).

Principal findings

In general, the tumours showed numerous gains and losses of large chromosomal regions. The most frequent minimal recurrent regions of gain were 1p33-p32.3, 1p31.3-p31.2 and 1p21.3 (all gained in 58% of the samples), as well as 1q21.2-q21.3 and 20q13.2 (both 55%). The most frequent minimal recurrent regions of loss were 10q25.3-q26.11, 13q13.3-q14.2 and 13q14.3-q21.1 (all lost in 64% of the samples), as well as 2q36.3-q37.2 (61%), 1q41 (55%) and 16q12.1-q12.2 (52%). Statistical analyses revealed that gain of 1p33-p32.3 and 1p21.3 was significantly associated with better patient survival (P = 0.021 and 0.046, respectively). Comparison with similar array CGH data from 44 leiomyosarcomas identified seven chromosomal regions; 1p36.32-p35.2, 1p21.3-p21.1, 1q32.1-q42.13, 2q14.1-q22.2, 4q33-q34.3, 6p25.1-p21.32 and 7p22.3-p13, which were significantly different in copy number between the MFHs and leiomyosarcomas.

Conclusions

A number of recurrent regions of gain and loss have been identified, some of which were associated with better patient survival. Several specific chromosomal regions with significant differences in copy number between MFHs and leiomyosarcomas were identified, and these aberrations may be used as additional tools for the differential diagnosis of MFHs and leiomyosarcomas.     

CEQer: A Graphical Tool for Copy Number and Allelic Imbalance Detection from Whole-Exome Sequencing Data

Copy number alterations (CNA) are common events occurring in leukaemias and solid tumors. Comparative Genome Hybridization (CGH) is actually the gold standard technique to analyze CNAs; however, CGH analysis requires dedicated instruments and is able to perform only low resolution Loss of Heterozygosity (LOH) analyses. Here we present CEQer (Comparative Exome Quantification analyzer), a new graphical, event-driven tool for CNA/allelic-imbalance (AI) coupled analysis of exome sequencing data. By using case-control matched exome data, CEQer performs a comparative digital exonic quantification to generate CNA data and couples this information with exome-wide LOH and allelic imbalance detection. This data is used to build mixed statistical/heuristic models allowing the identification of CNA/AI events. To test our tool, we initially used in silico generated data, then we performed whole-exome sequencing from 20 leukemic specimens and corresponding matched controls and we analyzed the results using CEQer. Taken globally, these analyses showed that the combined use of comparative digital exon quantification and LOH/AI allows generating very accurate CNA data. Therefore, we propose CEQer as an efficient, robust and user-friendly graphical tool for the identification of CNA/AI in the context of whole-exome sequencing data.     

Genomic changes and HPV type in cervical carcinoma

To identify chromosomal regions that may include the loci of abnormally expressed cellular genes and may be specifically altered depending on the histological subtype of the tumor, we studied primary cervical carcinoma using CGH and HPV genotyping. Eighty-seven percent of the primary tumors were positive for DNA of a "high-risk" HPV type (e.g., 16 or 18). In the cervical carcinomas, without reference to histologic subtype, overrepresentation of chromosome 3q was the most consistent chromosomal aberration with underrepresentation of chromosome 3p also a frequent finding. Chromosome arms 1q, 5p, 20q, and Xq were overrepresented in many tumors and 3p loss and 5p, 8q, and 16q gain were only associated with squamous cell carcinoma in this series.