Aneuploidy-dependent massive deregulation of the cellular transcriptome and apparent divergence of the Wnt/beta-catenin signaling pathway in human rectal carcinomas

Peripheral Blood Transcripts Predict Preoperative Obstructive Total Anomalous Pulmonary Venous Connection

The lack of accessible noninvasive tools to examine the molecular alterations limits our understanding of the causes of total anomalous pulmonary venous connection (TAPVC), as well as the identification of effective operational strategies. Here, we consecutively enrolled peripheral leukocyte transcripts of 26 preoperative obstructive and 22 non-obstructive patients with TAPVC. Two-hundred and fifty six differentially expressed mRNA and 27 differentially expressed long noncoding RNA transcripts were dysregulated. The up-regulated mRNA was enriched in the hydrogen peroxide catabolic process, response to mechanical stimulus, neutrophil degranulation, hemostasis, response to bacterium, and the NABA CORE MATRISOME pathway, all of which are associated with the development of fibrosis. Furthermore, we constructed predictive models using multiple machine-learning algorithms and tested the performance in the validation set. The mRNA NR3C2 and lncRNA MEG3 were screened based on multiple iterations. The random forest prediction model can predict preoperative obstruction patients in the validation set with high accuracy (area under curve = 1; sensitivity = 1). These data highlight the potential of peripheral leukocyte transcripts to evaluate obstructive-related pathophysiological alterations, leading to precision healthcare solutions that could improve patient survival after surgery. It also provides a novel direction for the study of preoperative obstructive TAPVC.

CircRNA expression profiling of PBMCs from patients with hepatocellular carcinoma by RNA-sequencing

Circular RNAs (circRNAs) are differentially expressed in various cancer types. The present study aimed to investigate the expression and clinical implication of circRNAs in hepatocellular carcinoma (HCC) and to evaluate the potential of circRNAs as diagnostic biomarkers for HCC. CircRNA expression was profiled in 19 patients with HCC and 19 normal controls using ribosomal RNA-depleted RNAs. Differentially expressed circRNAs (DE-circRNAs) between HCC and controls were identified using CIRI2 and distinct circRNA expression signatures were screened. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes were used to predict the potential functions of these DE-circRNAs and the circRNA-miRNA-mRNA regulatory networks were then constructed. Several DE-circRNAs were selected and confirmed by RT-qPCR. A total of 40 DE-circRNAs (27 upregulated and 13 downregulated) were identified between patients with HCC and controls. Functional annotation indicated that these DE-circRNAs were involved in cellular components, molecular functions and cancer-associated pathways related to HCC. These included pathways in cancer, TNF signaling pathway, hepatitis B, hepatitis C and hepatocyte differentiation. The circRNA-miRNA-mRNA regulatory network was generated based on 11 candidate circRNAs. Receiver operating characteristic curve analysis indicated that Homo sapiens (hsa)_circ_0073239, hsa_circ_007090, hsa_circ_0008304, hsa_circ_0017586, hsa_circ_0000369 and hsa_circ_0001181 may serve as potential biomarkers for HCC. Results from Cell Counting Kit-8 assay suggested that small interfering RNA targeting hsa_circ_0001181 reduced the proliferation of HepG2 cells, which implicated it as a potential therapeutic target for HCC. Therefore, in the present study, the differential expression pattern and important role of circRNAs in HCC were determined. The present results highlight the diagnostic potential of circRNAs in HCC and provide novel insight into the development of and treatment approaches for HCC.

Impact of storage conditions on peripheral leukocytes transcriptome

Leukocytes reflect the physiological and pathological states of each individual, and transcriptomic data of leukocytes have been used to reflect health conditions. Since the overall impact of ex vivo conditions on the leukocyte transcriptome before RNA stabilization remains unclear, we evaluated the influence of temporary storage conditions on the leukocyte transcriptome through RNA sequencing. We collected peripheral blood with EDTA tubes, which were processed immediately or stored either at 4 °C or room temperature (RT, 18-22 °C) for 2 h, 6 h and 24 h. Total cellular RNA was extracted from 42 leukocyte samples after red blood cells lysis for subsequent RNA sequencing. We applied weighted gene co-expression network analysis to construct co-expression networks of mRNA and lncRNA among the samples, and then performed gene ontology (GO) term enrichment to explore possible biological processes affected by storage conditions. Storage conditions change the gene expression of peripheral leukocytes. Comparing with fresh leukocytes, storage for 24 h at 4 °C and RT affected 1515 (1.51%) and 10,823 (10.82%) genes, respectively. Pathway enrichment analysis identified nucleosome assembly enriched in up-regulated genes at both conditions. When blood was stored at RT for 24 h, genes involved in apoptotic signaling pathway, negative regulation of cell cycle and lymphocyte activation were upregulated, while the relative proportion of neutrophils was significantly decreased. Temporary storage conditions profoundly affect the gene expression profiles of leukocytes and might further change cell viability and state. Storage of blood samples at 4 °C within 6 h largely maintains their original transcriptome.

Chromosomal Density of Cancer Up-Regulated Genes, Aberrant Enhancer Activity and Cancer Fitness Genes Are Associated with Transcriptional Cis-Effects of Broad Copy Number Gains in Colorectal Cancer

Broad Copy Number Gains (BCNGs) are copy-number increases of chromosomes or large segments of chromosomal arms. Publicly-available single-nucleotide polymorphism (SNP) array and RNA-Seq data of colon adenocarcinoma (COAD) samples from The Cancer Genome Atlas (TCGA) consortium allowed us to design better control groups in order to identify changes in expression due to highly recurrent BCNGs (in chromosomes 20, 8, 7, 13). We identified: (1) Overexpressed Transcripts (OverT), transcripts whose expression increases in "COAD groups bearing a specific BCNG" in comparison to "control COAD groups" not bearing it, and (2) up-regulated/down-regulated transcripts, transcripts whose expression increases/decreases in COAD groups in comparison to normal colon tissue. An analysis of gene expression reveals a correlation between the density of up-regulated genes per selected chromosome and the recurrence rate of their BCNGs. We report an enrichment of gained enhancer activity and of cancer fitness genes among OverT genes. These results support the hypothesis that the chromosomal density of overexpressed cancer fitness genes might play a significant role in the selection of gained chromosomes during cancer evolution. Analysis of functional pathways associated with OverT suggest that some multi-subunit protein complexes (eIF2, eIF3, CSTF and CPSF) are candidate targets for silencing transcriptional therapy.

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.

Single Chromosome Aneuploidy Induces Genome-Wide Perturbation of Nuclear Organization and Gene Expression

