Reconstructing a metazoan genetic pathway with transcriptome-wide epistasis measurements

DAF-18/PTEN protects LIN-35/Rb from CLP-1/CAPN-mediated cleavage to promote starvation resistance

Starvation resistance is a fundamental trait with profound influence on fitness and disease risk. DAF-18, the Caenorhabditis elegans ortholog of the tumor suppressor PTEN, promotes starvation resistance. PTEN is a dual phosphatase, and DAF-18 promotes starvation resistance as a lipid phosphatase by antagonizing insulin/IGF and PI3K signaling, activating the tumor suppressor DAF-16/FoxO. However, if or how DAF-18/PTEN protein-phosphatase activity promotes starvation resistance is unknown. Using genetic, genomic, bioinformatic, and biochemical approaches, we identified the C. elegans retinoblastoma/RB protein homolog, LIN-35/Rb, as a critical mediator of the effect of DAF-18/PTEN on starvation resistance. We show that DAF-18/PTEN protects LIN-35/Rb from cleavage by the μ-Calpain homolog CLP-1/CAPN, and that LIN-35/Rb together with the repressive DREAM complex promotes starvation resistance. We conclude that the tumor suppressors DAF-18/PTEN and LIN-35/Rb function in a linear pathway, with LIN-35/Rb and the rest of the DREAM complex functioning as a transcriptional effector of DAF-18/PTEN protein-phosphatase activity resulting in repression of germline gene expression. This work is significant for revealing a network of tumor suppressors that promote survival during cellular and developmental quiescence.

DAF-18/PTEN protects LIN-35/Rb from CLP-1/CAPN-mediated cleavage to promote starvation resistance

Starvation resistance is a fundamental trait with profound influence on fitness and disease risk. DAF-18, the C. elegans ortholog of the tumor suppressor PTEN, promotes starvation resistance. PTEN is a dual phosphatase, and DAF-18 promotes starvation resistance as a lipid phosphatase by antagonizing insulin/IGF and PI3K signaling, activating the tumor suppressor DAF-16/FoxO. However, if or how DAF-18/PTEN protein-phosphatase activity promotes starvation resistance is unknown. Using genetic, genomic, bioinformatic, and biochemical approaches, we identified the C. elegans retinoblastoma/RB protein homolog, LIN-35/Rb, as a critical mediator of the effect of DAF-18/PTEN on starvation resistance. We show that DAF-18/PTEN protects LIN-35/Rb from cleavage by the μ-Calpain homolog CLP-1/CAPN, and that LIN-35/Rb together with the repressive DREAM complex promote starvation resistance. We conclude that the tumor suppressors DAF-18/PTEN and LIN-35/Rb function in a linear pathway, with LIN-35/Rb and the rest of the DREAM complex functioning as a transcriptional effector of DAF-18/PTEN protein-phosphatase activity resulting in repression of germline gene expression. This work is significant for revealing a network of tumor suppressors that promote survival during cellular and developmental quiescence.

Diverse epistatic effects in barley-powdery mildew interactions localize to host chromosome hotspots

Barley Mildew locus a (Mla) encodes a multi-allelic series of nucleotide-binding leucine-rich repeat (NLR) receptors that specify recognition to diverse cereal diseases. We exploited time-course transcriptome dynamics of barley and derived immune mutants infected with the powdery mildew fungus, Blumeria hordei (Bh), to infer gene effects governed by Mla6 and two other loci significant to disease development, Blufensin1 (Bln1), and Required for Mla6 resistance3 (rar3 = Sgt1 ΔKL308-309 ). Interactions of Mla6 and Bln1 resulted in diverse epistatic effects on the Bh-induced barley transcriptome, differential immunity to Pseudomonas syringae expressing the effector protease AvrPphB, and reaction to Bh. From a total of 468 barley NLRs, 115 were grouped under different gene effect models; genes classified under these models localized to host chromosome hotspots. The corresponding Bh infection transcriptome was classified into nine co-expressed modules, linking differential expression with pathogen structures, signifying that disease is regulated by an inter-organismal network that diversifies the response.

Patterns of Fitness and Gene Expression Epistasis Generated by Beneficial Mutations in the rho and rpoB Genes of Escherichia coli during High-Temperature Adaptation

