Wild-type APC predicts poor prognosis in microsatellite-stable proximal colon cancer.

BACKGROUND: APC mutations (APC-mt) occur in ∼70% of colorectal cancers (CRCs), but their relationship to prognosis is unclear. METHODS: APC prognostic value was evaluated in 746 stage I-IV CRC patients, stratifying for tumour location and microsatellite instability (MSI). Microarrays were used to identify a gene signature that could classify APC mutation status, and classifier ability to predict prognosis was examined in an independent cohort. RESULTS: Wild-type APC microsatellite stable (APC-wt/MSS) tumours from the proximal colon showed poorer overall and recurrence-free survival (OS, RFS) than APC-mt/MSS proximal, APC-wt/MSS distal and APC-mt/MSS distal tumours (OS HR⩾1.79, P⩽0.015; RFS HR⩾1.88, P⩽0.026). APC was a stronger prognostic indicator than BRAF, KRAS, PIK3CA, TP53, CpG island methylator phenotype or chromosomal instability status (P⩽0.036). Microarray analysis similarly revealed poorer survival in MSS proximal cancers with an APC-wt-like signature (P=0.019). APC status did not affect outcomes in MSI tumours. In a validation on 206 patients with proximal colon cancer, APC-wt-like signature MSS cases showed poorer survival than APC-mt-like signature MSS or MSI cases (OS HR⩾2.50, P⩽0.010; RFS HR⩾2.14, P⩽0.025). Poor prognosis APC-wt/MSS proximal tumours exhibited features of the sessile serrated neoplasia pathway (P⩽0.016). CONCLUSIONS: APC-wt status is a marker of poor prognosis in MSS proximal colon cancer.

SISPA-Seq for rapid whole genome surveys of bacterial isolates.

Whole genome sequencing (WGS) of large isolate collections has many applications, yet sequencing costs are still significant. We sought to develop a rapid and cost efficient WGS method to address fundamental questions in clinical microbiology. We evaluated the performance of SISPA (Sequence-Independent, Single-Primer Amplification) combined with next-generation sequencing (SISPA-Seq) of 75 clinical isolates of Acinetobacter baumannii to establish whether SISPA-Seq resulted in sufficient coverage and quality to (1) determine strain phylogenetic placement and (2) and carriage of known antibiotic resistance (AbR) genes. Strains for which whole genome sequences were available were included for validation. Two libraries for each strain were constructed from separate SISPA reactions with different barcoded primers, using genomic DNA prepared from either high quality or rapid heat-lysis preparations. SISPA-Seq resulted in a median of 65× genome coverage when reads from both primer sets were combined. Coverage and quality were sufficient for detection of AbR genes by comparison of reads to the ARG-ANNOT database and were often sufficient to distinguish between different allelic variants of the same gene. kSNP and RAxML were used to construct a robust phylogeny based on single-nucleotide variants (SNVs) that showed that the SISPA-Seq data was sufficient for sensitive and accurate phylogenetic placement. Advantages of the SISPA-Seq method include inexpensive and rapid DNA preparation and a typical total cost less than one-half that of standard genome sequencing. In summary, SISPA-Seq can be used to survey whole genomes of a large strain collection and identify strains that should be targeted for additional sequencing.

PIK3CA and PTEN gene and exon mutation-specific clinicopathologic and molecular associations in colorectal cancer.

