CRISPR off-target analysis in genetically engineered rats and mice.

Despite widespread use of CRISPR, comprehensive data on the frequency and impact of Cas9-mediated off-targets in modified rodents are limited. Here we present deep-sequencing data from 81 genome-editing projects on mouse and rat genomes at 1,423 predicted off-target sites, 32 of which were confirmed, and show that high-fidelity Cas9 versions reduced off-target mutation rates in vivo. Using whole-genome sequencing data from ten mouse embryos, treated with a single guide RNA (sgRNA), and from their genetic parents, we found 43 off-targets, 30 of which were predicted by an adapted version of GUIDE-seq.

Mutational landscape of antibody variable domains reveals a switch modulating the interdomain conformational dynamics and antigen binding.

Somatic mutations within the antibody variable domains are critical to the immense capacity of the immune repertoire. Here, via a deep mutational scan, we dissect how mutations at all positions of the variable domains of a high-affinity anti-VEGF antibody G6.31 impact its antigen-binding function. The resulting mutational landscape demonstrates that large portions of antibody variable domain positions are open to mutation, and that beneficial mutations can be found throughout the variable domains. We determine the role of one antigen-distal light chain position 83, demonstrating that mutation at this site optimizes both antigen affinity and thermostability by modulating the interdomain conformational dynamics of the antigen-binding fragment. Furthermore, by analyzing a large number of human antibody sequences and structures, we demonstrate that somatic mutations occur frequently at position 83, with corresponding domain conformations observed for G6.31. Therefore, the modulation of interdomain dynamics represents an important mechanism during antibody maturation in vivo.

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.

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.

Disruption of PH-kinase domain interactions leads to oncogenic activation of AKT in human cancers.

The protein kinase v-akt murine thymoma viral oncogene homolog (AKT), a key regulator of cell survival and proliferation, is frequently hyperactivated in human cancers. Intramolecular pleckstrin homology (PH) domain-kinase domain (KD) interactions are important in maintaining AKT in an inactive state. AKT activation proceeds after a conformational change that dislodges the PH from the KD. To understand these autoinhibitory interactions, we generated mutations at the PH-KD interface and found that most of them lead to constitutive activation of AKT. Such mutations are likely another mechanism by which activation may occur in human cancers and other diseases. In support of this likelihood, we found somatic mutations in AKT1 at the PH-KD interface that have not been previously described in human cancers. Furthermore, we show that the AKT1 somatic mutants are constitutively active, leading to oncogenic signaling. Additionally, our studies show that the AKT1 mutants are not effectively inhibited by allosteric AKT inhibitors, consistent with the requirement for an intact PH-KD interface for allosteric inhibition. These results have important implications for therapeutic intervention in patients with AKT mutations at the PH-KD interface.

Regulation of Class IA PI 3-kinases: C2 domain-iSH2 domain contacts inhibit p85/p110alpha and are disrupted in oncogenic p85 mutants.

We previously proposed a model of Class IA PI3K regulation in which p85 inhibition of p110alpha requires (i) an inhibitory contact between the p85 nSH2 domain and the p110alpha helical domain, and (ii) a contact between the p85 nSH2 and iSH2 domains that orients the nSH2 so as to inhibit p110alpha. We proposed that oncogenic truncations of p85 fail to inhibit p110 due to a loss of the iSH2-nSH2 contact. However, we now find that within the context of a minimal regulatory fragment of p85 (the nSH2-iSH2 fragment, termed p85ni), the nSH2 domain rotates much more freely (tau(c) approximately 12.7 ns) than it could if it were interacting rigidly with the iSH2 domain. These data are not compatible with our previous model. We therefore tested an alternative model in which oncogenic p85 truncations destabilize an interface between the p110alpha C2 domain (residue N345) and the p85 iSH2 domain (residues D560 and N564). p85ni-D560K/N564K shows reduced inhibition of p110alpha, similar to the truncated p85ni-572(STOP). Conversely, wild-type p85ni poorly inhibits p110alphaN345K. Strikingly, the p110alphaN345K mutant is inhibited to the same extent by the wild-type or truncated p85ni, suggesting that mutation of p110alpha-N345 is not additive with the p85ni-572(STOP) mutation. Similarly, the D560K/N564K mutation is not additive with the p85ni-572(STOP) mutant for downstream signaling or cellular transformation. Thus, our data suggests that mutations at the C2-iSH2 domain contact and truncations of the iSH2 domain, which are found in human tumors, both act by disrupting the C2-iSH2 domain interface.

