Natural Variation in Gene Expression Modulates the Severity of Mutant Phenotypes.

Many mutations cause genetic disorders. However, two people inheriting the same mutation often have different severity of symptoms, and this is partly genetic. The effects of genetic background on mutant phenotypes are poorly understood, but predicting them is critical for personalized medicine. To study this phenomenon comprehensively and systematically, we used RNAi to compare loss-of-function phenotypes for ∼1,400 genes in two isolates of C. elegans and find that ∼20% of genes differ in the severity of phenotypes in these two genetic backgrounds. Crucially, this effect of genetic background on the severity of both RNAi and mutant phenotypes can be predicted from variation in the expression levels of the affected gene. This is also true in mammalian cells, suggesting it is a general property of genetic networks. We suggest that differences in the manifestation of mutant phenotypes between individuals are largely the result of natural variation in gene expression.

The majority of animal genes are required for wild-type fitness.

Almost all eukaryotic genes are conserved, suggesting that they have essential functions. However, only a minority of genes have detectable loss-of-function phenotypes in experimental assays, and multiple theories have been proposed to explain this discrepancy. Here, we use RNA-mediated interference in C. elegans to examine how knockdown of any gene affects the overall fitness of worm populations. Whereas previous studies typically assess phenotypes that are detectable by eye after a single generation, we monitored growth quantitatively over several generations. In contrast to previous estimates, we find that, in these multigeneration population assays, the majority of genes affect fitness, and this suggests that genetic networks are not robust to mutation. Our results demonstrate that, in a single environmental condition, most animal genes play essential roles. This is a higher proportion than for yeast genes, and we suggest that the source of negative selection is different in animals and in unicellular eukaryotes.

EPIC: software toolkit for elution profile-based inference of protein complexes.

