Epidemiological studies in the information and genomics era: experience of the Clinical Genome of Cancer Project in São Paulo, Brazil

Genomics is expanding the horizons of epidemiology, providing a new dimension for classical epidemiological studies and inspiring the development of large-scale multicenter studies with the statistical power necessary for the assessment of gene-gene and gene-environment interactions in cancer etiology and prognosis. This paper describes the methodology of the Clinical Genome of Cancer Project in São Paulo, Brazil (CGCP), which includes patients with nine types of tumors and controls. Three major epidemiological designs were used to reach specific objectives: cross-sectional studies to examine gene expression, case-control studies to evaluate etiological factors, and follow-up studies to analyze genetic profiles in prognosis. The clinical groups included patients' data in the electronic database through the Internet. Two approaches were used for data quality control: continuous data evaluation and data entry consistency. A total of 1749 cases and 1509 controls were entered into the CGCP database from the first trimester of 2002 to the end of 2004. Continuous evaluation showed that, for all tumors taken together, only 0.5% of the general form fields still included potential inconsistencies by the end of 2004. Regarding data entry consistency, the highest percentage of errors (11.8%) was observed for the follow-up form, followed by 6.7% for the clinical form, 4.0% for the general form, and only 1.1% for the pathology form. Good data quality is required for their transformation into useful information for clinical application and for preventive measures. The use of the Internet for communication among researchers and for data entry is perhaps the most innovative feature of the CGCP. The monitoring of patients' data guaranteed their quality.

Gene expression profiling of clinical stages II and III breast cancer

Clinical stage (CS) is an established indicator of breast cancer outcome. In the present study, a cDNA microarray platform containing 692 genes was used to identify molecular differences between CSII and CSIII disease. Tumor samples were collected from patients with CSII or CSIII breast cancer, and normal breast tissue was collected from women without invasive cancer. Seventy-eight genes were deregulated in CSIII tumors and 22 in CSII tumors when compared to normal tissue, and 20 of them were differentially expressed in both CSII and CSIII tumors. In addition, 58 genes were specifically altered in CSIII and expression of 6 of them was tested by real time RT-PCR in another cohort of patients with CSII or CSIII breast cancer and in women without cancer. Among these genes, MAX, KRT15 and S100A14, but not APOBEC3G or KRT19, were differentially expressed on both CSIII and CSII tumors as compared to normal tissue. Increased HMOX1 levels were detected only in CSIII tumors and may represent a molecular marker of this stage. A clear difference in gene expression pattern occurs at the normal-to-cancer transition; however, most of the differentially expressed genes are deregulated in tumors of both CS (II and III) compared to normal breast tissue.

CYP1A1 and GSTP1 polymorphisms in an oral cancer case-control study

CYP1A1 and GSTP1 polymorphisms have been associated with a higher risk to develop several cancers, including oral squamous cell carcinoma (OSCC), which is closely related to tobacco and alcohol consumption. Both genes code for enzymes that have an important role in activating or detoxifying carcinogenic elements found in tobacco and other compounds, and polymorphic variants of these genes may result in alterations of the enzymatic activity. The CYP1A1 gene codes for the enzyme aryl hydrocarbon hydroxylase, which is responsible for the metabolism of polycyclic aromatic hydrocarbons. The investigated polymorphism, Ile/Val, seems to increase the activity of the enzyme in homozygous individuals, leading to an accumulation of carcinogens. The Ile/Val polymorphism occurs because of an A->G transition at exon 7, resulting in the CYP1A1*2B allele. The GSTP1*B variant shows an A->G transition at exon 5, changing the amino acid Ile to Val, with a reduced catalytic activity of the enzyme. Due to this reduction, the carriers of mutant alleles lost the capability to metabolize carcinogens, which could be responsible for a higher susceptibility to cancer. We conducted a case-control study in a group of 72 cases with newly diagnosed OSCC and 60 healthy controls matched for age, gender, smoking habits, and ethnicity. We used PCR methods to identify the allelic variants CYP1A1*2B and GSTP1*B. The data obtained showed no statistically significant association of allelic or genotypic variants of CYP1A1*2B (OR = 1.06; 95% CI = 0.49-2.29) and GSTP1*B (OR = 1.40; 95% CI = 0.70-2.79) with OSCC.

Pipeline for macro- and microarray analyses