Epistasis is caused by genetic interactions among mutations that affect fitness. To characterize properties and potential mechanisms of epistasis, we engineered eight double mutants that combined mutations from the rho and rpoB genes of Escherichia coli. The two genes encode essential functions for transcription, and the mutations in each gene were chosen because they were beneficial for adaptation to thermal stress (42.2 °C). The double mutants exhibited patterns of fitness epistasis that included diminishing returns epistasis at 42.2 °C, stronger diminishing returns between mutations with larger beneficial effects and both negative and positive (sign) epistasis across environments (20.0 °C and 37.0 °C). By assessing gene expression between single and double mutants, we detected hundreds of genes with gene expression epistasis. Previous work postulated that highly connected hub genes in coexpression networks have low epistasis, but we found the opposite: hub genes had high epistasis values in both coexpression and protein-protein interaction networks. We hypothesized that elevated epistasis in hub genes reflected that they were enriched for targets of Rho termination but that was not the case. Altogether, gene expression and coexpression analyses revealed that thermal adaptation occurred in modules, through modulation of ribonucleotide biosynthetic processes and ribosome assembly, the attenuation of expression in genes related to heat shock and stress responses, and with an overall trend toward restoring gene expression toward the unstressed state.

Systematic mapping of organism-scale gene-regulatory networks in aging using population asynchrony

In aging, physiologic networks decline in function at rates that differ between individuals, producing a wide distribution of lifespan. Though 70% of human lifespan variance remains unexplained by heritable factors, little is known about the intrinsic sources of physiologic heterogeneity in aging. To understand how complex physiologic networks generate lifespan variation, new methods are needed. Here, we present Asynch-seq, an approach that uses gene-expression heterogeneity within isogenic populations to study the processes generating lifespan variation. By collecting thousands of single-individual transcriptomes, we capture the Caenorhabditis elegans "pan-transcriptome"-a highly resolved atlas of non-genetic variation. We use our atlas to guide a large-scale perturbation screen that identifies the decoupling of total mRNA content between germline and soma as the largest source of physiologic heterogeneity in aging, driven by pleiotropic genes whose knockdown dramatically reduces lifespan variance. Our work demonstrates how systematic mapping of physiologic heterogeneity can be applied to reduce inter-individual disparities in aging.

Combinatorial transcriptomic and genetic dissection of insulin/IGF-1 signaling-regulated longevity in Caenorhabditis elegans

Classical genetic analysis is invaluable for understanding the genetic interactions underlying specific phenotypes, but requires laborious and subjective experiments to characterize polygenic and quantitative traits. Contrarily, transcriptomic analysis enables the simultaneous and objective identification of multiple genes whose expression changes are associated with specific phenotypes. Here, we conducted transcriptomic analysis of genes crucial for longevity using datasets with daf-2/insulin/IGF-1 receptor mutant Caenorhabditis elegans. Our analysis unraveled multiple epistatic relationships at the transcriptomic level, in addition to verifying genetically established interactions. Our combinatorial analysis also revealed transcriptomic changes associated with longevity conferred by daf-2 mutations. In particular, we demonstrated that the extent of lifespan changes caused by various mutant alleles of the longevity transcription factor daf-16/FOXO matched their effects on transcriptomic changes in daf-2 mutants. We identified specific aging-regulating signaling pathways and subsets of structural and functional RNA elements altered by different genes in daf-2 mutants. Lastly, we elucidated the functional cooperation between several longevity regulators, based on the combination of transcriptomic and molecular genetic analysis. These data suggest that different biological processes coordinately exert their effects on longevity in biological networks. Together our work demonstrates the utility of transcriptomic dissection analysis for identifying important genetic interactions for physiological processes, including aging and longevity.

Transcriptome analyses describe the consequences of persistent HIF-1 over-activation in Caenorhabditis elegans

