Ancestral polyploidy in seed plants and angiosperms.

Whole-genome duplication (WGD), or polyploidy, followed by gene loss and diploidization has long been recognized as an important evolutionary force in animals, fungi and other organisms, especially plants. The success of angiosperms has been attributed, in part, to innovations associated with gene or whole-genome duplications, but evidence for proposed ancient genome duplications pre-dating the divergence of monocots and eudicots remains equivocal in analyses of conserved gene order. Here we use comprehensive phylogenomic analyses of sequenced plant genomes and more than 12.6 million new expressed-sequence-tag sequences from phylogenetically pivotal lineages to elucidate two groups of ancient gene duplications-one in the common ancestor of extant seed plants and the other in the common ancestor of extant angiosperms. Gene duplication events were intensely concentrated around 319 and 192 million years ago, implicating two WGDs in ancestral lineages shortly before the diversification of extant seed plants and extant angiosperms, respectively. Significantly, these ancestral WGDs resulted in the diversification of regulatory genes important to seed and flower development, suggesting that they were involved in major innovations that ultimately contributed to the rise and eventual dominance of seed plants and angiosperms.

Diverse somatic mutation patterns and pathway alterations in human cancers.

The systematic characterization of somatic mutations in cancer genomes is essential for understanding the disease and for developing targeted therapeutics. Here we report the identification of 2,576 somatic mutations across approximately 1,800 megabases of DNA representing 1,507 coding genes from 441 tumours comprising breast, lung, ovarian and prostate cancer types and subtypes. We found that mutation rates and the sets of mutated genes varied substantially across tumour types and subtypes. Statistical analysis identified 77 significantly mutated genes including protein kinases, G-protein-coupled receptors such as GRM8, BAI3, AGTRL1 (also called APLNR) and LPHN3, and other druggable targets. Integrated analysis of somatic mutations and copy number alterations identified another 35 significantly altered genes including GNAS, indicating an expanded role for galpha subunits in multiple cancer types. Furthermore, our experimental analyses demonstrate the functional roles of mutant GNAO1 (a Galpha subunit) and mutant MAP2K4 (a member of the JNK signalling pathway) in oncogenesis. Our study provides an overview of the mutational spectra across major human cancers and identifies several potential therapeutic targets.

Complete Khoisan and Bantu genomes from southern Africa.

The genetic structure of the indigenous hunter-gatherer peoples of southern Africa, the oldest known lineage of modern human, is important for understanding human diversity. Studies based on mitochondrial and small sets of nuclear markers have shown that these hunter-gatherers, known as Khoisan, San, or Bushmen, are genetically divergent from other humans. However, until now, fully sequenced human genomes have been limited to recently diverged populations. Here we present the complete genome sequences of an indigenous hunter-gatherer from the Kalahari Desert and a Bantu from southern Africa, as well as protein-coding regions from an additional three hunter-gatherers from disparate regions of the Kalahari. We characterize the extent of whole-genome and exome diversity among the five men, reporting 1.3 million novel DNA differences genome-wide, including 13,146 novel amino acid variants. In terms of nucleotide substitutions, the Bushmen seem to be, on average, more different from each other than, for example, a European and an Asian. Observed genomic differences between the hunter-gatherers and others may help to pinpoint genetic adaptations to an agricultural lifestyle. Adding the described variants to current databases will facilitate inclusion of southern Africans in medical research efforts, particularly when family and medical histories can be correlated with genome-wide data.

Comparative and population mitogenomic analyses of Madagascar's extinct, giant 'subfossil' lemurs.