PURPOSE: PIK3CA and PTEN mutations are prevalent in colorectal cancer and potential markers of response to mitogen-activated protein/extracellular signal-regulated kinase inhibitors and anti-EGF receptor antibody therapy. Relationships between phosphoinositide 3-kinase (PI3K) pathway mutation, clinicopathologic characteristics, molecular features, and prognosis remain controversial. EXPERIMENTAL DESIGN: A total of 1,093 stage I-IV colorectal cancers were screened for PIK3CA (exons 9 and 20), KRAS (codons 12-13), BRAF (codon 600) mutations, and microsatellite instability (MSI). PTEN (exons 3-8) and CpG island methylator phenotype (CIMP) status were determined in 744 and 489 cases. PIK3CA data were integrated with 17 previous reports (n = 5,594). RESULTS: PIK3CA and PTEN mutations were identified in 11.9% and 5.8% of colorectal cancers. PTEN mutation was associated with proximal tumors, mucinous histology, MSI-high (MSI-H), CIMP-high (CIMP-H), and BRAF mutation (P < 0.02). PIK3CA mutation was related to older age, proximal tumors, mucinous histology, and KRAS mutation (P < 0.04). In integrated cohort analysis, PIK3CA exon 9 and 20 mutations were overrepresented in proximal, CIMP-low (CIMP-L), and KRAS-mutated cancers (P ≤ 0.011). Comparing PIK3CA exonic mutants, exon 20 mutation was associated with MSI-H, CIMP-H, and BRAF mutation, and exon 9 mutation was associated with KRAS mutation (P ≤ 0.027). Disease-free survival for stage II/III colorectal cancers did not differ by PI3K pathway status. CONCLUSION: PI3K pathway mutation is prominent in proximal colon cancers, with PIK3CA exon 20 and PTEN mutations associated with features of the sessile-serrated pathway (MSI-H/CIMP-H/BRAF(mut)), and PIK3CA exon 9 (and to a lesser extent exon 20) mutation associated with features of the traditional serrated pathway (CIMP-L/KRAS(mut)) of tumorigenesis. Our data highlight the PI3K pathway as a therapeutic target in distinct colorectal cancer subtypes.

SMAD2, SMAD3 and SMAD4 mutations in colorectal cancer.

Activation of the canonical TGF-ß signaling pathway provides growth inhibitory signals in the normal intestinal epithelium. Colorectal cancers (CRCs) frequently harbor somatic mutations in the pathway members TGFBR2 and SMAD4, but to what extent mutations in SMAD2 or SMAD3 contribute to tumorigenesis is unclear. A cohort of 744 primary CRCs and 36 CRC cell lines were sequenced for SMAD4, SMAD2, and SMAD3 and analyzed for allelic loss by single-nucleotide polymorphism (SNP) microarray analysis. Mutation spectra were compared between the genes, the pathogenicity of mutations was assessed, and relationships with clinicopathologic features were examined. The prevalence of SMAD4, SMAD2, and SMAD3 mutations in sporadic CRCs was 8.6% (64 of 744), 3.4% (25 of 744), and 4.3% (32 of 744), respectively. A significant overrepresentation of two genetic hits was detected for SMAD4 and SMAD3, consistent with these genes acting as tumor suppressors. SMAD4 mutations were associated with mucinous histology. The mutation spectra of SMAD2 and SMAD3 were highly similar to that of SMAD4, both in mutation type and location within the encoded proteins. In silico analyses suggested the majority of the mutations were pathogenic, with most missense changes predicted to reduce protein stability or hinder SMAD complex formation. The latter altered interface residues or disrupted the phosphorylation-regulated Ser-Ser-X-Ser motifs within SMAD2 and SMAD3. Functional analyses of selected mutations showed reductions in SMAD3 transcriptional activity and SMAD2-SMAD4 complex formation. Joint biallelic hits in SMAD2 and SMAD3 were overrepresented and mutually exclusive to SMAD4 mutation, underlining the critical roles of these three proteins within the TGF-ß signaling pathway.

KRAS mutation is associated with lung metastasis in patients with curatively resected colorectal cancer.

PURPOSE: Oncogene mutations contribute to colorectal cancer development. We searched for differences in oncogene mutation profiles between colorectal cancer metastases from different sites and evaluated these as markers for site of relapse. EXPERIMENTAL DESIGN: One hundred colorectal cancer metastases were screened for mutations in 19 oncogenes, and further 61 metastases and 87 matched primary cancers were analyzed for genes with identified mutations. Mutation prevalence was compared between (a) metastases from liver (n = 65), lung (n = 50), and brain (n = 46), (b) metastases and matched primary cancers, and (c) metastases and an independent cohort of primary cancers (n = 604). Mutations differing between metastasis sites were evaluated as markers for site of relapse in 859 patients from the VICTOR trial. RESULTS: In colorectal cancer metastases, mutations were detected in 4 of 19 oncogenes: BRAF (3.1%), KRAS (48.4%), NRAS (6.2%), and PIK3CA (16.1%). KRAS mutation prevalence was significantly higher in lung (62.0%) and brain (56.5%) than in liver metastases (32.3%; P = 0.003). Mutation status was highly concordant between primary cancer and metastasis from the same individual. Compared with independent primary cancers, KRAS mutations were more common in lung and brain metastases (P < 0.005), but similar in liver metastases. Correspondingly, KRAS mutation was associated with lung relapse (HR = 2.1; 95% CI, 1.2 to 3.5, P = 0.007) but not liver relapse in patients from the VICTOR trial. CONCLUSIONS: KRAS mutation seems to be associated with metastasis in specific sites, lung and brain, in colorectal cancer patients. Our data highlight the potential of somatic mutations for informing surveillance strategies.

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

