HES-Mediated Repression of Pten in Caenorhabditis elegans.

The hairy/enhancer-of-split (HES) group of transcription factors controls embryonic development, often by acting downstream of the Notch signaling pathway; however, little is known about postembryonic roles of these proteins. In Caenorhabditis elegans, the six proteins that make up the REF-1 family are considered to be HES orthologs that act in both Notch-dependent and Notch-independent pathways to regulate embryonic events. To further our understanding of how the REF-1 family works to coordinate postembryonic cellular events, we performed a functional characterization of the REF-1 family member, HLH-25. We show that, after embryogenesis, hlh-25 expression persists throughout every developmental stage, including dauer, into adulthood. Like animals that carry loss-of-function alleles in genes required for normal cell-cycle progression, the phenotypes of hlh-25 animals include reduced brood size, unfertilized oocytes, and abnormal gonad morphology. Using gene expression microarray, we show that the HLH-25 transcriptional network correlates with the phenotypes of hlh-25 animals and that the C. elegans Pten ortholog, daf-18, is one major hub in the network. Finally, we show that HLH-25 regulates C. elegans lifespan and dauer recovery, which correlates with a role in the transcriptional repression of daf-18 activity. Collectively, these data provide the first genetic evidence that HLH-25 may be a functional ortholog of mammalian HES1, which represses PTEN activity in mice and human cells.

In Caenorhabditis elegans, HLH-25 is one of six members of the REF-1 family of basic helix-loop-helix (bHLH) transcription factors (Alper and Kenyon 2001). This family is characterized by the presence of two bHLH domains and are considered functional orthologs of the hairy/enhancer-of-split (HES) family, partly because of their roles as HES-like transducers of Notch signaling in early embryos (Neves and Priess 2005). REF-1, the first family member to be identified via genetic screens, acts downstream of both Notch and GATA signaling to direct endoderm specification (Hermann ; Neves ) and to direct neuronal lineage decisions (Lanjuin ), including the V-ray lineage in C. elegans males (Ross ). In addition to acting downstream of embryonic Notch signaling events, the family member HLH-29 is targeted by LIN-12/Notch during vulva muscle formation in L4 stage animals (Li ). HLH-29 also acts postembryonically to transcriptionally regulate a diverse set of genes, including those required for chemorepulsive behaviors (McMiller ), IP3-mediated movement of unfertilized oocytes and fertilized eggs through the spermatheca (White ), iron homeostasis, and oxidative stress response (Quach ).

Biological roles have not yet been identified for HLH-25. Notch-dependent, embryonic expression of the gene occurs in descendants of the AB blastomeres, cells that eventually give rise to the pharynx, the nervous system, and the hypodermis (Neves and Priess 2005; Schnabel and Priess 1997). Notch-independent, embryonic expression of is mediated by the GATA factor MED-1 and occurs in descendants of the MS blastomere, which also give rise to the pharynx and muscle and to the somatic gonad (Broitman-Maduro ; Schnabel and Priess 1997). Yeast two-hybrid analysis and protein binding microarrays suggest that the REF-1 family proteins fail to heterodimerize with other bHLH proteins and instead bind E-box-like sequences as homodimers (Grove ; Powell ). Of all the REF-1 family proteins, HLH-25 has the greatest degree of DNA-binding promiscuity and tightly binds to five E-box-like sequences. As a comparison, REF-1 binds tightly and HLH-26 binds weakly to the same E-box-like sequence, whereas HLH-29 only binds weakly to two sequences (Grove ). Gene Ontology (GO) of the genes containing HLH-25 binding sites predicts that HLH-25 plays critical roles in cell division and development, including embryonic development, larval development, growth, and specification of cell fate (Grove ).

HLH-25’s putative role in regulating cell division and embryonic and postembryonic development supports previous studies suggesting that the C. elegans REF-1 family proteins are orthologs of the mammalian HES protein family (Neves ). Like the other REF-1 family proteins, the sequence similarity between HLH-25 and the HES protein family lies within the bHLH domain. HLH-25 is more similar to HES1 than to the other mouse HES proteins but is most similar to the mouse Atonal homolog protein Math1 when amino acid residues flanking the bHLH domain are included in the analysis. Both Atonal1 and Hes1 play well-established roles in directing neuronal cell fate decisions (Akazawa ; Ben-Arie ; Ouji ; Tomita ; Wang ), as do the respective, previously characterized C. elegans orthologs LIN-32 and LIN-22 (Miller and Portman 2011; Wrischnik and Kenyon 1997). Other studies suggest that the human and mouse HES1 proteins are oncogenic (Axelson 2004; Gao ; Sang ), underscoring the roles that these proteins play in regulating cell division and cell proliferation.