CRISPR whole-genome screening identifies new necroptosis regulators and RIPK1 alternative splicing.

The necroptotic cell death pathway is a key component of human pathogen defense that can become aberrantly derepressed during tissue homeostasis to contribute to multiple types of tissue damage and disease. While formation of the necrosome kinase signaling complex containing RIPK1, RIPK3, and MLKL has been extensively characterized, additional mechanisms of its regulation and effector functions likely remain to be discovered. We screened 19,883 mouse protein-coding genes by CRISPR/Cas9-mediated gene knockout for resistance to cytokine-induced necroptosis and identified 112 regulators and mediators of necroptosis, including 59 new candidate pathway components with minimal or no effect on cell growth in the absence of necroptosis induction. Among these, we further characterized the function of PTBP1, an RNA binding protein whose activity is required to maintain RIPK1 protein abundance by regulating alternative splice-site selection.

Author Correction: Comprehensive analysis of single molecule sequencing-derived complete genome and whole transcriptome of Hyposidra talaca nuclear polyhedrosis virus.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has been fixed in the paper.

Comprehensive analysis of single molecule sequencing-derived complete genome and whole transcriptome of Hyposidra talaca nuclear polyhedrosis virus.

We sequenced the Hyposidra talaca NPV (HytaNPV) double stranded circular DNA genome using PacBio single molecule sequencing technology. We found that the HytaNPV genome is 139,089 bp long with a GC content of 39.6%. It encodes 141 open reading frames (ORFs) including the 37 baculovirus core genes, 25 genes conserved among lepidopteran baculoviruses, 72 genes known in baculovirus, and 7 genes unique to the HytaNPV genome. It is a group II alphabaculovirus that codes for the F protein and lacks the gp64 gene found in group I alphabaculovirus viruses. Using RNA-seq, we confirmed the expression of the ORFs identified in the HytaNPV genome. Phylogenetic analysis showed HytaNPV to be closest to BusuNPV, SujuNPV and EcobNPV that infect other tea pests, Buzura suppressaria, Sucra jujuba, and Ectropis oblique, respectively. We identified repeat elements and a conserved non-coding baculovirus element in the genome. Analysis of the putative promoter sequences identified motif consistent with the temporal expression of the genes observed in the RNA-seq data.

Silencing of retrotransposons by SETDB1 inhibits the interferon response in acute myeloid leukemia.

A propensity for rewiring genetic and epigenetic regulatory networks, thus enabling sustained cell proliferation, suppression of apoptosis, and the ability to evade the immune system, is vital to cancer cell propagation. An increased understanding of how this is achieved is critical for identifying or improving therapeutic interventions. In this study, using acute myeloid leukemia (AML) human cell lines and a custom CRISPR/Cas9 screening platform, we identify the H3K9 methyltransferase SETDB1 as a novel, negative regulator of innate immunity. SETDB1 is overexpressed in many cancers, and loss of this gene in AML cells triggers desilencing of retrotransposable elements that leads to the production of double-stranded RNAs (dsRNAs). This is coincident with induction of a type I interferon response and apoptosis through the dsRNA-sensing pathway. Collectively, our findings establish a unique gene regulatory axis that cancer cells can exploit to circumvent the immune system.

A rapid method for estimation of baculovirus titer based on viable cell size.

