Trends in socioeconomic inequalities in cancer mortality in Barcelona: 1992–2003

Abstract Background The objective of this study was to assess trends in cancer mortality by educational level in Barcelona from 1992 to 2003. Methods The study population comprised Barcelona inhabitants aged 20 years or older. Data on cancer deaths were supplied by the system of information on mortality. Educational level was obtained from the municipal census. Age-standardized rates by educational level were calculated. We also fitted Poisson regression models to estimate the relative index of inequality (RII) and the Slope Index of Inequalities (SII). All were calculated for each sex and period (1992–1994, 1995–1997, 1998–2000, and 2001–2003). Results Cancer mortality was higher in men and women with lower educational level throughout the study period. Less-schooled men had higher mortality by stomach, mouth and pharynx, oesophagus, larynx and lung cancer. In women, there were educational inequalities for cervix uteri, liver and colon cancer. Inequalities of overall and specific types of cancer mortality remained stable in Barcelona; although a slight reduction was observed for some cancers. Conclusion This study has identified those cancer types presenting the greatest inequalities between men and women in recent years and shown that in Barcelona there is a stable trend in inequalities in the burden of cancer.

Inequalities in mortality of men by oral and pharyngeal cancer in Barcelona, Spain and São Paulo, Brazil, 1995–2003

Abstract Background Large inequalities of mortality by most cancers in general, by mouth and pharynx cancer in particular, have been associated to behaviour and geopolitical factors. The assessment of socioeconomic covariates of cancer mortality may be relevant to a full comprehension of distal determinants of the disease, and to appraise opportune interventions. The objective of this study was to compare socioeconomic inequalities in male mortality by oral and pharyngeal cancer in two major cities of Europe and South America. Methods The official system of information on mortality provided data on deaths in each city; general censuses informed population data. Age-adjusted death rates by oral and pharyngeal cancer for men were independently assessed for neighbourhoods of Barcelona, Spain, and São Paulo, Brazil, from 1995 to 2003. Uniform methodological criteria instructed the comparative assessment of magnitude, trends and spatial distribution of mortality. General linear models assessed ecologic correlations between death rates and socioeconomic indices (unemployment, schooling levels and the human development index) at the inner-city area level. Results obtained for each city were subsequently compared. Results Mortality of men by oral and pharyngeal cancer ranked higher in Barcelona (9.45 yearly deaths per 100,000 male inhabitants) than in Spain and Europe as a whole; rates were on decrease. São Paulo presented a poorer profile, with higher magnitude (11.86) and stationary trend. The appraisal of ecologic correlations indicated an unequal and inequitably distributed burden of disease in both cities, with poorer areas tending to present higher mortality. Barcelona had a larger gradient of mortality than São Paulo, indicating a higher inequality of cancer deaths across its neighbourhoods. Conclusion The quantitative monitoring of inequalities in health may contribute to the formulation of redistributive policies aimed at the concurrent promotion of wellbeing and social justice. The assessment of groups experiencing a higher burden of disease can instruct health services to provide additional resources for expanding preventive actions and facilities aimed at early diagnosis, standardized treatments and rehabilitation.

Transcriptional effects of 1,25 dihydroxyvitamin D3 physiological and supra-physiological concentrations in breast cancer organotypic culture

