A multiparameter network reveals extensive divergence between C. elegans bHLH transcription factors

Simple nutrients bypass the requirement for HLH-30 in coupling lysosomal nutrient sensing to survival

Lysosomes are ubiquitous acidified organelles that degrade intracellular and extracellular material trafficked via multiple pathways. Lysosomes also sense cellular nutrient levels to regulate target of rapamycin (TOR) kinase, a signaling enzyme that drives growth and suppresses activity of the MiT/TFE family of transcription factors that control biogenesis of lysosomes. In this study, we subjected worms lacking basic helix–loop–helix transcription factor 30 (hlh-30), the Caenorhabditis elegans MiT/TFE ortholog, to starvation followed by refeeding to understand how this pathway regulates survival with variable nutrient supply. Loss of HLH-30 markedly impaired survival in starved larval worms and recovery upon refeeding bacteria. Remarkably, provision of simple nutrients in a completely defined medium (C. elegans maintenance medium [CeMM]), specifically glucose and linoleic acid, restored lysosomal acidification, TOR activation, and survival with refeeding despite the absence of HLH-30. Worms deficient in lysosomal lipase 2 (lipl-2), a lysosomal enzyme that is transcriptionally up-regulated in starvation in an HLH-30–dependent manner, also demonstrated increased mortality with starvation–refeeding that was partially rescued with glucose, suggesting a critical role for LIPL-2 in lipid metabolism under starvation. CeMM induced transcription of vacuolar proton pump subunits in hlh-30 mutant worms, and knockdown of vacuolar H⁺-ATPase 12 (vha-12) and its upstream regulator, nuclear hormone receptor 31 (nhr-31), abolished the rescue with CeMM. Loss of Ras-related GTP binding protein C homolog 1 RAGC-1, the ortholog for mammalian RagC/D GTPases, conferred starvation–refeeding lethality, and RAGC-1 overexpression was sufficient to rescue starved hlh-30 mutant worms, demonstrating a critical need for TOR activation with refeeding. These results show that HLH-30 activation is critical for sustaining survival during starvation–refeeding stress via regulating TOR. Glucose and linoleic acid bypass the requirement for HLH-30 in coupling lysosome nutrient sensing to survival.

Lysosomes play a central role in coupling the nutrient state of the cell to growth and survival decisions. This study uncovers a critical role for HLH-30, the nematode ortholog of the mammalian MiT/TFE family of master regulators of lysosome biogenesis, in survival under starvation and refeeding conditions. Refeeding simple nutrients bypasses the requirement for HLH-30 to permit survival.

Transcription factor Ptf1a in development, diseases and reprogramming

The transcription factor Ptf1a is a crucial helix-loop-helix (bHLH) protein selectively expressed in the pancreas, retina, spinal cord, brain, and enteric nervous system. Ptf1a is preferably assembled into a transcription trimeric complex PTF1 with an E protein and Rbpj (or Rbpjl). In pancreatic development, Ptf1a is indispensable in controlling the expansion of multipotent progenitor cells as well as the specification and maintenance of the acinar cells. In neural tissues, Ptf1a is transiently expressed in the post-mitotic cells and specifies the inhibitory neuronal cell fates, mostly mediated by downstream genes such as Tfap2a/b and Prdm13. Mutations in the coding and non-coding regulatory sequences resulting in Ptf1a gain- or loss-of-function are associated with genetic diseases such as pancreatic and cerebellar agenesis in the rodent and human. Surprisingly, Ptf1a alone is sufficient to reprogram mouse or human fibroblasts into tripotential neural stem cells. Its pleiotropic functions in many biological processes remain to be deciphered in the future.

Global transcriptional regulation of innate immunity by ATF-7 in C. elegans

The nematode Caenorhabditis elegans has emerged as a genetically tractable animal host in which to study evolutionarily conserved mechanisms of innate immune signaling. We previously showed that the PMK-1 p38 mitogen-activated protein kinase (MAPK) pathway regulates innate immunity of C. elegans through phosphorylation of the CREB/ATF bZIP transcription factor, ATF-7. Here, we have undertaken a genomic analysis of the transcriptional response of C. elegans to infection by Pseudomonas aeruginosa, combining genome-wide expression analysis by RNA-seq with ATF-7 chromatin immunoprecipitation followed by sequencing (ChIP-Seq). We observe that PMK-1-ATF-7 activity regulates a majority of all genes induced by pathogen infection, and observe ATF-7 occupancy in regulatory regions of pathogen-induced genes in a PMK-1-dependent manner. Moreover, functional analysis of a subset of these ATF-7-regulated pathogen-induced target genes supports a direct role for this transcriptional response in host defense. The genome-wide regulation through PMK-1- ATF-7 signaling reveals a striking level of control over the innate immune response to infection through a single transcriptional regulator.

The Human Transcription Factors

Transcription factors (TFs) recognize specific DNA sequences to control chromatin and transcription, forming a complex system that guides expression of the genome. Despite keen interest in understanding how TFs control gene expression, it remains challenging to determine how the precise genomic binding sites of TFs are specified and how TF binding ultimately relates to regulation of transcription. This review considers how TFs are identified and functionally characterized, principally through the lens of a catalog of over 1,600 likely human TFs and binding motifs for two-thirds of them. Major classes of human TFs differ markedly in their evolutionary trajectories and expression patterns, underscoring distinct functions. TFs likewise underlie many different aspects of human physiology, disease, and variation, highlighting the importance of continued effort to understand TF-mediated gene regulation.

A tissue-specific enhancer of the C. elegans nhr-67/tailless gene drives coordinated expression in uterine stem cells and the differentiated anchor cell

The nhr-67 nuclear receptor gene of Caenorhabditis elegans encodes the ortholog of the Drosophila tailless and vertebrate Tlx genes. In C. elegans, nhr-67 plays multiple roles in the development of the uterus during L2 and L3 larval stages. Four pre-VU cells are born in the L2 stage and form the precursor complement for the ventral surface of the mature uterus. One of the four pre-VU cells becomes the anchor cell (AC), which exits the cell cycle and differentiates, while the remaining three VU cells serve as stem cells that populate the ventral uterus. The nhr-67 gene functions in the development of both VU cell lineages and AC differentiation. Hypomorphic mutations in nhr-67 identify a 276bp region of the distal promoter that is sufficient to activate nhr-67 expression in pre-VU cells and the AC. The 276bp region includes 8 conserved potential cis-acting sites, including two E boxes and a nuclear receptor binding site. Mutational analysis demonstrates that the two E boxes are required for expression of nhr-67 in uterine precursor cells. The E/daughterless ortholog HLH-2 binds these sites as a homodimer, thus playing a central role in activating nhr-67 expression in the uterine precursors. At least two other binding activities, one of which may be the nhr-25/Ftz-F1 nuclear receptor transcription factor, also contribute to uterine precursor cell expression. The organization of the nhr-67 uterine precursor enhancer is compared to similar conserved enhancers in the egl-43, lag-2, and lin-3 genes, which contain the same HLH-2-binding E boxes and are similarly expressed in both pre-VU cells and the AC. This basic regulatory module allows the coordinated expression of at least four genes. Expression of genes in different cells that must coordinate to form a mature organ is driven by a shared set of promoter elements, which integrate multiple transcription factor inputs.

Gateway-Compatible Yeast One-Hybrid and Two-Hybrid Assays

In the first section of this introduction, we provide background information for yeast two-hybrid (Y2H) assays that provide a genetic method for the identification and analysis of binary protein-protein interactions and that are complementary to biochemical methods such as immunoprecipitation. In the second section, we discuss yeast one-hybrid (Y1H) assays that provide a "gene-centered" (DNA-to-protein) genetic method to identify and study protein-DNA interactions between cis-regulatory elements and transcription factors (TFs). This method is complementary to "TF-centered" (protein-to-DNA) biochemical methods such as chromatin immunoprecipitation.

Synergistic Binding of bHLH Transcription Factors to the Promoter of the Maize NADP-ME Gene Used in C4 Photosynthesis Is Based on an Ancient Code Found in the Ancestral C3 State

