Validation of predictive models for germline mutations in DNA mismatch repair genes in colorectal cancer.

Lynch syndrome is defined by the presence of germline mutations in mismatch repair (MMR) genes. Several models have been recently devised that predict mutation carrier status (Myriad Genetics, Wijnen, Barnetson, PREMM and MMRpro models). Families at moderate-high risk for harboring a Lynch-associated mutation, referred to the BC Cancer Agency (BCCA) Hereditary Cancer Program (HCP), underwent mutation analysis, immunohistochemistry and/or microsatellite testing. Seventy-two tested cases were included. Twenty-five patients were mutation positive (34.7%) and 47 were mutation negative (65.3%). Nineteen of 43 patients who were both microsatellite stable and normal on immunohistochemistry for MLH1 and MSH2 were also genotyped for mutations in these genes; all 19 were negative for MMR gene mutations. Model-derived probabilities of harboring a MMR gene mutation in the proband were calculated and compared to observed results. The area under the ROC curves were 0.75 (95%CI; 0.63-0.87), 0.86 (0.7-0.96), 0.89 (0.82-0.97), 0.89 (0.81-0.98) and 0.93 (0.86-0.99) for the Myriad, Barnetson, Wijnen, MMRpro and PREMM models, respectively. The Amsterdam II criteria had a sensitivity and specificity of 0.76 and 0.74, respectively, in this cohort. The PREMM model demonstrated the best performance for predicting carrier status based on the positive likelihood ratios at the >10%, >20% and >30% probability thresholds. In this referred cohort, the PREMM model had the most favorable concordance index and predictive performance for carrier status based on the positive LR. These prediction models (PREMM, MMRPro and Wijnen) may soon replace the Amsterdam II and revised Bethesda criteria as a prescreening tool for Lynch mutations.

De novo transcriptome assembly with ABySS.

MOTIVATION: Whole transcriptome shotgun sequencing data from non-normalized samples offer unique opportunities to study the metabolic states of organisms. One can deduce gene expression levels using sequence coverage as a surrogate, identify coding changes or discover novel isoforms or transcripts. Especially for discovery of novel events, de novo assembly of transcriptomes is desirable. RESULTS: Transcriptome from tumor tissue of a patient with follicular lymphoma was sequenced with 36 base pair (bp) single- and paired-end reads on the Illumina Genome Analyzer II platform. We assembled approximately 194 million reads using ABySS into 66 921 contigs 100 bp or longer, with a maximum contig length of 10 951 bp, representing over 30 million base pairs of unique transcriptome sequence, or roughly 1% of the genome. AVAILABILITY AND IMPLEMENTATION: Source code and binaries of ABySS are freely available for download at http://www.bcgsc.ca/platform/bioinfo/software/abyss. Assembler tool is implemented in C++. The parallel version uses Open MPI. ABySS-Explorer tool is implemented in Java using the Java universal network/graph framework. CONTACT: ibirol@bcgsc.ca.

Somatic mutations altering EZH2 (Tyr641) in follicular and diffuse large B-cell lymphomas of germinal-center origin.

Follicular lymphoma (FL) and the GCB subtype of diffuse large B-cell lymphoma (DLBCL) derive from germinal center B cells. Targeted resequencing studies have revealed mutations in various genes encoding proteins in the NF-kappaB pathway that contribute to the activated B-cell (ABC) DLBCL subtype, but thus far few GCB-specific mutations have been identified. Here we report recurrent somatic mutations affecting the polycomb-group oncogene EZH2, which encodes a histone methyltransferase responsible for trimethylating Lys27 of histone H3 (H3K27). After the recent discovery of mutations in KDM6A (UTX), which encodes the histone H3K27me3 demethylase UTX, in several cancer types, EZH2 is the second histone methyltransferase gene found to be mutated in cancer. These mutations, which result in the replacement of a single tyrosine in the SET domain of the EZH2 protein (Tyr641), occur in 21.7% of GCB DLBCLs and 7.2% of FLs and are absent from ABC DLBCLs. Our data are consistent with the notion that EZH2 proteins with mutant Tyr641 have reduced enzymatic activity in vitro.

Cytogenetically balanced translocations are associated with focal copy number alterations.

Current cytogenetic methods (e.g., G-banding and multicolor chromosomal painting) allow detection of translocation events but lack the resolution to (a) locate the breakpoints precisely at the chromosome band level or (b) discriminate balanced translocations from translocations with copy number alterations not previously reported, or imperfectly balanced translocations. In this study, we demonstrate that cytogenetically balanced translocations are in fact frequently associated with segmental gain or loss of DNA. The recent development of a whole genome tiling path BAC array has enabled tiling resolution analysis of genomic segmental copy number status. Combining tiling resolution BAC array comparative genomic hybridization (array CGH) with G-Banding analysis and multicolor chromosomal painting approaches such as spectral karyotyping (SKY) facilitates high-resolution mapping of genomic alterations associated with imperfectly balanced translocations. Using a refined version of our CGH array we have deduced the copy number status throughout the genomes of three cytogenetically well-characterized prostate cancer cell lines (PC3, DU145, LNCaP) to determine whether translocations are associated with focal gains and losses of DNA. At 78 kb tiling resolution we identified the boundaries of 170, 80, and 34 known and novel copy number alterations (CNA) in these cell line genomes, respectively. Thirty-three of the 36 known translocations (92%, P < 0.001) in DU145 were associated with segmental CNA. Likewise, 80% (P < 0.001) of the known translocations showed association in LNCaP. Although many translocation breakpoints exhibit segmental alteration in PC3, the pattern of chromosomal rearrangements is too complex for use in comprehensive association with CNA boundaries. Our results reveal that imperfectly balanced translocations in tumor genomes are a phenomenon that occurs at frequencies much higher than previously demonstrated.

A fluorescence in situ hybridization study of ETV6-NTRK3 fusion gene in secretory breast carcinoma.

The translocation t(12;15)(p13;q25), in which the ETV6 gene from chromosome 12 is rearranged with the NTRK3 gene from chromosome 15, has recently been identified in secretory breast carcinoma (SBC). This fusion gene was initially described in congenital fibrosarcoma and congenital mesoblastic nephroma. The biological consequence of this translocation is the expression of a chimeric protein tyrosine kinase with potent transforming activity. To assess the frequency of t(12;15)(p13;q25) in breast cancer, we developed complementary probe sets (fusion and split-apart probes) for the detection of this translocation by fluorescence in situ hybridization (FISH) in paraffin-embedded, formalin-fixed tissue sections. We tested four histologically confirmed cases of SBC for the presence of the ETV6-NTRK3 gene fusion and then applied the FISH assay to tissue microarrays (TMAs) in order to screen 481 cases of formalin-fixed, paraffin-embedded invasive breast carcinomas of various histologic subtypes. Three of the four cases of SBC revealed fusion signals. Of the 481 cases in the TMAs, 202 gave signals of sufficient quality for screening by FISH, and only one case showed fusion signals in most or all of the tumor cells. On review of the histology of this case, SBC was confirmed. On the other hand, none of the fusion-negative breast cancers revealed SBC histology. In all cases, the results from the fusion and split-apart FISH assays for the ETV6-NTRK3 fusion genes were concordant. Our data suggest that the ETV6-NTRK3 fusion gene is a specific genetic alteration in SBC.