Accessing the Phenotype Gap: Enabling Systematic Investigation of Paralog Functional Complexity with CRISPR

Evaluating sequence and structural similarity metrics for predicting shared paralog functions

Gene duplication is the primary source of new genes, resulting in most genes having identifiable paralogs. Over time, paralog pairs may diverge in some respects but many retain the ability to perform the same functional role. Protein sequence identity is often used as a proxy for functional similarity and can predict shared functions between paralogs as revealed by synthetic lethal experiments. However, the advent of alternative protein representations, including embeddings from protein language models (PLMs) and predicted structures from AlphaFold, raises the possibility that alternative similarity metrics could better capture functional similarity between paralogs. Here, using two species (budding yeast and human) and two different definitions of shared functionality (shared protein-protein interactions and synthetic lethality), we evaluated a variety of alternative similarity metrics. For some tasks, predicted structural similarity or PLM similarity outperform sequence identity, but more importantly these similarity metrics are not redundant with sequence identity, i.e. combining them with sequence identity leads to improved predictions of shared functionality. By adding contextual features, representing similarity to homologous proteins within and across species, we can significantly enhance our predictions of shared paralog functionality. Overall, our results suggest that alternative similarity metrics capture complementary aspects of functional similarity beyond sequence identity alone.

Incomplete paralog compensation generates selective dependency on TRA2A in cancer

Paralogs often exhibit functional redundancy, allowing them to effectively compensate for each other's loss. However, this buffering mechanism is frequently disrupted in cancer, exposing unique paralog-specific vulnerabilities. Here, we identify a selective dependency on the splicing factor TRA2A. We find that TRA2A and its paralog TRA2B are synthetic lethal partners that function as widespread and largely redundant activators of both alternative and constitutive splicing. While loss of TRA2A alone is typically neutral due to compensation by TRA2B, we discover that a subset of cancer cell lines are highly TRA2A-dependent. Upon TRA2A depletion, these cell lines exhibit a lack of paralog buffering specifically on shared splicing targets, leading to defects in mitosis and cell death. Notably, TRA2B overexpression rescues both the aberrant splicing and lethality associated with TRA2A loss, indicating that paralog compensation is dosage-sensitive. Together, these findings reveal a complex dosage-dependent relationship between paralogous splicing factors, and highlight how dysfunctional paralog buffering can create a selective dependency in cancer.

Engineered CRISPR-Cas12a for higher-order combinatorial chromatin perturbations

Multiplexed genetic perturbations are critical for testing functional interactions among coding or non-coding genetic elements. Compared to double-stranded DNA cutting, repressive chromatin formation using CRISPR interference (CRISPRi) avoids genotoxicity and is more effective for perturbing non-coding regulatory elements in pooled assays. However, current CRISPRi pooled screening approaches are limited to targeting one to three genomic sites per cell. We engineer an Acidaminococcus Cas12a (AsCas12a) variant, multiplexed transcriptional interference AsCas12a (multiAsCas12a), that incorporates R1226A, a mutation that stabilizes the ribonucleoprotein-DNA complex via DNA nicking. The multiAsCas12a-KRAB fusion improves CRISPRi activity over DNase-dead AsCas12a-KRAB fusions, often rescuing the activities of lentivirally delivered CRISPR RNAs (crRNA) that are inactive when used with the latter. multiAsCas12a-KRAB supports CRISPRi using 6-plex crRNA arrays in high-throughput pooled screens. Using multiAsCas12a-KRAB, we discover enhancer elements and dissect the combinatorial function of cis-regulatory elements in human cells. These results instantiate a group testing framework for efficiently surveying numerous combinations of chromatin perturbations for biological discovery and engineering.

A collection of split-Gal4 drivers targeting conserved signaling ligands in Drosophila

