Genetic dissection of developmental pathways

D-SPIN constructs gene regulatory network models from multiplexed scRNA-seq data revealing organizing principles of cellular perturbation response

Gene regulatory networks within cells modulate the expression of the genome in response to signals and changing environmental conditions. Reconstructions of gene regulatory networks can reveal the information processing and control principles used by cells to maintain homeostasis and execute cell-state transitions. Here, we introduce a computational framework, D-SPIN, that generates quantitative models of gene regulatory networks from single-cell mRNA-seq datasets collected across thousands of distinct perturbation conditions. D-SPIN models the cell as a collection of interacting gene-expression programs, and constructs a probabilistic model to infer regulatory interactions between gene-expression programs and external perturbations. Using large Perturb-seq and drug-response datasets, we demonstrate that D-SPIN models reveal the organization of cellular pathways, sub-functions of macromolecular complexes, and the logic of cellular regulation of transcription, translation, metabolism, and protein degradation in response to gene knockdown perturbations. D-SPIN can also be applied to dissect drug response mechanisms in heterogeneous cell populations, elucidating how combinations of immunomodulatory drugs can induce novel cell states through additive recruitment of gene expression programs. D-SPIN provides a computational framework for constructing interpretable models of gene-regulatory networks to reveal principles of cellular information processing and physiological control.

Innexin function dictates the spatial relationship between distal somatic cells in the Caenorhabditis elegans gonad without impacting the germline stem cell pool

Gap-junctional signaling mediates myriad cellular interactions in metazoans. Yet, how gap junctions control the positioning of cells in organs is not well understood. Innexins compose gap junctions in invertebrates and affect organ architecture. Here, we investigate the roles of gap-junctions in controlling distal somatic gonad architecture and its relationship to underlying germline stem cells in Caenorhabditis elegans. We show that a reduction of soma-germline gap-junctional activity causes displacement of distal sheath cells (Sh1) towards the distal end of the gonad. We confirm, by live imaging, transmission electron microscopy, and antibody staining, that bare regions-lacking somatic gonadal cell coverage of germ cells-are present between the distal tip cell (DTC) and Sh1, and we show that an innexin fusion protein used in a prior study encodes an antimorphic gap junction subunit that mispositions Sh1. We determine that, contrary to the model put forth in the prior study based on this fusion protein, Sh1 mispositioning does not markedly alter the position of the borders of the stem cell pool nor of the progenitor cell pool. Together, these results demonstrate that gap junctions can control the position of Sh1, but that Sh1 position is neither relevant for GLP-1/Notch signaling nor for the exit of germ cells from the stem cell pool.

Harnessing the power of genetics: fast forward genetics in Caenorhabditis elegans

Forward genetics is a powerful tool to unravel molecular mechanisms of diverse biological processes. The success of genetic screens primarily relies on the ease of genetic manipulation of an organism and the availability of a plethora of genetic tools. The roundworm Caenorhabditis elegans has been one of the favorite models for genetic studies due to its hermaphroditic lifestyle, ease of maintenance, and availability of various genetic manipulation tools. The strength of C. elegans genetics is highlighted by the leading role of this organism in the discovery of several conserved biological processes. In this review, the principles and strategies for forward genetics in C. elegans are discussed. Further, the recent advancements that have drastically accelerated the otherwise time-consuming process of mutation identification, making forward genetic screens a method of choice for understanding biological functions, are discussed. The emphasis of the review has been on providing practical and conceptual pointers for designing genetic screens that will identify mutations, specifically disrupting the biological processes of interest.

Gap Junctions and NCA Cation Channels Are Critical for Developmentally Timed Sleep and Arousal in Caenorhabditis elegans