Humans first arrived on Madagascar only a few thousand years ago. Subsequent habitat destruction and hunting activities have had significant impacts on the island's biodiversity, including the extinction of megafauna. For example, we know of 17 recently extinct 'subfossil' lemur species, all of which were substantially larger (body mass ∼11-160 kg) than any living population of the ∼100 extant lemur species (largest body mass ∼6.8 kg). We used ancient DNA and genomic methods to study subfossil lemur extinction biology and update our understanding of extant lemur conservation risk factors by i) reconstructing a comprehensive phylogeny of extinct and extant lemurs, and ii) testing whether low genetic diversity is associated with body size and extinction risk. We recovered complete or near-complete mitochondrial genomes from five subfossil lemur taxa, and generated sequence data from population samples of two extinct and eight extant lemur species. Phylogenetic comparisons resolved prior taxonomic uncertainties and confirmed that the extinct subfossil species did not comprise a single clade. Genetic diversity estimates for the two sampled extinct species were relatively low, suggesting small historical population sizes. Low genetic diversity and small population sizes are both risk factors that would have rendered giant lemurs especially susceptible to extinction. Surprisingly, among the extant lemurs, we did not observe a relationship between body size and genetic diversity. The decoupling of these variables suggests that risk factors other than body size may have as much or more meaning for establishing future lemur conservation priorities.

Polar and brown bear genomes reveal ancient admixture and demographic footprints of past climate change.

Polar bears (PBs) are superbly adapted to the extreme Arctic environment and have become emblematic of the threat to biodiversity from global climate change. Their divergence from the lower-latitude brown bear provides a textbook example of rapid evolution of distinct phenotypes. However, limited mitochondrial and nuclear DNA evidence conflicts in the timing of PB origin as well as placement of the species within versus sister to the brown bear lineage. We gathered extensive genomic sequence data from contemporary polar, brown, and American black bear samples, in addition to a 130,000- to 110,000-y old PB, to examine this problem from a genome-wide perspective. Nuclear DNA markers reflect a species tree consistent with expectation, showing polar and brown bears to be sister species. However, for the enigmatic brown bears native to Alaska's Alexander Archipelago, we estimate that not only their mitochondrial genome, but also 5-10% of their nuclear genome, is most closely related to PBs, indicating ancient admixture between the two species. Explicit admixture analyses are consistent with ancient splits among PBs, brown bears and black bears that were later followed by occasional admixture. We also provide paleodemographic estimates that suggest bear evolution has tracked key climate events, and that PB in particular experienced a prolonged and dramatic decline in its effective population size during the last ca. 500,000 years. We demonstrate that brown bears and PBs have had sufficiently independent evolutionary histories over the last 4-5 million years to leave imprints in the PB nuclear genome that likely are associated with ecological adaptation to the Arctic environment.

Nucleosome organization in the Drosophila genome.

Comparative genomics of nucleosome positions provides a powerful means for understanding how the organization of chromatin and the transcription machinery co-evolve. Here we produce a high-resolution reference map of H2A.Z and bulk nucleosome locations across the genome of the fly Drosophila melanogaster and compare it to that from the yeast Saccharomyces cerevisiae. Like Saccharomyces, Drosophila nucleosomes are organized around active transcription start sites in a canonical -1, nucleosome-free region, +1 arrangement. However, Drosophila does not incorporate H2A.Z into the -1 nucleosome and does not bury its transcriptional start site in the +1 nucleosome. At thousands of genes, RNA polymerase II engages the +1 nucleosome and pauses. How the transcription initiation machinery contends with the +1 nucleosome seems to be fundamentally different across major eukaryotic lines.

Sequencing the nuclear genome of the extinct woolly mammoth.

