A conserved chronobiological complex times C. elegans development.

The mammalian PAS-domain protein PERIOD (PER) and its C. elegans orthologue LIN-42 have been proposed to constitute an evolutionary link between two distinct, circadian and developmental, timing systems. However, while the function of PER in animal circadian rhythms is well understood molecularly and mechanistically, this is not true for the function of LIN-42 in timing rhythmic development. Here, using targeted deletions, we find that the LIN-42 PAS domains are dispensable for the protein's function in timing molts. Instead, we observe arrhythmic molts upon deletion of a distinct sequence element, conserved with PER. We show that this element mediates stable binding to KIN-20, the C. elegans CK1δ/ε orthologue. We demonstrate that CK1δ phosphorylates LIN-42 and define two conserved helical motifs, CK1δ-binding domain A (CK1BD-A) and CK1BD-B, that have distinct roles in controlling CK1δ-binding and kinase activity in vitro . KIN-20 and the LIN-42 CK1BD are required for proper molting timing in vivo . These interactions mirror the central role of a stable circadian PER-CK1 complex in setting a robust ∼24-hour period. Hence, our results establish LIN-42/PER - KIN-20/CK1δ/ε as a functionally conserved signaling module of two distinct chronobiological systems.

Dynamics of miRNA accumulation during C. elegans larval development.

Temporally and spatially controlled accumulation underlies the functions of microRNAs (miRNAs) in various developmental processes. In Caenorhabditis elegans, this is exemplified by the temporal patterning miRNAs lin-4 and let-7, but for most miRNAs, developmental expression patterns remain poorly resolved. Indeed, experimentally observed long half-lives may constrain possible dynamics. Here, we profile miRNA expression throughout C. elegans postembryonic development at high temporal resolution, which identifies dynamically expressed miRNAs. We use mathematical models to explore the underlying mechanisms. For let-7, we can explain, and experimentally confirm, a striking stepwise accumulation pattern through a combination of rhythmic transcription and stage-specific regulation of precursor processing by the RNA-binding protein LIN-28. By contrast, the dynamics of several other miRNAs cannot be explained by regulation of production rates alone. Specifically, we show that a combination of oscillatory transcription and rhythmic decay drive rhythmic accumulation of miR-235, orthologous to miR-92 in other animals. We demonstrate that decay of miR-235 and additional miRNAs depends on EBAX-1, previously implicated in target-directed miRNA degradation (TDMD). Taken together, our results provide insight into dynamic miRNA decay and establish a resource to studying both the developmental functions of, and the regulatory mechanisms acting on, miRNAs.

Convergence of multiple RNA-silencing pathways on GW182/TNRC6.

The RNA-binding protein TRIM71/LIN-41 is a phylogenetically conserved developmental regulator that functions in mammalian stem cell reprogramming, brain development, and cancer. TRIM71 recognizes target mRNAs through hairpin motifs and silences them through molecular mechanisms that await identification. Here, we uncover that TRIM71 represses its targets through RNA-supported interaction with TNRC6/GW182, a core component of the miRNA-induced silencing complex (miRISC). We demonstrate that AGO2, TRIM71, and UPF1 each recruit TNRC6 to specific sets of transcripts to silence them. As cellular TNRC6 levels are limiting, competition occurs among the silencing pathways, such that the loss of AGO proteins or of AGO binding to TNRC6 enhances the activities of the other pathways. We conclude that a miRNA-like silencing activity is shared among different mRNA silencing pathways and that the use of TNRC6 as a central hub provides a means to integrate their activities.

C. elegans molting requires rhythmic accumulation of the Grainyhead/LSF transcription factor GRH-1.

C. elegans develops through four larval stages that are rhythmically terminated by molts, that is, the synthesis and shedding of a cuticular exoskeleton. Each larval cycle involves rhythmic accumulation of thousands of transcripts, which we show here relies on rhythmic transcription. To uncover the responsible gene regulatory networks (GRNs), we screened for transcription factors that promote progression through the larval stages and identified GRH-1, BLMP-1, NHR-23, NHR-25, MYRF-1, and BED-3. We further characterize GRH-1, a Grainyhead/LSF transcription factor, whose orthologues in other animals are key epithelial cell-fate regulators. We find that GRH-1 depletion extends molt durations, impairs cuticle integrity and shedding, and causes larval death. GRH-1 is required for, and accumulates prior to, each molt, and preferentially binds to the promoters of genes expressed during this time window. Binding to the promoters of additional genes identified in our screen furthermore suggests that we have identified components of a core molting-clock GRN. Since the mammalian orthologues of GRH-1, BLMP-1 and NHR-23, have been implicated in rhythmic homeostatic skin regeneration in mouse, the mechanisms underlying rhythmic C. elegans molting may apply beyond nematodes.