C4 photosynthesis has evolved repeatedly from the ancestral C3 state to generate a carbon concentrating mechanism that increases photosynthetic efficiency. This specialized form of photosynthesis is particularly common in the PACMAD clade of grasses, and is used by many of the world's most productive crops. The C4 cycle is accomplished through cell-type-specific accumulation of enzymes but cis-elements and transcription factors controlling C4 photosynthesis remain largely unknown. Using the NADP-Malic Enzyme (NADP-ME) gene as a model we tested whether mechanisms impacting on transcription in C4 plants evolved from ancestral components found in C3 species. Two basic Helix-Loop-Helix (bHLH) transcription factors, ZmbHLH128 and ZmbHLH129, were shown to bind the C4NADP-ME promoter from maize. These proteins form heterodimers and ZmbHLH129 impairs trans-activation by ZmbHLH128. Electrophoretic mobility shift assays indicate that a pair of cis-elements separated by a seven base pair spacer synergistically bind either ZmbHLH128 or ZmbHLH129. This pair of cis-elements is found in both C3 and C4 Panicoid grass species of the PACMAD clade. Our analysis is consistent with this cis-element pair originating from a single motif present in the ancestral C3 state. We conclude that C4 photosynthesis has co-opted an ancient C3 regulatory code built on G-box recognition by bHLH to regulate the NADP-ME gene. More broadly, our findings also contribute to the understanding of gene regulatory networks controlling C4 photosynthesis.

Unconventional function of an Achaete-Scute homolog as a terminal selector of nociceptive neuron identity

Proneural genes are among the most early-acting genes in nervous system development, instructing blast cells to commit to a neuronal fate. Drosophila Atonal and Achaete-Scute complex (AS-C) genes, as well as their vertebrate orthologs, are basic helix-loop-helix (bHLH) transcription factors with such proneural activity. We show here that a C. elegans AS-C homolog, hlh-4, functions in a fundamentally different manner. In the embryonic, larval, and adult nervous systems, hlh-4 is expressed exclusively in a single nociceptive neuron class, ADL, and its expression in ADL is maintained via transcriptional autoregulation throughout the life of the animal. However, in hlh-4 null mutants, the ADL neuron is generated and still appears neuronal in overall morphology and expression of panneuronal and pansensory features. Rather than acting as a proneural gene, we find that hlh-4 is required for the ADL neuron to function properly, to adopt its correct morphology, to express its unusually large repertoire of olfactory receptor-encoding genes, and to express other known features of terminal ADL identity, including neurotransmitter phenotype, neuropeptides, ion channels, and electrical synapse proteins. hlh-4 is sufficient to induce ADL identity features upon ectopic expression in other neuron types. The expression of ADL terminal identity features is directly controlled by HLH-4 via a phylogenetically conserved E-box motif, which, through bioinformatic analysis, we find to constitute a predictive feature of ADL-expressed terminal identity markers. The lineage that produces the ADL neuron was previously shown to require the conventional, transient proneural activity of another AS-C homolog, hlh-14, demonstrating sequential activities of distinct AS-C-type bHLH genes in neuronal specification. Taken together, we have defined here an unconventional function of an AS-C-type bHLH gene as a terminal selector of neuronal identity and we speculate that such function could be reflective of an ancestral function of an "ur-" bHLH gene.

Intracellular Trafficking and Persistence of Acinetobacter baumannii Requires Transcription Factor EB

Adhesion is an initial and important step in Acinetobacter baumannii infections. However, the mechanism of entrance and persistence inside host cells is unclear and remains to be understood. In this study, we report that, in addition to its known role in host defense against Gram-positive bacterial infection, TFEB also plays an important role in the intracellular trafficking of A. baumannii in host cells. TFEB was activated shortly after A. baumannii infection and is required for its persistence within host cells. Additionally, using the C. elegans infection model by A. baumannii, the TFEB orthologue HLH-30 was required for survival of the nematode to infection, although nuclear translocation of HLH-30 was not required.

Acinetobacter baumannii is a significant human pathogen associated with hospital-acquired infections. While adhesion, an initial and important step in A. baumannii infection, is well characterized, the intracellular trafficking of this pathogen inside host cells remains poorly studied. Here, we demonstrate that transcription factor EB (TFEB) is activated after A. baumannii infection of human lung epithelial cells (A549). We also show that TFEB is required for the invasion and persistence inside A549 cells. Consequently, lysosomal biogenesis and autophagy activation were observed after TFEB activation which could increase the death of A549 cells. In addition, using the Caenorhabditis elegans infection model by A. baumannii, the TFEB orthologue HLH-30 was required for survival of the nematode to infection, although nuclear translocation of HLH-30 was not required. These results identify TFEB as a conserved key factor in the pathogenesis of A. baumannii.

IMPORTANCE Adhesion is an initial and important step in Acinetobacter baumannii infections. However, the mechanism of entrance and persistence inside host cells is unclear and remains to be understood. In this study, we report that, in addition to its known role in host defense against Gram-positive bacterial infection, TFEB also plays an important role in the intracellular trafficking of A. baumannii in host cells. TFEB was activated shortly after A. baumannii infection and is required for its persistence within host cells. Additionally, using the C. elegans infection model by A. baumannii, the TFEB orthologue HLH-30 was required for survival of the nematode to infection, although nuclear translocation of HLH-30 was not required.

Systems biology of embryonic development: Prospects for a complete understanding of the Caenorhabditis elegans embryo

The convergence of developmental biology and modern genomics tools brings the potential for a comprehensive understanding of developmental systems. This is especially true for the Caenorhabditis elegans embryo because its small size, invariant developmental lineage, and powerful genetic and genomic tools provide the prospect of a cellular resolution understanding of messenger RNA (mRNA) expression and regulation across the organism. We describe here how a systems biology framework might allow large-scale determination of the embryonic regulatory relationships encoded in the C. elegans genome. This framework consists of two broad steps: (a) defining the "parts list"-all genes expressed in all cells at each time during development and (b) iterative steps of computational modeling and refinement of these models by experimental perturbation. Substantial progress has been made towards defining the parts list through imaging methods such as large-scale green fluorescent protein (GFP) reporter analysis. Imaging results are now being augmented by high-resolution transcriptome methods such as single-cell RNA sequencing, and it is likely the complete expression patterns of all genes across the embryo will be known within the next few years. In contrast, the modeling and perturbation experiments performed so far have focused largely on individual cell types or genes, and improved methods will be needed to expand them to the full genome and organism. This emerging comprehensive map of embryonic expression and regulatory function will provide a powerful resource for developmental biologists, and would also allow scientists to ask questions not accessible without a comprehensive picture. This article is categorized under: Invertebrate Organogenesis > Worms Technologies > Analysis of the Transcriptome Gene Expression and Transcriptional Hierarchies > Gene Networks and Genomics.

Gateway Recombinational Cloning

The Gateway recombinatorial cloning system was developed for cloning multiple DNA fragments in parallel (e.g., in 96-well formats) in a standardized manner using the same enzymes. Gateway cloning is based on the highly specific integration and excision reactions of bacteriophage λ into and out of the Escherichia coli genome. Because the sites of recombination ("att" sites) are much longer (25-242 bp) than restriction sites, they are extremely unlikely to occur by chance in DNA fragments. Therefore, the same recombination enzyme can be used to robustly clone many different fragments of variable size in parallel reactions.

Zebrafish Transcription Factor ORFeome for Gene Discovery and Regulatory Network Elucidation

The completion of the zebrafish genome sequence and advances in miniaturization and multiplexing were essential to the creation of techniques such as RNA-seq, ChIP-seq, and high-throughput behavioral and chemical screens. Multiplexing was also instrumental in the recent enhancement of the classic yeast one-hybrid interaction techniques to provide unprecedented discovery capabilities for protein-DNA interactions. Unfortunately its use for zebrafish research is currently hampered by the lack of an open reading frame (ORF) clone collection. As a first step toward a complete collection, we describe a small library of transcriptional regulatory proteins comprising 142 ORFs and its potential applications.

The MXL-3/SBP-1 Axis Is Responsible for Glucose-Dependent Fat Accumulation in C. elegans