In 1994, two independent groups extracted DNA from several Pleistocene epoch mammoths and noted differences among individual specimens. Subsequently, DNA sequences have been published for a number of extinct species. However, such ancient DNA is often fragmented and damaged, and studies to date have typically focused on short mitochondrial sequences, never yielding more than a fraction of a per cent of any nuclear genome. Here we describe 4.17 billion bases (Gb) of sequence from several mammoth specimens, 3.3 billion (80%) of which are from the woolly mammoth (Mammuthus primigenius) genome and thus comprise an extensive set of genome-wide sequence from an extinct species. Our data support earlier reports that elephantid genomes exceed 4 Gb. The estimated divergence rate between mammoth and African elephant is half of that between human and chimpanzee. The observed number of nucleotide differences between two particular mammoths was approximately one-eighth of that between one of them and the African elephant, corresponding to a separation between the mammoths of 1.5-2.0 Myr. The estimated probability that orthologous elephant and mammoth amino acids differ is 0.002, corresponding to about one residue per protein. Differences were discovered between mammoth and African elephant in amino-acid positions that are otherwise invariant over several billion years of combined mammalian evolution. This study shows that nuclear genome sequencing of extinct species can reveal population differences not evident from the fossil record, and perhaps even discover genetic factors that affect extinction.

Genetic diversity and population structure of the endangered marsupial Sarcophilus harrisii (Tasmanian devil).

The Tasmanian devil (Sarcophilus harrisii) is threatened with extinction because of a contagious cancer known as Devil Facial Tumor Disease. The inability to mount an immune response and to reject these tumors might be caused by a lack of genetic diversity within a dwindling population. Here we report a whole-genome analysis of two animals originating from extreme northwest and southeast Tasmania, the maximal geographic spread, together with the genome from a tumor taken from one of them. A 3.3-Gb de novo assembly of the sequence data from two complementary next-generation sequencing platforms was used to identify 1 million polymorphic genomic positions, roughly one-quarter of the number observed between two genetically distant human genomes. Analysis of 14 complete mitochondrial genomes from current and museum specimens, as well as mitochondrial and nuclear SNP markers in 175 animals, suggests that the observed low genetic diversity in today's population preceded the Devil Facial Tumor Disease disease outbreak by at least 100 y. Using a genetically characterized breeding stock based on the genome sequence will enable preservation of the extant genetic diversity in future Tasmanian devil populations.

Translational and rotational settings of H2A.Z nucleosomes across the Saccharomyces cerevisiae genome.

The nucleosome is the fundamental building block of eukaryotic chromosomes. Access to genetic information encoded in chromosomes is dependent on the position of nucleosomes along the DNA. Alternative locations just a few nucleotides apart can have profound effects on gene expression. Yet the nucleosomal context in which chromosomal and gene regulatory elements reside remains ill-defined on a genomic scale. Here we sequence the DNA of 322,000 individual Saccharomyces cerevisiae nucleosomes, containing the histone variant H2A.Z, to provide a comprehensive map of H2A.Z nucleosomes in functionally important regions. With a median 4-base-pair resolution, we identify new and established signatures of nucleosome positioning. A single predominant rotational setting and multiple translational settings are evident. Chromosomal elements, ranging from telomeres to centromeres and transcriptional units, are found to possess characteristic nucleosomal architecture that may be important for their function. Promoter regulatory elements, including transcription factor binding sites and transcriptional start sites, show topological relationships with nucleosomes, such that transcription factor binding sites tend to be rotationally exposed on the nucleosome surface near its border. Transcriptional start sites tended to reside about one helical turn inside the nucleosome border. These findings reveal an intimate relationship between chromatin architecture and the underlying DNA sequence it regulates.

Complete mitochondrial genome of a Pleistocene jawbone unveils the origin of polar bear.

The polar bear has become the flagship species in the climate-change discussion. However, little is known about how past climate impacted its evolution and persistence, given an extremely poor fossil record. Although it is undisputed from analyses of mitochondrial (mt) DNA that polar bears constitute a lineage within the genetic diversity of brown bears, timing estimates of their divergence have differed considerably. Using next-generation sequencing technology, we have generated a complete, high-quality mt genome from a stratigraphically validated 130,000- to 110,000-year-old polar bear jawbone. In addition, six mt genomes were generated of extant polar bears from Alaska and brown bears from the Admiralty and Baranof islands of the Alexander Archipelago of southeastern Alaska and Kodiak Island. We show that the phylogenetic position of the ancient polar bear lies almost directly at the branching point between polar bears and brown bears, elucidating a unique morphologically and molecularly documented fossil link between living mammal species. Molecular dating and stable isotope analyses also show that by very early in their evolutionary history, polar bears were already inhabitants of the Artic sea ice and had adapted very rapidly to their current and unique ecology at the top of the Arctic marine food chain. As such, polar bears provide an excellent example of evolutionary opportunism within a widespread mammalian lineage.

