BACKGROUND

Non-small-cell lung cancer with sensitive mutations of the epidermal growth factor receptor (EGFR) is highly responsive to EGFR tyrosine kinase inhibitors such as gefitinib, but little is known about how its efficacy and safety profile compares with that of standard chemotherapy.

METHODS

We randomly assigned 230 patients with metastatic, non-small-cell lung cancer and EGFR mutations who had not previously received chemotherapy to receive gefitinib or carboplatin-paclitaxel. The primary end point was progression-free survival; secondary end points included overall survival, response rate, and toxic effects.

RESULTS

In the planned interim analysis of data for the first 200 patients, progression-free survival was significantly longer in the gefitinib group than in the standard-chemotherapy group (hazard ratio for death or disease progression with gefitinib, 0.36; P<0.001), resulting in early termination of the study. The gefitinib group had a significantly longer median progression-free survival (10.8 months, vs. 5.4 months in the chemotherapy group; hazard ratio, 0.30; 95% confidence interval, 0.22 to 0.41; P<0.001), as well as a higher response rate (73.7% vs. 30.7%, P<0.001). The median overall survival was 30.5 months in the gefitinib group and 23.6 months in the chemotherapy group (P=0.31). The most common adverse events in the gefitinib group were rash (71.1%) and elevated aminotransferase levels (55.3%), and in the chemotherapy group, neutropenia (77.0%), anemia (64.6%), appetite loss (56.6%), and sensory neuropathy (54.9%). One patient receiving gefitinib died from interstitial lung disease.

CONCLUSIONS

First-line gefitinib for patients with advanced non-small-cell lung cancer who were selected on the basis of EGFR mutations improved progression-free survival, with acceptable toxicity, as compared with standard chemotherapy. (UMIN-CTR number, C000000376.)

Carcinoma of the breast and ovary account for one-third of all cancers occurring in women and together are responsible for approximately one-quarter of cancer-related deaths in females. The HER-2/neu proto-oncogene is amplified in 25 to 30 percent of human primary breast cancers and this alteration is associated with disease behavior. In this report, several similarities were found in the biology of HER-2/neu in breast and ovarian cancer, including a similar incidence of amplification, a direct correlation between amplification and over-expression, evidence of tumors in which overexpression occurs without amplification, and the association between gene alteration and clinical outcome. A comprehensive study of the gene and its products (RNA and protein) was simultaneously performed on a large number of both tumor types. This analysis identified several potential shortcomings of the various methods used to evaluate HER-2/neu in these diseases (Southern, Northern, and Western blots, and immunohistochemistry) and provided information regarding considerations that should be addressed when studying a gene or gene product in human tissue. The data presented further support the concept that the HER-2/neu gene may be involved in the pathogenesis of some human cancers.

A novel microRNA (miRNA) quantification method has been developed using stem-loop RT followed by TaqMan PCR analysis. Stem-loop RT primers are better than conventional ones in terms of RT efficiency and specificity. TaqMan miRNA assays are specific for mature miRNAs and discriminate among related miRNAs that differ by as little as one nucleotide. Furthermore, they are not affected by genomic DNA contamination. Precise quantification is achieved routinely with as little as 25 pg of total RNA for most miRNAs. In fact, the high sensitivity, specificity and precision of this method allows for direct analysis of a single cell without nucleic acid purification. Like standard TaqMan gene expression assays, TaqMan miRNA assays exhibit a dynamic range of seven orders of magnitude. Quantification of five miRNAs in seven mouse tissues showed variation from less than 10 to more than 30,000 copies per cell. This method enables fast, accurate and sensitive miRNA expression profiling and can identify and monitor potential biomarkers specific to tissues or diseases. Stem-loop RT-PCR can be used for the quantification of other small RNA molecules such as short interfering RNAs (siRNAs). Furthermore, the concept of stem-loop RT primer design could be applied in small RNA cloning and multiplex assays for better specificity and efficiency.

OBJECTIVE

Most patients with advanced pancreas cancer experience pain and must limit their daily activities because of tumor-related symptoms. To date, no treatment has had a significant impact on the disease. In early studies with gemcitabine, patients with pancreas cancer experienced an improvement in disease-related symptoms. Based on those findings, a definitive trial was performed to assess the effectiveness of gemcitabine in patients with newly diagnosed advanced pancreas cancer.

METHODS