An essential characteristic of sleep is heightened arousal threshold, with decreased behavioral response to external stimuli. The molecular and cellular mechanisms underlying arousal threshold changes during sleep are not fully understood. We report that loss of UNC-7 or UNC-9 innexin function dramatically reduced sleep and decreased arousal threshold during developmentally timed sleep in Caenorhabditis elegans UNC-7 function was required in premotor interneurons and UNC-9 function was required in motor neurons in this paradigm. Simultaneous transient overexpression of UNC-7 and UNC-9 was sufficient to induce anachronistic sleep in adult animals. Moreover, loss of UNC-7 or UNC-9 suppressed the increased sleep of EGL-4 gain-of-function animals, which have increased cyclic-GMP-dependent protein kinase activity. These results suggest C. elegans gap junctions may act downstream of previously identified sleep regulators. In other paradigms, the NCA cation channels act upstream of gap junctions. Consistent with this, diminished NCA channel activity in C. elegans robustly increased arousal thresholds during sleep bouts in L4-to-adult developmentally timed sleep. Total time in sleep bouts was only modestly increased in animals lacking NCA channel auxiliary subunit UNC-79, whereas increased channel activity dramatically decreased sleep. Loss of EGL-4 or innexin proteins suppressed UNC-79 loss-of-function sleep and arousal defects. In Drosophila, the ion channel narrow abdomen, an ortholog of the C. elegans NCA channels, drive the pigment dispersing factor (PDF) neuropeptide release, regulating circadian behavior. However, in C. elegans, we found that loss of the PDF receptor PDFR-1 did not suppress gain-of-function sleep defects, suggesting an alternative downstream pathway. This study emphasizes the conservation and importance of neuronal activity modulation during sleep, and unequivocally demonstrates that gap junction function is critical for normal sleep.

RNA Sequencing Analysis of the msl2msl3, crl, and ggps1 Mutants Indicates that Diverse Sources of Plastid Dysfunction Do Not Alter Leaf Morphology Through a Common Signaling Pathway

Determining whether individual genes function in the same or in different pathways is an important aspect of genetic analysis. As an alternative to the construction of higher-order mutants, we used contemporary expression profiling methods to perform pathway analysis on several Arabidopsis thaliana mutants, including the mscS-like (msl)2msl3 double mutant. MSL2 and MSL3 are implicated in plastid ion homeostasis, and msl2msl3 double mutants exhibit leaves with a lobed periphery, a rumpled surface, and disturbed mesophyll cell organization. Similar developmental phenotypes are also observed in other mutants with defects in a range of other chloroplast or mitochondrial functions, including biogenesis, gene expression, and metabolism. We wished to test the hypothesis that the common leaf morphology phenotypes of these mutants are the result of a characteristic nuclear expression pattern that is generated in response to organelle dysfunction. RNA-Sequencing was performed on aerial tissue of msl2msl3 geranylgeranyl diphosphate synthase 1 (ggps1), and crumpled leaf (crl) mutants. While large groups of co-expressed genes were identified in pairwise comparisons between genotypes, we were only able to identify a small set of genes that showed similar expression profiles in all three genotypes. Subsequent comparison to the previously published gene expression profiles of two other mutants, yellow variegated 2 (var2) and scabra3 (sca3), failed to reveal a common pattern of gene expression associated with superficially similar leaf morphology defects. Nor did we observe overlap between genes differentially expressed in msl2msl3, crl, and ggps1 and a previously identified retrograde core response module. These data suggest that a common retrograde signaling pathway initiated by organelle dysfunction either does not exist in these mutants or cannot be identified through transcriptomic methods. Instead, the leaf developmental defects observed in these mutants may be achieved through a number of independent pathways.

Gene network inference by probabilistic scoring of relationships from a factorized model of interactions

Motivation: Epistasis analysis is an essential tool of classical genetics for inferring the order of function of genes in a common pathway. Typically, it considers single and double mutant phenotypes and for a pair of genes observes whether a change in the first gene masks the effects of the mutation in the second gene. Despite the recent emergence of biotechnology techniques that can provide gene interaction data on a large, possibly genomic scale, few methods are available for quantitative epistasis analysis and epistasis-based network reconstruction.

Results: We here propose a conceptually new probabilistic approach to gene network inference from quantitative interaction data. The approach is founded on epistasis analysis. Its features are joint treatment of the mutant phenotype data with a factorized model and probabilistic scoring of pairwise gene relationships that are inferred from the latent gene representation. The resulting gene network is assembled from scored pairwise relationships. In an experimental study, we show that the proposed approach can accurately reconstruct several known pathways and that it surpasses the accuracy of current approaches.

