Identification of genes interacting with rnt-1 through large-scale RNAi screening in Caenorhabditis elegans.

Although many critical roles of the RUNX family proteins have already been identified, little attention has been given to how these proteins interact with other factors. Elucidating RUNX protein interactions will help extend our understanding of their roles in normal development and tumorigenesis. In this study, we performed large-scale RNAi screening to identify genes that genetically interact with rnt-1, the sole homolog of RUNX protein in the nematode Caenorhabditis elegans. To this end, we took advantage of the fact that C. elegans can survive a severe loss of RNT-1 function with only mild phenotypes, and we looked for genes that caused a synthetic phenotype in the rnt-1 mutant background. We identified seven genes, three of which (cdk-8, cic-1, and sur-2) are involved in transcription, two of which (pgp-2 and cct-5) are involved in stress response, and two of which (D2045.7 and W09D10.4) are involved in signaling cascades, according to their functional gene ontology terms. We further confirmed that the CDK8-containing mediator complex genetically interacts with RNT-1 by showing that knockdown of each component of the CDK8 mediator complex caused a synthetic phenotype, that is, the exploded intestine through the vulva (Eiv) phenotype, in the rnt-1 mutant background. We also identified a putative target gene, acs-4, which is regulated by the RNT-1 and CDK8 mediator complex. Our results strengthen the notion that the CDK8 mediator complex may also act together with RUNX proteins in mammals.

The RUNX genes encode a family of transcription factors that contain a conserved DNA-binding Runt domain. Roles of RUNX genes in various aspects of animal development, such as the development of hematopoietic stem cells, bone, and neuronal cells, have been genetically elucidated through extensive studies of mutations in RUNX genes (Blyth ). For instance, translocations on RUNX1 loci were found to occur most commonly in acute leukemia (Osato ; Michaud ; Osato 2004). Mutations of RUNX1 in mice also result in abnormal development of cells in nervous system, endothelial cells, and immune cells (Wang ). Mutations of RUNX2 in humans induce incomplete development of osteoblasts, which results in cleidocranial dysplasia and osteosarcoma (Banerjee ; Werner ). Further studies of mutations in mice identified that RUNX2 is critical for maturation of osteoblasts in osteogenesis (Otto ; Aberg ). Elegant studies in mice revealed that RUNX3 is expressed in multiple tissues, such as endothelial cells in the gastrointestinal tract, T cells, dendritic cells, and neuronal cells, and that RUNX3 has roles in both cell proliferation and differentiation (Bangsow ; Fukamachi and Ito 2004). Interestingly, mutations in RUNX1 and RUNX2 showed the tendency of haploinsufficiency not only in mice but also in humans. Familial platelet disorder, which is caused by haploinsufficiency of RUNX1, leads to acute myelogenous leukemia (Song ). Cleidocranial dysplasia, which is induced by mutations in RUNX2, was found to be the result of a haploinsufficient mutation on RUNX2 (Mundlos ).

It is well-known that a major function of RUNX genes in cells is to regulate the balance between cell proliferation and differentiation (Coffman 2003). Because the decision between cell proliferation and differentiation must reflect the precise cellular environment, the roles and action mechanism of RUNX genes have to be exerted in a context-dependent manner (Braun and Woollard 2009). To reflect precise cellular status, multiple regulations including various splicing forms, overlapped expression of RUNX genes, regulation by upstream activators, numerous RUNX-targeted genes (Cohen 2009), and posttranscriptional regulations (Bae and Lee 2006) were identified for RUNX genes. In addition, RUNX proteins do not act alone, but in conjunction with cofactors; most RUNX proteins form heterodimers with evolutionarily conserved cofactors such as CBFβ. Although rich information on the players that act with RUNX proteins has been obtained in the past, the interacting partners that act transiently, weakly, or indirectly with RUNX proteins have, so far, been overlooked or ignored. Therefore, it would be worth using a genetic model system that allows one to identify genetic interactions of genes by examining the resulting phenotypes.

The nematode Caenorhabditis elegans is a good system for studying genetic interaction of RUNX proteins because there is a sole ortholog of RUNX, (Nam ; Lee ), and the existing mutations of are not lethal. For example, the genetic interaction of with , a cyclin-dependent kinase inhibitor, allowed identification of the role of in cell fate commitment in the hypodermal seam cell division, which is similar to the functions of mammalian RUNXs (Nimmo ). Also, , a cofactor of , interacts with and its mutation or knockdown affects hypoplasia of hypodermal seam cells (Kagoshima ; Shim and Lee 2008). In this study, we screened genes that genetically interact with through a large-scale RNAi screening in the viable mutant allele . We identified seven genes that genetically interact with by observation of the phenotype that is caused by the reduction of both genes’ functions but not by either one alone. We then focused on the components of the CDK8 mediator complex and showed that all of them genetically interact with . Disruption of both and a component of the CDK8 mediator complex showed the "exploded intestine through the vulva" phenotype. We demonstrate that is a putative target gene that is regulated simultaneously by and the CDK8-containing mediator complex.