Baculovirus protein expression system is a powerful tool for producing recombinant proteins. To optimize conditions for efficient recombinant protein expression, it is important to determine titer of virus stock for arriving at an optimal multiplicity of infection (MOI) that maximizes recombinant protein expression. Traditionally plaque assays have been used for titer determination. Other methods such as endpoint dilution, quantitative real-time polymerase chain reaction and flow cytometry have been developed to aid the determination of virus titers. However, most of these methods are time-consuming and labor intensive. In this regard, a simple and rapid method for determination of virus titers based on the cytopathic effects that lead to viable cell size increase following virus infection is presented in this paper. In this study, the Vi-CELL (Beckman Coulter) was used to measure cell-diameter over a range of virus dilutions, following infection. Applying statistical modeling techniques, the viable cell-diameter data was used to estimate the virus titer. The results indicated that the viable cell-diameter based titer estimation to be reliable and comparable to titers determined by the traditional plaque assay.

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.

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.

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.

Combined targeting of BRAF and CRAF or BRAF and PI3K effector pathways is required for efficacy in NRAS mutant tumors.

BACKGROUND: Oncogenic RAS is a highly validated cancer target. Attempts at targeting RAS directly have so far not succeeded in the clinic. Understanding downstream RAS-effectors that mediate oncogenesis in a RAS mutant setting will help tailor treatments that use RAS-effector inhibitors either alone or in combination to target RAS-driven tumors. METHODOLOGY/PRINCIPAL FINDINGS: In this study, we have investigated the sufficiency of targeting RAS-effectors, RAF, MEK and PI3-Kinase either alone or in combination in RAS mutant lines, using an inducible shRNA in vivo mouse model system. We find that in colon cancer cells harboring a KRAS(G13D) mutant allele, knocking down KRAS alone or the RAFs in combination or the RAF effectors, MEK1 and MEK2, together is effective in delaying tumor growth in vivo. In melanoma cells harboring an NRAS(Q61L) or NRAS(Q61K) mutant allele, we find that targeting NRAS alone or both BRAF and CRAF in combination or both BRAF and PIK3CA together showed efficacy. CONCLUSION/SIGNIFICANCE: Our data indicates that targeting oncogenic NRAS-driven melanomas require decrease in both pERK and pAKT downstream of RAS-effectors for efficacy. This can be achieved by either targeting both BRAF and CRAF or BRAF and PIK3CA simultaneously in NRAS mutant tumor cells.

Somatic mutations in p85alpha promote tumorigenesis through class IA PI3K activation.

Members of the mammalian phosphoinositide-3-OH kinase (PI3K) family of proteins are critical regulators of various cellular process including cell survival, growth, proliferation, and motility. Oncogenic activating mutations in the p110alpha catalytic subunit of the heterodimeric p110/p85 PI3K enzyme are frequent in human cancers. Here we show the presence of frequent mutations in p85alpha in colon cancer, a majority of which occurs in the inter-Src homology-2 (iSH2) domain. These mutations uncouple and retain p85alpha's p110-stabilizing activity, while abrogating its p110-inhibitory activity. The p85alpha mutants promote cell survival, AKT activation, anchorage-independent cell growth, and oncogenesis in a p110-dependent manner.

Estimation of baculovirus titer based on viable cell size.

In this paper, a simple and rapid protocol for determination of baculovirus titers based on increasing viable insect cell size/diameter following virus infection is presented. There are different methods available for determining virus titers such as plaque assays end-point dilution, quantitative real-time polymerase chain reaction and flow cytometry. However, most of these methods are time consuming and labor intensive. The titer estimation method presented here can be completed in approximately 28 h from start to finish. In this method, the Vi-CELL (Beckman Coulter) was used to measure cell diameter change over a range of virus dilutions, following infection. The cell diameter change data were used to compute the virus titer using a statistical method called the method of moments that we have described previously.