Protein complexes are key macromolecular machines of the cell, but their description remains incomplete. We and others previously reported an experimental strategy for global characterization of native protein assemblies based on chromatographic fractionation of biological extracts coupled to precision mass spectrometry analysis (chromatographic fractionation-mass spectrometry, CF-MS), but the resulting data are challenging to process and interpret. Here, we describe EPIC (elution profile-based inference of complexes), a software toolkit for automated scoring of large-scale CF-MS data to define high-confidence multi-component macromolecules from diverse biological specimens. As a case study, we used EPIC to map the global interactome of Caenorhabditis elegans, defining 612 putative worm protein complexes linked to diverse biological processes. These included novel subunits and assemblies unique to nematodes that we validated using orthogonal methods. The open source EPIC software is freely available as a Jupyter notebook packaged in a Docker container (https://hub.docker.com/r/baderlab/bio-epic/).

A pair of RNA-binding proteins controls networks of splicing events contributing to specialization of neural cell types.

Alternative splicing is important for the development and function of the nervous system, but little is known about the differences in alternative splicing between distinct types of neurons. Furthermore, the factors that control cell-type-specific splicing and the physiological roles of these alternative isoforms are unclear. By monitoring alternative splicing at single-cell resolution in Caenorhabditis elegans, we demonstrate that splicing patterns in different neurons are often distinct and highly regulated. We identify two conserved RNA-binding proteins, UNC-75/CELF and EXC-7/Hu/ELAV, which regulate overlapping networks of splicing events in GABAergic and cholinergic neurons. We use the UNC-75 exon network to discover regulators of synaptic transmission and to identify unique roles for isoforms of UNC-64/Syntaxin, a protein required for synaptic vesicle fusion. Our results indicate that combinatorial regulation of alternative splicing in distinct neurons provides a mechanism to specialize metazoan nervous systems.

The combinatorial control of alternative splicing in C. elegans.

Normal development requires the right splice variants to be made in the right tissues at the right time. The core splicing machinery is engaged in all splicing events, but which precise splice variant is made requires the choice between alternative splice sites-for this to occur, a set of splicing factors (SFs) must recognize and bind to short RNA motifs in the pre-mRNA. In C. elegans, there is known to be extensive variation in splicing patterns across development, but little is known about the targets of each SF or how multiple SFs combine to regulate splicing. Here we combine RNA-seq with in vitro binding assays to study how 4 different C. elegans SFs, ASD-1, FOX-1, MEC-8, and EXC-7, regulate splicing. The 4 SFs chosen all have well-characterised biology and well-studied loss-of-function genetic alleles, and all contain RRM domains. Intriguingly, while the SFs we examined have varied roles in C. elegans development, they show an unexpectedly high overlap in their targets. We also find that binding sites for these SFs occur on the same pre-mRNAs more frequently than expected suggesting extensive combinatorial control of splicing. We confirm that regulation of splicing by multiple SFs is often combinatorial and show that this is functionally significant. We also find that SFs appear to combine to affect splicing in two modes-they either bind in close proximity within the same intron or they appear to bind to separate regions of the intron in a conserved order. Finally, we find that the genes whose splicing are regulated by multiple SFs are highly enriched for genes involved in the cytoskeleton and in ion channels that are key for neurotransmission. Together, this shows that specific classes of genes have complex combinatorial regulation of splicing and that this combinatorial regulation is critical for normal development to occur.

A compendium of RNA-binding motifs for decoding gene regulation.

RNA-binding proteins are key regulators of gene expression, yet only a small fraction have been functionally characterized. Here we report a systematic analysis of the RNA motifs recognized by RNA-binding proteins, encompassing 205 distinct genes from 24 diverse eukaryotes. The sequence specificities of RNA-binding proteins display deep evolutionary conservation, and the recognition preferences for a large fraction of metazoan RNA-binding proteins can thus be inferred from their RNA-binding domain sequence. The motifs that we identify in vitro correlate well with in vivo RNA-binding data. Moreover, we can associate them with distinct functional roles in diverse types of post-transcriptional regulation, enabling new insights into the functions of RNA-binding proteins both in normal physiology and in human disease. These data provide an unprecedented overview of RNA-binding proteins and their targets, and constitute an invaluable resource for determining post-transcriptional regulatory mechanisms in eukaryotes.

C. elegans SUP-46, an HNRNPM family RNA-binding protein that prevents paternally-mediated epigenetic sterility.

BACKGROUND: In addition to DNA, gametes contribute epigenetic information in the form of histones and non-coding RNA. Epigenetic programs often respond to stressful environmental conditions and provide a heritable history of ancestral stress that allows for adaptation and propagation of the species. In the nematode C. elegans, defective epigenetic transmission often manifests as progressive germline mortality. We previously isolated sup-46 in a screen for suppressors of the hexosamine pathway gene mutant, gna-2(qa705). In this study, we examine the role of SUP-46 in stress resistance and progressive germline mortality. RESULTS: We identified SUP-46 as an HNRNPM family RNA-binding protein, and uncovered a highly novel role for SUP-46 in preventing paternally-mediated progressive germline mortality following mating. Proximity biotinylation profiling of human homologs (HNRNPM, MYEF2) identified proteins of ribonucleoprotein complexes previously shown to contain non-coding RNA. Like HNRNPM and MYEF2, SUP-46 was associated with multiple RNA granules, including stress granules, and also formed granules on active chromatin. SUP-46 depletion disrupted germ RNA granules and caused ectopic sperm, increased sperm transcripts, and chronic heat stress sensitivity. SUP-46 was also required for resistance to acute heat stress, and a conserved "MYEF2" motif was identified that was needed for stress resistance. CONCLUSIONS: In mammals, non-coding RNA from the sperm of stressed males has been shown to recapitulate paternal stress phenotypes in the offspring. Our results suggest that HNRNPM family proteins enable stress resistance and paternally-mediated epigenetic transmission that may be conserved across species.

Heritability and genetic basis of protein level variation in an outbred population.

The genetic basis of heritable traits has been studied for decades. Although recent mapping efforts have elucidated genetic determinants of transcript levels, mapping of protein abundance has lagged. Here, we analyze levels of 4084 GFP-tagged yeast proteins in the progeny of a cross between a laboratory and a wild strain using flow cytometry and high-content microscopy. The genotype of trans variants contributed little to protein level variation between individual cells but explained >50% of the variance in the population's average protein abundance for half of the GFP fusions tested. To map trans-acting factors responsible, we performed flow sorting and bulk segregant analysis of 25 proteins, finding a median of five protein quantitative trait loci (pQTLs) per GFP fusion. Further, we find that cis-acting variants predominate; the genotype of a gene and its surrounding region had a large effect on protein level six times more frequently than the rest of the genome combined. We present evidence for both shared and independent genetic control of transcript and protein abundance: More than half of the expression QTLs (eQTLs) contribute to changes in protein levels of regulated genes, but several pQTLs do not affect their cognate transcript levels. Allele replacements of genes known to underlie trans eQTL hotspots confirmed the correlation of effects on mRNA and protein levels. This study represents the first genome-scale measurement of genetic contribution to protein levels in single cells and populations, identifies more than a hundred trans pQTLs, and validates the propagation of effects associated with transcript variation to protein abundance.

Comparative RNAi screens in C. elegans and C. briggsae reveal the impact of developmental system drift on gene function.

Although two related species may have extremely similar phenotypes, the genetic networks underpinning this conserved biology may have diverged substantially since they last shared a common ancestor. This is termed Developmental System Drift (DSD) and reflects the plasticity of genetic networks. One consequence of DSD is that some orthologous genes will have evolved different in vivo functions in two such phenotypically similar, related species and will therefore have different loss of function phenotypes. Here we report an RNAi screen in C. elegans and C. briggsae to identify such cases. We screened 1333 genes in both species and identified 91 orthologues that have different RNAi phenotypes. Intriguingly, we find that recently evolved genes of unknown function have the fastest evolving in vivo functions and, in several cases, we identify the molecular events driving these changes. We thus find that DSD has a major impact on the evolution of gene function and we anticipate that the C. briggsae RNAi library reported here will drive future studies on comparative functional genomics screens in these nematodes.

Nuclear receptor binding protein 1 regulates intestinal progenitor cell homeostasis and tumour formation.

Genetic screens in simple model organisms have identified many of the key components of the conserved signal transduction pathways that are oncogenic when misregulated. Here, we identify H37N21.1 as a gene that regulates vulval induction in let-60(n1046gf), a strain with a gain-of-function mutation in the Caenorhabditis elegans Ras orthologue, and show that somatic deletion of Nrbp1, the mouse orthologue of this gene, results in an intestinal progenitor cell phenotype that leads to profound changes in the proliferation and differentiation of all intestinal cell lineages. We show that Nrbp1 interacts with key components of the ubiquitination machinery and that loss of Nrbp1 in the intestine results in the accumulation of Sall4, a key mediator of stem cell fate, and of Tsc22d2. We also reveal that somatic loss of Nrbp1 results in tumourigenesis, with haematological and intestinal tumours predominating, and that nuclear receptor binding protein 1 (NRBP1) is downregulated in a range of human tumours, where low expression correlates with a poor prognosis. Thus NRBP1 is a conserved regulator of cell fate, that plays an important role in tumour suppression.

Systematic mapping of genetic interactions in Caenorhabditis elegans identifies common modifiers of diverse signaling pathways.

Most heritable traits, including disease susceptibility, are affected by interactions between multiple genes. However, we understand little about how genes interact because very few possible genetic interactions have been explored experimentally. We have used RNA interference in Caenorhabditis elegans to systematically test approximately 65,000 pairs of genes for their ability to interact genetically. We identify approximately 350 genetic interactions between genes functioning in signaling pathways that are mutated in human diseases, including components of the EGF/Ras, Notch and Wnt pathways. Most notably, we identify a class of highly connected 'hub' genes: inactivation of these genes can enhance the phenotypic consequences of mutation of many different genes. These hub genes all encode chromatin regulators, and their activity as genetic hubs seems to be conserved across animals. We propose that these genes function as general buffers of genetic variation and that these hub genes may act as modifier genes in multiple, mechanistically unrelated genetic diseases in humans.

Global impact of RNA polymerase II elongation inhibition on alternative splicing regulation.

The rate of RNA polymerase II (Pol II) elongation can influence splice site selection in nascent transcripts, yet the extent and physiological relevance of this kinetic coupling between transcription and alternative splicing (AS) is not well understood. We performed experiments to perturb Pol II elongation and then globally compared AS patterns with genome-wide Pol II occupancy. RNA binding and RNA processing functions were significantly enriched among the genes with Pol II elongation inhibition-dependent changes in AS. Under conditions that interfere with Pol II elongation, including cell stress, increased Pol II occupancy was detected in the intronic regions flanking the alternative exons in these genes, and these exons generally became more included. A disproportionately high fraction of these exons introduced premature termination codons that elicited nonsense-mediated mRNA decay (NMD), thereby further reducing transcript levels. Our results provide evidence that kinetic coupling between transcription, AS, and NMD affords a rapid mechanism by which cells can respond to changes in growth conditions, including cell stress, to coordinate the levels of RNA processing factors with mRNA levels.

Genome-wide analysis of alternative splicing in Caenorhabditis elegans.

Alternative splicing (AS) plays a crucial role in the diversification of gene function and regulation. Consequently, the systematic identification and characterization of temporally regulated splice variants is of critical importance to understanding animal development. We have used high-throughput RNA sequencing and microarray profiling to analyze AS in C. elegans across various stages of development. This analysis identified thousands of novel splicing events, including hundreds of developmentally regulated AS events. To make these data easily accessible and informative, we constructed the C. elegans Splice Browser, a web resource in which researchers can mine AS events of interest and retrieve information about their relative levels and regulation across development. The data presented in this study, along with the Splice Browser, provide the most comprehensive set of annotated splice variants in C. elegans to date, and are therefore expected to facilitate focused, high resolution in vivo functional assays of AS function.

Identification of a family of species-selective complex I inhibitors as potential anthelmintics.

Soil-transmitted helminths (STHs) are major pathogens infecting over a billion people. There are few classes of anthelmintics and there is an urgent need for new drugs. Many STHs use an unusual form of anaerobic metabolism to survive the hypoxic conditions of the host gut. This requires rhodoquinone (RQ), a quinone electron carrier. RQ is not made or used by vertebrate hosts making it an excellent therapeutic target. Here we screen 480 structural families of natural products to find compounds that kill Caenorhabditis elegans specifically when they require RQ-dependent metabolism. We identify several classes of compounds including a family of species-selective inhibitors of mitochondrial respiratory complex I. These identified complex I inhibitors have a benzimidazole core and we determine key structural requirements for activity by screening 1,280 related compounds. Finally, we show several of these compounds kill adult STHs. We suggest these species-selective complex I inhibitors are potential anthelmintics.

Predicting genetic modifier loci using functional gene networks.

Most phenotypes are genetically complex, with contributions from mutations in many different genes. Mutations in more than one gene can combine synergistically to cause phenotypic change, and systematic studies in model organisms show that these genetic interactions are pervasive. However, in human association studies such nonadditive genetic interactions are very difficult to identify because of a lack of statistical power--simply put, the number of potential interactions is too vast. One approach to resolve this is to predict candidate modifier interactions between loci, and then to specifically test these for associations with the phenotype. Here, we describe a general method for predicting genetic interactions based on the use of integrated functional gene networks. We show that in both Saccharomyces cerevisiae and Caenorhabditis elegans a single high-coverage, high-quality functional network can successfully predict genetic modifiers for the majority of genes. For C. elegans we also describe the construction of a new, improved, and expanded functional network, WormNet 2. Using this network we demonstrate how it is possible to rapidly expand the number of modifier loci known for a gene, predicting and validating new genetic interactions for each of three signal transduction genes. We propose that this approach, termed network-guided modifier screening, provides a general strategy for predicting genetic interactions. This work thus suggests that a high-quality integrated human gene network will provide a powerful resource for modifier locus discovery in many different diseases.

A probabilistic view of gene function.

Cells are controlled by the complex and dynamic actions of thousands of genes. With the sequencing of many genomes, the key problem has shifted from identifying genes to knowing what the genes do; we need a framework for expressing that knowledge. Even the most rigorous attempts to construct ontological frameworks describing gene function (e.g., the Gene Ontology project) ultimately rely on manual curation and are thus labor-intensive and subjective. But an alternative exists: the field of functional genomics is piecing together networks of gene interactions, and although these data are currently incomplete and error-prone, they provide a glimpse of a new, probabilistic view of gene function. We outline such a framework, which revolves around a statistical description of gene interactions derived from large, systematically compiled data sets. In this probabilistic view, pleiotropy is implicit, all data have errors and the definition of gene function is an iterative process that ultimately converges on the correct functions. The relationships between the genes are defined by the data, not by hand. Even this comprehensive view fails to capture key aspects of gene function, not least their dynamics in time and space, showing that there are limitations to the model that must ultimately be addressed.

A systematic RNAi screen identifies a critical role for mitochondria in C. elegans longevity.

We report a systematic RNA interference (RNAi) screen of 5,690 Caenorhabditis elegans genes for gene inactivations that increase lifespan. We found that genes important for mitochondrial function stand out as a principal group of genes affecting C. elegans lifespan. A classical genetic screen identified a mutation in the mitochondrial leucyl-tRNA synthetase gene (lrs-2) that impaired mitochondrial function and was associated with longer-lifespan. The long-lived worms with impaired mitochondria had lower ATP content and oxygen consumption, but differential responses to free-radical and other stresses. These data suggest that the longer lifespan of C. elegans with compromised mitochrondria cannot simply be assigned to lower free radical production and suggest a more complex coupling of metabolism and longevity.

Encapsulation of Caenorhabditis elegans in Water-in-Water Microdroplets to Study the Worm Viability: Alternative Avenue to Manipulate Microdroplet Environment.

Due to the great biocompatibility of the aqueous two phase system (ATPS), biological cells have been widely encapsulated in ATPS microdroplets (diameter < 50 µm). However, the immobilization of relatively large multicellular organisms such as Caenorhabditis elegans in ATPS droplets remains challenging as the spontaneous generation of droplets greater than 200 µm is difficult without external perturbations. In this study, we utilize a microneedle-assisted coflow microfludic channel to passively form ATPS microdroplets larger than 200 µm and successfully entrap C. elegans in the microdroplets. We monitor the worm viability and its temporal stroke frequency up to 6 h. We study the effects of dextran (DEX)-to-polyethylene glycol (PEG) flow ratios and worm concentration on the droplet diameter, worm encapsulation efficiency, and the number of droplets containing individual worms. Larger ATPS microdroplets (>200 µm) form in the ranges of capillary number (Ca) between 0.020 to 0.20 and Weber number (We) between 10-5 and 10-3. An ATPS with the encapsulation ability and biocompatibility can offer an alternative immobilization tool for multicellular organisms to existing platforms such as water/oil droplets.

An intestinally secreted host factor promotes microsporidia invasion of C. elegans.

Microsporidia are ubiquitous obligate intracellular pathogens of animals. These parasites often infect hosts through an oral route, but little is known about the function of host intestinal proteins that facilitate microsporidia invasion. To identify such factors necessary for infection by Nematocida parisii, a natural microsporidian pathogen of Caenorhabditis elegans, we performed a forward genetic screen to identify mutant animals that have a Fitness Advantage with Nematocida (Fawn). We isolated four fawn mutants that are resistant to Nematocida infection and contain mutations in T14E8.4, which we renamed aaim-1 (Antibacterial and Aids invasion by Microsporidia). Expression of AAIM-1 in the intestine of aaim-1 animals restores N. parisii infectivity and this rescue of infectivity is dependent upon AAIM-1 secretion. N. parisii spores in aaim-1 animals are improperly oriented in the intestinal lumen, leading to reduced levels of parasite invasion. Conversely, aaim-1 mutants display both increased colonization and susceptibility to the bacterial pathogen Pseudomonas aeruginosa and overexpression ofaaim-1 reduces P. aeruginosa colonization. Competitive fitness assays show that aaim-1 mutants are favored in the presence of N. parisii but disadvantaged on P. aeruginosa compared to wild-type animals. Together, this work demonstrates how microsporidia exploits a secreted protein to promote host invasion. Our results also suggest evolutionary trade-offs may exist to optimizing host defense against multiple classes of pathogens.

Control of feeding behavior in C. elegans by human G protein-coupled receptors permits screening for agonist-expressing bacteria.

G protein-coupled receptors (GPCRs) have a key role in many biological processes and are important drug targets for many human diseases. Therefore, understanding the molecular interactions between GPCRs and their ligands would improve drug design. Here, we describe an approach that allows the rapid identification of functional agonists expressed in bacteria. Transgenic Caenorhabditis elegans expressing the human chemokine receptor 5 (CCR5) in nociceptive neurons show avoidance behavior on encounter with the ligand MIP-1alpha and avoid feeding on Escherichia coli expressing MIP-1alpha compared with control bacteria. This system allows a simple activity screen, based on the distribution of transgenic worms in a binary food-choice assay, without a requirement for protein purification or tagging. By using this approach, a library of 68 MIP-1alpha variants was screened, and 13 critical agonist residues involved in CCR5 activation were identified, four of which (T8, A9, N22, and A25) have not been described previously, to our knowledge. Identified residues were subsequently validated in receptor binding assays and by calcium flux assays in mammalian cells. This approach serves not only for structure/function studies as demonstrated, but may be used to facilitate the discovery of agonists within bacterial libraries.