Complete genome of Candidatus Chloracidobacterium thermophilum, a chlorophyll-based photoheterotroph belonging to the phylum Acidobacteria.

Candidatus Chloracidobacterium thermophilum, which naturally inhabits microbial mats of alkaline siliceous hot springs in Yellowstone National Park, is the only known chlorophototroph in the phylum Acidobacteria. The Ca. C. thermophilum genome was composed of two chromosomes (2,683,362 bp and 1,012,010 bp), and both encoded essential genes. The genome included genes to produce chlorosomes, the Fenna-Matthews-Olson protein, bacteriochlorophylls a and c as principal pigments, and type-1, homodimeric reaction centres. Ca. C. thermophilum is an aerobic photoheterotroph that lacks the ability to synthesize several essential nutrients. Key genes of all known carbon fixation pathways were absent, as were genes for assimilatory nitrate and sulfate reduction and vitamin B(12) synthesis. Genes for the synthesis of branched-chain amino acids (valine, isoleucine and leucine) were also absent, but genes for catabolism of these compounds were present. This observation suggested that Ca. C. thermophilum may synthesize branched-chain amino acids from an intermediate(s) of the catabolic pathway by reversing these reactions. The genome encoded an aerobic respiratory electron transport chain that included NADH dehydrogenase, alternative complex III and cytochrome oxidase. The chromosomes of the laboratory isolate were compared with assembled, metagenomic scaffolds from the major Ca. C. thermophilum population in hot-spring mats. The larger chromosomes of the two populations were highly syntenous but significantly divergent (~13%) in sequence. In contrast, the smaller chromosomes have undergone numerous rearrangements, contained many transposases, and might be less constrained by purifying selection than the large chromosomes. Some transposases were homologous to those of mat community members from other phyla.

Intraspecific phylogenetic analysis of Siberian woolly mammoths using complete mitochondrial genomes.

We report five new complete mitochondrial DNA (mtDNA) genomes of Siberian woolly mammoth (Mammuthus primigenius), sequenced with up to 73-fold coverage from DNA extracted from hair shaft material. Three of the sequences present the first complete mtDNA genomes of mammoth clade II. Analysis of these and 13 recently published mtDNA genomes demonstrates the existence of two apparently sympatric mtDNA clades that exhibit high interclade divergence. The analytical power afforded by the analysis of the complete mtDNA genomes reveals a surprisingly ancient coalescence age of the two clades, approximately 1-2 million years, depending on the calibration technique. Furthermore, statistical analysis of the temporal distribution of the (14)C ages of these and previously identified members of the two mammoth clades suggests that clade II went extinct before clade I. Modeling of protein structures failed to indicate any important functional difference between genomes belonging to the two clades, suggesting that the loss of clade II more likely is due to genetic drift than a selective sweep.

Analysis of complete mitochondrial genomes from extinct and extant rhinoceroses reveals lack of phylogenetic resolution.