Chronic exposure to elevated glucose levels leads to fatty acid accumulation, which promotes the development of metabolic diseases such as obesity and type 2 diabetes. MXL-3 is a conserved transcriptional factor that modulates the inhibition of lipolysis in Caenorhabditis elegans. However, the role of MXL-3 in lipid metabolism during nutrient excess remains unknown. We hypothesized that inhibition of MXL-3 prevents glucose-dependent fat accumulation. Nematodes from wild-type N2, MXL-3::GFP and sbp-1 or mxl-3 null strains were grown on standard, high glucose or high glucose plus metformin plates for 24 h. Using laser-scanning confocal microscopy, we monitored the glucose-induced activation of MXL-3 labeled with GFP (MXL-3::GFP). Lipid levels were determined by Oil Red O (ORO) staining and gas chromatography/mass spectrometry, and gene expression was assessed by qRT-PCR. We found that high glucose activated MXL-3 by increasing its rate of nuclear entry, which in turn increased lipid levels via sterol regulatory element-binding protein (SBP-1). This activated critical genes that synthesize long chain unsaturated fatty acids (MUFAs and PUFAs) and repress lipolytic genes. Interestingly, the anti-diabetic drug metformin inhibited MXL-3 activation and subsequently prevented glucose-dependent fat accumulation. These findings highlight the importance of the MXL-3/SBP-1 axis in the regulation of lipid metabolism during nutritional excess and provide new insight into the mechanism by which metformin prevents lipid accumulation. This study also suggests that inhibition of MXL-3 may serve as a potential target for the treatment of chronic metabolic diseases, including obesity, type 2 diabetes, and cardiovascular disease.

Identification of Human Lineage-Specific Transcriptional Coregulators Enabled by a Glossary of Binding Modules and Tunable Genomic Backgrounds

Transcription factors (TFs) control cellular processes by binding specific DNA motifs to modulate gene expression. Motif enrichment analysis of regulatory regions can identify direct and indirect TF binding sites. Here, we created a glossary of 108 non-redundant TF-8mer "modules" of shared specificity for 671 metazoan TFs from publicly available and new universal protein binding microarray data. Analysis of 239 ENCODE TF chromatin immunoprecipitation sequencing datasets and associated RNA sequencing profiles suggest the 8mer modules are more precise than position weight matrices in identifying indirect binding motifs and their associated tethering TFs. We also developed GENRE (genomically equivalent negative regions), a tunable tool for construction of matched genomic background sequences for analysis of regulatory regions. GENRE outperformed four state-of-the-art approaches to background sequence construction. We used our TF-8mer glossary and GENRE in the analysis of the indirect binding motifs for the co-occurrence of tethering factors, suggesting novel TF-TF interactions. We anticipate that these tools will aid in elucidating tissue-specific gene-regulatory programs.

Multi-omics Analyses of Starvation Responses Reveal a Central Role for Lipoprotein Metabolism in Acute Starvation Survival in C. elegans

Starvation causes comprehensive metabolic changes, which are still not fully understood. Here, we used quantitative proteomics and RNA sequencing to examine the temporal starvation responses in wild-type Caenorhabditis elegans and animals lacking the transcription factor HLH-30. Our findings show that starvation alters the abundance of hundreds of proteins and mRNAs in a temporal manner, many of which are involved in central metabolic pathways, including lipoprotein metabolism. We demonstrate that premature death of hlh-30 animals under starvation can be prevented by knockdown of either vit-1 or vit-5, encoding two different lipoproteins. We further show that the size and number of intestinal lipid droplets under starvation are altered in hlh-30 animals, which can be rescued by knockdown of vit-1. Taken together, this indicates that survival of hlh-30 animals under starvation is closely linked to regulation of intestinal lipid stores. We provide the most detailed poly-omic analysis of starvation responses to date, which serves as a resource for further mechanistic studies of starvation.

A bHLH Code for Sexually Dimorphic Form and Function of the C. elegans Somatic Gonad

How sexually dimorphic gonads are generated is a fundamental question at the interface of developmental and evolutionary biology [1-3]. In C. elegans, sexual dimorphism in gonad form and function largely originates in different apportionment of roles to three regulatory cells of the somatic gonad primordium in young larvae. Their essential roles include leading gonad arm outgrowth, serving as the germline niche, connecting to epithelial openings, and organizing reproductive organ development. The development and function of the regulatory cells in both sexes requires the basic-helix-loop-helix (bHLH) transcription factor HLH-2, the sole ortholog of the E proteins mammalian E2A and Drosophila Daughterless [4-8], yet how they adopt different fates to execute their different roles has been unknown. Here, we show that each regulatory cell expresses a distinct complement of bHLH-encoding genes-and therefore distinct HLH-2:bHLH dimers-and formulate a "bHLH code" hypothesis for regulatory cell identity. We support this hypothesis by showing that the bHLH gene complement is both necessary and sufficient to confer particular regulatory cell fates. Strikingly, prospective regulatory cells can be directly reprogrammed into other regulatory cell types simply by loss or ectopic expression of bHLH genes, and male-to-female and female-to-male transformations indicate that the code is instructive for sexual dimorphism. The bHLH code appears to be embedded in a bow-tie regulatory architecture [9, 10], wherein sexual, positional, temporal, and lineage inputs connect through bHLH genes to diverse outputs for terminal features and provides a plausible mechanism for the evolutionary plasticity of gonad form seen in nematodes [11-15].

Destabilization of the TWIST1/E12 complex dimerization following the R154P point-mutation of TWIST1: an in silico approach

Background

The bHLH transcription factor TWIST1 plays a key role in the embryonic development and in tumorigenesis. Some loss-of-function mutations of the TWIST1 gene have been shown to cause an autosomal dominant craniosynostosis, known as the Saethre-Chotzen syndrome (SCS). Although the functional impacts of many TWIST1 mutations have been experimentally reported, little is known on the molecular mechanisms underlying their loss-of-function. In a previous study, we highlighted the predictive value of in silico molecular dynamics (MD) simulations in deciphering the molecular function of TWIST1 residues.

Results

Here, since the substitution of the arginine 154 amino acid by a glycine residue (R154G) is responsible for the SCS phenotype and the substitution of arginine 154 by a proline experimentally decreases the dimerizing ability of TWIST1, we investigated the molecular impact of this point mutation using MD approaches. Consistently, MD simulations highlighted a clear decrease in the stability of the α-helix during the dimerization of the mutated R154P TWIST1/E12 dimer compared to the wild-type TE complex, which was further confirmed in vitro using immunoassays.

Conclusions

Our study demonstrates that MD simulations provide a structural explanation for the loss-of-function associated with the SCS TWIST1 mutation and provides a proof of concept of the predictive value of these MD simulations. This in silico methodology could be used to determine reliable pharmacophore sites, leading to the application of docking approaches in order to identify specific inhibitors of TWIST1 complexes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12900-017-0076-x) contains supplementary material, which is available to authorized users.

PRIMA: a gene-centered, RNA-to-protein method for mapping RNA-protein interactions

Interactions between RNA binding proteins (RBPs) and mRNAs are critical to post-transcriptional gene regulation. Eukaryotic genomes encode thousands of mRNAs and hundreds of RBPs. However, in contrast to interactions between transcription factors (TFs) and DNA, the interactome between RBPs and RNA has been explored for only a small number of proteins and RNAs. This is largely because the focus has been on using 'protein-centered' (RBP-to-RNA) interaction mapping methods that identify the RNAs with which an individual RBP interacts. While powerful, these methods cannot as of yet be applied to the entire RBPome. Moreover, it may be desirable for a researcher to identify the repertoire of RBPs that can interact with an mRNA of interest-in a 'gene-centered' manner-yet few such techniques are available. Here, we present Protein-RNA Interaction Mapping Assay (PRIMA) with which an RNA 'bait' can be tested versus multiple RBP 'preys' in a single experiment. PRIMA is a translation-based assay that examines interactions in the yeast cytoplasm, the cellular location of mRNA translation. We show that PRIMA can be used with small RNA elements, as well as with full-length Caenorhabditis elegans 3' UTRs. PRIMA faithfully recapitulated numerous well-characterized RNA-RBP interactions and also identified novel interactions, some of which were confirmed in vivo. We envision that PRIMA will provide a complementary tool to expand the depth and scale with which the RNA-RBP interactome can be explored.

The CDK8 Complex and Proneural Proteins Together Drive Neurogenesis from a Mesodermal Lineage