Materials and Methods

Nematode strains

The Bristol N2 strain was used as the wild-type C. elegans strain. The mutant strain was obtained from the Caenorhabditis Genetics Center (CGC, Minneapolis, MN), and the , , and mutant strains were kind gifts from National Bioresources Project (Tokyo, Japan). The double mutants were generated by using standard C. elegans techniques and confirmed by single-worm PCR to identify deletions. All strains were grown at 20° on standard Nematode growth media plates (Brenner 1974).

Feeding RNAi screening

The RNAi library by J. Ahringer (Cambridge, UK), which covers >80% of C. elegans open reading frames, was used in the screen. We screened genes on chromosomes I and III. All Escherichia coli strains were streaked and cultured with ampicillin on Luria broth. Before worms were fed with E. coli strains, 1 mM IPTG was treated to induce ds RNA transcription of the target genes. The phenotypes were observed in the progeny from the mothers that had been subject to RNAi from their L4 stage. We compared phenotypes of wild-type N2 and animals after each RNAi and further confirmed the phenotypes with other alleles, and , at least three times. HT115 bacteria carrying L4440, the plasmid of the empty vector pPD129.36, were used as a control RNAi.

Microscopy and imaging

To detect RNAi-induced phenotypes, worms were paralyzed with 1 mM levamisole and mounted on 5% agar pads. DIC images were observed using a Zeiss Axioplan 2 microscope.

Identification of candidate target genes from database

