Jorge Zamora

Exome sequencing identifies recurrent mutations of the splicing factor SF3B1 gene in chronic lymphocytic leukemia

Here we perform whole-exome sequencing of samples from 105 individuals with chronic lymphocytic leukemia (CLL), the most frequent leukemia in adults in Western countries. We found 1,246 somatic mutations potentially affecting gene function and identified 78 genes with predicted functional alterations in more than one tumor sample. Among these genes, SF3B1, encoding a subunit of the spliceosomal U2 small nuclear ribonucleoprotein (snRNP), is somatically mutated in 9.7% of affected individuals. Further analysis in 279 individuals with CLL showed that SF3B1 mutations were associated with faster disease progression and poor overall survival. This work provides the first comprehensive catalog of somatic mutations in CLL with relevant clinical correlates and defines a large set of new genes that may drive the development of this common form of leukemia. The results reinforce the idea that targeting several well-known genetic pathways, including mRNA splicing, could be useful in the treatment of CLL and other malignancies.

Analysis of the genetic phylogeny of multifocal prostate cancer identifies multiple independent clonal expansions in neoplastic and morphologically normal prostate tissue

Genome-wide DNA sequencing was used to decrypt the phylogeny of multiple samples from distinct areas of cancer and morphologically normal tissue taken from the prostates of three men. Mutations were present at high levels in morphologically normal tissue distant from the cancer, reflecting clonal expansions, and the underlying mutational processes at work in morphologically normal tissue were also at work in cancer. Our observations demonstrate the existence of ongoing abnormal mutational processes, consistent with field effects, underlying carcinogenesis. This mechanism gives rise to extensive branching evolution and cancer clone mixing, as exemplified by the coexistence of multiple cancer lineages harboring distinct ERG fusions within a single cancer nodule. Subsets of mutations were shared either by morphologically normal and malignant tissues or between different ERG lineages, indicating earlier or separate clonal cell expansions. Our observations inform on the origin of multifocal disease and have implications for prostate cancer therapy in individual cases.

Extensive transduction of nonrepetitive DNA mediated by L1 retrotransposition in cancer genomes

Introduction The human genome is peppered with mobile repetitive elements called long interspersed nuclear element–1 (L1) retrotransposons. Propagating through RNA and cDNA intermediates, these molecular parasites copy and insert themselves throughout the genome, with potentially disruptive effects on neighboring genes or regulatory sequences. In the germ line, unique sequence downstream of L1 elements can also be retrotransposed if transcription continues beyond the repeat, a process known as 3′ transduction. There has been growing interest in retrotransposition and 3′ transduction as a possible source of somatic mutations during tumorigenesis. The activity of individual L1 elements fluctuates during tumor evolution. In a lung tumor, hundreds of 3′ transductions arose from a small number of active L1 source elements (colored circles on outer rim of circle). As the tumor evolved from the preinvasive common ancestor to invasive cancer, individual elements exhibited variable activity over time. Rationale To explore whether 3′ transductions are frequent in cancer, we developed a bioinformatic algorithm for identifying somatically acquired retrotranspositions in cancer genomes. We applied our algorithm to 290 cancer samples from 244 patients across 12 tumor types. The unique downstream sequence mobilized with 3′ transductions effectively fingerprints the L1 source element, providing insights into the activity of individual L1 loci across the genome. Results Across the 290 samples, we identified 2756 somatic L1 retrotranspositions. Tumors from 53% of patients had at least one such event, with colorectal and lung cancers being most frequently affected (93% and 75% of patients, respectively). Somatic 3′ transductions comprised 24% of events, half of which represented mobilizations of unique sequence alone, without any accompanying L1 sequence. Overall, 95% of 3′ transductions identified derived from only 72 germline L1 source elements, with as few as four loci accounting for 50% of events. In a given sample, the same source element could generate 50 or more somatic transductions, scattered extensively across the genome. About 5% of somatic transductions arose from L1 source elements that were themselves somatic retrotranspositions. In three of the cases in which we sequenced more than one sample from a patient’s tumor, we were able to place 3′ transductions on the phylogenetic tree. We found that the activity of individual source elements fluctuated during tumor evolution, with different subclones exhibiting much variability in which elements were “on” and which were “off.” The ability to identify the individual L1 source elements active in a given tumor enabled us to study the promoter methylation of those elements specifically. We found that 3′ transduction activity in a patient’s tumor was always associated with hypomethylation of that element. Overall, 2.3% of transductions distributed exons or entire genes to other sites in the genome, and many more mobilized deoxyribonuclease I (DNAse-I) hypersensitive sites or transcription factor binding sites identified by the ENCODE project. Occasionally, somatic L1 insertions inserted near coding sequence and redistributed these exons elsewhere in the genome. However, we found no general effects of retrotranspositions on transcription levels of genes at the insertion points and no evidence for aberrant RNA species resulting from somatically acquired transposable elements. Indeed, as with germline retrotranspositions, somatic insertions exhibited a strong enrichment in heterochromatic, gene-poor regions of the genome. Conclusion Somatic 3′ transduction occurs frequently in human tumors, and in some cases transduction events can scatter exons, genes, and regulatory elements widely across the genome. Dissemination of these sequences appears to be due to a small number of highly active L1 elements, whose activity can wax and wane during tumor evolution. The majority of the retrotransposition events are likely to be harmless “passenger” mutations. Hitchhiking through the tumor genome Retrotransposons are DNA repeat sequences that are constantly on the move. By poaching certain cellular enzymes, they copy and insert themselves at new sites in the genome. Sometimes they carry along adjacent DNA sequences, a process called 3′ transduction. Tubio et al. found that 3′ transduction is a common event in human tumors. Because this process can scatter genes and regulatory sequences across the genome, it may represent yet another mechanism by which tumor cells acquire new mutations that help them survive and grow. Science, this issue p. 10.1126/science.1251343 Tumor genomes are peppered with mobile repeat sequences that carry along adjacent DNA when they insert into new genomic sites. Long interspersed nuclear element–1 (L1) retrotransposons are mobile repetitive elements that are abundant in the human genome. L1 elements propagate through RNA intermediates. In the germ line, neighboring, nonrepetitive sequences are occasionally mobilized by the L1 machinery, a process called 3′ transduction. Because 3′ transductions are potentially mutagenic, we explored the extent to which they occur somatically during tumorigenesis. Studying cancer genomes from 244 patients, we found that tumors from 53% of the patients had somatic retrotranspositions, of which 24% were 3′ transductions. Fingerprinting of donor L1s revealed that a handful of source L1 elements in a tumor can spawn from tens to hundreds of 3′ transductions, which can themselves seed further retrotranspositions. The activity of individual L1 elements fluctuated during tumor evolution and correlated with L1 promoter hypomethylation. The 3′ transductions disseminated genes, exons, and regulatory elements to new locations, most often to heterochromatic regions of the genome.

