Metabolic balance in colorectal cancer is maintained by optimal Wnt signaling levels

Wnt pathways are important for the modulation of tissue homeostasis, and their deregulation is linked to cancer development. Canonical Wnt signaling is hyperactivated in many human colorectal cancers due to genetic alterations of the negative Wnt regulator APC. However, the expression levels of Wnt‐dependent targets vary between tumors, and the mechanisms of carcinogenesis concomitant with this Wnt signaling dosage have not been understood. In this study, we integrate whole‐genome CRISPR/Cas9 screens with large‐scale multi‐omic data to delineate functional subtypes of cancer. We engineer APC loss‐of‐function mutations and thereby hyperactivate Wnt signaling in cells with low endogenous Wnt activity and find that the resulting engineered cells have an unfavorable metabolic equilibrium compared with cells which naturally acquired Wnt hyperactivation. We show that the dosage level of oncogenic Wnt hyperactivation impacts the metabolic equilibrium and the mitochondrial phenotype of a given cell type in a context‐dependent manner. These findings illustrate the impact of context‐dependent genetic interactions on cellular phenotypes of a central cancer driver mutation and expand our understanding of quantitative modulation of oncogenic signaling in tumorigenesis.     

Multiplexed CRISPR/CAS9‐mediated engineering of pre‐clinical mouse models bearing native human B cell receptors

B‐cell receptor (BCR) knock‐in (KI) mouse models play an important role in vaccine development and fundamental immunological studies. However, the time required to generate them poses a bottleneck. Here we report a one‐step CRISPR/Cas9 KI methodology to combine the insertion of human germline immunoglobulin heavy and light chains at their endogenous loci in mice. We validate this technology with the rapid generation of three BCR KI lines expressing native human precursors, instead of computationally inferred germline sequences, to HIV broadly neutralizing antibodies. We demonstrate that B cells from these mice are fully functional: upon transfer to congenic, wild type mice at controlled frequencies, such B cells can be primed by eOD‐GT8 60mer, a germline‐targeting immunogen currently in clinical trials, recruited to germinal centers, secrete class‐switched antibodies, undergo somatic hypermutation, and differentiate into memory B cells. KI mice expressing functional human BCRs promise to accelerate the development of vaccines for HIV and other infectious diseases.     

Minimized combinatorial CRISPR screens identify genetic interactions in autophagy

Combinatorial CRISPR-Cas screens have advanced the mapping of genetic interactions, but their experimental scale limits the number of targetable gene combinations. Here, we describe 3Cs multiplexing, a rapid and scalable method to generate highly diverse and uniformly distributed combinatorial CRISPR libraries. We demonstrate that the library distribution skew is the critical determinant of its required screening coverage. By circumventing iterative cloning of PCR-amplified oligonucleotides, 3Cs multiplexing facilitates the generation of combinatorial CRISPR libraries with low distribution skews. We show that combinatorial 3Cs libraries can be screened with minimal coverages, reducing associated efforts and costs at least 10-fold. We apply a 3Cs multiplexing library targeting 12,736 autophagy gene combinations with 247,032 paired gRNAs in viability and reporter-based enrichment screens. In the viability screen, we identify, among others, the synthetic lethal WDR45B-PIK3R4 and the proliferation-enhancing ATG7-KEAP1 genetic interactions. In the reporter-based screen, we identify over 1,570 essential genetic interactions for autophagy flux, including interactions among paralogous genes, namely ATG2A-ATG2B, GABARAP-MAP1LC3B and GABARAP-GABARAPL2. However, we only observe few genetic interactions within paralogous gene families of more than two members, indicating functional compensation between them. This work establishes 3Cs multiplexing as a platform for genetic interaction screens at scale.     

Behavioral fingerprints predict insecticide and anthelmintic mode of action

