Integrated analysis of copy number variation and genome-wide expression profiling in colorectal cancer tissues

Decoupling the dynamic mechanism revealed by FGFR2 mutation-induced population shift

The fibroblast growth factor receptor 2 (FGFR2) is a key component in cellular signaling networks, and its dysfunctional activation has been implicated in various diseases including cancer and developmental disorders. Mutations at the activation loop (A-loop) have been suggested to trigger an increased basal kinase activity. However, the molecular mechanism underlying this highly dynamic process has not been fully understood due to the limitation of static structural information. Here, we conducted multiple, large-scale Gaussian accelerated molecular dynamics simulations of five (K659E, K659N, K659M, K659Q, and K659T) FGFR2 mutants at the A-loop, and comprehensively analyzed the dynamic molecular basis of FGFR2 activation. The results quantified the population shift of each system, revealing that all mutants had a higher proportion of active-like states. Using Markov state models, we extracted the representative structure of different conformational states and identified key residues related to the increased kinase activity. Furthermore, community network analysis showed enhanced information connections in the mutants, highlighting the long-range allosteric communication between the A-loop and the hinge region. Our findings may provide insights into the dynamic mechanism for FGFR2 dysfunctional activation and allosteric drug discovery.Communicated by Ramaswamy H. Sarma.

The chromatin remodeler CHD6 promotes colorectal cancer development by regulating TMEM65-mediated mitochondrial dynamics via EGF and Wnt signaling

Chromodomain helicase DNA binding protein (CHD) family plays critical roles in regulating gene transcription. The family is linked to cancer disease, but the family member's role in tumorigenesis remains largely unknown. Here, we report that CHD6 is highly expressed in colorectal cancer (CRC). CHD6 knockdown inhibited cancer cell proliferation, migration, invasion, and tumorigenesis. Consistently, Villin-specific Chd6 knockout in mice attenuates cancer formation in AOM/DSS model. We found that aberrant EGF signals promoted the stability of CHD6 by diminishing ubiquitin-mediated degradation. EGF signal inhibits GSK3β activity, which in turn prevents phosphodegron formation of CHD6, thereby hindering E3 ligase FBXW7-mediated CHD6 ubiquitination and degradation. CHD6's chromatin remodeler activity engages in binding Wnt signaling transcription factor TCF4 to facilitate the transcriptional expression of TMEM65, a mitochondrial inner membrane protein involved in ATP production and mitochondrial dynamics. In addition, Wnt signaling is also an upstream regulator of CHD6. CHD6 promoter contains TCF4 and β-catenin binding site, and CHD6 can be transcriptionally activated by Wnt ligand to facilitate TMEM65 transcription. Thus CHD6-TMEM65 axis can be regulated by both EGF and Wnt signaling pathways through two different mechanisms. We further illustrate that CHD6-TMEM65 axis is deregulated in cancer and that co-administration of Wnt inhibitor LGK974 and the anti-EGFR monoclonal antibody cetuximab largely restricted the growth of patient-derived xenografts of CRC. Targeting CHD6-TMEM65 axis may be effective for cancer intervention.

CHD6 promotes broad nucleosome eviction for transcriptional activation in prostate cancer cells

Despite being a member of the chromodomain helicase DNA-binding protein family, little is known about the exact role of CHD6 in chromatin remodeling or cancer disease. Here we show that CHD6 binds to chromatin to promote broad nucleosome eviction for transcriptional activation of many cancer pathways. By integrating multiple patient cohorts for bioinformatics analysis of over a thousand prostate cancer datasets, we found CHD6 expression elevated in prostate cancer and associated with poor prognosis. Further comprehensive experiments demonstrated that CHD6 regulates oncogenicity of prostate cancer cells and tumor development in a murine xenograft model. ChIP-Seq for CHD6, along with MNase-Seq and RNA-Seq, revealed that CHD6 binds on chromatin to evict nucleosomes from promoters and gene bodies for transcriptional activation of oncogenic pathways. These results demonstrated a key function of CHD6 in evicting nucleosomes from chromatin for transcriptional activation of prostate cancer pathways.

Mechanistic Insights into the Protection Effect of Argonaute-RNA Complex on the HCV Genome

While host miRNA usually plays an antiviral role, the relentless tides of viral evolution have carved out a mechanism to recruit host miRNA as a viral protector. By complementing miR-122 at the 5' end of the genome, the hepatitis C virus (HCV) gene can form a complex with Argonaute 2 (Ago2) protein to protect the 5' end of HCV RNA from exonucleolytic attacks. Experiments showed that the disruption of the stem-loop 1(SL1) structure and the 9th nucleotide (T9) of HCV site 1 RNA could enhance the affinity of the Ago2 protein to the HCV site 1 RNA (target RNA). However, the underlying mechanism of how the conformation and dynamics of the Ago2: miRNA: target RNA complex is affected by the SL1 and T9 remains unclear. To address this, we performed large-scale molecular dynamics simulations on the AGO2-miRNA complex binding with the WT target, T9-abasic target and SL1-disruption target, respectively. The results revealed that the T9 and SL1 structures could induce the departing motion of the PAZ, PIWI and N domains, propping up the mouth of the central groove which accommodates the target RNA, causing the instability of the target RNA and disrupting the Ago2 binding. The coordinated motion among the PAZ, PIWI and N domains were also weakened by the T9 and SL1 structures. Moreover, we proposed a new model wherein the Ago2 protein could adopt a more constraint conformation with the proximity and more correlated motions of the PAZ, N and PIWI domains to protect the target RNA from dissociation. These findings reveal the mechanism of the Ago2-miRNA complex's protective effect on the HCV genome at the atomic level, which will offer guidance for the design of drugs to confront the protection effect and engineering of Ago2 as a gene-regulation tool.

CSMD1 Shows Complex Patterns of Somatic Copy Number Alterations and Expressions of mRNAs and Target Micro RNAs in Esophageal Squamous Cell Carcinoma

Simple Summary

Human Cub and Sushi Multiple Domains 1 (CSMD1) is a novel candidate tumor-suppressor gene. We investigated CSMD1 in esophageal squamous cell carcinoma (ESCC) by performing an integrated analysis of somatic DNA alterations (i.e., copy number alteration, allelic imbalance, and loss of heterozygosity) with RNA expressions (mRNA and target miRNAs) on specimens from the same ESCC patients, using data from SNP, miRNA, and RT-PCR arrays. Our results indicate that the CSMD1 gene may play a role in the development of ESCC through complex patterns involving somatic alterations and mRNA expression. Furthermore, somatic copy number alterations in SNPs located in non-coding regions of CSMD1 appear to influence expression of both this gene and its target miRNAs.

Abstract