Metazoan animals rely on oxygen for survival, but during normal development and homeostasis, animals are often challenged by hypoxia (low oxygen). In metazoans, many of the critical hypoxia responses are mediated by the evolutionarily conserved hypoxia-inducible transcription factors (HIFs). The stability and activity of HIF complexes are strictly regulated. In the model organism C. elegans, HIF-1 stability and activity are negatively regulated by VHL-1, EGL-9, RHY-1 and SWAN-1. Importantly, C. elegans mutants carrying strong loss-of-function mutations in these genes are viable, and this provides opportunities to interrogate the molecular consequences of persistent HIF-1 over-activation. We find that the genome-wide gene expression patterns are compellingly similar in these mutants, supporting models in which RHY-1, VHL-1 and EGL-9 function in common pathway(s) to regulate HIF-1 activity. These studies illuminate the diversified biological roles played by HIF-1, including metabolism and stress response. Genes regulated by persistent HIF-1 over-activation overlap with genes responsive to pathogens, and they overlap with genes regulated by DAF-16. As crucial stress regulators, HIF-1 and DAF-16 converge on key stress-responsive genes and function synergistically to enable hypoxia survival.

Transcriptome analyses describe the consequences of persistent HIF-1 over-activation in Caenorhabditis elegans

Metazoan animals rely on oxygen for survival, but during normal development and homeostasis, animals are often challenged by hypoxia (low oxygen). In metazoans, many of the critical hypoxia responses are mediated by the evolutionarily conserved hypoxia-inducible transcription factors (HIFs). The stability and activity of HIF complexes are strictly regulated. In the model organism C. elegans, HIF-1 stability and activity are negatively regulated by VHL-1, EGL-9, RHY-1 and SWAN-1. Importantly, C. elegans mutants carrying strong loss-of-function mutations in these genes are viable, and this provides opportunities to interrogate the molecular consequences of persistent HIF-1 over-activation. We find that the genome-wide gene expression patterns are compellingly similar in these mutants, supporting models in which RHY-1, SWAN-1 and EGL-9 function in common pathway(s) to regulate HIF-1 activity. These studies illuminate the diversified biological roles played by HIF-1, including metabolism, hypoxia and other stress responses, reproduction and development. Genes regulated by persistent HIF-1 over-activation overlap with genes responsive to pathogens, and they overlap with genes regulated by DAF-16. As crucial stress regulators, HIF-1 and DAF-16 converge on key stress-responsive genes and function synergistically to enable hypoxia survival.

A Caenorhabditis elegans Male Pheromone Feminizes Germline Gene Expression in Hermaphrodites and Imposes Life-History Costs

Sex pheromones not only improve the reproductive success of the recipients, but also impose costs, such as a reduced life span. The underlying mechanisms largely remain to be elucidated. Here, we show that even a brief exposure to physiological amounts of the dominant Caenorhabditis elegans male pheromone, ascr#10, alters the expression of thousands of genes in hermaphrodites. The most dramatic effect on the transcriptome is the upregulation of genes expressed during oogenesis and the downregulation of genes associated with male gametogenesis. This result reveals a way in which social signals help to resolve the inherent conflict between spermatogenesis and oogenesis in a simultaneous hermaphrodite, presumably to optimally align reproductive function with the presence of potential mating partners. We also found that exposure to ascr#10 increased the risk of persistent intestinal infections in hermaphrodites due to pathological pharyngeal hypertrophy. Thus, our study reveals ways in which the male pheromone can not only have beneficial effects on the recipients' reproduction, but also cause harmful consequences that reduce life span.

A C. elegans male pheromone feminizes germline gene expression in hermaphrodites and imposes life-history costs

Sex pheromones improve reproductive success, but also impose costs. Here we show that even brief exposure to physiological amounts of the dominant C. elegans male pheromone, ascr#10, alters the expression of thousands of genes in hermaphrodites. The most dramatic effect on the transcriptome was the upregulation of genes expressed during oogenesis and downregulation of genes associated with male gametogenesis. Among the detrimental effects of ascr#10 on hermaphrodites is the increased risk of persistent infections caused by pathological pharyngeal hypertrophy. Our results reveal a way in which social signals help to resolve the inherent conflict between spermatogenesis and oogenesis in a simultaneous hermaphrodite, presumably to optimally align reproductive function to the presence of potential mating partners. They also show that the beneficial effects of the pheromone are accompanied by harmful consequences that reduce lifespan.

Neuronal temperature perception induces specific defenses that enable C. elegans to cope with the enhanced reactivity of hydrogen peroxide at high temperature