Novel invertebrate‐killing compounds are required in agriculture and medicine to overcome resistance to existing treatments. Because insecticides and anthelmintics are discovered in phenotypic screens, a crucial step in the discovery process is determining the mode of action of hits. Visible whole‐organism symptoms are combined with molecular and physiological data to determine mode of action. However, manual symptomology is laborious and requires symptoms that are strong enough to see by eye. Here, we use high‐throughput imaging and quantitative phenotyping to measure Caenorhabditis elegans behavioral responses to compounds and train a classifier that predicts mode of action with an accuracy of 88% for a set of ten common modes of action. We also classify compounds within each mode of action to discover substructure that is not captured in broad mode‐of‐action labels. High‐throughput imaging and automated phenotyping could therefore accelerate mode‐of‐action discovery in invertebrate‐targeting compound development and help to refine mode‐of‐action categories.     

Sympatry leads to reduced body condition in chickadees that occasionally hybridize

Both abiotic and biotic drivers influence species distributions. Abiotic drivers such as climate have received considerable attention, even though biotic drivers such as hybridization often interact with abiotic drivers. We sought to explore the (1) costs of co‐occurrence for ecologically similar species that hybridize and (2) associations between ecological factors and condition to understand how abiotic and biotic factors influence species distributions. For two closely related and ecologically similar songbirds, black‐capped and mountain chickadees, we characterized body condition, as a proxy for fitness, using a 1358‐individual range‐wide dataset. We compared body condition in sympatry and allopatry with several abiotic and biotic factors using species‐specific generalized linear mixed models. We generated genomic data for a subset of 217 individuals to determine the extent of hybridization‐driven admixture in our dataset. Within this data subset, we found that ~11% of the chickadees had hybrid ancestry, and all hybrid individuals had typical black‐capped chickadee plumage. In the full dataset, we found that birds of both species, independent of demographic and abiotic factors, had significantly lower body condition when occurring in sympatry than birds in allopatry. This could be driven by either the inclusion of cryptic, likely poor condition, hybrids in our full dataset, competitive interactions in sympatry, or range edge effects. We are currently unable to discriminate between these mechanisms. Our findings have implications for mountain chickadees in particular, which will encounter more black‐capped chickadees as black‐capped chickadee ranges shift upslope and could lead to local declines in mountain chickadee populations.     

CRISPR/Cas12a-mediated CHO genome engineering can be effectively integrated at multiple stages of the cell line generation process for bioproduction

Most biopharmaceuticals produced today are generated using Chinese hamster ovary (CHO) cells, therefore significant attention is focused on methods to improve CHO cell productivity and product quality. The discovery of gene-editing tools, such as CRISPR/Cas9, offers new opportunities to improve CHO cell bioproduction through cell line engineering. Recently an additional CRISPR-associated protein, Cas12a (Cpf1), was shown to be effective for gene editing in eukaryotic cells, including CHO. In this study, we demonstrate the successful application of CRISPR/Cas12a for the generation of clonally derived CHO knockout (KO) cell lines with improved product quality attributes. While we found Cas12a efficiency to be highly dependent on the targeting RNA used, we were able to generate CHO KO cell lines using small screens of only 96–320 clonally derived cell lines. Additionally, we present a novel bulk culture analysis approach that can be used to quickly assess CRISPR RNA efficiency and determine ideal screen sizes for generating genetic KO cell lines. Most critically, we find that Cas12a can be directly integrated into the cell line generation process through cotransfection with no negative impact on titer or screen size. Overall, our results show CRISPR/Cas12a to be an efficient and effective CHO genome editing tool.     

Phosphoproteomics reveals novel modes of function and inter‐relationships among PIKKs in response to genotoxic stress