At least some animal species can generate neurons from mesoderm or endoderm, but the underlying mechanisms remain unknown. We screened for C. elegans mutants in which the presumptive mesoderm-derived I4 neuron adopts a muscle-like cell fate. From this screen, we identified HLH-3, the C. elegans homolog of a mammalian proneural protein (Ascl1) used for in vitro neuronal reprogramming, as required for efficient I4 neurogenesis. We discovered that the CDK-8 Mediator kinase module acts together with a second proneural protein, HLH-2, and in parallel to HLH-3 to promote I4 neurogenesis. Genetic analysis revealed that CDK-8 most likely promotes I4 neurogenesis by inhibiting the CDK-7/CYH-1 (CDK7/cyclin H) kinase module of the transcription initiation factor TFIIH. Ectopic expression of HLH-2 and HLH-3 together promoted expression of neuronal features in non-neuronal cells. These findings reveal that the Mediator CDK8 kinase module can promote non-ectodermal neurogenesis and suggest that inhibiting CDK7/cyclin H might similarly promote neurogenesis.

Systematic dissection of genomic features determining transcription factor binding and enhancer function

Enhancers regulate gene expression through the binding of sequence-specific transcription factors (TFs) to cognate motifs. Various features influence TF binding and enhancer function-including the chromatin state of the genomic locus, the affinities of the binding site, the activity of the bound TFs, and interactions among TFs. However, the precise nature and relative contributions of these features remain unclear. Here, we used massively parallel reporter assays (MPRAs) involving 32,115 natural and synthetic enhancers, together with high-throughput in vivo binding assays, to systematically dissect the contribution of each of these features to the binding and activity of genomic regulatory elements that contain motifs for PPARγ, a TF that serves as a key regulator of adipogenesis. We show that distinct sets of features govern PPARγ binding vs. enhancer activity. PPARγ binding is largely governed by the affinity of the specific motif site and higher-order features of the larger genomic locus, such as chromatin accessibility. In contrast, the enhancer activity of PPARγ binding sites depends on varying contributions from dozens of TFs in the immediate vicinity, including interactions between combinations of these TFs. Different pairs of motifs follow different interaction rules, including subadditive, additive, and superadditive interactions among specific classes of TFs, with both spatially constrained and flexible grammars. Our results provide a paradigm for the systematic characterization of the genomic features underlying regulatory elements, applicable to the design of synthetic regulatory elements or the interpretation of human genetic variation.

Combinatorial bZIP dimers display complex DNA-binding specificity landscapes

How transcription factor dimerization impacts DNA-binding specificity is poorly understood. Guided by protein dimerization properties, we examined DNA binding specificities of 270 human bZIP pairs. DNA interactomes of 80 heterodimers and 22 homodimers revealed that 72% of heterodimer motifs correspond to conjoined half-sites preferred by partnering monomers. Remarkably, the remaining motifs are composed of variably-spaced half-sites (12%) or 'emergent' sites (16%) that cannot be readily inferred from half-site preferences of partnering monomers. These binding sites were biochemically validated by EMSA-FRET analysis and validated in vivo by ChIP-seq data from human cell lines. Focusing on ATF3, we observed distinct cognate site preferences conferred by different bZIP partners, and demonstrated that genome-wide binding of ATF3 is best explained by considering many dimers in which it participates. Importantly, our compendium of bZIP-DNA interactomes predicted bZIP binding to 156 disease associated SNPs, of which only 20 were previously annotated with known bZIP motifs.

C. elegans and its bacterial diet as a model for systems-level understanding of host-microbiota interactions

Resident microbes of the human body, particularly the gut microbiota, provide essential functions for the host, and, therefore, have important roles in human health as well as mitigating disease. It is difficult to study the mechanisms by which the microbiota affect human health, especially at a systems-level, due to heterogeneity of human genomes, the complexity and heterogeneity of the gut microbiota, the challenge of growing these bacteria in the laboratory, and the lack of bacterial genetics in most microbiotal species. In the last few years, the interspecies model of the nematode Caenorhabditis elegans and its bacterial diet has proven powerful for studying host-microbiota interactions, as both the animal and its bacterial diet can be subjected to large-scale and high-throughput genetic screening. The high level of homology between many C. elegans and human genes, as well as extensive similarities between human and C. elegans metabolism, indicates that the findings obtained from this interspecies model may be broadly relevant to understanding how the human microbiota affects physiology and disease. In this review, we summarize recent systems studies on how bacteria interact with C. elegans and affect life history traits.

A lysosome-centered view of nutrient homeostasis

Lysosomes are highly acidic cellular organelles traditionally viewed as sacs of enzymes involved in digesting extracellular or intracellular macromolecules for the regeneration of basic building blocks, cellular housekeeping, or pathogen degradation. Bound by a single lipid bilayer, lysosomes receive their substrates by fusing with endosomes or autophagosomes, or through specialized translocation mechanisms such as chaperone-mediated autophagy or microautophagy. Lysosomes degrade their substrates using up to 60 different soluble hydrolases and release their products either to the cytosol through poorly defined exporting and efflux mechanisms or to the extracellular space by fusing with the plasma membrane. However, it is becoming evident that the role of the lysosome in nutrient homeostasis goes beyond the disposal of waste or the recycling of building blocks. The lysosome is emerging as a signaling hub that can integrate and relay external and internal nutritional information to promote cellular and organismal homeostasis, as well as a major contributor to the processing of energy-dense molecules like glycogen and triglycerides. Here we describe the current knowledge of the nutrient signaling pathways governing lysosomal function, the role of the lysosome in nutrient mobilization, and how lysosomes signal other organelles, distant tissues, and even themselves to ensure energy homeostasis in spite of fluctuations in energy intake. At the same time, we highlight the value of genomics approaches to the past and future discoveries of how the lysosome simultaneously executes and controls cellular homeostasis.

Multiparameter functional diversity of human C2H2 zinc finger proteins

C2H2 zinc finger proteins represent the largest and most enigmatic class of human transcription factors. Their C2H2-ZF arrays are highly variable, indicating that most will have unique DNA binding motifs. However, most of the binding motifs have not been directly determined. In addition, little is known about whether or how these proteins regulate transcription. Most of the ∼700 human C2H2-ZF proteins also contain at least one KRAB, SCAN, BTB, or SET domain, suggesting that they may have common interacting partners and/or effector functions. Here, we report a multifaceted functional analysis of 131 human C2H2-ZF proteins, encompassing DNA binding sites, interacting proteins, and transcriptional response to genetic perturbation. We confirm the expected diversity in DNA binding motifs and genomic binding sites, and provide motif models for 78 previously uncharacterized C2H2-ZF proteins, most of which are unique. Surprisingly, the diversity in protein-protein interactions is nearly as high as diversity in DNA binding motifs: Most C2H2-ZF proteins interact with a unique spectrum of co-activators and co-repressors. Thus, multiparameter diversification likely underlies the evolutionary success of this large class of human proteins.

A gene-centered C. elegans protein-DNA interaction network provides a framework for functional predictions

Transcription factors (TFs) play a central role in controlling spatiotemporal gene expression and the response to environmental cues. A comprehensive understanding of gene regulation requires integrating physical protein–DNA interactions (PDIs) with TF regulatory activity, expression patterns, and phenotypic data. Although great progress has been made in mapping PDIs using chromatin immunoprecipitation, these studies have only characterized ~10% of TFs in any metazoan species. The nematode C. elegans has been widely used to study gene regulation due to its compact genome with short regulatory sequences. Here, we delineated the largest gene‐centered metazoan PDI network to date by examining interactions between 90% of C. elegans TFs and 15% of gene promoters. We used this network as a backbone to predict TF binding sites for 77 TFs, two‐thirds of which are novel, as well as integrate gene expression, protein–protein interaction, and phenotypic data to predict regulatory and biological functions for multiple genes and TFs.

Microsporidia Intracellular Development Relies on Myc Interaction Network Transcription Factors in the Host