Hydrogen peroxide is the most common reactive chemical that organisms face on the microbial battlefield. The rate with which hydrogen peroxide damages biomolecules required for life increases with temperature, yet little is known about how organisms cope with this temperature-dependent threat. Here, we show that Caenorhabditis elegans nematodes use temperature information perceived by sensory neurons to cope with the temperature-dependent threat of hydrogen peroxide produced by the pathogenic bacterium Enterococcus faecium. These nematodes preemptively induce the expression of specific hydrogen peroxide defenses in response to perception of high temperature by a pair of sensory neurons. These neurons communicate temperature information to target tissues expressing those defenses via an insulin/IGF1 hormone. This is the first example of a multicellular organism inducing their defenses to a chemical when they sense an inherent enhancer of the reactivity of that chemical.

Dietary-derived vitamin B12 protects Caenorhabditis elegans from thiol-reducing agents

BACKGROUND

One-carbon metabolism, which includes the folate and methionine cycles, involves the transfer of methyl groups which are then utilised as a part of multiple physiological processes including redox defence. During the methionine cycle, the vitamin B12-dependent enzyme methionine synthetase converts homocysteine to methionine. The enzyme S-adenosylmethionine (SAM) synthetase then uses methionine in the production of the reactive methyl carrier SAM. SAM-binding methyltransferases then utilise SAM as a cofactor to methylate proteins, small molecules, lipids, and nucleic acids.

RESULTS

We describe a novel SAM methyltransferase, RIPS-1, which was the single gene identified from forward genetic screens in Caenorhabditis elegans looking for resistance to lethal concentrations of the thiol-reducing agent dithiothreitol (DTT). As well as RIPS-1 mutation, we show that in wild-type worms, DTT toxicity can be overcome by modulating vitamin B12 levels, either by using growth media and/or bacterial food that provide higher levels of vitamin B12 or by vitamin B12 supplementation. We show that active methionine synthetase is required for vitamin B12-mediated DTT resistance in wild types but is not required for resistance resulting from RIPS-1 mutation and that susceptibility to DTT is partially suppressed by methionine supplementation. A targeted RNAi modifier screen identified the mitochondrial enzyme methylmalonyl-CoA epimerase as a strong genetic enhancer of DTT resistance in a RIPS-1 mutant. We show that RIPS-1 is expressed in the intestinal and hypodermal tissues of the nematode and that treating with DTT, β-mercaptoethanol, or hydrogen sulfide induces RIPS-1 expression. We demonstrate that RIPS-1 expression is controlled by the hypoxia-inducible factor pathway and that homologues of RIPS-1 are found in a small subset of eukaryotes and bacteria, many of which can adapt to fluctuations in environmental oxygen levels.

CONCLUSIONS

This work highlights the central importance of dietary vitamin B12 in normal metabolic processes in C. elegans, defines a new role for this vitamin in countering reductive stress, and identifies RIPS-1 as a novel methyltransferase in the methionine cycle.

Defining characteristics and conservation of poorly annotated genes in Caenorhabditis elegans using WormCat 2.0

Omics tools provide broad datasets for biological discovery. However, the computational tools for identifying important genes or pathways in RNA-seq, proteomics, or GWAS (Genome-Wide Association Study) data depend on Gene Ontogeny annotations and are biased toward well-described pathways. This limits their utility as poorly annotated genes, which could have novel functions, are often passed over. Recently, we developed an annotation and category enrichment tool for Caenorhabditis elegans genomic data, WormCat, which provides an intuitive visualization output. Unlike Gene Ontogeny-based enrichment tools, which exclude genes with no annotation information, WormCat 2.0 retains these genes as a special UNASSIGNED category. Here, we show that the UNASSIGNED gene category enrichment exhibits tissue-specific expression patterns and can include genes with biological functions identified in published datasets. Poorly annotated genes are often considered to be potentially species-specific and thus, of reduced interest to the biomedical community. Instead, we find that around 3% of the UNASSIGNED genes have human orthologs, including some linked to human diseases. These human orthologs themselves have little annotation information. A recently developed method that incorporates lineage relationships (abSENSE) indicates that the failure of BLAST to detect homology explains the apparent lineage specificity for many UNASSIGNED genes. This suggests that a larger subset could be related to human genes. WormCat provides an annotation strategy that allows the association of UNASSIGNED genes with specific phenotypes and known pathways. Building these associations in C. elegans, with its robust genetic tools, provides a path to further functional study and insight into these understudied genes.