BACKGROUND: We are developing a cross-species comparison strategy to distinguish between cancer driver- and passenger gene alteration candidates, by utilizing the difference in genomic location of orthologous genes between the human and other mammals. As an initial test of this strategy, we conducted a pilot study with human colorectal cancer (CRC) and its mouse model C57BL/6J ApcMin/+, focusing on human 5q22.2 and 18q21.1-q21.2. METHODS: We first performed bioinformatics analysis on the evolution of 5q22.2 and 18q21.1-q21.2 regions. Then, we performed exon-targeted sequencing, real time quantitative polymerase chain reaction (qPCR), and real time quantitative reverse transcriptase PCR (qRT-PCR) analyses on a number of genes of both regions with both human and mouse colon tumors. RESULTS: These two regions (5q22.2 and 18q21.1-q21.2) are frequently deleted in human CRCs and encode genuine colorectal tumor suppressors APC and SMAD4. They also encode genes such as MCC (mutated in colorectal cancer) with their role in CRC etiology unknown. We have discovered that both regions are evolutionarily unstable, resulting in genes that are clustered in each human region being found scattered at several distinct loci in the genome of many other species. For instance, APC and MCC are within 200 kb apart in human 5q22.2 but are 10 Mb apart in the mouse genome. Importantly, our analyses revealed that, while known CRC driver genes APC and SMAD4 were disrupted in both human colorectal tumors and tumors from ApcMin/+ mice, the questionable MCC gene was disrupted in human tumors but appeared to be intact in mouse tumors. CONCLUSIONS: These results indicate that MCC may not actually play any causative role in early colorectal tumorigenesis. We also hypothesize that its disruption in human CRCs is likely a mere result of its close proximity to APC in the human genome. Expanding this pilot study to the entire genome may identify more questionable genes like MCC, facilitating the discovery of new CRC driver gene candidates.

Genome sequence of the necrotrophic plant pathogen Pythium ultimum reveals original pathogenicity mechanisms and effector repertoire.

BACKGROUND: Pythium ultimum is a ubiquitous oomycete plant pathogen responsible for a variety of diseases on a broad range of crop and ornamental species. RESULTS: The P. ultimum genome (42.8 Mb) encodes 15,290 genes and has extensive sequence similarity and synteny with related Phytophthora species, including the potato blight pathogen Phytophthora infestans. Whole transcriptome sequencing revealed expression of 86% of genes, with detectable differential expression of suites of genes under abiotic stress and in the presence of a host. The predicted proteome includes a large repertoire of proteins involved in plant pathogen interactions, although, surprisingly, the P. ultimum genome does not encode any classical RXLR effectors and relatively few Crinkler genes in comparison to related phytopathogenic oomycetes. A lower number of enzymes involved in carbohydrate metabolism were present compared to Phytophthora species, with the notable absence of cutinases, suggesting a significant difference in virulence mechanisms between P. ultimum and more host-specific oomycete species. Although we observed a high degree of orthology with Phytophthora genomes, there were novel features of the P. ultimum proteome, including an expansion of genes involved in proteolysis and genes unique to Pythium. We identified a small gene family of cadherins, proteins involved in cell adhesion, the first report of these in a genome outside the metazoans. CONCLUSIONS: Access to the P. ultimum genome has revealed not only core pathogenic mechanisms within the oomycetes but also lineage-specific genes associated with the alternative virulence and lifestyles found within the pythiaceous lineages compared to the Peronosporaceae.

The dynamic genome of Hydra.