Microsporidia are ubiquitous parasites that infect a wide range of animal hosts, and these fungal-related microbes undergo their entire replicative lifecycle inside of host cells. Despite being widespread in the environment and causing medical and agricultural harm, virtually nothing is known about the host factors important to facilitate their growth and development inside of host cells. Here, we perform a genetic screen to identify host transcription factors important for development of the microsporidian pathogen Nematocida parisii inside intestinal cells of its natural host, the nematode Caenorhabditis elegans. Through this screen, we identified the C. elegans Myc family of transcription factors as key host regulators of microsporidia growth and development. The Mad-like transcription factor MDL-1, and the Max-like transcription factors MXL-1 and MXL-2 promote pathogen levels, while the Myc-Mondo-like transcription factor MML-1 inhibits pathogen levels. We used epistasis analysis to show that MDL-1 and MXL-1, which are thought to function as a heterodimer, appear to be acting canonically. In contrast, MXL-2 and MML-1, which are also thought to function as a heterodimer, appear to be acting in separate pathways (noncanonically) in the context of pathogen infection. We also found that both MDL-1::GFP and MML-1::GFP are expressed in intestinal cells during infection. These findings provide novel insight into the host transcription factors that regulate microsporidia development.

Evolution of New cis-Regulatory Motifs Required for Cell-Specific Gene Expression in Caenorhabditis

Patterning of C. elegans vulval cell fates relies on inductive signaling. In this induction event, a single cell, the gonadal anchor cell, secretes LIN-3/EGF and induces three out of six competent precursor cells to acquire a vulval fate. We previously showed that this developmental system is robust to a four-fold variation in lin-3/EGF genetic dose. Here using single-molecule FISH, we find that the mean level of expression of lin-3 in the anchor cell is remarkably conserved. No change in lin-3 expression level could be detected among C. elegans wild isolates and only a low level of change—less than 30%—in the Caenorhabditis genus and in Oscheius tipulae. In C. elegans, lin-3 expression in the anchor cell is known to require three transcription factor binding sites, specifically two E-boxes and a nuclear-hormone-receptor (NHR) binding site. Mutation of any of these three elements in C. elegans results in a dramatic decrease in lin-3 expression. Yet only a single E-box is found in the Drosophilae supergroup of Caenorhabditis species, including C. angaria, while the NHR-binding site likely only evolved at the base of the Elegans group. We find that a transgene from C. angaria bearing a single E-box is sufficient for normal expression in C. elegans. Even a short 58 bp cis-regulatory fragment from C. angaria with this single E-box is able to replace the three transcription factor binding sites at the endogenous C. elegans lin-3 locus, resulting in the wild-type expression level. Thus, regulatory evolution occurring in cis within a 58 bp lin-3 fragment, results in a strict requirement for the NHR binding site and a second E-box in C. elegans. This single-cell, single-molecule, quantitative and functional evo-devo study demonstrates that conserved expression levels can hide extensive change in cis-regulatory site requirements and highlights the evolution of new cis-regulatory elements required for cell-specific gene expression.

Diversification of mechanisms regulating gene expression of key developmental factors is a major force in the evolution of development. However, in the past, comparisons of gene expression across different species have often been qualitative (i.e. ‘expression is on versus off’ in a certain cell) without precise quantification. New experimental methods now allow us to quantitatively compare the expression of gene homologs across species, with single cell resolution. Moreover, the development of genome editing tools enables the dissection of regulatory DNA sequences that drive gene expression. We use here a well-established “textbook” example of animal organogenesis in the microscopic nematode, Caenorhabditis elegans, focusing on the expression of lin-3, coding for the main inducer of the vulva, in a single cell called the anchor cell. We find that the lin-3 expression level is remarkably conserved, with 20–25 messenger RNAs per anchor cell, in species that are molecularly as distant as fish and mammals. This conservation occurs despite substantial changes and compensation in the regulatory elements required for cell-specific gene expression.

Cell-Autonomous and Non-Cell-Autonomous Regulation of a Feeding State-Dependent Chemoreceptor Gene via MEF-2 and bHLH Transcription Factors

Food and feeding-state dependent changes in chemoreceptor gene expression may allow Caenorhabditis elegans to modify their chemosensory behavior, but the mechanisms essential for these expression changes remain poorly characterized. We had previously shown that expression of a feeding state-dependent chemoreceptor gene, srh-234, in the ADL sensory neuron of C. elegans is regulated via the MEF-2 transcription factor. Here, we show that MEF-2 acts together with basic helix-loop-helix (bHLH) transcription factors to regulate srh-234 expression as a function of feeding state. We identify a cis-regulatory MEF2 binding site that is necessary and sufficient for the starvation-induced down regulation of srh-234 expression, while an E-box site known to bind bHLH factors is required to drive srh-234 expression in ADL. We show that HLH-2 (E/Daughterless), HLH-3 and HLH-4 (Achaete-scute homologs) act in ADL neurons to regulate srh-234 expression. We further demonstrate that the expression levels of srh-234 in ADL neurons are regulated remotely by MXL-3 (Max-like 3 homolog) and HLH-30 (TFEB ortholog) acting in the intestine, which is dependent on insulin signaling functioning specifically in ADL neurons. We also show that this intestine-to-neuron feeding-state regulation of srh-234 involves a subset of insulin-like peptides. These results combined suggest that chemoreceptor gene expression is regulated by both cell-autonomous and non-cell-autonomous transcriptional mechanisms mediated by MEF2 and bHLH factors, which may allow animals to fine-tune their chemosensory responses in response to changes in their feeding state.

Plasticity in chemoreceptor gene expression may be a simple strategy by which an animal can modulate its chemosensory responses in changing external and internal state conditions. However, the transcriptional mechanisms required for these chemoreceptor gene expression changes are poorly understood. Here, we describe the identification of a transcriptional module(s) consisting of MEF-2 and basic helix-loop-helix (bHLH) transcription factors and their cognate binding sites in Caenorhabditis elegans that act together in ADL sensory neurons to properly regulate expression of a feeding-state dependent chemoreceptor gene. We also showed that chemoreceptor gene expression in ADL neurons are regulated remotely by bHLH factors acting in the intestine through an insulin-mediated signaling pathway, implying a sensory neuron-gut interaction for modulating chemoreceptor gene expression as a function of feeding state. This work describes transcriptional mechanisms mediated by MEF-2 and bHLH factors by which the expression of individual chemoreceptor genes in C. elegans are changed in response to changes in feeding state conditions.

Anti-aging treatments slow propagation of synucleinopathy by restoring lysosomal function

Aging is the major risk factor for neurodegenerative diseases that are also associated with impaired proteostasis, resulting in abnormal accumulation of protein aggregates. However, the role of aging in development and progression of disease remains elusive. Here, we used Caenorhabditis elegans models to show that aging-promoting genetic variations accelerated the rate of cell-to-cell transmission of SNCA/α-synuclein aggregates, hallmarks of Parkinson disease, and the progression of disease phenotypes, such as nerve degeneration, behavioral deficits, and reduced life span. Genetic and pharmacological anti-aging manipulations slowed the spread of aggregates and the associated phenotypes. Lysosomal degradation was significantly impaired in aging models, while anti-aging treatments reduced the impairment. Transgenic expression of hlh-30p::hlh-30, the master controller of lysosomal biogenesis, alleviated intercellular transmission of aggregates in the aging model. Our results demonstrate that the rate of aging closely correlates with the rate of aggregate propagation and that general anti-aging treatments can slow aggregate propagation and associated disease progression by restoring lysosomal function.

Dry and wet approaches for genome-wide functional annotation of conventional and unconventional transcriptional activators

Transcription factors (TFs) are master gene products that regulate gene expression in response to a variety of stimuli. They interact with DNA in a sequence-specific manner using a variety of DNA-binding domain (DBD) modules. This allows to properly position their second domain, called "effector domain", to directly or indirectly recruit positively or negatively acting co-regulators including chromatin modifiers, thus modulating preinitiation complex formation as well as transcription elongation. At variance with the DBDs, which are comprised of well-defined and easily recognizable DNA binding motifs, effector domains are usually much less conserved and thus considerably more difficult to predict. Also not so easy to identify are the DNA-binding sites of TFs, especially on a genome-wide basis and in the case of overlapping binding regions. Another emerging issue, with many potential regulatory implications, is that of so-called "moonlighting" transcription factors, i.e., proteins with an annotated function unrelated to transcription and lacking any recognizable DBD or effector domain, that play a role in gene regulation as their second job. Starting from bioinformatic and experimental high-throughput tools for an unbiased, genome-wide identification and functional characterization of TFs (especially transcriptional activators), we describe both established (and usually well affordable) as well as newly developed platforms for DNA-binding site identification. Selected combinations of these search tools, some of which rely on next-generation sequencing approaches, allow delineating the entire repertoire of TFs and unconventional regulators encoded by the any sequenced genome.