The DNA damage response (DDR) is a complex signaling network that relies on cascades of protein phosphorylation, which are initiated by three protein kinases of the family of PI3‐kinase‐related protein kinases (PIKKs): ATM, ATR, and DNA‐PK. ATM is missing or inactivated in the genome instability syndrome, ataxia‐telangiectasia (A‐T). The relative shares of these PIKKs in the response to genotoxic stress and the functional relationships among them are central questions in the genome stability field. We conducted a comprehensive phosphoproteomic analysis in human wild‐type and A‐T cells treated with the double‐strand break‐inducing chemical, neocarzinostatin, and validated the results with the targeted proteomic technique, selected reaction monitoring. We also matched our results with 34 published screens for DDR factors, creating a valuable resource for identifying strong candidates for novel DDR players. We uncovered fine‐tuned dynamics between the PIKKs following genotoxic stress, such as DNA‐PK‐dependent attenuation of ATM. In A‐T cells, partial compensation for ATM absence was provided by ATR and DNA‐PK, with distinct roles and kinetics. The results highlight intricate relationships between these PIKKs in the DDR.     

Immunomodulatory activity of extracellular Hsp70 mediated via paired receptors Siglec‐5 and Siglec‐14

The intracellular chaperone heat‐shock protein 70 (Hsp70) can be secreted from cells, but its extracellular role is unclear, as the protein has been reported to both activate and suppress the innate immune response. Potential immunomodulatory receptors on myelomonocytic lineage cells that bind extracellular Hsp70 are not well defined. Siglecs are Ig‐superfamily lectins on mammalian leukocytes that recognize sialic acid‐bearing glycans and thereby modulate immune responses. Siglec‐5 and Siglec‐14, expressed on monocytes and neutrophils, share identical ligand‐binding domains but have opposing signaling functions. Based on phylogenetic analyses of these receptors, we predicted that endogenous sialic acid‐independent ligands should exist. An unbiased screen revealed Hsp70 as a ligand for Siglec‐5 and Siglec‐14. Hsp70 stimulation through Siglec‐5 delivers an anti‐inflammatory signal, while stimulation through Siglec‐14 is pro‐inflammatory. The functional consequences of this interaction are also addressed in relation to a SIGLEC14 polymorphism found in humans. Our results demonstrate that an endogenous non‐sialic acid‐bearing molecule can be either a danger‐associated or self‐associated signal through paired Siglecs, and may explain seemingly contradictory prior reports on extracellular Hsp70 action.     

Receptor‐mediated clustering of FIP200 bypasses the role of LC3 lipidation in autophagy

Autophagosome formation requires multiple autophagy‐related (ATG) factors. However, we find that a subset of autophagy substrates remains robustly targeted to the lysosome in the absence of several core ATGs, including the LC3 lipidation machinery. To address this unexpected result, we performed genome‐wide CRISPR screens identifying genes required for NBR1 flux in ATG7^KO cells. We find that ATG7‐independent autophagy still requires canonical ATG factors including FIP200. However, in the absence of LC3 lipidation, additional factors are required including TAX1BP1 and TBK1. TAX1BP1''s ability to cluster FIP200 around NBR1 cargo and induce local autophagosome formation enforces cargo specificity and replaces the requirement for lipidated LC3. In support of this model, we define a ubiquitin‐independent mode of TAX1BP1 recruitment to NBR1 puncta, highlighting that TAX1BP1 recruitment and clustering, rather than ubiquitin binding per se, is critical for function. Collectively, our data provide a mechanistic basis for reports of selective autophagy in cells lacking the lipidation machinery, wherein receptor‐mediated clustering of upstream autophagy factors drives continued autophagosome formation.     

Toward an integrated map of genetic interactions in cancer cells

Cancer genomes often harbor hundreds of molecular aberrations. Such genetic variants can be drivers or passengers of tumorigenesis and create vulnerabilities for potential therapeutic exploitation. To identify genotype‐dependent vulnerabilities, forward genetic screens in different genetic backgrounds have been conducted. We devised MINGLE, a computational framework to integrate CRISPR/Cas9 screens originating from different libraries building on approaches pioneered for genetic network discovery in model organisms. We applied this method to integrate and analyze data from 85 CRISPR/Cas9 screens in human cancer cells combining functional data with information on genetic variants to explore more than 2.1 million gene‐background relationships. In addition to known dependencies, we identified new genotype‐specific vulnerabilities of cancer cells. Experimental validation of predicted vulnerabilities identified GANAB and PRKCSH as new positive regulators of Wnt/β‐catenin signaling. By clustering genes with similar genetic interaction profiles, we drew the largest genetic network in cancer cells to date. Our scalable approach highlights how diverse genetic screens can be integrated to systematically build informative maps of genetic interactions in cancer, which can grow dynamically as more data are included.     