Background: Human Cub and Sushi Multiple Domains 1 (CSMD1) is a novel candidate tumor-suppressor gene that codes for multiple domains, including complement regulatory and adhesion proteins, and has recently been shown to have alterations in multiple cancers. We investigated CSMD1 in esophageal squamous cell carcinoma (ESCC) by performing an integrated analysis on somatic copy number alterations (CNAs), including copy-number gain or loss, allelic imbalance (AI), loss of heterozygosity (LOH), and the expressions of mRNA and its target miRNAs on specimens from the same patients with ESCC. Results: (i) Two-thirds of ESCC patients had all three types of alterations studied—somatic DNA alterations in 70%, and abnormal expressions of CSMD1 RNA in 69% and in target miRNAs in 66%; patterns among these alterations were complex. (ii) In total, 97% of 888 CSMD1 SNPs studied showed somatic DNA alterations, with most located near exons 4–11, 24–25, 39–40, 55–56, and 69–70. (iii) In total, 68% of SNPs with a CNA were correlated with expression of CSMD1. (iv) A total of 33 correlations between non-coding SNPs and expression of CSMD1 target miRs were found. Conclusions: Our results indicate that the CSMD1 gene may play a role in ESCC through complex patterns of DNA alterations and RNA and miRNA expressions. Alterations in some somatic SNPs in non-coding regions of CSMD1 appear to influence expression of this gene and its target miRNAs.

Mechanistic insights into the clinical Y96D mutation with acquired resistance to AMG510 in the KRASG12C

Special oncogenic mutations in the RAS proteins lead to the aberrant activation of RAS and its downstream signaling pathways. AMG510, the first approval drug for KRAS, covalently binds to the mutated cysteine 12 of KRASG12C protein and has shown promising antitumor activity in clinical trials. Recent studies have reported that the clinically acquired Y96D mutation could severely affect the effectiveness of AMG510. However, the underlying mechanism of the drug-resistance remains unclear. To address this, we performed multiple microsecond molecular dynamics simulations on the KRASG12C−AMG510 and KRASG12C/Y96D−AMG510 complexes at the atomic level. The direct interaction between the residue 96 and AMG510 was impaired owing to the Y96D mutation. Moreover, the mutation yielded higher flexibility and more coupled motion of the switch II and α3-helix, which led to the departing motion of the switch II and α3-helix. The resulting departing motion impaired the interaction between the switch II and α3-helix and subsequently induced the opening and loosening of the AMG510 binding pocket, which further disrupted the interaction between the key residues in the pocket and AMG510 and induced an increased solvent exposure of AMG510. These findings reveal the resistance mechanism of AMG510 to KRASG12C/Y96D, which will help to offer guidance for the development of KRAS targeted drugs to overcome acquired resistance.

Oncogenic Role of miR-200c-3p in High-Grade Serous Ovarian Cancer Progression via Targeting the 3'-Untranslated Region of DLC1

High-grade serous ovarian cancer (HGSC) is the most common ovarian cancer with highly metastatic properties. A small non-coding RNA, microRNA (miRNA) was discovered to be a major regulator in many types of cancers through binding at the 3′-untranslated region (3′UTR), leading to degradation of the mRNA. In this study, we sought to investigate the underlying mechanisms involved in the dysregulation of miR-200c-3p in HGSC progression and metastasis. We identified the upregulation of miR-200c-3p expression in different stages of HGSC clinical samples and the downregulation of the tumor suppressor gene, Deleted in Liver Cancer 1 (DLC1), expression. Over expression of miR-200c-3p in HGSC cell lines downregulated DLC1 but upregulated the epithelial marker, E-cadherin (CDH1). Based on in silico analysis, two putative binding sites were found within the 3′UTR of DLC1, and we confirmed the direct binding of miR-200c-3p to the target binding motif at position 1488–1495 bp of 3′UTR of DLC1 by luciferase reporter assay in a SKOV3 cell line co-transfected with vectors and miR-200c-3p mimic. These data showed that miR-200c-3p regulated the progression of HGSC by regulating DLC1 expression post-transcription and can be considered as a promising target for therapeutic purposes.

Tribbles Pseudokinase 3 Regulation and Contribution to Cancer

Simple Summary

Accumulating evidence supports a key function for Tribbles proteins in oncogenesis, both in leukemia and solid tumors. However, the exact role of these proteins is hard to define since in a context-dependent manner they can function as both oncogenes and tumor suppressors. Their complex role arises from the capacity to interact with a wide range of target molecules thereby acting as molecular scaffolds and signaling regulators of multiple pathways. This review focuses on one particular Tribbles family member, namely, TRIB3, addressing its gene and protein expression, as well as its role in cancer development and progression.

Abstract

The first Tribbles protein was identified as critical for the coordination of morphogenesis in Drosophila melanogaster. Three mammalian homologs were subsequently identified, with a structure similar to classic serine/threonine kinases, but lacking crucial amino acids for the catalytic activity. Thereby, the very weak ATP affinity classifies TRIB proteins as pseudokinases. In this review, we provide an overview of the regulation of TRIB3 gene expression at both transcriptional and post-translational levels. Despite the absence of kinase activity, TRIB3 interferes with a broad range of cellular processes through protein–protein interactions. In fact, TRIB3 acts as an adaptor/scaffold protein for many other proteins such as kinase-dependent proteins, transcription factors, ubiquitin ligases, or even components of the spliceosome machinery. We then state the contribution of TRIB3 to cancer development, progression, and metastasis. TRIB3 dysregulation can be associated with good or bad prognosis. Indeed, as TRIB3 interacts with and regulates the activity of many key signaling components, it can act as a tumor-suppressor or oncogene in a context-dependent manner.

Reactivation of oncogenes involved in G1/S transcription and apoptosis pathways by low dose decitabine promotes HT29 human colon cancer cell growth in vitro

BACKGROUND

To examine the effects of low-dose decitabine (DAC) on the proliferation of HT-29 cell lines, and to explore the central mechanism by which low-dose DAC affects HT-29 cell proliferation using a systematic biological approach.

METHODS

First, we examined the global effects of DAC on cell proliferation, the cell cycle, and apoptosis in HT29 colon cancer cells. Then, a series test of cluster (STC) analysis and weighted gene coexpression network analysis (WGCNA) were employed to identify critical pathways involved in the response to DAC treatment using 3 datasets from the GEO database. Finally, the expression changes and promoter methylation levels of hub genes were further confirmed by in vitro experiments.

RESULTS

Low-dose DAC (less than 1 µM) promoted the proliferation and colony formation ability of HT-29 cell lines. The results of the system-level analysis, including STC analysis, WGCNA, and Gene set variation analysis (GSVA), showed that DAC modulated 3 critical pathways: G1/S-specific transcription involved in E2F-mediated regulation of Cyclin E-associated events, apoptosis pathways, and EMT pathways. Subsequent in vitro experiments showed that low-dose DAC (0.1 µM) promoted G1/S-specific transcription and decreased apoptosis rates. Then, several regulatory hub oncogenes in these 3 pathways, CCNE1, E2F1, BCL2, PCNA, FOXC1, VIM, CXCL1, and VCAM1, were further confirmed to be activated by DAC at either the mRNA or protein level. We chose the oncogene BCL2 as an example and detected its methylation status and the effect of low-dose DAC on BCL2 expression. Data from TCGA and Oncomine databases demonstrated that BCL2 was decreased in colon cancer compared with normal mucosa. Further analysis showed that BCL2 had an increased degree of promoter methylation in 12 methylated sites in colon cancer compared with normal colon tissues. Bisulfite sequencing PCR showed that low-dose DAC decreased the methylation rate at the BCL2 promoter region.