A specific type of Argonaute phosphorylation regulates binding to microRNAs during C. elegans development.

Argonaute proteins are at the core of the microRNA-mediated gene silencing pathway essential for animals. In C. elegans, the microRNA-specific Argonautes ALG-1 and ALG-2 regulate multiple processes required for proper animal developmental timing and viability. Here we identified a phosphorylation site on ALG-1 that modulates microRNA association. Mutating ALG-1 serine 642 into a phospho-mimicking residue impairs microRNA binding and causes embryonic lethality and post-embryonic phenotypes that are consistent with alteration of microRNA functions. Monitoring microRNA levels in alg-1 phosphorylation mutant animals shows that microRNA passenger strands increase in abundance but are not preferentially loaded into ALG-1, indicating that the miRNA binding defects could lead to microRNA duplex accumulation. Our genetic and biochemical experiments support protein kinase A (PKA) KIN-1 as the putative kinase that phosphorylates ALG-1 serine 642. Our data indicate that PKA triggers ALG-1 phosphorylation to regulate its microRNA association during C. elegans development.

Coupling of growth rate and developmental tempo reduces body size heterogeneity in C. elegans.

Animals increase by orders of magnitude in volume during development. Therefore, small variations in growth rates among individuals could amplify to a large heterogeneity in size. By live imaging of C. elegans, we show that amplification of size heterogeneity is prevented by an inverse coupling of the volume growth rate to the duration of larval stages and does not involve strict size thresholds for larval moulting. We perturb this coupling by changing the developmental tempo through manipulation of a transcriptional oscillator that controls the duration of larval development. As predicted by a mathematical model, this perturbation alters the body volume. Model analysis shows that an inverse relation between the period length and the growth rate is an intrinsic property of genetic oscillators and can occur independently of additional complex regulation. This property of genetic oscillators suggests a parsimonious mechanism that counteracts the amplification of size differences among individuals during development.

H3K9me selectively blocks transcription factor activity and ensures differentiated tissue integrity.

The developmental role of histone H3K9 methylation (H3K9me), which typifies heterochromatin, remains unclear. In Caenorhabditis elegans, loss of H3K9me leads to a highly divergent upregulation of genes with tissue and developmental-stage specificity. During development H3K9me is lost from differentiated cell type-specific genes and gained at genes expressed in earlier developmental stages or other tissues. The continuous deposition of H3K9me2 by the SETDB1 homolog MET-2 after terminal differentiation is necessary to maintain repression. In differentiated tissues, H3K9me ensures silencing by restricting the activity of a defined set of transcription factors at promoters and enhancers. Increased chromatin accessibility following the loss of H3K9me is neither sufficient nor necessary to drive transcription. Increased ATAC-seq signal and gene expression correlate at a subset of loci positioned away from the nuclear envelope, while derepressed genes at the nuclear periphery remain poorly accessible despite being transcribed. In conclusion, H3K9me deposition can confer tissue-specific gene expression and maintain the integrity of terminally differentiated muscle by restricting transcription factor activity.

Gene expression oscillations in C. elegans underlie a new developmental clock.

During C. elegans larval development, thousands of genes, accounting for >20% of the transcriptome, exhibit oscillatory expression with large amplitudes. The time of peaking varies for different genes, but expression generally peaks once per larval stage, with both the oscillation period and larval stage duration varying in concert with temperature. This and other evidence support the existence of a gene expression oscillator that functions as a developmental clock. In this article, we review what is known about the biology, architecture and possible mechanisms of this clock. We compare it to other oscillators, and highlight tools and approaches suited to its study. Finally, we point out implications of these wide-spread and dynamic changes of gene expression on any type of gene expression profiling experiment in C. elegans larvae and how such experiments need to be controlled.

Protease-mediated processing of Argonaute proteins controls small RNA association.

Small RNA pathways defend the germlines of animals against selfish genetic elements, yet pathway activities need to be contained to prevent silencing of self genes. Here, we reveal a proteolytic mechanism that controls endogenous small interfering (22G) RNA activity in the Caenorhabditis elegans germline to protect genome integrity and maintain fertility. We find that DPF-3, a P-granule-localized N-terminal dipeptidase orthologous to mammalian dipeptidyl peptidase (DPP) 8/9, processes the unusually proline-rich N termini of WAGO-1 and WAGO-3 Argonaute (Ago) proteins. Without DPF-3 activity, these WAGO proteins lose their proper complement of 22G RNAs. Desilencing of repeat-containing and transposon-derived transcripts, DNA damage, and acute sterility ensue. These phenotypes are recapitulated when WAGO-1 and WAGO-3 are rendered resistant to DPF-3-mediated processing, identifying them as critical substrates of DPF-3. We conclude that N-terminal processing of Ago proteins regulates their activity and promotes silencing of selfish genetic elements by ensuring Ago association with appropriate small RNAs.