Structural and functional relationship of mammalian and nematode ferritins

Ferritin is a unique buffering protein in iron metabolism. By storing or releasing iron in a tightly controlled manner, it prevents the negative effects of free ferrous ions on biomolecules in all domains of life - from bacteria to mammals. This review focuses on the structural features and activity of the ferritin protein family with an emphasis on nematode ferritins and the similarities in their biological roles with mammalian ferritins. The conservative characteristic of the ferritin family across the species originates from the ferroxidase activity against redox-active iron. The antioxidative function of these proteins translates into their involvement in a wide range of important biological processes, e.g., aging, fat metabolism, immunity, anticancer activity, and antipathogenic activity. Moreover, disturbances in ferritin expression lead to severe iron-associated diseases. Research on the Caenorhabditis elegans model organism may allow us to better understand the wide spectrum of mechanisms involving ferritin activity.

Crosstalk in oxygen homeostasis networks: SKN-1/NRF inhibits the HIF-1 hypoxia-inducible factor in Caenorhabditis elegans

During development, homeostasis, and disease, organisms must balance responses that allow adaptation to low oxygen (hypoxia) with those that protect cells from oxidative stress. The evolutionarily conserved hypoxia-inducible factors are central to these processes, as they orchestrate transcriptional responses to oxygen deprivation. Here, we employ genetic strategies in C. elegans to identify stress-responsive genes and pathways that modulate the HIF-1 hypoxia-inducible factor and facilitate oxygen homeostasis. Through a genome-wide RNAi screen, we show that RNAi-mediated mitochondrial or proteasomal dysfunction increases the expression of hypoxia-responsive reporter Pnhr-57::GFP in C. elegans. Interestingly, only a subset of these effects requires hif-1. Of particular importance, we found that skn-1 RNAi increases the expression of hypoxia-responsive reporter Pnhr-57::GFP and elevates HIF-1 protein levels. The SKN-1/NRF transcription factor has been shown to promote oxidative stress resistance. We present evidence that the crosstalk between HIF-1 and SKN-1 is mediated by EGL-9, the prolyl hydroxylase that targets HIF-1 for oxygen-dependent degradation. Treatment that induces SKN-1, such as heat or gsk-3 RNAi, increases expression of a Pegl-9::GFP reporter, and this effect requires skn-1 function and a putative SKN-1 binding site in egl-9 regulatory sequences. Collectively, these data support a model in which SKN-1 promotes egl-9 transcription, thereby inhibiting HIF-1. We propose that this interaction enables animals to adapt quickly to changes in cellular oxygenation and to better survive accompanying oxidative stress.

Protein Complexes Form a Basis for Complex Hybrid Incompatibility

Proteins are the workhorses of the cell and execute many of their functions by interacting with other proteins forming protein complexes. Multi-protein complexes are an admixture of subunits, change their interaction partners, and modulate their functions and cellular physiology in response to environmental changes. When two species mate, the hybrid offspring are usually inviable or sterile because of large-scale differences in the genetic makeup between the two parents causing incompatible genetic interactions. Such reciprocal-sign epistasis between inter-specific alleles is not limited to incompatible interactions between just one gene pair; and, usually involves multiple genes. Many of these multi-locus incompatibilities show visible defects, only in the presence of all the interactions, making it hard to characterize. Understanding the dynamics of protein-protein interactions (PPIs) leading to multi-protein complexes is better suited to characterize multi-locus incompatibilities, compared to studying them with traditional approaches of genetics and molecular biology. The advances in omics technologies, which includes genomics, transcriptomics, and proteomics can help achieve this end. This is especially relevant when studying non-model organisms. Here, we discuss the recent progress in the understanding of hybrid genetic incompatibility; omics technologies, and how together they have helped in characterizing protein complexes and in turn multi-locus incompatibilities. We also review advances in bioinformatic techniques suitable for this purpose and propose directions for leveraging the knowledge gained from model-organisms to identify genetic incompatibilities in non-model organisms.