CONCLUSIONS

We concluded that low-dose DAC treatment resulted in a cancer-promoting effect in HT29 cell lines. Mechanistically, high methylation levels at the promoter region of oncogenes with dominant effects in CRC, such as BCL2 in HT29, might play a role in suppressing CRC by inhibiting oncogene expression. Low-dose DAC treatment triggered BCL2 expression by decreasing its promoter methylation level, thereby resulting in cancer promotion.

Acquired resistance to DZNep-mediated apoptosis is associated with copy number gains of AHCY in a B-cell lymphoma model

BACKGROUND

Enhancer of zeste homolog 2 (EZH2) is considered an important driver of tumor development and progression by its histone modifying capabilities. Inhibition of EZH2 activity is thought to be a potent treatment option for eligible cancer patients with an aberrant EZH2 expression profile, thus the indirect EZH2 inhibitor 3-Deazaneplanocin A (DZNep) is currently under evaluation for its clinical utility. Although DZNep blocks proliferation and induces apoptosis in different tumor types including lymphomas, acquired resistance to DZNep may limit its clinical application.

METHODS

To investigate possible mechanisms of acquired DZNep resistance in B-cell lymphomas, we generated a DZNep-resistant clone from a previously DZNep-sensitive B-cell lymphoma cell line by long-term treatment with increasing concentrations of DZNep (ranging from 200 to 2000 nM) and compared the molecular profiles of resistant and wild-type clones. This comparison was done using molecular techniques such as flow cytometry, copy number variation assay (OncoScan and TaqMan assays), fluorescence in situ hybridization, Western blot, immunohistochemistry and metabolomics analysis.

RESULTS

Whole exome sequencing did not indicate the acquisition of biologically meaningful single nucleotide variants. Analysis of copy number alterations, however, demonstrated among other acquired imbalances an amplification (about 30 times) of the S-adenosyl-L-homocysteine hydrolase (AHCY) gene in the resistant clone. AHCY is a direct target of DZNep and is critically involved in the biological methylation process, where it catalyzes the reversible hydrolysis of S-adenosyl-L-homocysteine to L-homocysteine and adenosine. The amplification of the AHCY gene is paralleled by strong overexpression of AHCY at both the transcriptional and protein level, and persists upon culturing the resistant clone in a DZNep-free medium.

CONCLUSIONS

This study reveals one possible molecular mechanism how B-cell lymphomas can acquire resistance to DZNep, and proposes AHCY as a potential biomarker for investigation during the administration of EZH2-targeted therapy with DZNep.

Multi-omics Analysis of Primary Cell Culture Models Reveals Genetic and Epigenetic Basis of Intratumoral Phenotypic Diversity

Uncovering the functionally essential variations related to tumorigenesis and tumor progression from cancer genomics data is still challenging due to the genetic diversity among patients, and extensive inter- and intra-tumoral heterogeneity at different levels of gene expression regulation, including but not limited to the genomic, epigenomic, and transcriptional levels. To minimize the impact of germline genetic heterogeneities, in this study, we establish multiple primary cultures from the primary and recurrent tumors of a single patient with hepatocellular carcinoma (HCC). Multi-omics sequencing was performed for these cultures that encompass the diversity of tumor cells from the same patient. Variations in the genome sequence, epigenetic modification, and gene expression are used to infer the phylogenetic relationships of these cell cultures. We find the discrepancy among the relationships revealed by single nucleotide variations (SNVs) and transcriptional/epigenomic profiles from the cell cultures. We fail to find overlap between sample-specific mutated genes and differentially expressed genes (DEGs), suggesting that most of the heterogeneous SNVs among tumor stages or lineages of the patient are functionally insignificant. Moreover, copy number alterations (CNAs) and DNA methylation variation within gene bodies, rather than promoters, are significantly correlated with gene expression variability among these cell cultures. Pathway analysis of CNA/DNA methylation-related genes indicates that a single cell clone from the recurrent tumor exhibits distinct cellular characteristics and tumorigenicity, and such an observation is further confirmed by cellular experiments both in vitro and in vivo. Our systematic analysis reveals that CNAs and epigenomic changes, rather than SNVs, are more likely to contribute to the phenotypic diversity among subpopulations in the tumor. These findings suggest that new therapeutic strategies targeting gene dosage and epigenetic modification should be considered in personalized cancer medicine. This culture model may be applied to the further identification of plausible determinants of cancer metastasis and relapse.

Deregulation of CSMD1 targeted by microRNA-10b drives gastric cancer progression through the NF-κB pathway

Aim: This study aimed to investigate the oncogenic activity of microRNA-10b by targeting CUB and sushi multiple domains protein 1 (CSMD1) in human gastric cancer (GC) and the underlying mechanisms. Methods: The expression of CSMD1 in human GC tissues was evaluated by real-time reverse transcription polymerase chain reaction (RT-PCR), immunoblotting, and immunohistochemical analysis. The expressive abundance of microRNA-10b was detected by stem-loop RT-PCR. Molecular and cellular techniques, including lentiviral vector-mediated knockdown or overexpression, were used to elucidate the effect of microRNA-10b on the expression of CSMD1. Results: CSMD1 was targeted and downregulated by microRNA-10b in human GC tissues and cells, and the down-regulated expression of CSMD1 contributed to poor survival. The knockdown of microRNA-10b expression inhibited cell proliferation in GC cells in vitro and tumor growth in vivo. The inhibition of microRNA-10b expression repressed invasion and migration of HGC27 cells and retarded GC cells metastasis to the liver in Balb/c nude mice. The up-regulated expression of microRNA-10b promoted the proliferation and metastasis of MKN74 cell in vitro. Intratumoral injection of microRNA-10b mimic also promoted the growth and metastasis of tumor xenografts in Balb/c nude mice. Mechanistically, microRNA-10b promoted the invasion and metastasis of human GC cells through inhibiting the expression of CSMD1, leading to the activation of the nuclear factor-κB (NF-κB) pathway that links inflammation to carcinogenesis, subsequently resulting in the upregulation of c-Myc, cyclin D1 (CCND1), and epithelial-mesenchymal transition (EMT) markers. Conclusions: The findings established that microRNA-10b is an oncomiR that drives metastasis. Moreover, a set of critical tumor suppressor mechanisms was defined that microRNA-10b overcame to drive human GC progression.

The landscape of genomic copy number alterations in colorectal cancer and their consequences on gene expression levels and disease outcome