Availability and implementation: Source code is available at http://github.com/biolab/red.

Contact:blaz.zupan@fri.uni-lj.si

Supplementary information:Supplementary data are available at Bioinformatics online.

From genes to drugs: target validation in Caenorhabditis elegans

A central challenge in post-genomic drug discovery is selection of relevant therapeutic targets from a large pool of candidates, so that resources are invested productively. Key to meeting this challenge is demonstrating target function in disease-related pathways in vivo. The nematode Caenorhabditis elegans, a well-characterized experimental system for genetic analysis of biological regulatory pathways, provides a powerful means of assessing the impact of modulating target function on biological processes, thus facilitating selection of high-value targets.:

A statistical procedure to map high-order epistasis for complex traits

Genetic interactions or epistasis have been thought to play a pivotal role in shaping the formation, development and evolution of life. Previous work focused on lower-order interactions between a pair of genes, but it is obviously inadequate to explain a complex network of genetic interactions and pathways. We review and assess a statistical model for characterizing high-order epistasis among more than two genes or quantitative trait loci (QTLs) that control a complex trait. The model includes a series of start-of-the-art standard procedures for estimating and testing the nature and magnitude of QTL interactions. Results from simulation studies and real data analysis warrant the statistical properties of the model and its usefulness in practice. High-order epistatic mapping will provide a routine procedure for charting a detailed picture of the genetic regulation mechanisms underlying the phenotypic variation of complex traits.

Analysis of aging in Caenorhabditis elegans

This chapter is dedicated to the study of aging in Caenorhabditis elegans (C. elegans). The assays are divided into two sections. In the first section, we describe detailed protocols for performing life span analysis in solid and liquid medium. In the second section, we describe various assays for measuring age-related changes. Our laboratory has been involved in several fruitful collaborations with non-C. elegans researchers keen on testing a role for their favorite gene in modulating aging (Carrano et al., 2009; Dong et al., 2007; Raices et al., 2008; Wolff et al., 2006). But even with the guidance of trained worm biologists, this undertaking can be daunting. We hope that this chapter will serve as a worthy compendium for those researchers who may or may not have immediate access to laboratories studying C. elegans.

Genetic analysis of IP3 and calcium signalling pathways in C. elegans

BACKGROUND

The nematode, Caenorhabditis elegans is an established model system that is particularly well suited to genetic analysis. C. elegans is easily manipulated and we have an in depth knowledge of many aspects of its biology. Thus, it is an attractive system in which to pursue integrated studies of signalling pathways. C. elegans has a complement of calcium signalling molecules similar to that of other animals.

SCOPE OF REVIEW

We focus on IP3 signalling. We describe how forward and reverse genetic approaches, including RNAi, have resulted in a tool kit which enables the analysis of IP3/Ca2+ signalling pathways. The importance of cell and tissue specific manipulation of signalling pathways and the use of epistasis analysis are highlighted. We discuss how these tools have increased our understanding of IP3 signalling in specific developmental, physiological and behavioural roles. Approaches to imaging calcium signals in C. elegans are considered.

MAJOR CONCLUSIONS

A wide selection of tools is available for the analysis of IP3/Ca2+ signalling in C. elegans. This has resulted in detailed descriptions of the function of IP3/Ca2+ signalling in the animal's biology. Nevertheless many questions about how IP3 signalling regulates specific processes remain.

GENERAL SIGNIFICANCE

Many of the approaches described may be applied to other calcium signalling systems. C. elegans offers the opportunity to dissect pathways, perform integrated studies and to test the importance of the properties of calcium signalling molecules to whole animal function, thus illuminating the function of calcium signalling in animals. This article is part of a Special Issue entitled Biochemical, biophysical and genetic approaches to intracellular calcium signalling.

Dissection of genetic pathways in C. elegans