Chromosomal aneuploidy is a defining feature of carcinomas and results in tumor-entity specific genomic imbalances. For instance, most sporadic colorectal carcinomas carry extra copies of chromosome 7, an aneuploidy that emerges already in premalignant adenomas, and is maintained throughout tumor progression and in derived cell lines. A comprehensive understanding on how chromosomal aneuploidy affects nuclear organization and gene expression, i.e., the nucleome, remains elusive. We now analyzed a cell line established from healthy colon mucosa with a normal karyotype (46,XY) and its isogenic derived cell line that acquired an extra copy of chromosome 7 as its sole anomaly (47,XY,+7). We studied structure/function relationships consequent to aneuploidization using genome-wide chromosome conformation capture (Hi-C), RNA sequencing and protein profiling. The gain of chromosome 7 resulted in an increase of transcript levels of resident genes as well as genome-wide gene and protein expression changes. The Hi-C analysis showed that the extra copy of chromosome 7 is reflected in more interchromosomal contacts between the triploid chromosomes. Chromatin organization changes are observed genome-wide, as determined by changes in A/B compartmentalization and topologically associating domain (TAD) boundaries. Most notably, chromosome 4 shows a profound loss of chromatin organization, and chromosome 14 contains a large A/B compartment switch region, concurrent with resident gene expression changes. No changes to the nuclear position of the additional chromosome 7 territory were observed when measuring distances of chromosome painting probes by interphase FISH. Genome and protein data showed enrichment in signaling pathways crucial for malignant transformation, such as the HGF/MET-axis. We conclude that a specific chromosomal aneuploidy has profound impact on nuclear structure and function, both locally and genome-wide. Our study provides a benchmark for the analysis of cancer nucleomes with complex karyotypes.

Impact of RNA integrity and blood sample storage conditions on the gene expression analysis

Background

The reliability of RNA sequencing (RNA-seq) output is affected by the quality of RNAs, which is in turn dependent on the quality of samples. Therefore, the purposes of this study were to reconsider the threshold of the RNA integrity number (RIN) and propose a simple and efficient storage scheme of blood samples for RNA-seq.

Patients and methods

The RNAs were extracted from blood samples that were stored at different conditions and used for sequencing. The bioinformatic analyses were performed to evaluate the impact of RNA integrity and blood sample storage conditions on the gene expression analysis.

Results

Our outcomes showed that the samples with RIN values more than 5.3 scarcely affected the quantitative results of RNA-seq, and the influence of inherent cellular physiological processes on RNA-seq output could be negligible.

Conclusion

The blood samples stored at 4°C within 7 days with RIN values more than 5.3 were available for RNA-seq.

Screening effective differential expression genes for hepatic carcinoma with metastasis in the peripheral blood mononuclear cells by RNA-seq

Tumor metastasis is a multistep process involving a number of genetic alterations so that the genetic diagnosis is got increasingly attentions today. The aim of this study was to use RNA-seq to screen the effective differential expression genes in the peripheral blood mononuclear cells for the hepatic carcinoma with metastasis. The results showed that hepatic carcinoma samples gathered according to different metastasis. CCL3, CCL3L1, JUN, IL8, and IL1B were identified in inflammation mediated by chemokine and cytokine signaling pathway (P00031) in the hepatic carcinoma samples with metastasis, and subsequently confirmed by quantitative real-time polymerase chain reaction. In conclusions, CCL3, CCL3L1, JUN, IL8, and IL1B have the potential to be considered as candidates for future molecular diagnosis of the hepatic carcinoma with metastasis. This work may provide us with new visions into the metastasis process and potential efficient clinical diagnosis in the future.

"Vessels in the Storm": Searching for Prognostic and Predictive Angiogenic Factors in Colorectal Cancer

High expectations are placed upon anti-angiogenic compounds for metastatic colorectal cancer (mCRC), the first malignancy for which such type of treatment has been approved. Indeed, clinical trials have confirmed that targeting the formation of new vessels can improve in many cases clinical outcomes of mCRC patients. However, current anti-angiogenic drugs are far from obtaining the desirable or expected curative results. Many are the factors probably involved in such disappointing results, but particular attention is currently focused on the validation of biomarkers able to improve the direction of treatment protocols. Because clinical studies have clearly demonstrated that serum or tissue concentration of some angiogenic factors is associated with the evolution of the disease of mCRC patients, they are currently explored as potential biomarkers of prognosis and of tumor response to therapy. However, the complex biology underlying CRC -induced angiogenesis is a hurdle in finding rapid solutions. The aim of this review was to explore molecular mechanisms that determine the formation of tumor-associated vessels during CRC progression, and to discuss the potential role of angiogenic factors as diagnostic, prognostic and predictive biomarkers in CRC.

Patterns of Chromosomal Aberrations in Solid Tumors

Chromosomal abnormalities are a defining feature of solid tumors. Such cytogenetic alterations are mainly classified into structural chromosomal aberrations and copy number alterations, giving rise to aneuploid karyotypes. The increasing detection of these genetic changes allowed the description of specific tumor entities and the associated patterns of gene expression. In fact, tumor-specific landscapes of gross genomic copy number changes, including aneuploidies of entire chromosome arms and chromosomes result in a global deregulation of the transcriptome of cancer cells. Furthermore, the molecular characterization of cytogenetic abnormalities has provided insights into the mechanisms of tumorigenesis and has, in a few instances, led to the clinical implementation of effective diagnostic and prognostic tools, as well as treatment strategies that target a specific genetic abnormality.

Consequences of Aneuploidy in Cancer: Transcriptome and Beyond

Cancer cells differ from normal healthy cells in multiple aspects ranging from altered cellular signaling through metabolic changes to aberrant chromosome content, so called aneuploidy. The large-scale changes in copy numbers of chromosomes or large chromosomal regions due to aneuploidy alter significantly the gene expression, as several hundreds of genes are gained or lost. Comparison of quantitative genome, transcriptome and proteome data enables dissection of the molecular causes that underlie the gene expression changes observed in cancer cells and provides a new perspective on the molecular consequences of aneuploidy. Here, we will map to what degree aneuploidy affects the expression of genes located on the affected chromosomes. We will also address the effects of aneuploidy on global gene expression in cancer cells as well as whether and how it may contribute to the physiology of cancer cells.

Detection of candidate tumor driver genes using a fully integrated Bayesian approach

DNA copy number alterations (CNAs), including amplifications and deletions, can result in significant changes in gene expression and are closely related to the development and progression of many diseases, especially cancer. For example, CNA-associated expression changes in certain genes (called candidate tumor driver genes) can alter the expression levels of many downstream genes through transcription regulation and cause cancer. Identification of such candidate tumor driver genes leads to discovery of novel therapeutic targets for personalized treatment of cancers. Several approaches have been developed for this purpose by using both copy number and gene expression data. In this study, we propose a Bayesian approach to identify candidate tumor driver genes, in which the copy number and gene expression data are modeled together, and the dependency between the two data types is modeled through conditional probabilities. The proposed joint modeling approach can identify CNA and differentially expressed genes simultaneously, leading to improved detection of candidate tumor driver genes and comprehensive understanding of underlying biological processes. We evaluated the proposed method in simulation studies, and then applied to a head and neck squamous cell carcinoma data set. Both simulation studies and data application show that the joint modeling approach can significantly improve the performance in identifying candidate tumor driver genes, when compared with other existing approaches.