Developmental function and state transitions of a gene expression oscillator in Caenorhabditis elegans.

Gene expression oscillators can structure biological events temporally and spatially. Different biological functions benefit from distinct oscillator properties. Thus, finite developmental processes rely on oscillators that start and stop at specific times, a poorly understood behavior. Here, we have characterized a massive gene expression oscillator comprising > 3,700 genes in Caenorhabditis elegans larvae. We report that oscillations initiate in embryos, arrest transiently after hatching and in response to perturbation, and cease in adults. Experimental observation of the transitions between oscillatory and non-oscillatory states at high temporal resolution reveals an oscillator operating near a Saddle Node on Invariant Cycle (SNIC) bifurcation. These findings constrain the architecture and mathematical models that can represent this oscillator. They also reveal that oscillator arrests occur reproducibly in a specific phase. Since we find oscillations to be coupled to developmental processes, including molting, this characteristic of SNIC bifurcations endows the oscillator with the potential to halt larval development at defined intervals, and thereby execute a developmental checkpoint function.

A branched heterochronic pathway directs juvenile-to-adult transition through two LIN-29 isoforms.

Robust organismal development relies on temporal coordination of disparate physiological processes. In Caenorhabditis elegans, the heterochronic pathway controls a timely juvenile-to-adult (J/A) transition. This regulatory cascade of conserved proteins and small RNAs culminates in accumulation of the transcription factor LIN-29, which triggers coordinated execution of transition events. We report that two LIN-29 isoforms fulfill distinct functions. Functional specialization is a consequence of distinct isoform expression patterns, not protein sequence, and we propose that distinct LIN-29 dose sensitivities of the individual J/A transition events help to ensure their temporal ordering. We demonstrate that unique isoform expression patterns are generated by the activities of LIN-41 for lin-29a, and of HBL-1 for lin-29b, whereas the RNA-binding protein LIN-28 coordinates LIN-29 isoform activity, in part by regulating both hbl-1 and lin-41. Our findings reveal that coordinated transition from juvenile to adult involves branching of a linear pathway to achieve timely control of multiple events.

The RNA hairpin binder TRIM71 modulates alternative splicing by repressing MBNL1.

TRIM71/LIN-41, a phylogenetically conserved regulator of development, controls stem cell fates. Mammalian TRIM71 exhibits both RNA-binding and protein ubiquitylation activities, but the functional contribution of either activity and relevant primary targets remain poorly understood. Here, we demonstrate that TRIM71 shapes the transcriptome of mouse embryonic stem cells (mESCs) predominantly through its RNA-binding activity. We reveal that TRIM71 binds targets through 3' untranslated region (UTR) hairpin motifs and that it acts predominantly by target degradation. TRIM71 mutations implicated in etiogenesis of human congenital hydrocephalus impair target silencing. We identify a set of primary targets consistently regulated in various human and mouse cell lines, including MBNL1 (Muscleblind-like protein 1). MBNL1 promotes cell differentiation through regulation of alternative splicing, and we demonstrate that TRIM71 promotes embryonic splicing patterns through MBNL1 repression. Hence, repression of MBNL1-dependent alternative splicing may contribute to TRIM71's function in regulating stem cell fates.

let-7 coordinates the transition to adulthood through a single primary and four secondary targets.

The juvenile-to-adult (J/A) transition, or puberty, is a period of extensive changes of animal body morphology and function. The onset of puberty is genetically controlled, and the let-7 miRNA temporally regulates J/A transition events in nematodes and mammals. Here, we uncover the targets and downstream pathways through which Caenorhabditis elegans let-7 controls male and female sexual organ morphogenesis and skin progenitor cell fates. We find that let-7 directs all three processes by silencing a single target, the post-transcriptional regulator lin-41 In turn, the RNA-binding protein LIN41/TRIM71 regulates these processes by silencing only four target mRNAs. Thus, by silencing LIN41, let-7 activates LIN-29a and MAB-10 (an early growth response-type transcription factor and its NAB1/2-orthologous cofactor, respectively) to terminate progenitor cell self-renewal and to promote vulval integrity. By contrast, let-7 promotes development of the male sexual organ by up-regulating DMD-3 and MAB-3, two Doublesex/MAB-3 domain-containing transcription factors. Our results provide mechanistic insight into how a linear chain of post-transcriptional regulators diverges in the control of a small set of transcriptional regulators to achieve a coordinated J/A transition.

Timing mechanism of sexually dimorphic nervous system differentiation.