The pipeline for macro- and microarray analyses (PMmA) is a set of scripts with a web interface developed to analyze DNA array data generated by array image quantification software. PMmA is designed for use with single- or double-color array data and to work as a pipeline in five classes (data format, normalization, data analysis, clustering, and array maps). It can also be used as a plugin in the BioArray Software Environment, an open-source database for array analysis, or used in a local version of the web service. All scripts in PMmA were developed in the PERL programming language and statistical analysis functions were implemented in the R statistical language. Consequently, our package is a platform-independent software. Our algorithms can correctly select almost 90% of the differentially expressed genes, showing a superior performance compared to other methods of analysis. The pipeline software has been applied to 1536 expressed sequence tags macroarray public data of sugarcane exposed to cold for 3 to 48 h. PMmA identified thirty cold-responsive genes previously unidentified in this public dataset. Fourteen genes were up-regulated, two had a variable expression and the other fourteen were down-regulated in the treatments. These new findings certainly were a consequence of using a superior statistical analysis approach, since the original study did not take into account the dependence of data variability on the average signal intensity of each gene. The web interface, supplementary information, and the package source code are available, free, to non-commercial users at http://ipe.cbmeg.unicamp.br/pub/PMmA.

Comparison of blood neoangiogenesis and lymphatic vascularization in colorectal adenomas from patients with and without concomitant colorectal cancer

Blood and lymphatic vessel proliferation is essential for tumor growth and progression. Most colorectal carcinomas develop from adenomas (adenoma-carcinoma sequence) in a process due to accumulation of molecular genetic alterations. About 5% of adenomatous polyps are expected to become malignant, but data on the differential angiogenic patterns of these lesions in patients with and without concomitant cancer are missing. The aim of the present study is to compare the angiogenic and lymphatic patterns of adenomatous polyps from patients with and without sporadic cancer. Thirty adenomatous polyps (15 from patients with another principal malignant lesion, and 15 from patients without cancer) were submitted to immunohistochemical staining for CD105 (marker for neoangiogenesis) and D2-40 (marker for lymphatic endothelium). Microvessel density and total vascular area were determined by computer image analysis to quantify the immunostained and total areas, and to assess the number of microvessels. Adenomas from patients with carcinoma showed significantly higher values of total vascular area determined by immunostaining for CD105 (cutoff value = 4386 µm²; P = 0.019) and of lymphatic microvessel density determined by immunostaining with D2-40 (cutoff value = 11.5; P = 0.041) when compared with those from patients without cancer. The present data indicate a significant increase in blood microvascular area and in lymphatic microvascular counts in adenomas removed from patients with cancer.

Growth inhibitory effect and Chk1-dependent signaling involved in G2/M arrest on human gastric cancer cells induced by diallyl disulfide

Diallyl disulfide (DADS) inhibits growth and induces cell cycle G2/M arrest in human gastric cancer MGC803 cells. In this study, 15 mg/L DADS exerted similar effects on growth and cell cycle arrest in human gastric cancer BGC823 cells. Due to the importance of cell cycle redistribution in DADS-mediated anti-carcinogenic effects, we investigated the role of checkpoint kinases (Chk1 and Chk2) during DADS-induced cell cycle arrest. We hypothesized that DADS could mediate G2/M phase arrest through either Chk1 or Chk2 signal transduction pathways. We demonstrated that DADS induced the accumulation of phosphorylated Chk1, but not of Chk2, and that DADS down-regulated Cdc25C and cyclin B1. The expression of mRNA and total protein for Chkl and Chk2 was unchanged. Chk1 is specifically phosphorylated by ATR (ATM-RAD3-related gene). Western blot analysis showed that phospho-ATR was activated by DADS. Taken together, these data suggest that cell cycle G2/M arrest, which was associated with accumulation of the phosphorylated forms of Chk1, but not of Chk2, was involved in the growth inhibition induced by DADS in the human gastric cancer cell line BGC823. Furthermore, the DADS-induced G2/M checkpoint response is mediated by Chk1 signaling through ATR/Chk1/Cdc25C/cyclin B1, and is independent of Chk2.

5-Methyltetrahydrofolate-homocysteine methyltransferase gene polymorphism (MTR) and risk of head and neck cancer