Background Vitamin D transcriptional effects were linked to tumor growth control, however, the hormone targets were determined in cell cultures exposed to supra physiological concentrations of 1,25(OH)2D3 (50-100nM). Our aim was to evaluate the transcriptional effects of 1,25(OH)2D3 in a more physiological model of breast cancer, consisting of fresh tumor slices exposed to 1,25(OH)2D3 at concentrations that can be attained in vivo. Methods Tumor samples from post-menopausal breast cancer patients were sliced and cultured for 24 hours with or without 1,25(OH)2D3 0.5nM or 100nM. Gene expression was analyzed by microarray (SAM paired analysis, FDR≤0.1) or RT-qPCR (p≤0.05, Friedman/Wilcoxon test). Expression of candidate genes was then evaluated in mammary epithelial/breast cancer lineages and cancer associated fibroblasts (CAFs), exposed or not to 1,25(OH)2D3 0.5nM, using RT-qPCR, western blot or immunocytochemistry. Results 1,25(OH)2D3 0.5nM or 100nM effects were evaluated in five tumor samples by microarray and seven and 136 genes, respectively, were up-regulated. There was an enrichment of genes containing transcription factor binding sites for the vitamin D receptor (VDR) in samples exposed to 1,25(OH)2D3 near physiological concentration. Genes up-modulated by both 1,25(OH)2D3 concentrations were CYP24A1, DPP4, CA2, EFTUD1, TKTL1, KCNK3. Expression of candidate genes was subsequently evaluated in another 16 samples by RT-qPCR and up-regulation of CYP24A1, DPP4 and CA2 by 1,25(OH)2D3 was confirmed. To evaluate whether the transcripitonal targets of 1,25(OH)2D3 0.5nM were restricted to the epithelial or stromal compartments, gene expression was examined in HB4A, C5.4, SKBR3, MDA-MB231, MCF-7 lineages and CAFs, using RT-qPCR. In epithelial cells, there was a clear induction of CYP24A1, CA2, CD14 and IL1RL1. In fibroblasts, in addition to CYP24A1 induction, there was a trend towards up-regulation of CA2, IL1RL1, and DPP4. A higher protein expression of CD14 in epithelial cells and CA2 and DPP4 in CAFs exposed to 1,25(OH)2D3 0.5nM was detected. Conclusions In breast cancer specimens a short period of 1,25(OH)2D3 exposure at near physiological concentration modestly activates the hormone transcriptional pathway. Induction of CYP24A1, CA2, DPP4, IL1RL1 expression appears to reflect 1,25(OH)2D3 effects in epithelial as well as stromal cells, however, induction of CD14 expression is likely restricted to the epithelial compartment.

Independent of ErbB1 gene copy number, EGF stimulates migration but is not associated with cell proliferation in non-small cell lung cancer

Abstract Background Lung cancer often exhibits molecular changes, such as the overexpression of the ErbB1 gene. ErbB1 encodes epidermal growth factor receptor (EGFR), a tyrosine kinase receptor, involved mainly in cell proliferation and survival. EGFR overexpression has been associated with more aggressive disease, poor prognosis, low survival rate and low response to therapy. ErbB1 amplification and mutation are associated with tumor development and are implicated in ineffective treatment. The aim of the present study was to investigate whether the ErbB1 copy number affects EGFR expression, cell proliferation or cell migration by comparing two different cell lines. Methods The copies of ErbB1 gene was evaluated by FISH. Immunofluorescence and Western blotting were performed to determine location and expression of proteins mentioned in the present study. Proliferation was studied by flow cytometry and cell migration by wound healing assay and time lapse. Results We investigated the activation and function of EGFR in the A549 and HK2 lung cancer cell lines, which contain 3 and 6 copies of ErbB1, respectively. The expression of EGFR was lower in the HK2 cell line. EGFR was activated after stimulation with EGF in both cell lines, but this activation did not promote differences in cellular proliferation when compared to control cells. Inhibiting EGFR with AG1478 did not modify cellular proliferation, confirming previous data. However, we observed morphological alterations, changes in microfilament organization and increased cell migration upon EGF stimulation. However, these effects did not seem to be consequence of an epithelial-mesenchymal transition. Conclusion EGFR expression did not appear to be associated to the ErbB1 gene copy number, and neither of these aspects appeared to affect cell proliferation. However, EGFR activation by EGF resulted in cell migration stimulation in both cell lines.

Prognostic factors for survival in patients with colorectal liver metastases: experience of a single brazilian cancer center