The freshwater cnidarian Hydra was first described in 1702 and has been the object of study for 300 years. Experimental studies of Hydra between 1736 and 1744 culminated in the discovery of asexual reproduction of an animal by budding, the first description of regeneration in an animal, and successful transplantation of tissue between animals. Today, Hydra is an important model for studies of axial patterning, stem cell biology and regeneration. Here we report the genome of Hydra magnipapillata and compare it to the genomes of the anthozoan Nematostella vectensis and other animals. The Hydra genome has been shaped by bursts of transposable element expansion, horizontal gene transfer, trans-splicing, and simplification of gene structure and gene content that parallel simplification of the Hydra life cycle. We also report the sequence of the genome of a novel bacterium stably associated with H. magnipapillata. Comparisons of the Hydra genome to the genomes of other animals shed light on the evolution of epithelia, contractile tissues, developmentally regulated transcription factors, the Spemann-Mangold organizer, pluripotency genes and the neuromuscular junction.

Transcriptome-guided characterization of genomic rearrangements in a breast cancer cell line.

We have identified new genomic alterations in the breast cancer cell line HCC1954, using high-throughput transcriptome sequencing. With 120 Mb of cDNA sequences, we were able to identify genomic rearrangement events leading to fusions or truncations of genes including MRE11 and NSD1, genes already implicated in oncogenesis, and 7 rearrangements involving other additional genes. This approach demonstrates that high-throughput transcriptome sequencing is an effective strategy for the characterization of genomic rearrangements in cancers.

An integrated genomic analysis of human glioblastoma multiforme.

Glioblastoma multiforme (GBM) is the most common and lethal type of brain cancer. To identify the genetic alterations in GBMs, we sequenced 20,661 protein coding genes, determined the presence of amplifications and deletions using high-density oligonucleotide arrays, and performed gene expression analyses using next-generation sequencing technologies in 22 human tumor samples. This comprehensive analysis led to the discovery of a variety of genes that were not known to be altered in GBMs. Most notably, we found recurrent mutations in the active site of isocitrate dehydrogenase 1 (IDH1) in 12% of GBM patients. Mutations in IDH1 occurred in a large fraction of young patients and in most patients with secondary GBMs and were associated with an increase in overall survival. These studies demonstrate the value of unbiased genomic analyses in the characterization of human brain cancer and identify a potentially useful genetic alteration for the classification and targeted therapy of GBMs.

Genetic variation in an individual human exome.

There is much interest in characterizing the variation in a human individual, because this may elucidate what contributes significantly to a person's phenotype, thereby enabling personalized genomics. We focus here on the variants in a person's 'exome,' which is the set of exons in a genome, because the exome is believed to harbor much of the functional variation. We provide an analysis of the approximately 12,500 variants that affect the protein coding portion of an individual's genome. We identified approximately 10,400 nonsynonymous single nucleotide polymorphisms (nsSNPs) in this individual, of which approximately 15-20% are rare in the human population. We predict approximately 1,500 nsSNPs affect protein function and these tend be heterozygous, rare, or novel. Of the approximately 700 coding indels, approximately half tend to have lengths that are a multiple of three, which causes insertions/deletions of amino acids in the corresponding protein, rather than introducing frameshifts. Coding indels also occur frequently at the termini of genes, so even if an indel causes a frameshift, an alternative start or stop site in the gene can still be used to make a functional protein. In summary, we reduced the set of approximately 12,500 nonsilent coding variants by approximately 8-fold to a set of variants that are most likely to have major effects on their proteins' functions. This is our first glimpse of an individual's exome and a snapshot of the current state of personalized genomics. The majority of coding variants in this individual are common and appear to be functionally neutral. Our results also indicate that some variants can be used to improve the current NCBI human reference genome. As more genomes are sequenced, many rare variants and non-SNP variants will be discovered. We present an approach to analyze the coding variation in humans by proposing multiple bioinformatic methods to hone in on possible functional variation.

Novel computational methods for increasing PCR primer design effectiveness in directed sequencing.

