pgMAP: a pipeline to enable guide RNA read mapping from dual-targeting CRISPR screens.

We developed pgMAP, an analysis pipeline to map gRNA sequencing reads from dual-targeting CRISPR screens. pgMAP output includes a dual gRNA read counts table and quality control metrics including the proportion of correctly-paired reads and CRISPR library sequencing coverage across all time points and samples. pgMAP is implemented using Snakemake and is available open-source under the MIT license at https://github.com/fredhutch/pgmap_pipeline.

Discovery of synthetic lethal and tumor suppressor paralog pairs in the human genome.

CRISPR screens have accelerated the discovery of important cancer vulnerabilities. However, single-gene knockout phenotypes can be masked by redundancy among related genes. Paralogs constitute two-thirds of the human protein-coding genome, so existing methods are likely inadequate for assaying a large portion of gene function. Here, we develop paired guide RNAs for paralog genetic interaction mapping (pgPEN), a pooled CRISPR-Cas9 single- and double-knockout approach targeting more than 2,000 human paralogs. We apply pgPEN to two cell types and discover that 12% of human paralogs exhibit synthetic lethality in at least one context. We recover known synthetic lethal paralogs MEK1/MEK2, important drug targets CDK4/CDK6, and other synthetic lethal pairs including CCNL1/CCNL2. Additionally, we identify ten tumor suppressor paralog pairs whose compound loss promotes cell proliferation. These findings nominate drug targets and suggest that paralog genetic interactions could shape the landscape of positive and negative selection in cancer.

Integrative oncogene-dependency mapping identifies RIT1 vulnerabilities and synergies in lung cancer.

CRISPR-based cancer dependency maps are accelerating advances in cancer precision medicine, but adequate functional maps are limited to the most common oncogenes. To identify opportunities for therapeutic intervention in other rarer subsets of cancer, we investigate the oncogene-specific dependencies conferred by the lung cancer oncogene, RIT1. Here, genome-wide CRISPR screening in KRAS, EGFR, and RIT1-mutant isogenic lung cancer cells identifies shared and unique vulnerabilities of each oncogene. Combining this genetic data with small-molecule sensitivity profiling, we identify a unique vulnerability of RIT1-mutant cells to loss of spindle assembly checkpoint regulators. Oncogenic RIT1M90I weakens the spindle assembly checkpoint and perturbs mitotic timing, resulting in sensitivity to Aurora A inhibition. In addition, we observe synergy between mutant RIT1 and activation of YAP1 in multiple models and frequent nuclear overexpression of YAP1 in human primary RIT1-mutant lung tumors. These results provide a genome-wide atlas of oncogenic RIT1 functional interactions and identify components of the RAS pathway, spindle assembly checkpoint, and Hippo/YAP1 network as candidate therapeutic targets in RIT1-mutant lung cancer.

CRISPR base editor screens identify variant function at scale.

Cuella-Martin et al. (2021) and Hanna et al. (2021) showcase CRISPR base editing in large-scale pooled screens in human cells to discover both loss- and gain-of-function variants, enabling protein structure/function insights and clinical variant interpretation.

Whole exome sequencing in patients with Williams-Beuren syndrome followed by disease modeling in mice points to four novel pathways that may modify stenosis risk.

Supravalvular aortic stenosis (SVAS) is a narrowing of the aorta caused by elastin (ELN) haploinsufficiency. SVAS severity varies among patients with Williams-Beuren syndrome (WBS), a rare disorder that removes one copy of ELN and 25-27 other genes. Twenty percent of children with WBS require one or more invasive and often risky procedures to correct the defect while 30% have no appreciable stenosis, despite sharing the same basic genetic lesion. There is no known medical therapy. Consequently, identifying genes that modify SVAS offers the potential for novel modifier-based therapeutics. To improve statistical power in our rare-disease cohort (N = 104 exomes), we utilized extreme-phenotype cohorting, functional variant filtration and pathway-based analysis. Gene set enrichment analysis of exome-wide association data identified increased adaptive immune system variant burden among genes associated with SVAS severity. Additional enrichment, using only potentially pathogenic variants known to differ in frequency between the extreme phenotype subsets, identified significant association of SVAS severity with not only immune pathway genes, but also genes involved with the extracellular matrix, G protein-coupled receptor signaling and lipid metabolism using both SKAT-O and RQTest. Complementary studies in Eln+/-; Rag1-/- mice, which lack a functional adaptive immune system, showed improvement in cardiovascular features of ELN insufficiency. Similarly, studies in mixed background Eln+/- mice confirmed that variations in genes that increase elastic fiber deposition also had positive impact on aortic caliber. By using tools to improve statistical power in combination with orthogonal analyses in mice, we detected four main pathways that contribute to SVAS risk.