CONTEXT: Liver metastases are a common event in the clinical outcome of patients with colorectal cancer and account for 2/3 of deaths from this disease. There is considerable controversy among the data in the literature regarding the results of surgical treatment and prognostic factors of survival, and no analysis have been done in a large cohort of patients in Brazil. OBJECTIVES: To characterize the results of surgical treatment of patients with colorectal liver metastases, and to establish prognostic factors of survival in a Brazilian population. METHOD: This was a retrospective study of patients undergoing liver resection for colorectal metastases in a tertiary cancer hospital from 1998 to 2009. We analyzed epidemiologic variables and the clinical characteristics of primary tumors, metastatic disease and its treatment, surgical procedures and follow-up, and survival results. Survival analyzes were done by the Kaplan-Meier method and the log-rank test was applied to determine the influence of variables on overall and disease-free survival. All variables associated with survival with P<0.20 in univariate analysis, were included in multivariate analysis using a Cox proportional hazard regression model. RESULTS: During the period analyzed, 209 procedures were performed on 170 patients. Postope-rative mortality in 90 days was 2.9% and 5-year overall survival was 64.9%. Its independent prognostic factors were the presence of extrahepatic disease at diagnosis of liver metastases, bilateral nodules and the occurrence of major complications after liver surgery. The estimated 5-year disease-free survival was 39.1% and its prognostic factors included R1 resection, extrahepatic disease, bilateral nodules, lymph node involvement in the primary tumor and primary tumors located in the rectum. CONCLUSION: Liver resection for colorectal metastases is safe and effective and the analysis of prognostic factors of survival in a large cohort of Brazilian patients showed similar results to those pointed in international series. The occurrence of major postoperative complications appears to be able to compromise overall survival and further investigation in needed in this topic.

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.

Reorganization of secondary and tertiary health care levels: impact on the outcomes of oral cancer screening in the São Paulo State, Brazil

This study presents the strategies for prevention and early detection of oral cancer by means of screening in the elderly population of São Paulo, the richest and the most populous state of Brazil. This research was a retrospective longitudinal study based on the analysis of secondary data. The variables - number of participating cities, coverage of screening, and number of suspicious and confirmed cases of oral cancer - were divided into two periods: 2001-2004 and 2005-2008. Data were analyzed statistically by the chi-square test at 5% significance level. The implementation of a nationwide public oral health policy in 2004 and the reorganization of the secondary and tertiary health care were evaluated as mediator factors able to interfere in the achieved outcomes. From 2001 to 2008, 2,229,273 oral examinations were performed. There was an addition of 205 participating cities by the end of the studied period (p<0.0001). The coverage of oral cancer screening increased from 4.1% to 16% (p<0.0001). There was a decrease in the number of suspicious lesions (from 9% in 2005 to 5% in 2008) (p<0.0001) and in the rate of confirmed oral cancer cases per 100,000 examinations (from 20.89 in 2001 to 10.40 in 2008) (p<0.0001). After 8 years of screening, there was a decrease in the number of suspicious lesions and confirmed cases of oral cancer in the population. The reorganization of secondary and tertiary health care levels of oral care seems to have contributed to modify these numbers, having a positive impact on the outcomes of oral cancer screening in the São Paulo State.

Prevention of skin reactions due to teletherapy in women with breast cancer: a comprehensive review

One of the possible courses of cancer treatment is teletherapy, and one of the most important adverse side effects are skin reactions, an ailment more commonly called radiodermatitis. The main purpose of this study is to analyze knowledge of the evidence about topical products used in the prevention of radiodermatitis, to support care delivery to women with breast cancer during teletherapy. The research method used here is the comprehensive literature review. Four databases were used to select the bibliography. The sample consists of 15 articles. The data shows that, among the topical products analyzed here, Calendula, corticosteroids and Xclair have shown significant protective effects, underlining their actions. The lack of articles published in Brazil highlights the need for further research in this area, seeking better care quality through the use of products with scientifically proven efficiency.

High-throughput sequencing of small RNA transcriptomes reveals critical biological features targeted by microRNAs in cell models used for squamous cell cancer research