One hundred twenty-six patients with advanced symptomatic pancreas cancer completed a lead-in period to characterize and stabilize pain and were randomized to receive either gemcitabine 1,000 mg/m2 weekly x 7 followed by 1 week of rest, then weekly x 3 every 4 weeks thereafter (63 patients), or to fluorouracil (5-FU) 600 mg/m2 once weekly (63 patients). The primary efficacy measure was clinical benefit response, which was a composite of measurements of pain (analgesic consumption and pain intensity), Karnofsky performance status, and weight. Clinical benefit required a sustained >> or = 4 weeks) improvement in at least one parameter without worsening in any others. Other measures of efficacy included response rate, time to progressive disease, and survival.

RESULTS

Clinical benefit response was experienced by 23.8% of gemcitabine-treated patients compared with 4.8% of 5-FU-treated patients (P = .0022). The median survival durations were 5.65 and 4.41 months for gemcitabine-treated and 5-FU-treated patients, respectively (P = .0025). The survival rate at 12 months was 18% for gemcitabine patients and 2% for 5-FU patients. Treatment was well tolerated.

CONCLUSIONS

This study demonstrates that gemcitabine is more effective than 5-FU in alleviation of some disease-related symptoms in patients with advanced, symptomatic pancreas cancer. Gemcitabine also confers a modest survival advantage over treatment with 5-FU.

Immune checkpoint inhibitors, which unleash a patient's own T cells to kill tumors, are revolutionizing cancer treatment. To unravel the genomic determinants of response to this therapy, we used whole-exome sequencing of non-small cell lung cancers treated with pembrolizumab, an antibody targeting programmed cell death-1 (PD-1). In two independent cohorts, higher nonsynonymous mutation burden in tumors was associated with improved objective response, durable clinical benefit, and progression-free survival. Efficacy also correlated with the molecular smoking signature, higher neoantigen burden, and DNA repair pathway mutations; each factor was also associated with mutation burden. In one responder, neoantigen-specific CD8+ T cell responses paralleled tumor regression, suggesting that anti-PD-1 therapy enhances neoantigen-specific T cell reactivity. Our results suggest that the genomic landscape of lung cancers shapes response to anti-PD-1 therapy.