Quantification of Cooperativity in Heterodimer-DNA Binding Improves the Accuracy of Binding Specificity Models

Many transcription factors (TFs) have the ability to cooperate on DNA elements as heterodimers. Despite the significance of TF heterodimerization for gene regulation, a quantitative understanding of cooperativity between various TF dimer partners and its impact on heterodimer DNA binding specificity models is still lacking. Here, we used a novel integrative approach, combining microfluidics-steered measurements of dimer-DNA assembly with mechanistic modeling of the implicated protein-protein-DNA interactions to quantitatively interrogate the cooperative DNA binding behavior of the adipogenic peroxisome proliferator-activated receptor γ (PPARγ):retinoid X receptor α (RXRα) heterodimer. Using the high throughput MITOMI (mechanically induced trapping of molecular interactions) platform, we derived equilibrium DNA binding data for PPARγ, RXRα, as well as the PPARγ:RXRα heterodimer to more than 300 target DNA sites and variants thereof. We then quantified cooperativity underlying heterodimer-DNA binding and derived an integrative heterodimer DNA binding constant. Using this cooperativity-inclusive constant, we were able to build a heterodimer-DNA binding specificity model that has superior predictive power than the one based on a regular one-site equilibrium. Our data further revealed that individual nucleotide substitutions within the target site affect the extent of cooperativity in PPARγ:RXRα-DNA binding. Our study therefore emphasizes the importance of assessing cooperativity when generating DNA binding specificity models for heterodimers.

Deciphering the molecular mechanisms underlying the binding of the TWIST1/E12 complex to regulatory E-box sequences

The TWIST1 bHLH transcription factor controls embryonic development and cancer processes. Although molecular and genetic analyses have provided a wealth of data on the role of bHLH transcription factors, very little is known on the molecular mechanisms underlying their binding affinity to the E-box sequence of the promoter. Here, we used an in silico model of the TWIST1/E12 (TE) heterocomplex and performed molecular dynamics (MD) simulations of its binding to specific (TE-box) and modified E-box sequences. We focused on (i) active E-box and inactive E-box sequences, on (ii) modified active E-box sequences, as well as on (iii) two box sequences with modified adjacent bases the AT- and TA-boxes. Our in silico models were supported by functional in vitro binding assays. This exploration highlighted the predominant role of protein side-chain residues, close to the heart of the complex, at anchoring the dimer to DNA sequences, and unveiled a shift towards adjacent ((-1) and (-1*)) bases and conserved bases of modified E-box sequences. In conclusion, our study provides proof of the predictive value of these MD simulations, which may contribute to the characterization of specific inhibitors by docking approaches, and their use in pharmacological therapies by blocking the tumoral TWIST1/E12 function in cancers.

Identification of plant transcription factor target sequences

Regulation of gene expression depends on specific cis-regulatory sequences located in the gene promoter regions. These DNA sequences are recognized by transcription factors (TFs) in a sequence-specific manner, and their identification could help to elucidate the regulatory networks that underlie plant physiological responses to developmental programs or to environmental adaptation. Here we review recent advances in high throughput methodologies for the identification of plant TF binding sites. Several approaches offer a map of the TF binding locations in vivo and of the dynamics of the gene regulatory networks. As an alternative, high throughput in vitro methods provide comprehensive determination of the DNA sequences recognized by TFs. These advances are helping to decipher the regulatory lexicon and to elucidate transcriptional network hierarchies in plants in response to internal or external cues. This article is part of a Special Issue entitled: Plant Gene Regulatory Mechanisms and Networks, edited by Dr. Erich Grotewold and Dr. Nathan Springer.

Transcriptional and epigenetic regulation of autophagy in aging

Macroautophagy is a major intracellular degradation process recognized as playing a central role in cell survival and longevity. This multistep process is extensively regulated at several levels, including post-translationally through the action of conserved longevity factors such as the nutrient sensor TOR. More recently, transcriptional regulation of autophagy genes has emerged as an important mechanism for ensuring the somatic maintenance and homeostasis necessary for a long life span. Autophagy is increased in many long-lived model organisms and contributes significantly to their longevity. In turn, conserved transcription factors, particularly the helix-loop-helix transcription factor TFEB and the forkhead transcription factor FOXO, control the expression of many autophagy-related genes and are important for life-span extension. In this review, we discuss recent progress in understanding the contribution of these transcription factors to macroautophagy regulation in the context of aging. We also review current research on epigenetic changes, such as histone modification by the deacetylase SIRT1, that influence autophagy-related gene expression and additionally affect aging. Understanding the molecular regulation of macroautophagy in relation to aging may offer new avenues for the treatment of age-related diseases.

An inter-species protein-protein interaction network across vast evolutionary distance

In cellular systems, biophysical interactions between macromolecules underlie a complex web of functional interactions. How biophysical and functional networks are coordinated, whether all biophysical interactions correspond to functional interactions, and how such biophysical-versus-functional network coordination is shaped by evolutionary forces are all largely unanswered questions. Here, we investigate these questions using an "inter-interactome" approach. We systematically probed the yeast and human proteomes for interactions between proteins from these two species and functionally characterized the resulting inter-interactome network. After a billion years of evolutionary divergence, the yeast and human proteomes are still capable of forming a biophysical network with properties that resemble those of intra-species networks. Although substantially reduced relative to intra-species networks, the levels of functional overlap in the yeast-human inter-interactome network uncover significant remnants of co-functionality widely preserved in the two proteomes beyond human-yeast homologs. Our data support evolutionary selection against biophysical interactions between proteins with little or no co-functionality. Such non-functional interactions, however, represent a reservoir from which nascent functional interactions may arise.

Contrasting invertebrate immune defense behaviors caused by a single gene, the Caenorhabditis elegans neuropeptide receptor gene npr-1

Background

The invertebrate immune system comprises physiological mechanisms, physical barriers and also behavioral responses. It is generally related to the vertebrate innate immune system and widely believed to provide nonspecific defense against pathogens, whereby the response to different pathogen types is usually mediated by distinct signalling cascades. Recent work suggests that invertebrate immune defense can be more specific at least at the phenotypic level. The underlying genetic mechanisms are as yet poorly understood.

Results

We demonstrate in the model invertebrate Caenorhabditis elegans that a single gene, a homolog of the mammalian neuropeptide Y receptor gene, npr-1, mediates contrasting defense phenotypes towards two distinct pathogens, the Gram-positive Bacillus thuringiensis and the Gram-negative Pseudomonas aeruginosa. Our findings are based on combining quantitative trait loci (QTLs) analysis with functional genetic analysis and RNAseq-based transcriptomics. The QTL analysis focused on behavioral immune defense against B. thuringiensis, using recombinant inbred lines (RILs) and introgression lines (ILs). It revealed several defense QTLs, including one on chromosome X comprising the npr-1 gene. The wildtype N2 allele for the latter QTL was associated with reduced defense against B. thuringiensis and thus produced an opposite phenotype to that previously reported for the N2 npr-1 allele against P. aeruginosa. Analysis of npr-1 mutants confirmed these contrasting immune phenotypes for both avoidance behavior and nematode survival. Subsequent transcriptional profiling of C. elegans wildtype and npr-1 mutant suggested that npr-1 mediates defense against both pathogens through p38 MAPK signaling, insulin-like signaling, and C-type lectins. Importantly, increased defense towards P. aeruginosa seems to be additionally influenced through the induction of oxidative stress genes and activation of GATA transcription factors, while the repression of oxidative stress genes combined with activation of Ebox transcription factors appears to enhance susceptibility to B. thuringiensis.

Conclusions

Our findings highlight the role of a single gene, npr-1, in fine-tuning nematode immune defense, showing the ability of the invertebrate immune system to produce highly specialized and potentially opposing immune responses via single regulatory genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2603-8) contains supplementary material, which is available to authorized users.

Mondo complexes regulate TFEB via TOR inhibition to promote longevity in response to gonadal signals