Abstract Background The implication of post-transcriptional regulation by microRNAs in molecular mechanisms underlying cancer disease is well documented. However, their interference at the cellular level is not fully explored. Functional in vitro studies are fundamental for the comprehension of their role; nevertheless results are highly dependable on the adopted cellular model. Next generation small RNA transcriptomic sequencing data of a tumor cell line and keratinocytes derived from primary culture was generated in order to characterize the microRNA content of these systems, thus helping in their understanding. Both constitute cell models for functional studies of microRNAs in head and neck squamous cell carcinoma (HNSCC), a smoking-related cancer. Known microRNAs were quantified and analyzed in the context of gene regulation. New microRNAs were investigated using similarity and structural search, ab initio classification, and prediction of the location of mature microRNAs within would-be precursor sequences. Results were compared with small RNA transcriptomic sequences from HNSCC samples in order to access the applicability of these cell models for cancer phenotype comprehension and for novel molecule discovery. Results Ten miRNAs represented over 70% of the mature molecules present in each of the cell types. The most expressed molecules were miR-21, miR-24 and miR-205, Accordingly; miR-21 and miR-205 have been previously shown to play a role in epithelial cell biology. Although miR-21 has been implicated in cancer development, and evaluated as a biomarker in HNSCC progression, no significant expression differences were seen between cell types. We demonstrate that differentially expressed mature miRNAs target cell differentiation and apoptosis related biological processes, indicating that they might represent, with acceptable accuracy, the genetic context from which they derive. Most miRNAs identified in the cancer cell line and in keratinocytes were present in tumor samples and cancer-free samples, respectively, with miR-21, miR-24 and miR-205 still among the most prevalent molecules at all instances. Thirteen miRNA-like structures, containing reads identified by the deep sequencing, were predicted from putative miRNA precursor sequences. Strong evidences suggest that one of them could be a new miRNA. This molecule was mostly expressed in the tumor cell line and HNSCC samples indicating a possible biological function in cancer. Conclusions Critical biological features of cells must be fully understood before they can be chosen as models for functional studies. Expression levels of miRNAs relate to cell type and tissue context. This study provides insights on miRNA content of two cell models used for cancer research. Pathways commonly deregulated in HNSCC might be targeted by most expressed and also by differentially expressed miRNAs. Results indicate that the use of cell models for cancer research demands careful assessment of underlying molecular characteristics for proper data interpretation. Additionally, one new miRNA-like molecule with a potential role in cancer was identified in the cell lines and clinical samples.

MicroRNA expression profile in head and neck cancer: HOX-cluster embedded microRNA-196a and microRNA-10b dysregulation implicated in cell proliferation

Abstract Background Current evidence implicates aberrant microRNA expression patterns in human malignancies; measurement of microRNA expression may have diagnostic and prognostic applications. Roles for microRNAs in head and neck squamous cell carcinomas (HNSCC) are largely unknown. HNSCC, a smoking-related cancer, is one of the most common malignancies worldwide but reliable diagnostic and prognostic markers have not been discovered so far. Some studies have evaluated the potential use of microRNA as biomarkers with clinical application in HNSCC. Methods MicroRNA expression profile of oral squamous cell carcinoma samples was determined by means of DNA microarrays. We also performed gain-of-function assays for two differentially expressed microRNA using two squamous cell carcinoma cell lines and normal oral keratinocytes. The effect of the over-expression of these molecules was evaluated by means of global gene expression profiling and cell proliferation assessment. Results Altered microRNA expression was detected for a total of 72 microRNAs. Among these we found well studied molecules, such as the miR-17-92 cluster, comprising potent oncogenic microRNA, and miR-34, recently found to interact with p53. HOX-cluster embedded miR-196a/b and miR-10b were up- and down-regulated, respectively, in tumor samples. Since validated HOX gene targets for these microRNAs are not consistently deregulated in HNSCC, we performed gain-of-function experiments, in an attempt to outline their possible role. Our results suggest that both molecules interfere in cell proliferation through distinct processes, possibly targeting a small set of genes involved in cell cycle progression. Conclusions Functional data on miRNAs in HNSCC is still scarce. Our data corroborate current literature and brings new insights into the role of microRNAs in HNSCC. We also show that miR-196a and miR-10b, not previously associated with HNSCC, may play an oncogenic role in this disease through the deregulation of cell proliferation. The study of microRNA alterations in HNSCC is an essential step to the mechanistic understanding of tumor formation and could lead to the discovery of clinically relevant biomarkers.

Complex organochlorine pesticide mixtures as determinant factor for breast cancer risk: a population-based case-control study in the Canary Islands (Spain)