The functional effect of the A>G transition at position 2756 on the MTR gene (5-methyltetrahydrofolate-homocysteine methyltransferase), involved in folate metabolism, may be a risk factor for head and neck squamous cell carcinoma (HNSCC). The frequency of MTR A2756G (rs1805087) polymorphism was compared between HNSCC patients and individuals without history of neoplasias. The association of this polymorphism with clinical histopathological parameters was evaluated. A total of 705 individuals were included in the study. The polymerase chain reaction-restriction fragment length polymorphism technique was used to genotype the polymorphism. For statistical analysis, the chi-square test (univariate analysis) was used for comparisons between groups and multiple logistic regression (multivariate analysis) was used for interactions between the polymorphism and risk factors and clinical histopathological parameters. Using univariate analysis, the results did not show significant differences in allelic or genotypic distributions. Multivariable analysis showed that tobacco and alcohol consumption (P < 0.05), AG genotype (P = 0.019) and G allele (P = 0.028) may be predictors of the disease and a higher frequency of the G polymorphic allele was detected in men with HNSCC compared to male controls (P = 0.008). The analysis of polymorphism regarding clinical histopathological parameters did not show any association with the primary site, aggressiveness, lymph node involvement or extension of the tumor. In conclusion, our data provide evidence that supports an association between the polymorphism and the risk of HNSCC.

In vitro modulation of Bcl-2 levels in small cell lung cancer cells: effects on cell viability

Small cell lung cancer (SCLC) is an aggressive disease, representing 15% of all cases of lung cancer, has high metastatic potential and low prognosis that urgently demands the development of novel therapeutic approaches. One of the proposed approaches has been the down-regulation of BCL2, with poorly clarified and controversial therapeutic value regarding SCLC. The use of anti-BCL2 small interfering RNA (siRNA) in SCLC has never been reported. The aim of the present study was to select and test the in vitro efficacy of anti-BCL2 siRNA sequences against the protein and mRNA levels of SCLC cells, and their effects on cytotoxicity and chemosensitization. Two anti-BCL2 siRNAs and the anti-BCL2 G3139 oligodeoxynucleotide (ODN) were evaluated in SCLC cells by the simultaneous determination of Bcl-2 and viability using a flow cytometry method recently developed by us in addition to Western blot, real-time reverse-transcription PCR, and cell growth after single and combined treatment with cisplatin. In contrast to previous reports about the use of ODN, a heterogeneous and up to 80% sequence-specific Bcl-2 protein knockdown was observed in the SW2, H2171 and H69 SCLC cell lines, although without significant sequence-specific reduction of cell viability, cell growth, or sensitization to cisplatin. Our results question previous data generated with antisense ODN and supporting the present concept of the therapeutic interest in BCL2 silencing per se in SCLC, and support the growing notion of the necessity of a multitargeting molecular approach for the treatment of cancer.

Prognosis of colorectal cancer with liver metastasis: value of a prognostic index

The objective of the present study was to explore the factors related to the prognosis of colorectal cancer (CRC) and to establish a prognostic model for the selection of patients who might benefit from hepatic resection for metastatic CRC. A total of 293 patients undergoing liver resection for metastatic CRC (172 males and 80 females ranging in age from 26 to 80 years) were selected and clinical, pathological and outcome data were examined in this retrospective study. The prognostic index (PI) of the patients was calculated on the basis of results of multivariate analysis. Patients were stratified into different groups, with survival curves projected according to PI. The 1-, 3-, and 5-year overall survival rates were 58.3, 26.4, and 11.3%, respectively. Univariate analysis indicated that degree of primary tumor differentiation, resection margin, preoperative carcinoembryonic antigen (CEA) level, number of liver metastases, and resection of liver metastases were associated with prognosis (P < 0.05). In multivariate analysis, the last three factors were found to be independent prognostic factors. The resection of liver metastases was a favorable factor. Patients were classified into three groups according to PI, which differed significantly in survival rate (P < 0.05). The individual survival rate was evaluated based on PI. Resection of hepatic colorectal metastases may produce long-term survival and cure. The proposed PI was easy to use, was highly predictive of patient outcome, and permitted categorization of patients into treatment groups.

Expression of TERT in precancerous gastric lesions compared to gastric cancer