With unique genetic and cell biological strengths, C. elegans has emerged as a powerful model system for studying many biological processes. These processes are typically regulated by complex genetic networks consisting of genes. Identifying those genes and organizing them into genetic pathways are two major steps toward understanding the mechanisms that regulate biological events. Forward genetic screens with various designs are a traditional approach for identifying candidate genes. The completion of the genome sequencing in C. elegans and the advent of high-throughput experimental techniques have led to the development of two additional powerful approaches: functional genomics and systems biology. Genes that are discovered by these approaches can be ordered into interacting pathways through a variety of strategies, involving genetics, cell biology, biochemistry, and functional genomics, to gain a more complete understanding of how gene regulatory networks control a particular biological process. The aim of this review is to provide an overview of the approaches available to identify and construct the genetic pathways using C. elegans.

A general model for multilocus epistatic interactions in case-control studies

Background

Epistasis, i.e., the interaction of alleles at different loci, is thought to play a central role in the formation and progression of complex diseases. The complexity of disease expression should arise from a complex network of epistatic interactions involving multiple genes.

Methodology

We develop a general model for testing high-order epistatic interactions for a complex disease in a case-control study. We incorporate the quantitative genetic theory of high-order epistasis into the setting of cases and controls sampled from a natural population. The new model allows the identification and testing of epistasis and its various genetic components.

Conclusions

Simulation studies were used to examine the power and false positive rates of the model under different sampling strategies. The model was used to detect epistasis in a case-control study of inflammatory bowel disease, in which five SNPs at a candidate gene were typed, leading to the identification of a significant three-locus epistasis.

Insights into the molecular determinants of proton inhibition in an acid-inactivated degenerins and mammalian epithelial Na(+) channel

Mammalian ASIC channels of the DEG/ENaC superfamily are gated by extracellular protons and function to mediate touch and pain sensitivity, learning and memory, and fear conditioning. The recently solved crystal structure of chicken ASIC1a and preliminary functional studies suggested that a highly negatively charged pocket in the extracellular domain of the channel might be the primary proton binding domain. However, more recent extensive mutagenesis analysis paints a more complex mechanism of channel gating, involving binding of protons at sites immediately after the first transmembrane domain (TM1) and displacement of inhibitory Ca(2+) ions from the acidic pocket in the extracellular domain and from another Ca(2+) binding site at the mouth of the pore. We recently identified and functionally characterized Caenorhabditis elegans ACD-1, the first acid-inactivated DEG/ENaC channel. ACD-1 is expressed in C. elegans amphid glia and functions with neuronal DEG/ENaC channel DEG-1 to mediate acid avoidance and chemotaxis to the amino acid lysine. The post-TM1 residues that were proposed to bind protons in ASIC1a are not conserved in ACD-1, but some of the amino acids constituting the acidic pocket are. However, ACD-1 proton sensitivity is completely independent from extracellular Ca(2+), and protons appear to bind the channel in a less cooperative manner. We thus wondered if residues in the acidic pocket might contribute to ACD-1 acid sensitivity. We show here that while ACD-1 sensitivity to extracellular protons is influenced by mutations in the acidic pocket, other sites are likely to participate. We also report that one histidine at the base of the thumb and residues in the channel pore influence proton inhibition in a voltage-independent manner, suggesting that they affect the coupling of proton binding with the gating rather than proton binding itself. We conclude that ACD-1 inhibition by protons is likely mediated by binding of proton ions to multiple sites throughout the extracellular domain of the channel. Our data also support a model in which residues in the acidic pocket contribute to determining the channel state perhaps by changing the strength of the interaction between adjacent thumb and finger domains.

Characterizing the transcriptional regulation of let-721, a Caenorhabditis elegans homolog of human electron flavoprotein dehydrogenase