[EN] Background: All the relevant risk factors contributing to breast cancer etiology are not fully known. Exposure to organochlorine pesticides has been linked to an increased incidence of the disease, although not all data have been consistent. Most published studies evaluated the exposure to organochlorines individually, ignoring the potential effects exerted by the mixtures of chemicals. Methods: This population-based study was designed to evaluate the profile of mixtures of organochlorines detected in 103 healthy women and 121 women diagnosed with breast cancer from Gran Canaria Island, and the relation between the exposure to these compounds and breast cancer risk.Results: The most prevalent mixture of organochlorines among healthy women was the combination of lindane and endrin, and this mixture was not detected in any affected women. Breast cancer patients presented more frequently a combination of aldrin, dichlorodiphenyldichloroethylene (DDE) and dichlorodiphenyldichloroethane (DDD), and this mixture was not found in any healthy woman. After adjusting for covariables, the risk of breast cancer was moderately associated with DDD (OR = 1.008, confidence interval 95% 1.001-1.015, p = 0.024).Conclusions: This study indicates that healthy women show a very different profile of organochlorine pesticide mixtures than breast cancer patients, suggesting that organochlorine pesticide mixtures could play a relevant role in breast cancer risk.

Radiation induced apoptosis and initial DNA damage are inversely related in locally advanced breast cancer patients

[EN] Background: DNA-damage assays, quantifying the initial number of DNA double-strand breaks induced by radiation, have been proposed as a predictive test for radiation-induced toxicity. Determination of radiation-induced apoptosis in peripheral blood lymphocytes by flow cytometry analysis has also been proposed as an approach for predicting normal tissue responses following radiotherapy. The aim of the present study was to explore the association between initial DNA damage, estimated by the number of double-strand breaks induced by a given radiation dose, and the radio-induced apoptosis rates observed. Methods: Peripheral blood lymphocytes were taken from 26 consecutive patients with locally advanced breast carcinoma. Radiosensitivity of lymphocytes was quantified as the initial number of DNA double-strand breaks induced per Gy and per DNA unit (200 Mbp). Radio-induced apoptosis at 1, 2 and 8 Gy was measured by flow cytometry using annexin V/propidium iodide. Results: Radiation-induced apoptosis increased in order to radiation dose and data fitted to a semi logarithmic mathematical model. A positive correlation was found among radio-induced apoptosis values at different radiation doses: 1, 2 and 8 Gy (p < 0.0001 in all cases). Mean DSB/Gy/DNA unit obtained was 1.70 ± 0.83 (range 0.63-4.08; median, 1.46). A statistically significant inverse correlation was found between initial damage to DNA and radio-induced apoptosis at 1 Gy (p = 0.034). A trend toward 2 Gy (p = 0.057) and 8 Gy (p = 0.067) was observed after 24 hours of incubation. Conclusions: An inverse association was observed for the first time between these variables, both considered as predictive factors to radiation toxicity.

Prediction of clinical toxicity in locally advanced head and neck cancer patients by radio-induced apoptosis in peripheral blood lymphocytes (PBLs)

[EN] Head and neck cancer is treated mainly by surgery and radiotherapy. Normal tissue toxicity due to x-ray exposure is a limiting factor for treatment success. Many efforts have been employed to develop predictive tests applied to clinical practice. Determination of lymphocyte radio-sensitivity by radio-induced apoptosis arises as a possible method to predict tissue toxicity due to radiotherapy. The aim of the present study was to analyze radio-induced apoptosis of peripheral blood lymphocytes in head and neck cancer patients and to explore their role in predicting radiation induced toxicity. Seventy nine consecutive patients suffering from head and neck cancer, diagnosed and treated in our institution, were included in the study. Toxicity was evaluated using the Radiation Therapy Oncology Group scale. Peripheral blood lymphocytes were isolated and irradiated at 0, 1, 2 and 8 Gy during 24 hours. Apoptosis was measured by flow cytometry using annexin V/propidium iodide. Lymphocytes were marked with CD45 APC-conjugated monoclonal antibody. Radiation-induced apoptosis increased in order to radiation dose and fitted to a semi logarithmic model defined by two constants: α and β. α, as the origin of the curve in the Y axis determining the percentage of spontaneous cell death, and β, as the slope of the curve determining the percentage of cell death induced at a determined radiation dose, were obtained. β value was statistically associated to normal tissue toxicity in terms of severe xerostomia, as higher levels of apoptosis were observed in patients with low toxicity (p = 0.035; Exp(B) 0.224, I.C.95% (0.060-0.904)). These data agree with our previous results and suggest that it is possible to estimate the radiosensitivity of peripheral blood lymphocytes from patients determining the radiation induced apoptosis with annexin V/propidium iodide staining. β values observed define an individual radiosensitivity profile that could predict late toxicity due to radiotherapy in locally advanced head and neck cancer patients. Anyhow, prospective studies with different cancer types and higher number of patients are needed to validate these results.