Identification of a microRNA expression signature for chemoradiosensitivity of colorectal cancer cells, involving miRNAs-320a, -224, -132 and let7g

BACKGROUND AND PURPOSE

Preoperative chemoradiotherapy (CRT) represents the standard treatment for locally advanced rectal cancer. Tumor response and progression vary considerably. MicroRNAs represent master regulators of gene expression, and may therefore contribute to this diversity.

MATERIAL AND METHODS

Genome-wide microRNA (miRNA) profiling was performed for 12 colorectal cancer (CRC) cell lines and an individual in vitro signature of chemoradiosensitivity was established. Functional relevance of selected miRNAs was established by transfecting miRNA-mimics into SW480 and SW837 cells. The prognostic value of selected miRNAs was assessed in 128 pretherapeutic patient biopsies.

RESULTS

Thirty-six miRNAs were identified to significantly correlate with sensitivity to CRT (Q < 0.05) including miR-320a and other miRNAs involved in the MAPK-, TGF- and Wnt-pathway. Transfection of selected miRNAs (let-7g, miR-132, miR-224, miR-320a) each induced a shift of sensitivity. High expression of let-7 g was associated with a good prognosis in rectal cancer patients (P = 0.03).

CONCLUSIONS

This is the first report of a miRNA expression signature for in vitro chemoradiosensitivity of CRC cell lines. Many of the identified miRNAs have not been linked to the response to CRT and may represent potential molecular targets to sensitize resistant cancers. If further validated, let7g expression may serve as predictive biomarker.

The role of aberrant promoter hypermethylation of DACT1 in bladder urothelial carcinoma

The purpose of this study was to determine the relationship between hypermethylation of DACT1 gene promoter and lower mRNA expression in bladder urothelial carcinoma tissue. The methylation status of 29 urothelial carcinoma samples and 29 normal tissue samples were examined by methylation-specific polymerase chain reaction (MSP). The DACT1 mRNA transcript levels and DACT1 protein levels in all samples were then evaluated to define the relationship between the methylation status of the DACT1 promoter and its expression at the transcriptional and translational levels. Decreased expression of DACT1 was detected in 89.66% of urothelial carcinomas (26/29; P < 0.005). Promoter hypermethylation was found in 58.62% (17/29) urothelial carcinomas and 25% (7/29) normal tissues, respectively (P < 0.05). DACT1 expression was lower in tissues where the DACT1 gene promoter was hypermethylated than in unmethylated tissues (0.25±0.17 vs 0.69±0.30, P < 0.05). DACT1 gene hypermethylation was closely related to tumor size, grade and stage (P < 0.05). Our results indicate that silencing and downregulation of DACT1 mRNA may be implicated in carcinogenesis and the progression of bladder urothelial carcinoma, and may be a potential prognostic factor.

Integrated analysis of genome-wide copy number alterations and gene expression in microsatellite stable, CpG island methylator phenotype-negative colon cancer

Microsatellite stable (MSS), CpG island methylator phenotype (CIMP)-negative colorectal tumors, the most prevalent molecular subtype of colorectal cancer, are associated with extensive copy number alteration (CNA) events and aneuploidy. We report on the identification of characteristic recurrent CNA (with frequency >25%) events and associated gene expression profiles for a total of 40 paired tumor and adjacent normal colon tissues using genome-wide microarrays. We observed recurrent CNAs, namely gains at 1q, 7p, 7q, 8p12-11, 8q, 12p13, 13q, 20p, 20q, Xp, and Xq and losses at 1p36, 1p31, 1p21, 4p15-12, 4q12-35, 5q21-22, 6q26, 8p, 14q, 15q11-12, 17p, 18p, 18q, 21q21-22, and 22q. Within these genomic regions we identified 356 genes with significant differential expression (P < 0.0001 and ±1.5-fold change) in the tumor compared to adjacent normal tissue. Gene ontology and pathway analyses indicated that many of these genes were involved in functional mechanisms that regulate cell cycle, cell death, and metabolism. An amplicon present in >70% of the tumor samples at 20q11-20q13 contained several cancer-related genes (AHCY, POFUT1, RPN2, TH1L, and PRPF6) that were upregulated and demonstrated a significant linear correlation (P < 0.05) for gene dosage and gene expression. Copy number loss at 8p, a CNA associated with adenocarcinoma and poor prognosis, was observed in >50% of the tumor samples and demonstrated a significant linear correlation for gene dosage and gene expression for two potential tumor suppressor genes, MTUS1 (8p22) and PPP2CB (8p12). The results from our integration analysis illustrate the complex relationship between genomic alterations and gene expression in colon cancer.

The consequences of chromosomal aneuploidy on the transcriptome of cancer cells

Chromosomal aneuploidies are a defining feature of carcinomas, i.e., tumors of epithelial origin. Such aneuploidies result in tumor specific genomic copy number alterations. The patterns of genomic imbalances are tumor specific, and to a certain extent specific for defined stages of tumor development. Genomic imbalances occur already in premalignant precursor lesions, i.e., before the transition to invasive disease, and their distribution is maintained in metastases, and in cell lines derived from primary tumors. These observations are consistent with the interpretation that tumor specific genomic imbalances are drivers of malignant transformation. Naturally, this precipitates the question of how such imbalances influence the expression of resident genes. A number of laboratories have systematically integrated copy number alterations with gene expression changes in primary tumors and metastases, cell lines, and experimental models of aneuploidy to address the question as to whether genomic imbalances deregulate the expression of one or few key genes, or rather affect the cancer transcriptome more globally. The majority of these studies showed that gene expression levels follow genomic copy number. Therefore, gross genomic copy number changes, including aneuploidies of entire chromosome arms and chromosomes, result in a massive deregulation of the transcriptome of cancer cells. This article is part of a Special Issue entitled: Chromatin in time and space.

The molecular basis of chemoradiosensitivity in rectal cancer: implications for personalized therapies

Introduction

Preoperative chemoradiotherapy represents the standard treatment for patients with locally advanced rectal cancer. Unfortunately, the response of individual tumors to multimodal treatment is not uniform and ranges from complete response to complete resistance. This poses a particular problem for patients with a priori resistant tumors because they may be exposed to irradiation and chemotherapy, treatment regimens that are both expensive and at times toxic, without benefit. Accordingly, there is a strong need to establish molecular biomarkers that predict the response of an individual patient’s tumor to multimodal treatment and that indicate treatment-associated toxicities prior to therapy. Such biomarkers may guide clinicians in choosing the best possible treatment for each individual patient. In addition, these biomarkers could be used to identify novel molecular targets and thereby assist in implementing novel strategies to sensitize a priori resistant tumors to multimodal treatment regimens.

Objective

