Large-scale RNAi screens identify novel genes that interact with the C. elegans retinoblastoma pathway as well as splicing-related components with synMuv B activity
BMC Developmental Biology volume 7, Article number: 30 (2007)
The retinoblastoma tumor suppressor (Rb) acts in a conserved pathway that is deregulated in most human cancers. Inactivation of the single Rb-related gene in Caenorhabditis elegans, lin-35, has only limited effects on viability and fertility, yet causes changes in cell-fate and cell-cycle regulation when combined with inactivation of specific other genes. For instance, lin-35 Rb is a synthetic multivulva (synMuv) class B gene, which causes a multivulva phenotype when inactivated simultaneously with a class A or C synMuv gene.
We used the ORFeome RNAi library to identify genes that interact with C. elegans lin-35 Rb and identified 57 genes that showed synthetic or enhanced RNAi phenotypes in lin-35 mutants as compared to rrf-3 and eri-1 RNAi hypersensitive mutants. Based on characterizations of a deletion allele, the synthetic lin-35 interactor zfp-2 was found to suppress RNAi and to cooperate with lin-35 Rb in somatic gonad development. Interestingly, ten splicing-related genes were found to function similar to lin-35 Rb, as synMuv B genes that prevent inappropriate vulval induction. Partial inactivation of specific spliceosome components revealed further similarities with lin-35 Rb functions in cell-cycle control, transgene expression and restricted expression of germline granules.
We identified an extensive series of candidate lin-35 Rb interacting genes and validated zfp-2 as a novel lin-35 synthetic lethal gene. In addition, we observed a novel role for a subset of splicing components in lin-35 Rb-controlled processes. Our data support novel hypotheses about possibilities for anti-cancer therapies and multilevel regulation of gene expression.
The retinoblastoma gene (Rb) was the first tumor suppressor gene to be genetically identified and cloned, based on germline mutations in familial cases of retinoblastoma . Since, deregulation of the Rb pathway has been found to be a general aspect of a wide variety of human cancers, with ~80% of all sporadic cancers containing alterations in Rb or its regulatory components [2, 3]. The frequency of mutations in the Rb gene itself varies greatly between different tumor types, which likely relates to expression levels of other members of the Rb family (p130 and p107) that are closely related in structure and show partly redundant functions in vivo .
Proteins of the pRb family function as transcriptional repressors to control cell division, differentiation and death [, and citations herein]. To control cell proliferation, pRb family members inhibit transcription of S phase promoting genes by binding and blocking activating E2F transcription factors (in mammals: E2F-1, E2F-2 or E2F-3 in association with DP-1 or DP-2). Another group of E2F transcription factors acts together with the pRb protein family to repress transcription, including heterodimers of mammalian E2F-4 or E2F-5 and one of the two DP proteins . Besides E2F family members, pRb-related proteins have been reported to physically interact with a large number of other cellular proteins, including components of chromatin remodeling complexes [e.g.: [7, 8]]. Although the physiological relevance of many interactions remains to be determined, the functions of pRb likely involve a variety of different protein complexes that either block activation of transcription, or actively repress transcription through recruitment of histone-modification and chromatin-remodeling complexes.
The gene lin-35 encodes the single Caenorhabditis elegans homolog of the pRb family. Animals with putative null mutations in lin-35 are viable and fertile , although homozygous mutants show a reduced brood size, delayed development, and larval arrest with low penetrance [[10, 11] and our unpublished observations]. Previous studies revealed redundant functions of lin-35 in a variety of processes such as vulva formation, cell-cycle control, development of the pharynx and somatic gonad, larval development, distal tip cell (DTC) migration and meiotic progression [9–21].
The first and best-characterized function of lin-35 is inhibition of the vulval cell fate [ and references herein]. In this process, lin-35 acts as a synthetic Multivulva (synMuv) class B gene redundantly with class A and C genes . Only animals that contain mutations in two different classes, e.g. A and B, display a Multivulva phenotype (a phenotype that depends on two simultaneous mutations is referred to as "synthetic", hence synMuv). Several synMuv B genes act also in cell-cycle control [12, 22] and many encode homologs of well-known partners of pRb, including efl-1 E2F and dpl-1 DP . Biochemical purification of Rb complexes revealed extensive overlap between components of the Drosophila myb/dREAM complex and C. elegans synMuv B proteins [6, 24]. Together, these data highlight the evolutionary conservation of the Rb pathway as well as the valuable contributions of synMuv B mutants in identifying functional partners of Rb family members.
Since lin-35 is a non-essential gene and its functions are redundant, we decided to perform a synthetic screen to identify additional genes that interact with lin-35 Rb. This screen aimed to identify genes whose functions are essential specifically when lin-35 Rb is inactive. As the Rb pathway is generally compromised in tumor cells, this strategy has the potential of identifying targets for anticancer drugs, because inactivating a gene that is synthetic lethal with loss of Rb activity could preferentially reduce the viability of cancer cells. In a large-scale RNAi approach based on the ORFeome RNAi library , we identified 57 lin-35 interacting genes. In addition, we found a novel group of synMuv B genes that encode components of the splicing machinery. Thus, our screens have revealed many novel genes that may act in parallel or in a linear pathway with lin-35 Rb, some of which may be followed as candidate targets for cancer therapies.
Results and discussion
To expand on previous studies of synthetic interactions with lin-35 Rb, we used RNAi by feeding to inhibit the function of 10,953 genes present in the ORFeome RNAi library v1.1. We compared the effects of feeding RNAi in wild-type (WT) animals and lin-35(n2239) Rb candidate "null" mutants . In the initial screen, we identified 523 genes that showed a variety of synthetic or enhanced phenotypes in lin-35(n2239) as compared to WT (N2) animals (see: Methods; Additional file 1, Supplementary Table 1A; for legend, see Additional file 6).
lin-35 Rb and several other genes that act with Rb in transcriptional repression and chromatin remodeling in other species have been identified in C. elegans based on their inhibitory role in vulval cell-fate specification. Several strains deficient for these so-called synMuv B genes, including lin-35 mutants, were recently found to be RNAi hypersensitive . To exclude candidates that show a non-specific stronger RNAi phenotype in lin-35(n2239), we included the non-Rb related RNAi hypersensitive strain rrf-3(pk1426) in a secondary screen . This assay removed 279 initially selected genes, as they showed similar RNAi phenotypes in rrf-3 and lin-35 mutants. However, for 244 candidate genes, we reproducibly observed a stronger feeding RNAi phenotype in lin-35(n2239) than in rrf-3(pk1426) and WT animals (Supplementary Table 1A). The majority of genes (50/69) for which no RNAi phenotype had been previously reported belonged to this group (according to wormbase version WS153 , see Additional file 1, Supplementary Table 1A).
To further validate the specificity of the interactions with lin-35 Rb, we tested the 244 candidate genes in a more extensive panel of RNAi hypersensitive strains, including rrf-3(pk1426), eri-1(mg366), the synMuv B mutants hpl-2(ok917), lin-15B(n744), lin-35(n2339) as well as lin-37(n758), and double mutant eri-1(mg366); lin-15B(n744) animals. In agreement with Cui et al  but not Wang et al , we observed that lin-37(n758) synMuv B animals are RNAi hypersensitive (Table 1). The strain eri-1; lin-15B was included as an extreme RNAi hypersensitive control, with the caveat that this strain is somewhat sick and was found to display a variety of defects with low penetrance (data not shown).
In the combined screens, we identified 57 genes whose RNAi phenotype was reproducibly stronger in lin-35, or in multiple synMuv B mutants, than in WT, rrf-3 and eri-1. These 57 genes can be classified in three groups: 5 genes showed an RNAi phenotype substantially stronger in lin-35(n2239) than in any other background, 31 genes displayed strong RNAi phenotypes in lin-35(n2239) and in one or two other synMuv B mutants, and 21 genes showed strong RNAi phenotypes in all synMuv B mutants examined (Figure 1) (see Additional file 1, Supplementary Table 1B). We consider all three classes of interest, as synthetic interactions could either be specific for lin-35 or take place with multiple Rb-pathway components. The group of 31 candidates includes pha-1, which was previously shown to interact with lin-35  and provides a first validation that our approach identified specific interactions. In addition, the various classes of genes identified are in agreement with previous genetic studies that revealed redundant functions of lin-35 in various developmental processes [9–21]. A high percentage of the genes identified are predicted to fulfill roles in transcription, DNA replication and chromatin structure, as was expected for genes that act redundantly with lin-35 Rb (Figure 1C–D). Genes involved in protein degradation were also expected to interact with the lin-35 transcriptional repressor. Such genes, as well as genes involved in signaling, cell cycle and RNA metabolism, were all well represented in our screens as in previous screens (Fig. 1C–D, see Additional file 3).
Recently, Lehner et al. described the existence of "hub" genes. This group of genes encodes chromatin regulators that modify diverse signaling pathways and therefore genetically interact with many other genes . As lin-35 is expected to function in chromatin remodeling and interacted with a large number of genes in our screen, we suggest that Rb is one of such "hub" genes.
Importantly, redundant (synthetic) interactions with lin-35 Rb have often been revealed through partial loss-of-function mutations, while complete loss of function of the same genes caused lethality (e.g.: pha-1 , psa-1 , him-17  and gon-14 ; the class A synMuv genes are a clear exception). Consistent with this notion, we injected dsRNA to achieve strong RNAi phenotypes for a sample of 12 candidate genes (including pha-1, the known lin-35 interactor), and all but zfp-2 caused lethality/sterility with or without co-injection of lin-35 dsRNA (see Additional file 1, Supplementary Table 1B). Therefore, synthetic interaction with lin-35 Rb may be detected with partial loss-of-function alleles (or by feeding dsRNA) but not in a setting approaching strong loss-of-function or null conditions (as may be achieved by injecting dsRNA). Thus, independent validation of these genetic interactions will require appropriate alleles.
Similarly, a feeding RNAi phenotype in an eri-1 or rrf-3 mutant background does not rule out genetic interaction with lin-35 Rb. However, genes that show an apparent feeding RNAi phenotype only in lin-35 mutants, and not even in eri-1; lin-15B RNAi hypersensitive mutants, are the best candidates for a functional interaction with lin-35 Rb. Only five genes showed this level of specificity: zfp-2, F25H5.5, rsr-2, Y17G7B.2 and tir-1.
F25H5.5 encodes the putative C. elegans homolog of Claspin, which in vertebrates mediates activation of the Chk1 kinase as part of a DNA replication checkpoint response . As pRb and Claspin both negatively regulate cell-cycle progression, a synthetic lethal interaction is conceivable. It will be attractive to examine whether reducing human Claspin activity specifically increases the killing of Rb-minus tumor cells, if not alone then possibly in response to e.g. gamma-irradiation (experiments in progress).
The predicted Y17G7B.2 product is similar to the trithorax group protein Ash2, a subunit of a histone H3 (Lys4) methyltransferase complex . Genetic interaction between the Ash2 and pRb family members is plausible, as both form part of protein complexes that regulate gene expression through histone modification, and physical interactions between components of these complexes have been reported . Indeed, others recently observed a genetic interaction between mutant alleles of Y17G7B.2/ash-2 and lin-35 (G. Soete and HC Korswagen, manuscript in preparation), which validates the identification of this gene in our screen.
The tir-1 gene encodes a protein with a conserved Toll/IL-1 resistance (TIR) domain and has been shown to act in a conserved MAP kinase pathway in innate immunity [33, 34] as well as lateral signaling between specific neurons . The basis of a genetic interaction between tir-1 and lin-35 Rb is not immediately obvious. However, if confirmed it should be of considerable interest, as TIR-1 contains a single human ortholog, SARM. Like TIR-1, SARM contains a TIR domain, heat-armadillo repeats and two sterile alpha motifs (SAM) .
The zfp-2 gene encodes a zinc-finger transcription factor with six C2H2 zinc finger domains and a KRAB A domain, which is associated with transcriptional repression  (see Additional file 4). Human pRb has been reported to physically interact with two zinc finger proteins named RIZ1 and RBak. RIZ1 is commonly inactivated in human cancers, and RBak contains a KRAB domain and has homology with zfp-2 [37–39]. Finally, rsr-2 encodes a serine (S) and arginine (R)-rich protein with homology to the human splicing co-activator Srm300. The zfp-2 and rsr-2 genes and their interactions with lin-35 are discussed below.
The zinc finger zfp-2 cooperates with lin-35in somatic gonad development
Of the five candidates that interacted specifically with lin-35 (Figure 1A, top), viable mutations were only available for zfp-2 and tir-1. We did not observe a synthetic interaction between lin-35 and the tir-1 deletion alleles ok1052 and tm1111. However, these alleles likely do not cause tir-1 inactivation, as the deletions affect specific splice forms while two of the confirmed tir-1 mRNAs (F13B10.1b and d)  should be expressed normally. In contrast, we observed a strong synthetic interaction between the zfp-2 allele tm557 and the lin-35 candidate null mutations n745 and n2239 as well as lin-35(RNAi). The zfp-2(tm557) mutation is a 233 bp deletion which removes 130 bp of coding sequences and changes the reading frame shortly after the predicted translation-initiation codon.
Inactivation of zfp-2 alone by RNAi (feeding or injection) or the zfp-2(tm557) homozygous deletion did not cause any apparent abnormalities (Figure 2A,C). Although lin-35 mutants are also viable, their brood size is reduced (Figure 2A and ). We observed a variety of defects in lin-35(n2239) mutants with low penetrance, including reduced integrity of the proximal gonad (uterus and spermatheca), a low percentage of endomitotic oocytes (Emo), DTC migration defects (8.4%, n = 107) and abnormalities in vulval morphology (3.7 % everted vulva, n = 107). Loss-of-lin-35 function also allowed continued nuclear division of intestinal nuclei during larval development (variable and partly dependent on developmental stage and temperature, ranging between normal (32–34 nuclei) and 55 nuclei/animal [[17, 19], and our unpublished observations].
Importantly, lin-35(n2239); zfp-2(tm557) double mutants were almost completely sterile (Figure 2A). In addition, defects in gonad migration (18% n = 188) and vulval morphology (17% n = 188) were enhanced compared to lin-35 alone and polyploid nuclei were seen in the soma (Figure 2H). Such sterile double mutant animals showed dramatic morphological defects in the proximal somatic gonad, with absence of a normal uterine cavity and, based on misplaced sperm precursors and spermatids, dysfunctional spermatheca (Figure 2D,2F and 2H). Oocytes were present but in general did not ovulate, were not fertilized and often became polyploid (endomitotic oocytes: Emo phenotype ) (Figure 2H). These defects resemble phenotypes observed after ablation of somatic cell precursors in the germ line, which indicated that sheath and spermathecal cells are required for ovulation and fertilization .
Examination of the zfp-2 expression pattern suggested that the affected tissues express both lin-35 and zfp-2. While lin-35 Rb seems to be ubiquitously expressed [9, 42], zfp-2 expression is more restricted. Transgenic animals expressing the Green Fluorescent Protein (GFP) under the control of the zfp-2 promoter showed fluorescence in vulval cells and all somatic gonad structures such as spermatheca, sheath cells, uterine cells and distal tip cells (DTCs) (Figure 3). This expression pattern correlates well with the observed synthetic phenotype. Together, these observations indicate that lin-35 and zfp-2 function redundantly in development of the somatic gonad lineages derived from Z1 and Z4 (Figure 2B: DTCs, spermatheca, sheath cells and uterine cells). This may explain the abnormal migration of DTCs, morphological defects of the uterus and spermatheca, and synthetic sterile phenotype of lin-35(n2239); zfp-2(tm557) double mutants.
Although gonadal defects were apparent, the developmental program is partially completed and cell specification, based on reporter gene expression and structure, does occur (Figure 2F). Some of the sheath cells were apparently missing in the double mutants, but quantitative analysis was hampered by silencing of the sheath cell marker (lim-7::GFP ) upon injection of lin-35 and zfp-2 dsRNA (Supplementary Figure 3). Knock down of lin-35 and zfp-2 after completion of somatic gonad development did not cause sterility. Specifically, when lin-35 and zfp-2 dsRNA were injected in the gonad of L4 larvae or young adult animals, the injected animals remained fertile while their progeny were sterile. Although such experiments do not provide strong evidence, the results are consistent with redundant functions for lin-35 and zfp-2 in somatic gonad development, rather than oocyte maturation.
F35H8.3/zfp-2 was previously found in an RNAi screen as one of 59 genes required for cosuppression . Cosuppression refers to the coincident silencing of a repetitive transgene and corresponding endogenous gene, which involves chromatin remodeling and RNAi pathway components . Remarkably, we observed that zfp-2(tm557) mutants are not RNAi resistant but hypersensitive to RNAi for certain genes (Table 1). Thus, zfp-2 and lin-35 both suppress RNAi, which indicates a more extensive functional relationship. As yet another similarity with lin-35, zfp-2 RNAi also caused silencing of the rol-6 transgene (see Additional file 5). Based on the similarities and genetic interaction with lin-35, we favor the model that ZFP-2 acts as a zinc-finger transcription factor in transcriptional suppression.
In conclusion, zfp-2 and lin-35 show a synthetic lethal interaction and cooperate in development of the somatic gonad, likely acting in parallel or together in transcription repression and/or chromatin remodeling.
Knock down of specific splicing components resembles the lin-35Rb phenotype
The synMuv B gene lin-35 Rb represses ectopic induction of the vulval cell fate redundantly with synMuv A and C genes. We did not observe synthetic Multivulva phenotypes in our screen, not even for the two known class A genes or single class C gene present in the library. This fits with our observations that class A genes tend to score false negative in feeding RNAi screens (our unpublished results).
Although synthetic interaction between null alleles often points to functions in parallel pathways, several observations indicate that class A and B synMuv proteins could in fact act in conjunction [46, 47]. Similarly, genetic interaction with lin-35 Rb in our screen could reveal genes that act in parallel pathways, but possibly also identify additional components of a LIN-35 protein complex or pathway. As such, it was of interest that a known synMuv B gene, mep-1, was found in the set of 244 candidates. To examine whether additional synMuv B genes were present among the putative lin-35 interactors, we performed feeding RNAi for the 244 candidates in the synMuv A mutant lin-15A(n767). This screen identified rsr-2, one of the five genes that showed a feeding RNAi phenotype specifically in lin-35 mutants (Figure 1A), as a novel synMuv B gene.
Despite the lin-35 specific feeding RNAi phenotype, rsr-2 dsRNA injection caused severe abnormalities such as sterility, embryonic lethality and larval arrest [, and this report]. Interestingly, rsr-2 encodes a serine (S) and arginine (R)-rich protein with homology to the human splicing co-activator Srm300, which acts with other splicing factors to perform its functions. To test whether additional splicing factors belong to the synMuv B family, we collected feeding RNAi clones corresponding to 135 putative C. elegans splicing-related genes (see Additional file 2, Supplementary Table 2A; legend, see Additional file 6). We assayed these 135 genes by feeding RNAi in L1 and young adult lin-15A mutants (see Additional file 2, Supplementary Table 2B). An additional nine genes showed a Muv phenotype in lin-15A and other class A mutants (Table 2), but not in wild-type or class B animals, classifying these genes as synMuv B. These nine genes include the seven members of the Sm family in C. elegans (snr-1, 2, 3, 4, 5, 6 and 7)  as well as two Sm-like genes (lsm-2/gut-2 and lsm-4), all of which are well conserved through evolution. Only lsm-2 and lsm-4 RNAi produced viable animals at 25°C, at which temperature the synMuv phenotype is most highly penetrant. Combination of lsm-2 or lsm-4 RNAi with a mutation in any one of four different Class A genes caused a Muv phenotype, ranging from 20% to 96% of the animals at 25°C (Table 2). Genetic epistasis analysis revealed that this splice factor associated synMuv B phenotype depends on a functional let-60 Ras signaling pathway and is reduced in lin-3 EGF mutants, in agreement with genetic epistasis analysis of other synMuv B genes (Table 2) [9, 50].
In addition to their role in vulval-fate determination, synMuv B genes show specific combinations of other loss-of-function characteristics that include: deregulation of cell-cycle entry, ectopic expression or silencing of transgenes, RNAi hypersensitivity and expression of germline P granules in the soma [26, 51]. Although the ten splicing-related synMuv B genes are essential for viability, partial inactivation by feeding RNAi induced other properties of synMuv B genes (Table 3). For example, like lin-35 inactivation (see above), rsr-2 RNAi led to formation of extra intestinal nuclei (41.3 ± 3.6 (n = 37) versus control: 33.3 ± 1.9 (n = 27) in the elt-2::GFP strain) (Figure 3). In addition, RNAi of several synMuv B genes induces ectopic expression of a lag-2::GFP reporter in the intestine  and we observed the same phenotype after RNAi of lsm-2 and lsm-4, a weaker phenotype after RNAi of snr genes, but no effect after rsr-2 RNAi (Table 3). Moreover, like lin-35 loss of function, feeding RNAi of snr-2 and snr-3 induced ectopic P-granules in the soma (Figure 4). lin-35 RNAi also causes silencing of the scm::GFP transgene and partial rescue of cyd-1 larval arrest. These latter effects were not observed after RNAi of the ten novel synMuv B genes. Together, the observed effects suggest that components of the general splicing machinery contribute surprisingly specific functions to some, but not all, processes mediated by lin-35 Rb and other synMuv B genes (Table 3).
How can spliceosome components behave as synMuv B genes? Evidently, complete inactivation of splicing should cause loss-of-function of all intron-containing genes, including synMuv B genes. However, we used partial loss-of-function conditions that should not completely block splicing, and three of the phenotypic readouts involve increased, rather than reduced, expression levels. Specifically, the synMuv phenotype likely depends on increased expression of lin-3 EGF in the hyp7 hypodermal syncytium , extra nuclear divisions require induction of S phase genes and ectopic lag-2 expression needs de novo expression of the transgene reporter. Importantly, RNAi for many other critical splicing components did not cause a Multivulva phenotype (see Additional file 2, Supplementary Table 2B), indicating that this phenotype is not a consequence of a general shutdown of the splicing process. A more specific defect in splicing of synMuv B pre-mRNAs cannot be excluded. However, it is unclear why splicing of one or more synMuv B messages would depend strongly on a subset of general splicing factors, and not on many others. Moreover, because the splicing related genes display non-overlapping synMuv B features, a single synMuv B pre-mRNA that is particularly sensitive to (alternative) splicing does not provide an explanation for the observed phenotypes (Table 3).
Functions independent from splicing have been reported for Sm and Sm-like (LSM) proteins and should be considered as alternative explanations [49, 52–54]. Studies in C. elegans have shown that Sm proteins are components of germline specific P granules and regulate multiple aspects of the germ-cell fate [49, 52]. While LSM proteins were not examined, these studies indicated that Sm ribonucleoprotein complexes can act distinct from pre-mRNA splicing, possibly in the posttranscriptional control of maternal mRNAs. Recently, cytoplasmic structures known as P bodies or GW bodies have been linked to mRNA degradation and storage . Yeast and human LSM 1–7 proteins have been shown to localize to these cytoplasmic bodies and have been implicated in the degradation of mRNA . Thus, the synMuv B activity of specific splicing related proteins could relate to alternative functions in mRNA metabolism and translation. Such functions could inhibit gene expression post-transcriptionally, in concert with transcriptional repression mediated by LIN-35 Rb complexes. If this model is correct, it could indicate a surprising overlap between repression of gene expression at the transcriptional level in the nucleus and post-transcriptional level in the cytoplasm.
Of great interest, a link between mRNA processing bodies and the microRNA (miRNA) and RNAi pathways has recently been established [55–57]. In addition, RNA-interference related pathways have been observed to cooperate with lin-35 Rb , and Sm proteins have also been found to interact with the RNAi machinery in C. elegans [58, 59]. Specifically: while snr-7 is required for RNAi and transgene expression , SNR-3 physically interacts with DCR-1, a key enzyme in the RNAi and miRNA pathways . Like lin-35 and rsr-2 (this study), dcr-1 also negatively regulates nuclear divisions in the intestine . Thus, cooperation between LIN-35 Rb and Sm/LSm protein functions might involve miRNA/siRNA pathways.
Further studies will be needed to examine potential coupling between transcription, splicing and RNAi-related processes. Based on the overlapping loss-of-function phenotypes, we propose a model in which specific splicing components cooperate with lin-35 synMuv B complexes in repression of gene expression, either by promoting mRNA turnover/translation inhibition or by taking part in multi-protein complexes that regulate transcription, chromatin remodeling, and/or miRNA/RNAi mediated post-transcriptional gene silencing (Figure 5).
Our results provide an extensive series of candidate lin-35 interacting genes, validate zfp-2 as a novel lin-35 synthetic lethal gene, and implicate a subset of spliceosome components in gene regulation in conjunction with the lin-35 Rb pathway.
We used the wild-type strain Bristol N2 and the following mutations:
LG I: lin-35(n745, n2239). LG II: rrf-3(pk1426), zfp-2(tm557), let-23(sy97). LG III: lin-37(n758), hpl-2(ok917). LG IV: eri-1(mg366), let-60(n1876), unc-22(s7), lin-3(n378). LG X: lin-15(n744), lin-15A(n767)
Integrated arrays: rtIs14 [elt-2::GFP; osm-10::HT150Q] IV (a gift from P. W. Faber and A. Hart). ezIs2 [fkh-6::GFP + unc-119(+)] . In addition, we used the integrated array ccIs4251 [myo-3::Ngfp-lacZ, pSAK4(myo-3::Mtgfp)], mapped to LGI (+2.5), to balance lin-35 (n2239) in strain lin-35(n2239)/ccIs4251 I, zfp-2(tm557) II.
RNAi experiments were performed as detailed in Rual et al. . In the screen, dsRNA was delivered by feeding L1 synchronized animals on agar plates containing IPTG (6 mM) plus ampicillin and tetracycline. 10,953 RNAi clones were tested in 24-well plates. In all RNAi experiments, wild-type animals and lin-35 mutants were simultaneously examined in duplicate (thus, 6 different genes were tested per 24-well plate). RNAi phenotypes were scored visually on days 4 and 6 after L1 worms were put to culture. Because of a substantial variation in feeding RNAi phenotypes, the large number of genes tested and range of phenotypes scored, quantitative determination of phenotypes was unpractical. Therefore, candidates were selected based on obvious and reproducibly stronger phenotypes in lin-35 mutants, as compared to other backgrounds described in the text. For all genes of interest, inserts of feeding RNAi vectors were checked by PCR amplification and DNA sequence analysis.
For many genes, a low penetrant RNAi phenotype was observed that became much stronger when combined with inactivation of lin-35 Rb. We refer to such phenotypes as "enhanced", rather than "synthetic", in combination with lin-35. Lethal (Let) is defined as animals unable to produce viable progeny. When possible, a more informative description is used: Emb (embryonic lethal), Lvl (larval lethal), Lva (larval arrest), or Ste (sterile). Other phenotypes scored were: Unc (uncoordinated), Prl (paralyzed), Dpy (dumpy), Bmd (body morphology defective), Sck (sick), Bli (blistered), Mlt (molting defective), Slm (slim), Him (high incidence of males), Pvl (protruding vulva), Muv (multivulva), Clr (clear), Slu (sluggish), Lon (long), Sma (small), Gro (growth rate abnormal), Egl (egg laying defective), Rbs (reduced brood size), Ooc (oocyte formation abnormal), Stp (sterile progeny), and Rup (ruptured). See Additional file 1, Supplementary Table 1, for further information and details.
To screen for synMuv B genes among splicing-related genes, we performed feeding RNAi starting from L4 and L1 stage animals.
Preparation of dsRNA for injection: DNA templates were amplified by PCR from L4440 plasmids using the following primers:
pL440-dest-RNAi-FOR (5'GTTTTCCCAGTCACGACGTT3') and pL440-dest-RNAi-REV (5'TGGATAACCGTATTACCGCC3'). PCR products were used as the templates for synthesis of dsRNA with an in vitro transcription kit (T3 T7 Ambion). The RNA were incubated at 72°C for 10 minutes followed by 30 minutes at 37°C to allow annealing. The quality and amount of dsRNA was examined by electrophoresis. For injection, we used dsRNA at a concentration of 0.5 to 2 μg/μl.
Worm larvae were permeabilized on coverslips using freeze/cracking. Samples were fixed in ice-cold methanol (10 min.) followed by ice-cold acetone (10 min.). K76 mouse monoclonal antibodies were used (1:2) to detect P granules. Rhodamine-conjugated anti-mouse antibodies were used as secondary antibodies (1:200).
For staining of DNA, animals were fixed in Carnoys fixative at 4°C for 24 hours. After rehydration and washes, samples were incubated for 10 minutes in a solution of 1 μg/ml DAPI.
Analysis of zfp-2expression
zfp-2::GFP transcriptional fusions were generated by in vitrorecombination using the promoterome clones and the destination vector pdest-DD03 as described in Dupuy et al. . Clone and primers sequences can be retrieved from Promoterome database .
Friend SH, Bernards R, Rogelj S, Weinberg RA, Rapaport JM, Albert DM, Dryja TP: A human DNA segment with properties of the gene that predisposes to retinoblastoma and osteosarcoma. Nature. 1986, 323 (6089): 643-646. 10.1038/323643a0.
Weinberg RA: The retinoblastoma protein and cell cycle control. Cell. 1995, 81 (3): 323-330. 10.1016/0092-8674(95)90385-2.
Sherr CJ: Cancer cell cycles. Science. 1996, 274 (5293): 1672-1677. 10.1126/science.274.5293.1672.
Classon M, Harlow E: The retinoblastoma tumour suppressor in development and cancer. Nat Rev Cancer. 2002, 2 (12): 910-917. 10.1038/nrc950.
Dimova DK, Stevaux O, Frolov MV, Dyson NJ: Cell cycle-dependent and cell cycle-independent control of transcription by the Drosophila E2F/RB pathway. Genes Dev. 2003, 17 (18): 2308-2320. 10.1101/gad.1116703.
Frolov MV, Dyson NJ: Molecular mechanisms of E2F-dependent activation and pRB-mediated repression. J Cell Sci. 2004, 117 (Pt 11): 2173-81. 10.1242/jcs.01227.
Korenjak M, Taylor-Harding B, Binne UK, Satterlee JS, Stevaux O, Aasland R, White-Cooper H, Dyson N, Brehm A: Native E2F/RBF complexes contain Myb-interacting proteins and repress transcription of developmentally controlled E2F target genes. Cell. 2004, 119 (2): 181-193. 10.1016/j.cell.2004.09.034.
Morris EJ, Dyson NJ: Retinoblastoma protein partners. Adv Cancer Res. 2001, 82: 1-54.
Lu X, Horvitz HR: lin-35 and lin-53, two genes that antagonize a C. elegans Ras pathway, encode proteins similar to Rb and its binding protein RbAp48. Cell. 1998, 95 (7): 981-991. 10.1016/S0092-8674(00)81722-5.
Fay DS, Keenan S, Han M: fzr-1 and lin-35/Rb function redundantly to control cell proliferation in C. elegans as revealed by a nonbiased synthetic screen. Genes Dev. 2002, 16 (4): 503-517. 10.1101/gad.952302.
Cardoso C, Couillault C, Mignon-Ravix C, Millet A, Ewbank JJ, Fontes M, Pujol N: XNP-1/ATR-X acts with RB, HP1 and the NuRD complex during larval development in C. elegans. Dev Biol. 2005, 278 (1): 49-59. 10.1016/j.ydbio.2004.10.014.
Boxem M, van den Heuvel S: lin-35 Rb and cki-1 Cip/Kip cooperate in developmental regulation of G1 progression in C. elegans. Development. 2001, 128 (21): 4349-4359.
Fay DS, Large E, Han M, Darland M: lin-35/Rb and ubc-18, an E2 ubiquitin-conjugating enzyme, function redundantly to control pharyngeal morphogenesis in C. elegans. Development. 2003, 130 (14): 3319-3330. 10.1242/dev.00561.
Fay DS, Qiu X, Large E, Smith CP, Mango S, Johanson BL: The coordinate regulation of pharyngeal development in C. elegans by lin-35/Rb, pha-1, and ubc-18. Dev Biol. 2004, 271 (1): 11-25. 10.1016/j.ydbio.2004.03.022.
Bender AM, Wells O, Fay DS: lin-35/Rb and xnp-1/ATR-X function redundantly to control somatic gonad development in C. elegans. Dev Biol. 2004, 273 (2): 335-349. 10.1016/j.ydbio.2004.06.009.
Cui M, Fay DS, Han M: lin-35/Rb cooperates with the SWI/SNF complex to control Caenorhabditis elegans larval development. Genetics. 2004, 167 (3): 1177-1185. 10.1534/genetics.103.024554.
Saito RM, Perreault A, Peach B, Satterlee JS, van den Heuvel S: The CDC-14 phosphatase controls developmental cell-cycle arrest in C. elegans. Nat Cell Biol. 2004, 6 (8): 777-783. 10.1038/ncb1154.
Reddy KC, Villeneuve AM: C. elegans HIM-17 links chromatin modification and competence for initiation of meiotic recombination. Cell. 2004, 118 (4): 439-452. 10.1016/j.cell.2004.07.026.
Grishok A, Sharp PA: Negative regulation of nuclear divisions in Caenorhabditis elegans by retinoblastoma and RNA interference-related genes. Proc Natl Acad Sci USA. 2005, 102 (48): 17360-17365. 10.1073/pnas.0508989102.
Chesney MA, Kidd AR, Kimble J: gon-14 functions with class B and class C synthetic multivulva genes to control larval growth in Caenorhabditis elegans. Genetics. 2006, 172 (2): 915-928. 10.1534/genetics.105.048751.
Ceol CJ, Horvitz HR: A new class of C. elegans synMuv genes implicates a Tip60/NuA4-like HAT complex as a negative regulator of Ras signaling. Dev Cell. 2004, 6 (4): 563-576. 10.1016/S1534-5807(04)00065-6.
Boxem M, van den Heuvel S: C. elegans class B synthetic multivulva genes act in G(1) regulation. Curr Biol. 2002, 12 (11): 906-911. 10.1016/S0960-9822(02)00844-8.
Ceol CJ, Horvitz HR: dpl-1 DP and efl-1 E2F act with lin-35 Rb to antagonize Ras signaling in C. elegans vulval development. Mol Cell. 2001, 7 (3): 461-73. 10.1016/S1097-2765(01)00194-0.
Lewis PW, Beall EL, Fleischer TC, Georlette D, Link AJ, Botchan MR: Identification of a Drosophila Myb-E2F2/RBF transcriptional repressor complex. Genes Dev. 2004, 18 (23): 2929-40. 10.1101/gad.1255204.
Rual JF, Ceron J, Koreth J, Hao T, Nicot AS, Hirozane-Kishikawa T, Vandenhaute J, Orkin SH, Hill DE, van den Heuvel S, et al: Toward improving Caenorhabditis elegans phenome mapping with an ORFeome-based RNAi library. Genome Res. 2004, 14 (10B): 2162-2168. 10.1101/gr.2505604.
Wang D, Kennedy S, Conte D, Kim JK, Gabel HW, Kamath RS, Mello CC, Ruvkun G: Somatic misexpression of germline P granules and enhanced RNA interference in retinoblastoma pathway mutants. Nature. 2005, 436 (7050): 593-597. 10.1038/nature04010.
Simmer F, Tijsterman M, Parrish S, Koushika SP, Nonet ML, Fire A, Ahringer J, Plasterk RH: Loss of the putative RNA-directed RNA polymerase RRF-3 makes C. elegans hypersensitive to RNAi. Curr Biol. 2002, 12 (15): 1317-1319. 10.1016/S0960-9822(02)01041-2.
Cui M, Kim EB, Han M: Diverse chromatin remodeling genes antagonize the RB-involved SynMuv pathways in C. elegans. Plos Genet. 2006, 2 (5): e74-10.1371/journal.pgen.0020074.
Lehner B, Crombie C, Tischler J, Fortunato A, Fraser AG: Systematic mapping of genetic interactions in Caenorhabditis elegans identifies common modifiers of diverse signaling pathways. Nat Genet. 2006, 38 (8): 896-903. 10.1038/ng1844.
Freire Raimundo, van Vugt Marcel, Mamely Ivan, Medema Rene: Claspin: Timing the cell cycle arrest when the genome is damaged. Cell Cycle. 2006, 5 (24): 2831-2834.
Ringrose Leonie, Paro Renato: Epigenetic regulation of cellular memory by the Polycomb and Trithorax group proteins. Annu Rev Genet. 2004, 38: 413-43. 10.1146/annurev.genet.38.072902.091907.
Wysocka J, Myers MP, Laherty CD, Eisenman RN, Herr W: Human Sin3 deacetylase and trithorax-related Set1/Ash2 histone H3-K4 methyltransferase are tethered together selectively by the cell-proliferation factor HCF-1. Genes Dev. 2003, 17 (7): 896-911. 10.1101/gad.252103.
Couillault Carole, Pujol Nathalie, Reboul Jérôme, Sabatier Laurence, Guichou Jean-François, Kohara Yuji, Ewbank Jonathan: TLR-independent control of innate immunity in Caenorhabditis elegans by the TIR domain adaptor protein TIR-1, an ortholog of human SARM. Nature Immunology. 2004, 5: 488-494. 10.1038/ni1060.
Liberati Nicole, Fitzgerald Katherine, Kim Dennis, Feinbaum Rhonda, Golenbock Douglas, Ausubel Frederick: Requirement for a conserved Toll/interleukin-1 resistance domain protein in the Caenorhabditis elegans immune response. Proc Natl Acad Sci USA. 2004, 101 (17): 6593-6598. 10.1073/pnas.0308625101.
Chuang Chiou-Fen, Bargmann Cornelia: A Toll-interleukin 1 repeat protein at the synapse specifies asymmetric odorant receptor expression via ASK1 MAPKKK signaling. Genes & Dev. 2005, 19: 270-281. 10.1101/gad.1276505.
Urrutia R: KRAB-containing zinc-finger repressor proteins. Genome Biol. 2003, 4 (10): 231-10.1186/gb-2003-4-10-231.
Buyse IM, Shao G, Huang S: The retinoblastoma protein binds to RIZ, a zinc-finger protein that shares an epitope with the adenovirus E1A protein. Proc Natl Acad Sci USA. 1995, 92 (10): 4467-4471. 10.1073/pnas.92.10.4467.
Canote R, Du Y, Carling T, Tian F, Peng Z, Huang S: The tumor suppressor gene RIZ in cancer gene therapy (review). Oncol Rep. 2002, 9 (1): 57-60.
Skapek SX, Jansen D, Wei TF, McDermott T, Huang W, Olson EN, Lee EY: Cloning and characterization of a novel Kruppel-associated box family transcriptional repressor that interacts with the retinoblastoma gene product, RB. J Biol Chem. 2000, 275 (10): 7212-7223. 10.1074/jbc.275.10.7212.
Greenstein D, Hird S, Plasterk RH, Andachi Y, Kohara Y, Wang B, Finney M, Ruvkun G: Targeted mutations in the Caenorhabditis elegans POU homeo box gene ceh-18 cause defects in oocyte cell cycle arrest, gonad migration, and epidermal differentiation. Genes Dev. 1994, 8 (16): 1935-1948.
McCarter J, Bartlett B, Dang T, Schedl T: Soma-Germ Cell Interactions in Caenorhabditis elegans: Multiple events of hermaphrodite germline development require the somatic sheat and spermathecal lineages. Dev Biol. 1997, 181: 121-143. 10.1006/dbio.1996.8429.
Myers TR, Greenwald I: lin-35 Rb acts in the major hypodermis to oppose ras-mediated vulval induction in C. elegans. Dev Cell. 2005, 8 (1): 117-123. 10.1016/j.devcel.2004.11.015.
Kimble J, Hirsh D: The postembryonic cell lineages of the hermaphrodite and male gonads in Caenorhabditis elegans. Dev Biol. 1979, 70 (2): 396-417. 10.1016/0012-1606(79)90035-6.
Hall DH, Winfrey VP, Blaeuer G, Hoffman LH, Furuta T, Rose KL, Hobert O, Greenstein D: Ultrastructural features of the adult hermaphrodite gonad of Caenorhabditis elegans: relations between the germ line and soma. Dev Biol. 1999, 212 (1): 101-123. 10.1006/dbio.1999.9356.
Robert VJ, Sijen T, van Wolfswinkel J, Plasterk RH: Chromatin and RNAi factors protect the C. elegans germline against repetitive sequences. Genes Dev. 2005, 19 (7): 782-787. 10.1101/gad.332305.
Walhout AJ, Sordella R, Lu X, Hartley JL, Temple GF, Brasch MA, Thierry-Mieg N, Vidal M: Protein interaction mapping in C. elegans using proteins involved in vulval development. Science. 2000, 287 (5450): 116-22. 10.1126/science.287.5450.116.
Davison EM, Harrison MM, Walhout AJ, Vidal M, Horvitz HR: lin-8, which antagonizes Caenorhabditis elegans Ras-mediated vulval induction, encodes a novel nuclear protein that interacts with the LIN-35 Rb protein. Genetics. 2005, 171 (3): 1017-31. 10.1534/genetics.104.034173.
Longman D, McGarvey T, McCracken S, Johnstone IL, Blencowe BJ, Caceres JF: Multiple interactions between SRm160 and SR family proteins in enhancer-dependent splicing and development of C. elegans. Curr Biol. 2001, 11 (24): 1923-1933. 10.1016/S0960-9822(01)00589-9.
Barbee SA, Lublin AL, Evans TC: A novel function for the Sm proteins in germ granule localization during C. elegans embryogenesis. Curr Biol. 2002, 12 (17): 1502-1506. 10.1016/S0960-9822(02)01111-9.
Cui M, Chen J, Myers TR, Hwang BJ, Sternberg PW, Greenwald I, Han M: SynMuv genes redundantly inhibit lin-3/EGF expression to prevent inappropriate vulval induction in C. elegans. Dev Cell. 2006, 10 (5): 667-672. 10.1016/j.devcel.2006.04.001.
Poulin G, Dong Y, Fraser AG, Hopper NA, Ahringer J: Chromatin regulation and sumoylation in the inhibition of Ras-induced vulval development in Caenorhabditis elegans. Embo J. 2005, 24 (14): 2613-2623. 10.1038/sj.emboj.7600726.
Barbee SA, Evans TC: The Sm proteins regulate germ cell specification during early C. elegans embryogenesis. Dev Biol. 2006, 291 (1): 132-143. 10.1016/j.ydbio.2005.12.011.
Shet U, Parker R: Decapping and decay of messenger RNA occur in cytoplasmic processing bodies. Science. 2003, 300 (5620): 805-8. 10.1126/science.1082320.
Ingelfinger D, Arndt-Jovin DJ, Luhrmann R, Achsel T: The human LSm1-7 proteins colocalize with the mRNA-degrading enzymes Dcp1/2 and Xrnl in distinct cytoplasmic foci. Rna. 2002, 8 (12): 1489-1501.
Ding L, Spencer A, Morita K, Han M: The developmental timing regulator AIN-1 interacts with miRISCs and may target the argonaute protein ALG-1 to cytoplasmic P bodies in C. elegans. Mol Cell. 2005, 19 (4): 437-47. 10.1016/j.molcel.2005.07.013.
Liu J, Rivas FV, Wohlschlegel J, Yates JR, Parker R, Hannon GJ: A role for the P-body component GW182 in microRNA function. Nat Cell Biol. 2005, 7 (12): 1261-6. 10.1038/ncb1333.
Pillai RS, Bhattacharyya SN, Artus CG, Zoller T, Cougot N, Basyuk E, Bertrand E, Filipowicz W: Inhibition of translational initiation by Let-7 MicroRNA in human cells. Science. 2005, 309 (5740): 1573-6. 10.1126/science.1115079.
Kim JK, Gabel HW, Kamath RS, Tewari M, Pasquinelli A, Rual JF, Kennedy S, Dybbs M, Bertin N, Kaplan JM, et al: Functional genomic analysis of RNA interference in C. elegans. Science. 2005, 308 (5725): 1164-1167. 10.1126/science.1109267.
Duchaine TF, Wohlschlegel JA, Kennedy S, Bei Y, Conte D, Pang K, Brownell DR, Harding S, Mitani S, Ruvkun G, et al: Functional proteomics reveals the biochemical niche of C. elegans DCR-1 in multiple small-RNA-mediated pathways. Cell. 2006, 124 (2): 343-354. 10.1016/j.cell.2005.11.036.
Chang W, Tilmann C, Thoemke K, Markussen FH, Mathies LD, Kimble J, Zarkower D: A forkhead protein controls sexual identity of the C. elegans male somatic gonad. Development. 2004, 131 (6): 1425-1436. 10.1242/dev.01012.
Dupuy D, Li Q, Deplancke B, Boxem M, Hao T, Lamesch P, Sequerra R, Bosak S, Doucette-Stamm L, IA H, et al: A First version of the Caenorhabditis elegans Promoterome. Genome Res. 2004, 14: 2169-2175. 10.1101/gr.2497604.
Kennedy S, Wang D, Ruvkun G: A conserved siRNA-degrading RNase negatively regulates RNA interference in C. elegans. Nature. 2004, 427 (6975): 645-649. 10.1038/nature02302.
Simmer F, Moorman C, van der Linden AM, Kuijk E, van den Berghe PV, Kamath RS, Fraser AG, Ahringer J, Plasterk RH: Genome-wide RNAi of C. elegans using the hypersensitive rrf-3 strain reveals novel gene functions. PLoS Biol. 2003, 1 (1): E12-10.1371/journal.pbio.0000012.
Kamath RS, Fraser AG, Dong Y, Poulin G, Durbin R, Gotta M, Kanapin A, Le Bot N, Moreno S, Sohrmann M, et al: Systematic functional analysis of the Caenorhabditis elegans genome using RNAi. Nature. 2003, 421 (6920): 231-237. 10.1038/nature01278.
Lehner B, Calixto A, Crombie C, Tischler J, Fortunato A, Chalfie M, Fraser A: Loss of LIN-35, the Caenorhabditis elegans ortholog of the tumor suppressor p105Rb, results in enhanced RNA interference. Genome Biol. 2006, 7 (1): R4-10.1186/gb-2006-7-1-r4.
The biology and geneome of C. elegans database. [http://www.wormbase.org]
Promoterome database. [http://vidal.dfci.harvard.edu/promoteromedb/]
We thank Mike Boxem and Erick Morris for critical reading of the manuscript, Kursheed Wani for technical assistance, John Satterlee and members of the Hart and van den Heuvel laboratories for helpful discussions. We are grateful to the Hart, Ruvkun and Horvitz laboratories, the Caenorhabditis Genetics Center and the National Bioresource Project for C. elegans (Japan) for strains and information. This research was supported by "AstraZeneca's Academic Sponsorship Program" and the National Institutes of Health (P01 CA95281). JC was a recipient of a postdoctoral fellowship from the Spanish Government (MEC).
JC conducted most of the experiments. JFR and MV generated and provided the RNAi library. AC participated in some of the RNAi screens. DD generated the zfp-2 expression vector. SV conceived the study, and participated in its design and coordination. JC and SV wrote the manuscript. All authors have read and approved the final manuscript.
Electronic supplementary material
Additional File 1: Supplementary Table 1. Results from RNAi screens and functional categories of candidates genes. (additional legend : Additional file 6) (XLS 242 KB)
About this article
Cite this article
Ceron, J., Rual, JF., Chandra, A. et al. Large-scale RNAi screens identify novel genes that interact with the C. elegans retinoblastoma pathway as well as splicing-related components with synMuv B activity. BMC Dev Biol 7, 30 (2007). https://doi.org/10.1186/1471-213X-7-30