Aneuploidy, the unbalanced state of the chromosome content, represents a hallmark of most solid tumors, including colorectal cancer. Such aneuploidies result in tumor specific genomic imbalances, which emerge in premalignant precursor lesions. Moreover, increasing levels of chromosomal instability have been observed in adenocarcinomas and are maintained in distant metastases. A number of studies have systematically integrated copy number alterations with gene expression changes in primary carcinomas, cell lines, and experimental models of aneuploidy. In fact, chromosomal aneuploidies target a number of genes conferring a selective advantage for the metabolism of the cancer cell. Copy number alterations not only have a positive correlation with expression changes of the majority of genes on the altered genomic segment, but also have effects on the transcriptional levels of genes genome-wide. Finally, copy number alterations have been associated with disease outcome; nevertheless, the translational applicability in clinical practice requires further studies. Here, we (i) review the spectrum of genetic alterations that lead to colorectal cancer, (ii) describe the most frequent copy number alterations at different stages of colorectal carcinogenesis, (iii) exemplify their positive correlation with gene expression levels, and (iv) discuss copy number alterations that are potentially involved in disease outcome of individual patients.

Genetic Aberration Analysis in Thai Colorectal Adenoma and Early-Stage Adenocarcinoma Patients by Whole-Exome Sequencing

Colorectal adenomas are precursor lesions of colorectal adenocarcinoma. The transition from adenoma to carcinoma in patients with colorectal cancer (CRC) has been associated with an accumulation of genetic aberrations. However, criteria that can screen adenoma progression to adenocarcinoma are still lacking. This present study is the first attempt to identify genetic aberrations, such as the somatic mutations, copy number variations (CNVs), and high-frequency mutated genes, found in Thai patients. In this study, we identified the genomic abnormality of two sample groups. In the first group, five cases matched normal-colorectal adenoma-colorectal adenocarcinoma. In the second group, six cases matched normal-colorectal adenomas. For both groups, whole-exome sequencing was performed. We compared the genetic aberration of the two sample groups. In both normal tissues compared with colorectal adenoma and colorectal adenocarcinoma analyses, somatic mutations were observed in the tumor suppressor gene APC (Adenomatous polyposis coli) in eight out of ten patients. In the group of normal tissue comparison with colorectal adenoma tissue, somatic mutations were also detected in Catenin Beta 1 (CTNNB1), Family With Sequence Similarity 123B (FAM123B), F-Box And WD Repeat Domain Containing 7 (FBXW7), Sex-Determining Region Y-Box 9 (SOX9), Low-Density Lipoprotein Receptor-Related Protein 5 (LRP5), Frizzled Class Receptor 10 (FZD10), and AT-Rich Interaction Domain 1A (ARID1A) genes, which are involved in the Wingless-related integration site (Wnt) signaling pathway. In the normal tissue comparison with colorectal adenocarcinoma tissue, Kirsten retrovirus-associated DNA sequences (KRAS), Tumor Protein 53 (TP53), and Ataxia-Telangiectasia Mutated (ATM) genes are found in the receptor tyrosine kinase-RAS (RTK-RAS) signaling pathway and p53 signaling pathway, respectively. These results suggest that APC and TP53 may act as a potential screening marker for colorectal adenoma and early-stage CRC. This preliminary study may help identify patients with adenoma and early-stage CRC and may aid in establishing prevention and surveillance strategies to reduce the incidence of CRC.

PLAGL2 and POFUT1 are regulated by an evolutionarily conserved bidirectional promoter and are collaboratively involved in colorectal cancer by maintaining stemness

Background

Our previous study revealed that PLAGL2 or POFUT1 can promote tumorigenesis and maintain significant positive correlations in colorectal cancer (CRC). However, the mechanism leading to the co-expression and the underlying functional and biological implications remain unclear.

Methods

Clinical tumor tissues and TCGA dataset were utilized to analyze the co-expression of PLAGL2 and POFUT1. Luciferase reporter assays, specially made bidirectional promoter vectors and ectopic expression of 3’UTR were employed to study the mechanisms of co-expression. In vitro and in vivo assays were performed to further confirm the oncogenic function of both. The sphere formation assay, immunofluorescence, Western blot and qRT-PCR were performed to investigate the effect of both genes in colorectal cancer stem cells (CSCs).

Findings

PLAGL2 and POFUT1 maintained co-expression in CRC (r = 0.91, p < .0001). An evolutionarily conserved bidirectional promoter, rather than post-transcriptional regulation by competing endogenous RNAs, caused the co-expression of PLAGL2 and POFUT1 in CRC. The bidirectional gene pair PLAGL2/POFUT1 was subverted in CRC and acted synergistically to promote colorectal tumorigenesis by maintaining stemness of colorectal cancer stem cells through the Wnt and Notch pathways. Finally, PLAGL2 and POFUT1 share transcription factor binding sites, and introducing mutations into promoter regions with shared transcription regulatory elements led to a decrease in the PLAGL2/POFUT1 promoter activity in both directions.

Interpretation

Our team identified for the first time a bidirectional promoter pair oncogene, PLAGL2-POFUT1, in CRC. The two genes synergistically promote the progression of CRC and affect the characteristics of CSCs, which can offer promising intervention targets for clinicians and researchers.

Fund

National Nature Science Foundation of China, the Hunan province projects of Postgraduate Independent Exploration and Innovation of Central South University.

aCGH Analysis of Predictive Biomarkers for Response to Bevacizumab plus Oxaliplatin- or Irinotecan-Based Chemotherapy in Patients with Metastatic Colorectal Cancer

BACKGROUND

The randomized phase III study (WJOG4407G) showed equivalent efficacy between FOLFOX and FOLFIRI in combination with bevacizumab as the first-line treatment for metastatic colorectal cancer (mCRC). We studied whole genome copy number profiles using array-based comparative genomic hybridization (aCGH) analysis of tumor tissue samples obtained in this study. The aim of this study was to identify gene copy number alterations that could aid in selecting either FOLFOX or FOLFIRI in combination with bevacizumab for patients with mCRC.

MATERIALS AND METHODS

DNA was purified from 154 pretreatment formalin-fixed paraffin-embedded tissue samples (75 from the FOLFOX arm and 79 from the FOLFIRI arm) of 395 patients enrolled in the WJOG4407G trial and analyzed by aCGH. Genomic regions greater than 1.2-fold were regarded as copy number gain (CNG).

RESULTS

Patient characteristics between the treatment arms were well balanced except for tumor laterality (left side; 64% in FOLFOX arm and 80% in FOLFIRI arm, p = .07). FOLFIRI showed a trend toward better response rate (RR), progression-free survival (PFS) and overall survival (OS) than FOLFOX in the patients with CNG of chromosome 8q24.1 (Fisher's exact test, p = .134 for RR; interaction test, p = .102 for PFS and p = .003 for OS) and 8q24.2 (Fisher's exact test, p = .179 for RR; interaction test, p = .144 for PFS and p = .002 for OS).

CONCLUSION

Chromosome 8q24.1-q24.2 may contain genes that could potentially serve as predictive markers for selecting either FOLFOX or FOLFIRI in combination with bevacizumab for treatment of patients with mCRC.

IMPLICATIONS FOR PRACTICE