Large‐scale survey and database of high affinity ligands for peptide recognition modules

Many proteins involved in signal transduction contain peptide recognition modules (PRMs) that recognize short linear motifs (SLiMs) within their interaction partners. Here, we used large‐scale peptide‐phage display methods to derive optimal ligands for 163 unique PRMs representing 79 distinct structural families. We combined the new data with previous data that we collected for the large SH3, PDZ, and WW domain families to assemble a database containing 7,984 unique peptide ligands for 500 PRMs representing 82 structural families. For 74 PRMs, we acquired enough new data to map the specificity profiles in detail and derived position weight matrices and binding specificity logos based on multiple peptide ligands. These analyses showed that optimal peptide ligands resembled peptides observed in existing structures of PRM‐ligand complexes, indicating that a large majority of the phage‐derived peptides are likely to target natural peptide‐binding sites and could thus act as inhibitors of natural protein–protein interactions. The complete dataset has been assembled in an online database (http://www.prm‐db.org) that will enable many structural, functional, and biological studies of PRMs and SLiMs.     

How low can you go? Introducing SeXY: sex identification from low‐quantity sequencing data despite lacking assembled sex chromosomes

Accurate sex identification is crucial for elucidating the biology of a species. In the absence of directly observable sexual characteristics, sex identification of wild fauna can be challenging, if not impossible. Molecular sexing offers a powerful alternative to morphological sexing approaches. Here, we present SeXY, a novel sex‐identification pipeline, for very low‐coverage shotgun sequencing data from a single individual. SeXY was designed to utilize low‐effort screening data for sex identification and does not require a conspecific sex‐chromosome assembly as reference. We assess the accuracy of our pipeline to data quantity by downsampling sequencing data from 100,000 to 1000 mapped reads and to reference genome selection by mapping to a variety of reference genomes of various qualities and phylogenetic distance. We show that our method is 100% accurate when mapping to a high‐quality (highly contiguous N50 > 30 Mb) conspecific genome, even down to 1000 mapped reads. For lower‐quality reference assemblies (N50 < 30 Mb), our method is 100% accurate with 50,000 mapped reads, regardless of reference assembly quality or phylogenetic distance. The SeXY pipeline provides several advantages over previously implemented methods; SeXY (i) requires sequencing data from only a single individual, (ii) does not require assembled conspecific sex chromosomes, or even a conspecific reference assembly, (iii) takes into account variation in coverage across the genome, and (iv) is accurate with only 1000 mapped reads in many cases.     

High-content imaging-based pooled CRISPR screens in mammalian cells

CRISPR (clustered regularly interspaced short palindromic repeats)-based gene inactivation provides a powerful means for linking genes to particular cellular phenotypes. CRISPR-based screening typically uses large genomic pools of single guide RNAs (sgRNAs). However, this approach is limited to phenotypes that can be enriched by chemical selection or FACS sorting. Here, we developed a microscopy-based approach, which we name optical enrichment, to select cells displaying a particular CRISPR-induced phenotype by automated imaging-based computation, mark them by photoactivation of an expressed photoactivatable fluorescent protein, and then isolate the fluorescent cells using fluorescence-activated cell sorting (FACS). A plugin was developed for the open source software μManager to automate the phenotypic identification and photoactivation of cells, allowing ∼1.5 million individual cells to be screened in 8 h. We used this approach to screen 6,092 sgRNAs targeting 544 genes for their effects on nuclear size regulation and identified 14 bona fide hits. These results present a scalable approach to facilitate imaging-based pooled CRISPR screens.