To find genes whose knockdown causes the exploded intestine through the vulva phenotype, we identified candidate target genes from the WormMart at WormBase homepage using "exploded through vulva" as keywords (http://www.wormbase.org). To select genes that harbor the RNT-1 binding sites, R(G/T)ACCRCA, in their promoter regions, we identified the numbers and locations of the RNT-1 binding sites in the 2-kb region of the promoters by visual inspection using the CLC sequence viewer 6 (http://www.clcbio.com).

RNA preparation, reverse-transcription, and quantitative real-time PCR analysis

Total RNA of mixed-stage worms was isolated with TRIzol reagent (Invitrogen) by freeze–thaw method given in the standard protocol served with the manufacturer’s manual. cDNA was synthesized with RevertAid M-MuLV reverse-transcriptase (Fermentas Canada) using oligo-dT primer and subjected to PCR amplification. Real-time quantitative PCR was performed with BIO-RAD iQ SYBR Green Supermix in iQ5 as described in the manufacturer’s manual. Primers used in quantitative real-time PCR were generated by Primer3 (http://frodo.wi.mit.edu) with default setting. All the results of quantitative real-time PCR were triplicated in each of three independent samples. The gene was used for normalization. Student t test was used for statistical analyses.

Results and Discussion

The rationale of this study

Although RUNX genes in mammals such as RUNX1 and RUNX2 are known to be haploinsufficient (Blyth ), this seems not to be the case for in C. elegans because even homozygous animals that possess any of the three existing in C. elegans survive and do not show severe phenotypes. It is possible that deletion mutant alleles may not be null mutations, although they delete considerable portions of the coding regions. Larger deletions have not been isolated despite extensive random mutagenesis followed by PCR-based screening (J. H. Ahnn, personal communication), and the phenotype caused by RNAi includes embryonic lethality, which is obviously more severe than that of existing deletion mutations. In addition, there is no chromosomal deficiency available that uncovers the locus, suggesting that in C. elegans may be haploinsufficient and that the three existing deletion mutants are not null, but rather reduction-of-function mutations. Consistently, we were able to detect transcripts RT-PCR even though they were produced as truncated transcripts (Figure 1A, bottom panel). Reduction-of-function mutations are not the best resources for genetics because they make it difficult to determine the exact function of a gene. However, they can sometimes serve as precious tools for elucidating genetic interactions, which otherwise could not be identified because of lethality. We reasoned that the existing deletion mutations of provide a sensitized genetic background so that one could identify interacting genes by looking for synthetic phenotypes.

Figure 1

RNAi screening to identify genes that show a synthetic phenotype with rnt-1. (A) The gene structure of rnt-1 and its deletion mutant alleles. The white bars in the diagram indicate the regions deleted by the rnt-1 mutant alleles. Arrowheads indicate a pair of primers that target full-length of rnt-1 transcripts. The bottom panel shows the rnt-1 cDNA products, which were amplified from wild-type N2, ok351, tm388, and tm491 strains. (B) The scheme of RNAi screening. We performed unbiased RNAi screening with knockdown of genes on chromosomes I and III. We screened for the genes that showed severe morphological defects in the rnt-1 mutant background, but not in the wild-type background. (C) Images of synthetic phenotypes observed in rnt-1 mutant animals after RNAi of each gene listed. Arrows indicate the sites where rupture occurred. Arrowheads indicate dead embryos. *An arrested worm at its larval stage. The scale bar is 200 μm.

Identification of genes that interact with rnt-1 by RNAi screening for synthetic phenotypes

Although the allele has the largest deletion within the coding region, it causes the mildest phenotype among the deletion mutants, that is, a slightly short body length and delayed development. We therefore decided to use allele as our sensitized genetic background for synthetic phenotype screening. To isolate genes that genetically interact with , we performed RNAi screening for genes whose RNAi caused an obvious degenerative phenotype only in the background, but not in the wild-type background (Figure 1B). Among approximately 5000 genes on chromosomes I and III that we screened, we found that seven genes caused an –dependent synthetic phenotype (Figure 1C, summarized in Table 1). The synthetic phenotypes, such as ruptured (Thastrup ), exploded intestine through vulva (Eiv), larval arrest (Lva), and embryonic lethality (Emb), were observed. RNAi of and W09D10.4 showed ruptured phenotype. However, , , and exhibited an Eiv phenotype. Knockdown of caused ruptured phenotype, but Lva also occurred in this case. We also found that RNAi of D2045.7 in the mutant background induced complete embryonic lethality. Among them, , , and can be grouped as components of the transcriptional machinery, and can be grouped as stress-related genes, and D2045.7 and W09D10.4 can be grouped as components of signal cascade, as categorized by the functional GO terms. and are orthologs of the components of the CDK8 mediator complex, CDK8 and CycC, respectively. is the ortholog of MED23, which is a member of the mediator complex. The mediator complex is a core regulatory integrator that links DNA binding factors and transcriptional machinery (Casamassimi and Napoli 2007). These results suggest that the function of RNT-1 as a transcription factor is closely related with specific transcriptional machinery components in C. elegans. was cloned as an ortholog of multidrug resistance gene (Lincke ) and is required for the formation of gut granules and lysosome-related fat storage (Schroeder ). is a component of cytosolic chaperonin (Leroux and Candido 1995). It was reported that knockdown of might increase the expression of through increased oxidative stress evoked by unidentified mechanisms (Kahn ). The fact that responds to environmental stress (Lee ) led us to speculate that a severe defect in the response to stresses may have induced the synthetic phenotype observed. It would be of interest to dissect the molecular aspects of the synthetic phenotypes caused by the reduction of both and genes. In addition, the mechanism of RNT-1 interaction with signaling components such as D2045.7 and W09D10.4 may shed light on the new action mode of RUNX proteins.

Table 1

Genes that were positively selected with severe morphological defects

Sequence Name	Gene Name	LG	Description	Phenotype
C34G6.4	pgp-2	I	ABC transporter family, multidrug resistance family	Rup
F39B2.4	sur-2	I	Mediator complex component, MED23	Eiv
F39H11.3	cdk-8	I	Cyclin-dependent serine/threonine protein kinase, CDK8	Eiv
C07G2.3	cct-5	III	Cytosolic chaperonin	Rup, Lva
D2045.7		III	Serine/threonine protein kinase/TGFβ-stimulated factor	Emb
H14E04.5	cic-1	III	CDK8 kinase–activating protein, Cyclin C	Eiv
W09D10.4		III	Serine/threonine protein phosphatase	Rup

Rup, ruptured; Eiv, exploded intestine through vulva; Lva, larval arrest; Emb, embryonic lethality.

The gene, located on chromosome I, is a cofactor of RNT-1, and double mutations in and are known to cause a similar phenotype as the exploded intestine through the vulva (Kagoshima ). However, because RNAi clone was not included in the Ahringer library, we have not tested the gene in our screening. Likewise, we may have missed additional genes that interact with because we did not screen genes on chromosomes other than chromosomes I and III. Screening of more chromosomes may result in the identification of new players in RNT-1 action.

CDK8 mediator complex genetically interacts with rnt-1

RNAi of and in the mutant background caused a similar phenotype to each other, that is, the Eiv phenotype (Figure 2A). We examined the possibility that other components of the CDK8 mediator complex such as and also genetically interact with . When we knocked down or in the mutant background, we obtained the same Eiv phenotype, but not in the wild-type background (Figure 2A). We tested three deletion mutations, , , and , and we observed similar effects (Figure 2A). Interestingly, RNAi of and greatly increased the incidence of Eiv phenotype compared to RNAi of other components. To further clarify the genetic interaction between and , we generated and examined the double mutants. The double mutants were maintained only as heterozygous for because of the lethality attributable to protruded intestine through the vulva phenotype. We found that 95.4% of the double mutants, which were confirmed by single-worm PCR, exhibited the Eiv phenotype (n = 22; Figure 2C), whereas single-mutant animals and single-mutant animals showed 0% (n >100) and 5% Eiv phenotype (n = 78), respectively.

Figure 2

CDK8 mediator complex genetically interacts with rnt-1. (A) DIC images of animals with the genotype of wild-type, rnt-1(ok351), rnt-1(tm388), or rnt-1(tm491). Each worm was grown in a control or CDK8 mediator complex RNAi plate, and young adult worms were observed. Note that the intestine, not the gonad, was leaked out through the vulva. (B) The percentage of exploded intestine through vulva (Eiv) animals in wild-type and rnt-1 mutant animals treated with RNAi of CDK8 mediator complex components. Total n represents total number of worms examined. (C) DIC images of the vulva region of rnt-1(ok351) cdk-8(tm1238) double mutants. White bar indicates a scale bar (50 μm).

An interesting observation was that the Eiv phenotype was such that the intestine, not the egg-laying apparatus, was spilled through the exploded vulva. One explanation for this phenotype could be that the cell cycles in the hypodermis or vulval cells were deregulated by and defects. To test this possibility, we examined the numbers of vulval precursor cells (VPCs) and of hypodermal seam cells. However, we could not identify significant differences (data not shown), suggesting that there might be other mechanisms that are regulated by both and the CDK8 mediator complex near VPCs and hypodermal seam cells. It would be of interest to examine whether RNT-1 and the CDK8-containing mediator complex are involved in the regulation of signal transduction pathways such as TGFβ and WNT pathways, which may be responsible for the differentiation of hypodermal cells.

acs-4 is a putative target gene of RNT-1 and CDK8 mediator complex

To identify genes that are regulated by RNT-1 and the CDK8 mediator complex, we searched for genes that are reported to cause similar phenotype to our exploded intestine through the vulva phenotype. The keywords to extract the candidate genes from the WormBase were "exploded through vulva." Among them, we selected those with the binding site(s) in their promoter regions (Supporting Information, Table S1). We then performed RNAi experiments of these candidate genes to determine whether their RNAi caused the same phenotype as RNAi in mutants. We also measured the amount of transcripts of putative target genes by real-time quantitative PCR to determine whether their transcription was dependent on RNT-1 and CDK8. We found that was a gene from the list that showed the Eiv phenotype and was coregulated by and (Figure 3A). The transcript of was decreased only when both and were knocked down. We confirmed that the Eiv phenotype was induced by RNAi of alone in the wild-type background (Figure 3, B and C) at a lower frequency than the RNAi in the mutant background, suggesting that is not the only target gene. encodes a fatty acyl CoA synthase and its role has been studied in terms of metabolism. Metabolic problems cause a variety of defects in development. Deprivation of fatty acids as a result of mutation in lipid metabolism causes developmental defects such as embryonic lethality and larval arrest in C. elegans, which are triggered by incomplete cytoskeletal components, reduction of function in signaling pathway, and abnormal biophysical property of membrane (Vrablik and Watts 2012). It has been recently reported that fatty-acid biosynthesis affects the polarity of intestine cells in C. elegans (Zhang ). The intestine of C. elegans stores lots of fatty acids. It is conceivable that deregulation of , which is regulated by and CDK8 mediator complex, may have led to abnormal lipid status of the worms, which has resulted in the Eiv phenotype in general.

Figure 3

acs-4 is a putative target gene of rnt-1 and CDK8 mediator complex. (A) Relative mRNA level of acs-4 detected by quantitative real-time PCR in the mixed stage of wild-type and rnt-1 mutant animals with or without cdk-8 RNAi. It was triplicated and normalized by the act-1 gene. *p value <0.05. NS, not significant. (B) The percentage of exploded intestine through vulva (Eiv) animals when acs-4 was knocked down by RNAi. (C) DIC images of the vulva region of the animals in which acs-4 was knocked down by RNAi. White bar indicate a scale bar (50 μm).

In summary, our results suggest that may physically interact with the CDK8 mediator complex in C. elegans. Consistent with our results, physical and genetic interactions between Lozenge, a RUNX gene, and CDK8 mediator complex have recently been reported in flies (Gobert ). Taken together, we propose that the action mechanism of RUNX family proteins may be evolutionarily conserved, and that it would be of a great interest to establish the role of CDK8 mediator complex in the function of RUNX proteins in mammals.