The objective of this study was to determine the levels of TERT mRNA and TERT protein expression in stomach precancerous lesions such as intestinal metaplasia (IM) and gastric ulcer (GU) and compare them to gastric cancer (GC). Real-time PCR was performed to detect TERT mRNA expression levels in 35 biopsies of IM, 30 of GU, and 22 of GC and their respective normal mucosas. TERT protein was detected by immunohistochemistry in 68 samples, 34 of IM, 23 of GU, and 11 of GC. Increased TERT mRNA expression levels were observed in a significant number of cases, i.e., 46% of IM, 50% of GU, and 79% of GC. The relative mean level of TERT mRNA after normalization with the β-actin reference gene and comparison with the respective adjacent normal mucosa was slightly increased in the IM and GU groups, 2.008 ± 2.605 and 2.730 ± 4.120, respectively, but high TERT mRNA expression was observed in the GC group (17.271 ± 33.852). However, there were no statistically significant differences between the three groups. TERT protein-positive immunostaining was observed in 38% of IM, 39% of GU, and 55% of GC. No association of TERT mRNA and protein expression with Helicobacter pylori infection or other clinicopathological variables was demonstrable, except for the incomplete type vs the complete type of IM. This study confirms previous data of the high expression of both TERT mRNA and protein in gastric cancer and also demonstrates this type of changed expression in IM and GU, thus suggesting that TERT expression may be deregulated in precursor lesions that participate in the early stages of gastric carcinogenesis.

Transcriptional regulation of bidirectional gene pairs by 17-β-estradiol in MCF-7 breast cancer cells

Using cDNA microarray analysis, we previously identified a set of differentially expressed genes in primary breast tumors based on the status of estrogen and progesterone receptors. In the present study, we performed an integrated computer-assisted and manual search of potential estrogen response element (ERE) binding sites in the promoter region of these genes to characterize their potential to be regulated by estrogen receptors (ER). Publicly available databases were used to annotate the position of these genes in the genome and to extract a 5’flanking region 2 kb upstream to 2 kb downstream of the transcription start site for transcription binding site analysis. The search for EREs and other binding sites was performed using several publicly available programs. Overall, approximately 40% of the genes analyzed were potentially able to be regulated by estrogen via ER. In addition, 17% of these genes are located very close to other genes organized in a head-to-head orientation with less than 1.0 kb between their transcript units, sharing a bidirectional promoter, and could be classified as bidirectional gene pairs. Using quantitative real-time PCR, we further investigated the effects of 17β-estradiol and antiestrogens on the expression of the bidirectional gene pairs in MCF-7 breast cancer cells. Our results showed that some of these gene pairs, such as TXNDC9/EIF5B, GALNS/TRAPPC2L, and SERINC1/PKIB, are modulated by 17β-estradiol via ER in MCF-7 breast cancer cells. Here, we also characterize the promoter region of potential ER-regulated genes and provide new information on the transcriptional regulation of bidirectional gene pairs.

Serum metabolic profiling of human gastric cancer based on gas chromatography/mass spectrometry

Research on molecular mechanisms of carcinogenesis plays an important role in diagnosing and treating gastric cancer. Metabolic profiling may offer the opportunity to understand the molecular mechanism of carcinogenesis and help to non-invasively identify the potential biomarkers for the early diagnosis of human gastric cancer. The aims of this study were to explore the underlying metabolic mechanisms of gastric cancer and to identify biomarkers associated with morbidity. Gas chromatography/mass spectrometry (GC/MS) was used to analyze the serum metabolites of 30 Chinese gastric cancer patients and 30 healthy controls. Diagnostic models for gastric cancer were constructed using orthogonal partial least squares discriminant analysis (OPLS-DA). Acquired metabolomic data were analyzed by the nonparametric Wilcoxon test to find serum metabolic biomarkers for gastric cancer. The OPLS-DA model showed adequate discrimination between cancer and non-cancer cohorts while the model failed to discriminate different pathological stages (I-IV) of gastric cancer patients. A total of 44 endogenous metabolites such as amino acids, organic acids, carbohydrates, fatty acids, and steroids were detected, of which 18 differential metabolites were identified with significant differences. A total of 13 variables were obtained for their greatest contribution in the discriminating OPLS-DA model [variable importance in the projection (VIP) value >1.0], among which 11 metabolites were identified using both VIP values (VIP >1) and the Wilcoxon test. These metabolites potentially revealed perturbations of glycolysis and of amino acid, fatty acid, cholesterol, and nucleotide metabolism of gastric cancer patients. These results suggest that gastric cancer serum metabolic profiling has great potential in detecting this disease and helping to understand its metabolic mechanisms.

Efficacy of the dietary histone deacetylase inhibitor butyrate alone or in combination with vitamin A against proliferation of MCF-7 human breast cancer cells