The aim of this review is to summarize recent findings about the molecular basis of treatment resistance and treatment toxicity in patients with rectal cancer. Whole-genome, as well as single-biomarker or multibiomarker, analyses and their potential implications will be highlighted. At the end, we will outline a future vision of rectal cancer treatment in the era of personalized medicine.

Genomic alterations in rectal tumors and response to neoadjuvant chemoradiotherapy: an exploratory study

Background

Neoadjuvant chemoradiotherapy is the treatment of choice in advanced rectal cancer, even though there are many patients who will not benefit from it. There are still no effective methods for predicting which patients will respond or not. The present study aimed to define the genomic profile of rectal tumors and to identify alterations that are predictive of response in order to optimize therapeutic strategies.

Methods

Forty-eight candidates for neoadjuvant chemoradiotherapy were recruited and their pretherapy biopsies analyzed by array Comparative Genomic Hybridization (aCGH). Pathologic response was evaluated by tumor regression grade.

Results

Both Hidden Markov Model and Smoothing approaches identified similar alterations, with a prevalence of DNA gains. Non responsive patients had a different alteration profile from responsive ones, with a higher number of genome changes mainly located on 2q21, 3q29, 7p22-21, 7q21, 7q36, 8q23-24, 10p14-13, 13q12, 13q31-34, 16p13, 17p13-12 and 18q23 chromosomal regions.

Conclusions

This exploratory study suggests that an in depth characterization of chromosomal alterations by aCGH would provide useful predictive information on response to neoadjuvant chemoradiotherapy and could help to optimize therapy in rectal cancer patients.

The data discussed in this study are available on the NCBI Gene Expression Omnibus [GEO: GSE25885].

A genomic strategy for the functional validation of colorectal cancer genes identifies potential therapeutic targets

Genes that are highly overexpressed in tumor cells can be required for tumor cell survival and have the potential to be selective therapeutic targets. In an attempt to identify such targets, we combined a functional genomics and a systems biology approach to assess the consequences of RNAi-mediated silencing of overexpressed genes that were selected from 140 gene expression profiles from colorectal cancers (CRCs) and matched normal mucosa. In order to identify credible models for in-depth functional analysis, we first confirmed the overexpression of these genes in 25 different CRC cell lines. We then identified five candidate genes that profoundly reduced the viability of CRC cell lines when silenced with either siRNAs or short-hairpin RNAs (shRNAs), i.e., HMGA1, TACSTD2, RRM2, RPS2 and NOL5A. These genes were further studied by systematic analysis of comprehensive gene expression profiles generated following siRNA-mediated silencing. Exploration of these RNAi-specific gene expression signatures allowed the identification of the functional space in which the five genes operate and showed enrichment for cancer-specific signaling pathways, some known to be involved in CRC. By comparing the expression of the RNAi signature genes with their respective expression levels in an independent set of primary rectal carcinomas, we could recapitulate these defined RNAi signatures, therefore, establishing the biological relevance of our observations. This strategy identified the signaling pathways that are affected by the prominent oncogenes HMGA1 and TACSTD2, established a yet unknown link between RRM2 and PLK1 and identified RPS2 and NOL5A as promising potential therapeutic targets in CRC.

A genomic strategy for the functional validation of colorectal cancer genes identifies potential therapeutic targets

Genes that are highly overexpressed in tumor cells can be required for tumor cell survival and have the potential to be selective therapeutic targets. In an attempt to identify such targets, we combined a functional genomics and a systems biology approach to assess the consequences of RNAi-mediated silencing of overexpressed genes that were selected from 140 gene expression profiles from colorectal cancers (CRCs) and matched normal mucosa. In order to identify credible models for in-depth functional analysis, we first confirmed the overexpression of these genes in 25 different CRC cell lines. We then identified five candidate genes that profoundly reduced the viability of CRC cell lines when silenced with either siRNAs or short-hairpin RNAs (shRNAs), i.e., HMGA1, TACSTD2, RRM2, RPS2 and NOL5A. These genes were further studied by systematic analysis of comprehensive gene expression profiles generated following siRNA-mediated silencing. Exploration of these RNAi-specific gene expression signatures allowed the identification of the functional space in which the five genes operate and showed enrichment for cancer-specific signaling pathways, some known to be involved in CRC. By comparing the expression of the RNAi signature genes with their respective expression levels in an independent set of primary rectal carcinomas, we could recapitulate these defined RNAi signatures, therefore, establishing the biological relevance of our observations. This strategy identified the signaling pathways that are affected by the prominent oncogenes HMGA1 and TACSTD2, established a yet unknown link between RRM2 and PLK1 and identified RPS2 and NOL5A as promising potential therapeutic targets in CRC.

Mutated KRAS results in overexpression of DUSP4, a MAP-kinase phosphatase, and SMYD3, a histone methyltransferase, in rectal carcinomas

Mutations of the KRAS oncogene are predictive for resistance to treatment with antibodies against the epithelial growth factor receptor in patients with colorectal cancer. Overcoming this therapeutic dilemma could potentially be achieved by the introduction of drugs that inhibit signaling pathways that are activated by KRAS mutations. To identify comprehensively such signaling pathways, we profiled pretreatment biopsies and normal mucosa from 65 patients with locally advanced rectal cancer-30 of which carried mutated KRAS-using global gene expression microarrays. By comparing all tumor tissues exclusively to matched normal mucosa, we could improve assay sensitivity, and identified a total of 22,297 features that were differentially expressed (adjusted P-value <0.05) between normal mucosa and cancer, including several novel potential rectal cancer genes. We then used this comprehensive description of the rectal cancer transcriptome as the baseline for identifying KRAS-dependent alterations. The presence of activating KRAS mutations is significantly correlated to an upregulation of 13 genes (adjusted P-value <0.05), among them DUSP4, a MAP-kinase phosphatase, and SMYD3, a histone methyltransferase. Inhibition of the expression of both genes has previously been shown using the MEK1-inhibitor PD98059 and the antibacterial compound Novobiocin, respectively. These findings suggest a potential approach to overcome resistance to treatment with antibodies against the epithelial growth factor receptor in patients with KRAS-mutant rectal carcinomas.

-8p12-23 and +20q are predictors of subtypes and metastatic pathways in colorectal cancer: construction of tree models using comparative genomic hybridization data

