|
|
|
RESEARCH PAPER
1 Program in Biochemistry and Molecular Biology and 2 Department of Cellular and Molecular Pharmacology, University of California, San Francisco, California 94143-2280, USA; 3 Huffington Center on Aging, Baylor College of Medicine, Houston, Texas 77030, USA
| Abstract |
|---|
|
|
|---|
[Keywords: DAF-12; In Vitro Genomic Selection; intracellular receptor; response element; target genes]
Received May 6, 2004; revised version accepted August 23, 2004.
The nematode C. elegans develops in 2 d from embryos, through four larval stages L1L4, to adults. During the 23-wk C. elegans life cycle, IR DAF-12 (Antebi et al. 2000
; Snow and Larsen 2000
) functions as a developmental regulator with at least three broad phenotypic effects. First, in response to unfavorable environmental conditions, including lack of food, overcrowding, and elevated temperature, DAF-12 drives formation of dauer larvae (Riddle et al. 1981
). Dauers are long-lived, nonfeeding, stress-resistant specialized L3 larvae that resume development on return to favorable conditions (Riddle and Albert 1997
). Second, DAF-12 functions as a heterochronic gene to regulate C. elegans developmental age (Antebi et al. 1998
). Finally, DAF-12 affects C. elegans adult lifespan (Larsen et al. 1995
; Gems et al. 1998
; Hsin and Kenyon 1999
; Gerisch et al. 2001
; Jia et al. 2002
).
Multiple alleles of daf-12 with distinct protein sequence alterations partially uncouple its phenotypic effects (Antebi et al. 2000
). Genetic analyses place DAF-12, as dauer regulator, at the convergence of TGF-
(transforming growth factor
), cGMP (cycle guanosine 3'-5'-monophosphate), and insulin/IGF-I (insulin growth factor I)-like signaling pathways (Riddle and Albert 1997
). This implies that DAF-12 and other C. elegans receptors, like mammalian IRs, control broad and diverse gene networks (Sengupta et al. 1994
; Larsen et al. 1995
; Antebi et al. 1998
; Carmi et al. 1998
; Kostrouchova et al. 1998
; Hsin and Kenyon 1999
; Asahina et al. 2000
; Gissendanner and Sluder 2000
; Much et al. 2000
; Lindblom et al. 2001
). However, as in mammals, the downstream effectors and the mechanistic basis for the multiple outputs of DAF-12 are unknown. Moreover, unlike mammals, little is known about the molecular biology and molecular physiology of C. elegans IR action, that is, ligands, genomic binding sites, response elements, or target genes have not been identified. Thus, to exploit the advantages of C. elegans toward deciphering DAF-12 gene networks, we must first investigate the molecular biology of DAF-12 action.
As one approach, random selection and amplification methodologies (Blackwell and Weintraub 1990
) might be used to identify high-affinity binding sites for DAF-12. However, such sequences would not necessarily represent bona fide response elements, nor would they identify genes regulated by DAF-12. In this study, we sought to develop an approach in which, starting with the C. elegans genome, we could search biochemically for genomic fragments bound by DAF-12, identify binding sites, test their activities as response elements, and identify DAF-12 target genes.
| Results |
|---|
|
|
|---|
Drawing from concepts of exponential enrichment methodology (Blackwell and Weintraub 1990
; Ellington and Szostak 1990
; Tuerk and Gold 1990
; Shtatland et al. 2000
), we developed an in vitro selection and amplification method, In Vitro Genomic Selection, in which immobilized recombinant DAF-12 DNA-binding domain (DAF-12 DBD) was used to screen total C. elegans genomic DNA digested with restriction endonuclease Sau3AI. The fragments that emerge contain DNA sequences bound directly by DAF-12, and in contrast to synthetic optimized sequences, the binding sites reside in their native C. elegans genomic context. This method is distinct from those in which the bound fragments are incorporated into libraries (Sompayrac and Danna 1990
) or those that use indirect readouts of DNA:protein interactions (Wilson et al. 1991
). Kinzler and Vogelstein (1989
) described Whole Genome PCR to identify genomic fragments containing previously known transcription factor-binding sites. We used In Vitro Genomic Selection to identify unknown DAF-12-binding sites, response elements, and DAF-12-regulated target genes.
After four rounds of selection and PCR amplification (Fig. 1AE), we isolated 47 clones containing C. elegans genomic DNA fragments that selectively bound to DAF-12; these clones represented 26 distinct genomic fragments (Table 1; Supplementary Table S1). Ten of the 26 fragments were overrepresented; three fragments denoted 4.1, 4.2, and 4.3 were recovered five, seven, and six times, respectively; one fragment was recovered three times, and six fragments were recovered twice (Supplementary Table S1).
|
|
|
We compared by quantitative real-time (qRT) PCR the relative abundance of 15 recovered genomic fragments (from the fourth, third, and second rounds of selection, normalized to total genomic DNA) and six neighboring nonspecific control genomic fragments, revealing up to 105-fold enrichment of specific fragments, such as 4.1 and 4.2, and up to 103-fold loss of neighboring nonspecific fragments, including 4.1control and 4.2control (Fig. 1F). The 4.1 and 4.2 fragments were enriched more than 107-fold over the corresponding nonspecific control fragments. During successive rounds of selection and amplification, six of 10 qRT PCR-analyzed fragments recovered in the third and second rounds of the selection showed decreased enrichment after the second or third round; after the fourth round, the abundance levels of 4 fragments, 2.4, 2.14, 2.27, and 3.3, fell below the starting levels (setting as 1 the starting abundance of each fragment of total C. elegans genomic DNA), suggesting that those fragments are likely nonspecific. In contrast, each of the five examined fragments from the fourth round (4.1, 4.2, 4.6, 4.13, and 4.19) was enriched throughout the procedure (Fig. 1F). We conclude that during In Vitro Genomic Selection, particular C. elegans fragments, likely those with the highest relative affinity for DAF-12, were progressively enriched with successive rounds of selection and amplification, whereas nonspecific genomic fragments were lost.
Selected C. elegans genomic fragments display DAF-12 response element activity in yeast
We carried out a functional assay in Saccharomyces cerevisiae to determine if the three most overrepresented C. elegans DNA fragments isolated by In Vitro Genomic Selection conferred DAF-12 response element activity. Although yeast lacks IRs and metazoan-specific cofactors, it has proved useful as a simplified physiological setting for studies of mammalian IR activities (Schena and Yamamoto 1988
). Curiously, response elements that confer negative or positive regulation in a given mammalian context typically display positive regulatory activity in yeast; thus, yeast identifies many functional response elements but does not reliably indicate the direction of regulation.
We found that full-length DAF-12 weakly activated transcription from yeast promoter-reporter constructs bearing linked 4.1 or 4.2 sequences and failed to activate from the 4.3 sequence (Fig. 2AC). Such weak activity in yeast might be expected if DAF-12 action in C. elegans is normally dependent on binding a hormonal ligand. Indeed, mammalian glucocorticoid and androgen receptors are nearly inactive in the absence of their cognate ligands; however, derivatives of these receptors lacking their C-terminal ligand-binding domains (LBDs) behave as constitutive activators (Godowski et al. 1987
; Rundlett et al. 1990
). Correspondingly, we found that a truncated DAF-12 lacking its putative LBD, DAF-12 N500, activated transcription robustly from all three reporter constructs (Fig. 2DF); as a control, the DAF-12 DBD alone was inactive. Thus, these experiments show that 4.1, 4.2, and 4.3 fragments can serve as DAF-12 response elements in yeast and suggest that DAF-12 function may be ligand dependent.
|
DR5 elements function as DAF-12-binding sites and response elements
We next conducted in vitro binding experiments using fluorescence anisotropy with labeled oligonucleotides. DAF-12 bound with high affinity to an oligonucleotide bearing the DR5 sequence from 4.2 genomic fragment (DR5gt, Fig. 3A). As DR5gt contained two distinct hexamers (AGGACA and AGTGCA), we investigated their relative affinities. A synthetic DR5 element bearing two AGTGCA hexamers (DR5tt) bound to DAF-12 with a higher affinity than the natural DR5gt element (Fig. 3A), suggesting that AGTGCA might serve as a high-affinity binding half-site for DAF-12.
|
We next determined whether the DR5-bearing oligonucleotides were sufficient to mediate DAF-12 N500-dependent activation in yeast. As shown in Figure 3D, a single DR5gt sequence was indeed competent to mediate a DAF-12 response. Furthermore, mutant elements with decreased DAF-12 binding affinity (Fig. 3B) conferred little or no transcriptional activity (Fig. 3D), whereas mutations outside the DAF-12 recognition sequences (DR5 m0) did not affect activity. Finally, 5-bp spacing between half-sites in a direct repeat orientation displayed strong response element activity, and increased spacing (to seven bases, DR5 + sp2, Fig. 3C), decreased spacing (DR04, Fig. 3C), or an inverted repeat configuration with various spacing (IR05, 7, Fig. 3C) reduced or eliminated the activity (Fig. 3D,E). Interestingly, DAF-12 N500 activated at moderate levels from DR4, DR3, and IR2 elements (Fig. 3E). Indeed, one of the recovered fragments, 4.7, which was active in yeast (Fig. 2G), contained an IR2 element (Table 2). Thus, DR5 is sufficient and an element containing two half-sites appears to be necessary for DAF-12 response element activity in yeast.
We found that two other C. elegans IRs, CHR3 (Kostrouchova et al. 1998
) and NHR-25 (Gissendanner and Sluder 2000
), failed to activate transcription from the DR5gt DAF-12 response element, whereas DAF-12 was inactive at response elements (identified by M. Van Gilst, unpubl.) specific to CHR3 or NHR-25; DAF-12 also failed to activate transcription from a vertebrate IR response element-containing consensus AGGTCA half-sites (data not shown). Furthermore, transcriptional activation from the 4.2 fragment required functional DAF-12, as a DAF-12 DBD mutant (Antebi et al. 2000
) that is inactive in C. elegans was similarly devoid of activity in the yeast reporter assay (data not shown). Thus, although further study will likely reveal context effects that complicate this picture (Yamamoto et al. 1998
), it is apparent that DAF-12 action is selective and response element specific.
Multiple DAF-12-binding half-sites contribute to overall response element activity from the selected C. elegans genomic fragments
We found that mutations of "solo" half-sites (those not in the DR5 configuration; Supplementary Fig. S2A) within 4.2 modestly reduced DAF-12-mediated activation in yeast and that combinations of mutated half-sites greatly decreased or abolished activity (Supplementary Fig. S2BD). Double and triple mutations within AGGACA half-sites (m01, m41, m51; Supplementary Fig. S2D) had less pronounced effects than did double mutants within the AGTGCA and AGTGCG elements (m1, m21, m31), consistent with our in vitro binding competition assay results. Similar results were obtained with the 4.1 fragment, where mutations in a solo half-site (m21) or in a DR5 element (m11) decreased slightly the activation by DAF-12 but in combination abolished activation (Supplementary Fig. S2EF). We conclude that the 4.2 and 4.1 genomic fragments carry composite response elements, each bearing multiple functional hexameric DAF-12-binding sites; it is likely that In Vitro Genomic Selection biases toward recovery of DNA fragments carrying multiple binding sites for the selecting protein. Such complex organization seems to be characteristic of many normal genomic enhancers (Payvar et al. 1982
; Huang et al. 1999
).
DAF-12 regulates LIT-1 gene expression in C. elegans
We next tested whether DAF-12 could activate transcription in C. elegans from a green fluorescent protein (GFP) reporter gene linked to 4.2. As seen in Figure 4A, the microinjected reporter was strongly expressed in the larval and adult pharynx of wild-type N2 C. elegans, whereas expression was dramatically reduced in an isogenic strain (rh61rh412) carrying a loss-of-function daf-12 null-like allele (Antebi et al. 2000
) (Fig. 4C). Furthermore, injection of a reporter bearing the 4.2 m11m31 mutant element into N2 C. elegans yielded little GFP expression (Fig. 4B). Thus, a genomic DNA fragment bound specifically by DAF-12 in vitro displayed DAF-12 response element activity in vivo both in yeast and in C. elegans. Moreover, DAF-12-activated transcription from the 4.2 sequences in C. elegans was tissue specific.
|
Strikingly, the LIT-1GFP construct with the 4.2 m11m31 mutant response element produced increased GFP expression in seam and vulva cells at larval stage 4 (L4) in wild-type animals (Fig. 4H,K), compared with the wild-type 4.2 construct (Fig. 4G,J). We also observed elevated expression when the wild-type LIT-1GFP construct was microinjected into the daf-12 (rh61rh412) mutant strain (Fig. 4I,L). Thus, DAF-12 may repress LIT-1 expression in seam and vulva cells in L4 animals.
We next used qRT PCR to measure endogenous LIT-1 mRNA levels in whole animals at L4. We observed a greater than twofold increased accumulation of both LIT-1 mRNA isoforms, a and b, in daf-12 (m583) mutant C. elegans (bearing a daf-12 null-like allele; Antebi et al. 2000
), compared with N2 worms (Fig. 4M). These results strongly suggest that LIT-1 is a bona fide DAF-12 target gene.
Recovered C. elegans fragments are linked to DAF-12-regulated genes and gene clusters
Given our findings with 4.2, we determined whether other isolated fragments were linked to DAF-12-regulated genes. We measured by qRT PCR the mRNA levels of 75 C. elegans genes proximal to 28 recovered genomic fragments that exhibited DAF-12-dependent transcriptional activity in yeast. We examined relative gene expression in wild-type N2 and daf-12 (m583) mutant C. elegans at the L4 stage, as LIT-1 gene regulation and aspects of DAF-12 function, for example, proper gonadal migration, occur in L4 (Antebi et al. 1998
, 2000
). We found (Fig. 5A,B) that the expression of various genes (including GST family member gst-37, Skp1 ubiquitin ligase-like skr-8 gene, Snf7 homolog-like C37C3.3, hedgehog-like grd-8, Notch pathway factor sel-8a/lag-3a, and C. elegans caveolin cav-1) was decreased by greater than twofold in the daf-12 mutant, whereas the expression of others (including furrowed-like T07H6.5, soluble guanylyl cyclase gcy-36, putative mitochondrial protein C37C3.8b, prion-like Q/N rich pqn-89, intermediate filament ifa-2, adt-2 protease, phosphocholine cytidylyltransferase-like F08C6.2, cuticle collagen-like col-130, UDP-glucuronosyltransferase-like F08G5.5, and ornithine decarboxylase odc-1) was increased by greater than twofold in the daf-12 mutant (see Supplemental Materials for gene references). The expression of various control genes, including CAP-1, CAP-2, RGR-1, ACT-5, RPL-19, DPY-26, and others, did not change (Fig. 4M; data not shown).
|
The 3.16 fragment mapped to an intron of the K11C4.2 ORF. In L4, we found that although K11C4.2 did not display DAF-12-dependent expression, the nearby K11C4.1 (serine/threonine kinase-like) gene transcript was down-regulated in daf-12 mutant C. elegans, whereas the odc-1 message was up-regulated (Fig. 5B, center panel).
The 2.11 C. elegans genomic fragment mapped to two locations 110 kb apart on chromosome V, an intron/coding region junction in the nearly identical Y19D10A.5 and C01B4.8 ORFs (Supplementary Table S3), which themselves reside within two regions of high similarity spanning >50 kb each with an identical organization of predicted ORFs. In L4, we found that the expression of six of seven genes examined near the 2.11 fragment Y19D10A.7/F56A4.9 (EGF and L-domain-containing protein), Y19D10A.9/F56A4.2 (lectin C-type domain-containing protein), Y19D10A.10/F56A4.11 (anion, sugar transporter-like), Y19D10A.12/C01B4.9 (monocarboxylate transporter-like), Y19D10A.5/C01B4.8 (substrate transporter-like), and Y19D10A.4/C01B4.7 (aldose 1-epimerase family, sugar transporter-like)was dependent on DAF-12, as daf-12 mutant worms exhibited significantly lower abundance of the corresponding mRNAs compared with N2 animals (Fig. 5B, right panel). The accumulation of the Y19D10A.8/F56A4.10 (anionic transporter-like) transcript differed by less than twofold (we have not determined whether these analyzed transcripts originate from the Y19D10A ORFs or C01B4/F56A4 ORFs).
Stage-specific expression patterns of DAF-12 dependent genes
We measured by qRT PCR the transcripts from the 75 genes proximal to DAF-12 response element-containing fragments as a function of developmental stage and DAF-12 activity. Specifically, we examined transcript abundance in N2 C. elegans relative to that in daf-12, daf-7, and daf-7/daf-12 double mutant worms. Under 25°C growth conditions, all of these animals developed normally for 24 h, reaching L2; by 36 h, however, the N2, daf-12, and daf-7/daf-12 C. elegans were in L3, whereas the daf-7 mutants were in L2d, an extended pre-dauer stage; daf-7 mutants molted into dauer larvae after 50 h.
We found that UDP-glucuronosyltransferase-like F08G5.5 gene expression was down-regulated in daf-12 relative to N2 (Fig. 6A), and that the down-regulation was more pronounced during dauer induction in daf-7; expression was partially restored in daf-7/daf-12 mutants in both the L2 and L3/L2d stages. F08G5.5 was unique in displaying this specific pattern of regulation. A second subset of genes displayed up-regulation during dauer induction (daf-7) in a manner that was dependent on DAF-12: the 3.13 genomic fragment-linked gst-37 gene and 3.10-linked Y38H6C.9 (similar/associated with F-box proteins) gene; Furrowed-like T07H6.5, acyltransferase-like E03H4.7, kinase-like ZC404.9, and F28E10.1 (Fig. 6A). LIT-1b kinase expression was down-regulated in mutant L2 animals relative to wild type, but slightly up-regulated in L3/L2d with strong dependence on DAF-12. Furthermore, at least four genes in a gene cluster linked to the 2.11 response elementY19D10A.7/F56A4.9, Y19D10A.8/F56A4.10, Y19D10A.9/F56A4.2, Y19D10A.12/C01B4.9were strongly expressed in N2 L2 animals relative to daf-12, daf-7, and daf-7/daf-12 mutant worms (Fig. 6B).
|
| Discussion |
|---|
|
|
|---|
By identifying genomic sites bound directly by a regulatory factor of interest, In Vitro Genomic Selection differs from a chromatin immunoprecipitation/microarray-based method (Weinmann et al. 2002
), which yields DNA fragments occupied in vivo either directly or indirectly by the factor. As such, this "ChIP on chip" approach is not as useful as a screening method for discovery of the preferred DNA-binding sites for a given regulatory factor, and to our knowledge, the method had not been applied for that purpose. ChIP on chip also requires technical resources not required by In Vitro Genomic Selection, such as an antibody directed against an epitope accessible when the target protein is assembled into many different complexes, and a whole genome microarray.
We analyzed the genomic fragments identified by our method for response element activity in yeast. Most of the active fragments contained at least one DR5 element, which was sufficient to confer DAF-12-dependent transcriptional activation. Whereas this element architecture implies that DAF-12 may act as a homodimer, a series of additional functional configurations, including heterotypic DR5 half-sites (e.g., high-affinity AGTGCA and low-affinity AGGACA), variations in half-site spacing, and direct rather than inverted repeat orientations invite speculation that DAF-12 might also heterodimerize, perhaps acting in concert with other C. elegans IRs. Such behavior is known for certain vertebrate IRs and could increase substantially the combinatorial diversity of DAF-12 target gene regulation.
Our experiments demonstrated that DAF-12 regulates LIT-1 gene expression in C. elegans, serving either as transcriptional activator or repressor, depending on developmental stage and cell type. Many potential mechanisms might explain the differential activities of DAF-12. For example, the transport or metabolism of a DAF-12 hormone-like ligand might be stage or cell type specific. Alternatively, the expression or activity of a DAF-12 coregulator might be selectively controlled. A third possibility is that other transcriptional regulators may bind at a lit-1 response element and differentially affect the activity of DAF-12. Indeed, examples of all three of these mechanisms have been described for mammalian IRs (Funder et al. 1988
; Yamamoto et al. 1992
; Yoon et al. 2001
).
LIT-1 encodes a MAP kinase that acts in the WNT signaling pathway to regulate cell fate during early embryogenesis and larval development (see Supplemental Materials for references). DAF-12 and LIT-1 are commonly expressed in the same cells, including pharynx, seam, and vulva, and at similar times during development. It will be interesting in future studies to pursue the relationship of DAF-12 and LIT-1. Indeed, in experiments with a temperature-sensitive lit-1 mutant (lit-1 [or131]) or with LIT-1 RNAi, we have observed protruding vulva phenotypes and defects in gonadal migration and development (see Supplementary Fig. S3; Supplementary Table S6), similar to certain daf-12 mutant alleles such as rh61 (Antebi et al. 1998
, 2000
). Although LIT-1 has not been shown previously to affect vulval development, WNT signaling influences that process (Eisenmann et al. 1998
). Moreover, some of the lit-1 mutant phenotypes, such as defects in gonadal migration, vulval morphogenesis, and vacuoles, are reminiscent of aspects of the daf-9 mutant phenotype (Gerisch et al. 2001
; Jia et al. 2002
). DAF-9 is a cytochrome P450-related gene that regulates dauer formation and reproductive development and has been proposed to function in the metabolism of a potential DAF-12 ligand.
We defined candidate DAF-12 response elements as those selected genomic fragments that conferred DAF-12-activated expression in a yeast reporter assay. This criterion likely identifies only a subset of bona fide DAF-12 response elements, as fragments lacking activity in yeast might well be functional in C. elegans. It is intriguing in this regard that one of the selected fragments inactive in yeast maps to the coding region of sdf-9, which affects dauer formation (Ohkura et al. 2003
); it will be interesting to determine if sdf-9 is indeed regulated by DAF-12.
In this study, we surveyed candidate target gene expression by measuring steady-state mRNA levels in whole animals. This is a simple but insensitive assay, as DAF-12-dependent transcriptional regulation of target genes is likely to be cell and stage specific. Moreover, a target gene might be up-regulated in one cell and down-regulated in another, as seen with the LIT-1 GFP reporter, such that even robust regulation might be masked entirely. Nevertheless, our whole-animal, real-time PCR analyses allowed us rapidly to screen many candidates and to identify 40 genes exhibiting DAF-12-dependent regulation (see Tables 1, 3).
|
Identification of DAF-12 target genes, together with a glimpse at their developmental expression patterns, invites speculation on their physiological roles. For example, the UDP-glucuronosyltransferase-like F08G5.5 gene might serve as a checkpoint whose expression and function ensures reproductive development, whereas its absence provokes dauer formation. Similarly, the DAF-12 target genes on the Y19D10A cosmid may provide a signaling function during reproductive development together with their likely roles in metabolic regulation. A group of genes, including F28E10.1, T07H6.5, E03H4.7, Y38H6C.9, gst-37, and ZC404.9, displays DAF-12-dependent up-regulation during dauer induction (Fig. 6A), conceivably triggering early aspects of dauer formation. Indeed, we found that reducing ZC404.9 function by RNAi compromised the ability of N2 worms to form transient dauers at high (27°C) temperature (see Supplementary Table S6).
Numerous DAF-12 targets display putative metabolic roles and dauer program-dependent regulation, perhaps consistent with the shift from aerobic to anaerobic metabolism during dauer formation (Holt and Riddle 2003
). During reproductive development in L4, we observed DAF-12-dependent expression of the Notch pathway factor sel-8/lag-3 (see Supplemental Materials for reference), which functions in vulval and gonadal development, potentially mediating DAF-12 function in these processes. DAF-12-dependent regulation of Skp1 ubiquitin ligase-like skr-8 gene expression may reflect a role for ubiquitin-mediated protein degradation in DAF-12 action. Indeed, skr-8 affects C. elegans lifespan and is regulated by insulin signaling (Murphy et al. 2003
), implying that skr-8 may be downstream of a convergence of the DAF-12 and insulin pathways, where the insulin signaling mediator DAF-16 interacts physically with DAF-12 (Dowell et al. 2003
). Finally, DAF-12-dependent regulation of caveolin expression is intriguing, as caveolin is critical for cellular function and developmental processes in C. elegans and other organisms (Sternberg and Schmid 1999
; van Deurs et al. 2003
), regulating, for example, TGF-
signaling in mammalian cells (Razani et al. 2001
; Di Guglielmo et al. 2003
) and participating in cholesterol transport; cav-1, similarly to skr-8, is a target of insulin signaling in C. elegans (Murphy et al. 2003
). Thus, numerous DAF-12 target genes can be fitted logically into regulatory networks governing aspects of C. elegans reproductive development and dauer formation.
In metazoans such as Drosophila and humans, cis-acting elements that modulate transcription can reside upstream, within, or downstream of the regulated genes, and can be very remote, up to hundreds of kilobases away (Levine and Tjian 2003
). In C. elegans, few such elements have been characterized (Zhang and Emmons 2000
). In this study, we found DAF-12 response elements upstream, in introns and coding regions, and downstream of genes under DAF-12 control (see Supplementary Tables S1S3; Figs. 4, 5, 6). The substantially higher density of genes within the C. elegans genome, relative to Drosophila or mammals, implies that C. elegans harbors mechanisms that tightly couple DAF-12 response elements, for example, to cognate target genes, and exclude nearby nontarget genes from influence by DAF-12. Conceivably, C. elegans uses strategies proposed to function in other metazoans, such as insulators, proximal promoter-enhancer tethering, and gene competition (Levine and Tjian 2003
), or instead has evolved other mechanisms.
We found that C. elegans genomic fragments displaying DAF-12 response element activity in yeast were commonly tightly linked to DAF-12-regulated gene clusters, ranging from at least two to at least seven genes per cluster (the absolute dimensions of the clusters, and their developmental stage specificity, have not been fully defined). Interestingly, clustering of C. elegans muscle-specific genes has been observed, prompting the suggestion that an "active chromatin" domain or a common regulatory element might explain the findings (Roy et al. 2002
). Our finding of clusters in which individual genes appear to be regulated differentially (i.e., DAF-12 activated or repressed within the same cluster) seems to argue against the simplest versions of these mechanisms. This degree of complexity is not unique, however, because bacteriophage studies established more than 30 years ago that a single regulator, acting from a single response element, could regulate one gene positively and another negatively (Reichardt and Kaiser 1971
; Johnson et al. 1979
). The fact that we were unable to identify additional DAF-12-binding sites by computational analysis (FUZZNUC [EMBOSS]: http://bioweb.pasteur.fr/seqanal/interfaces/fuzznuc.html) of sequences within the clusters may indicate that a single response element exerts primary control over multiple genes. It will be interesting to determine if any of the genomic fragments harbor "locus control regions", which have been shown in Drosophila and vertebrates to impose global regulation over clusters of related genes (Grosveld 1999
; Pavlova et al. 1999
; Lee et al. 2003
).
Although response elements and target genes for IRs have been identified in other metazoans, our findings are significant for several reasons: First, we established proof of principle for the notion that biologically functional genetic elements can be isolated from metazoan genomic DNA using a purified DBD. Second, C. elegans encodes 284 IR genes (Sluder and Maina 2001
), far more than any other known organism (c.f., 48 in the human genome [Maglich et al. 2001
], 21 in D. melanogaster [Rubin et al. 2000
]); within C. elegans, DAF-12 is closely related to two other receptors. It will be illuminating to determine roles for the many C. elegans regulators and to assess functional relationships within the DAF-12 subfamily. Finally, it may prove possible in C. elegans, with its sequenced genome, defined cell lineage, powerful genetics, and other experimental advantages, to describe the complete regulatory network for a global effector such as DAF-12; a full appraisal of the "physiological spectrum" of a regulator would enrich our understanding of regulatory evolution and the principles that define regulatory networks.
In Vitro Genomic Selection seems potentially applicable to any DNA-binding protein for which the functional DBD can be purified, yielding naturally occurring DNA-binding sites in their genomic context. Although minimal DNA binding sequences can be derived by random selection (Blackwell and Weintraub 1990
), the biological significance of such sequences is unknown. Similarly, response element and target gene identification through sequential random selection and computational genomic analysis is a useful approach for some regulators (Berman et al. 2002
) but may be less so for IRs, which typically bind to response elements containing relatively few binding sites and high half-site degeneracy. Clearly, the genomic context of factor binding is a critical component of factor activity (Lefstin and Yamamoto 1998
; Yamamoto et al. 1998
; Rogatsky et al. 2002
), and factor-binding sites typically reside within "composite response elements" (Yamamoto et al. 1992
), which include sites for multiple factors, and which nucleate assembly of complexes whose composition and action depend strongly on the precise sites within the element. Fragments isolated by In Vitro Genomic Selection allow protein binding and activity to be assessed either in isolation or in a functional combinatorial context. Finally, the response elements elucidated by In Vitro Genomic Selection provide "tags" to genes that are regulated by the factor of interest, thus leading to identification of direct transcriptional target genes.
| Materials and methods |
|---|
|
|
|---|
Two micrograms of purified (see Methods in Supplemental Materials) recombinant Flag-tagged DAF-12 DBD was immobilized on 30 µL of anti-Flag M2 affinity matrix (Sigma) by incubation for 90 min at 4°C in 150 µL of TBS-2 buffer (50 mM Tris at pH 7.4, 150 mM NaCl, 10 µM Zn(C2H3O2)2, 0.1 mM EDTA). The matrix was then washed with TBS-2 buffer containing 1 M NaCl, followed by two washes with the binding buffer (20 mM HEPES at pH 7.9, 8% glycerol, 10 mM MgCl2, 100 mM KCl, 10 µM Zn(C2H3O2)2). All washes were performed at 4°C for 5 min on a nutator. Immobilized DAF-12 DBD was then incubated with 5 µg of C. elegans genomic DNA digested with 16 units of Sau3AI restriction endonuclease (NEB) overnight, in 150 µL of binding buffer for 25 min at room temperature on a nutator. ProteinDNA bead complexes were washed with 150 µL of the binding buffer for 5 min on nutator at 4°C, followed by a wash with binding buffer containing 250 mM KCl and a third wash containing 500 mM KCl. DAF-12-bound DNA was eluted from the beads with 250 µL of the binding buffer containing 1 M KCl (10 min with 150 µL followed by 15 min with 100 µL; incubations on a nutator).
TE (250 µL) was added to the eluate and the DNA was phenol/chloroform-extracted and precipitated in the presence of 5 µg of glycogen; the pellet was redissolved in 12 µL of H2O. Either 1 or 0.1 pmole of double-stranded NotISau3A oligonucleotide (hybridized 5'-GGGGCGGCCGC-3' to 5'-P-GATCGCGGCCGC CCC-3') was ligated to 10 µL of DNA with 0.5 µL of T4 DNA Ligase (NEB) in 20 µL total volume. Two microliters of the ligation mixture was used as template in a 100-µL first PCR reaction (1.4 µL of Vent DNA Polymerase [NEB], 400 pmole of NotISau3A PCR+ primer [5'-GGGGCGGCCGCGATC-3'], 1x of Thermal Pol Buffer [NEB], 250 µM dNTPs, 0.5 mM extra MgSO4; PCR program: [1] 92°C for 4 min, [2] 92°C for 1 min, [3] 57°C for 45 sec, [4] 72°C for 2 min, [5] GoTo (2) 28 times, [6] 72°C for 5 min, [7] 4°C).
Amplified DNA was purified using Qiaquick columns (Qiagen) and 2 µg of DNA was used for subsequent rounds of binding, performed as described earlier, except a 250-mM KCl wash was used instead of the 500-mM KCl wash. After each binding/elution round, PCR was performed for a total of four rounds. After the second, third, or fourth PCR, purified DNA was digested with NotI restriction endonuclease (NEB) and subcloned into the NotI site of pBluescriptIIKS+ vector (Stratagene) for sequencing and analysis.
qRT PCR analysis of C. elegans genomic fragments
After each round of In Vitro Genomic Selection and amplification, the product mixture was analyzed for relative abundance of individual C. elegans genomic fragments using qRT PCR. Each 50 µL qRT PCR reaction contained 0.5 ng of Sau3A-digested product, 1.25 U of Taq DNA polymerase (Invitrogen), 1.5 mM MgCl2, 125 µM dNTP mix, 0.6 µM of each primer (listed in Supplementary Table S4 for each genomic fragment), and 0.2x SYBR Green I dye (Molecular Probes) in 1x Taq buffer (Invitrogen). qRT PCR was performed in triplicate with an Opticon-2 DNA Engine (MJ Research) and analyzed using the 
Ct method. After 37 cycles, a melting curve of the product was generated between 70°C and 94°C, reading every 0.2°C. The data were normalized to the values of amplification from Sau3A-digested total C. elegans genomic DNA for each genomic selection and control fragment.
Fluorescence anisotropy
Fluorescein-labeled double-stranded oligonucleotides, either 1 or 10 nM of DR5tt and DR5gt (Fig. 3C), in 2 mL of 20 mM Tris (pH 7.4), 100 mM KCl, and 5 mM MgCl2, were used for fluorescence anisotropy measurements with GST-DAF-12 DBD (see Methods in Supplemental Materials). Excitation was performed at 485-nm wavelength and emission was monitored at 515 nm using a Photon Technologies International fluorometer. Anisotropy calculations and Kd determinations by curve fitting were performed as previously described (Gill et al. 1991
; Jagath et al. 1998
) using KaleidaGraph 3.51 (Synergy Software).
Transcriptional reporter and target gene expression assays
See Methods in Supplemental Materials for the yeast transcriptional reporter assay. For the C. elegans transcriptional reporter assay, the wild-type or m11m31 mutant 4.2 fragment was subcloned upstream of a minimal pes-10 promoter driving GFP expression in pPD97.78 vector (plasmids pYSCR0001 and pYSCR0002, respectively). Fifty nanograms per microliter of the constructs were microinjected into wild-type N2 or daf-12 mutant rh61rh412 C. elegans with 50 ng/µL rol-6 marker. daf-12 (rh61rh412) is a loss-of-function class 3 mutant allele characterized by dauer-defective phenotype and impenetrant extragonadal heterochronic phenotype (Antebi et al. 2000
).
For target gene GFP expression assays, the wild-type or m11m31 mutant LIT-1b genomic fragment (nucleotides 45278561 from W06F12 cosmid) was subcloned into pPD95.67 vector to make an in-frame LIT-1GFP fusion (plasmids pYSCR0003 and pYSCR0004, respectively). Ten nanograms per microliter of the constructs were microinjected in wild-type N2 or daf-12 mutant rh61rh412 C. elegans with 50 ng/µL rol-6 marker. pPD97.78 and pPD95.67 vectors were generously provided by Andy Fire.
qRT PCR analysis of putative DAF-12 target gene expression
Wild-type N2 and daf-12(m583) mutant C. elegans strains were grown at 20°C, and daf-7(e1372) and daf-7(e1372)/daf-12(m20) mutant strains were grown at 15°C or 20°C on normal NGM-light or 8x tryptone plates pregrown with OP50 bacteria. Similar to daf-12 (rh61rh412), daf-12(m583) and daf-12(m20) are loss-of-function class 3 mutant alleles characterized by dauer-defective phenotype and impenetrant extragonadal heterochronic phenotype (Antebi et al. 2000
). daf-7(e1372) form constitutive dauers at nonpermissive temperatures (
20°C), and daf-7(e1372)/daf-12(m20) double mutants are dauer defective at all temperatures. daf-7(e1372) and daf-7(e1372)/daf-12(m20) C. elegans stains were generously provided by Pamela Larsen (University of Texas, San Antonio). Worms were growth synchronized by bleaching followed by harvest of embryos. Embryos were plated on NGM-light plates with OP50 bacterial layer and grown at 25°C; larvae were harvested at 24 h for L2, 36 h for L3, and 43 h for L4, and washed three times in 1x M9 buffer before RNA isolation. Total RNA from >50,000 worms was isolated using TRI Reagent (MRC, Inc.) with subsequent purification and DNase I treatment using RNeasy-Mini kit (Qiagen). Random-primed cDNA was prepared from 3 µg of total RNA using the ProtoScript first-strand cDNA synthesis kit (NEB). Each 50 µL qRT PCR reaction was as described earlier, except it contained 120 ng of cDNA and 1 µM of each primer (listed in Supplementary Table S5 for each target gene). qRT PCR was performed in duplicate as described earlier. RNA polymerase II large subunit AMA-1 gene was used as an internal control for data normalization. After 37 or 45 cycles, a melting curve of a product was generated as described earlier. AMA-1 and RGR-1 primers were designed by Wally Wang. The data shown are representative of two to five independent worm growths and RNA isolations.
| Acknowledgments |
|---|
|
|
|---|
| Footnotes |
|---|
Article and publication are at http://www.genesdev.org/cgi/doi/10.1101/gad.1218504.
4 Corresponding author. E-MAIL yamamoto{at}cgl.ucsf.edu; FAX (415) 476-6129. ![]()
| References |
|---|
|
|
|---|
Antebi, A., Yeh, W.H., Tait, D., Hedgecock, E.M., and Riddle, D.L. 2000. daf-12 encodes a nuclear receptor that regulates the dauer diapause and developmental age in C. elegans. Genes & Dev. 14: 1512-1527.
Asahina, M., Ishihara, T., Jindra, M., Kohara, Y., Katsura, I., and Hirose, S. 2000. The conserved nuclear receptor Ftz-F1 is required for embryogenesis, moulting and reproduction in Caenorhabditis elegans. Genes Cells 5: 711-723.[Abstract]
Ashburner, M. 1990. Puffs, genes, and hormones revisited. Cell 61: 1-3.[CrossRef][Medline]
Berman, B.P., Nibu, Y., Pfeiffer, B.D., Tomancak, P., Celniker, S.E., Levine, M., Rubin, G.M., and Eisen, M.B. 2002. Exploiting transcription factor binding site clustering to identify cis-regulatory modules involved in pattern formation in the Drosophila genome. Proc. Natl. Acad. Sci. 99: 757-762.
Blackwell, T.K. and Weintraub, H. 1990. Differences and similarities in DNA-binding preferences of MyoD and E2A protein complexes revealed by binding site selection. Science 250: 1104-1110.
Carmi, I., Kopczynski, J.B., and Meyer, B.J. 1998. The nuclear hormone receptor SEX-1 is an X-chromosome signal that determines nematode sex. Nature 396: 168-173.[CrossRef][Medline]
Chawla, A., Repa, J.J., Evans, R.M., and Mangelsdorf, D.J. 2001. Nuclear receptors and lipid physiology: Opening the X-files. Science 294: 1866-1870.
Di Guglielmo, G.M., Le Roy, C., Goodfellow, A.F., and Wrana, J.L. 2003. Distinct endocytic pathways regulate TGF-
receptor signalling and turnover. Nat. Cell Biol. 5: 410-421.[CrossRef][Medline]
Dowell, P., Otto, T.C., Adi, S., and Lane, M.D. 2003. Convergence of PPARg and Foxo1 signaling pathways. J. Biol. Chem. 278: 45485-45491.
Eisenmann, D.M., Maloof, J.N., Simske, J.S., Kenyon, C., and Kim, S.K. 1998. The
-catenin homolog BAR-1 and LET-60 Ras coordinately regulate the Hox gene lin-39 during Caenorhabditis elegans vulval development. Development 125: 3667-3680.[Abstract]
Ellington, A.D. and Szostak, J.W. 1990. In vitro selection of RNA molecules that bind specific ligands. Nature 346: 818-822.[CrossRef][Medline]
Fisk, G.J. and Thummel, C.S. 1995. Isolation, regulation, and DNA-binding properties of three Drosophila nuclear hormone receptor superfamily members. Proc. Natl. Acad. Sci. 92: 10604-10608.
Freedman, L.P., Arce, V., and Perez Fernandez, R. 1994. DNA sequences that act as high affinity targets for the vitamin D3 receptor in the absence of the retinoid X receptor. Mol. Endocrinol. 8: 265-273.[Abstract]
Funder, J.W., Pearce, P.T., Smith, R., and Smith, A.I. 1988. Mineralocorticoid action: Target tissue specificity is enzyme, not receptor, mediated. Science 242: 583-585.
Gaudet, J. and Mango, S.E. 2002. Regulation of organogenesis by the Caenorhabditis elegans FoxA protein PHA-4. Science 295: 821-825.
Gems, D., Sutton, A.J., Sundermeyer, M.L., Albert, P.S., King, K.V., Edgley, M.L., Larsen, P.L., and Riddle, D.L. 1998. Two pleiotropic classes of daf-2 mutation affect larval arrest, adult behavior, reproduction and longevity in Caenorhabditis elegans. Genetics 150: 129-155.
Gerisch, B., Weitzel, C., Kober-Eisermann, C., Rottiers, V., and Antebi, A. 2001. A hormonal signaling pathway influencing C. elegans metabolism, reproductive development, and life span. Dev. Cell 1: 841-851.[CrossRef][Medline]
Gill, S.C., Weitzel, S.E., and von Hippel, P.H. 1991. Escherichia coli
70 and NusA proteins. I. Binding interactions with core RNA polymerase in solution and within the transcription complex. J. Mol. Biol. 220: 307-324.[CrossRef][Medline]
Gissendanner, C.R. and Sluder, A.E. 2000. nhr-25, the Caenorhabditis elegans ortholog of ftz-f1, is required for epidermal and somatic gonad development. Dev. Biol. 221: 259-272.[CrossRef][Medline]
Godowski, P.J., Rusconi, S., Miesfeld, R., and Yamamoto, K.R. 1987. Glucocorticoid receptor mutants that are constitutive activators of transcriptional enhancement. Nature 325: 365-368.[CrossR