In this study, we describe a functional characterization of HLH-25 by using reverse genetics and molecular approaches. We show that is expressed postembryonically in body wall muscles and in neurons and that the phenotypes of animals include reduced brood size, unfertilized oocytes, and abnormal gonad morphology, all of which are consistent with defective cell-cycle or cell division events during embryogenesis. We show that the predicted HLH-25 transcriptional network, based on gene expression microarray, includes the Pten ortholog DAF-18 and other genes that mediate cell division, embryonic development, and larval development. Finally, we show that the ability of HLH-25 to facilitate dauer recovery and lifespan correlates with our mRNA studies suggesting that HLH-25 transcriptionally represses /Pten expression.

Materials and Methods

C. elegans growth and culture conditions

The following strains were used: N2 Bristol wild-type (Brenner 1974); VC1220, II; RB712, (ok480) IV; IC748, quIs18; CMJ3001, cmjEx31(Phlh-25::GFP::unc-54 ’UTR in pPD95.67 + pRF4); CMJ3002, II; cmjEx32[(Phlh-25::hlh-25) + (Pmyo2::mCherry::unc-54 3’UTR)]; and CMJ3003, II; . The extrachromosomal array cmjEx31 includes 1540 nucleotides upstream of the predicted initiation codon, plus the first 30 nucleotides of the coding region for , fused in-frame to sequences coding green fluorescent protein. The extrachromosomal array cmjEx32 includes the genomic DNA sequences spanning from 1540 nucleotides upstream of the initiator codon to 573 nucleotides downstream of the predicted stop codon for . VC1220 animals were outcrossed 10 times to N2 animals, and final homozygosity was confirmed by polymerase chain reaction (PCR). Animals were maintained at 20° on nematode growth media (NGM) agar plates seeded with the Escherichia coli strain OP50 and were synchronized by alkaline hypochlorite treatment as previously described (Lewis and Fleming 1995). Bacteria-mediated RNA interference was performed as previously described (Kamath ; Quach ) via HT115-producing dsRNA for either the control gene or for . The genes and are predicted to produce identical mRNAs and it is impossible to selectively silence .

Microscopy

For imaging, animals were anesthetized either with 0.2% levamisole (Govindan ) or with 10 mM NaN3, mounted on 2% agarose pads containing 10 mM NaN3, and imaged using either a 40× or a 60× water-immersion objective with the Nikon Eclipse 90i microscope equipped with a Nikon Coolsnap CCD camera. Z projections (sums of intensity) were obtained using NIS-Elements AR, version 4.0. For expression profiling, animals were cultured for at least two generations at 16° and then were synchronized and allowed to grow to the appropriate stage (L1−L4). Forty animals per larval stage were examined for expression. For dauer stage expression, 10 starvation-induced dauer were imaged.

Egg-laying assays

L4-stage animals were selected to individual NGM plates and allowed to molt to adulthood at 20°. At the start of the egg-laying period, each animal was moved daily to a new plate, already seeded with OP50, until 3 d after eggs or unfertilized oocytes were no longer detected. Eggs and oocytes were counted on the day that the hermaphrodite was removed. The plate was then incubated at 20° for 24 hr, and the numbers of live progeny and unhatched eggs were counted. Eggs that did not hatch within the 24-hr period were scored as dead, whereas the number of live progeny was checked again in the next 24 hr and then removed. For rescue assays, the ratio of unfertilized oocytes to fertilized eggs was calculated daily.

Lifespan assays

Lifespan assays were performed as previously reported (Quach ). Lifespan survival comparisons were done individually for each of three biological replicates using the log rank test in GraphPad Prism 6 (GraphPad Software, Inc) and were then repeated with the data from all biological replicates combined into one experiment. Combined data are presented as a graph in Figure 5. Individual data are shown in Table S2.

Figure 5

HLH-25 regulates lifespan and dauer recovery. (A) hlh-25(ok1710) animals (n = 282) have longer than lifespans than wild-type animals (n = 264, P-value < 0.0001). (B) Representative experiment showing the percentage of animals with pharyngeal pumping during dauer recovery at 20°. P-value < 0.0001 and P-value = 0.938 when hlh-25(ok1710) animals and daf-18(ok480) animals, respectively, were compared with wild-type. (C) Representative experiment showing the percentage of animals with pharyngeal pumping during dauer recovery at 20°. P-value < 0.0001 and P-value = 0.2771 when hlh-25(ok1710) and hlh-25(ok1710;cmjEx32) animals, respectively, were compared with wild-type.