Constitutive gene expression profile segregates toxicity in locally advanced breast cancer patients treated with high-dose hyperfractionated radical radiotherapy

[EN] Breast cancer patients show a wide variation in normal tissue reactions after radiotherapy. The individual sensitivity to x-rays limits the efficiency of the therapy. Prediction of individual sensitivity to radiotherapy could help to select the radiation protocol and to improve treatment results. The aim of this study was to assess the relationship between gene expression profiles of ex vivo un-irradiated and irradiated lymphocytes and the development of toxicity due to high-dose hyperfractionated radiotherapy in patients with locally advanced breast cancer. Raw data from microarray experiments were uploaded to the Gene Expression Omnibus Database http://www.ncbi.nlm.nih.gov/geo/ (GEO accession GSE15341). We obtained a small group of 81 genes significantly regulated by radiotherapy, lumped in 50 relevant pathways. Using ANOVA and t-test statistical tools we found 20 and 26 constitutive genes (0 Gy) that segregate patients with and without acute and late toxicity, respectively. Non-supervised hierarchical clustering was used for the visualization of results. Six and 9 pathways were significantly regulated respectively. Concerning to irradiated lymphocytes (2 Gy), we founded 29 genes that separate patients with acute toxicity and without it. Those genes were gathered in 4 significant pathways. We could not identify a set of genes that segregates patients with and without late toxicity. In conclusion, we have found an association between the constitutive gene expression profile of peripheral blood lymphocytes and the development of acute and late toxicity in consecutive, unselected patients. These observations suggest the possibility of predicting normal tissue response to irradiation in high-dose non-conventional radiation therapy regimens. Prospective studies with higher number of patients are needed to validate these preliminary results.

Analyzing the dependence structure of microarray data: a copula–based approach

The main aim of this Ph.D. dissertation is the study of clustering dependent data by means of copula functions with particular emphasis on microarray data. Copula functions are a popular multivariate modeling tool in each field where the multivariate dependence is of great interest and their use in clustering has not been still investigated. The first part of this work contains the review of the literature of clustering methods, copula functions and microarray experiments. The attention focuses on the K–means (Hartigan, 1975; Hartigan and Wong, 1979), the hierarchical (Everitt, 1974) and the model–based (Fraley and Raftery, 1998, 1999, 2000, 2007) clustering techniques because their performance is compared. Then, the probabilistic interpretation of the Sklar’s theorem (Sklar’s, 1959), the estimation methods for copulas like the Inference for Margins (Joe and Xu, 1996) and the Archimedean and Elliptical copula families are presented. In the end, applications of clustering methods and copulas to the genetic and microarray experiments are highlighted. The second part contains the original contribution proposed. A simulation study is performed in order to evaluate the performance of the K–means and the hierarchical bottom–up clustering methods in identifying clusters according to the dependence structure of the data generating process. Different simulations are performed by varying different conditions (e.g., the kind of margins (distinct, overlapping and nested) and the value of the dependence parameter ) and the results are evaluated by means of different measures of performance. In light of the simulation results and of the limits of the two investigated clustering methods, a new clustering algorithm based on copula functions (‘CoClust’ in brief) is proposed. The basic idea, the iterative procedure of the CoClust and the description of the written R functions with their output are given. The CoClust algorithm is tested on simulated data (by varying the number of clusters, the copula models, the dependence parameter value and the degree of overlap of margins) and is compared with the performance of model–based clustering by using different measures of performance, like the percentage of well–identified number of clusters and the not rejection percentage of H0 on . It is shown that the CoClust algorithm allows to overcome all observed limits of the other investigated clustering techniques and is able to identify clusters according to the dependence structure of the data independently of the degree of overlap of margins and the strength of the dependence. The CoClust uses a criterion based on the maximized log–likelihood function of the copula and can virtually account for any possible dependence relationship between observations. Many peculiar characteristics are shown for the CoClust, e.g. its capability of identifying the true number of clusters and the fact that it does not require a starting classification. Finally, the CoClust algorithm is applied to the real microarray data of Hedenfalk et al. (2001) both to the gene expressions observed in three different cancer samples and to the columns (tumor samples) of the whole data matrix.