A substantial body of evidence suggests the genetic heterogeneous pattern and multiple pathways in colorectal cancer initiation and progression. In this study, we construct a branching tree and multiple distance-based tree models to elucidate these genetic patterns and pathways in colorectal cancer by using a data set comprised of 244 cases of comparative genomic hybridization. We identify the six most common gains of chromosomal regions of 7p (37.0%), 7q11-32 (34.8%), 8q (48.3%), 13q (49.1%), 20p (36.1%), and 20q (50.4%), and the nine most common losses of 1p13-36 (30.9%), 4p15 (24.3%), 4q33-34 (24.3%), 8p12-23 (50.9%), 15q13-14 (23.5%), 15q24-25 (24.3%), 17p (34.8%), 18p (36.5%), and 18q (61.7%) in colorectal cancer. We classify colorectal cancer into two distinct groups: one preceding with -8p12-23, and the other with +20q. The sample-based classification tree also demonstrates that colorectal cancer can be classified into multiple subtypes marked by -8p12-23 and +20q. By comparing chromosomal abnormalities between primary and metastatic colorectal cancer, we identify five potential metastatic pathways: (-18q, -18p), (-8p12-23, -4p15, -4q33-34), (+20q, +20p), (+20q, +7p, +7q11-32), and +8q. -8p12-23 and +20q are inferred to be the two marker events of colorectal cancer metastasis. The current oncogenetic tree models may contribute to our understanding towards molecular genetics in colorectal cancer. Particularly, the metastatic pathways we describe may provide pivotal clues for metastatic candidate genes, and thus impact on the prediction and intervention of metastatic colorectal cancer.

A gene expression signature for chemoradiosensitivity of colorectal cancer cells

PURPOSE

The standard treatment of patients with locally advanced rectal cancers comprises preoperative 5-fluorouracil-based chemoradiotherapy followed by standardized surgery. However, tumor response to multimodal treatment has varied greatly, ranging from complete resistance to complete pathologic regression. The prediction of the response is, therefore, an important clinical need.

METHODS AND MATERIALS

To establish in vitro models for studying the molecular basis of this heterogeneous tumor response, we exposed 12 colorectal cancer cell lines to 3 μM of 5-fluorouracil and 2 Gy of radiation. The differences in treatment sensitivity were then correlated with the pretherapeutic gene expression profiles of these cell lines.

RESULTS

We observed a heterogeneous response, with surviving fractions ranging from 0.28 to 0.81, closely recapitulating clinical reality. Using a linear model analysis, we identified 4,796 features whose expression levels correlated significantly with the sensitivity to chemoradiotherapy (Q <.05), including many genes involved in the mitogen-activated protein kinase signaling pathway or cell cycle genes. These data have suggested a potential relevance of the insulin and Wnt signaling pathways for treatment response, and we identified STAT3, RASSF1, DOK3, and ERBB2 as potential therapeutic targets. The microarray measurements were independently validated for a subset of these genes using real-time polymerase chain reactions.

CONCLUSION

We are the first to report a gene expression signature for the in vitro chemoradiosensitivity of colorectal cancer cells. We anticipate that this analysis will unveil molecular biomarkers predictive of the response of rectal cancers to chemoradiotherapy and enable the identification of genes that could serve as targets to sensitize a priori resistant primary tumors.

Distinguishing between cancer driver and passenger gene alteration candidates via cross-species comparison: a pilot study

Background

We are developing a cross-species comparison strategy to distinguish between cancer driver- and passenger gene alteration candidates, by utilizing the difference in genomic location of orthologous genes between the human and other mammals. As an initial test of this strategy, we conducted a pilot study with human colorectal cancer (CRC) and its mouse model C57BL/6J Apc^Min/+, focusing on human 5q22.2 and 18q21.1-q21.2.

Methods

We first performed bioinformatics analysis on the evolution of 5q22.2 and 18q21.1-q21.2 regions. Then, we performed exon-targeted sequencing, real time quantitative polymerase chain reaction (qPCR), and real time quantitative reverse transcriptase PCR (qRT-PCR) analyses on a number of genes of both regions with both human and mouse colon tumors.

Results

These two regions (5q22.2 and 18q21.1-q21.2) are frequently deleted in human CRCs and encode genuine colorectal tumor suppressors APC and SMAD4. They also encode genes such as MCC (mutated in colorectal cancer) with their role in CRC etiology unknown. We have discovered that both regions are evolutionarily unstable, resulting in genes that are clustered in each human region being found scattered at several distinct loci in the genome of many other species. For instance, APC and MCC are within 200 kb apart in human 5q22.2 but are 10 Mb apart in the mouse genome. Importantly, our analyses revealed that, while known CRC driver genes APC and SMAD4 were disrupted in both human colorectal tumors and tumors from Apc^Min/+ mice, the questionable MCC gene was disrupted in human tumors but appeared to be intact in mouse tumors.

Conclusions

These results indicate that MCC may not actually play any causative role in early colorectal tumorigenesis. We also hypothesize that its disruption in human CRCs is likely a mere result of its close proximity to APC in the human genome. Expanding this pilot study to the entire genome may identify more questionable genes like MCC, facilitating the discovery of new CRC driver gene candidates.

Impact of RNA degradation on gene expression profiling

BACKGROUND

Gene expression profiling is a highly sensitive technique which is used for profiling tumor samples for medical prognosis. RNA quality and degradation influence the analysis results of gene expression profiles. The impact of this influence on the profiles and its medical impact is not fully understood. As patient samples are very valuable for clinical studies, it is necessary to establish criteria for the RNA quality to be able to use these samples in later analysis.

METHODS

To investigate the effects of RNA integrity on gene expression profiling, whole genome expression arrays were used. We used tumor biopsies from patients diagnosed with locally advanced rectal cancer. To simulate degradation, the isolated total RNA of all patients was subjected to heat-induced degradation in a time-dependent manner. Expression profiling was then performed and data were analyzed bioinformatically to assess the differences.

RESULTS

The differences introduced by RNA degradation were largely outweighed by the biological differences between the patients. Only a relatively small number of probes (275 out of 41,000) show a significant effect due to degradation. The genes that show the strongest effect due to RNA degradation were, especially, those with short mRNAs and probe positions near the 5' end.

CONCLUSIONS

Degraded RNA from tumor samples (RIN > 5) can still be used to perform gene expression analysis. A much higher biological variance between patients is observed compared to the effect that is imposed by degradation of RNA. Nevertheless there are genes, very short ones and those with the probe binding side close to the 5' end that should be excluded from gene expression analysis when working with degraded RNA. These results are limited to the Agilent 44 k microarray platform and should be carefully interpreted when transferring to other settings.

Integrative genomics reveals mechanisms of copy number alterations responsible for transcriptional deregulation in colorectal cancer

To evaluate the mechanisms and consequences of chromosomal aberrations in colorectal cancer (CRC), we used a combination of spectral karyotyping, array comparative genomic hybridization (aCGH), and array-based global gene expression profiling on 31 primary carcinomas and 15 established cell lines. Importantly, aCGH showed that the genomic profiles of primary tumors are recapitulated in the cell lines. We revealed a preponderance of chromosome breakpoints at sites of copy number variants (CNVs) in the CRC cell lines, a novel mechanism of DNA breakage in cancer. The integration of gene expression and aCGH led to the identification of 157 genes localized within high-level copy number changes whose transcriptional deregulation was significantly affected across all of the samples, thereby suggesting that these genes play a functional role in CRC. Genomic amplification at 8q24 was the most recurrent event and led to the overexpression of MYC and FAM84B. Copy number dependent gene expression resulted in deregulation of known cancer genes such as APC, FGFR2, and ERBB2. The identification of only 36 genes whose localization near a breakpoint could account for their observed deregulated expression demonstrates that the major mechanism for transcriptional deregulation in CRC is genomic copy number changes resulting from chromosomal aberrations.