The molecular mechanisms that control the timing of sexual differentiation in the brain are poorly understood. We found that the timing of sexually dimorphic differentiation of postmitotic, sex-shared neurons in the nervous system of the Caenorhabditis elegans male is controlled by the temporally regulated miRNA let-7 and its target lin-41, a translational regulator. lin-41 acts through lin-29a, an isoform of a conserved Zn finger transcription factor, expressed in a subset of sex-shared neurons only in the male. Ectopic lin-29a is sufficient to impose male-specific features at earlier stages of development and in the opposite sex. The temporal, sexual and spatial specificity of lin-29a expression is controlled intersectionally through the lin-28/let-7/lin-41 heterochronic pathway, sex chromosome configuration and neuron-type-specific terminal selector transcription factors. Two Doublesex-like transcription factors represent additional sex- and neuron-type specific targets of LIN-41 and are regulated in a similar intersectional manner.

An interplay of miRNA abundance and target site architecture determines miRNA activity and specificity.

MicroRNAs often occur in families whose members share an identical 5' terminal 'seed' sequence. The seed is a major determinant of miRNA activity, and family members are thought to act redundantly on target mRNAs with perfect seed matches, i.e. sequences complementary to the seed. However, recently sequences outside the seed were reported to promote silencing by individual miRNA family members. Here, we examine this concept and the importance of miRNA specificity for the robustness of developmental gene control. Using the let-7 miRNA family in Caenorhabditis elegans, we find that seed match imperfections can increase specificity by requiring extensive pairing outside the miRNA seed region for efficient silencing and that such specificity is needed for faithful worm development. In addition, for some target site architectures, elevated miRNA levels can compensate for a lack of complementarity outside the seed. Thus, some target sites require higher miRNA concentration for silencing than others, contrasting with a traditional binary distinction between functional and non-functional sites. We conclude that changing miRNA concentrations can alter cellular miRNA target repertoires. This diversifies possible biological outcomes of miRNA-mediated gene regulation and stresses the importance of target validation under physiological conditions to understand miRNA functions in vivo.

Evolutionary plasticity of the NHL domain underlies distinct solutions to RNA recognition.

RNA-binding proteins regulate all aspects of RNA metabolism. Their association with RNA is mediated by RNA-binding domains, of which many remain uncharacterized. A recently reported example is the NHL domain, found in prominent regulators of cellular plasticity like the C. elegans LIN-41. Here we employ an integrative approach to dissect the RNA specificity of LIN-41. Using computational analysis, structural biology, and in vivo studies in worms and human cells, we find that a positively charged pocket, specific to the NHL domain of LIN-41 and its homologs (collectively LIN41), recognizes a stem-loop RNA element, whose shape determines the binding specificity. Surprisingly, the mechanism of RNA recognition by LIN41 is drastically different from that of its more distant relative, the fly Brat. Our phylogenetic analysis suggests that this reflects a rapid evolution of the domain, presenting an interesting example of a conserved protein fold that acquired completely different solutions to RNA recognition.

Two distinct transcription termination modes dictated by promoters.

Transcription termination determines the ends of transcriptional units and thereby ensures the integrity of the transcriptome and faithful gene regulation. Studies in yeast and human cells have identified the exoribonuclease XRN2 as a key termination factor for protein-coding genes. Here we performed a genome-wide investigation of RNA polymerase II (Pol II) transcription termination in XRN2-deficient Caenorhabditis elegans and observed two distinct modes of termination. Although a subset of genes requires XRN2, termination of other genes appears both independent of, and refractory to, XRN2. XRN2 independence is not merely a consequence of failure to recruit XRN2, since XRN2 is present on-and promotes Pol II accumulation near the polyadenylation sites of-both gene classes. Unexpectedly, promoters instruct the choice of termination mode, but XRN2-independent termination additionally requires a compatible region downstream from the 3' end cleavage site. Hence, different termination mechanisms may work with different configurations of Pol II complexes dictated by promoters.

Bayesian prediction of RNA translation from ribosome profiling.

Ribosome profiling via high-throughput sequencing (ribo-seq) is a promising new technique for characterizing the occupancy of ribosomes on messenger RNA (mRNA) at base-pair resolution. The ribosome is responsible for translating mRNA into proteins, so information about its occupancy offers a detailed view of ribosome density and position which could be used to discover new translated open reading frames (ORFs), among other things. In this work, we propose Rp-Bp, an unsupervised Bayesian approach to predict translated ORFs from ribosome profiles. We use state-of-the-art Markov chain Monte Carlo techniques to estimate posterior distributions of the likelihood of translation of each ORF. Hence, an important feature of Rp-Bp is its ability to incorporate and propagate uncertainty in the prediction process. A second novel contribution is automatic Bayesian selection of read lengths and ribosome P-site offsets (BPPS). We empirically demonstrate that our read length selection technique modestly improves sensitivity by identifying more canonical and non-canonical ORFs. Proteomics- and quantitative translation initiation sequencing-based validation verifies the high quality of all of the predictions. Experimental comparison shows that Rp-Bp results in more peptide identifications and proteomics-validated ORF predictions compared to another recent tool for translation prediction.