Method | Open | Published:
Synthetic lethal analysis of Caenorhabditis elegans posterior embryonic patterning genes identifies conserved genetic interactions
Genome Biologyvolume 6, Article number: R45 (2005)
Phenotypic robustness is evidenced when single-gene mutations do not result in an obvious phenotype. It has been suggested that such phenotypic stability results from 'buffering' activities of homologous genes as well as non-homologous genes acting in parallel pathways. One approach to characterizing mechanisms of phenotypic robustness is to identify genetic interactions, specifically, double mutants where buffering is compromised. To identify interactions among genes implicated in posterior patterning of the Caenorhabditis elegans embryo, we measured synthetic lethality following RNA interference of 22 genes in 15 mutant strains. A pair of homologous T-box transcription factors (tbx-8 and tbx-9) is found to interact in both C. elegans and C. briggsae, indicating that their compensatory function is conserved. Furthermore, a muscle module is defined by transitive interactions between the MyoD homolog hlh-1, another basic helix-loop-helix transcription factor, hnd-1, and the MADS-box transcription factor unc-120. Genetic interactions within a homologous set of genes involved in vertebrate myogenesis indicate broad conservation of the muscle module and suggest that other genetic modules identified in C. elegans will be conserved.
Forward and reverse genetic screens in flies and worms indicate that most genes are not essential to the development or viability of the organism [1–4]. Two primary explanations for such phenotypic robustness to mutation have been offered: homologous gene products may directly compensate for one another's function, or indirect compensation of function may occur through regulatory networks via non-homologous genes acting in alternative pathways or feedback mechanisms [5, 6]. In both Saccharomyces cerevisiae and C. elegans, genes with at least one homolog are less likely than unique genes to have a loss-of-function phenotype [7, 8]. However, homology accounts for no more than two thirds of the observed phenotypic robustness to mutation in S. cerevisiae and even less in C. elegans, indicating a significant role of the regulatory network in genetic buffering [7, 8]. By identifying gene disruptions that are viable in wild type but lethal in a specific mutant background, synthetic lethal screens can shed light on how regulatory networks buffer gene function . Although genetic interactions are being mapped on a genome-wide scale in yeast , no such efforts have been reported for an animal system. We report here the use of existing mutants and RNA interference (RNAi) to assemble a synthetic lethal matrix in C. elegans. Our aim was to build on prior knowledge by focusing on a characterized set of genes likely to have interactions among them.
The C blastomere is one of five somatic founder blastomeres in the C. elegans embryo. Each founder blastomere produces a characteristic set of cell types via an invariant lineage; the C blastomere predominantly gives rise to muscle and epidermal cells . Maternal and zygotic activities of the homeodomain protein PAL-1 specify the identity and maintain the development of the C-blastomere lineage [11, 12]. We have identified genes whose expression depends either directly or indirectly on PAL-1 function ('PAL-1 targets') and shown that they affect specification, differentiation, and morphogenesis of C-lineage cells. PAL-1 targets include many transcription factors and based on their temporal and spatial expression patterns in wild type, we have proposed a model for their regulatory relationships . However, RNAi of most PAL-1 targets does not result in a detectable phenotype . Because the regulatory network specified by PAL-1 appears to function autonomously while controlling multiple, discrete developmental processes (that is, specification and morphogenesis of muscle and epidermal cells) [11, 13, 14], we reasoned that the network is organized in a modular fashion and that PAL-1 targets may have overlapping or compensatory functions. We set out to measure synthetic lethality for as many pairs of PAL-1 targets as possible, and we included three additional genes encoding transcriptional regulators (lin-26, unc-62, and ceh-40) whose expression patterns and phenotypes suggest that they may have genetic interactions with PAL-1 or its targets [15, 16].
Results and discussion
A synthetic lethal matrix
We began our analysis with 26 genes, 23 PAL-1 targets and three genes likely to have genetic interactions with PAL-1 targets. Mutations in 15 of these genes have been isolated, and appropriate strains were obtained from the Caenorhabditis Genetics Center (Table 1). To screen for synthetic lethal interactions, RNAi of the 13 nonessential genes was performed and embryonic lethality scored in each of the mutant strains as well as wild type to generate a synthetic lethal matrix (Figure 1 and Additional data file 1). However, RNAi of a given gene could not be used if it results in 100% lethality in wild type (pal-1, elt-1, unc-62, and lin-26), and so the matrix is not entirely symmetric. To reduce the penetrance of lethality following RNAi for these four genes we used a lower concentration of double-stranded RNA (dsRNA), but the results were variable (data not shown), and we decided that it was an unreliable approach. We also tried to assemble the matrix using exclusively RNAi, but we found that soaking worms in two different dsRNA molecules targeting two different genes at once to be unacceptably inefficient (Additional data file 2).
Given the effects of mutation alone without RNAi, and the effects of RNAi in wild type, it may not be obvious from the raw data which disruptions result in significantly more lethality when paired than when alone. Measurements of lethality were thus duplicated and a t-test was used to assign statistical significance to each of the interactions tested (Figure 1 and Additional data file 1). Seven genetic interactions were identified (P < 0.001) among six genes out of 195 total interactions tested. In addition to the genes shown in Figure 1, RNAi of the remaining genes for which no mutant was available was performed in all 16 strains. Because synthetic lethality was measured only a single time for each of these possible interactions, the data could not be subjected to statistical analysis. However, RNAi of these additional genes did not appear to elevate lethality in any of the genetic backgrounds tested (data not shown).
Conservation of a genetic interaction
The T-box family of proteins regulates diverse aspects of animal development, including cell fate specification, differentiation, and morphogenesis . The C. elegans genome contains 21 of these animal-specific genes, compared to 32 in C. briggsae and eight in D. melanogaster [18, 19], suggesting duplication or retention within nematodes. tbx-37 and tbx-38 have been shown to redundantly control embryonic patterning of the C. elegans AB lineage , and our systematic approach identified tbx-8 and tbx-9 as a synthetic lethal pair. tbx-8 and tbx-9 are expressed in nearly identical temporal and spatial patterns in the early embryo, with the strongest expression detected in dorsal epidermal cells and posterior expression depending on pal-1 function [13, 21, 22]. tbx-8(ok656), tbx-8(RNAi), and tbx-9(RNAi) each result in a low frequency of embryonic lethality (less than 5%) with 1-10% of hatching larvae displaying posterior morphological defects (Figure 2) [21, 22]. However, 50% of tbx-8(ok656); tbx-9(RNAi) embryos fail to hatch and those that do display severe morphological defects and die as larvae (Figure 2). Phenotypic analysis demonstrates that tbx-8 and tbx-9 are required for intercalation of dorsal epidermal cells thus explaining the morphological defects observed following their disruption .
The overlapping function of tbx-8 and tbx-9 is not the residual result of a recent gene duplication event. The C. briggsae genome contains a pair of syntenic orthologs (by best reciprocal blastp match) , and as in C. elegans, the C. briggsae orthologs have overlapping function during morphogenesis (Table 2), indicating that the functional relationship between this pair of genes has been selectively maintained for over 100 million years. This observation suggests that the overlapping functions of other pairs of T-box genes, which are likely to be common throughout the Metazoa, may also be conserved.
A muscle differentiation module
It has been suggested that modules of genes comprise units of biological function . Because they operate in concerted fashion, genes in a module are expected to have increased likelihood of being co-expressed and sharing physical and genetic interactions relative to a random set of genes. Clustering a large synthetic lethal matrix in S. cerevisiae identified multiple functional modules , and we identify a module in our matrix around the hlh-1 gene (Figure 1). We detect five genetic interactions (P < 0.001) out of six tested among hlh-1, hnd-1, and unc-120. All three genes encode transcription factors expressed exclusively in the muscle progenitor cells of the early embryo [13, 24, 25]. Disruption of function by mutation or RNAi of any one of the three genes results in a low frequency of embryonic lethality with embryos arresting paralyzed at the two-fold stage (Pat) (Figure 3). Embryos require muscle contraction to elongate and pass from the two-fold to three-fold embryonic stage, and the Pat phenotype results from a complete failure of muscle differentiation [26, 27]. A large fraction of animals that do hatch after disruption of hlh-1, unc-120, or hnd-1 show a less severe phenotype (uncoordinated, dumpy larvae) likely reflecting partial muscle function [26, 27] (Figure 3). However, disruption of function of any two of the three genes significantly elevates the frequency of Pat embryos, with disruption of hlh-1 and unc-120 resulting in 100% Pat embryos (Figure 3). The fact that hlh-1, unc-120, and hnd-1 have identical individual and synthetic phenotypes, in addition to their identical early embryonic expression patterns, indicate that they comprise a muscle-differentiation module. Although fewer genes define this module than those identified in larger screens in yeast , the specificity of its defining phenotype indicates that hlh-1, unc-120, and hnd-1 function in concert, and it is likely that future synthetic screens will identify additional members of the module.
If both RNAi and genetic disruptions of gene function are null, the synthetic lethal matrix should be symmetric. However, we can only test the pair of reciprocal interactions when mutations of both genes are available and RNAi of neither gene results in 100% lethality in wild type. Although there is some symmetry in the matrix around hlh-1 (Figure 1), the interactions detected are not all reciprocal (Figure 3). The efficiency of RNAi varies from gene to gene and often fails to result in null phenotypes, suggesting a significant number of false negatives among the interactions tested and explaining the imperfect symmetry among the interactions detected. The pattern of interactions defining the muscle module suggests that hlh-1 is the most essential of the three genes, or that its function is the most potent, and hnd-1 the least (Figure 3), a model that accounts for the reported lack of genetic interaction between heterozygous hlh-1 and homozygous hnd-1 mutations .
The functional importance of the muscle module we identified is underscored by its conservation. Well characterized genetic interactions exist among the vertebrate homologs of hlh-1 - the myogenic basic helix-loop-helix (bHLH) regulatory factors MyoD, myf5, myogenin, and MRF4 . In addition, interactions between vertebrate myogenic bHLH genes and homologs of unc-120 (the MEF2 group of MADS-box regulators) have been described . Interactions among and between vertebrate myogenic bHLH and MADS-box regulators is likely to be the result of cooperative activation of common target genes as well as their ability to activate each other transcriptionally [28–30]. It is possible that hlh-1, hnd-1, and unc-120, like their vertebrate counterparts, display modular genetic interactions because they share common target genes, and they may also activate each other transcriptionally.
We show that the facility of RNAi in C. elegans can be used to screen a large number of genes for synthetic phenotypes. Using an unbiased experimental approach we identified both easy-to-predict genetic interactions between similar co-expressed genes and difficult-to-predict genetic interactions between non-homologous genes. We also used reverse genetics in C. briggsae to show that genetic interaction between a pair of co-expressed homologous transcription factors (tbx-8 and tbx-9) is conserved, suggesting an adaptive role of overlapping function as opposed to the transient result of gene duplication. Furthermore, we identify a genetic module defined by reciprocal interactions among three transcription factors that is required for differentiation of muscle. A similar network among homologous genes is conserved in vertebrates, suggesting that a single master regulator cannot robustly control muscle development. Given the wealth of existing mutants and the ease of RNAi in C. elegans, our approach can be scaled up to identify the core set of metazoan developmental genetic interactions.
Materials and methods
DNA templates for in vitro transcription were prepared by amplification of 100-1000 base pair (bp) sequences from cDNA by PCR. PCR primers were designed to amplify unique sequences in order to ensure specificity of RNAi (Additional data file 3), and a minimal T7 promoter sequence was added to the 5' end of both primers. PCR DNA was purified on a DNA Clean and Concentrator-5 column (Zymo Research), and a high-yield in vitro transcription kit (Promega) was used to produce dsRNA. The dsRNA was reannealed by heating to 90°C and then cooling 6°C/min to 25°C before organic extraction with a 1:1 mixture of phenol and chloroform. Following ethanol precipitation, dsRNA was dissolved in Soaking Buffer  at high concentration, which was then measured by spectrophotometer and diluted to 5 mg/ml.
Eight to twelve larval stage 4 (L4) worms were soaked in 0.75 μl 5 mg/ml dsRNA for 24 h at 20°C and then transferred to a recovery plate with Escherichia coli strain OP50 for 12 h at 25°C and then transferred to a score plate with OP50. After 12 h at 25°C on the score plate, adults were removed, and after another 24 h at 25°C, embryos and larvae were counted under a dissecting microscope (with the exception of unc-120(st364) which is temperature-sensitive and was therefore scored at 15°C). For strains carrying viable alleles, embryonic lethality was scored for the homozygous progeny of homozygous parents. For strains carrying 100% lethal alleles, embryonic lethality was scored for the mixed progeny of either heterozygous or hemizygous parents. For RNAi in C. briggsae, 1 mg/ml dsRNA was injected into the germline of young wild-type adults (strain AF16), and the worms were rescued and scored in the same way as the soaked C. elegans.
For statistical analysis of the data, percent embryonic lethality was converted to survival (1 - % embryonic lethality). All survival values were normalized by the survival of N2 (wild type) in Soaking Buffer (no dsRNA). For each genetic interaction tested, the null hypothesis used is that the survival of the double disruption (mutation and RNAi) equals the product of survivals for each single disruption (mutation alone × RNAi alone). The pairs of duplicate survival measurements for each single disruption are multiplied to generate four values for the expected survival of the double disruption, and a t-test is used to determine if the observed survival (n = 2) is significantly different from expected (n = 4). For example, if two genes function independently, and the survival following disruption of each alone is 90%, then our model expects 81% survival following double disruption since the 90% of individuals surviving one disruption have a 10% chance of dying from the other disruption. If survival is significantly different from the expected 81% then we conclude that the two genes interact.
For microscopy, embryos were mounted on 2% agar pads under coverslips and allowed to develop in a humid box for 12-15 h at 25°C. A Zeiss Axiophot with 100× differential interference contrast optics was used and digital images captured.
Additional data files
The following additional data are available with the online version of this paper. Additional data file 1 contains replicate embryonic lethality measurements and corresponding P-values for synthetic lethality. Additional data file 2 contains data showing that double RNAi by soaking worms in dsRNA for two genes at once is an inefficient means of synthetic genetic analysis. Additional data file 3 lists primer sequences used in this study to amplify coding sequences from cDNA for use as templates for dsRNA.
Nusslein-Volhard C: Of flies and fishes. Science. 1994, 266: 572-574.
Riddle DL: C. elegans II. 1997, Plainview, NY: Cold Spring Harbor Laboratory Press
Boutros M, Kiger AA, Armknecht S, Kerr K, Hild M, Koch B, Haas SA, Consortium HF, Paro R, Perrimon N: Genome-wide RNAi analysis of growth and viability in Drosophila cells. Science. 2004, 303: 832-835. 10.1126/science.1091266.
Kamath RS, Ahringer J: Genome-wide RNAi screening in Caenorhabditis elegans. Methods. 2003, 30: 313-321. 10.1016/S1046-2023(03)00050-1.
Hartman JLt, Garvik B, Hartwell L: Principles for the buffering of genetic variation. Science. 2001, 291: 1001-1004. 10.1126/science.291.5506.1001.
Wagner A: Robustness against mutations in genetic networks of yeast. Nat Genet. 2000, 24: 355-361. 10.1038/74174.
Gu Z, Steinmetz LM, Gu X, Scharfe C, Davis RW, Li WH: Role of duplicate genes in genetic robustness against null mutations. Nature. 2003, 421: 63-66. 10.1038/nature01198.
Conant GC, Wagner A: Duplicate genes and robustness to transient gene knock-downs in Caenorhabditis elegans. Proc R Soc Lond B Biol Sci. 2004, 271: 89-96. 10.1098/rspb.2003.2560.
Tong AH, Lesage G, Bader GD, Ding H, Xu H, Xin X, Young J, Berriz GF, Brost RL, Chang M, et al: Global mapping of the yeast genetic interaction network. Science. 2004, 303: 808-813. 10.1126/science.1091317.
Sulston JE, Schierenberg E, White JG, Thomson JN: The embryonic cell lineage of the nematode Caenorhabditis elegans. Dev Biol. 1983, 100: 64-119. 10.1016/0012-1606(83)90201-4.
Hunter CP, Kenyon C: Spatial and temporal controls target pal-1 blastomere-specification activity to a single blastomere lineage in C. elegans embryos. Cell. 1996, 87: 217-226. 10.1016/S0092-8674(00)81340-9.
Edgar LG, Carr S, Wang H, Wood WB: Zygotic expression of the caudal homolog pal-1 is required for posterior patterning in Caenorhabditis elegans embryogenesis. Dev Biol. 2001, 229: 71-88. 10.1006/dbio.2000.9977.
Baugh LR, Hill AA, Claggett JM, Hill-Harfe K, Wen JC, Slonim DK, Brown EL, Hunter CP: The homeodomain protein PAL-1 specifies a lineage-specific regulatory network in the C. elegans embryo. Development. 2005, 132: 1843-1845. 10.1242/dev.01782.
Draper BW, Mello CC, Bowerman B, Hardin J, Priess JR: MEX-3 is a KH domain protein that regulates blastomere identity in early C. elegans embryos. Cell. 1996, 87: 205-216. 10.1016/S0092-8674(00)81339-2.
Labouesse M, Hartwieg E, Horvitz HR: The Caenorhabditis elegans LIN-26 protein is required to specify and/or maintain all non-neuronal ectodermal cell fates. Development. 1996, 122: 2579-2588.
Van Auken K, Weaver D, Robertson B, Sundaram M, Saldi T, Edgar L, Elling U, Lee M, Boese Q, Wood WB: Roles of the Homothorax/Meis/Prep homolog UNC-62 and the Exd/Pbx homologs CEH-20 and CEH-40 in C. elegans embryogenesis. Development. 2002, 129: 5255-5268.
Showell C, Binder O, Conlon FL: T-box genes in early embryogenesis. Dev Dyn. 2004, 229: 201-218. 10.1002/dvdy.10480.
Stein LD, Bao Z, Blasiar D, Blumenthal T, Brent MR, Chen N, Chinwalla A, Clarke L, Clee C, Coghlan A, et al: The genome sequence of Caenorhabditis briggsae: a platform for comparative genomics. PLoS Biol. 2003, 1: E45-10.1371/journal.pbio.0000045.
Riechmann JL, Heard J, Martin G, Reuber L, Jiang C, Keddie J, Adam L, Pineda O, Ratcliffe OJ, Samaha RR, et al: Arabidropsis transcription factors: genome-wide comparative analysis among eukaryotes. Science. 2000, 290: 2105-2110. 10.1126/science.290.5499.2105.
Good K, Ciosk R, Nance J, Neves A, Hill RJ, Priess JR: The T-box transcription factors TBX-37 and TBX-38 link GLP-1/Notch signaling to mesoderm induction in C. elegans embryos. Development. 2004, 131: 1967-1978. 10.1242/dev.01088.
Andachi Y: Caenorhabditis elegans T-box genes tbx-9 and tbx-8 are required for formation of hypodermis and body-wall muscle in embryogenesis. Genes Cells. 2004, 9: 331-344. 10.1111/j.1356-9597.2004.00725.x.
Pocock R, Ahringer J, Mitsch M, Maxwell S, Woollard A: A regulatory network of T-box genes and the even-skipped homologue vab-7 controls patterning and morphogenesis in C. elegans. Development. 2004, 131: 2373-2385. 10.1242/dev.01110.
Hartwell LH, Hopfield JJ, Leibler S, Murray AW: From molecular to modular cell biology. Nature. 1999, 402 (Suppl): C47-C52. 10.1038/35011540.
Krause M, Fire A, Harrison SW, Priess J, Weintraub H: CeMyoD accumulation defines the body wall muscle cell fate during C. elegans embryogenesis. Cell. 1990, 63: 907-919. 10.1016/0092-8674(90)90494-Y.
Mathies LD, Henderson ST, Kimble J: The C. elegans Hand gene controls embryogenesis and early gonadogenesis. Development. 2003, 130: 2881-2892. 10.1242/dev.00483.
Hresko MC, Williams BD, Waterston RH: Assembly of body wall muscle and muscle cell attachment structures in Caenorhabditis elegans. J Cell Biol. 1994, 124: 491-506. 10.1083/jcb.124.4.491.
Williams BD, Waterston RH: Genes critical for muscle development and function in Caenorhabditis elegans identified through lethal mutations. J Cell Biol. 1994, 124: 475-490. 10.1083/jcb.124.4.475.
Weintraub H: The MyoD family and myogenesis: redundancy, networks, and thresholds. Cell. 1993, 75: 1241-1244. 10.1016/0092-8674(93)90610-3.
Yun K, Wold B: Skeletal muscle determination and differentiation: story of a core regulatory network and its context. Curr Opin Cell Biol. 1996, 8: 877-889. 10.1016/S0955-0674(96)80091-3.
Molkentin JD, Olson EN: Combinatorial control of muscle development by basic helix-loop-helix and MADS-box transcription factors. Proc Natl Acad Sci USA. 1996, 93: 9366-9373. 10.1073/pnas.93.18.9366.
Maeda I, Kohara Y, Yamamoto M, Sugimoto A: Large-scale analysis of gene function in Caenorhabditis elegans by high-throughput RNAi. Curr Biol. 2001, 11: 171-176. 10.1016/S0960-9822(01)00052-5.
We thank the Kimble lab for providing strain JK3276, and the Caenorhabditis Genetics Center for providing the other strains used in this work. This work was funded by an NIH grant GM64429 to C.P.H.