KRAS and BRAF mutations in patients with rectal cancer treated with preoperative chemoradiotherapy

BACKGROUND AND PURPOSE

KRAS and BRAF are mutated in 35% and 10% of colorectal cancers, respectively. However, data specifically for locally advanced rectal cancers are scarce, and the frequency of KRAS mutations in codons 61 and 146 remains to be established.

MATERIALS AND METHODS

DNA was isolated from pre-therapeutic biopsies of 94 patients who were treated within two phase-III clinical trials receiving preoperative chemoradiotherapy. Mutation status of KRAS exons 1-3 and BRAF exon 15 was established using the ABI PRISM Big Dye Sequencing Kit and subsequently correlated with clinical parameters.

RESULTS

Overall, KRAS was mutated in 45 patients (48%). Twenty-nine mutations (64%) were located in codon 12, 10 mutations (22%) in codon 13, and 3 mutations (7%) in codons 61 and 146. No V600E BRAF mutation was detected. The presence of KRAS mutations was correlated neither with tumor response or lymph node status after preoperative chemoradiotherapy nor with overall survival or disease-free survival. When KRAS exon 1 mutations were separated based on the amino-acid exchange, we again failed to detect significant correlations (p=0.052). However, G12V mutations appeared to be associated with higher rates of tumor regression than G13D mutations (p=0.012).

CONCLUSION

We are the first to report the mutation status of KRAS and BRAF in pre-therapeutic biopsies from locally advanced rectal cancers. The high number of KRAS mutations in codons 61 and 146 emphasizes the importance to expand current mutation analyses, whereas BRAF mutations are not relevant for rectal carcinogenesis. Although the KRAS mutation status was not correlated with response, the subtle difference between G12V and G13D mutations warrants analysis of a larger patient population.

[Individualizing treatment for locally advanced rectal cancer]

Based on results of the German Rectal Cancer Study Group CAO/ARO/AIO-94 trial, long-term chemoradiotherapy (RT/CTx) is recommended as standard treatment for locally advanced rectal cancer (UICC stages II/III) in the lower two thirds of the rectum (0-12 cm from the anocutaneous verge). Tumor response to neoadjuvant therapy is very heterogeneous, ranging from complete remission to total resistance to RT/CTx. To fulfill the clinical requirement of individual and risk-adapted multimodal treatment, distinct progress in translational research has been achieved (e.g. gene profiling). However, in clinical reality "individualization" of the therapy of rectal cancer patients has not actually been realized. This can be achieved only on the basis of successful randomized clinical trials (e.g. the CAO/ARO/AIO-04 and GAST-05 trials) translationally combined with basic scientific approaches. One simple first step toward individualizing rectal cancer therapy is being made with the ongoing GAST-05 trial. This investigator initiated phase II trial funded by the German Research Foundation (Deutsche Forschungsgemeinschaft) excludes preoperative RT/CTx for patients with rectal cancer localized in the upper third of the rectum, using only quality controlled principles of radical surgery (partial vs total mesorectal excision) followed by adjuvant chemotherapy.

Can survival prediction be improved by merging gene expression data sets?

BACKGROUND

High-throughput gene expression profiling technologies generating a wealth of data, are increasingly used for characterization of tumor biopsies for clinical trials. By applying machine learning algorithms to such clinically documented data sets, one hopes to improve tumor diagnosis, prognosis, as well as prediction of treatment response. However, the limited number of patients enrolled in a single trial study limits the power of machine learning approaches due to over-fitting. One could partially overcome this limitation by merging data from different studies. Nevertheless, such data sets differ from each other with regard to technical biases, patient selection criteria and follow-up treatment. It is therefore not clear at all whether the advantage of increased sample size outweighs the disadvantage of higher heterogeneity of merged data sets. Here, we present a systematic study to answer this question specifically for breast cancer data sets. We use survival prediction based on Cox regression as an assay to measure the added value of merged data sets.

RESULTS

Using time-dependent Receiver Operating Characteristic-Area Under the Curve (ROC-AUC) and hazard ratio as performance measures, we see in overall no significant improvement or deterioration of survival prediction with merged data sets as compared to individual data sets. This apparently was due to the fact that a few genes with strong prognostic power were not available on all microarray platforms and thus were not retained in the merged data sets. Surprisingly, we found that the overall best performance was achieved with a single-gene predictor consisting of CYB5D1.

CONCLUSIONS

Merging did not deteriorate performance on average despite (a) The diversity of microarray platforms used. (b) The heterogeneity of patients cohorts. (c) The heterogeneity of breast cancer disease. (d) Substantial variation of time to death or relapse. (e) The reduced number of genes in the merged data sets. Predictors derived from the merged data sets were more robust, consistent and reproducible across microarray platforms. Moreover, merging data sets from different studies helps to better understand the biases of individual studies and can lead to the identification of strong survival factors like CYB5D1 expression.

DATE analysis: A general theory of biological change applied to microarray data

In contrast to conventional data mining, which searches for specific subsets of genes (extensive variables) to correlate with specific phenotypes, DATE analysis correlates intensive state variables calculated from the same datasets. At the heart of DATE analysis are two biological equations of state not dependent on genetic pathways. This result distinguishes DATE analysis from other bioinformatics approaches. The dimensionless state variable F quantifies the relative overall cellular activity of test cells compared to well-chosen reference cells. The variable pi(i) is the fold-change in the expression of the ith gene of test cells relative to reference. It is the fraction phi of the genome undergoing differential expression-not the magnitude pi-that controls biological change. The state variable phi is equivalent to the control strength of metabolic control analysis. For tractability, DATE analysis assumes a linear system of enzyme-connected networks and exploits the small average contribution of each cellular component. This approach was validated by reproducible values of the state variables F, RNA index, and phi calculated from random subsets of transcript microarray data. Using published microarray data, F, RNA index, and phi were correlated with: (1) the blood-feeding cycle of the malaria parasite, (2) embryonic development of the fruit fly, (3) temperature adaptation of Killifish, (4) exponential growth of cultured S. pneumoniae, and (5) human cancers. DATE analysis was applied to aCGH data from the great apes. A good example of the power of DATE analysis is its application to genomically unstable cancers, which have been refractory to data mining strategies.

Chromosomal copy number changes of locally advanced rectal cancers treated with preoperative chemoradiotherapy