Bevacizumab has been used as a standard first-line treatment for patients with metastatic colorectal cancer (mCRC) in combination with either oxaliplatin-based or irinotecan-based chemotherapy. Until now, there has been no predictive marker to choose between the two combination chemotherapies. This array-based comparative genomic hybridization analysis revealed that the difference in therapeutic effect between the two combination chemotherapies is prominent in patients with mCRC with gene copy number gain in chromosome 8p24.1-p24.2. Such patients showed more favorable response and survival when treated with irinotecan-based combination chemotherapy. Overlapping genes commonly found in this region may be predictive biomarkers of the efficacy of the combination chemotherapy with bevacizumab.

Positive Caricature Transcriptomic Effects Associated with Broad Genomic Aberrations in Colorectal Cancer

We re-examined the correlation between Broad Genomic Aberrations (BGAs) and transcriptomic profiles in Colorectal Cancer (CRC). Two types of BGAs have been examined: Broad Copy-Number Abnormal regions (BCNAs), distinguished in gain- and loss-type, and Copy-Neutral Loss of Heterozygosities (CNLOHs). Transcripts are classified as “OverT” or “UnderT” if overexpressed or underexpressed comparing CRCs bearing a specific BGA to CRCs not bearing it and as “UpT” or “DownT” if upregulated or downregulated in cancer compared to normal tissue. BGA-associated effects were evaluated by changes in the “Chromosomal Distribution Index” (CDI) of different transcript classes. Data show that UpT are more sensitive than DownT to BCNA-associated gene dosage effects. “Over-UpT” genes are upregulated in cancer and further overexpressed by gene dosage, defining the so called “positive caricature transcriptomic effect”. When Over-UpT genes are ranked according to overexpression, top positions are occupied by genes implicated at the functional and therapeutic level in CRC. We show that cancer-upregulated transcripts are sensitive markers of BCNA-induced effects and suggest that analysis of positive caricature transcriptomic effects can provide clues toward the identification of BCNA-associated cancer driver genes.

Comprehensive bioinformatics analysis of the characterization and determination underlying mechanisms of over-expression and co-expression of genes residing on 20q in colorectal cancer

The Long arm of chromosome 20 (20q) is closely related to the development of colorectal cancer, so identifying the expression profile of genes on 20q through a comprehensive overview is indispensable. In this article, preliminar experimental data, several available databases and bioinformatics tools such as the Cancer Genome Atlas, the Encyclopedia of DNA Elements, the JASPAR database and starBase were combined to analyze the correlation between genes and chromosomal aberrations, microRNA and transcription factors, as well as to explore the expression feature and potential regulative mechanism. The results showed that the most frequently unregulated genes in colorectal cancer arelocated on chromosome 20q, present a significant CNA–mRNA correlation.Furthermore, the genes with mRNA overexpression showed co-expression features and tended to be clustered within the same genomic neighborhoods. Then, several genes were selected to carry out further analysis and demonstrated that shared transcription factors, a conserved bidirectional promoter, and competition for a limited pool of microRNAin the 3’UTR of mRNA may be the underlying mechanisms behind the co-expression of physically adjacent genes.Finally, the databases, Lentivirus shRNA, and qPCR were used to find that these adjacent genes with co-expression cooperatively participated in the same biological pathways associated with the pathogenesis and development of colorectal cancer.

Comprehensive analysis of copy number aberrations in microsatellite stable colon cancer in view of stromal component

BACKGROUND

Somatic copy number aberrations (CNAs) are common acquired changes in cancer cells having an important role in the progression of colon cancer (colorectal cancer, CRC). This study aimed to perform a characterisation of CNA and their impact in gene expression.

METHODS

Copy number aberrations were inferred from SNP array data in a series of 99 CRC. Copy number aberration events were calculated and used to assess the association between copy number dosage, clinical and molecular characteristics of the tumours, and gene expression changes. All analyses were adjusted for the quantity of stroma in each sample, which was inferred from gene expression data.

RESULTS

High heterogeneity among samples was observed; the proportion of altered genome ranged between 0.04 and 26.6%. Recurrent CNA regions with gains were frequent in chromosomes 7p, 8q, 13q, and 20, whereas 8p, 17p, and 18 cumulated losses. A significant positive correlation was observed between the number of somatic mutations and total CNA (Spearman's r=0.42, P=0.006). Approximately 37% of genes located in CNA regions changed their level of expression and the average partial correlation (adjusted for stromal content) with copy number was 0.54 (interquartile range 0.20 to 0.81). Altered genes showed enrichment in pathways relevant for CRC. Tumours classified as CMS2 and CMS4 by the consensus molecular subtyping showed higher frequency of CNA. Losses of one small region in 1p36.33, with gene CDK11B, were associated with poor prognosis. More than 66% of the recurrent CNA were validated in the The Cancer Genome Atlas (TCGA) data when analysed with the same procedure. Furthermore, 79% of the genes with altered expression in our data were validated in the TCGA.

CONCLUSIONS

Although CNA are frequent events in microsatellite stable CRC, few focal recurrent regions were found. These aberrations have strong effects on gene expression and contribute to deregulate relevant cancer pathways. Owing to the diploid nature of stromal cells, it is important to consider the purity of tumour samples to accurately calculate CNA events in CRC.

The Chromodomain Helicase DNA-Binding Chromatin Remodelers: Family Traits that Protect from and Promote Cancer

A plethora of mutations in chromatin regulators in diverse human cancers is emerging, attesting to the pivotal role of chromatin dynamics in tumorigenesis. A recurrent theme is inactivation of the chromodomain helicase DNA-binding (CHD) family of proteins-ATP-dependent chromatin remodelers that govern the cellular machinery's access to DNA, thereby controlling fundamental processes, including transcription, proliferation, and DNA damage repair. This review highlights what is currently known about how genetic and epigenetic perturbation of CHD proteins and the pathways that they regulate set the stage for cancer, providing new insight for designing more effective anti-cancer therapies.

DNA alterations in Cd133+ and Cd133- tumour cells enriched from intra-operative human colon tumour biopsies

Background

Tumour stem cells are considered important to promote disease progression, recurrence and treatment resistance following chemotherapy in colon cancer. However, genomic analyses of colorectal cancer have mainly been performed on integrated tumour tissue consisting of several different cell types in addition to differentiated tumour cells. The purpose of the present study was to compare genomic alterations in two cell fractions enriched of CD133+ and CD133−/EpCAM+ cells, respectively, obtained from fresh intraoperative human tumour biopsies.

Methods

The tumour biopsies were fractionated into CD133+ and CD133−/EpCAM+ cells by immunomagnetic separation, confirmed by immunocytochemistry and Q-PCR. DNA were extracted and used for array comparative genome hybridization (aCGH) after whole genome amplification. Frozen tumour tissue biopsies were used for DNA/RNA extraction and Q-PCR analyses to check for DNA alterations detected in the cell fractions.

Results

The number and size of DNA alterations were equally distributed across the cell fractions; however, large deletions were detected on chromosome 1, 7 and 19 in CD133−/EpCAM+ cells. Deletions were frequent in both cell fractions and a deletion on chromosome 19p was confirmed in 90% of the patients.

Conclusion

Isolation of enriched cells derived from tumour tissue revealed mainly genomic deletions, which were not observed in tumour tissue DNA analyses. CD133+ cells were genetically heterogeneous among patients without any defined profile compared to CD133−/EpCAM+ cells.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-017-3206-8) contains supplementary material, which is available to authorized users.