LET-721 is the Caenorhabditis elegans ortholog of electron-transferring flavoprotein dehydrogenase (ETFDH). We are studying this protein in C. elegans in order to establish a tractable model system for further exploration of ETFDH structure and function. ETFDH is an inner mitochondrial membrane localized enzyme that plays a key role in the beta-oxidation of fatty acids and catabolism of amino acids and choline. ETFDH accepts electrons from at least twelve mitochondrial matrix flavoprotein dehydrogenases via an intermediate dimer protein and transfers the electrons to ubiquinone. In humans, ETFDH mutations result in the autosomal recessive metabolic disorder, multiple acyl-CoA dehydrogenase deficiency. Mutants of let-721 in C. elegans are either maternal effect lethals or semi-sterile. let-721 is transcribed in the pharynx, body wall muscle, hypoderm, intestine and somatic gonad. In addition, the subcellular localization of LET-721 agrees with predictions that it is localized to mitochondria. We identified and confirmed three cis-regulatory sequences (pha-site, rep-site, and act-site). Phylogenetic footprinting of each site indicates that they are conserved between four Caenorhabditis species. The pha-site mapped roughly 1,300 bp upstream of let-721's translational start site and is necessary for expression in pharyngeal tissues. The rep-site mapped roughly 830 bp upstream of the translational start site and represses expression of LET-721 within pharyngeal tissues. The act-site mapped roughly 800 bp upstream of the translational start site and is required for expression within spermatheca, body wall muscle, pharynx, and intestine. Taken together, we find that LET-721 is a mitochondrially expressed protein that is under complex transcriptional controls.

Raw mediates antagonism of AP-1 activity in Drosophila

High baselines of transcription factor activities represent fundamental obstacles to regulated signaling. Here we show that in Drosophila, quenching of basal activator protein 1 (AP-1) transcription factor activity serves as a prerequisite to its tight spatial and temporal control by the JNK (Jun N-terminal kinase) signaling cascade. Our studies indicate that the novel raw gene product is required to limit AP-1 activity to leading edge epidermal cells during embryonic dorsal closure. In addition, we provide the first evidence that the epidermis has a Basket JNK-independent capacity to activate AP-1 targets and that raw function is required broadly throughout the epidermis to antagonize this activity. Finally, our mechanistic studies of the three dorsal-open group genes [raw, ribbon (rib), and puckered (puc)] indicate that these gene products provide at least two tiers of JNK/AP-1 regulation. In addition to Puckered phosphatase function in leading edge epidermal cells as a negative-feedback regulator of JNK signaling, the three dorsal-open group gene products (Raw, Ribbon, and Puckered) are required more broadly in the dorsolateral epidermis to quench a basal, signaling-independent activity of the AP-1 transcription factor.

A novel gain-of-function mutant of the cyclic GMP-dependent protein kinase egl-4 affects multiple physiological processes in Caenorhabditis elegans

cGMP-dependent protein kinases are key intracellular transducers of cell signaling. We identified a novel dominant mutation in the C. elegans egl-4 cGMP-dependent protein kinase (PKG) and show that this mutation causes increased normal gene activity although it is associated with a reduced EGL-4 protein level. Prior phenotypic analyses of this gain-of-function mutant demonstrated a reduced longevity and a reduced feeding behavior when the animals were left unperturbed. We characterize several additional phenotypes caused by increased gene activity of egl-4. These phenotypes include a small body size, reduced locomotion in the presence of food, a pale intestine, increased intestinal fat storage, and a decreased propensity to form dauer larvae. The multiple phenotypes of egl-4 dominant mutants are consistent with an instructive signaling role of PKG to control many aspects of animal physiology. This is among the first reported gain-of-function mutations in this enzyme of central physiological importance. In a genetic screen we have identified extragenic suppressors of this gain-of-function mutant. Thus, this mutant promises to be a useful tool for identifying downstream targets of PKG.

RNA interference of peroxisome-related genes in C. elegans: a new model for human peroxisomal disorders

RNA-mediated interference (RNAi) for the posttranscriptional silencing of genes was used to evaluate the importance of various peroxisomal enzymes and peroxins for the development of Caenorhabditis elegans and to compare the roles of these proteins in the nematode to their roles in yeasts and humans. The nematode counterparts of the human ATP-binding cassette half-transporters, the enzymes alkyldihydroxyacetonephosphate synthase and Delta(3,5)-Delta (2,4)-dienoyl-CoA isomerase, the receptors for peroxisomal membrane and matrix proteins (Pex19p and Pex5p), and components of the docking and translocation machineries for matrix proteins (Pex13p and Pex12p) are essential for the development of C. elegans. Unexpectedly, RNAi silencing of the acyl-CoA synthetase-mediated activation of fatty acids, the alpha- and beta-oxidation of fatty acids, the intraperoxisomal decomposition of hydrogen peroxide, and the peroxins Pex1p, Pex2p, and Pex6p had no apparent effect on C. elegans development. The described analysis of functional gene knockouts through RNAi provides a basis for the use of C. elegans as a valuable model system with which to study the molecular and physiological defects underlying the human peroxisomal disorders.