Standard treatment of rectal cancer patients comprises preoperative chemoradiotherapy followed by radical surgery. However, clinicians are faced with the problem that response rates vary from one individual to another. Predictive biomarkers would therefore be helpful. To identify genomic imbalances that might assist in stratifying tumors into responsive or nonresponsive categories, we used metaphase comparative genomic hybridization to prospectively analyze pretherapeutic biopsies from 42 patients with locally advanced rectal cancers. These patients were subsequently treated with 5-fluorouracil-based preoperative chemoradiotherapy. Based on downsizing of the T-category, 21 rectal cancers were later classified as responsive, while the other 21 were nonresponsive. Comparing these two groups, we could show that gains of chromosomal regions 7q32 approximately q36 and 7q11 approximately q31, as well as amplifications of 20q11 approximately q13, were significantly associated with responsiveness to preoperative chemoradiotherapy (P<0.05). However, the probability of detecting these copy number changes by chance is high (P=0.21). Our primary results suggest that pretherapeutic evaluation of chromosomal copy number changes may be of value for response prediction of rectal cancers to preoperative chemoradiotherapy. This will require validation in a larger cohort of patients.

Preoperative chemoradiotherapy in locally advanced rectal cancer: correlation of a gene expression-based response signature with recurrence

Preoperative chemoradiotherapy is recommended for locally advanced rectal cancer (UICC stage II/III). We recently demonstrated that responsive and nonresponsive tumors showed differential expression levels of 54 genes. In this follow-up study, we investigated the relationship between this gene set and disease-free (DFS) and overall survival (OS). Pretherapeutic biopsies from 30 participants in the CAO/ARO/AIO-94 trial of the German Rectal Cancer Study Group were analyzed using gene expression microarrays. Statistical analysis was performed to identify differentially expressed genes between recurrent and nonrecurrent tumors and to correlate these changes with disease recurrence and outcome. After a median follow-up of 59 months, seven of eight patients with recurrent disease was a nonresponder, and one responsive tumor recurred. Response to chemoradiotherapy was significantly correlated with an improved DFS (log rank P=0.028), whereas OS did not differ significantly (P=0.11). Applying a class comparison analysis, we identified 20 genes that were differentially expressed between recurrent and nonrecurrent tumors (P<0.001). Analyzing the first two principal components of the 54 genes previously identified to predict response, we observed that this response signature correlated with an increased risk of cancer recurrence. These data suggest that the genetic basis of local response also affects the genetic basis of tumor recurrence. Genes that are indicative of nonresponse to preoperative chemoradiotherapy might also be linked to an increased risk of tumor recurrence.

A computational procedure to identify significant overlap of differentially expressed and genomic imbalanced regions in cancer datasets

The integration of high-throughput genomic data represents an opportunity for deciphering the interplay between structural and functional organization of genomes and for discovering novel biomarkers. However, the development of integrative approaches to complement gene expression (GE) data with other types of gene information, such as copy number (CN) and chromosomal localization, still represents a computational challenge in the genomic arena. This work presents a computational procedure that directly integrates CN and GE profiles at genome-wide level. When applied to DNA/RNA paired data, this approach leads to the identification of Significant Overlaps of Differentially Expressed and Genomic Imbalanced Regions (SODEGIR). This goal is accomplished in three steps. The first step extends to CN a method for detecting regional imbalances in GE. The second part provides the integration of CN and GE data and identifies chromosomal regions with concordantly altered genomic and transcriptional status in a tumor sample. The last step elevates the single-sample analysis to an entire dataset of tumor specimens. When applied to study chromosomal aberrations in a collection of astrocytoma and renal carcinoma samples, the procedure proved to be effective in identifying discrete chromosomal regions of coordinated CN alterations and changes in transcriptional levels.

Integrating chromosomal aberrations and gene expression profiles to dissect rectal tumorigenesis

Background

Accurate staging of rectal tumors is essential for making the correct treatment choice. In a previous study, we found that loss of 17p, 18q and gain of 8q, 13q and 20q could distinguish adenoma from carcinoma tissue and that gain of 1q was related to lymph node metastasis. In order to find markers for tumor staging, we searched for candidate genes on these specific chromosomes.

Methods

We performed gene expression microarray analysis on 79 rectal tumors and integrated these data with genomic data from the same sample series. We performed supervised analysis to find candidate genes on affected chromosomes and validated the results with qRT-PCR and immunohistochemistry.

Results

Integration of gene expression and chromosomal instability data revealed similarity between these two data types. Supervised analysis identified up-regulation of EFNA1 in cases with 1q gain, and EFNA1 expression was correlated with the expression of a target gene (VEGF). The BOP1 gene, involved in ribosome biogenesis and related to chromosomal instability, was over-expressed in cases with 8q gain. SMAD2 was the most down-regulated gene on 18q, and on 20q, STMN3 and TGIF2 were highly up-regulated. Immunohistochemistry for SMAD4 correlated with SMAD2 gene expression and 18q loss.

Conclusion

On basis of integrative analysis this study identified one well known CRC gene (SMAD2) and several other genes (EFNA1, BOP1, TGIF2 and STMN3) that possibly could be used for rectal cancer characterization.

Chromosomal breakpoints in primary colon cancer cluster at sites of structural variants in the genome

Genomic aberrations on chromosome 8 are common in colon cancer, and are associated with lymph node and distant metastases as well as with disease susceptibility. This prompted us to generate a high-resolution map of genomic imbalances of chromosome 8 in 51 primary colon carcinomas using a custom-designed genomic array consisting of a tiling path of BAC clones. This analysis confirmed the dominant role of this chromosome. Unexpectedly, the position of the breakpoints suggested colocalization with structural variants in the human genome. In order to map these sites with increased resolution and to extend the analysis to the entire genome, we analyzed a subset of these tumors (n = 32) by comparative genomic hybridization on a 185K oligonucleotide array platform. Our comprehensive map of the colon cancer genome confirmed recurrent and specific low-level copy number changes of chromosomes 7, 8, 13, 18, and 20, and unveiled additional, novel sites of genomic imbalances including amplification of a histone gene cluster on chromosome 6p21.1-21.33 and deletions on chromosome 4q34-35. The systematic comparison of segments of copy number change with gene expression profiles showed that genomic imbalances directly affect average expression levels. Strikingly, we observed a significant association of chromosomal breakpoints with structural variants in the human genome: 41% of all copy number changes occurred at sites of such copy number variants (P < 2.2e(-16)). Such an association has not been previously described and reveals a yet underappreciated plasticity of the colon cancer genome; it also points to potential mechanisms for the induction of chromosomal breakage in cancer cells.

Beta-catenin expression results in p53-independent DNA damage and oncogene-induced senescence in prelymphomagenic thymocytes in vivo

The expression of beta-catenin, a potent oncogene, is causally linked to tumorigenesis. Therefore, it was surprising that the transgenic expression of oncogenic beta-catenin in thymocytes resulted in thymic involution instead of lymphomagenesis. In this report, we demonstrate that this is because the expression of oncogenic beta-catenin induces DNA damage, growth arrest, oncogene-induced senescence (OIS), and apoptosis of immature thymocytes. In p53-deficient mice, the expression of oncogenic beta-catenin still results in DNA damage and OIS, but the thymocytes survive and eventually progress to thymic lymphoma. beta-Catenin-induced thymic lymphomas are distinct from lymphomas that arise in p53(-/-) mice. They are CD4(-) CD8(-), while p53-dependent lymphomas are largely CD4(+) CD8(+), and they develop at an earlier age and in the absence of c-Myc expression or Notch1 signaling. Thus, we report that oncogenic beta-catenin-induced, p53-independent growth arrest and OIS and p53-dependent apoptosis protect developing thymocytes from transformation by oncogenic beta-catenin.