Hypoxic response regulators RHY-1 and EGL-9/PHD promote longevity through a VHL-1-independent transcriptional response

HIF-1-mediated adaptation to changes in oxygen availability is a critical aspect of healthy physiology. HIF is regulated by a conserved mechanism whereby EGLN/PHD family members hydroxylate HIF in an oxygen-dependent manner, targeting it for ubiquitination by Von-Hippel-Lindau (VHL) family members, leading to its proteasomal degradation. The activity of the only C. elegans PHD family member, EGL-9, is also regulated by a hydrogen sulfide sensing cysteine-synthetase-like protein, CYSL-1, which is, in turn, regulated by RHY-1/acyltransferase. Over the last decade, multiple seminal studies have established a role for the hypoxic response in regulating longevity, with mutations in vhl-1 substantially extending C. elegans lifespan through a HIF-1-dependent mechanism. However, studies on other components of the hypoxic signaling pathway that similarly stabilize HIF-1 have shown more mixed results, suggesting that mutations in egl-9 and rhy-1 frequently fail to extend lifespan. Here, we show that egl-9 and rhy-1 mutants suppress the long-lived phenotype of vhl-1 mutants. We also show that RNAi of rhy-1 extends lifespan of wild-type worms while decreasing lifespan of vhl-1 mutant worms. We further identify VHL-1-independent gene expression changes mediated by EGL-9 and RHY-1 and find that a subset of these genes contributes to longevity regulation. The resulting data suggest that changes in HIF-1 activity derived by interactions with EGL-9 likely contribute greatly to its role in regulation of longevity.

Ferritin is regulated by a neuro-intestinal axis in the nematode Caenorhabditis elegans

Iron is vital for the life of most organisms. However, when dysregulated, iron can catalyze the formation of oxygen (O₂) radicals that can destroy any biological molecule and thus lead to oxidative injury and death. Therefore, iron metabolism must be tightly regulated at all times, as well as coordinated with the metabolism of O₂. However, how is this achieved at the whole animal level is not well understood. Here, we explore this question using the nematode Caenorhabditis elegans. Exposure of worms to O₂ starvation conditions (i.e. hypoxia) induces a major upregulation in levels of the conserved iron-cage protein ferritin 1 (ftn-1) in the intestine, while exposure to 21% O₂ decreases ftn-1 level. This O₂-dependent inhibition is mediated by O₂-sensing neurons that communicate with the intestine through neurotransmitter and neuropeptide signalling, and requires the activity of hydroxylated HIF-1. By contrast, the induction of ftn-1 in hypoxia appears to be HIF-1-independent. This upregulation provides protection against Pseudomonas aeruginosa bacteria and oxidative injury. Taken together, our studies uncover a neuro-intestine axis that coordinates O₂ and iron responses at the whole animal level.

•

The expression of ferritin 1 (ftn-1) is tightly regulated by O₂ tension.

•

O₂-sensing neurons inhibit the expression of ftn-1 in the intestine at 21% O₂.

•

Hydroxylated–HIF–1 inhibits the expression of ftn-1 at 21% O₂.

•

ftn-1 is important for protecting against Pseudomonas aeruginosa bacteria.

Highly Combinatorial Genetic Interaction Analysis Reveals a Multi-Drug Transporter Influence Network

Many traits are complex, depending non-additively on variant combinations. Even in model systems, such as the yeast S. cerevisiae, carrying out the high-order variant-combination testing needed to dissect complex traits remains a daunting challenge. Here, we describe “X-gene” genetic analysis (XGA), a strategy for engineering and profiling highly combinatorial gene perturbations. We demonstrate XGA on yeast ABC transporters by engineering 5,353 strains, each deleted for a random subset of 16 transporters, and profiling each strain’s resistance to 16 compounds. XGA yielded 85,648 genotype-to-resistance observations, revealing high-order genetic interactions for 13 of the 16 transporters studied. Neural networks yielded intuitive functional models and guided exploration of fluconazole resistance, which was influenced non-additively by five genes. Together, our results showed that highly combinatorial genetic perturbation can functionally dissect complex traits, supporting pursuit of analogous strategies in human cells and other model systems.