High-throughput and quantitative genome-wide messenger RNA sequencing for molecular phenotyping.

BackgroundWe present a genome-wide messenger RNA (mRNA) sequencing technique that converts small amounts of RNA from many samples into molecular phenotypes. It encompasses all steps from sample preparation to sequence analysis and is applicable to baseline profiling or perturbation measurements.ResultsMultiplex sequencing of transcript 3′ ends identifies differential transcript abundance independent of gene annotation. We show that increasing biological replicate number while maintaining the total amount of sequencing identifies more differentially abundant transcripts.ConclusionsThis method can be implemented on polyadenylated RNA from any organism with an annotated reference genome and in any laboratory with access to Illumina sequencing.

Appraising the relevance of DNA copy number loss and gain in prostate cancer using whole genome DNA sequence data.

A variety of models have been proposed to explain regions of recurrent somatic copy number alteration (SCNA) in human cancer. Our study employs Whole Genome DNA Sequence (WGS) data from tumor samples (n = 103) to comprehensively assess the role of the Knudson two hit genetic model in SCNA generation in prostate cancer. 64 recurrent regions of loss and gain were detected, of which 28 were novel, including regions of loss with more than 15% frequency at Chr4p15.2-p15.1 (15.53%), Chr6q27 (16.50%) and Chr18q12.3 (17.48%). Comprehensive mutation screens of genes, lincRNA encoding sequences, control regions and conserved domains within SCNAs demonstrated that a two-hit genetic model was supported in only a minor proportion of recurrent SCNA losses examined (15/40). We found that recurrent breakpoints and regions of inversion often occur within Knudson model SCNAs, leading to the identification of ZNF292 as a target gene for the deletion at 6q14.3-q15 and NKX3.1 as a two-hit target at 8p21.3-p21.2. The importance of alterations of lincRNA sequences was illustrated by the identification of a novel mutational hotspot at the KCCAT42, FENDRR, CAT1886 and STCAT2 loci at the 16q23.1-q24.3 loss. Our data confirm that the burden of SCNAs is predictive of biochemical recurrence, define nine individual regions that are associated with relapse, and highlight the possible importance of ion channel and G-protein coupled-receptor (GPCR) pathways in cancer development. We concluded that a two-hit genetic model accounts for about one third of SCNA indicating that mechanisms, such haploinsufficiency and epigenetic inactivation, account for the remaining SCNA losses.

Biological and prognostic impact of APOBEC-induced mutations in the spectrum of plasma cell dyscrasias and multiple myeloma cell lines

Biological and prognostic impact of APOBEC-induced mutations in the spectrum of plasma cell dyscrasias and multiple myeloma cell lines

Erratum: Analysis of the genetic phylogeny of multifocal prostate cancer identifies multiple independent clonal expansions in neoplastic and morphologically normal prostate tissue (Nature Genetics (2015) 47 (367-372))

In the version of this article initially published, the following two sentences were omitted from the Acknowledgments: “Sequencing was carried out at the Millard and Muriel Jacobs Genome Facility at the California Institute of Technology. This work was supported by US National Institutes of Health grants to P.W.S. (GM084389) and to R.V.A. (AI056189), by Cornell University salary and start-up funds to E.M.S. and by the Howard Hughes Medical Institute to P.W.S.” The error has been corrected in the HTML and PDF versions of the article.