Eric Stawiski

Recurrent R-spondin fusions in colon cancer

Identifying and understanding changes in cancer genomes is essential for the development of targeted therapeutics. Here we analyse systematically more than 70 pairs of primary human colon tumours by applying next-generation sequencing to characterize their exomes, transcriptomes and copy-number alterations. We have identified 36,303 protein-altering somatic changes that include several new recurrent mutations in the Wnt pathway gene TCF7L2, chromatin-remodelling genes such as TET2 and TET3 and receptor tyrosine kinases including ERBB3. Our analysis for significantly mutated cancer genes identified 23 candidates, including the cell cycle checkpoint kinase ATM. Copy-number and RNA-seq data analysis identified amplifications and corresponding overexpression of IGF2 in a subset of colon tumours. Furthermore, using RNA-seq data we identified multiple fusion transcripts including recurrent gene fusions involving R-spondin family members RSPO2 and RSPO3 that together occur in 10% of colon tumours. The RSPO fusions were mutually exclusive with APC mutations, indicating that they probably have a role in the activation of Wnt signalling and tumorigenesis. Consistent with this we show that the RSPO fusion proteins were capable of potentiating Wnt signalling. The R-spondin gene fusions and several other gene mutations identified in this study provide new potential opportunities for therapeutic intervention in colon cancer.

Comprehensive genomic analysis identifies SOX2 as a frequently amplified gene in small-cell lung cancer

Small-cell lung cancer (SCLC) is an exceptionally aggressive disease with poor prognosis. Here, we obtained exome, transcriptome and copy-number alteration data from approximately 53 samples consisting of 36 primary human SCLC and normal tissue pairs and 17 matched SCLC and lymphoblastoid cell lines. We also obtained data for 4 primary tumors and 23 SCLC cell lines. We identified 22 significantly mutated genes in SCLC, including genes encoding kinases, G protein–coupled receptors and chromatin-modifying proteins. We found that several members of the SOX family of genes were mutated in SCLC. We also found SOX2 amplification in ∼27% of the samples. Suppression of SOX2 using shRNAs blocked proliferation of SOX2-amplified SCLC lines. RNA sequencing identified multiple fusion transcripts and a recurrent RLF-MYCL1 fusion. Silencing of MYCL1 in SCLC cell lines that had the RLF-MYCL1 fusion decreased cell proliferation. These data provide an in-depth view of the spectrum of genomic alterations in SCLC and identify several potential targets for therapeutic intervention.

An Evolutionary Model-Based Algorithm for Accurate Phylogenetic Breakpoint Mapping and Subtype Prediction in HIV-1

Genetically diverse pathogens (such as Human Immunodeficiency virus type 1, HIV-1) are frequently stratified into phylogenetically or immunologically defined subtypes for classification purposes. Computational identification of such subtypes is helpful in surveillance, epidemiological analysis and detection of novel variants, e.g., circulating recombinant forms in HIV-1. A number of conceptually and technically different techniques have been proposed for determining the subtype of a query sequence, but there is not a universally optimal approach. We present a model-based phylogenetic method for automatically subtyping an HIV-1 (or other viral or bacterial) sequence, mapping the location of breakpoints and assigning parental sequences in recombinant strains as well as computing confidence levels for the inferred quantities. Our Subtype Classification Using Evolutionary ALgorithms (SCUEAL) procedure is shown to perform very well in a variety of simulation scenarios, runs in parallel when multiple sequences are being screened, and matches or exceeds the performance of existing approaches on typical empirical cases. We applied SCUEAL to all available polymerase (pol) sequences from two large databases, the Stanford Drug Resistance database and the UK HIV Drug Resistance Database. Comparing with subtypes which had previously been assigned revealed that a minor but substantial (≈5%) fraction of pure subtype sequences may in fact be within- or inter-subtype recombinants. A free implementation of SCUEAL is provided as a module for the HyPhy package and the Datamonkey web server. Our method is especially useful when an accurate automatic classification of an unknown strain is desired, and is positioned to complement and extend faster but less accurate methods. Given the increasingly frequent use of HIV subtype information in studies focusing on the effect of subtype on treatment, clinical outcome, pathogenicity and vaccine design, the importance of accurate, robust and extensible subtyping procedures is clear.