BACKGROUND: The scientific literature contains many examples where DNA sequence analyses have been used to provide definitive answers to phylogenetic problems that traditional (non-DNA based) approaches alone have failed to resolve. One notable example concerns the rhinoceroses, a group for which several contradictory phylogenies were proposed on the basis of morphology, then apparently resolved using mitochondrial DNA fragments. RESULTS: In this study we report the first complete mitochondrial genome sequences of the extinct ice-age woolly rhinoceros (Coelodonta antiquitatis), and the threatened Javan (Rhinoceros sondaicus), Sumatran (Dicerorhinus sumatrensis), and black (Diceros bicornis) rhinoceroses. In combination with the previously published mitochondrial genomes of the white (Ceratotherium simum) and Indian (Rhinoceros unicornis) rhinoceroses, this data set putatively enables reconstruction of the rhinoceros phylogeny. While the six species cluster into three strongly supported sister-pairings: (i) The black/white, (ii) the woolly/Sumatran, and (iii) the Javan/Indian, resolution of the higher-level relationships has no statistical support. The phylogenetic signal from individual genes is highly diffuse, with mixed topological support from different genes. Furthermore, the choice of outgroup (horse vs tapir) has considerable effect on reconstruction of the phylogeny. The lack of resolution is suggestive of a hard polytomy at the base of crown-group Rhinocerotidae, and this is supported by an investigation of the relative branch lengths. CONCLUSION: Satisfactory resolution of the rhinoceros phylogeny may not be achievable without additional analyses of substantial amounts of nuclear DNA. This study provides a compelling demonstration that, in spite of substantial sequence length, there are significant limitations with single-locus phylogenetics. We expect further examples of this to appear as next-generation, large-scale sequencing of complete mitochondrial genomes becomes commonplace in evolutionary studies. "The human factor in classification is nowhere more evident than in dealing with this superfamily (Rhinocerotoidea)." G. G. Simpson (1945).

Characterization of meiotic crossovers and gene conversion by whole-genome sequencing in Saccharomyces cerevisiae.

BACKGROUND: Meiotic recombination alters frequency and distribution of genetic variation, impacting genetics and evolution. In the budding yeast, DNA double strand breaks (DSBs) and D loops form either crossovers (COs) or non-crossovers (NCOs), which occur at many sites in the genome. Differences at the nucleotide level associated with COs and NCOs enable us to detect these recombination events and their distributions. RESULTS: We used high throughput sequencing to uncover over 46 thousand single nucleotide polymorphisms (SNPs) between two budding yeast strains and investigated meiotic recombinational events. We provided a detailed analysis of CO and NCO events, including number, size range, and distribution on chromosomes. We have detected 91 COs, very close to the average number from previous genetic studies, as well as 21 NCO events and mapped the positions of these events with high resolution. We have obtained DNA sequence-level evidence for a wide range of sizes of chromosomal regions involved in CO and NCO events. We show that a large fraction of the COs are accompanied by gene conversion (GC), indicating that meiotic recombination changes allelic frequencies, in addition to redistributing existing genetic variations. CONCLUSION: This work is the first reported study of meiotic recombination using high throughput sequencing technologies. Our results show that high-throughput sequencing is a sensitive method to uncover at single-base resolution details of CO and NCO events, including some complex patterns, providing new clues about the mechanism of this fundamental process.

Comparison of next generation sequencing technologies for transcriptome characterization.