iGC-an integrated analysis package of gene expression and copy number alteration

BACKGROUND

With the advancement in high-throughput technologies, researchers can simultaneously investigate gene expression and copy number alteration (CNA) data from individual patients at a lower cost. Traditional analysis methods analyze each type of data individually and integrate their results using Venn diagrams. Challenges arise, however, when the results are irreproducible and inconsistent across multiple platforms. To address these issues, one possible approach is to concurrently analyze both gene expression profiling and CNAs in the same individual.

RESULTS

We have developed an open-source R/Bioconductor package (iGC). Multiple input formats are supported and users can define their own criteria for identifying differentially expressed genes driven by CNAs. The analysis of two real microarray datasets demonstrated that the CNA-driven genes identified by the iGC package showed significantly higher Pearson correlation coefficients with their gene expression levels and copy numbers than those genes located in a genomic region with CNA. Compared with the Venn diagram approach, the iGC package showed better performance.

CONCLUSION

The iGC package is effective and useful for identifying CNA-driven genes. By simultaneously considering both comparative genomic and transcriptomic data, it can provide better understanding of biological and medical questions. The iGC package's source code and manual are freely available at https://www.bioconductor.org/packages/release/bioc/html/iGC.html .

Integrated Analysis of Genome-Wide Copy Number Alterations and Gene Expression Profiling of Lung Cancer in Xuanwei, China

Objectives

Lung cancer in Xuanwei (LCXW), China, is known throughout the world for its distinctive characteristics, but little is known about its pathogenesis. The purpose of this study was to screen potential novel “driver genes” in LCXW.

Methods

Genome-wide DNA copy number alterations (CNAs) were detected by array-based comparative genomic hybridization and differentially expressed genes (DEGs) by gene expression microarrays in 8 paired LCXW and non-cancerous lung tissues. Candidate driver genes were screened by integrated analysis of CNAs and DEGs. The candidate genes were further validated by real-time quantitative polymerase chain reaction.

Results

Large numbers of CNAs and DEGs were detected, respectively. Some of the most frequently occurring CNAs included gains at 5p15.33-p15.32, 5p15.1-p14.3, and 5p14.3-p14.2 and losses at 11q24.3, 21q21.1, 21q22.12-q22.13, and 21q22.2. Integrated analysis of CNAs and DEGs identified 24 candidate genes with frequent copy number gains and concordant upregulation, which were considered potential oncogenes, including CREB3L4, TRIP13, and CCNE2. In addition, the analysis identified 19 candidate genes with a negative association between copy number change and expression change, considered potential tumor suppressor genes, including AHRR, NKD2, and KLF10. One of the most studied oncogenes, MYC, may not play a carcinogenic role in LCXW.

Conclusions

This integrated analysis of CNAs and DEGs identified several potential novel LCXW-related genes, laying an important foundation for further research on the pathogenesis of LCXW and identification of novel biomarkers or therapeutic targets.

Chromothripsis and progression-free survival in metastatic colorectal cancer

Metastatic dissemination of the primary tumor is the major cause of death in colorectal cancer (CRC) patients. Multiple chromosomal breaks and chromothripsis, a phenomenon involving multiple chromosomal fragmentations occurring in a single catastrophic event, are associated with cancer genesis, progression and developing of metastases. The aim of this study was to evaluate the effect of chromothripsis and total breakpoint count (breakpoint instability index) on progression-free survival (PFS). A total of 19 patients with metastatic CRC (mCRC) receiving FOLFOX first-line palliative chemotherapy between August, 2011 and October, 2012 were selected for this study. The results indicated that the highest breakpoint count was observed in chromosomes 1, 2 and 6. Chromothripsis was detected in 52.6% of the study patients. Furthermore, chromothripsis was associated with an increased median PFS (mPFS; 14 vs. 8 months, respectively; P=0.03), but an association with overall survival was not identified. The present study demonstrated that chromothripsis affected CRC patient survival, suggesting a role for this event as a prognostic and predictive marker in mCRC treatment.

DNA copy number alterations, gene expression changes and disease-free survival in patients with colorectal cancer: a 10 year follow-up

BACKGROUND

DNA copy number alterations (CNAs) and gene expression changes have amply been encountered in colorectal cancers (CRCs), but the extent at which CNAs affect gene expression, as well as their relevance for tumor development, are still poorly defined. Here we aimed at assessing the clinical relevance of these parameters in a 10 year follow-up study.

METHODS

Tumors and normal adjacent colon mucosa, obtained at primary surgery from 21 CRC patients, were subjected to (i) high-resolution array CGH (a-CGH) for the detection of CNAs and (ii) microarray-based transcriptome profiling for the detection of gene expression (GE) changes. Correlations between these genomic and transcriptomic changes and their associations with clinical and histopathological parameters were assessed with the aim to identify molecular signatures associated with disease-free survival of the CRC patients during a 10 year follow-up.

RESULTS

DNA copy number gains were frequently detected in chromosomes 7, 8q, 13, 19, 20q and X, whereas DNA copy number losses were frequently detected in chromosomes 1p, 4, 8p, 15, 17p, 18, 19 and 22q. None of these alterations were observed in all samples. In addition, we found that 2,498 genes were up- and that 1,094 genes were down-regulated in the tumor samples compared to their corresponding normal mucosa (p < 0.01). The expression of 65 genes was found to be significantly associated with prognosis (p < 0.01). Specifically, we found that up-regulation of the IL17RA, IGF2BP2 and ABCC2 genes, and of genes acting in the mTOR and cytokine receptor pathways, were strongly associated with a poor survival. Subsequent integrated analyses revealed that increased expression levels of the MMP9, BMP7, UBE2C, I-CAM, NOTCH3, NOTCH1, PTGES2, HMGB1 and ERBB3 genes were associated with copy number gains, whereas decreased expression levels of the MUC1, E2F2, HRAS and SIRT3 genes were associated with copy number losses. Pathways related to cell cycle progression, eicosanoid metabolism, and TGF-β and apoptosis signaling, were found to be most significantly affected.

CONCLUSIONS

Our results suggest that CNAs in CRC tumor tissues are associated with concomitant changes in the expression of cancer-related genes. In other genes epigenetic mechanism may be at work. Up-regulation of the IL17RA, IGF2BP2 and ABCC2 genes, and of genes acting in the mTOR and cytokine receptor pathways, appear to be associated with a poor survival. These alterations may, in addition to Dukes' staging, be employed as new prognostic biomarkers for the prediction of clinical outcome in CRC patients.

Novel Somatic Copy Number Alteration Identified for Cervical Cancer in the Mexican American Population