Isolation and solubilization of proteins after TRIzol extraction of RNA and DNA from patient material following prolonged storage

A systems approach is being applied in many areas of the biological sciences, particularly in cancer research. The coordinated, simultaneous extraction of DNA, RNA, and proteins from a single sample is crucial for accurate correlations between genomic aberrations and their consequences on the transcriptome and proteome. We present an approach to extract and completely solubilize up to 98% of the total protein recovered from archived samples following TRIzoL isolation of RNA and DNA. We also demonstrate using polyacrylamide gel electrophoresis (PAGE) and Western blot analysis that the proteins, representing both a wide molecular weight range and some posttranslational modifications, such as protein phosphorylation, remain stable in phenol-ethanol for up to 3 years at -20 degrees C.

Single nucleotide polymorphism array analysis of chromosomal instability patterns discriminates rectal adenomas from carcinomas

Total mesorectal excision (TME) is the standard treatment for rectal cancer, while transanal endoscopic microsurgery (TEM) is a recently introduced surgical approach for the treatment of rectal adenomas. Incorrect preoperative staging before TEM is a problem. To identify genetic changes that might correlate with tumour stage and could lead to optimized treatment selection we performed a genome-wide chromosomal instability search in a homogeneous, clinical cohort of rectal tumours. 78 rectal tumours during different clinical stages were analysed with 10K single nucleotide polymorphism (SNP) arrays. Logistic regression was performed to build a quantitative model of specific chromosomal aberrations. Overall, most cases (95%) had one or more chromosomal aberrations. We observed a clear correlation between the total number of aberrations and the different tumour stages. Specifically, the chromosomal events: gain of 8q22-24, 13q and 20q, and loss of 17p and 18q12-22, were far more abundant in carcinoma than in adenoma. In adenoma fractions from cases with a carcinoma (infiltrating at least in the submucosa), twice the amount of such 'malignant aberrations' was observed, compared to pure adenomas. Furthermore, combined aberrations such as gain of 13q and loss of 18q were only found in adenomatous fractions of carcinomas and not in benign lesions. Based on these five genomic events associated with carcinoma, a clear distinction between adenoma and carcinoma tissue could be made. These data should be validated further in order that they may be used in preoperative staging of rectal tumours.

Gene expression profiling reveals a massive, aneuploidy-dependent transcriptional deregulation and distinct differences between lymph node-negative and lymph node-positive colon carcinomas

To characterize patterns of global transcriptional deregulation in primary colon carcinomas, we did gene expression profiling of 73 tumors [Unio Internationale Contra Cancrum stage II (n = 33) and stage III (n = 40)] using oligonucleotide microarrays. For 30 of the tumors, expression profiles were compared with those from matched normal mucosa samples. We identified a set of 1,950 genes with highly significant deregulation between tumors and mucosa samples (P < 1e-7). A significant proportion of these genes mapped to chromosome 20 (P = 0.01). Seventeen genes had a >5-fold average expression difference between normal colon mucosa and carcinomas, including up-regulation of MYC and of HMGA1, a putative oncogene. Furthermore, we identified 68 genes that were significantly differentially expressed between lymph node-negative and lymph node-positive tumors (P < 0.001), the functional annotation of which revealed a preponderance of genes that play a role in cellular immune response and surveillance. The microarray-derived gene expression levels of 20 deregulated genes were validated using quantitative real-time reverse transcription-PCR in >40 tumor and normal mucosa samples with good concordance between the techniques. Finally, we established a relationship between specific genomic imbalances, which were mapped for 32 of the analyzed colon tumors by comparative genomic hybridization, and alterations of global transcriptional activity. Previously, we had conducted a similar analysis of primary rectal carcinomas. The systematic comparison of colon and rectal carcinomas revealed a significant overlap of genomic imbalances and transcriptional deregulation, including activation of the Wnt/beta-catenin signaling cascade, suggesting similar pathogenic pathways.

Stage-specific alterations of the genome, transcriptome, and proteome during colorectal carcinogenesis

To identify sequential alterations of the genome, transcriptome, and proteome during colorectal cancer progression, we have analyzed tissue samples from 36 patients, including the complete mucosa-adenoma-carcinoma sequence from 8 patients. Comparative genomic hybridization (CGH) revealed patterns of stage specific, recurrent genomic imbalances. Gene expression analysis on 9K cDNA arrays identified 58 genes differentially expressed between normal mucosa and adenoma, 116 genes between adenoma and carcinoma, and 158 genes between primary carcinoma and liver metastasis (P < 0.001). Parallel analysis of our samples by CGH and expression profiling revealed a direct correlation of chromosomal copy number changes with chromosome-specific average gene expression levels. Protein expression was analyzed by two-dimensional gel electrophoresis and subsequent mass spectrometry. Although there was no direct match of differentially expressed proteins and genes, the majority of them belonged to identical pathways or networks. In conclusion, increasing genomic instability and a recurrent pattern of chromosomal imbalances as well as specific gene and protein expression changes correlate with distinct stages of colorectal cancer progression. Chromosomal aneuploidies directly affect average resident gene expression levels, thereby contributing to a massive deregulation of the cellular transcriptome. The identification of novel genes and proteins might deliver molecular targets for diagnostic and therapeutic interventions.

Gene expression in 16q is associated with survival and differs between Sørlie breast cancer subtypes

We have investigated the relationship between gene expression and chromosomal positions in 402 breast cancer patients. Using an overrepresentation approach based on Fisher's exact test, we identified disproportionate contributions of specific chromosomal positions to genes associated with survival. Our major finding is that the gene expression in the long arm of chromosome 16 stands out in its relationship to survival. This arm contributes 36 (18%) and 55 (11%) genes to lists negatively associated with recurrence-free survival (set to sizes 200 and 500). This is a highly disproportionate contribution from the 313 (2%) genes in this arm represented on the used Affymetrix U133A and B microarray platforms (Bonferroni corrected Fisher test: P < 2.2 x 10(-16)). We also demonstrate differential expression in 16q across tumor subtypes, which suggests that the ERBB2, basal, and luminal B tumors progress along a high grade-poor prognosis path, while luminal A and normal-like tumors progress along a low grade-good prognosis path, in accordance with a previously proposed model of tumor progression. We conclude that important biological information can be extracted from gene expression data in breast cancer by studying non-random connections between chromosomal positions and gene expression. This article contains Supplementary Material available at http://www.interscience.wiley.com/jpages/1045-2257/suppmat.