Communication between cells in metazoan organisms is mediated by a remarkably small number of highly conserved signaling pathways. Given this small number of signaling pathways, the existence of multiple related ligands for many of these pathways represents a key evolutionary innovation for encoding complexity into cell-cell signaling. Relatedly, crosstalk between pathways is another critical feature, which allows a modest number of pathways to ultimately generate an enormously diverse range of outcomes. It would thus be useful to have genetic tools to identify and manipulate not only those cells that express a given signaling ligand but also those cells that specifically coexpress pairs of signaling ligands. We present a collection of split-Gal4 knock-in lines targeting many of the ligands for highly conserved signaling pathways in Drosophila (Notch, Hedgehog, fibroblast growth factor (FGF), epidermal growth factor (EGF), transforming growth factor β (TGFβ), Janus kinase/signal transducer and activator of transcription (JAK/STAT), Jun kinase (JNK), and platelet-derived growth factor (PDGF)/vascular endothelial growth factor (VEGF)-related receptor (PVR). We demonstrate that these lines faithfully recapitulate the endogenous expression pattern of their targets and that they can be used to identify cells and tissues that coexpress pairs of ligands. As a proof of principle, we demonstrate that the 4th chromosome TGFβ ligands myoglianin and maverick are broadly coexpressed in muscles and other tissues of both larva and adults and that the JAK/STAT ligands upd2 and upd3 are partially coexpressed from cells of the midgut following gut damage. Together with our previously collection of split-Gal4 lines targeting the 7 Wnt ligands, this resource allows Drosophila researchers to identify and genetically manipulate cells that specifically express pairs of conserved ligands from nearly all the major intercellular signaling pathways.

A collection of split-Gal4 drivers targeting conserved signaling ligands in Drosophila

Communication between cells in metazoan organisms is mediated by a remarkably small number of highly conserved signaling pathways. Given the relatively small number of signaling pathways, the existence of multiple related ligands for many of these pathways is thought to represent a key evolutionary innovation for encoding complexity into cell-cell signaling. Relatedly, crosstalk and other interactions between pathways is another critical feature which allows a modest number pathways to ultimately generate an enormously diverse range of outcomes. It would thus be useful to have genetic tools to identify and manipulate not only those cells which express a given signaling ligand, but also those cells that specifically co-express pairs of signaling ligands. Here, we present a collection of split-Gal4 knock-in lines targeting many of the ligands for highly conserved signaling pathways in Drosophila (Notch, Hedgehog, FGF, EGF, TGFβ, JAK/STAT, JNK, and PVR). We demonstrate that these lines faithfully recapitulate the endogenous expression pattern of their targets, and that they can be used to specifically identify the cells and tissues that co-express pairs of signaling ligands. As a proof of principle, we demonstrate that the 4th chromosome TGFβ ligands myoglianin and maverick are broadly co-expressed in muscles and other tissues of both larva and adults, and that the JAK/STAT ligands upd2 and upd3 are partially co-expressed from cells of the midgut following gut damage. Together with our previously collection of split-Gal4 lines targeting the seven Wnt ligands, this resource allows Drosophila researchers to identify and genetically manipulate cells that specifically express pairs of conserved ligands from nearly all the major intercellular signaling pathways.

Higher resolution pooled genome-wide CRISPR knockout screening in Drosophila cells using Integration and Anti-CRISPR (IntAC)

CRISPR screens enable systematic, scalable genotype-to-phenotype mapping. We previously developed a pooled CRISPR screening method for Drosophila melanogaster and mosquito cell lines using plasmid transfection and site-specific integration to introduce single guide (sgRNA) libraries, followed by PCR and sequencing of integrated sgRNAs. While effective, the method relies on early constitutive Cas9 activity that potentially can lead to discrepancies between genome edits and sgRNAs detected by PCR, reducing screen accuracy. To address this issue, we introduce a new method to co-transfect a plasmid expressing the anti-CRISPR protein AcrIIa4 to suppress Cas9 activity during early sgRNA expression, which we term "IntAC" (integrase with anti-CRISPR). IntAC allowed us to construct a new CRISPR screening approach driven by the high strength dU6:3 promoter. This new library dramatically improved precision-recall of fitness genes across the genome, retrieving 90-95% of essential gene groups within 5% error, allowing us to generate the most comprehensive list of cell fitness genes yet assembled for Drosophila. Our analysis determined that elevated sgRNA levels, made permissible by the IntAC approach, drove much of the improvement. The Drosophila fitness genes show strong correlation with human fitness genes and underscore the effects of paralogs on gene essentiality. We further demonstrate that IntAC combined with a targeted sgRNA sub-library enabled precise positive selection of a transporter under solute overload. IntAC represents a straightforward enhancement to existing Drosophila CRISPR screening methods, dramatically increasing accuracy, and might also be broadly applicable to virus-free CRISPR screens in other cell types, including mosquito, lepidopteran, tick, and mammalian cells.

A genetic screen of transcription factors in the Drosophila melanogaster abdomen identifies novel pigmentation genes

Gene regulatory networks specify the gene expression patterns needed for traits to develop. Differences in these networks can result in phenotypic differences between organisms. Although loss-of-function genetic screens can identify genes necessary for trait formation, gain-of-function screens can overcome genetic redundancy and identify loci whose expression is sufficient to alter trait formation. Here, we leveraged transgenic lines from the Transgenic RNAi Project at Harvard Medical School to perform both gain- and loss-of-function CRISPR/Cas9 screens for abdominal pigmentation phenotypes. We identified measurable effects on pigmentation patterns in the Drosophila melanogaster abdomen for 21 of 55 transcription factors in gain-of-function experiments and 7 of 16 tested by loss-of-function experiments. These included well-characterized pigmentation genes, such as bab1 and dsx, and transcription factors that had no known role in pigmentation, such as slp2. Finally, this screen was partially conducted by undergraduate students in a Genetics Laboratory course during the spring semesters of 2021 and 2022. We found this screen to be a successful model for student engagement in research in an undergraduate laboratory course that can be readily adapted to evaluate the effect of hundreds of genes on many different Drosophila traits, with minimal resources.

The TSC22D, WNK, and NRBP gene families exhibit functional buffering and evolved with Metazoa for cell volume regulation

The ability to sense and respond to osmotic fluctuations is critical for the maintenance of cellular integrity. We used gene co-essentiality analysis to identify an unappreciated relationship between TSC22D2, WNK1, and NRBP1 in regulating cell volume homeostasis. All of these genes have paralogs and are functionally buffered for osmo-sensing and cell volume control. Within seconds of hyperosmotic stress, TSC22D, WNK, and NRBP family members physically associate into biomolecular condensates, a process that is dependent on intrinsically disordered regions (IDRs). A close examination of these protein families across metazoans revealed that TSC22D genes evolved alongside a domain in NRBPs that specifically binds to TSC22D proteins, which we have termed NbrT (NRBP binding region with TSC22D), and this co-evolution is accompanied by rapid IDR length expansion in WNK-family kinases. Our study reveals that TSC22D, WNK, and NRBP genes evolved in metazoans to co-regulate rapid cell volume changes in response to osmolarity.

Lhx4 surpasses its paralog Lhx3 in promoting the differentiation of spinal V2a interneurons

Paralog factors are considered to ensure the robustness of biological processes by providing redundant activity in cells where they are co-expressed. However, the specific contribution of each factor is frequently underestimated. In the developing spinal cord, multiple families of transcription factors successively contribute to differentiate an initially homogenous population of neural progenitors into a myriad of neuronal subsets with distinct molecular, morphological, and functional characteristics. The LIM-homeodomain transcription factors Lhx3, Lhx4, Isl1 and Isl2 promote the segregation and differentiation of spinal motor neurons and V2 interneurons. Based on their high sequence identity and their similar distribution, the Lhx3 and Lhx4 paralogs are considered to contribute similarly to these processes. However, the specific contribution of Lhx4 has never been studied. Here, we provide evidence that Lhx3 and Lhx4 are present in the same cell populations during spinal cord development. Similarly to Lhx3, Lhx4 can form multiproteic complexes with Isl1 or Isl2 and the nuclear LIM interactor NLI. Lhx4 can stimulate a V2-specific enhancer more efficiently than Lhx3 and surpasses Lhx3 in promoting the differentiation of V2a interneurons in chicken embryo electroporation experiments. Finally, Lhx4 inactivation in mice results in alterations of differentiation of the V2a subpopulation, but not of motor neuron production, suggesting that Lhx4 plays unique roles in V2a differentiation that are not compensated by the presence of Lhx3. Thus, Lhx4 could be the major LIM-HD factor involved in V2a interneuron differentiation during spinal cord development and should be considered for in vitro differentiation of spinal neuronal populations.

Efficient gene knockout and genetic interaction screening using the in4mer CRISPR/Cas12a multiplex knockout platform

Genetic interactions mediate the emergence of phenotype from genotype, but technologies for combinatorial genetic perturbation in mammalian cells are challenging to scale. Here, we identify background-independent paralog synthetic lethals from previous CRISPR genetic interaction screens, and find that the Cas12a platform provides superior sensitivity and assay replicability. We develop the in4mer Cas12a platform that uses arrays of four independent guide RNAs targeting the same or different genes. We construct a genome-scale library, Inzolia, that is ~30% smaller than a typical CRISPR/Cas9 library while also targeting ~4000 paralog pairs. Screens in cancer cells demonstrate discrimination of core and context-dependent essential genes similar to that of CRISPR/Cas9 libraries, as well as detection of synthetic lethal and masking/buffering genetic interactions between paralogs of various family sizes. Importantly, the in4mer platform offers a fivefold reduction in library size compared to other genetic interaction methods, substantially reducing the cost and effort required for these assays.

Higher-order combinatorial chromatin perturbations by engineered CRISPR-Cas12a for functional genomics

Multiplexed genetic perturbations are critical for testing functional interactions among coding or non-coding genetic elements. Compared to double-stranded DNA cutting, repressive chromatin formation using CRISPR interference (CRISPRi) avoids genotoxicity and is more effective for perturbing non-coding regulatory elements in pooled assays. However, current CRISPRi pooled screening approaches are limited to targeting 1-3 genomic sites per cell. To develop a tool for higher-order ( > 3) combinatorial targeting of genomic sites with CRISPRi in functional genomics screens, we engineered an Acidaminococcus Cas12a variant -- referred to as mul tiplexed transcriptional interference AsCas12a (multiAsCas12a). multiAsCas12a incorporates a key mutation, R1226A, motivated by the hypothesis of nicking-induced stabilization of the ribonucleoprotein:DNA complex for improving CRISPRi activity. multiAsCas12a significantly outperforms prior state-of-the-art Cas12a variants in combinatorial CRISPRi targeting using high-order multiplexed arrays of lentivirally transduced CRISPR RNAs (crRNA), including in high-throughput pooled screens using 6-plex crRNA array libraries. Using multiAsCas12a CRISPRi, we discover new enhancer elements and dissect the combinatorial function of cis-regulatory elements. These results instantiate a group testing framework for efficiently surveying potentially numerous combinations of chromatin perturbations for biological discovery and engineering.

Iterative editing of multiple genes using CRISPR/Cas9 in C. elegans

Certain sets of genes are derived from gene duplication and share substantial sequence similarity in C. elegans , presenting a significant challenge in determining the specific roles of each gene and their collective impact on cellular processes. Here, we show that a collection of genes can be disrupted in a single animal via multiple rounds of CRISPR/Cas9 mediated genome editing. We found that up to three genes can be simultaneously disrupted in a single editing event with high efficiency. Our approach offers an opportunity to explore the genetic interaction and molecular underpinning of gene clusters with redundant function.

Tackling redundancy: genetic mechanisms underlying paralog compensation in plants

Gene duplication is a powerful source of biological innovation giving rise to paralogous genes that undergo diverse fates. Redundancy between paralogous genes is an intriguing outcome of duplicate gene evolution, and its maintenance over evolutionary time has long been considered a paradox. Redundancy can also be dubbed 'a geneticist's nightmare': It hinders the predictability of genome editing outcomes and limits our ability to link genotypes to phenotypes. Genetic studies in yeast and plants have suggested that the ability of ancient redundant duplicates to compensate for dosage perturbations resulting from a loss of function depends on the reprogramming of gene expression, a phenomenon known as active compensation. Starting from considerations on the stoichiometric constraints that drive the evolutionary stability of redundancy, this review aims to provide insights into the mechanisms of active compensation between duplicates that could be targeted for breaking paralog dependencies - the next frontier in plant functional studies.

Pooled genome-wide CRISPR activation screening for rapamycin resistance genes in Drosophila cells

Loss-of-function and gain-of-function genetic perturbations provide valuable insights into gene function. In Drosophila cells, while genome-wide loss-of-function screens have been extensively used to reveal mechanisms of a variety of biological processes, approaches for performing genome-wide gain-of-function screens are still lacking. Here, we describe a pooled CRISPR activation (CRISPRa) screening platform in Drosophila cells and apply this method to both focused and genome-wide screens to identify rapamycin resistance genes. The screens identified three genes as novel rapamycin resistance genes: a member of the SLC16 family of monocarboxylate transporters (CG8468), a member of the lipocalin protein family (CG5399), and a zinc finger C2H2 transcription factor (CG9932). Mechanistically, we demonstrate that CG5399 overexpression activates the RTK-Akt-mTOR signaling pathway and that activation of insulin receptor (InR) by CG5399 requires cholesterol and clathrin-coated pits at the cell membrane. This study establishes a novel platform for functional genetic studies in Drosophila cells.

Paralog Explorer: A resource for mining information about paralogs in common research organisms

Paralogs are genes which arose via gene duplication, and when such paralogs retain overlapping or redundant function, this poses a challenge to functional genetics research. Recent technological advancements have made it possible to systematically probe gene function for redundant genes using dual or multiplex gene perturbation, and there is a need for a simple bioinformatic tool to identify putative paralogs of a gene(s) of interest. We have developed Paralog Explorer (https://www.flyrnai.org/tools/paralogs/), an online resource that allows researchers to quickly and accurately identify candidate paralogous genes in the genomes of the model organisms D. melanogaster, C. elegans, D. rerio, M. musculus, and H. sapiens. Paralog Explorer deploys an effective between-species ortholog prediction software, DIOPT, to analyze within-species paralogs. Paralog Explorer allows users to identify candidate paralogs, and to navigate relevant databases regarding gene co-expression, protein-protein and genetic interaction, as well as gene ontology and phenotype annotations. Altogether, this tool extends the value of current ortholog prediction resources by providing sophisticated features useful for identification and study of paralogous genes.

Trans-omics analysis of insulin action reveals a cell growth subnetwork which co-regulates anabolic processes

Insulin signaling promotes anabolic metabolism to regulate cell growth through multi-omic interactions. To obtain a comprehensive view of the cellular responses to insulin, we constructed a trans-omic network of insulin action in Drosophila cells that involves the integration of multi-omic data sets. In this network, 14 transcription factors, including Myc, coordinately upregulate the gene expression of anabolic processes such as nucleotide synthesis, transcription, and translation, consistent with decreases in metabolites such as nucleotide triphosphates and proteinogenic amino acids required for transcription and translation. Next, as cell growth is required for cell proliferation and insulin can stimulate proliferation in a context-dependent manner, we integrated the trans-omic network with results from a CRISPR functional screen for cell proliferation. This analysis validates the role of a Myc-mediated subnetwork that coordinates the activation of genes involved in anabolic processes required for cell growth.

•

A trans-omic network of insulin action in Drosophila cells was constructed

•

Insulin co-regulates various anabolic processes in a time-dependent manner

•

The trans-omic network and a CRISPR screen for cell proliferation were integrated

•

A Myc-mediated subnetwork promoting anabolic processes is required for cell growth

Cross-species prediction of essential genes in insects

MOTIVATION

Insects possess a vast phenotypic diversity and key ecological roles. Several insect species also have medical, agricultural and veterinary importance as parasites and disease vectors. Therefore, strategies to identify potential essential genes in insects may reduce the resources needed to find molecular players in central processes of insect biology. However, most predictors of essential genes in multicellular eukaryotes using machine learning rely on expensive and laborious experimental data to be used as gene features, such as gene expression profiles or protein-protein interactions, even though some of this information may not be available for the majority of insect species with genomic sequences available.

RESULTS

Here, we present and validate a machine learning strategy to predict essential genes in insects using sequence-based intrinsic attributes (statistical and physicochemical data) together with the predictions of subcellular location and transcriptomic data, if available. We gathered information available in public databases describing essential and non-essential genes for Drosophila melanogaster (fruit fly, Diptera) and Tribolium castaneum (red flour beetle, Coleoptera). We proceeded by computing intrinsic and extrinsic attributes that were used to train statistical models in one species and tested by their capability of predicting essential genes in the other. Even models trained using only intrinsic attributes are capable of predicting genes in the other insect species, including the prediction of lineage-specific essential genes. Furthermore, the inclusion of RNA-Seq data is a major factor to increase classifier performance.

AVAILABILITY AND IMPLEMENTATION

The code, data and final models produced in this study are freely available at https://github.com/g1o/GeneEssentiality/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The Arabidopsis AAC Proteins CIL and CIA2 Are Sub-functionalized Paralogs Involved in Chloroplast Development

The Arabidopsis gene Chloroplast Import Apparatus 2 (CIA2) encodes a transcription factor that positively affects the activity of nuclear genes for chloroplast ribosomal proteins and chloroplast protein import machineries. CIA2-like (CIL) is the paralogous gene of CIA2. We generated a cil mutant by site-directed mutagenesis and compared it with cia2 and cia2cil double mutant. Phenotype of the cil mutant did not differ from the wild type under our growth conditions, except faster growth and earlier time to flowering. Compared to cia2, the cia2cil mutant showed more impaired chloroplast functions and reduced amounts of plastid ribosomal RNAs. In silico analyses predict for CIA2 and CIL a C-terminal CCT domain and an N-terminal chloroplast transit peptide (cTP). Chloroplast (and potentially nuclear) localization was previously shown for HvCMF3 and HvCMF7, the homologs of CIA2 and CIL in barley. We observed nuclear localization of CIL after transient expression in Arabidopsis protoplasts. Surprisingly, transformation of cia2 with HvCMF3, HvCMF7, or with a truncated CIA2 lacking the predicted cTP could partially rescue the pale-green phenotype of cia2. These data are discussed with respect to potentially overlapping functions between CIA2, CIL, and their barley homologs and to the function of the putative cTPs of CIA2 and CIL.

Using CRISPR-Cas9-based genome engineering tools in Drosophila melanogaster

Drosophila melanogaster has been used as a model organism for over a century. Mutant-based analyses have been used extensively to understand the genetic basis of different cellular processes, including development, neuronal function and diseases. Most of the earlier genetic mutants and specific tools were generated by random insertions and deletion strategies and then mapped to specific genomic loci. Since all genomic regions are not equally accessible to random mutations and insertions, many genes still remain uncharacterized. Low efficiency of targeted genomic manipulation approaches that rely on homologous recombination, and difficulty in generating resources for sequence-specific endonucleases, such as ZFNs (Zinc Finger Nucleases) and TALENs (Transcription Activator-Like Effector Nucleases), could not make these gene targeting techniques very popular. However, recently RNA directed DNA endonucleases, such as CRISPR-Cas, have transformed genome engineering owing to their comparative ease, versatility, and low expense. With the added advantage of preexisting genetic tools, CRISPR-Cas-based manipulations are being extensively used in Drosophila melanogaster and simultaneously being fine-tuned for specific experimental requirements. In this chapter, I will discuss various uses of CRISPR-Cas-based genetic engineering and specific design methods in Drosophila melanogaster. I will summarize various already available tools that are being utilized in conjunction with CRISPR-Cas technology to generate specific genetic manipulation and are being optimized to address specific questions. Finally, I will discuss the future directions of Drosophila genetics research and how CRISPR-Cas can be utilized to target specific questions, addressing which has not been possible thus far.

Functional interrogation of a SARS-CoV-2 host protein interactome identifies unique and shared coronavirus host factors

The ongoing severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic has devastated the global economy and claimed more than 1.7 million lives, presenting an urgent global health crisis. To identify host factors required for infection by SARS-CoV-2 and seasonal coronaviruses, we designed a focused high-coverage CRISPR-Cas9 library targeting 332 members of a recently published SARS-CoV-2 protein interactome. We leveraged the compact nature of this library to systematically screen SARS-CoV-2 at two physiologically relevant temperatures along with three related coronaviruses (human coronavirus 229E [HCoV-229E], HCoV-NL63, and HCoV-OC43), allowing us to probe this interactome at a much higher resolution than genome-scale studies. This approach yielded several insights, including potential virus-specific differences in Rab GTPase requirements and glycosylphosphatidylinositol (GPI) anchor biosynthesis, as well as identification of multiple pan-coronavirus factors involved in cholesterol homeostasis. This coronavirus essentiality catalog could inform ongoing drug development efforts aimed at intercepting and treating coronavirus disease 2019 (COVID-19) and help prepare for future coronavirus outbreaks.

Cross-species identification of PIP5K1-, splicing- and ubiquitin-related pathways as potential targets for RB1-deficient cells

The RB1 tumor suppressor is recurrently mutated in a variety of cancers including retinoblastomas, small cell lung cancers, triple-negative breast cancers, prostate cancers, and osteosarcomas. Finding new synthetic lethal (SL) interactions with RB1 could lead to new approaches to treating cancers with inactivated RB1. We identified 95 SL partners of RB1 based on a Drosophila screen for genetic modifiers of the eye phenotype caused by defects in the RB1 ortholog, Rbf1. We validated 38 mammalian orthologs of Rbf1 modifiers as RB1 SL partners in human cancer cell lines with defective RB1 alleles. We further show that for many of the RB1 SL genes validated in human cancer cell lines, low activity of the SL gene in human tumors, when concurrent with low levels of RB1 was associated with improved patient survival. We investigated higher order combinatorial gene interactions by creating a novel Drosophila cancer model with co-occurring Rbf1, Pten and Ras mutations, and found that targeting RB1 SL genes in this background suppressed the dramatic tumor growth and rescued fly survival whilst having minimal effects on wild-type cells. Finally, we found that drugs targeting the identified RB1 interacting genes/pathways, such as UNC3230, PYR-41, TAK-243, isoginkgetin, madrasin, and celastrol also elicit SL in human cancer cell lines. In summary, we identified several high confidence, evolutionarily conserved, novel targets for RB1-deficient cells that may be further adapted for the treatment of human cancer.

Investigating rare and ultrarare epilepsy syndromes with Drosophila models

One in three epilepsy cases is drug resistant, and seizures often begin in infancy, when they are life-threatening and when therapeutic options are highly limited. An important tool for prioritizing and validating genes associated with epileptic conditions, which is suitable for large-scale screening, is disease modeling in Drosophila. Approximately two-thirds of disease genes are conserved in Drosophila, and gene-specific fly models exhibit behavioral changes that are related to symptoms of epilepsy. Models are based on behavior readouts, seizure-like attacks and paralysis following stimulation, and neuronal, cell-biological readouts that are in the majority based on changes in nerve cell activity or morphology. In this review, we focus on behavioral phenotypes. Importantly, Drosophila modeling is independent of, and complementary to, other approaches that are computational and based on systems analysis. The large number of known epilepsy-associated gene variants indicates a need for efficient research strategies. We will discuss the status quo of epilepsy disease modelling in Drosophila and describe promising steps towards the development of new drugs to reduce seizure rates and alleviate other epileptic symptoms.

Genome-Scale Identification of SARS-CoV-2 and Pan-coronavirus Host Factor Networks

The coronavirus disease 2019 (COVID-19) pandemic has claimed the lives of over one million people worldwide. The causative agent, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is a member of the Coronaviridae family of viruses that can cause respiratory infections of varying severity. The cellular host factors and pathways co-opted during SARS-CoV-2 and related coronavirus life cycles remain ill defined. To address this gap, we performed genome-scale CRISPR knockout screens during infection by SARS-CoV-2 and three seasonal coronaviruses (HCoV-OC43, HCoV-NL63, and HCoV-229E). These screens uncovered host factors and pathways with pan-coronavirus and virus-specific functional roles, including major dependency on glycosaminoglycan biosynthesis, sterol regulatory element-binding protein (SREBP) signaling, bone morphogenetic protein (BMP) signaling, and glycosylphosphatidylinositol biosynthesis, as well as a requirement for several poorly characterized proteins. We identified an absolute requirement for the VMP1, TMEM41, and TMEM64 (VTT) domain-containing protein transmembrane protein 41B (TMEM41B) for infection by SARS-CoV-2 and three seasonal coronaviruses. This human coronavirus host factor compendium represents a rich resource to develop new therapeutic strategies for acute COVID-19 and potential future coronavirus pandemics.

Genome-scale identification of SARS-CoV-2 and pan-coronavirus host factor networks

The COVID-19 pandemic has claimed the lives of more than one million people worldwide. The causative agent, SARS-CoV-2, is a member of the Coronaviridae family, which are viruses that cause respiratory infections of varying severity. The cellular host factors and pathways co-opted by SARS-CoV-2 and other coronaviruses in the execution of their life cycles remain ill-defined. To develop an extensive compendium of host factors required for infection by SARS-CoV-2 and three seasonal coronaviruses (HCoV-OC43, HCoV-NL63, and HCoV-229E), we performed parallel genome-scale CRISPR knockout screens. These screens uncovered multiple host factors and pathways with pan-coronavirus and virus-specific functional roles, including major dependency on glycosaminoglycan biosynthesis, SREBP signaling, and glycosylphosphatidylinositol biosynthesis, as well as an unexpected requirement for several poorly characterized proteins. We identified an absolute requirement for the VTT-domain containing protein TMEM41B for infection by SARS-CoV-2 and all other coronaviruses. This human Coronaviridae host factor compendium represents a rich resource to develop new therapeutic strategies for acute COVID-19 and potential future coronavirus spillover events.

HIGHLIGHTS

Genome-wide CRISPR screens for SARS-CoV-2, HCoV-OC43, HCoV-NL63, and HCoV-229E coronavirus host factors.Parallel genome-wide CRISPR screening uncovered host factors and pathways with pan-coronavirus and virus-specific functional roles.Coronaviruses co-opt multiple biological pathways, including glycosaminoglycan biosynthesis, SREBP signaling, and glycosylphosphatidylinositol biosynthesis and anchoring, among others.TMEM41B - a poorly understood factor with roles in autophagy and lipid mobilization - is a critical pan-coronavirus host factor.

Vsx1 and Chx10 paralogs sequentially secure V2 interneuron identity during spinal cord development

Paralog factors are usually described as consolidating biological systems by displaying redundant functionality in the same cells. Here, we report that paralogs can also cooperate in distinct cell populations at successive stages of differentiation. In mouse embryonic spinal cord, motor neurons and V2 interneurons differentiate from adjacent progenitor domains that share identical developmental determinants. Therefore, additional strategies secure respective cell fate. In particular, Hb9 promotes motor neuron identity while inhibiting V2 differentiation, whereas Chx10 stimulates V2a differentiation while repressing motor neuron fate. However, Chx10 is not present at the onset of V2 differentiation and in other V2 populations. In the present study, we show that Vsx1, the single paralog of Chx10, which is produced earlier than Chx10 in V2 precursors, can inhibit motor neuron differentiation and promote V2 interneuron production. However, the single absence of Vsx1 does not impact on V2 fate consolidation, suggesting that lack of Vsx1 may be compensated by other factors. Nevertheless, Vsx1 cooperates with Chx10 to prevent motor neuron differentiation in early V2 precursors although these two paralog factors are not produced in the same cells. Hence, this study uncovers an original situation, namely labor division, wherein paralog genes cooperate at successive steps of neuronal development.

Functional interrogation of a SARS-CoV-2 host protein interactome identifies unique and shared coronavirus host factors

The ongoing SARS-CoV-2 pandemic has devastated the global economy and claimed nearly one million lives, presenting an urgent global health crisis. To identify host factors required for infection by SARS-CoV-2 and seasonal coronaviruses, we designed a focused high-coverage CRISPR-Cas9 library targeting 332 members of a recently published SARS-CoV-2 protein interactome. We leveraged the compact nature of this library to systematically screen four related coronaviruses (HCoV-229E, HCoV-NL63, HCoV-OC43 and SARS-CoV-2) at two physiologically relevant temperatures (33 °C and 37 °C), allowing us to probe this interactome at a much higher resolution relative to genome scale studies. This approach yielded several new insights, including unexpected virus and temperature specific differences in Rab GTPase requirements and GPI anchor biosynthesis, as well as identification of multiple pan-coronavirus factors involved in cholesterol homeostasis. This coronavirus essentiality catalog could inform ongoing drug development efforts aimed at intercepting and treating COVID-19, and help prepare for future coronavirus outbreaks.

HIGHLIGHTS

Focused CRISPR screens targeting host factors in the SARS-CoV-2 interactome were performed for SARS-CoV-2, HCoV-229E, HCoV-NL63, and HCoV-OC43 coronaviruses.Focused interactome CRISPR screens achieve higher resolution compared to genome-wide screens, leading to the identification of critical factors missed by the latter.Parallel CRISPR screens against multiple coronaviruses uncover host factors and pathways with pan-coronavirus and virus-specific functional roles.The number of host proteins that interact with a viral bait protein is not proportional to the number of functional interactors.Novel SARS-CoV-2 host factors are expressed in relevant cell types in the human airway.

Combined Proteomic and Genetic Interaction Mapping Reveals New RAS Effector Pathways and Susceptibilities

Activating mutations in RAS GTPases drive many cancers, but limited understanding of less-studied RAS interactors, and of the specific roles of different RAS interactor paralogs, continues to limit target discovery. We developed a multistage discovery and screening process to systematically identify genes conferring RAS-related susceptibilities in lung adenocarcinoma. Using affinity purification mass spectrometry, we generated a protein-protein interaction map of RAS interactors and pathway components containing hundreds of interactions. From this network, we constructed a CRISPR dual knockout library targeting 119 RAS-related genes that we screened for KRAS-dependent genetic interactions (GI). This approach identified new RAS effectors, including the adhesion controller RADIL and the endocytosis regulator RIN1, and >250 synthetic lethal GIs, including a potent KRAS-dependent interaction between RAP1GDS1 and RHOA. Many GIs link specific paralogs within and between gene families. These findings illustrate the power of multiomic approaches to uncover synthetic lethal combinations specific for hitherto untreatable cancer genotypes. SIGNIFICANCE: We establish a deep network of protein-protein and genetic interactions in the RAS pathway. Many interactions validated here demonstrate important specificities and redundancies among paralogous RAS regulators and effectors. By comparing synthetic lethal interactions across KRAS-dependent and KRAS-independent cell lines, we identify several new combination therapy targets for RAS-driven cancers.This article is highlighted in the In This Issue feature, p. 1775.

The FOXO's Advantages of Being a Family: Considerations on Function and Evolution

The nematode Caenorhabditis elegans possesses a unique (with various isoforms) FOXO transcription factor DAF-16, which is notorious for its role in aging and its regulation by the insulin-PI3K-AKT pathway. In humans, five genes (including a protein-coding pseudogene) encode for FOXO transcription factors that are targeted by the PI3K-AKT axis, such as in C. elegans. This common regulation and highly conserved DNA-binding domain are the pillars of this family. In this review, I will discuss the possible meaning of possessing a group of very similar proteins and how it can generate additional functionality to more complex organisms. I frame this discussion in relation to the much larger super family of Forkhead proteins to which they belong. FOXO members are very often co-expressed in the same cell type. The overlap of function and expression creates a certain redundancy that might be a safeguard against the accidental loss of FOXO function, which could otherwise lead to disease, particularly, cancer. This is one of the points that will be examined in this "family affair" report.

Genetic interaction mapping and exon-resolution functional genomics with a hybrid Cas9-Cas12a platform

Systematic mapping of genetic interactions (GIs) and interrogation of the functions of sizable genomic segments in mammalian cells represent important goals of biomedical research. To advance these goals, we present a CRISPR (clustered regularly interspaced short palindromic repeats)-based screening system for combinatorial genetic manipulation that employs coexpression of CRISPR-associated nucleases 9 and 12a (Cas9 and Cas12a) and machine-learning-optimized libraries of hybrid Cas9-Cas12a guide RNAs. This system, named Cas Hybrid for Multiplexed Editing and screening Applications (CHyMErA), outperforms genetic screens using Cas9 or Cas12a editing alone. Application of CHyMErA to the ablation of mammalian paralog gene pairs reveals extensive GIs and uncovers phenotypes normally masked by functional redundancy. Application of CHyMErA in a chemogenetic interaction screen identifies genes that impact cell growth in response to mTOR pathway inhibition. Moreover, by systematically targeting thousands of alternative splicing events, CHyMErA identifies exons underlying human cell line fitness. CHyMErA thus represents an effective screening approach for GI mapping and the functional analysis of sizable genomic regions, such as alternative exons.

Large-Scale Transgenic Drosophila Resource Collections for Loss- and Gain-of-Function Studies

The Transgenic RNAi Project (TRiP), a Drosophila melanogaster functional genomics platform at Harvard Medical School, was initiated in 2008 to generate and distribute a genome-scale collection of RNA interference (RNAi) fly stocks. To date, it has generated >15,000 RNAi fly stocks. As this covers most Drosophila genes, we have largely transitioned to development of new resources based on CRISPR technology. Here, we present an update on our libraries of publicly available RNAi and CRISPR fly stocks, and focus on the TRiP-CRISPR overexpression (TRiP-OE) and TRiP-CRISPR knockout (TRiP-KO) collections. TRiP-OE stocks express single guide RNAs targeting upstream of a gene transcription start site. Gene activation is triggered by coexpression of catalytically dead Cas9 fused to an activator domain, either VP64-p65-Rta or Synergistic Activation Mediator. TRiP-KO stocks express one or two single guide RNAs targeting the coding sequence of a gene or genes. Cutting is triggered by coexpression of Cas9, allowing for generation of indels in both germline and somatic tissue. To date, we have generated >5000 TRiP-OE or TRiP-KO stocks for the community. These resources provide versatile, transformative tools for gene activation, gene repression, and genome engineering.

CRISPR-induced double-strand breaks trigger recombination between homologous chromosome arms

CRISPR-Cas9-based genome editing has transformed the life sciences, enabling virtually unlimited genetic manipulation of genomes: The RNA-guided Cas9 endonuclease cuts DNA at a specific target sequence and the resulting double-strand breaks are mended by one of the intrinsic cellular repair pathways. Imprecise double-strand repair will introduce random mutations such as indels or point mutations, whereas precise editing will restore or specifically edit the locus as mandated by an endogenous or exogenously provided template. Recent studies indicate that CRISPR-induced DNA cuts may also result in the exchange of genetic information between homologous chromosome arms. However, conclusive data of such recombination events in higher eukaryotes are lacking. Here, we show that in Drosophila, the detected Cas9-mediated editing events frequently resulted in germline-transmitted exchange of chromosome arms-often without indels. These findings demonstrate the feasibility of using the system for generating recombinants and also highlight an unforeseen risk of using CRISPR-Cas9 for therapeutic intervention.

Evolutionary history of the medaka long-wavelength sensitive genes and effects of artificial regression by gene loss on behavioural photosensitivity

Tandem gene duplication has led to an expansion of cone-opsin repertoires in many fish, but the resulting functional advantages have only been conjectured without empirical demonstration. Medaka (Oryzias latipes and O. sakaizumii) have eight (two red, three green, two blue, and one violet) cone opsin genes. Absorbance maxima (λ_max) of the proteins vary from 356 nm to 562 nm, but those of the red opsins (long-wavelength sensitive; LWS) are nearly identical, obscuring the necessity of their coexistence. Here, we compared the LWSa and LWSb loci of these sister species and found that the gene duplication occurred long before the latipes–sakaizumii speciation (4–18 million years ago), and the high sequence similarity between the paralogues is the result of at least two events of gene conversion. These repetitive gene conversions would indicate the importance for medaka of retaining two identical LWSs in the genome. However, a newly established medaka mutant with a single LWS showed no defect in LWS expression or behavioural red-light sensitivity, demonstrating functional redundancy of the paralogs. Thus, as with many other genes after whole-genome duplication, the redundant LWS might be on the way to being lost from the current cone opsin repertoire. Thus, non-allelic gene conversion may temporarily provide an easier and more frequent solution than gene loss for reducing genetic diversity, which should be considered when assessing history of gene evolution by phylogenetic analyses.

Ca2+ concentration-dependent premature death of igfbp5a-/- fish reveals a critical role of IGF signaling in adaptive epithelial growth

The phenotype gap is a challenge for genetically dissecting redundant endocrine signaling pathways, such as the six isoforms in the insulin-like growth factor binding protein (IGFBP) family. Although overexpressed IGFBPs can inhibit or potentiate IGF actions or have IGF-independent actions, mutant mice lacking IGFBP-encoding genes do not exhibit major phenotypes. We found that although zebrafish deficient in igfbp5a did not show overt phenotypes when raised in Ca²⁺-rich solutions, they died prematurely in low Ca²⁺ conditions. A group of epithelial cells expressing igfbp5a take up Ca²⁺ and proliferate under low Ca²⁺ conditions because of activation of IGF signaling. Deletion of igfbp5a blunted low Ca²⁺ stress-induced IGF signaling and impaired adaptive proliferation. Reintroducing zebrafish Igfbp5a, but not its ligand binding-deficient mutant, restored adaptive proliferation. Similarly, adaptive proliferation was restored in zebrafish lacking igfbp5a by expression of human IGFBP5, but not two cancer-associated IGFBP5 mutants. Knockdown of IGFBP5 in human colon carcinoma cells resulted in reduced IGF-stimulated cell proliferation. These results reveal a conserved mechanism by which a locally expressed Igfbp regulates organismal Ca²⁺ homeostasis and survival by activating IGF signaling in epithelial cells and promoting their proliferation in Ca²⁺-deficient states. These findings underscore the importance of physiological context when analyzing loss-of-function phenotypes of endocrine factors.

Discovering Health Benefits of Phytochemicals with Integrated Analysis of the Molecular Network, Chemical Properties and Ethnopharmacological Evidence

Identifying the health benefits of phytochemicals is an essential step in drug and functional food development. While many in vitro screening methods have been developed to identify the health effects of phytochemicals, there is still room for improvement because of high cost and low productivity. Therefore, researchers have alternatively proposed in silico methods, primarily based on three types of approaches; utilizing molecular, chemical or ethnopharmacological information. Although each approach has its own strength in analyzing the characteristics of phytochemicals, previous studies have not considered them all together. Here, we apply an integrated in silico analysis to identify the potential health benefits of phytochemicals based on molecular analysis and chemical properties as well as ethnopharmacological evidence. From the molecular analysis, we found an average of 415.6 health effects for 591 phytochemicals. We further investigated ethnopharmacological evidence of phytochemicals and found that on average 129.1 (31%) of the predicted health effects had ethnopharmacological evidence. Lastly, we investigated chemical properties to confirm whether they are orally bio-available, drug available or effective on certain tissues. The evaluation results indicate that the health effects can be predicted more accurately by cooperatively considering the molecular analysis, chemical properties and ethnopharmacological evidence.

Regulation of ciliary membrane protein trafficking and signalling by kinesin motor proteins

Primary cilia are antenna-like sensory organelles that regulate a substantial number of cellular signalling pathways in vertebrates, both during embryonic development as well as in adulthood, and mutations in genes coding for ciliary proteins are causative of an expanding group of pleiotropic diseases known as ciliopathies. Cilia consist of a microtubule-based axoneme core, which is subtended by a basal body and covered by a bilayer lipid membrane of unique protein and lipid composition. Cilia are dynamic organelles, and the ability of cells to regulate ciliary protein and lipid content in response to specific cellular and environmental cues is crucial for balancing ciliary signalling output. Here we discuss mechanisms involved in regulation of ciliary membrane protein trafficking and signalling, with main focus on kinesin-2 and kinesin-3 family members.

Pooled genome-wide CRISPR screening for basal and context-specific fitness gene essentiality in Drosophila cells

Genome-wide screens in Drosophila cells have offered numerous insights into gene function, yet a major limitation has been the inability to stably deliver large multiplexed DNA libraries to cultured cells allowing barcoded pooled screens. Here, we developed a site-specific integration strategy for library delivery and performed a genome-wide CRISPR knockout screen in Drosophila S2R+ cells. Under basal growth conditions, 1235 genes were essential for cell fitness at a false-discovery rate of 5%, representing the highest-resolution fitness gene set yet assembled for Drosophila, including 407 genes which likely duplicated along the vertebrate lineage and whose orthologs were underrepresented in human CRISPR screens. We additionally performed context-specific fitness screens for resistance to or synergy with trametinib, a Ras/ERK/ETS inhibitor, or rapamycin, an mTOR inhibitor, and identified key regulators of each pathway. The results present a novel, scalable, and versatile platform for functional genomic screens in invertebrate cells.

IGF-Binding Proteins: Why Do They Exist and Why Are There So Many?

Insulin-like growth factors (IGFs) are key growth-promoting peptides that act as both endocrine hormones and autocrine/paracrine growth factors. In the bloodstream and in local tissues, most IGF molecules are bound by one of the members of the IGF-binding protein (IGFBP) family, of which six distinct types exist. These proteins bind to IGF with an equal or greater affinity than the IGF1 receptor and are thus in a key position to regulate IGF signaling globally and locally. Binding to an IGFBP increases the half-life of IGF in the circulation and blocks its potential binding to the insulin receptor. In addition to these classical roles, IGFBPs have been shown to modulate IGF signaling locally under various conditions. Although members of the IGFBP family share significant sequence homology, they each have unique structural features and play distinct roles. These IGFBP genes also have different modes of regulation and distinct expression patterns. Some IGFBPs have been found to bind to their own receptors or to translocate into the interior compartments of cells where they may execute IGF-independent actions. In spite of this functional and regulatory diversity, it has been puzzling that loss-of-function studies have yielded relatively little information about the physiological functions of IGFBPs. In this review, we suggest that evolution has tended to retain an array of IGFBPs in order to facilitate fine-tuning of IGF signaling. We explore the emerging explanation that many IGFBP functions have evolved to allow the targeted adjustment of IGF signaling under stressful or irregular conditions, which would likely not be revealed in a standard laboratory setting.

Engineering the Drosophila Genome for Developmental Biology

The recent development of transposon and CRISPR-Cas9-based tools for manipulating the fly genome in vivo promises tremendous progress in our ability to study developmental processes. Tools for introducing tags into genes at their endogenous genomic loci facilitate imaging or biochemistry approaches at the cellular or subcellular levels. Similarly, the ability to make specific alterations to the genome sequence allows much more precise genetic control to address questions of gene function.