Germline removal provokes longevity in several species and shifts resources towards survival and repair. Several Caenorhabditis elegans transcription factors regulate longevity arising from germline removal; yet, how they work together is unknown. Here we identify a Myc-like HLH transcription factor network comprised of Mondo/Max-like complex (MML-1/MXL-2) to be required for longevity induced by germline removal, as well as by reduced TOR, insulin/IGF signalling and mitochondrial function. Germline removal increases MML-1 nuclear accumulation and activity. Surprisingly, MML-1 regulates nuclear localization and activity of HLH-30/TFEB, a convergent regulator of autophagy, lysosome biogenesis and longevity, by downregulating TOR signalling via LARS-1/leucyl-transfer RNA synthase. HLH-30 also upregulates MML-1 upon germline removal. Mammalian MondoA/B and TFEB show similar mutual regulation. MML-1/MXL-2 and HLH-30 transcriptomes show both shared and preferential outputs including MDL-1/MAD-like HLH factor required for longevity. These studies reveal how an extensive interdependent HLH transcription factor network distributes responsibility and mutually enforces states geared towards reproduction or survival.

Overlapping cell population expression profiling and regulatory inference in C. elegans

Background

Understanding gene expression across the diverse metazoan cell types during development is critical to understanding their function and regulation. However, most cell types have not been assayed for expression genome-wide.

Results

We applied a novel approach we term “Profiling of Overlapping Populations of cells (POP-Seq)” to assay differential expression across all embryonic cells in the nematode Caenorhabditis elegans. In this approach, we use RNA-seq to define the transcriptome of diverse partially overlapping FACS-sorted cell populations. This identified thousands of transcripts differentially expressed across embryonic cells. Hierarchical clustering analysis identified over 100 sets of coexpressed genes corresponding to distinct patterns of cell type specific expression. We identified thousands of candidate regulators of these clusters based on enrichment of transcription factor motifs and experimentally determined binding sites.

Conclusions

Our analysis provides new insight into embryonic gene regulation, and provides a resource for improving our knowledge of tissue-specific expression and its regulation throughout C. elegans development.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2482-z) contains supplementary material, which is available to authorized users.

Control of lysosomal biogenesis and Notch-dependent tissue patterning by components of the TFEB-V-ATPase axis in Drosophila melanogaster

In vertebrates, TFEB (transcription factor EB) and MITF (microphthalmia-associated transcription factor) family of basic Helix-Loop-Helix (bHLH) transcription factors regulates both lysosomal function and organ development. However, it is not clear whether these 2 processes are interconnected. Here, we show that Mitf, the single TFEB and MITF ortholog in Drosophila, controls expression of vacuolar-type H(+)-ATPase pump (V-ATPase) subunits. Remarkably, we also find that expression of Vha16-1 and Vha13, encoding 2 key components of V-ATPase, is patterned in the wing imaginal disc. In particular, Vha16-1 expression follows differentiation of proneural regions of the disc. These regions, which will form sensory organs in the adult, appear to possess a distinctive endolysosomal compartment and Notch (N) localization. Modulation of Mitf activity in the disc in vivo alters endolysosomal function and disrupts proneural patterning. Similar to our findings in Drosophila, in human breast epithelial cells we observe that impairment of the Vha16-1 human ortholog ATP6V0C changes the size and function of the endolysosomal compartment and that depletion of TFEB reduces ligand-independent N signaling activity. Our data suggest that lysosomal-associated functions regulated by the TFEB-V-ATPase axis might play a conserved role in shaping cell fate.

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.

A Glimpse to Background and Characteristics of Major Molecular Biological Networks

Recently, biology has become a data intensive science because of huge data sets produced by high throughput molecular biological experiments in diverse areas including the fields of genomics, transcriptomics, proteomics, and metabolomics. These huge datasets have paved the way for system-level analysis of the processes and subprocesses of the cell. For system-level understanding, initially the elements of a system are connected based on their mutual relations and a network is formed. Among omics researchers, construction and analysis of biological networks have become highly popular. In this review, we briefly discuss both the biological background and topological properties of major types of omics networks to facilitate a comprehensive understanding and to conceptualize the foundation of network biology.

Dimerization-driven degradation of C. elegans and human E proteins

In this study, Sallee et al. demonstrate that E-protein dimer formation can promote C. elegans and human bHLH protein instability. By investigating HLH-2, the sole C. elegans E protein, the authors show that HLH-2 functions as a homodimer for sequential roles in AC specification and differentiation and that the functional dimer is targeted for degradation in VUs, the “opposite” fate. The findings indicate that dimerization-driven regulation of bHLH protein stability may be a conserved mechanism for differential regulation in specific cell contexts.

E proteins are conserved regulators of growth and development. We show that the Caenorhabditis elegans E-protein helix–loop–helix-2 (HLH-2) functions as a homodimer in directing development and function of the anchor cell (AC) of the gonad, the critical organizer of uterine and vulval development. Our structure–function analysis of HLH-2 indicates that dimerization drives its degradation in other uterine cells (ventral uterine precursor cells [VUs]) that initially have potential to be the AC. We also provide evidence that this mode of dimerization-driven down-regulation can target other basic HLH (bHLH) dimers as well. Remarkably, human E proteins can functionally substitute for C. elegans HLH-2 in regulating AC development and also display dimerization-dependent degradation in VUs. Our results suggest that dimerization-driven regulation of bHLH protein stability may be a conserved mechanism for differential regulation in specific cell contexts.

Nonconsensus Protein Binding to Repetitive DNA Sequence Elements Significantly Affects Eukaryotic Genomes

Recent genome-wide experiments in different eukaryotic genomes provide an unprecedented view of transcription factor (TF) binding locations and of nucleosome occupancy. These experiments revealed that a large fraction of TF binding events occur in regions where only a small number of specific TF binding sites (TFBSs) have been detected. Furthermore, in vitro protein-DNA binding measurements performed for hundreds of TFs indicate that TFs are bound with wide range of affinities to different DNA sequences that lack known consensus motifs. These observations have thus challenged the classical picture of specific protein-DNA binding and strongly suggest the existence of additional recognition mechanisms that affect protein-DNA binding preferences. We have previously demonstrated that repetitive DNA sequence elements characterized by certain symmetries statistically affect protein-DNA binding preferences. We call this binding mechanism nonconsensus protein-DNA binding in order to emphasize the point that specific consensus TFBSs do not contribute to this effect. In this paper, using the simple statistical mechanics model developed previously, we calculate the nonconsensus protein-DNA binding free energy for the entire C. elegans and D. melanogaster genomes. Using the available chromatin immunoprecipitation followed by sequencing (ChIP-seq) results on TF-DNA binding preferences for ~100 TFs, we show that DNA sequences characterized by low predicted free energy of nonconsensus binding have statistically higher experimental TF occupancy and lower nucleosome occupancy than sequences characterized by high free energy of nonconsensus binding. This is in agreement with our previous analysis performed for the yeast genome. We suggest therefore that nonconsensus protein-DNA binding assists the formation of nucleosome-free regions, as TFs outcompete nucleosomes at genomic locations with enhanced nonconsensus binding. In addition, here we perform a new, large-scale analysis using in vitro TF-DNA preferences obtained from the universal protein binding microarrays (PBM) for ~90 eukaryotic TFs belonging to 22 different DNA-binding domain types. As a result of this new analysis, we conclude that nonconsensus protein-DNA binding is a widespread phenomenon that significantly affects protein-DNA binding preferences and need not require the presence of consensus (specific) TFBSs in order to achieve genome-wide TF-DNA binding specificity.

Interactions between proteins and DNA trigger many important biological processes. Therefore, to fully understand how the information encoded on the DNA transcribes into RNA, which in turn translates into proteins in the cell, we need to unravel the molecular design principles of protein-DNA interactions. It is known that many interactions occur when a protein is attracted to a specific short segment on the DNA called a specific protein-DNA binding motif. Strikingly, recent experiments revealed that many regulatory proteins reproducibly bind to different regions on the DNA lacking such specific motifs. This suggests that fundamental molecular mechanisms responsible for protein-DNA recognition specificity are not fully understood. Here, using high-throughput protein-DNA binding data obtained by two entirely different methods for ~100 TFs in each case, we show that DNA regions possessing certain repetitive sequence elements exert the statistical attractive potential on DNA-binding proteins, and as a result, such DNA regions are enriched in bound proteins. This is in agreement with our previous analysis performed for the yeast genome. We use the term nonconsensus protein-DNA binding in order to describe protein-DNA interactions that occur in the absence of specific protein-DNA binding motifs. Here we demonstrate that the identified nonconsensus effect is highly significant for a variety of organismal genomes and it affects protein-DNA binding preferences and nucleosome occupancy at the genome-wide level.