The Emerging Mutational Landscape of G-proteins and G-protein Coupled Receptors in Cancer

Aberrant expression and activity of G proteins and G-protein-coupled receptors (GPCRs) are frequently associated with tumorigenesis. Deep sequencing studies show that 4.2% of tumours carry activating mutations in GNAS (encoding Gαs), and that oncogenic activating mutations in genes encoding Gαq family members (GNAQ or GNA11) are present in ∼66% and ∼6% of melanomas arising in the eye and skin, respectively. Furthermore, nearly 20% of human tumours harbour mutations in GPCRs. Many human cancer-associated viruses also express constitutively active viral GPCRs. These studies indicate that G proteins, GPCRs and their linked signalling circuitry represent novel therapeutic targets for cancer prevention and treatment.

A comprehensive transcriptional portrait of human cancer cell lines

Tumor-derived cell lines have served as vital models to advance our understanding of oncogene function and therapeutic responses. Although substantial effort has been made to define the genomic constitution of cancer cell line panels, the transcriptome remains understudied. Here we describe RNA sequencing and single-nucleotide polymorphism (SNP) array analysis of 675 human cancer cell lines. We report comprehensive analyses of transcriptome features including gene expression, mutations, gene fusions and expression of non-human sequences. Of the 2,200 gene fusions catalogued, 1,435 consist of genes not previously found in fusions, providing many leads for further investigation. We combine multiple genome and transcriptome features in a pathway-based approach to enhance prediction of response to targeted therapeutics. Our results provide a valuable resource for studies that use cancer cell lines.

Oncogenic ERBB3 Mutations in Human Cancers

The human epidermal growth factor receptor (HER) family of tyrosine kinases is deregulated in multiple cancers either through amplification, overexpression, or mutation. ERBB3/HER3, the only member with an impaired kinase domain, although amplified or overexpressed in some cancers, has not been reported to carry oncogenic mutations. Here, we report the identification of ERBB3 somatic mutations in ~11% of colon and gastric cancers. We found that the ERBB3 mutants transformed colonic and breast epithelial cells in a ligand-independent manner. However, the mutant ERBB3 oncogenic activity was dependent on kinase-active ERBB2. Furthermore, we found that anti-ERBB antibodies and small molecule inhibitors effectively blocked mutant ERBB3-mediated oncogenic signaling and disease progression in vivo.

Coreceptor Tropism in Human Immunodeficiency Virus Type 1 Subtype D: High Prevalence of CXCR4 Tropism and Heterogeneous Composition of Viral Populations

ABSTRACT In human immunodeficiency virus type 1 (HIV-1) subtype B, CXCR4 coreceptor use ranges from ∼20% in early infection to ∼50% in advanced disease. Coreceptor use by non-subtype B HIV is less well characterized. We studied coreceptor tropism of subtype A and D HIV-1 collected from 68 pregnant, antiretroviral drug-naive Ugandan women (HIVNET 012 trial). None of 33 subtype A or 10 A/D-recombinant viruses used the CXCR4 coreceptor. In contrast, nine (36%) of 25 subtype D viruses used both CXCR4 and CCR5 coreceptors. Clonal analyses of the nine subtype D samples with dual or mixed tropism revealed heterogeneous viral populations comprised of X4-, R5-, and dual-tropic HIV-1 variants. In five of the six samples with dual-tropic strains, V3 loop sequences of dual-tropic clones were identical to those of cocirculating R5-tropic clones, indicating the presence of CXCR4 tropism determinants outside of the V3 loop. These dual-tropic variants with R5-tropic-like V3 loops, which we designated “dual-R,” use CCR5 much more efficiently than CXCR4, in contrast to dual-tropic clones with X4-tropic-like V3 loops (“dual-X”). These observations have implications for pathogenesis and treatment of subtype D-infected individuals, for the association between V3 sequence and coreceptor tropism phenotype, and for understanding potential mechanisms of evolution from exclusive CCR5 use to efficient CXCR4 use by subtype D HIV-1.