BACKGROUND: Polymerase chain reaction (PCR) is used in directed sequencing for the discovery of novel polymorphisms. As the first step in PCR directed sequencing, effective PCR primer design is crucial for obtaining high-quality sequence data for target regions. Since current computational primer design tools are not fully tuned with stable underlying laboratory protocols, researchers may still be forced to iteratively optimize protocols for failed amplifications after the primers have been ordered. Furthermore, potentially identifiable factors which contribute to PCR failures have yet to be elucidated. This inefficient approach to primer design is further intensified in a high-throughput laboratory, where hundreds of genes may be targeted in one experiment. RESULTS: We have developed a fully integrated computational PCR primer design pipeline that plays a key role in our high-throughput directed sequencing pipeline. Investigators may specify target regions defined through a rich set of descriptors, such as Ensembl accessions and arbitrary genomic coordinates. Primer pairs are then selected computationally to produce a minimal amplicon set capable of tiling across the specified target regions. As part of the tiling process, primer pairs are computationally screened to meet the criteria for success with one of two PCR amplification protocols. In the process of improving our sequencing success rate, which currently exceeds 95% for exons, we have discovered novel and accurate computational methods capable of identifying primers that may lead to PCR failures. We reveal the laboratory protocols and their associated, empirically determined computational parameters, as well as describe the novel computational methods which may benefit others in future primer design research. CONCLUSION: The high-throughput PCR primer design pipeline has been very successful in providing the basis for high-quality directed sequencing results and for minimizing costs associated with labor and reprocessing. The modular architecture of the primer design software has made it possible to readily integrate additional primer critique tests based on iterative feedback from the laboratory. As a result, the primer design software, coupled with the laboratory protocols, serves as a powerful tool for low and high-throughput primer design to enable successful directed sequencing.

The diploid genome sequence of an individual human.

Presented here is a genome sequence of an individual human. It was produced from approximately 32 million random DNA fragments, sequenced by Sanger dideoxy technology and assembled into 4,528 scaffolds, comprising 2,810 million bases (Mb) of contiguous sequence with approximately 7.5-fold coverage for any given region. We developed a modified version of the Celera assembler to facilitate the identification and comparison of alternate alleles within this individual diploid genome. Comparison of this genome and the National Center for Biotechnology Information human reference assembly revealed more than 4.1 million DNA variants, encompassing 12.3 Mb. These variants (of which 1,288,319 were novel) included 3,213,401 single nucleotide polymorphisms (SNPs), 53,823 block substitutions (2-206 bp), 292,102 heterozygous insertion/deletion events (indels)(1-571 bp), 559,473 homozygous indels (1-82,711 bp), 90 inversions, as well as numerous segmental duplications and copy number variation regions. Non-SNP DNA variation accounts for 22% of all events identified in the donor, however they involve 74% of all variant bases. This suggests an important role for non-SNP genetic alterations in defining the diploid genome structure. Moreover, 44% of genes were heterozygous for one or more variants. Using a novel haplotype assembly strategy, we were able to span 1.5 Gb of genome sequence in segments >200 kb, providing further precision to the diploid nature of the genome. These data depict a definitive molecular portrait of a diploid human genome that provides a starting point for future genome comparisons and enables an era of individualized genomic information.

A Sanger/pyrosequencing hybrid approach for the generation of high-quality draft assemblies of marine microbial genomes.

Since its introduction a decade ago, whole-genome shotgun sequencing (WGS) has been the main approach for producing cost-effective and high-quality genome sequence data. Until now, the Sanger sequencing technology that has served as a platform for WGS has not been truly challenged by emerging technologies. The recent introduction of the pyrosequencing-based 454 sequencing platform (454 Life Sciences, Branford, CT) offers a very promising sequencing technology alternative for incorporation in WGS. In this study, we evaluated the utility and cost-effectiveness of a hybrid sequencing approach using 3730xl Sanger data and 454 data to generate higher-quality lower-cost assemblies of microbial genomes compared to current Sanger sequencing strategies alone.

Sequence survey of receptor tyrosine kinases reveals mutations in glioblastomas.

It is now clear that tyrosine kinases represent attractive targets for therapeutic intervention in cancer. Recent advances in DNA sequencing technology now provide the opportunity to survey mutational changes in cancer in a high-throughput and comprehensive manner. Here we report on the sequence analysis of members of the receptor tyrosine kinase (RTK) gene family in the genomes of glioblastoma brain tumors. Previous studies have identified a number of molecular alterations in glioblastoma, including amplification of the RTK epidermal growth factor receptor. We have identified mutations in two other RTKs: (i) fibroblast growth receptor 1, including the first mutations in the kinase domain in this gene observed in any cancer, and (ii) a frameshift mutation in the platelet-derived growth factor receptor-alpha gene. Fibroblast growth receptor 1, platelet-derived growth factor receptor-alpha, and epidermal growth factor receptor are all potential entry points to the phosphatidylinositol 3-kinase and mitogen-activated protein kinase intracellular signaling pathways already known to be important for neoplasia. Our results demonstrate the utility of applying DNA sequencing technology to systematically assess the coding sequence of genes within cancer genomes.