Cervical cancer affects millions of Americans, but the rate for cervical cancer in the Mexican American is approximately twice that for non-Mexican Americans. The etiologies of cervical cancer are still not fully understood. A number of somatic mutations, including several copy number alterations (CNAs), have been identified in the pathogenesis of cervical carcinomas in non-Mexican Americans. Thus, the purpose of this study was to investigate CNAs in association with cervical cancer in the Mexican American population. We conducted a pilot study of genome-wide CNA analysis using 2.5 million markers in four diagnostic groups: reference (n = 125), low grade dysplasia (cervical intraepithelial neoplasia (CIN)-I, n = 4), high grade dysplasia (CIN-II and -III, n = 5) and invasive carcinoma (squamous cell carcinoma (SCC), n = 5) followed by data analyses using Partek. We observed a statistically-significant difference of CNA burden between case and reference groups of different sizes (>100 kb, 10-100 kb and 1-10 kb) of CNAs that included deletions and amplifications, e.g., a statistically-significant difference of >100 kb deletions was observed between the reference (6.6%) and pre-cancer and cancer (91.3%) groups. Recurrent aberrations of 98 CNA regions were also identified in cases only. However, none of the CNAs have an impact on cancer progression. A total of 32 CNA regions identified contained tumor suppressor genes and oncogenes. Moreover, the pathway analysis revealed endometrial cancer and estrogen signaling pathways associated with this cancer (p < 0.05) using Kyoto Encyclopedia of Genes and Genomes (KEGG). This is the first report of CNAs identified for cervical cancer in the U.S. Latino population using high density markers. We are aware of the small sample size in the study. Thus, additional studies with a larger sample are needed to confirm the current findings.

Strategies for Integrated Analysis of Genetic, Epigenetic, and Gene Expression Variation in Cancer: Addressing the Challenges

The development and progression of cancer, a collection of diseases with complex genetic architectures, is facilitated by the interplay of multiple etiological factors. This complexity challenges the traditional single-platform study design and calls for an integrated approach to data analysis. However, integration of heterogeneous measurements of biological variation is a non-trivial exercise due to the diversity of the human genome and the variety of output data formats and genome coverage obtained from the commonly used molecular platforms. This review article will provide an introduction to integration strategies used for analyzing genetic risk factors for cancer. We critically examine the ability of these strategies to handle the complexity of the human genome and also accommodate information about the biological and functional interactions between the elements that have been measured-making the assessment of disease risk against a composite genomic factor possible. The focus of this review is to provide an overview and introduction to the main strategies and to discuss where there is a need for further development.

Influenza Virus and Chromatin: Role of the CHD1 Chromatin Remodeler in the Virus Life Cycle

Influenza A virus requires ongoing cellular transcription to carry out the cap-snatching process. Chromatin remodelers modify chromatin structure to produce an active or inactive conformation, which enables or prevents the recruitment of transcriptional complexes to specific genes; viral transcription thus depends on chromatin dynamics. Influenza virus polymerase associates with chromatin components of the infected cell, such as RNA polymerase II (RNAP II) or the CHD6 chromatin remodeler. Here we show that another CHD family member, CHD1 protein, also interacts with the influenza virus polymerase complex. CHD1 recognizes the H3K4me3 (histone 3 with a trimethyl group in lysine 4) histone modification, a hallmark of active chromatin. Downregulation of CHD1 causes a reduction in viral polymerase activity, viral RNA transcription, and the production of infectious particles. Despite the dependence of influenza virus on cellular transcription, RNAP II is degraded when viral transcription is complete, and recombinant viruses unable to degrade RNAP II show decreased pathogenicity in the murine model. We describe the CHD1-RNAP II association, as well as the parallel degradation of both proteins during infection with viruses showing full or reduced induction of degradation. The H3K4me3 histone mark also decreased during influenza virus infection, whereas a histone mark of inactive chromatin, H3K27me3, remained unchanged. Our results indicate that CHD1 is a positive regulator of influenza virus multiplication and suggest a role for chromatin remodeling in the control of the influenza virus life cycle.

IMPORTANCE

Although influenza virus is not integrated into the genome of the infected cell, it needs continuous cellular transcription to synthesize viral mRNA. This mechanism implies functional association with host genome expression and thus depends on chromatin dynamics. Influenza virus polymerase associates with transcription-related factors, such as RNA polymerase II, and with chromatin remodelers, such as CHD6. We identified the association of viral polymerase with another chromatin remodeler, the CHD1 protein, which positively modulated viral polymerase activity, viral RNA transcription, and virus multiplication. Once viral transcription is complete, RNAP II is degraded in infected cells, probably as a virus-induced mechanism to reduce the antiviral response. CHD1 associated with RNAP II and paralleled its degradation during infection with viruses that induce full or reduced degradation. These findings suggest that RNAP II degradation and CHD1 degradation cooperate to reduce the antiviral response.

Integration of genome scale data for identifying new players in colorectal cancer

Colorectal cancers (CRCs) display a wide variety of genomic aberrations that may be either causally linked to their development and progression, or might serve as biomarkers for their presence. Recent advances in rapid high-throughput genetic and genomic analysis have helped to identify a plethora of alterations that can potentially serve as new cancer biomarkers, and thus help to improve CRC diagnosis, prognosis, and treatment. Each distinct data type (copy number variations, gene and microRNAs expression, CpG island methylation) provides an investigator with a different, partially independent, and complementary view of the entire genome. However, elucidation of gene function will require more information than can be provided by analyzing a single type of data. The integration of knowledge obtained from different sources is becoming increasingly essential for obtaining an interdisciplinary view of large amounts of information, and also for cross-validating experimental results. The integration of numerous types of genetic and genomic data derived from public sources, and via the use of ad-hoc bioinformatics tools and statistical methods facilitates the discovery and validation of novel, informative biomarkers. This combinatory approach will also enable researchers to more accurately and comprehensively understand the associations between different biologic pathways, mechanisms, and phenomena, and gain new insights into the etiology of CRC.

Developmental genes significantly afflicted by aberrant promoter methylation and somatic mutation predict overall survival of late-stage colorectal cancer

Carcinogenesis is an exceedingly complicated process, which involves multi-level dysregulations, including genomics (majorly caused by somatic mutation and copy number variation), DNA methylomics, and transcriptomics. Therefore, only looking into one molecular level of cancer is not sufficient to uncover the intricate underlying mechanisms. With the abundant resources of public available data in the Cancer Genome Atlas (TCGA) database, an integrative strategy was conducted to systematically analyze the aberrant patterns of colorectal cancer on the basis of DNA copy number, promoter methylation, somatic mutation and gene expression. In this study, paired samples in each genomic level were retrieved to identify differentially expressed genes with corresponding genetic or epigenetic dysregulations. Notably, the result of gene ontology enrichment analysis indicated that the differentially expressed genes with corresponding aberrant promoter methylation or somatic mutation were both functionally concentrated upon developmental process, suggesting the intimate association between development and carcinogenesis. Thus, by means of random walk with restart, 37 significant development-related genes were retrieved from a priori-knowledge based biological network. In five independent microarray datasets, Kaplan-Meier survival and Cox regression analyses both confirmed that the expression of these genes was significantly associated with overall survival of Stage III/IV colorectal cancer patients.

Somatic gene copy number alterations in colorectal cancer: new quest for cancer drivers and biomarkers