Dauer recovery assays

To induce dauer formation, embryos were collected by hypochlorite treatment, washed extensively in water, and then approximately 200 embryos were placed onto an NGM plate containing a sparse lawn of OP50 (∼10 µL of saturated OP50, 24 hr before use). Embryos were allowed to hatch, and the dauer larva to develop, undisturbed at 27° for at least 96 hr. Cultures were visually examined to ensure that at least 90% of the larvae were arrested as dauer. Failure to include OP50 on the plate resulted in animals that arrested as starved L1 stage larvae. For dauer recovery, 25 dauers per strain were individually moved to a plate with a thick lawn of OP50 and incubated at 20°. Dauers were examined every 3 hr during a 24-hr period for fat accumulation, pharyngeal pumping, crescent formation, and then at 8-hr intervals for another 24 hr for oocyte/embryo formation, and onset of egg-laying. The rates of onset of pharyngeal pumping were plotted as survival from dauer, and were compared by log-rank test using GraphPad Prism 6. The timing of dauer recovery varied with the day of the assay, with the onset of pumping occurring in 50% of wild-type animals as early as 3 hr or as late as 8.5 hr. In all assays the trend was the same, and the data shown are representative of one experiment.

Total RNA isolation, cDNA synthesis, reverse transcription, and quantitative PCR (qPCR)

Total RNA extraction and cDNA synthesis were carried out as previously described (Quach ) except total RNA was extracted from synchronized L4 stage animals. After cDNA synthesis, qPCR assays were performed with Taqman Gene Expression Assays (Invitrogen) for detection of amplicon, using relative quantitation with normalization against the endogenous control gene (Hoogewijs ). Values shown represent the averages of at least three different experiments. The ID numbers for the Taqman Assays used in this study are available upon request.

Gene expression microarray