The combined treatment with histone deacetylase inhibitors (HDACi) and retinoids has been suggested as a potential epigenetic strategy for the control of cancer. In the present study, we investigated the effects of treatment with butyrate, a dietary HDACi, combined with vitamin A on MCF-7 human breast cancer cells. Cell proliferation was evaluated by the crystal violet staining method. MCF-7 cells were plated at 5 x 10(4) cells/mL and treated with butyrate (1 mM) alone or combined with vitamin A (10 µM) for 24 to 120 h. Cell proliferation inhibition was 34, 10 and 46% following treatment with butyrate, vitamin A and their combination, respectively, suggesting that vitamin A potentiated the inhibitory activities of butyrate. Furthermore, exposure to this short-chain fatty acid increased the level of histone H3K9 acetylation by 9.5-fold (Western blot), but not of H4K16, and increased the expression levels of p21WAF1 by 2.7-fold (Western blot) and of RARβ by 2.0-fold (quantitative real-time PCR). Our data show that RARβ may represent a molecular target for butyrate in breast cancer cells. Due to its effectiveness as a dietary HDACi, butyrate should be considered for use in combinatorial strategies with more active retinoids, especially in breast cancers in which RARβ is epigenetically altered.

Leptin receptor (Ob-R) mRNA expression and serum leptin concentration in patients with colorectal and metastatic colorectal cancer

The objective of the present study was to investigate the effect of leptin on the progression of colorectal carcinoma to metastatic disease by analyzing the serum leptin concentration and Ob-R gene expression in colon cancer tissues. Tissue samples were obtained from 31 patients who underwent surgical resection for colon (18 cases) and metastatic colon (13 cases) cancer. Serum leptin concentration was determined by an enzyme-linked immunosorbent assay (ELISA) and Ob-R mRNA expression by real-time polymerase chain reaction (RT-PCR) for both groups. ELISA data were analyzed by the Student t-test and RT-PCR data were analyzed by the Mann-Whitney U-test. RT-PCR results demonstrated that mRNA expression of Ob-R in human metastatic colorectal cancer was higher than in local colorectal cancer tissues. On the other hand, mean serum leptin concentration was significantly higher in local colorectal cancer patients compared to patients with metastatic colorectal cancer. The results of the present study suggest a role for leptin in the progression of colon cancer to metastatic disease without weight loss. In other words, significantly increased Ob-R mRNA expression and decreased serum leptin concentration in patients with metastatic colon cancer indicate that sensitization to leptin activity may be a major indicator of metastasis to the colon tissue and the determination of leptin concentration and leptin gene expression may be used to aid the diagnosis.

AKT inhibitor suppresses hyperthermia-induced Ndrg2 phosphorylation in gastric cancer cells

Hyperthermia is one of the most effective adjuvant treatments for various cancers with few side effects. However, the underlying molecular mechanisms still are not known. N-myc downstream-regulated gene 2 (NDRG2), a tumor suppressor, has been shown to be involved in diverse cellular stresses including hypoxia, lipotoxicity, etc. In addition, Ndrg2 has been reported to be related to progression of gastric cancer. In the current study, our data showed that the apoptosis rate of MKN28 cells increased relatively rapidly to 13.4% by 24 h after treatment with hyperthermia (42°C for 1 h) compared to 5.1% in control cells (P < 0.05). Nevertheless, there was no obvious change in the expression level of total Ndrg2 during this process. Further investigation demonstrated that the relative phosphorylation levels of Ndrg2 at Ser332, Thr348 increased up to 3.2- and 1.9-fold (hyperthermia groupvs control group) at 3 h in MKN28 cells, respectively (P < 0.05). We also found that heat treatment significantly increased AKT phosphorylation. AKT inhibitor VIII (10 µM) decreased the phosphorylation level of Ndrg2 induced by hyperthermia. Accordingly, the apoptosis rate rose significantly in MKN28 cells (16.4%) treated with a combination of AKT inhibitor VIII and hyperthermia compared to that (6.8%) of cells treated with hyperthermia alone (P < 0.05). Taken together, these data demonstrated that Ndrg2 phosphorylation could be induced by hyperthermia in an AKT-dependent manner in gastric cancer cells. Furthermore, AKT inhibitor VIII suppressed Ndrg2 phosphorylation and rendered gastric cancer cells susceptible to apoptosis induced by hyperthermia.