Mapping and analysis of Caenorhabditis elegans transcription factor sequence specificities

Caenorhabditis elegans is a powerful model for studying gene regulation, as it has a compact genome and a wealth of genomic tools. However, identification of regulatory elements has been limited, as DNA-binding motifs are known for only 71 of the estimated 763 sequence-specific transcription factors (TFs). To address this problem, we performed protein binding microarray experiments on representatives of canonical TF families in C. elegans, obtaining motifs for 129 TFs. Additionally, we predict motifs for many TFs that have DNA-binding domains similar to those already characterized, increasing coverage of binding specificities to 292 C. elegans TFs (∼40%). These data highlight the diversification of binding motifs for the nuclear hormone receptor and C2H2 zinc finger families and reveal unexpected diversity of motifs for T-box and DM families. Motif enrichment in promoters of functionally related genes is consistent with known biology and also identifies putative regulatory roles for unstudied TFs.

DOI:http://dx.doi.org/10.7554/eLife.06967.001

Many scientists use ‘model’ species—such as the fruit fly or a nematode worm called Caenorhabditis elegans—in their research because these organisms have useful features that make it easier to carry out many experiments. For example, C. elegans has a smaller genome compared to many other animals, which is useful for studying the roles of individual genes or stretches of DNA.

Transcription factors are a type of protein that can bind to specific stretches of DNA and help to switch certain genes on or off. These ‘motifs’ may be close to the gene or further away in the genome, and therefore, must stand out amongst the rest of the DNA, like lights on a landing strip. However, the motifs for only 10% of the estimated 763 transcription factors in C. elegans have been identified so far.

In this study, Narasimhan, Lambert, Yang et al. used a technique called a ‘protein binding microarray’ to identify the motifs for many more of the C. elegans transcription factors. These findings were then used to predict motifs for other transcription factors. Together, these methods increased the proportion of C. elegans transcription factors with known DNA-binding motifs from 10% to around 40%.

Now that more DNA motifs have been identified, it is possible to look for similarities and differences between them. For example, Narasimhan, Lambert, Yang et al. found that transcription factors with similar sequences can bind to very varied motifs. On the other hand, some transcription factors that are very different are able to recognize very similar motifs. The experiments also indicate that motifs found very close to genes—in sequences known as ‘promoters’—may be able to interact with many proteins to influence the activity of genes.

Narasimhan, Lambert, Yang et al.'s findings increase the number of C. elegans transcription factors with a motif, bringing the knowledge of these proteins more in line with the better-studied transcription factors of humans and fruit flies. The next challenge is to identify DNA motifs for the remaining 60% of transcription factors.

DOI:http://dx.doi.org/10.7554/eLife.06967.002

Human gene-centered transcription factor networks for enhancers and disease variants

Gene regulatory networks (GRNs) comprising interactions between transcription factors (TFs) and regulatory loci control development and physiology. Numerous disease-associated mutations have been identified, the vast majority residing in non-coding regions of the genome. As current GRN mapping methods test one TF at a time and require the use of cells harboring the mutation(s) of interest, they are not suitable to identify TFs that bind to wild-type and mutant loci. Here, we use gene-centered yeast one-hybrid (eY1H) assays to interrogate binding of 1,086 human TFs to 246 enhancers, as well as to 109 non-coding disease mutations. We detect both loss and gain of TF interactions with mutant loci that are concordant with target gene expression changes. This work establishes eY1H assays as a powerful addition to the toolkit of mapping human GRNs and for the high-throughput characterization of genomic variants that are rapidly being identified by genome-wide association studies.

Plasticity of chemoreceptor gene expression: Sensory and circuit inputs modulate state-dependent chemoreceptors

Animals dramatically modify their chemosensory behaviors when starved, which could allow them to alter and optimize their food-search strategies. Dynamic changes in the gene expression of chemoreceptors may be a general mechanism underlying food and state-dependent changes in chemosensory behaviors. In our recent study,(1) we identified chemoreceptors in the ADL sensory neuron type of C. elegans that are modulated by feeding state and food availability. Here, we highllight our recent findings by which sensory inputs into ADL, neuronal outputs from ADL, and circuit inputs from the RMG interneuron, which is electrically connected to ADL, are required to regulate an ADL-expressed chemoreceptor. This sensory and circuit-mediated regulation of chemoreceptor gene expression is dependent on cell-autonomous pathways acting in ADL, e.g. KIN-29, DAF-2, OCR-2 and calcium signaling, and circuit inputs from RMG mediated by NPR-1. Based on these findings, we propose an intriguing but speculative feedback modulatory circuit mechanism by which sensory perception of food and internal state signals may be coupled to regulate ADL-expressed chemoreceptors, which may allow animals to precisely regulate and fine-tune their chemosensory neuron responses as a function of feeding state.

The Bicoid class homeodomain factors ceh-36/OTX and unc-30/PITX cooperate in C. elegans embryonic progenitor cells to regulate robust development

While many transcriptional regulators of pluripotent and terminally differentiated states have been identified, regulation of intermediate progenitor states is less well understood. Previous high throughput cellular resolution expression studies identified dozens of transcription factors with lineage-specific expression patterns in C. elegans embryos that could regulate progenitor identity. In this study we identified a broad embryonic role for the C. elegans OTX transcription factor ceh-36, which was previously shown to be required for the terminal specification of four neurons. ceh-36 is expressed in progenitors of over 30% of embryonic cells, yet is not required for embryonic viability. Quantitative phenotyping by computational analysis of time-lapse movies of ceh-36 mutant embryos identified cell cycle or cell migration defects in over 100 of these cells, but most defects were low-penetrance, suggesting redundancy. Expression of ceh-36 partially overlaps with that of the PITX transcription factor unc-30. unc-30 single mutants are viable but loss of both ceh-36 and unc-30 causes 100% lethality, and double mutants have significantly higher frequencies of cellular developmental defects in the cells where their expression normally overlaps. These factors are also required for robust expression of the downstream developmental regulator mls-2/HMX. This work provides the first example of genetic redundancy between the related yet evolutionarily distant OTX and PITX families of bicoid class homeodomain factors and demonstrates the power of quantitative developmental phenotyping in C. elegans to identify developmental regulators acting in progenitor cells.

How do regulatory networks evolve and expand throughout evolution?

Throughout evolution, regulatory networks need to expand and adapt to accommodate novel genes and gene functions. However, the molecular details explaining how gene networks evolve remain largely unknown. Recent studies demonstrate that changes in transcription factors contribute to the evolution of regulatory networks. In particular, duplication of transcription factors followed by specific mutations in their DNA-binding or interaction domains propels the divergence and emergence of new networks. The innate promiscuity and modularity of regulatory networks contributes to their evolvability: duplicated promiscuous regulators and their target promoters can acquire mutations that lead to gradual increases in specificity, allowing neofunctionalization or subfunctionalization.

Evolution of transcription factor function as a mechanism for changing metazoan developmental gene regulatory networks

The form that an animal takes during development is directed by gene regulatory networks (GRNs). Developmental GRNs interpret maternally deposited molecules and externally supplied signals to direct cell-fate decisions, which ultimately leads to the arrangements of organs and tissues in the organism. Genetically encoded modifications to these networks have generated the wide range of metazoan diversity that exists today. Most studies of GRN evolution focus on changes to cis-regulatory DNA, and it was historically theorized that changes to the transcription factors that bind to these cis-regulatory modules (CRMs) contribute to this process only rarely. A growing body of evidence suggests that changes to the coding regions of transcription factors play a much larger role in the evolution of developmental gene regulatory networks than originally imagined. Just as cis-regulatory changes make use of modular binding site composition and tissue-specific modules to avoid pleiotropy, transcription factor coding regions also predominantly evolve in ways that limit the context of functional effects. Here, we review the recent works that have led to this unexpected change in the field of Evolution and Development (Evo-Devo) and consider the implications these studies have had on our understanding of the evolution of developmental processes.