Gene expression microarray (GEO accession #GSE65417), including probe preparation, hybridization, fluidics run and chip scan, was performed by Georgia State University DNA/Protein Core Facility. Global gene expression in synchronized populations of VC1220 animals was compared to expression in N2 (wild-type) animals using GeneChip C. elegans Genome Array (Affymetrix). Data collection was carried out using GCOS 1.4 software (Affymetrix). Data analysis was performed using GeneSpring GX 11 Software (Agilent, Palo Alto, CA), and probe intensity values were normalized using robust multichip average algorithm. The quality controls on samples and on probe sets were performed stepwise to detect the outlying samples and the poor probe sets. The principal components analysis score plot and hybridization controls plot were applied for sample detection. Those probe sets passing minimal detection cutoffs and quality control measurements subsequently were used for statistical analysis. The Student’s t-test was performed to find the candidates for differential expression, and genes with significant signal level between different conditions (P < 0.05) were collected. In addition, fold change analysis were performed on the genes with significant expression, and all genes showing greater than two-fold change were considered putative targets.

Functional analysis by GO

GO analysis was performed via the Database for Annotation, Visualization and Integrated Discovery (i.e., DAVID), version 6.7, to cluster related target genes (Huang da ,b). Additional GO terms and functional information for putative targets were assigned based on Wormbase annotations (Harris ). Interactions and connectivity between target gene products were identified using the Search Tool for the Retrieval of Interacting Genes/Proteins (Szklarczyk ) and redrawn manually.

Mobility assays

For locomotion assays, individual L4 stage animals were placed by platinum wire onto an unseeded NGM plate, allowed to recover for 5 min, and then examined after the start of forward locomotion for 30 sec. For thrashing assays, individual animals were placed into 50 µL of M9 medium. Each lateral movement of the head was counted as one thrash, and movement was monitored for one 30-sec interval per animal. Each biological replica included at least 25 animals per strain, and each assay was repeated for three biological replicas.

Data availability

The Gene expression microarray (GEO accession #GSE65417) will release for public access on January 1, 2016.

Results

Phenotypes of hlh-25(ok1710) animals and Phlh-25::GFP activity suggest roles for HLH-25 in reproduction

Previous studies indicate that is expressed in eight-cell stage embryos in response to Notch signalling, and in later stage embryos in response to GATA signalling (Broitman-Maduro ; Neves and Priess 2005). Using an hlh-25::GFP transcriptional reporter (supporting information, Figure S1), we detected postembryonic expression of at all larval and adult stages, including dauer stage, in unidentified neurons of the head and tail, the anal depressor muscle, anterior and posterior intestinal cells, and the head and body wall muscles of transgenic animals (Figure 1). We detected expression in the PDE sensory neuron pairs and rare, transient expression in migrating distal tip cells of L3/L4 stage animals. We also detected stage-specific expression in multiple neurons of the ventral nerve cord in starvation-induced dauer stage animals. We were able to use cell body location and commissure morphology to identify three of the motor neurons as VD1, DA3, and DA7. Expression in the body wall and head muscles and in the PDE neuron pair was not evident in dauer larvae.

Figure 1

Postembryonic expression of Phlh-25::GFP. (A−D) Nomarski images (left) accompanied by green fluorescent protein fluorescent images (right). (A) Expression of cmjEx31 in multiple head neurons (arrowheads). This pattern was detected in animals at all stages, with a few additional neurons at dauer stage. (B) L3 expression in distal tip cell (arrow, red dashed circle), unidentified head neuron (arrowhead) and head muscles (asterisk). (C) L3 stage expression in PDE (arrow). Arrowheads indicate PDE processes. (D) Dauer stage expression in ventral nerve cord neurons. In all images anterior is to the right, dorsal is up. Scale bar represents 50 µm.

The allele is a null allele that was generated by ethyl methanesulfonate mutagenesis. It is a 1550-bp deletion that spans from 292 bp upstream of the initiator codon to 348 bp downstream of the terminator codon (Figure S1). We outcrossed these animals 10 times, confirmed homozygosity for the allele by PCR, and then examined them for mutant phenotypes. animals produced fewer live progeny at 20° than wild-type animals (Figure 2A, P-value = 0.0027) whereas laying more unfertilized oocytes (P-value = 0.00014) throughout the egg-laying period.

Figure 2

Brood size and unfertilized oocytes hlh-25(ok1710) animals. (A) hlh-25(ok1710) animals have fewer live progeny and lay more unfertilized eggs than wild-type animals. Expression of the extrachromosomal array cmjEx32, which carries hlh-25 genomic sequences, rescues this phenotype in hlh-25(ok1710) animals. Graph shows mean of total numbers of live progeny (dark bars) and unfertilized eggs (light bars) produced during the entire lifespans of wild-type (n = 7), hlh-25(ok1710) (n = 12), and hlh-25(ok1710; cmjEx32) animals (n = 9). *P-value < 0.05, **P-value < 0.005. (B) Average numbers of live progeny (solid lines) and unfertilized eggs (dashed lines) produced per day of egg laying by the same hlh-25(ok1710) animals (triangles) and wild-type animals (circles) as those shown in (A).

Adult hermaphrodites produce more oocytes than sperm and will continue to lay unfertilized eggs for a short period after depleting their sperm supply. Under conditions where food is plentiful, and when cultured at 20°, wild-type hermaphrodites lay eggs over the first 5 d of adulthood (Ward and Carrel 1979; White ). As previously described, we found that wild-type animals lay an increasing number of unfertilized eggs as the sperm supply diminishes. In our assays, wild-type animals continued to lay a small number of unfertilized eggs for one additional day after the active egg-laying period ended (egg-laying day 6), to produce an average total of 49 unfertilized oocytes (Figure 2B, Table S1). animals continued to produce unfertilized oocytes for 4 days after the active egg-laying period ended (egg-laying day 9), and produced an average total of 232 unfertilized oocytes. This phenotype has been recently named uno-o (Riesen ), and is also seen in animals that carry null alleles of the REF-1 family protein HLH-29 (White ). We were able to partially rescue the uno-o phenotype in animals by expressing from an extrachromosomal array, , so that transgenic animals produced an average of 38% fewer unfertilized oocytes than animals (Figure 2A).

We also examined the gonad arms of adult hermaphrodites. In wild-type animals, oocytes production occurs as progressive movement through two U-shaped gonad arms, starting with the production of germline nuclei at the distal end of each gonad arm, and ending with the ovulation of the oocyte into the sperm-filled spermatheca at the ends proximal to the uterus (Hubbard and Greenstein 2000). In the distal, syncytial region of each arm, the germline nuclei go through mitotic and meiotic divisions while being surrounded by a commonly shared cytoplasmic core. Nuclei start to compartmentalize into oocytes as they move progressively closer to a bend in the gonad arm, and are fully compartmentalized, or budded, as they leave the bend and proceed towards the spermatheca (Green ). We found that organization of the gonad arm and germline progression were abnormal in 63% of animals. As shown in Figure 3, and in Figure S2 and Figure S3, gonad architecture was variably abnormal and included constrictions in the distal and proximal gonad arms, increased apoptosis, increased oocyte number, delayed oocyte compartmentalization, abnormal chromatin condensation, and arrested embryos in the uterus.

Figure 3

Gonad and oocyte morphologies of hlh-25(ok1710) animals. Nomarski (A, C) and merged fluorescent (B, D) images of representative gonad architecture and oocyte morphology in wild-type animals (A, B) as analyzed in a strain coexpressing a fluorescent chromosome marker (mCherry-Histone H2b) and a marker targeting green fluorescent protein to the plasma membrane (PI4,5P2) (Green ), (B). Representative gonad architecture and oocyte morphology in hlh-25(ok1710) animals (C, D). Phenotypes include increased oocyte number, distal gonad arm constriction (arrow), increased apoptosis (^), and arrested embryos in uterus (#). Insert in (C) shows magnified view of distal gonad arm constriction. Inserts in (B) and (D) show embryonic nuclei. In all images, ventral is up, anterior is to the right, scale bar = 50 µm.

Other phenotypes of animals included decreased movement as monitored by thrashing and locomotion assays (Figure S4), leading us to classify mutants as locomotion variants. This phenotype correlates with the expression of in the body wall muscles and in head neurons, but we could not accurately assay for rescue of this phenotype by transgenic expression of because transgenic animals also carried the marker, which causes a dominant locomotor defect. We did not pursue the locomotion phenotype further in this study.

DAF-18 is one major hub in the hlh-25 transcriptional network

We used gene expression microarray to compare the expression profiles of late L4/young-adult stage animals to profiles of age-matched, wild-type animals. We found that the expression of 630 genes, 506 of which were up-regulated and 124 of which were down-regulated, was altered by at least 2.0-fold (Table S3), and confirmed by reverse transcription (RT)-qPCR that the expression of 21 genes was significantly affected in L4-stage animals (Table 1). As shown in Table 1, the direction of the change in expression was the same for 18 of the 21 genes. It is possible that the opposite, but significant change in the three remaining genes was due to differences in the developmental age of animals at the times of RNA extraction. Interestingly, 152 of these genes, indicated by asterisks in Table 1 and Table S3, contain predicted HLH-25 binding sites within their putative promoters (Grove ). Hypergeometric analysis, indicated that this overlap of 152 genes was significantly lower than expected when compared to an overlap between the list of genes identified by protein binding microarray to a randomly chosen gene set of the same size (P-value <2.5 × 10-05).

Table 1

Validation of randomly selected HLH-25 targets by RT-qPCR

WormBase ID	Sequence ID	Gene Name	log2 FC (Microarray)	log2 FC (RT-qPCR)
WBGene00000230	F27C1.7	atp-3	1.14405	1.203510
WBGene00015102	B0280.5	cpg-2	1.29866	2.277368
WBGene00000913	T07A9.6	daf-18a	1.57550	1.181618
WBGene00001263	K04H4.1	emb-9	1.77821	0.688850
WBGene00010305	F59A2.5	F59A2.5a	−1.17632	−1.630390
WBGene00001394	W02A2.1	fat-2a	1.03562	1.718458
WBGene00014095	ZK829.4	gdh-1a	2.48027	1.338850
WBGene00001758	Y45G12C.2	gst-10a	1.36737	0.774000
WBGene00003022	ZK418.4	lin-37	1.58496	1.489140
WBGene00003037	JC8.8	lin-54	0.59567	0.425533
WBGene00003230	W02A2.7	mex-5	1.04264	0.574521
WBGene00003231	AH6.5	mex-6	1.11770	0.656000
WBGene00003242	C37C3.6	mig-6	1.71370	0.877350
WBGene00003473	K11G9.6	mtl-1a	−4.08804	−1.578000
WBGene00004078	F52E1.1	pos-1	1.91456	1.017858
WBGene00004302	K01G5.4	ran-1a	1.49570	0.989226
WBGene00004736	K11D9.2	sca-1	1.78660	1.270413
WBGene00004984	ZC404.8	spn-4	2.06005	0.947934
WBGene00009221	F28F8.2	acs-2	−1.71450	1.107159
WBGene00003421	F09E8.3	msh-5	1.70929	−1.120290
WBGene00011733	T12D8.5	T12D8.5	−3.24793	0.441060

RT-qPCR, reverse transcription quantitative polymerase chain reaction.

Target previously predicted by protein binding microarray (Grove ).

We clustered related target genes based on enriched GO terms. A total of 328 genes were organized into 10 annotation clusters that all fit within the broad category of “Development” (top eight clusters are shown in Figure 4A). Using the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING 9.0) to identify predicted and known interactions between proteins encoded by the target genes (Szklarczyk ), we found that the regulatory network has a high degree of overlap and interconnectivity. For example, 38 genes clustered under the GO term “Sexual Reproduction.” Of those, 22 genes (58%) formed a single, highly interconnected, interaction network (Figure 4B) and also were associated with at least two other clusters. This trend of high connectivity among genes within a cluster and significant overlap among genes between clusters was evident for the entire network. We also identified six genes that we considered to be major hubs (Table 2). These genes were associated with at least four of the GO clusters and formed predicted interactions with at least 35 other targets in the network. Importantly, the hub genes underscore the developmental influence of the regulatory network, encoding helicases, cyclins, and RNA-binding proteins with established roles in embryogenesis (Table 2). The C. elegans homolog of the mammalian Pten tumor suppressor, DAF-18 (Masse ), is one of the major hubs, forming connections with 62 proteins distributed across six annotation clusters (Figure 4C, bold genes in Table S3).

Figure 4

Functional annotation of HLH-25 targets. (A) Distribution of HLH-25 targets into top eight GO clusters. (B) Predicted and known interactions between genes in the ‘sexual reproduction’ cluster. The hub genes cgh-1, cyb-2.1, cyb-2.2, cyb-3, and puf-3 are outlined in red. The hub gene daf-18 is not a part of the ‘sexual reproduction’ cluster. (C) Representative genes from the 62 HLH-25 targets that form predicted or known interactions with daf-18. All genes, except daf-18, clustered under “embryonic development ending in birth (EDB)” Most genes also clustered under at least one other top GO category, including “nucleoside/nucleotide transport” (green spheres).

Table 2

Hub genes in the HLH-25 network

Wormbase ID	Sequence ID	Gene Name	Functional Summarya	GO Cluster Association
				ED	S	P	M	G	C
WBGene00000479	C07H6.5.2	cgh-1	ATP-dependent helicase; prevents apoptosis of developing embryos	X	X	X	X	X
WBGene00000866	Y43E12A.1	cyb-2.1	Cyclin B family; required for oocyte maturation	X	X	X		X	X
WBGene00000868	T06E6.2a.2	cyb-3	G2/mitotic-specific, Cyclin B3	X	X	X			X
WBGene00000913	T07A9.6	daf-18	Phosphatase; negative regulator of insulin-like signaling	X			X	X
WBGene00003155	C25D7.6.2	mcm-3	DNA replication licensing factor	X			X		X
WBGene00004239	Y45F10A.2.1	puf-3	RNA binding protein; required for cell-cycle timing, spindle positioning, pronuclei formation	X	X				X

GO, Gene Ontology; ED, embryonic development ending in birth; S, sexual reproduction; P, nucleoside, proton transport; M, meiosis; G, regulation of growth; C, cell cycle.

Functional summaries are based on Wormbase annotations (www.wormbase.org).

HLH-25 functions upstream of DAF-18 to regulate dauer recovery

It is not clear whether the REF-1 proteins and the HES proteins arose from the same parent gene; however, conservation within the helix-loop-helix domains and the ability to transduce Notch signals suggest that REF-1 proteins and HES proteins are members of the same functional family. Importantly, the identification of as a major hub in the HLH-25 transcriptional network, together with the identification of predicted HLH-25 binding sites in the promoter (Grove ), underscores the functional homology between the HLH-25 and vertebrate HES1 protein. In the mouse thymus, HES1 acts downstream of Notch signalling to repress Pten activity during T-cell differentiation (Wong ), and the resistance of T-cell acute lymphoblastic leukemia cells to Notch inhibition depends on functional HES1 binding sites in the Pten promoter (Palomero ). Therefore, we further characterized the genetic interaction between the genes encoding the HES homolog HLH-25 and the PTEN homolog DAF-18.

Using RT-qPCR, we found that the levels of mRNA increase greater than two-fold in embryos (data not shown) and in L4 stage animals (Table 1) compared with age-matched controls. The levels of mRNA also increased in wild-type animals that were subjected to RNA interference (Figure S5), suggesting that the increased activity is due to loss of HLH-25 in animals. We failed to detect differences in transcriptional activity in wild-type and animals expressing Pdaf-18::GFP via a previously described, integrated transgene. These animals were generated in a loss-of-function background, which increases the transcriptional activity of the promoter (Solari ). It is possible that loss of masks the effects of HLH-25. Unfortunately, we have not yet been successful in generating stable transgenic lines expressing Pdaf-18::GFP in animals that produce functional DAF-2.

In C. elegans, longevity and entry into dauer stage are controlled through the insulin signalling pathway, and depend on the insulin receptor DAF-2. Decreased insulin signalling, via loss-of-function mutations in , increases C. elegans lifespan and causes a dauer-constitutive phenotype. These phenotypes are suppressed in animals: animals have reduced lifespan (Solari ) and animals that carry strong loss-of-function alleles are unable to undergo dauer development (Gil ; Mihaylova ; Rouault ). The reduced lifespan of animals can be rescued, and slightly extended beyond the wild-type lifespan, by the overexpression of from the extra-chromosomal array (Brisbin ; Solari ). We reasoned that since animals have increased activity, they should have an extended lifespan. As expected, the lifespan of animals (median LS = 19 d) is significantly longer than the lifespan of wild-animals (median LS = 16 d, P-value <0.0001) grown at 20° (Figure 5A, Table S2). This result is comparable with the extended lifespan previously reported for C. elegans females (Brisbin ).

Normally when wild-type, starvation-induced dauer larvae are moved to favorable growth conditions, they take approximately 60 minutes to commit to exiting the dauer stage at 25°, and they molt into the L4 stage at approximately 10 hr (Golden and Riddle 1982; Golden and Riddle 1984a,b,c). As animals begin to exit dauer stage and enter into the L4 molt, fat accumulation and pharyngeal pumping begin after approximately 60 min and 3 hr, respectively, of making the dauer recovery/exit decision (Proudfoot ). The onset of vulva formation is detectable within 12 hr of the exit decision, and recovering dauers enter late L4/early adult stage within 24 hr of the decision to exit dauer stage.

Because animals do not form dauers under starvation conditions, we used heat stress to induce dauer formation via an IIS-independent pathway in wild-type, , , and animals, and compared their dauer recovery phenotypes. We found that at least 80% of wild-type animals (n = 1199) reached early L4 stage within 12 hr after the switch to 20°. As expected, dauer recovery was similar to wild-type in animals (n = 485) but was longer in animals (n = 955). At 12 hr, approximately 65% of the animals, but only 10% of the animals, reached early L4. Interestingly, the dauer recovery timing of the animals (n = 738) was very similar to the animals: only 21% of animals reached L4 stage within 12 hr, and only 45% reached L4 stage within 24 hr. This delayed dauer recovery phenotype in animals was partially rescued by transgenic expression of (): 70% of transgenic animals (n = 300) reached L4 within 24 hr. For a more quantitative approach, we monitored the onset of pharyngeal pumping during dauer recovery. Dauer recovery in animals was not significantly different from recovery in wild-type animals; however, it took 1.5 times longer for the onset of pharyngeal pumping to occur in 50% of () animals (P-value <0.0001) than in wild-type animals (Figure 5B). This delay was rescued by by transgenic expression of (Figure 5C). These results, together with our gene expression data, suggest that HLH-25 acts upstream of to negatively regulate lifespan and dauer recovery.

Discussion

This study is a functional characterization of HLH-25, a member of the C. elegans REF-1 family of transcriptional regulators. Previous studies have shown that the gene is transcriptionally activated through Notch-signaling in cells that eventually give rise to the pharynx and nervous system, and through the Notch-independent GATA family of transcription factors in tissues that give rise to the pharynx, muscle, and somatic gonad. Here we show that expression of persists throughout the C. elegans life cycle. Although we identified a number of HLH-25 targets that are germline expressed and whose mRNAs are stored in oocytes before fertilization, we were unable to recover transgenic lines that report germline expression of . Future studies may further define the role of HLH-25 in embryogenesis, which includes determining whether the maternal contribution of HLH-25 is required for proper embryogenesis.

Our results show that the phenotypes of animals include reduced brood size, unfertilized oocytes in the uterus, abnormal gonad morphology/architecture, extended lifespan, and delayed dauer recovery. In previous studies, functional annotation of candidate target genes suggested overlapping roles in cell signaling events for all of the REF-1 family proteins, and HES-like roles for HLH-25 and HLH-29 in regulating embryonic and larval development, growth, cell fate specification, and reproductive behaviors (Grove ). Interestingly, this functional overlap is suggested, but not clearly evident, in the phenotypes of REF-1 family mutants that have been characterized to date. Both and animals, for example, show varying degrees of embryonic lethality that would ultimately lead to reduced brood size (Alper and Kenyon 2001; McMiller ; White ); however, the reduced brood size in animals in not due to lethality during embryogenesis, but to abnormalities that occur before fertilization. Likewise, mutations in , , and cause morphologic defects that affect sexual reproduction. Loss of affects cell fusion events during vulva morphogenesis (Alper and Kenyon 2001) and ray lineage decisions during morphogenesis of the male tail (Ross ). Loss of either or results in unfertilized oocytes, abnormal gonad morphologies, and the exploded through vulva phenotype; however, the underlying mechanisms of these phenotypes appears to differ in and animals. The unfertilized oocyte phenotype in animals, for example, is caused in part by altered expression of genes required for inositol 1, 4, 5-triphosphate (IP3) signaling, which affects the ability of oocytes and eggs to move through the spermatheca (White ). On the basis of our observations of independent ovulation events, however, spermatheca function in animals appears normal, and genes required for IP3 signaling are not affected by loss of HLH-25. These data suggest that the REF-1 family proteins do not act within a single transcriptional cascade to regulate specific events, and that each family member plays unique, stage-specific roles in C. elegans development.

Our gene expression microarray analysis suggests that there are six major hubs in the HLH-25 transcriptional network, including the C. elegans PTEN homolog, . We demonstrated by RT-qPCR that expression is normally repressed by HLH-25. This result is the first direct correlation between a HES protein target gene in mammals and a REF-1 family target gene in C. elegans. Hes1 directs cell division, in part, by repressing Pten expression in mouse (Wong ), in human large cell neuroendocrine carcinoma (Nasgashio ), and in regenerating Clara cells of wild-type lungs (Xing ). In C. elegans, DAF-18 is required during embryogenesis to maintain G2 arrest in germline precursor cells (Fukuyama ) via inhibition of Tor complex 1 signaling (Fukuyama ). We did not examine embryonic or L1-stage phenotypes of animals; however, we show that, like animals that overexpress (Brisbin ; Liu and Chin-Sang 2015), animals have an extended lifespan and take longer to exit dauer stage. Further examination of both the embryonic and postembryonic regulation of by HLH-25 will further our understanding of the mechanisms and biological significance of HES-dependent activation and repression of PTEN.

C. elegans enter the dauer stage in response to adverse growth conditions, maintain this state while conditions remain undesirable, and then exit, or recover from dauer stage when environmental conditions are more favorable (Hu 2007; Ouellet ). The decision and ability to enter into dauer stage requires the integration of dietary and environmental signals through multiple molecular pathways, including insulin/insulin-like growth factor, transforming growth factor-β, and guanylyl cyclase, and the molecular details of dauer entry in response to starvation, heat stress, and overcrowding have been well established (reviewed in Fielenbach and Antebi 2008). Less is known about the molecular requirements for dauer maintenance and dauer recovery. Studies suggest that GLP-1−dependent Notch activation in chemosensory neurons is required for dauer maintenance (Ouellet ) and that metabotropic acetylcholine signaling is required for the induction of an insulin-like signal during dauer recovery (Tissenbaum ). Previous studies have also shown that dauer recovery requires LIN-12−dependent Notch activation, though the downstream target of this signaling event has not been identified (Ouellet ). Future studies will determine whether HLH-25 acts downstream of Notch to modulate dauer recovery by repressing activity.

Finally, GO and protein interaction analysis of the HLH-25 network suggests that HLH-25 affects genes required for three highly connected processes: cell division, cell-cycle regulation, and reproduction, including gamete formation. The connectivity between the genes in these processes suggests that the HLH-25 regulatory network is amendable to change without the corresponding disruption in overall functionality. This feature is illustrated by the hub gene , and the hub genes , , , , and , which encode proteins required for early embryonic cytokinesis, meiosis II, G2/mitosis transition, licensing of DNA replication, and spindle positioning, respectively (Audhya ; Boag ; Hubstenberger ; Kaitna ; Sonneville ; van der Voet ). Although mutations in any of the hub genes, or in any of the nodes connected to them, result in severe phenotypes, including sterility, mutant animals are viable and have functional cell cycles. Future work will include further defining and validating the HLH-25 transcriptional network and correlating this network with the postembryonic functions of HLH-25.

Taken together, our results demonstrate that the C. elegans HES homolog, HLH-25, functions postembryonically to transcriptionally regulate genes required for cell division, cell-cycle regulation, and sexual reproduction, including the gene encoding the C. elegans Pten homolog, DAF-18.