Social, neurodevelopmental, endocrine, and head size differences associated with atypical deletions in Williams-Beuren syndrome.

Williams-Beuren syndrome (WBS) is a multisystem disorder caused by a hemizygous deletion on 7q11.23 encompassing 26-28 genes. An estimated 2-5% of patients have "atypical" deletions, which extend in the centromeric and/or telomeric direction from the WBS critical region. To elucidate clinical differentiators among these deletion types, we evaluated 10 individuals with atypical deletions in our cohort and 17 individuals with similarly classified deletions previously described in the literature. Larger deletions in either direction often led to more severe developmental delays, while deletions containing MAGI2 were associated with infantile spasms and seizures in patients. In addition, head size was notably smaller in those with centromeric deletions including AUTS2. Because children with atypical deletions were noted to be less socially engaged, we additionally sought to determine how atypical deletions relate to social phenotypes. Using the Social Responsiveness Scale-2, raters scored individuals with atypical deletions as having different social characteristics to those with typical WBS deletions (p = .001), with higher (more impaired) scores for social motivation (p = .005) in the atypical deletion group. In recognizing these distinctions, physicians can better identify patients, including those who may already carry a clinical or FISH WBS diagnosis, who may benefit from additional molecular evaluation, screening, and therapy. In addition to the clinical findings, we note mild endocrine findings distinct from those typically seen in WBS in several patients with telomeric deletions that included POR. Further study in additional telomeric deletion cases will be needed to confirm this observation.

Exome sequencing of 85 Williams-Beuren syndrome cases rules out coding variation as a major contributor to remaining variance in social behavior.

BACKGROUND: Large, multigenic deletions at chromosome 7q11.23 result in a highly penetrant constellation of physical and behavioral symptoms known as Williams-Beuren syndrome (WS). Of particular interest is the unusual social-cognitive profile evidenced by deficits in social cognition and communication reminiscent of autism spectrum disorders (ASD) that are juxtaposed with normal or even relatively enhanced social motivation. Interestingly, duplications in the same region also result in ASD-like phenotypes as well as social phobias. Thus, the region clearly regulates human social motivation and behavior, yet the relevant gene(s) have not been definitively identified. METHOD: Here, we deeply phenotyped 85 individuals with WS and used exome sequencing to analyze common and rare variation for association with the remaining variance in social behavior as assessed by the Social Responsiveness Scale. RESULTS: We replicated the previously reported unusual juxtaposition of behavioral symptoms in this new patient collection, but we did not find any new alleles of large effect in the targeted analysis of the remaining copy of genes in the Williams syndrome critical region. However, we report on two nominally significant SNPs in two genes that have been implicated in the cognitive and social phenotypes of Williams syndrome, BAZ1B and GTF2IRD1. Secondary discovery driven explorations focusing on known ASD genes and an exome wide scan do not highlight any variants of a large effect. CONCLUSIONS: Whole exome sequencing of 85 individuals with WS did not support the hypothesis that there are variants of large effect within the remaining Williams syndrome critical region that contribute to the social phenotype. This deeply phenotyped and genotyped patient cohort with a defined mutation provides the opportunity for similar analyses focusing on noncoding variation and/or other phenotypic domains.

Programmed evolution for optimization of orthogonal metabolic output in bacteria.

Current use of microbes for metabolic engineering suffers from loss of metabolic output due to natural selection. Rather than combat the evolution of bacterial populations, we chose to embrace what makes biological engineering unique among engineering fields - evolving materials. We harnessed bacteria to compute solutions to the biological problem of metabolic pathway optimization. Our approach is called Programmed Evolution to capture two concepts. First, a population of cells is programmed with DNA code to enable it to compute solutions to a chosen optimization problem. As analog computers, bacteria process known and unknown inputs and direct the output of their biochemical hardware. Second, the system employs the evolution of bacteria toward an optimal metabolic solution by imposing fitness defined by metabolic output. The current study is a proof-of-concept for Programmed Evolution applied to the optimization of a metabolic pathway for the conversion of caffeine to theophylline in E. coli. Introduced genotype variations included strength of the promoter and ribosome binding site, plasmid copy number, and chaperone proteins. We constructed 24 strains using all combinations of the genetic variables. We used a theophylline riboswitch and a tetracycline resistance gene to link theophylline production to fitness. After subjecting the mixed population to selection, we measured a change in the distribution of genotypes in the population and an increased conversion of caffeine to theophylline among the most fit strains, demonstrating Programmed Evolution. Programmed Evolution inverts the standard paradigm in metabolic engineering by harnessing evolution instead of fighting it. Our modular system enables researchers to program bacteria and use evolution to determine the combination of genetic control elements that optimizes catabolic or anabolic output and to maintain it in a population of cells. Programmed Evolution could be used for applications in energy, pharmaceuticals, chemical commodities, biomining, and bioremediation.