Comprehensive genomic analysis of malignant pleural mesothelioma identifies recurrent mutations, gene fusions and splicing alterations

We analyzed transcriptomes (n = 211), whole exomes (n = 99) and targeted exomes (n = 103) from 216 malignant pleural mesothelioma (MPM) tumors. Using RNA-seq data, we identified four distinct molecular subtypes: sarcomatoid, epithelioid, biphasic-epithelioid (biphasic-E) and biphasic-sarcomatoid (biphasic-S). Through exome analysis, we found BAP1, NF2, TP53, SETD2, DDX3X, ULK2, RYR2, CFAP45, SETDB1 and DDX51 to be significantly mutated (q-score ≥ 0.8) in MPMs. We identified recurrent mutations in several genes, including SF3B1 (∼2%; 4/216) and TRAF7 (∼2%; 5/216). SF3B1-mutant samples showed a splicing profile distinct from that of wild-type tumors. TRAF7 alterations occurred primarily in the WD40 domain and were, except in one case, mutually exclusive with NF2 alterations. We found recurrent gene fusions and splice alterations to be frequent mechanisms for inactivation of NF2, BAP1 and SETD2. Through integrated analyses, we identified alterations in Hippo, mTOR, histone methylation, RNA helicase and p53 signaling pathways in MPMs.

A systems analysis of mutational effects in HIV-1 protease and reverse transcriptase

The development of a quantitative understanding of viral evolution and the fitness landscape in HIV-1 drug resistance is a formidable challenge given the large number of available drugs and drug resistance mutations. We analyzed a dataset measuring the in vitro fitness of 70,081 virus samples isolated from HIV-1 subtype B infected individuals undergoing routine drug resistance testing. We assayed virus samples for in vitro replicative capacity in the absence of drugs as well as in the presence of 15 individual drugs. We employed a generalized kernel ridge regression to estimate main fitness effects and epistatic interactions of 1,859 single amino acid variants found within the HIV-1 protease and reverse transcriptase sequences. Models including epistatic interactions predict an average of 54.8% of the variance in replicative capacity across the 16 different environments and substantially outperform models based on main fitness effects only. We find that the fitness landscape of HIV-1 protease and reverse transcriptase is characterized by strong epistasis.

Spectrum of diverse genomic alterations define non–clear cell renal carcinoma subtypes

To further understand the molecular distinctions between kidney cancer subtypes, we analyzed exome, transcriptome and copy number alteration data from 167 primary human tumors that included renal oncocytomas and non–clear cell renal cell carcinomas (nccRCCs), consisting of papillary (pRCC), chromophobe (chRCC) and translocation (tRCC) subtypes. We identified ten significantly mutated genes in pRCC, including MET, NF2, SLC5A3, PNKD and CPQ. MET mutations occurred in 15% (10/65) of pRCC samples and included previously unreported recurrent activating mutations. In chRCC, we found TP53, PTEN, FAAH2, PDHB, PDXDC1 and ZNF765 to be significantly mutated. Gene expression analysis identified a five-gene set that enabled the molecular classification of chRCC, renal oncocytoma and pRCC. Using RNA sequencing, we identified previously unreported gene fusions, including ACTG1-MITF fusion. Ectopic expression of the ACTG1-MITF fusion led to cellular transformation and induced the expression of downstream target genes. Finally, we observed upregulation of the anti-apoptotic factor BIRC7 in MiTF-high RCC tumors, suggesting a potential therapeutic role for BIRC7 inhibitors.