BACKGROUND: We have developed a simulation approach to help determine the optimal mixture of sequencing methods for most complete and cost effective transcriptome sequencing. We compared simulation results for traditional capillary sequencing with "Next Generation" (NG) ultra high-throughput technologies. The simulation model was parameterized using mappings of 130,000 cDNA sequence reads to the Arabidopsis genome (NCBI Accession SRA008180.19). We also generated 454-GS20 sequences and de novo assemblies for the basal eudicot California poppy (Eschscholzia californica) and the magnoliid avocado (Persea americana) using a variety of methods for cDNA synthesis. RESULTS: The Arabidopsis reads tagged more than 15,000 genes, including new splice variants and extended UTR regions. Of the total 134,791 reads (13.8 MB), 119,518 (88.7%) mapped exactly to known exons, while 1,117 (0.8%) mapped to introns, 11,524 (8.6%) spanned annotated intron/exon boundaries, and 3,066 (2.3%) extended beyond the end of annotated UTRs. Sequence-based inference of relative gene expression levels correlated significantly with microarray data. As expected, NG sequencing of normalized libraries tagged more genes than non-normalized libraries, although non-normalized libraries yielded more full-length cDNA sequences. The Arabidopsis data were used to simulate additional rounds of NG and traditional EST sequencing, and various combinations of each. Our simulations suggest a combination of FLX and Solexa sequencing for optimal transcriptome coverage at modest cost. We have also developed ESTcalc http://fgp.huck.psu.edu/NG_Sims/ngsim.pl, an online webtool, which allows users to explore the results of this study by specifying individualized costs and sequencing characteristics. CONCLUSION: NG sequencing technologies are a highly flexible set of platforms that can be scaled to suit different project goals. In terms of sequence coverage alone, the NG sequencing is a dramatic advance over capillary-based sequencing, but NG sequencing also presents significant challenges in assembly and sequence accuracy due to short read lengths, method-specific sequencing errors, and the absence of physical clones. These problems may be overcome by hybrid sequencing strategies using a mixture of sequencing methodologies, by new assemblers, and by sequencing more deeply. Sequencing and microarray outcomes from multiple experiments suggest that our simulator will be useful for guiding NG transcriptome sequencing projects in a wide range of organisms.

Complete Genome Sequence of the Photoautotrophic and Bacteriochlorophyll e-Synthesizing Green Sulfur Bacterium Chlorobaculum limnaeum DSM 1677T.

Chlorobaculum limnaeum DSM 1677T is a mesophilic, brown-colored, chlorophototrophic green sulfur bacterium that produces bacteriochlorophyll e and the carotenoid isorenieratene as major pigments. This bacterium serves as a model organism in molecular research on photosynthesis, sulfur metabolism, and bacteriochlorophyll biosynthesis. We report here the complete genome sequence.

Draft Genome Sequence of Anoxybacillusayderensis Strain MT-Cab (Firmicutes).

The draft genome of the Gram-positive spore-forming Anoxybacillus ayderensis strain MT-Cab (Firmicutes), isolated from an enrichment culture of Chloracidobacterium thermophilum, was sequenced and comprises 2,577,015 bp in 92 contigs. The draft genome is predicted to consist of 2,699 protein-coding genes, 73 tRNA-coding genes, and an estimated 8 rRNA operons.

APC I1307K and the risk of prostate cancer.

The kin-cohort design has been proposed as an alternative to traditional case-control and cohort measures to evaluate inherited susceptibility to cancer in population-based studies. Here, we used this design to evaluate inherited susceptibility to prostate cancer associated with APC I1307K using data from the Molecular Epidemiology of Colorectal Cancer study. Two techniques were used to compare the incidence of prostate cancer in APC I1307K carriers. First, we compared the incidence of prostate cancer in relatives of mutation carriers and noncarriers using standard techniques for survival analysis. Second, we used the marginal maximum likelihood method for kin-cohort analysis to infer the genotypes in the relatives. We also evaluated APC I1307K in 75 Ashkenazi Jewish individuals with prostate cancer from 27 families enrolled in the University of Michigan Prostate Cancer Genetic Study. We observed a slightly increased risk of prostate cancer in relatives of APC I1307K carriers, however, this difference was not statistically significant (hazard ratio, 1.6; 95% confidence intervals, 0.7-3.4). Similar conclusions were drawn using both techniques for kin-cohort analysis. APC I1307K was found in 7.4% of families genotyped, which is slightly higher than the allele prevalence reported in Ashkenazi Jews in the general population. In addition, we did not observe loss of heterozygosity at APC or a somatic mutation near APC I1307K using microdissected tumor DNA from mutation carriers enrolled in the Prostate Cancer Genetic Study. Overall, the evidence for an association between APC I1307K and prostate cancer is not compelling. APC I1307K is unlikely to play a clinically meaningful role in susceptibility to prostate cancer.