Complete knowledge of all direct and indirect interactions between proteins in a given cell would represent an important milestone towards a comprehensive description of cellular mechanisms and functions. Although this goal is still elusive, considerable progress has been made-particularly for certain model organisms and functional systems. Currently, protein interactions and associations are annotated at various levels of detail in online resources, ranging from raw data repositories to highly formalized pathway databases. For many applications, a global view of all the available interaction data is desirable, including lower-quality data and/or computational predictions. The STRING database (http://string-db.org/) aims to provide such a global perspective for as many organisms as feasible. Known and predicted associations are scored and integrated, resulting in comprehensive protein networks covering >1100 organisms. Here, we describe the update to version 9.1 of STRING, introducing several improvements: (i) we extend the automated mining of scientific texts for interaction information, to now also include full-text articles; (ii) we entirely re-designed the algorithm for transferring interactions from one model organism to the other; and (iii) we provide users with statistical information on any functional enrichment observed in their networks.

BACKGROUND

Treatment with peginterferon alfa-2a alone produces significantly higher sustained virologic responses than treatment with interferon alfa-2a alone in patients with chronic hepatitis C virus (HCV) infection. We compared the efficacy and safety of peginterferon alfa-2a plus ribavirin, interferon alfa-2b plus ribavirin, and peginterferon alfa-2a alone in the initial treatment of chronic hepatitis C.

METHODS

A total of 1121 patients were randomly assigned to treatment and received at least one dose of study medication, consisting of 180 microg of peginterferon alfa-2a once weekly plus daily ribavirin (1000 or 1200 mg, depending on body weight), weekly peginterferon alfa-2a plus daily placebo, or 3 million units of interferon alfa-2b thrice weekly plus daily ribavirin for 48 weeks.

RESULTS

A significantly higher proportion of patients who received peginterferon alfa-2a plus ribavirin had a sustained virologic response (defined as the absence of detectable HCV RNA 24 weeks after cessation of therapy) than of patients who received interferon alfa-2b plus ribavirin (56 percent vs. 44 percent, P<0.001) or peginterferon alfa-2a alone (56 percent vs. 29 percent, P<0.001). The proportions of patients with HCV genotype 1 who had sustained virologic responses were 46 percent, 36 percent, and 21 percent, respectively, for the three regimens. Among patients with HCV genotype 1 and high base-line levels of HCV RNA, the proportions of those with sustained virologic responses were 41 percent, 33 percent, and 13 percent, respectively. The overall safety profiles of the three treatment regimens were similar; the incidence of influenza-like symptoms and depression was lower in the groups receiving peginterferon alfa-2a than in the group receiving interferon alfa-2b plus ribavirin.

CONCLUSIONS

In patients with chronic hepatitis C, once-weekly peginterferon alfa-2a plus ribavirin was tolerated as well as interferon alfa-2b plus ribavirin and produced significant improvements in the rate of sustained virologic response, as compared with interferon alfa-2b plus ribavirin or peginterferon alfa-2a alone.

Methods for alignment of protein sequences typically measure similarity by using a substitution matrix with scores for all possible exchanges of one amino acid with another. The most widely used matrices are based on the Dayhoff model of evolutionary rates. Using a different approach, we have derived substitution matrices from about 2000 blocks of aligned sequence segments characterizing more than 500 groups of related proteins. This led to marked improvements in alignments and in searches using queries from each of the groups.

Drug therapy for hypercholesterolaemia has remained controversial mainly because of insufficient clinical trial evidence for improved survival. The present trial was designed to evaluate the effect of cholesterol lowering with simvastatin on mortality and morbidity in patients with coronary heart disease (CHD). 4444 patients with angina pectoris or previous myocardial infarction and serum cholesterol 5.5-8.0 mmol/L on a lipid-lowering diet were randomised to double-blind treatment with simvastatin or placebo. Over the 5.4 years median follow-up period, simvastatin produced mean changes in total cholesterol, low-density-lipoprotein cholesterol, and high-density-lipoprotein cholesterol of -25%, -35%, and +8%, respectively, with few adverse effects. 256 patients (12%) in the placebo group died, compared with 182 (8%) in the simvastatin group. The relative risk of death in the simvastatin group was 0.70 (95% CI 0.58-0.85, p = 0.0003). The 6-year probabilities of survival in the placebo and simvastatin groups were 87.6% and 91.3%, respectively. There were 189 coronary deaths in the placebo group and 111 in the simvastatin group (relative risk 0.58, 95% CI 0.46-0.73), while noncardiovascular causes accounted for 49 and 46 deaths, respectively. 622 patients (28%) in the placebo group and 431 (19%) in the simvastatin group had one or more major coronary events. The relative risk was 0.66 (95% CI 0.59-0.75, p < 0.00001), and the respective probabilities of escaping such events were 70.5% and 79.6%. This risk was also significantly reduced in subgroups consisting of women and patients of both sexes aged 60 or more. Other benefits of treatment included a 37% reduction (p < 0.00001) in the risk of undergoing myocardial revascularisation procedures. This study shows that long-term treatment with simvastatin is safe and improves survival in CHD patients.

pRK212.2, a derivative of the broad host range plasmid RK2, contains two EcoRI cleavage fragments, A and B, neither of which can replicate by itself in Escherichia coli. Fragment A (41.7 kilobases), but not fragment B (14.4 kilobases), can be cloned by insertion into the unrelated plasmids mini-F and ColE1. Fragment B contains the origin of replication and the ampicillin-resistance determinant of RK2. Transformation of E. coli cells containing the mini-F-fragment A hybrid plasmid with fragment B DNA results in the recircularization and replication of fragment B as a nonmobilizable plasmid (pRK2067) with the copy number and incompatibility properties of RK2. Fragment B cannot be cloned in the absence of fragment A because the latter fragment suppresses a function, specified by fragment B, that results in loss of host cell viability. A small segment (2.4 kilobases) of fragment B that contains the RK2 origin of replication but no longer affects host cell growth in the absence of fragment A had been cloned previously by insertion into a ColE1 plasmid. This hybrid plasmid, designated pRK256, will replicate in E. coli polA mutants only when a fragment A-bearing helper plasmid is present. These results demonstrate that the potentially lethal function specified by fragment B of RK2 is not necessary for replication and that at least one trans-acting function is directly involved in RK2 replication.

An endothelium-derived 21-residue vasoconstrictor peptide, endothelin, has been isolated, and shown to be one of the most potent vasoconstrictors known. Cloning and sequencing of preproendothelin complementary DNA shows that mature endothelin is generated through an unusual proteolytic processing, and regional homologies to a group of neurotoxins suggest that endothelin is an endogenous modulator of voltage-dependent ion channels. Expression of the endothelin gene is regulated by several vasoactive agents, indicating the existence of a novel cardiovascular control system.

The human genome contains many thousands of long noncoding RNAs (lncRNAs). While several studies have demonstrated compelling biological and disease roles for individual examples, analytical and experimental approaches to investigate these genes have been hampered by the lack of comprehensive lncRNA annotation. Here, we present and analyze the most complete human lncRNA annotation to date, produced by the GENCODE consortium within the framework of the ENCODE project and comprising 9277 manually annotated genes producing 14,880 transcripts. Our analyses indicate that lncRNAs are generated through pathways similar to that of protein-coding genes, with similar histone-modification profiles, splicing signals, and exon/intron lengths. In contrast to protein-coding genes, however, lncRNAs display a striking bias toward two-exon transcripts, they are predominantly localized in the chromatin and nucleus, and a fraction appear to be preferentially processed into small RNAs. They are under stronger selective pressure than neutrally evolving sequences-particularly in their promoter regions, which display levels of selection comparable to protein-coding genes. Importantly, about one-third seem to have arisen within the primate lineage. Comprehensive analysis of their expression in multiple human organs and brain regions shows that lncRNAs are generally lower expressed than protein-coding genes, and display more tissue-specific expression patterns, with a large fraction of tissue-specific lncRNAs expressed in the brain. Expression correlation analysis indicates that lncRNAs show particularly striking positive correlation with the expression of antisense coding genes. This GENCODE annotation represents a valuable resource for future studies of lncRNAs.

Genotype imputation methods are now being widely used in the analysis of genome-wide association studies. Most imputation analyses to date have used the HapMap as a reference dataset, but new reference panels (such as controls genotyped on multiple SNP chips and densely typed samples from the 1,000 Genomes Project) will soon allow a broader range of SNPs to be imputed with higher accuracy, thereby increasing power. We describe a genotype imputation method (IMPUTE version 2) that is designed to address the challenges presented by these new datasets. The main innovation of our approach is a flexible modelling framework that increases accuracy and combines information across multiple reference panels while remaining computationally feasible. We find that IMPUTE v2 attains higher accuracy than other methods when the HapMap provides the sole reference panel, but that the size of the panel constrains the improvements that can be made. We also find that imputation accuracy can be greatly enhanced by expanding the reference panel to contain thousands of chromosomes and that IMPUTE v2 outperforms other methods in this setting at both rare and common SNPs, with overall error rates that are 15%-20% lower than those of the closest competing method. One particularly challenging aspect of next-generation association studies is to integrate information across multiple reference panels genotyped on different sets of SNPs; we show that our approach to this problem has practical advantages over other suggested solutions.

We have previously shown that a brief episode of ischemia slows the rate of ATP depletion during subsequent ischemic episodes. Additionally, intermittent reperfusion may be beneficial to the myocardium by washing out catabolites that have accumulated during ischemia. Thus, we proposed that multiple brief ischemic episodes might actually protect the heart from a subsequent sustained ischemic insult. To test this hypothesis, two sets of experiments were performed. In the first set, one group of dogs (n = 7) was preconditioned with four 5 min circumflex occlusions, each separated by 5 min of reperfusion, followed by a sustained 40 min occlusion. The control group (n = 5) received a single 40 min occlusion. In the second study, an identical preconditioning protocol was followed, and animals (n = 9) then received a sustained 3 hr occlusion. Control animals (n = 7) received a single 3 hr occlusion. Animals were allowed 4 days of reperfusion thereafter. Histologic infarct size then was measured and was related to the major baseline predictors of infarct size, including the anatomic area at risk and collateral blood flow. In the 40 min study, preconditioning with ischemia paradoxically limited infarct size to 25% of that seen in the control group (p less than .001). Collateral blood flows were not significantly different in the two groups. In the 3 hr study, there was no difference between infarct size in the preconditioned and control groups. The protective effect of preconditioning in the 40 min study may have been due to reduced ATP depletion and/or to reduced catabolite accumulation during the sustained occlusion. These results suggest that the multiple anginal episodes that often precede myocardial infarction in man may delay cell death after coronary occlusion, and thereby allow for greater salvage of myocardium through reperfusion therapy.

In 2008 we published the first set of guidelines for standardizing research in autophagy. Since then, research on this topic has continued to accelerate, and many new scientists have entered the field. Our knowledge base and relevant new technologies have also been expanding. Accordingly, it is important to update these guidelines for monitoring autophagy in different organisms. Various reviews have described the range of assays that have been used for this purpose. Nevertheless, there continues to be confusion regarding acceptable methods to measure autophagy, especially in multicellular eukaryotes. A key point that needs to be emphasized is that there is a difference between measurements that monitor the numbers or volume of autophagic elements (e.g., autophagosomes or autolysosomes) at any stage of the autophagic process vs. those that measure flux through the autophagy pathway (i.e., the complete process); thus, a block in macroautophagy that results in autophagosome accumulation needs to be differentiated from stimuli that result in increased autophagic activity, defined as increased autophagy induction coupled with increased delivery to, and degradation within, lysosomes (in most higher eukaryotes and some protists such as Dictyostelium) or the vacuole (in plants and fungi). In other words, it is especially important that investigators new to the field understand that the appearance of more autophagosomes does not necessarily equate with more autophagy. In fact, in many cases, autophagosomes accumulate because of a block in trafficking to lysosomes without a concomitant change in autophagosome biogenesis, whereas an increase in autolysosomes may reflect a reduction in degradative activity. Here, we present a set of guidelines for the selection and interpretation of methods for use by investigators who aim to examine macroautophagy and related processes, as well as for reviewers who need to provide realistic and reasonable critiques of papers that are focused on these processes. These guidelines are not meant to be a formulaic set of rules, because the appropriate assays depend in part on the question being asked and the system being used. In addition, we emphasize that no individual assay is guaranteed to be the most appropriate one in every situation, and we strongly recommend the use of multiple assays to monitor autophagy. In these guidelines, we consider these various methods of assessing autophagy and what information can, or cannot, be obtained from them. Finally, by discussing the merits and limits of particular autophagy assays, we hope to encourage technical innovation in the field.

Finding objective and effective thresholds for voxelwise statistics derived from neuroimaging data has been a long-standing problem. With at least one test performed for every voxel in an image, some correction of the thresholds is needed to control the error rates, but standard procedures for multiple hypothesis testing (e.g., Bonferroni) tend to not be sensitive enough to be useful in this context. This paper introduces to the neuroscience literature statistical procedures for controlling the false discovery rate (FDR). Recent theoretical work in statistics suggests that FDR-controlling procedures will be effective for the analysis of neuroimaging data. These procedures operate simultaneously on all voxelwise test statistics to determine which tests should be considered statistically significant. The innovation of the procedures is that they control the expected proportion of the rejected hypotheses that are falsely rejected. We demonstrate this approach using both simulations and functional magnetic resonance imaging data from two simple experiments.

CD4+CD25+ regulatory T cells (Treg) are instrumental in the maintenance of immunological tolerance. One critical question is whether Treg can only be generated in the thymus or can differentiate from peripheral CD4+CD25- naive T cells. In this paper, we present novel evidence that conversion of naive peripheral CD4+CD25- T cells into anergic/suppressor cells that are CD25+, CD45RB-/low and intracellular CTLA-4+ can be achieved through costimulation with T cell receptors (TCRs) and transforming growth factor beta (TGF-beta). Although transcription factor Foxp3 has been shown recently to be associated with the development of Treg, the physiological inducers for Foxp3 gene expression remain a mystery. TGF-beta induced Foxp3 gene expression in TCR-challenged CD4+CD25- naive T cells, which mediated their transition toward a regulatory T cell phenotype with potent immunosuppressive potential. These converted anergic/suppressor cells are not only unresponsive to TCR stimulation and produce neither T helper cell 1 nor T helper cell 2 cytokines but they also express TGF-beta and inhibit normal T cell proliferation in vitro. More importantly, in an ovalbumin peptide TCR transgenic adoptive transfer model, TGF-beta-converted transgenic CD4+CD25+ suppressor cells proliferated in response to immunization and inhibited antigen-specific naive CD4+ T cell expansion in vivo. Finally, in a murine asthma model, coadministration of these TGF-beta-induced suppressor T cells prevented house dust mite-induced allergic pathogenesis in lungs.

BCL-2 family proteins, which have either pro- or anti-apoptotic activities, have been studied intensively for the past decade owing to their importance in the regulation of apoptosis, tumorigenesis and cellular responses to anti-cancer therapy. They control the point of no return for clonogenic cell survival and thereby affect tumorigenesis and host-pathogen interactions and regulate animal development. Recent structural, phylogenetic and biological analyses, however, suggest the need for some reconsideration of the accepted organizational principles of the family and how the family members interact with one another during programmed cell death. Although these insights into interactions among BCL-2 family proteins reveal how these proteins are regulated, a unifying hypothesis for the mechanisms they use to activate caspases remains elusive.

BACKGROUND

Prediction of 3-dimensional protein structures from amino acid sequences represents one of the most important problems in computational structural biology. The community-wide Critical Assessment of Structure Prediction (CASP) experiments have been designed to obtain an objective assessment of the state-of-the-art of the field, where I-TASSER was ranked as the best method in the server section of the recent 7th CASP experiment. Our laboratory has since then received numerous requests about the public availability of the I-TASSER algorithm and the usage of the I-TASSER predictions.

RESULTS

An on-line version of I-TASSER is developed at the KU Center for Bioinformatics which has generated protein structure predictions for thousands of modeling requests from more than 35 countries. A scoring function (C-score) based on the relative clustering structural density and the consensus significance score of multiple threading templates is introduced to estimate the accuracy of the I-TASSER predictions. A large-scale benchmark test demonstrates a strong correlation between the C-score and the TM-score (a structural similarity measurement with values in [0, 1]) of the first models with a correlation coefficient of 0.91. Using a C-score cutoff>> -1.5 for the models of correct topology, both false positive and false negative rates are below 0.1. Combining C-score and protein length, the accuracy of the I-TASSER models can be predicted with an average error of 0.08 for TM-score and 2 A for RMSD.

CONCLUSIONS

The I-TASSER server has been developed to generate automated full-length 3D protein structural predictions where the benchmarked scoring system helps users to obtain quantitative assessments of the I-TASSER models. The output of the I-TASSER server for each query includes up to five full-length models, the confidence score, the estimated TM-score and RMSD, and the standard deviation of the estimations. The I-TASSER server is freely available to the academic community at http://zhang.bioinformatics.ku.edu/I-TASSER.

An algorithm is presented for the multiple alignment of sequences, either proteins or nucleic acids, that is both accurate and easy to use on microcomputers. The approach is based on the conventional dynamic-programming method of pairwise alignment. Initially, a hierarchical clustering of the sequences is performed using the matrix of the pairwise alignment scores. The closest sequences are aligned creating groups of aligned sequences. Then close groups are aligned until all sequences are aligned in one group. The pairwise alignments included in the multiple alignment form a new matrix that is used to produce a hierarchical clustering. If it is different from the first one, iteration of the process can be performed. The method is illustrated by an example: a global alignment of 39 sequences of cytochrome c.

Understanding how the immune system affects cancer development and progression has been one of the most challenging questions in immunology. Research over the past two decades has helped explain why the answer to this question has evaded us for so long. We now appreciate that the immune system plays a dual role in cancer: It can not only suppress tumor growth by destroying cancer cells or inhibiting their outgrowth but also promote tumor progression either by selecting for tumor cells that are more fit to survive in an immunocompetent host or by establishing conditions within the tumor microenvironment that facilitate tumor outgrowth. Here, we discuss a unifying conceptual framework called "cancer immunoediting," which integrates the immune system's dual host-protective and tumor-promoting roles.

The anterior insular cortex (AIC) is implicated in a wide range of conditions and behaviours, from bowel distension and orgasm, to cigarette craving and maternal love, to decision making and sudden insight. Its function in the re-representation of interoception offers one possible basis for its involvement in all subjective feelings. New findings suggest a fundamental role for the AIC (and the von Economo neurons it contains) in awareness, and thus it needs to be considered as a potential neural correlate of consciousness.

Colon carcinoma is the second most common cause of death from cancer. The isolation and characterization of tumorigenic colon cancer cells may help to devise novel diagnostic and therapeutic procedures. Although there is increasing evidence that a rare population of undifferentiated cells is responsible for tumour formation and maintenance, this has not been explored for colorectal cancer. Here, we show that tumorigenic cells in colon cancer are included in the high-density CD133+ population, which accounts for about 2.5% of the tumour cells. Subcutaneous injection of colon cancer CD133+ cells readily reproduced the original tumour in immunodeficient mice, whereas CD133- cells did not form tumours. Such tumours were serially transplanted for several generations, in each of which we observed progressively faster tumour growth without significant phenotypic alterations. Unlike CD133- cells, CD133+ colon cancer cells grew exponentially for more than one year in vitro as undifferentiated tumour spheres in serum-free medium, maintaining the ability to engraft and reproduce the same morphological and antigenic pattern of the original tumour. We conclude that colorectal cancer is created and propagated by a small number of undifferentiated tumorigenic CD133+ cells, which should therefore be the target of future therapies.