Celaj et al. introduce “X-gene” genetic analysis (XGA), a strategy for modeling complex systems by engineering and profiling highly combinatorial genetic perturbations. They apply XGA to 16 yeast ABC transporters, revealing many high-order genetic interactions. Neural network models yielded intuitive functional models and illuminated an ABC transporter influence network, supporting application of XGA to other organisms and processes.

Common and distinct transcriptional signatures of mammalian embryonic lethality

The Deciphering the Mechanisms of Developmental Disorders programme has analysed the morphological and molecular phenotypes of embryonic and perinatal lethal mouse mutant lines in order to investigate the causes of embryonic lethality. Here we show that individual whole-embryo RNA-seq of 73 mouse mutant lines (>1000 transcriptomes) identifies transcriptional events underlying embryonic lethality and associates previously uncharacterised genes with specific pathways and tissues. For example, our data suggest that Hmgxb3 is involved in DNA-damage repair and cell-cycle regulation. Further, we separate embryonic delay signatures from mutant line-specific transcriptional changes by developing a baseline mRNA expression catalogue of wild-type mice during early embryogenesis (4-36 somites). Analysis of transcription outside coding sequence identifies deregulation of repetitive elements in Morc2a mutants and a gene involved in gene-specific splicing. Collectively, this work provides a large scale resource to further our understanding of early embryonic developmental disorders.

Using Transcriptomes as Mutant Phenotypes Reveals Functional Regions of a Mediator Subunit in Caenorhabditis elegans

Although transcriptomes have recently been used as phenotypes with which to perform epistasis analyses, they are not yet used to study intragenic function/structure relationships. We developed a theoretical framework to study allelic series using transcriptomic phenotypes. As a proof-of-concept, we apply our methods to an allelic series of dpy-22, a highly pleiotropic Caenorhabditis elegans gene orthologous to the human gene MED12, which encodes a subunit of the Mediator complex. Our methods identify functional units within dpy-22 that modulate Mediator activity upon various genetic programs, including the Wnt and Ras modules.

Hypoxia-inducible factor cell non-autonomously regulates C. elegans stress responses and behavior via a nuclear receptor

The HIF (hypoxia-inducible factor) transcription factor is the master regulator of the metazoan response to chronic hypoxia. In addition to promoting adaptations to low oxygen, HIF drives cytoprotective mechanisms in response to stresses and modulates neural circuit function. How most HIF targets act in the control of the diverse aspects of HIF-regulated biology remains unknown. We discovered that a HIF target, the C. elegans gene cyp-36A1, is required for numerous HIF-dependent processes, including modulation of gene expression, stress resistance, and behavior. cyp-36A1 encodes a cytochrome P450 enzyme that we show controls expression of more than a third of HIF-induced genes. CYP-36A1 acts cell non-autonomously by regulating the activity of the nuclear hormone receptor NHR-46, suggesting that CYP-36A1 functions as a biosynthetic enzyme for a hormone ligand of this receptor. We propose that regulation of HIF effectors through activation of cytochrome P450 enzyme/nuclear receptor signaling pathways could similarly occur in humans.

The Caenorhabditis elegans Female-Like State: Decoupling the Transcriptomic Effects of Aging and Sperm Status

Understanding genome and gene function in a whole organism requires us to fully comprehend the life cycle and the physiology of the organism in question. Caenorhabditis elegans XX animals are hermaphrodites that exhaust their sperm after 3 d of egg-laying. Even though C. elegans can live for many days after cessation of egg-laying, the molecular physiology of this state has not been as intensely studied as other parts of the life cycle, despite documented changes in behavior and metabolism. To study the effects of sperm depletion and aging of C. elegans during the first 6 d of adulthood, we measured the transcriptomes of first-day adult hermaphrodites and sixth-day sperm-depleted adults, and, at the same time points, mutant fog-2(lf) worms that have a feminized germline phenotype. We found that we could separate the effects of biological aging from sperm depletion. For a large subset of genes, young adult fog-2(lf) animals had the same gene expression changes as sperm-depleted sixth-day wild-type hermaphrodites, and these genes did not change expression when fog-2(lf) females reached the sixth day of adulthood. Taken together, this indicates that changing sperm status causes a change in the internal state of the worm, which we call the female-like state. Our data provide a high-quality picture of the changes that happen in global gene expression throughout the period of early aging in the worm.