EGF gene polymorphism and the risk of incident primary melanoma.

Overexpression of the epidermal growth factor (EGF) pathway has been implicated in melanoma pathogenesis, and a recent case-control study identified a single nucleotide polymorphism (G to A) in the EGF gene where the G allele was associated with increased EGF expression and an increased risk of melanoma. To further evaluate this association, we conducted a case-control analysis from the Genes, Environment, and Melanoma study at the University of Michigan site using two different study designs. Incident cases of histopathologically confirmed first primary melanoma that were diagnosed between January 1, 2000 and December 31, 2000 from the University of Michigan Melanoma Clinic (n = 330) were compared with the following two different sources of nonmelanoma controls: spouse/friend controls (n = 84) and healthy volunteer controls from a case-control study of psoriasis (n = 148). Using a second analytic design, comparisons between multiple primary melanoma cases (n = 62) and single primary melanoma cases (n = 330) were also evaluated to estimate odds ratios (ORs). Genotyping for the single nucleotide substitution (G to A) at position 61 in the 5' untranslated region of the EGF gene was performed from genomic DNA, and epidemiological risk factors were assessed through a telephone interview. When EGF genotypes were compared between incident primary melanoma cases and the nonmelanoma controls, the risk associated with the homozygous G/G genotype was not statistically significantly associated with an increased risk for incident primary melanoma compared with the homozygous A/A genotype [OR, 1.09; 95% confidence interval (CI); 0.65-1.85]. No strong associations with EGF G/G genotype were observed in comparisons of multiple primary and single primary melanoma cases (OR, 0.66; 95% CI; 0.25-1.73). Case subjects with tumors >/=3.5 mm compared with those <3.5 mm were not significantly associated with the G/G genotype (OR, 0.54; 95% CI; 0.12-2.35). Our data do not support a significant association between melanoma and the EGF 61*G allele or the homozygous G/G genotype. The EGF polymorphism is not a reproducible risk factor for melanoma or thick melanoma in our data. The two analytic approaches used in the study provide evidence against a strong association between EGF 61*G and melanoma and demonstrate the potential utility of case-case designs for evaluating the role of single nucleotide polymorphisms and cancer. Additional independent studies will be required to elucidate relationships between genetic variation in the EGF gene and risk of melanoma.

Tumor-Infiltrating Lymphocytes, Crohn's-Like Lymphoid Reaction, and Survival From Colorectal Cancer.

BACKGROUND: While clinical outcomes from colorectal cancer (CRC) are influenced by stage at diagnosis and treatment, mounting evidence suggests that an enhanced lymphocytic reaction to a tumor may also be an informative prognostic indicator. METHODS: The roles of intratumoral T lymphocyte infiltration (TIL), peritumoral Crohn's-like lymphoid reaction (CLR), microsatellite instability (MSI), and clinicopathological characteristics in survival from CRC were examined using 2369 incident CRCs from a population-based case-control study in northern Israel. Cox proportional hazards regression was used to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) for CRC-specific and all-cause mortality in multivariable models adjusted for age, sex, ethnicity, grade, stage, and MSI. All statistical tests were two-sided. RESULTS: Tumors with TIL/high-powered field (HPF) of 2 or greater were associated with a statistically significant increase in CRC-specific (P < .001) and overall survival (P < .001) compared with tumors with TIL/HPF of less than 2. Similarly, tumors with a prominent CLR experienced better CRC-specific (P < .001) and overall survival (P < .001) as compared with those with no response. High TILs (HR = 0.76, 95% CI = 0.64 to 0.89, P < .001) and a prominent CLR (HR = 0.71, 95% CI = 0.62 to 0.80, P < .001), but not MSI, were associated with a statistically significant reduction in all-cause mortality after adjustment for established prognostic factors. CONCLUSIONS: TILs and CLR are both prognostic indicators for CRC after adjusting for traditional prognostic indicators.