Goalpha regulates volatile anesthetic action in Caenorhabditis elegans

To identify genes controlling volatile anesthetic (VA) action, we have screened through existing Caenorhabditis elegans mutants and found that strains with a reduction in Go signaling are VA resistant. Loss-of-function mutants of the gene goa-1, which codes for the alpha-subunit of Go, have EC(50)s for the VA isoflurane of 1.7- to 2.4-fold that of wild type. Strains overexpressing egl-10, which codes for an RGS protein negatively regulating goa-1, are also isoflurane resistant. However, sensitivity to halothane, a structurally distinct VA, is differentially affected by Go pathway mutants. The RGS overexpressing strains, a goa-1 missense mutant found to carry a novel mutation near the GTP-binding domain, and eat-16(rf) mutants, which suppress goa-1(gf) mutations, are all halothane resistant; goa-1(null) mutants have wild-type sensitivities. Double mutant strains carrying mutations in both goa-1 and unc-64, which codes for a neuronal syntaxin previously found to regulate VA sensitivity, show that the syntaxin mutant phenotypes depend in part on goa-1 expression. Pharmacological assays using the cholinesterase inhibitor aldicarb suggest that VAs and GOA-1 similarly downregulate cholinergic neurotransmitter release in C. elegans. Thus, the mechanism of action of VAs in C. elegans is regulated by Goalpha, and presynaptic Goalpha-effectors are candidate VA molecular targets.

Molecular signals versus the Loi de Balancement

Life history tradeoffs are often thought to be caused by the allocation of limited resources among competing traits such as reproduction, somatic growth and maintenance. One line of evidence supporting this comes from eliminating reproduction, for example, by surgically removing gonads. However, recent evidence from the nematode Caenorhabditis elegans suggests that the apparent tradeoffs it shows might not be due to resource allocation at all but rather to the effects of a molecular signal originating in the germ line that represses longevity. These results should cause us to rethink the interpretation of many classic experiments in life history evolution.

Suppression mechanisms: themes from variations

Suppressor analysis is a commonly used strategy to identify functional relationships between genes that might not have been revealed through other genetic or biochemical means. Many mechanisms that explain the phenomenon of genetic suppression have been described, but the wide variety of possible mechanisms can present a challenge to defining the relationship between a suppressor and the original gene. This article provides a broad framework for classifying suppression mechanisms and describes a series of genetic tests that can be applied to determine the most likely mechanism of suppression.

Functional genomics in the mouse: phenotype-based mutagenesis screens

Significant progress has been made in sequencing the genomes of several model organisms, and efforts are now underway to complete the sequencing of the human genome. In parallel with this effort, new approaches are being developed for the elucidation of the functional content of the human genome. The mouse will have an important role in this phase of the genome project as a model system. In this review we discuss and compare classical genetic approaches to gene function-phenotype-based mutagenesis screens aimed at the establishment of a large collection of single gene mutations affecting a wide range of phenotypic traits in the mouse. Whereas large scale genome-wide screens that are directed at the identification of all loci contributing to a specific phenotype may be impractical, region-specific saturation screens that provide mutations within a delimited chromosomal region are a feasible alternative. Region-specific screens in the mouse can be performed in only two generations by combining high-efficiency chemical mutagenesis with deletion complexes generated using embryonic stem (ES) cells. The ability to create and analyze deletion complexes rapidly, as well as to map novel chemically-induced mutations within these complexes, will facilitate systematic functional analysis of the mouse genome and corresponding gene sequences in humans. Furthermore, as the extent of the mouse genome sequencing effort is still uncertain, we underscore a necessity to direct sequencing efforts to those chromosomal regions that are targets for extensive mutagenesis screens.