Colorectal cancer (CRC) results from the accumulation of genetic alterations, and somatic copy number alterations (CNAs) are crucial for the development of CRC. Genome-wide survey of CNAs provides opportunities for identifying cancer driver genes in an unbiased manner. The detection of aberrant CNAs may provide novel markers for the early diagnosis and personalized treatment of CRC. A major challenge in array-based profiling of CNAs is to distinguish the alterations that play causative roles from the random alterations that accumulate during colorectal carcinogenesis. In this view, we systematically discuss the frequent CNAs in CRC, focusing on functional genes that have potential diagnostic, prognostic and therapeutic significance.

Coiled-coil domain containing 68 (CCDC68) demonstrates a tumor-suppressive role in pancreatic ductal adenocarcinoma

Using integrative genomics and functional screening, we identified coiled-coil domain containing 68 (CCDC68) as a novel putative tumor suppressor gene (TSG) in pancreatic ductal adenocarcinoma (PDAC). CCDC68 allelic losses were documented in 48% of primary PDAC patient tumors, 50% of PDAC cell lines and 30% of primary patient derived xenografts. We also discovered a single nucleotide polymorphism (SNP) variant (SNP rs1344011) that leads to exon skipping and generation of an unstable protein isoform CCDC68Δ(69-114) in 31% of PDAC patients. Overexpression of full length CCDC68 (CCDC68(wt)) in PANC-1 and Hs.766T PDAC cell lines lacking CDCC68 expression decreased proliferation and tumorigenicity in scid mice. In contrast, the downregulation of endogenous CCDC68 in MIAPaca-2 cells increased tumor growth rate. These effects were not observed with the deletion-containing isoform, CCDC68Δ(69-114).

Architects of the genome: CHD dysfunction in cancer, developmental disorders and neurological syndromes

Chromatin is vital to normal cells, and its deregulation contributes to a spectrum of human ailments. An emerging concept is that aberrant chromatin regulation culminates in gene expression programs that set the stage for the seemingly diverse pathologies of cancer, developmental disorders and neurological syndromes. However, the mechanisms responsible for such common etiology have been elusive. Recent evidence has implicated lesions affecting chromatin-remodeling proteins in cancer, developmental disorders and neurological syndromes, suggesting a common source for these different pathologies. Here, we focus on the chromodomain helicase DNA binding chromatin-remodeling family and the recent evidence for its deregulation in diverse pathological conditions, providing a new perspective on the underlying mechanisms and their implications for these prevalent human diseases.

Identification of diagnostic markers in colorectal cancer via integrative epigenomics and genomics data

Apart from genetic mutations, epigenetic alteration is a common phenomenon that contributes to neoplastic transformation in colorectal cancer. Transcriptional silencing of tumor-suppressor genes without changes in the DNA sequence is explained by the existence of promoter hypermethylation. To test this hypothesis, we integrated the epigenome and transcriptome data from a similar set of colorectal tissue samples. Methylation profiling was performed using the Illumina InfiniumHumanMethylation27 BeadChip on 55 paired cancer and adjacent normal epithelial cells. Fifteen of the 55 paired tissues were used for gene expression profiling using the Affymetrix GeneChip Human Gene 1.0 ST array. Validation was carried out on 150 colorectal tissues using the methylation-specific multiplex ligation-dependent probe amplification (MS-MLPA) technique. PCA and supervised hierarchical clustering in the two microarray datasets showed good separation between cancer and normal samples. Significant genes from the two analyses were obtained based on a ≥2-fold change and a false discovery rate (FDR) p-value of <0.05. We identified 1,081 differentially hypermethylated CpG sites and 36 hypomethylated CpG sites. We also found 709 upregulated and 699 downregulated genes from the gene expression profiling. A comparison of the two datasets revealed 32 overlapping genes with 27 being hypermethylated with downregulated expression and 4 hypermethylated with upregulated expression. One gene was found to be hypomethylated and downregulated. The most enriched molecular pathway identified was cell adhesion molecules that involved 4 overlapped genes, JAM2, NCAM1, ITGA8 and CNTN1. In the present study, we successfully identified a group of genes that showed methylation and gene expression changes in well-defined colorectal cancer tissues with high purity. The integrated analysis gives additional insight regarding the regulation of colorectal cancer-associated genes and their underlying mechanisms that contribute to colorectal carcinogenesis.

Evaluation of copy number variation and gene expression in neurofibromatosis type-1-associated malignant peripheral nerve sheath tumours

BACKGROUND

Neurofibromatosis type-1 (NF1) is a complex neurogenetic disorder characterised by the development of benign and malignant tumours of the peripheral nerve sheath (MPNSTs). Whilst biallelic NF1 gene inactivation contributes to benign tumour formation, additional cellular changes in gene structure and/or expression are required to induce malignant transformation. Although few molecular profiling studies have been performed on the process of progression of pre-existing plexiform neurofibromas to MPNSTs, the integrated analysis of copy number alterations (CNAs) and gene expression is likely to be key to understanding the molecular mechanisms underlying NF1-MPNST tumorigenesis. In a pilot study, we employed this approach to identify genes differentially expressed between benign and malignant NF1 tumours.

RESULTS

SPP1 (osteopontin) was the most differentially expressed gene (85-fold increase in expression), compared to benign plexiform neurofibromas. Short hairpin RNA (shRNA) knockdown of SPP1 in NF1-MPNST cells reduced tumour spheroid size, wound healing and invasion in four different MPNST cell lines. Seventy-six genes were found to exhibit concordance between CNA and gene expression level.

CONCLUSIONS

Pathway analysis of these genes suggested that glutathione metabolism and Wnt signalling may be specifically involved in NF1-MPNST development. SPP1 is associated with malignant transformation in NF1-associated MPNSTs and could prove to be an important target for therapeutic intervention.

Integrated DNA Copy Number and Gene Expression Regulatory Network Analysis of Non-small Cell Lung Cancer Metastasis

Integrative analysis of multi-level molecular profiles can distinguish interactions that cannot be revealed based on one kind of data in the analysis of cancer susceptibility and metastasis. DNA copy number variations (CNVs) are common in cancer cells, and their role in cell behaviors and relationship to gene expression (GE) is poorly understood. An integrative analysis of CNV and genome-wide mRNA expression can discover copy number alterations and their possible regulatory effects on GE. This study presents a novel framework to identify important genes and construct potential regulatory networks based on these genes. Using this approach, DNA copy number aberrations and their effects on GE in lung cancer progression were revealed. Specifically, this approach contains the following steps: (1) select a pool of candidate driver genes, which have significant CNV in lung cancer patient tumors or have a significant association with the clinical outcome at the transcriptional level; (2) rank important driver genes in lung cancer patients with good prognosis and poor prognosis, respectively, and use top-ranked driver genes to construct regulatory networks with the COpy Number and EXpression In Cancer (CONEXIC) method; (3) identify experimentally confirmed molecular interactions in the constructed regulatory networks using Ingenuity Pathway Analysis (IPA); and (4) visualize the refined regulatory networks with the software package Genatomy. The constructed CNV/mRNA regulatory networks provide important insights into potential CNV-regulated transcriptional mechanisms in lung cancer metastasis.