Identification of a novel cis-regulatory element involved in the heat shock response in Caenorhabditis elegans using microarray gene expression and computational methods

EnhancerTracker: Comparing cell-type-specific enhancer activity of DNA sequence triplets via an ensemble of deep convolutional neural networks

Motivation

Transcriptional enhancers - unlike promoters - are unrestrained by distance or strand orientation with respect to their target genes, making their computational identification a challenge. Further, there are insufficient numbers of confirmed enhancers for many cell types, preventing robust training of machine-learning-based models for enhancer prediction for such cell types.

Results

We present EnhancerTracker , a novel tool that leverages an ensemble of deep separable convolutional neural networks to identify cell-type-specific enhancers with the need of only two confirmed enhancers. EnhancerTracker is trained, validated, and tested on 52,789 putative enhancers obtained from the FANTOM5 Project and control sequences derived from the human genome. Unlike available tools, which accept one sequence at a time, the input to our tool is three sequences; the first two are enhancers active in the same cell type. EnhancerTracker outputs 1 if the third sequence is an enhancer active in the same cell type(s) where the first two enhancers are active. It outputs 0 otherwise. On a held-out set (15%), EnhancerTracker achieved an accuracy of 64%, a specificity of 93%, a recall of 35%, a precision of 84%, and an F1 score of 49%.

Availability and implementation

https://github.com/BioinformaticsToolsmith/EnhancerTracker.

Contact

hani.girgis@tamuk.edu.

Identification of a pluripotency-inducing small compound, PLU, that induces callus formation via Heat Shock Protein 90-mediated activation of auxin signaling

Plants retain the ability to generate a pluripotent tissue called callus by dedifferentiating somatic cells. A pluripotent callus can also be artificially induced by culturing explants with hormone mixtures of auxin and cytokinin, and an entire body can then be regenerated from the callus. Here we identified a pluripotency-inducing small compound, PLU, that induces the formation of callus with tissue regeneration potency without the external application of either auxin or cytokinin. The PLU-induced callus expressed several marker genes related to pluripotency acquisition via lateral root initiation processes. PLU-induced callus formation required activation of the auxin signaling pathway though the amount of active auxin was reduced by PLU treatment. RNA-seq analysis and subsequent experiments revealed that Heat Shock Protein 90 (HSP90) mediates a significant part of the PLU-initiated early events. We also showed that HSP90-dependent induction of TRANSPORT INHIBITOR RESPONSE 1, an auxin receptor gene, is required for the callus formation by PLU. Collectively, this study provides a new tool for manipulating and investigating the induction of plant pluripotency from a different angle from the conventional method with the external application of hormone mixtures.

The emergent role of mitochondrial surveillance in cellular health

Mitochondrial dysfunction is one of the primary causatives for many pathologies, including neurodegenerative diseases, cancer, metabolic disorders, and aging. Decline in mitochondrial functions leads to the loss of proteostasis, accumulation of ROS, and mitochondrial DNA damage, which further exacerbates mitochondrial deterioration in a vicious cycle. Surveillance mechanisms, in which mitochondrial functions are closely monitored for any sign of perturbations, exist to anticipate possible havoc within these multifunctional organelles with primitive origin. Various indicators of unhealthy mitochondria, including halted protein import, dissipated membrane potential, and increased loads of oxidative damage, are on the top of the lists for close monitoring. Recent research also indicates a possibility of reductive stress being monitored as part of a mitochondrial surveillance program. Upon detection of mitochondrial stress, multiple mitochondrial stress-responsive pathways are activated to promote the transcription of numerous nuclear genes to ameliorate mitochondrial damage and restore compromised cellular functions. Co-expression occurs through functionalization of transcription factors, allowing their binding to promoter elements to initiate transcription of target genes. This review provides a comprehensive summary of the intricacy of mitochondrial surveillance programs and highlights their roles in our cellular life. Ultimately, a better understanding of these surveillance mechanisms is expected to improve healthspan.

Digital Microfluidics Supported Microproteomics for Quantitative Proteome Analysis of Single Caenorhabditis elegans Nematodes

Miniaturization of sample preparation, including omissible manual sample handling steps, is key for reproducible nanoproteomics, as material is often restricted to only hundreds of cells or single model organisms. Here, we demonstrate a highly sensitive digital microfluidics (DMF)-based sample preparation workflow making use of single-pot solid-phase enhanced sample preparation (SP3) in combination with high-field asymmetric-waveform ion mobility spectrometry (FAIMS), and fast and sensitive ion trap detection on an Orbitrap tribrid MS system. Compared to a manual in-tube SP3-supported sample preparation, the numbers of identified peptides and proteins were markedly increased, while lower standard deviations between replicates were observed. We repeatedly identified up to 5000 proteins from single nematodes. Moreover, label-free quantification of protein changes in single Caenorhabditis elegans treated with a heat stimulus yielded 45 differentially abundant proteins when compared to the untreated control, highlighting the potential of this technology for low-input proteomics studies. LC-MS data have been deposited to the ProteomeXchange Consortium with the data set identifier PXD033143.

Binding of the HSF-1 DNA-binding domain to multimeric C. elegans consensus HSEs is guided by cooperative interactions

The protein HSF-1 is the controlling transcription factor of the heat-shock response (HSR). Its binding to the heat-shock elements (HSEs) induces the strong upregulation of conserved heat-shock proteins, including Hsp70s, Hsp40s and small HSPs. Next to these commonly known HSPs, more than 4000 other HSEs are found in the promoter regions of C. elegans genes. In microarray experiments, few of the HSE-containing genes are specifically upregulated during the heat-shock response. Most of the 4000 HSE-containing genes instead are unaffected by elevated temperatures and coexpress with genes unrelated to the HSR. This is also the case for several genes related to the HSP chaperone system, like dnj-12, dnj-13, and hsp-1. Interestingly, several promoters of the dedicated HSR-genes, like F44E5.4p, hsp-16.48p or hsp-16.2p, contain extended HSEs in their promoter region, composed of four or five HSE-elements instead of the common trimeric HSEs. We here aim at understanding how HSF-1 interacts with the different promoter regions. To this end we purify the nematode HSF-1 DBD and investigate the interaction with DNA sequences containing these regions. EMSA assays suggest that the HSF-1 DBD interacts with most of these HSE-containing dsDNAs, but with different characteristics. We employ sedimentation analytical ultracentrifugation (SV-AUC) to determine stoichiometry, affinity, and cooperativity of HSF-1 DBD binding to these HSEs. Interestingly, most HSEs show cooperative binding of the HSF-1 DBD with up to five DBDs being bound. In most cases binding to the HSEs of inducible promoters is stronger, even though the consensus scores are not always higher. The observed high affinity of HSF-1 DBD to the non-inducible HSEs of dnj-12, suggests that constitutive expression may be supported from some promoter regions, a fact that is evident for this transcription factor, that is essential also under non-stress conditions.

HSF-1 mediated combined ginsenosides ameliorating Alzheimer's disease like symptoms in Caernorhabditis elegans

There are few effective medications to treat Alzheimer's disease (AD). It has been suggested that several ginsenosides possess mild or moderate anti-AD activity. In our present work, a preferred combined ginsenosides was shown to have a more significant benefit effect on AD-like symptoms of worm paralysis and hypersensitivity to exogenous 5-HT in C. elegans. The combined ginsenosides can suppress Aβ deposits and Aβ oligomers, alleviating the toxicity induced by Aβ overexpression more effectively than used alone. Its anti-AD effect was partially abolished by hsf-1 RNAi knocked down or hsf-1 inactivation by point mutation, but not by daf-16 or skn-1 RNAi knocked down. Furthermore, it markedly activated hsp-16.2 gene expression downstream of HSF-1. Our results demonstrated that HSF-1 signaling pathway exerts an important role in mediating the therapeutic effect of combined ginsenosides on AD worms. These results provided powerful evidences and theoretical foundation for reshaping medicinal products of ginsenosides and ginseng on prevention of neurodegenerative diseases.

Quantitative Proteomic Analysis Revealed Broad Roles of N6-Methyladenosine in Heat Shock Response

As optimum temperature is essential for all living organisms, heat shock represents a challenging problem for their survival. Therefore, cellular response to heat shock is among the most extensively investigated stress response pathways; however, how the human proteome responds to heat shock has not been comprehensively investigated. In this study, we employed stable isotope labeling by amino acids in cell culture (SILAC), together with liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis, to fulfill an in-depth analysis of the alterations in the human proteome in M14 human melanoma cells in response to heat shock stress. We found that, after heat shock, 284 and 278 out of the 4319 quantified proteins were with substantially diminished and elevated expressions, respectively. We also examined the alterations in human kinome after heat shock by using our recently developed targeted proteomic method relying on parallel-reaction monitoring. Our results showed that the expression levels of 11 and 22 kinase proteins were increased and decreased, respectively, by at least 1.5-fold upon heat shock. By interrogating publicly available RNA-seq and m⁶A sequencing data, we observed that the elevated expression of more than 30 proteins, including CHEK1 and CCND3 kinases, could occur via an m⁶A-mediated mechanism. Furthermore, our results from single-base elongation and ligation-based quantitative polymerase chain reaction (qPCR) amplification (SELECT) and luciferase reporter assays revealed that heat shock gave rise to elevated m⁶A levels at A280 and A286 sites in the 5'-untranslated region of HSPH1 mRNA, thereby leading to increased translation of HSPH1 protein. Together, our discovery and targeted proteomic methods revealed the reprogramming of human proteome and kinome upon heat shock stress and provided insights into cellular responses toward heat shock stress.

Chronic temperature stress inhibits reproduction and disrupts endocytosis via chaperone titration in Caenorhabditis elegans

Background

Temperature influences biology at all levels, from altering rates of biochemical reactions to determining sustainability of entire ecosystems. Although extended exposure to elevated temperatures influences organismal phenotypes important for human health, agriculture, and ecology, the molecular mechanisms that drive these responses remain largely unexplored. Prolonged, mild temperature stress (48 h at 28 °C) has been shown to inhibit reproduction in Caenorhabditis elegans without significantly impacting motility or viability.

Results

Analysis of molecular responses to chronic stress using RNA-seq uncovers dramatic effects on the transcriptome that are fundamentally distinct from the well-characterized, acute heat shock response (HSR). While a large portion of the genome is differentially expressed ≥ 4-fold after 48 h at 28 °C, the only major class of oogenesis-associated genes affected is the vitellogenin gene family that encodes for yolk proteins (YPs). Whereas YP mRNAs decrease, the proteins accumulate and mislocalize in the pseudocoelomic space as early as 6 h, well before reproduction declines. A trafficking defect in a second, unrelated fluorescent reporter and a decrease in pre-synaptic neuronal signaling indicate that the YP mislocalization is caused by a generalized defect in endocytosis. Molecular chaperones are involved in both endocytosis and refolding damaged proteins. Decreasing levels of the major HSP70 chaperone, HSP-1, causes similar YP trafficking defects in the absence of stress. Conversely, increasing chaperone levels through overexpression of the transcription factor HSF-1 rescues YP trafficking and restores neuronal signaling.

Conclusions

These data implicate chaperone titration during chronic stress as a molecular mechanism contributing to endocytic defects that influence multiple aspects of organismal physiology. Notably, HSF-1 overexpression improves recovery of viable offspring after exposure to stress. These findings provide important molecular insights into understanding organismal responses to temperature stress as well as phenotypes associated with chronic protein misfolding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01008-1.

Transcriptional Enhancers in Drosophila

Key discoveries in Drosophila have shaped our understanding of cellular "enhancers." With a special focus on the fly, this chapter surveys properties of these adaptable cis-regulatory elements, whose actions are critical for the complex spatial/temporal transcriptional regulation of gene expression in metazoa. The powerful combination of genetics, molecular biology, and genomics available in Drosophila has provided an arena in which the developmental role of enhancers can be explored. Enhancers are characterized by diverse low- or high-throughput assays, which are challenging to interpret, as not all of these methods of identifying enhancers produce concordant results. As a model metazoan, the fly offers important advantages to comprehensive analysis of the central functions that enhancers play in gene expression, and their critical role in mediating the production of phenotypes from genotype and environmental inputs. A major challenge moving forward will be obtaining a quantitative understanding of how these cis-regulatory elements operate in development and disease.

The evolutionarily conserved ESRE stress response network is activated by ROS and mitochondrial damage

Background

Mitochondrial dysfunction causes or contributes to a wide variety of pathologies, including neurodegenerative diseases, cancer, metabolic diseases, and aging. Cells actively surveil a number of mitochondrial readouts to ensure that cellular homeostasis is maintained.

Results

In this article, we characterize the role of the ethanol and stress response element (ESRE) pathway in mitochondrial surveillance and show that it is robustly activated when the concentration of reactive oxygen species (ROS) in the cell increases. While experiments were mostly performed in Caenorhabditis elegans, we observed similar gene activation profile in human cell lines. The linear relationship between ROS and ESRE activation differentiates ESRE from known mitochondrial surveillance pathways, such as the mitochondrial unfolded protein response (UPR^mt), which monitor mitochondrial protein import. The ability of the ESRE network to be activated by increased ROS allows the cell to respond to oxidative and reductive stresses. The ESRE network works in tandem with other mitochondrial surveillance mechanisms as well, in a fashion that suggests a partially redundant hierarchy. For example, mutation of the UPR^mt pathway results in earlier and more robust activation of the ESRE pathway. Interestingly, full expression of ATFS-1, a key transcription factor for the UPR^mt, requires the presence of an ESRE motif in its promoter region.

Conclusion

The ESRE pathway responds to mitochondrial damage by monitoring ROS levels. This response is conserved in humans. The ESRE pathway is activated earlier when other mitochondrial surveillance pathways are unavailable during mitochondrial crises, potentially to mitigate stress and restore health. However, the exact mechanisms of pathway activation and crosstalk remain to be elucidated. Ultimately, a better understanding of this network, and its role in the constellation of mitochondrial and cellular stress networks, will improve healthspan.

Global Transcriptome Changes That Accompany Alterations in Serotonin Levels in Caenorhabditis elegans

Serotonin (5-hydroxytryptamine, 5-HT), is a phylogenetically ancient molecule best characterized as a neurotransmitter that modulates multiple aspects of mood and social cognition. The roles that 5-HT plays in normal and abnormal behavior are not fully understood but have been posited to be due to its common function as a 'defense signal'. However, 5-HT levels also systemically impact cell physiology, modulating cell division, migration, apoptosis, mitochondrial biogenesis, cellular metabolism and differentiation. Whether these diverse cellular effects of 5-HT also share a common basis is unclear. C. elegans provides an ideal system to interrogate the systemic effects of 5-HT, since lacking a blood-brain barrier, 5-HT synthesized and released by neurons permeates the organism to modulate neuronal as well as non-neuronal cells throughout the body. Here we used RNA-Seq to characterize the systemic changes in gene expression that occur in C. elegans upon altering 5-HT levels, and compared the transcriptomes to published datasets. We find that an acute increase in 5-HT is accompanied by a global decrease in gene expression levels, upregulation of genes involved in stress pathways, changes that significantly correlate with the published transcriptomes of animals that have activated defense and immune responses, and an increase in levels of phosphorylated eukaryotic initiation factor, eIF2α. In 5-HT deficient animals lacking tryptophan hydroxylase (tph-1(mg280)II) there is a net increase in gene expression, with an overrepresentation of genes related to development and chromatin. Surprisingly, the transcriptomes of animals with acute increases in 5-HT levels, and 5-HT deficiency do not overlap with transcriptomes of mutants with whom they share striking physiological resemblance. These studies are the first to catalog systemic transcriptome changes that occur upon alterations in 5-HT levels. They further show that in C. elegans changes in gene expression upon altering 5-HT levels, and changes in physiology, are not directly correlated.

A portable, low-cost device for precise control of specimen temperature under stereomicroscopes

To facilitate precise and convenient control of biological sample temperature, we developed a low-cost device that can be used independently or with any stereomicroscope. The purpose of the device is to control the thermal environment during experimental intervals in which a specimen must be manipulated outside of an incubator, e.g. for dissection or slide-mounting in preparation for imaging. Sample temperatures can be both cooled to below and heated to above room temperatures, and stably maintained at a precision of +/- 0.1˚C. To demonstrate the utility of this device, we report improved characterization of the penetrance of a short-acting temperature-sensitive allele in C. elegans embryos, and identification of the upper temperature threshold for embryonic viability for six Caenorhabditis species. By controlling the temperature environment even as a specimen is manipulated, this device offers consistency and flexibility, reduces environmental noise, and enables precision timing in experiments requiring temperature shifts.

Multivariable regulation of gene expression plasticity in metazoans

An important capacity of genes is the rapid change of expression levels to cope with the environment, known as expression responsiveness or plasticity. Elucidating the genomic mechanisms determining expression plasticity is critical for understanding the molecular basis of phenotypic plasticity, fitness and adaptation. In this study, we systematically quantified gene expression plasticity in four metazoan species by integrating changes of expression levels under a large number of genetic and environmental conditions. From this, we demonstrated that expression plasticity measures a distinct feature of gene expression that is orthogonal to other well-studied features, including gene expression level and tissue specificity/broadness. Expression plasticity is conserved across species with important physiological implications. The magnitude of expression plasticity is highly correlated with gene function and genes with high plasticity are implicated in disease susceptibility. Genome-wide analysis identified many conserved promoter cis-elements, trans-acting factors (such as CTCF), and gene body histone modifications (H3K36me3, H3K79me2 and H4K20me1) that are significantly associated with expression plasticity. Analysis of expression changes in perturbation experiments further validated a causal role of specific transcription factors and histone modifications. Collectively, this work reveals the general properties, physiological implications and multivariable regulation of gene expression plasticity in metazoans, extending the mechanistic understanding of gene regulation.

A secreted microRNA disrupts autophagy in distinct tissues of Caenorhabditis elegans upon ageing

Macroautophagy, a key player in protein quality control, is proposed to be systematically impaired in distinct tissues and causes coordinated disruption of protein homeostasis and ageing throughout the body. Although tissue-specific changes in autophagy and ageing have been extensively explored, the mechanism underlying the inter-tissue regulation of autophagy with ageing is poorly understood. Here, we show that a secreted microRNA, mir-83/miR-29, controls the age-related decrease in macroautophagy across tissues in Caenorhabditis elegans. Upregulated in the intestine by hsf-1/HSF1 with age, mir-83 is transported across tissues potentially via extracellular vesicles and disrupts macroautophagy by suppressing CUP-5/MCOLN, a vital autophagy regulator, autonomously in the intestine as well as non-autonomously in body wall muscle. Mutating mir-83 thereby enhances macroautophagy in different tissues, promoting protein homeostasis and longevity. These findings thus identify a microRNA-based mechanism to coordinate the decreasing macroautophagy in various tissues with age.

Decreased autophagy is a hallmark of ageing, but its inter-tissue regulation is poorly understood. Here, Zhou et al. identify mir-83 in C. elegans, which is transported across tissues and suppresses autophagy, contributing to age-related decline.

Ilex paraguariensis extract provides increased resistance against oxidative stress and protection against Amyloid beta-induced toxicity compared to caffeine in Caenorhabditis elegans

Ilex paraguariensis is a plant from South America, used to prepare a tea-like beverage rich in caffeine and polyphenols with antioxidant proprieties. Caffeine consumption is associated with a lower risk of age-associated neuropathologies, besides several extracts that have antioxidant proprieties are known to be neuroprotective, and oxidative stress strongly correlates with Aβ-toxicity. This study aims to investigate the neuroprotective effects of the Ilex paraguariensis hydroalcoholic extract (IPHE) and to evaluate if caffeine agent present in IPHE exerts neuroprotective effects in an amyloid beta-peptide (Aβ)-induced toxicity in Caenorhabditis elegans. The wild-type and CL2006 worms were treated with IPHE (2 and 4 mg/mL) or caffeine (200 and 400 μM) since larval stage 1 (L1) until they achieved the required age for each assay. IPHE and caffeine increased the lifespan and appeared to act directly by reactive oxygen species (ROS) scavenger in both wild-type and CL2006 worms, also conferred resistance against oxidative stress in wild-type animals. Furthermore, both treatments delayed Aβ-induced paralysis and decreased AChE activity in CL2006. The protective effect of IPHE against Aβ-induced paralysis was found to be dependent on heat shock factor hsf-1 and FOXO-family transcription factor daf-16, which are respectively involved in aging-related processes and chaperone synthesis, while that of caffeine was dependent only on daf-16. Mechanistically, IPHE and caffeine decreased the levels of Aβ mRNA in the CL2006 worms; however, only IPHE induced expression of the heat shock chaperonin hsp-16.2, involved in protein homeostasis. The results were overall better when treated with IPHE than with caffeine.

Male Sterility in Maize after Transient Heat Stress during the Tetrad Stage of Pollen Development

Shifts in the duration and intensity of ambient temperature impair plant development and reproduction, particularly male gametogenesis. Stress exposure causes meiotic defects or premature spore abortion in male reproductive organs, leading to male sterility. However, little is known about the mechanisms underlying stress and male sterility. To elucidate these mechanisms, we imposed a moderate transient heat stress on maize (Zea mays) plants at the tetrad stage of pollen development. After completion of pollen development at optimal conditions, stress responses were assessed in mature pollen. Transient heat stress resulted in reduced starch content, decreased enzymatic activity, and reduced pollen germination, resulting in sterility. A transcriptomic comparison pointed toward misregulation of starch, lipid, and energy biosynthesis-related genes. Metabolomic studies showed an increase of Suc and its monosaccharide components, as well as a reduction in pyruvate. Lipidomic analysis showed increased levels of unsaturated fatty acids and decreased levels of saturated fatty acids. In contrast, the majority of genes involved in developmental processes such as those required for auxin and unfolded protein responses, signaling, and cell wall biosynthesis remained unaltered. It is noteworthy that changes in the regulation of transcriptional and metabolic pathway genes, as well as heat stress proteins, remained altered even though pollen could recover during further development at optimal conditions. In conclusion, our findings demonstrate that a short moderate heat stress during the highly susceptible tetrad stage strongly affects basic metabolic pathways and thus generates germination-defective pollen, ultimately leading to severe yield losses in maize.

Genome structure predicts modular transcriptome responses to genetic and environmental conditions

Understanding the plasticity, robustness and modularity of transcriptome expression to genetic and environmental conditions is crucial to deciphering how organisms adapt in nature. To test how genome architecture influences transcriptome profiles, we quantified expression responses for distinct temperature-adapted genotypes of the nematode Caenorhabditis briggsae when exposed to chronic temperature stresses throughout development. We found that 56% of the 8,795 differentially expressed genes show genotype-specific changes in expression in response to temperature (genotype-by-environment interactions, GxE). Most genotype-specific responses occur under heat stress, indicating that cold vs. heat stress responses involve distinct genomic architectures. The 22 co-expression modules that we identified differ in their enrichment of genes with genetic vs. environmental vs. interaction effects, as well as their genomic spatial distributions, functional attributes and rates of molecular evolution at the sequence level. Genes in modules enriched for simple effects of either genotype or temperature alone tend to evolve especially rapidly, consistent with disproportionate influence of adaptation or weaker constraint on these subsets of loci. Chromosome-scale heterogeneity in nucleotide polymorphism, however, rather than the scale of individual genes predominates as the source of genetic differences among expression profiles, and natural selection regimes are largely decoupled between coding sequences and noncoding flanking sequences that contain cis-regulatory elements. These results illustrate how the form of transcriptome modularity and genome structure contribute to predictable profiles of evolutionary change.

RNA processing errors triggered by cadmium and integrator complex disruption are signals for environmental stress

BACKGROUND

Adaptive responses to stress are essential for cell and organismal survival. In metazoans, little is known about the impact of environmental stress on RNA homeostasis.

RESULTS

By studying the regulation of a cadmium-induced gene named numr-1 in Caenorhabditis elegans, we discovered that disruption of RNA processing acts as a signal for environmental stress. We find that NUMR-1 contains motifs common to RNA splicing factors and influences RNA splicing in vivo. A genome-wide screen reveals that numr-1 is strongly and specifically induced by silencing of genes that function in basal RNA metabolism including subunits of the metazoan integrator complex. Human integrator processes snRNAs for functioning with splicing factors, and we find that silencing of C. elegans integrator subunits disrupts snRNA processing, causes aberrant pre-mRNA splicing, and induces the heat shock response. Cadmium, which also strongly induces numr-1, has similar effects on RNA and the heat shock response. Lastly, we find that heat shock factor-1 is required for full numr-1 induction by cadmium.

CONCLUSION

Our results are consistent with a model in which disruption of integrator processing of RNA acts as a molecular damage signal initiating an adaptive stress response mediated by heat shock factor-1. When numr-1 is induced via this pathway in C. elegans, its function in RNA metabolism may allow it to mitigate further damage and thereby promote tolerance to cadmium.

Duplication of hsp-110 Is Implicated in Differential Success of Globodera Species under Climate Change

Managing the emergence and spread of crop pests and pathogens is essential for global food security. Understanding how organisms have adapted to their native climate is key to predicting the impact of climate change. The potato cyst nematodes Globodera pallida and G. rostochiensis are economically important plant pathogens that cause yield losses of up to 50% in potato. The two species have different thermal optima that may relate to differences in the altitude of their regions of origin in the Andes. Here, we demonstrate that juveniles of G. pallida are less able to recover from heat stress than those of G. rostochiensis. Genome-wide analysis revealed that while both Globodera species respond to heat stress by induction of various protective heat-inducible genes, G. pallida experiences heat stress at lower temperatures. We use C. elegans as a model to demonstrate the dependence of the heat stress response on expression of Heat Shock Factor-1 (HSF-1). Moreover, we show that hsp-110 is induced by heat stress in G. rostochiensis, but not in the less thermotolerant G. pallida. Sequence analysis revealed that this gene and its promoter was duplicated in G. rostochiensis and acquired thermoregulatory properties. We show that hsp-110 is required for recovery from acute thermal stress in both C. elegans and in G. rostochiensis. Our findings point towards an underlying molecular mechanism that allows the differential expansion of one species relative to another closely related species under current climate change scenarios. Similar mechanisms may be true of other invertebrate species with pest status.

HSB-1 Inhibition and HSF-1 Overexpression Trigger Overlapping Transcriptional Changes To Promote Longevity in Caenorhabditis elegans

Heat shock factor 1 (HSF-1) is a component of the heat shock response pathway that is induced by cytoplasmic proteotoxic stress. In addition to its role in stress response, HSF-1 also acts as a key regulator of the rate of organismal aging. Overexpression of HSF-1 promotes longevity in C. elegans via mechanisms that remain less understood. Moreover, genetic ablation of a negative regulator of HSF-1, termed as heat shock factor binding protein 1 (HSB-1), results in hsf-1-dependent life span extension in animals. Here we show that in the absence of HSB-1, HSF-1 acquires increased DNA binding activity to its genomic target sequence. Using RNA-Seq to compare the gene expression profiles of the hsb-1 mutant and hsf-1 overexpression strains, we found that while more than 1,500 transcripts show ≥1.5-fold upregulation due to HSF-1 overexpression, HSB-1 inhibition alters the expression of less than 500 genes in C. elegans. Roughly half of the differentially regulated transcripts in the hsb-1 mutant have altered expression also in hsf-1 overexpressing animals, with a strongly correlated fold-expression pattern between the two strains. In addition, genes that are upregulated via both HSB-1 inhibition and HSF-1 overexpression include numerous DAF-16 targets that have known functions in longevity regulation. This study identifies how HSB-1 acts as a specific regulator of the transactivation potential of HSF-1 in non-stressed conditions, thus providing a detailed understanding of the role of HSB-1/HSF-1 signaling pathway in transcriptional regulation and longevity in C. elegans.

Studying Transcriptional Enhancers: The Founder Fallacy, Validation Creep, and Other Biases

Transcriptional enhancers play a major role in regulating metazoan gene expression. Recent developments in genomics and next-generation sequencing have accelerated and revitalized the study of this important class of sequence elements. Increased interest and attention, however, has also led to troubling trends in the enhancer literature. In this Opinion, I describe some of these issues and show how they arise from shifting and nonuniform enhancer definitions, and genome-era biases. I discuss how they can lead to interpretative errors and an unduly narrow focus on certain aspects of enhancer biology to the potential exclusion of others.

Prefoldin 6 mediates longevity response from heat shock factor 1 to FOXO in C. elegans

Heat shock factor 1 (HSF-1) and forkhead box O (FOXO) are key transcription factors that protect cells from various stresses. In Caenorhabditis elegans, HSF-1 and FOXO together promote a long life span when insulin/IGF-1 signaling (IIS) is reduced. However, it remains poorly understood how HSF-1 and FOXO cooperate to confer IIS-mediated longevity. Here, we show that prefoldin 6 (PFD-6), a component of the molecular chaperone prefoldin-like complex, relays longevity response from HSF-1 to FOXO under reduced IIS. We found that PFD-6 was specifically required for reduced IIS-mediated longevity by acting in the intestine and hypodermis. We showed that HSF-1 increased the levels of PFD-6 proteins, which in turn directly bound FOXO and enhanced its transcriptional activity. Our work suggests that the prefoldin-like chaperone complex mediates longevity response from HSF-1 to FOXO to increase the life span in animals with reduced IIS.

Pyoverdine, a siderophore from Pseudomonas aeruginosa, translocates into C. elegans, removes iron, and activates a distinct host response

Pseudomonas aeruginosa, a re-emerging, opportunistic human pathogen, encodes a variety of virulence determinants. Pyoverdine, a siderophore produced by this bacterium, is essential for pathogenesis in mammalian infections. This observation is generally attributed to its roles in acquiring iron and/or regulating other virulence factors. Here we report that pyoverdine translocates into the host, where it binds and extracts iron. Pyoverdine-mediated iron extraction damages host mitochondria, disrupting their function and triggering mitochondrial turnover via autophagy. The host detects this damage via a conserved mitochondrial surveillance pathway mediated by the ESRE network. Our findings illuminate the pathogenic mechanisms of pyoverdine and highlight the importance of this bacterial product in host-pathogen interactions.

Temporal dynamics of gene expression in heat-stressed Caenorhabditis elegans

There is considerable insight into pathways and genes associated with heat-stress conditions. Most genes involved in stress response have been identified using mutant screens or gene knockdowns. Yet, there is limited understanding of the temporal dynamics of global gene expression in stressful environments. Here, we studied global gene expression profiles during 12 hours of heat stress in the nematode C. elegans. Using a high-resolution time series of increasing stress exposures, we found a distinct shift in gene expression patterns between 3–4 hours into the stress response, separating an initially highly dynamic phase from a later relatively stagnant phase. This turning point in expression dynamics coincided with a phenotypic turning point, as shown by a strong decrease in movement, survival and, progeny count in the days following the stress. Both detectable at transcriptional and phenotypic level, this study pin-points a relatively small time frame during heat stress at which enough damage is accumulated, making it impossible to recover the next few days.

Mechanism of biofilm-mediated stress resistance and lifespan extension in C. elegans

Bacteria naturally form communities of cells known as biofilms. However the physiological roles of biofilms produced by non-pathogenic microbiota remain largely unknown. To assess the impact of a biofilm on host physiology we explored the effect of several non-pathogenic biofilm-forming bacteria on Caenorhabditis elegans. We show that biofilm formation by Bacillus subtilis, Lactobacillus rhamnosus and Pseudomonas fluorescens induces C. elegans stress resistance. Biofilm also protects against pathogenic infection and prolongs lifespan. Total mRNA analysis identified a set of host genes that are upregulated in response to biofilm formation by B. subtilis. We further demonstrate that mtl-1 is responsible for the biofilm-mediated increase in oxidative stress resistance and lifespan extension. Induction of mtl-1 and hsp-70 promotes biofilm-mediated thermotolerance. ilys-2 activity accounts for biofilm-mediated resistance to Pseudomonas aeruginosa killing. These results reveal the importance of non-pathogenic biofilms for host physiology and provide a framework to study commensal biofilms in higher organisms.

Conserved signaling mechanisms in Drosophila heart development

Signal transduction through multiple distinct pathways regulates and orchestrates the numerous biological processes comprising heart development. This review outlines the roles of the FGFR, EGFR, Wnt, BMP, Notch, Hedgehog, Slit/Robo, and other signaling pathways during four sequential phases of Drosophila cardiogenesis-mesoderm migration, cardiac mesoderm establishment, differentiation of the cardiac mesoderm into distinct cardiac cell types, and morphogenesis of the heart and its lumen based on the proper positioning and cell shape changes of these differentiated cardiac cells-and illustrates how these same cardiogenic roles are conserved in vertebrates. Mechanisms bringing about the regulation and combinatorial integration of these diverse signaling pathways in Drosophila are also described. This synopsis of our present state of knowledge of conserved signaling pathways in Drosophila cardiogenesis and the means by which it was acquired should facilitate our understanding of and investigations into related processes in vertebrates. Developmental Dynamics 246:641-656, 2017. © 2017 Wiley Periodicals, Inc.

Contribution of trans regulatory eQTL to cryptic genetic variation in C. elegans

Background

Cryptic genetic variation (CGV) is the hidden genetic variation that can be unlocked by perturbing normal conditions. CGV can drive the emergence of novel complex phenotypes through changes in gene expression. Although our theoretical understanding of CGV has thoroughly increased over the past decade, insight into polymorphic gene expression regulation underlying CGV is scarce. Here we investigated the transcriptional architecture of CGV in response to rapid temperature changes in the nematode Caenorhabditis elegans. We analyzed regulatory variation in gene expression (and mapped eQTL) across the course of a heat stress and recovery response in a recombinant inbred population.

Results

We measured gene expression over three temperature treatments: i) control, ii) heat stress, and iii) recovery from heat stress. Compared to control, exposure to heat stress affected the transcription of 3305 genes, whereas 942 were affected in recovering animals. These affected genes were mainly involved in metabolism and reproduction. The gene expression pattern in recovering animals resembled both the control and the heat-stress treatment. We mapped eQTL using the genetic variation of the recombinant inbred population and detected 2626 genes with an eQTL in the heat-stress treatment, 1797 in the control, and 1880 in the recovery. The cis-eQTL were highly shared across treatments. A considerable fraction of the trans-eQTL (40–57%) mapped to 19 treatment specific trans-bands. In contrast to cis-eQTL, trans-eQTL were highly environment specific and thus cryptic. Approximately 67% of the trans-eQTL were only induced in a single treatment, with heat-stress showing the most unique trans-eQTL.

Conclusions

These results illustrate the highly dynamic pattern of CGV across three different environmental conditions that can be evoked by a stress response over a relatively short time-span (2 h) and that CGV is mainly determined by response related trans regulatory eQTL.

Electronic supplementary material

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

A conserved mitochondrial surveillance pathway is required for defense against Pseudomonas aeruginosa

All living organisms exist in a precarious state of homeostasis that requires constant maintenance. A wide variety of stresses, including hypoxia, heat, and infection by pathogens perpetually threaten to imbalance this state. Organisms use a battery of defenses to mitigate damage and restore normal function. Previously, we described a Caenorhabditis elegans-Pseudomonas aeruginosa assay (Liquid Killing) in which toxicity to the host is dependent upon the secreted bacterial siderophore pyoverdine. Although pyoverdine is also indispensable for virulence in mammals, its cytological effects are unclear. We used genetics, transcriptomics, and a variety of pathogen and chemical exposure assays to study the interactions between P. aeruginosa and C. elegans. Although P. aeruginosa can kill C. elegans through at least 5 different mechanisms, the defense responses activated by Liquid Killing are specific and selective and have little in common with innate defense mechanisms against intestinal colonization. Intriguingly, the defense response utilizes the phylogenetically-conserved ESRE (Ethanol and Stress Response Element) network, which we and others have previously shown to mitigate damage from a variety of abiotic stresses. This is the first report of this networks involvement in innate immunity, and indicates that host innate immune responses overlap with responses to abiotic stresses. The upregulation of the ESRE network in C. elegans is mediated in part by a family of bZIP proteins (including ZIP-2, ZIP-4, CEBP-1, and CEBP-2) that have overlapping and unique functions. Our data convincingly show that, following exposure to P. aeruginosa, the ESRE defense network is activated by mitochondrial damage, and that mitochondrial damage also leads to ESRE activation in mammals. This establishes a role for ESRE in a phylogenetically-conserved mitochondrial surveillance system important for stress response and innate immunity.

Conserved Transcription Factors Steer Growth-Related Genomic Programs in Daphnia

Ecological genomics aims to understand the functional association between environmental gradients and the genes underlying adaptive traits. Many genes that are identified by genome-wide screening in ecologically relevant species lack functional annotations. Although gene functions can be inferred from sequence homology, such approaches have limited power. Here, we introduce ecological regulatory genomics by presenting an ontology-free gene prioritization method. Specifically, our method combines transcriptome profiling with high-throughput cis-regulatory sequence analysis in the water fleas Daphnia pulex and Daphnia magna. It screens coexpressed genes for overrepresented DNA motifs that serve as transcription factor binding sites, thereby providing insight into conserved transcription factors and gene regulatory networks shaping the expression profile. We first validated our method, called Daphnia-cisTarget, on a D. pulex heat shock data set, which revealed a network driven by the heat shock factor. Next, we performed RNA-Seq in D. magna exposed to the cyanobacterium Microcystis aeruginosa. Daphnia-cisTarget identified coregulated gene networks that associate with the moulting cycle and potentially regulate life history changes in growth rate and age at maturity. These networks are predicted to be regulated by evolutionary conserved transcription factors such as the homologues of Drosophila Shavenbaby and Grainyhead, nuclear receptors, and a GATA family member. In conclusion, our approach allows prioritising candidate genes in Daphnia without bias towards prior knowledge about functional gene annotation and represents an important step towards exploring the molecular mechanisms of ecological responses in organisms with poorly annotated genomes.

An Enhancer's Length and Composition Are Shaped by Its Regulatory Task

Enhancers drive the gene expression patterns required for virtually every process in metazoans. We propose that enhancer length and transcription factor (TF) binding site composition—the number and identity of TF binding sites—reflect the complexity of the enhancer's regulatory task. In development, we define regulatory task complexity as the number of fates specified in a set of cells at once. We hypothesize that enhancers with more complex regulatory tasks will be longer, with more, but less specific, TF binding sites. Larger numbers of binding sites can be arranged in more ways, allowing enhancers to drive many distinct expression patterns, and therefore cell fates, using a finite number of TF inputs. We compare ~100 enhancers patterning the more complex anterior-posterior (AP) axis and the simpler dorsal-ventral (DV) axis in Drosophila and find that the AP enhancers are longer with more, but less specific binding sites than the (DV) enhancers. Using a set of ~3,500 enhancers, we find enhancer length and TF binding site number again increase with increasing regulatory task complexity. Therefore, to be broadly applicable, computational tools to study enhancers must account for differences in regulatory task.

Identification of β-catenin-interacting proteins in nuclear fractions of native rat collecting duct cells

The gene encoding the aquaporin-2 water channel is regulated transcriptionally in response to vasopressin. In the renal collecting duct, vasopressin stimulates the nuclear translocation and phosphorylation (at Ser⁵⁵²) of β-catenin, a multifunctional protein that acts as a transcriptional coregulator in the nucleus. The purpose of this study was to identify β-catenin-interacting proteins that might be involved in transcriptional regulation in rat inner medullary collecting duct (IMCD) cells, using experimental and computational approaches. We used a standard chromatin immunoprecipitation procedure coupled to mass spectrometry (ChIP-MS) in a nuclear fraction isolated from rat IMCD suspensions. Over four biological replicates, we reproducibly identified 43 β-catenin-binding proteins, including several known β-catenin-binding partners as well as novel interacting proteins. Multiple proteins involved in transcriptional regulation were identified (Taf1, Jup, Tdrd3, Cdh1, Cenpj, and several histones). Many of the identified β-catenin-binding partners were found in prior studies to translocate to the nucleus in response to vasopressin. There was only one DNA-binding transcription factor (TF), specifically Taf1, part of the RNA-polymerase II preinitiation complex. To identify sequence-specific TFs that might interact with β-catenin, Bayes' theorem was used to integrate data from several information sources. The analysis identified several TFs with potential binding sites in the Aqp2 gene promoter that could interact with β-catenin in the regulation of Aqp2 gene transcription, specifically Jun, Junb, Jund, Atf1, Atf2, Mef2d, Usf1, Max, Pou2f1, and Rxra. The findings provide information necessary for modeling the transcriptional response to vasopressin.

Hormetic heat stress and HSF-1 induce autophagy to improve survival and proteostasis in C. elegans

Stress-response pathways have evolved to maintain cellular homeostasis and to ensure the survival of organisms under changing environmental conditions. Whereas severe stress is detrimental, mild stress can be beneficial for health and survival, known as hormesis. Although the universally conserved heat-shock response regulated by transcription factor HSF-1 has been implicated as an effector mechanism, the role and possible interplay with other cellular processes, such as autophagy, remains poorly understood. Here we show that autophagy is induced in multiple tissues of Caenorhabditis elegans following hormetic heat stress or HSF-1 overexpression. Autophagy-related genes are required for the thermoresistance and longevity of animals exposed to hormetic heat shock or HSF-1 overexpression. Hormetic heat shock also reduces the progressive accumulation of PolyQ aggregates in an autophagy-dependent manner. These findings demonstrate that autophagy contributes to stress resistance and hormesis, and reveal a requirement for autophagy in HSF-1-regulated functions in the heat-shock response, proteostasis and ageing.

Mild heat stress has beneficial effects on organismal health and survival. Here, Kumsta et al. show that a mild heat shock and HSF-1 overexpression induce autophagy in multiple tissues of C. elegans and autophagy-related genes are essential for both heat shock-induced and HSF-1–mediated stress resistance and longevity.

Involvement of heat shock proteins on Mn-induced toxicity in Caenorhabditis elegans

BACKGROUND

All living cells display a rapid molecular response to adverse environmental conditions, and the heat shock protein family reflects one such example. Hence, failing to activate heat shock proteins can impair the cellular response. In the present study, we evaluated whether the loss of different isoforms of heat shock protein (hsp) genes in Caenorhabditis elegans would affect their vulnerability to Manganese (Mn) toxicity.

METHODS

We exposed wild type and selected hsp mutant worms to Mn (30 min) and next evaluated further the most susceptible strains. We analyzed survival, protein carbonylation (as a marker of oxidative stress) and Parkinson's disease related gene expression immediately after Mn exposure. Lastly, we observed dopaminergic neurons in wild type worms and in hsp-70 mutants following Mn treatment. Analysis of the data was performed by one-way or two way ANOVA, depending on the case, followed by post-hoc Bonferroni test if the overall p value was less than 0.05.

RESULTS

We verified that the loss of hsp-70, hsp-3 and chn-1 increased the vulnerability to Mn, as exposed mutant worms showed lower survival rate and increased protein oxidation. The importance of hsp-70 against Mn toxicity was then corroborated in dopaminergic neurons, where Mn neurotoxicity was aggravated. The lack of hsp-70 also blocked the transcriptional upregulation of pink1, a gene that has been linked to Parkinson's disease.

CONCLUSIONS

Taken together, our data suggest that Mn exposure modulates heat shock protein expression, particularly HSP-70, in C. elegans. Furthermore, loss of hsp-70 increases protein oxidation and dopaminergic neuronal degeneration following manganese exposure, which is associated with the inhibition of pink1 increased expression, thus potentially exacerbating the vulnerability to this metal.

Reactive oxygen species (ROS) and the heat stress response of Daphnia pulex: ROS-mediated activation of hypoxia-inducible factor 1 (HIF-1) and heat shock factor 1 (HSF-1) and the clustered expression of stress genes

BACKGROUND INFORMATION

Heat stress in ectotherms involves direct (e.g. protein damage) and/or indirect effects (temperature-induced hypoxia and ROS formation), which cause activation of the transcription factors (TF) heat shock factor 1 (HSF-1) and/or hypoxia-inducible factor 1 (HIF-1). The present study focused on the links between stress (ROS) signals, nuclear (n) and cytoplasmic (c) HSF-1/HIF-1 levels, and stress gene expression on mRNA and protein levels (e.g. heat-shock protein 90, HSP90) upon acute heat and ROS (H₂ O₂ ) stress.

RESULTS

Acute heat stress (30°C) evoked fluctuations in ROS level. Different feeding regimens, which affected the glutathione (GSH) level, allowed altering the frequency of ROS fluctuations. Other data showed fluctuation frequency to depend also on ROS production rate. The heat-induced slow or fast ROS fluctuations (at high or low GSH levels) evoked slow or fast fluctuations in the levels of nHIF-1α, nHSF-1 and gene products (mRNAs and protein), albeit after different time delays. Time delays to ROS fluctuations were, for example,shorter for nHIF-1α than for nHSF-1 fluctuations, and nHIF-1α fluctuations preceded and nHSF-1 fluctuations followed fluctuations in HSP90 mRNA level. Cytoplasmic TF levels either changed little (cHIF-1α) or showed a steady increase (cHSF-1). Applying acute H₂ O₂ stress (at 20°C) revealed effects on nHIF-1α and mRNA levels, but no significant effects on nHSF-1 level. Transcriptome data additionally showed coordinated fluctuations of mRNA levels upon acute heat stress, involving mRNAs for HSPs and other stress proteins, with all corresponding genes carrying DNA binding motifs for HIF-1 and HSF-1.

CONCLUSIONS

This study provided evidence for promoting effects of ROS and HIF-1 on early haemoglobin, HIF-1α and HSP90 mRNA expressions upon heat or ROS stress. The increasing cHSF-1 level likely affected nHSF-1 level and later HSP90 mRNA expression.

SIGNIFICANCE

Heat stress evoked ROS fluctuations, with this stress signal forwarded via nHIF-1 and nHSF-1 fluctuations to stress gene expression. The frequency of ROS fluctuations seemed to integrate information about ROS productionrate and GSH antioxidant buffer capacity, resulting in stress protein expression of different speed. Results of this study suggest ROS as early (pre-damage) and protein defects as later (post-damage) stress signals to trigger heat stress responses.

Structuring evolution: biochemical networks and metabolic diversification in birds

Background

Recurrence and predictability of evolution are thought to reflect the correspondence between genomic and phenotypic dimensions of organisms, and the connectivity in deterministic networks within these dimensions. Direct examination of the correspondence between opportunities for diversification imbedded in such networks and realized diversity is illuminating, but is empirically challenging because both the deterministic networks and phenotypic diversity are modified in the course of evolution. Here we overcome this problem by directly comparing the structure of a “global” carotenoid network – comprising of all known enzymatic reactions among naturally occurring carotenoids – with the patterns of evolutionary diversification in carotenoid-producing metabolic networks utilized by birds.

Results

We found that phenotypic diversification in carotenoid networks across 250 species was closely associated with enzymatic connectivity of the underlying biochemical network – compounds with greater connectivity occurred the most frequently across species and were the hotspots of metabolic pathway diversification. In contrast, we found no evidence for diversification along the metabolic pathways, corroborating findings that the utilization of the global carotenoid network was not strongly influenced by history in avian evolution.

Conclusions

The finding that the diversification in species-specific carotenoid networks is qualitatively predictable from the connectivity of the underlying enzymatic network points to significant structural determinism in phenotypic evolution.

Electronic supplementary material

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

Transcription factors GAF and HSF act at distinct regulatory steps to modulate stress-induced gene activation

The coordinated regulation of gene expression at the transcriptional level is fundamental to development and homeostasis. Inducible systems are invaluable when studying transcription because the regulatory process can be triggered instantaneously, allowing the tracking of ordered mechanistic events. Here, we use precision run-on sequencing (PRO-seq) to examine the genome-wide heat shock (HS) response in Drosophila and the function of two key transcription factors on the immediate transcription activation or repression of all genes regulated by HS. We identify the primary HS response genes and the rate-limiting steps in the transcription cycle that GAGA-associated factor (GAF) and HS factor (HSF) regulate. We demonstrate that GAF acts upstream of promoter-proximally paused RNA polymerase II (Pol II) formation (likely at the step of chromatin opening) and that GAF-facilitated Pol II pausing is critical for HS activation. In contrast, HSF is dispensable for establishing or maintaining Pol II pausing but is critical for the release of paused Pol II into the gene body at a subset of highly activated genes. Additionally, HSF has no detectable role in the rapid HS repression of thousands of genes.

HSF-1 activates the ubiquitin proteasome system to promote non-apoptotic developmental cell death in C. elegans

Apoptosis is a prominent metazoan cell death form. Yet, mutations in apoptosis regulators cause only minor defects in vertebrate development, suggesting that another developmental cell death mechanism exists. While some non-apoptotic programs have been molecularly characterized, none appear to control developmental cell culling. Linker-cell-type death (LCD) is a morphologically conserved non-apoptotic cell death process operating in Caenorhabditis elegans and vertebrate development, and is therefore a compelling candidate process complementing apoptosis. However, the details of LCD execution are not known. Here we delineate a molecular-genetic pathway governing LCD in C. elegans. Redundant activities of antagonistic Wnt signals, a temporal control pathway, and mitogen-activated protein kinase kinase signaling control heat shock factor 1 (HSF-1), a conserved stress-activated transcription factor. Rather than protecting cells, HSF-1 promotes their demise by activating components of the ubiquitin proteasome system, including the E2 ligase LET-70/UBE2D2 functioning with E3 components CUL-3, RBX-1, BTBD-2, and SIAH-1. Our studies uncover design similarities between LCD and developmental apoptosis, and provide testable predictions for analyzing LCD in vertebrates.

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

Embryos make numerous new cells as they develop, but also destroy many cells to remove the faulty ones and to ensure that tissues grow to the right size and shape. This deliberate form of cell death must be precisely regulated to prevent too many cells or healthy cells, from being destroyed. Understanding the molecular mechanisms that govern cell death is therefore important for understanding normal development and also human disease.

One well-studied process that leads to cell death is called apoptosis. This process carefully dismantles and breaks down the components of a cell, but does not seem to account for all cell death that occurs during animal development. Recently another developmental cell-death pathway, called the linker-cell-type death, was discovered in a small roundworm called Caenorhabditis elegans. This pathway appears to work in mammalian cells as well, and may help to break down nerve fibers that are not needed. However, many of this pathway’s component parts remained unknown.

Kinet, Malin et al. have now used a combination of genetics and cell biology in C. elegans to uncover the components of linker-cell-type death and to investigate how they interact. The results of these studies revealed a hierarchy of genetic interactions that governs this pathway in C. elegans. One protein called HSF-1 plays a particularly important role. This protein is a transcription factor and it binds to, and regulates, the activities of various genes. HSF-1 usually works in cells to protect them from stress, but Kinet, Malin et al. showed that it instead promotes linker-cell-type death by activating a molecular machine, called the proteasome, that breaks down proteins. The experiments also revealed two proteins (called BTBD-2 and SIAH-1) that may be important for shuttling specific proteins for degradation by the proteasome.

Three signalling pathways that regulate important developmental processes also regulate the activation of linker-cell-type death. Kinet, Malin et al. propose that these signalling pathways do so by working together to activate HSF-1, which in turn activates the genes that lead to the destruction of cells by the proteasome.

A future challenge is to understand in more detail how the more recently discovered cell death pathway actually kills cells. Further work could also explore how HSF-1, a protein that normally protects cells, is transformed into a cell-killing protein.

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

Interplay among Drosophila transcription factors Ets21c, Fos and Ftz-F1 drives JNK-mediated tumor malignancy

Cancer initiation and maintenance of the transformed cell state depend on altered cellular signaling and aberrant activities of transcription factors (TFs) that drive pathological gene expression in response to cooperating genetic lesions. Deciphering the roles of interacting TFs is therefore central to understanding carcinogenesis and for designing cancer therapies. Here, we use an unbiased genomic approach to define a TF network that triggers an abnormal gene expression program promoting malignancy of clonal tumors, generated in Drosophila imaginal disc epithelium by gain of oncogenic Ras (Ras^V12) and loss of the tumor suppressor Scribble (scrib¹). We show that malignant transformation of the ras^V12scrib¹ tumors requires TFs of distinct families, namely the bZIP protein Fos, the ETS-domain factor Ets21c and the nuclear receptor Ftz-F1, all acting downstream of Jun-N-terminal kinase (JNK). Depleting any of the three TFs improves viability of tumor-bearing larvae, and this positive effect can be enhanced further by their combined removal. Although both Fos and Ftz-F1 synergistically contribute to ras^V12scrib¹ tumor invasiveness, only Fos is required for JNK-induced differentiation defects and Matrix metalloprotease (MMP1) upregulation. In contrast, the Fos-dimerizing partner Jun is dispensable for JNK to exert its effects in ras^V12scrib¹ tumors. Interestingly, Ets21c and Ftz-F1 are transcriptionally induced in these tumors in a JNK- and Fos-dependent manner, thereby demonstrating a hierarchy within the tripartite TF network, with Fos acting as the most upstream JNK effector. Of the three TFs, only Ets21c can efficiently substitute for loss of polarity and cooperate with Ras^V12 in inducing malignant clones that, like ras^V12scrib¹ tumors, invade other tissues and overexpress MMP1 and the Drosophila insulin-like peptide 8 (Dilp8). While ras^V12ets21c tumors require JNK for invasiveness, the JNK activity is dispensable for their growth. In conclusion, our study delineates both unique and overlapping functions of distinct TFs that cooperatively promote aberrant expression of target genes, leading to malignant tumor phenotypes.

Summary: This study provides genetic evidence that malignancy driven by oncogenic Ras and loss of polarity requires transcription factors of three distinct protein families, acting in synergy downstream of JNK signaling.

RNA helicase HEL-1 promotes longevity by specifically activating DAF-16/FOXO transcription factor signaling in Caenorhabditis elegans

The homeostatic maintenance of the genomic DNA is crucial for regulating aging processes. However, the role of RNA homeostasis in aging processes remains unknown. RNA helicases are a large family of enzymes that regulate the biogenesis and homeostasis of RNA. However, the functional significance of RNA helicases in aging has not been explored. Here, we report that a large fraction of RNA helicases regulate the lifespan of Caenorhabditis elegans. In particular, we show that a DEAD-box RNA helicase, helicase 1 (HEL-1), promotes longevity by specifically activating the DAF-16/forkhead box O (FOXO) transcription factor signaling pathway. We find that HEL-1 is required for the longevity conferred by reduced insulin/insulin-like growth factor 1 (IGF-1) signaling (IIS) and is sufficient for extending lifespan. We further show that the expression of HEL-1 in the intestine and neurons contributes to longevity. HEL-1 enhances the induction of a large fraction of DAF-16 target genes. Thus, the RNA helicase HEL-1 appears to promote longevity in response to decreased IIS as a transcription coregulator of DAF-16. Because HEL-1 and IIS are evolutionarily well conserved, a similar mechanism for longevity regulation via an RNA helicase-dependent regulation of FOXO signaling may operate in mammals, including humans.

Identifying transcriptional cis-regulatory modules in animal genomes

Gene expression is regulated through the activity of transcription factors (TFs) and chromatin-modifying proteins acting on specific DNA sequences, referred to as cis-regulatory elements. These include promoters, located at the transcription initiation sites of genes, and a variety of distal cis-regulatory modules (CRMs), the most common of which are transcriptional enhancers. Because regulated gene expression is fundamental to cell differentiation and acquisition of new cell fates, identifying, characterizing, and understanding the mechanisms of action of CRMs is critical for understanding development. CRM discovery has historically been challenging, as CRMs can be located far from the genes they regulate, have few readily identifiable sequence characteristics, and for many years were not amenable to high-throughput discovery methods. However, the recent availability of complete genome sequences and the development of next-generation sequencing methods have led to an explosion of both computational and empirical methods for CRM discovery in model and nonmodel organisms alike. Experimentally, CRMs can be identified through chromatin immunoprecipitation directed against TFs or histone post-translational modifications, identification of nucleosome-depleted 'open' chromatin regions, or sequencing-based high-throughput functional screening. Computational methods include comparative genomics, clustering of known or predicted TF-binding sites, and supervised machine-learning approaches trained on known CRMs. All of these methods have proven effective for CRM discovery, but each has its own considerations and limitations, and each is subject to a greater or lesser number of false-positive identifications. Experimental confirmation of predictions is essential, although shortcomings in current methods suggest that additional means of validation need to be developed. For further resources related to this article, please visit the WIREs website.

CONFLICT OF INTEREST

The authors have declared no conflicts of interest for this article.

Enhancer modeling uncovers transcriptional signatures of individual cardiac cell states in Drosophila

Here we used discriminative training methods to uncover the chromatin, transcription factor (TF) binding and sequence features of enhancers underlying gene expression in individual cardiac cells. We used machine learning with TF motifs and ChIP data for a core set of cardiogenic TFs and histone modifications to classify Drosophila cell-type-specific cardiac enhancer activity. We show that the classifier models can be used to predict cardiac cell subtype cis-regulatory activities. Associating the predicted enhancers with an expression atlas of cardiac genes further uncovered clusters of genes with transcription and function limited to individual cardiac cell subtypes. Further, the cell-specific enhancer models revealed chromatin, TF binding and sequence features that distinguish enhancer activities in distinct subsets of heart cells. Collectively, our results show that computational modeling combined with empirical testing provides a powerful platform to uncover the enhancers, TF motifs and gene expression profiles which characterize individual cardiac cell fates.

The balanced regulation of Hsc70 by DNJ-13 and UNC-23 is required for muscle functionality

The molecular chaperone Hsc70 assists in the folding of non-native proteins together with its J domain- and BAG domain-containing cofactors. In Caenorhabditis elegans, two BAG domain-containing proteins can be identified, one of them being UNC-23, whose mutation induces severe motility dysfunctions. Using reporter strains, we find that the full-length UNC-23, in contrast to C-terminal fragments, localizes specifically to the muscular attachment sites. C-terminal fragments of UNC-23 instead perform all Hsc70-related functions, like ATPase stimulation and regulation of folding activity, albeit with lower affinity than BAG-1. Interestingly, overexpression of CFP-Hsc70 can induce muscular defects in wild-type nematodes that phenocopy the knockout of its cofactor UNC-23. Strikingly, the motility dysfunction in the unc-23 mutated strain can be cured specifically by down-regulation of the antagonistic Hsc70 cochaperone DNJ-13, implying that the severe phenotype is caused by misregulation of the Hsc70 cycle. These findings point out that the balanced action of cofactors in the ATP-driven cycle of Hsc70 is crucial for the contribution of Hsc70 to muscle functionality.

Ubiquitin-mediated response to microsporidia and virus infection in C. elegans

Microsporidia comprise a phylum of over 1400 species of obligate intracellular pathogens that can infect almost all animals, but little is known about the host response to these parasites. Here we use the whole-animal host C. elegans to show an in vivo role for ubiquitin-mediated response to the microsporidian species Nematocida parisii, as well to the Orsay virus, another natural intracellular pathogen of C. elegans. We analyze gene expression of C. elegans in response to N. parisii, and find that it is similar to response to viral infection. Notably, we find an upregulation of SCF ubiquitin ligase components, such as the cullin ortholog cul-6, which we show is important for ubiquitin targeting of N. parisii cells in the intestine. We show that ubiquitylation components, the proteasome, and the autophagy pathway are all important for defense against N. parisii infection. We also find that SCF ligase components like cul-6 promote defense against viral infection, where they have a more robust role than against N. parisii infection. This difference may be due to suppression of the host ubiquitylation system by N. parisii: when N. parisii is crippled by anti-microsporidia drugs, the host can more effectively target pathogen cells for ubiquitylation. Intriguingly, inhibition of the ubiquitin-proteasome system (UPS) increases expression of infection-upregulated SCF ligase components, indicating that a trigger for transcriptional response to intracellular infection by N. parisii and virus may be perturbation of the UPS. Altogether, our results demonstrate an in vivo role for ubiquitin-mediated defense against microsporidian and viral infections in C. elegans.

Integrin-linked kinase modulates longevity and thermotolerance in C. elegans through neuronal control of HSF-1

Integrin-signaling complexes play important roles in cytoskeletal organization and cell adhesion in many species. Components of the integrin-signaling complex have been linked to aging in both Caenorhabditis elegans and Drosophila melanogaster, but the mechanism underlying this function is unknown. Here, we investigated the role of integrin-linked kinase (ILK), a key component of the integrin-signaling complex, in lifespan determination. We report that genetic reduction of ILK in both C. elegans and Drosophila increased resistance to heat stress, and led to lifespan extension in C. elegans without majorly affecting cytoskeletal integrity. In C. elegans, longevity and thermotolerance induced by ILK depletion was mediated by heat-shock factor-1 (HSF-1), a major transcriptional regulator of the heat-shock response (HSR). Reduction in ILK levels increased hsf-1 transcription and activation, and led to enhanced expression of a subset of genes with roles in the HSR. Moreover, induction of HSR-related genes, longevity and thermotolerance caused by ILK reduction required the thermosensory neurons AFD and interneurons AIY, which are known to play a critical role in the canonical HSR. Notably, ILK was expressed in neighboring neurons, but not in AFD or AIY, implying that ILK reduction initiates cell nonautonomous signaling through thermosensory neurons to elicit a noncanonical HSR. Our results thus identify HSF-1 as a novel effector of the organismal response to reduced ILK levels and show that ILK inhibition regulates HSF-1 in a cell nonautonomous fashion to enhance stress resistance and lifespan in C. elegans.

Detecting clustering and ordering binding patterns among transcription factors via point process models

MOTIVATION

Recent development in ChIP-Seq technology has generated binding data for many transcription factors (TFs) in various cell types and cellular conditions. This opens great opportunities for studying combinatorial binding patterns among a set of TFs active in a particular cellular condition, which is a key component for understanding the interaction between TFs in gene regulation.

RESULTS

As a first step to the identification of combinatorial binding patterns, we develop statistical methods to detect clustering and ordering patterns among binding sites (BSs) of a pair of TFs. Testing procedures based on Ripley's K-function and its generalizations are developed to identify binding patterns from large collections of BSs in ChIP-Seq data. We have applied our methods to the ChIP-Seq data of 91 pairs of TFs in mouse embryonic stem cells. Our methods have detected clustering binding patterns between most TF pairs, which is consistent with the findings in the literature, and have identified significant ordering preferences, relative to the direction of target gene transcription, among the BSs of seven TFs. More interestingly, our results demonstrate that the identified clustering and ordering binding patterns between TFs are associated with the expression of the target genes. These findings provide new insights into co-regulation between TFs.

AVAILABILITY AND IMPLEMENTATION

See 'www.stat.ucla.edu/∼zhou/TFKFunctions/' for source code.

Hormone-responsive enhancer-activity maps reveal predictive motifs, indirect repression, and targeting of closed chromatin

Steroid hormones act as important developmental switches, and their nuclear receptors regulate many genes. However, few hormone-dependent enhancers have been characterized, and important aspects of their sequence architecture, cell-type-specific activating and repressing functions, or the regulatory roles of their chromatin structure have remained unclear. We used STARR-seq, a recently developed enhancer-screening assay, and ecdysone signaling in two different Drosophila cell types to derive genome-wide hormone-dependent enhancer-activity maps. We demonstrate that enhancer activation depends on cis-regulatory motif combinations that differ between cell types and can predict cell-type-specific ecdysone targeting. Activated enhancers are often not accessible prior to induction. Enhancer repression following hormone treatment seems independent of receptor motifs and receptor binding to the enhancer, as we show using ChIP-seq, but appears to rely on motifs for other factors, including Eip74. Our strategy is applicable to study signal-dependent enhancers for different pathways and across organisms.

Imogene: identification of motifs and cis-regulatory modules underlying gene co-regulation

Cis-regulatory modules (CRMs) and motifs play a central role in tissue and condition-specific gene expression. Here we present Imogene, an ensemble of statistical tools that we have developed to facilitate their identification and implemented in a publicly available software. Starting from a small training set of mammalian or fly CRMs that drive similar gene expression profiles, Imogene determines de novocis-regulatory motifs that underlie this co-expression. It can then predict on a genome-wide scale other CRMs with a regulatory potential similar to the training set. Imogene bypasses the need of large datasets for statistical analyses by making central use of the information provided by the sequenced genomes of multiple species, based on the developed statistical tools and explicit models for transcription factor binding site evolution. We test Imogene on characterized tissue-specific mouse developmental CRMs. Its ability to identify CRMs with the same specificity based on its de novo created motifs is comparable to that of previously evaluated ‘motif-blind’ methods. We further show, both in flies and in mammals, that Imogene de novo generated motifs are sufficient to discriminate CRMs related to different developmental programs. Notably, purely relying on sequence data, Imogene performs as well in this discrimination task as a previously reported learning algorithm based on Chromatin Immunoprecipitation (ChIP) data for multiple transcription factors at multiple developmental stages.

Quantitative modeling of a gene's expression from its intergenic sequence

Modeling a gene's expression from its intergenic locus and trans-regulatory context is a fundamental goal in computational biology. Owing to the distributed nature of cis-regulatory information and the poorly understood mechanisms that integrate such information, gene locus modeling is a more challenging task than modeling individual enhancers. Here we report the first quantitative model of a gene's expression pattern as a function of its locus. We model the expression readout of a locus in two tiers: 1) combinatorial regulation by transcription factors bound to each enhancer is predicted by a thermodynamics-based model and 2) independent contributions from multiple enhancers are linearly combined to fit the gene expression pattern. The model does not require any prior knowledge about enhancers contributing toward a gene's expression. We demonstrate that the model captures the complex multi-domain expression patterns of anterior-posterior patterning genes in the early Drosophila embryo. Altogether, we model the expression patterns of 27 genes; these include several gap genes, pair-rule genes, and anterior, posterior, trunk, and terminal genes. We find that the model-selected enhancers for each gene overlap strongly with its experimentally characterized enhancers. Our findings also suggest the presence of sequence-segments in the locus that would contribute ectopic expression patterns and hence were "shut down" by the model. We applied our model to identify the transcription factors responsible for forming the stripe boundaries of the studied genes. The resulting network of regulatory interactions exhibits a high level of agreement with known regulatory influences on the target genes. Finally, we analyzed whether and why our assumption of enhancer independence was necessary for the genes we studied. We found a deterioration of expression when binding sites in one enhancer were allowed to influence the readout of another enhancer. Thus, interference between enhancer activities was a possible factor necessitating enhancer independence in our model.

The conserved PBAF nucleosome-remodeling complex mediates the response to stress in Caenorhabditis elegans

To adapt to stress, cells must undergo major changes in their gene expression profiles. We have previously described a largely uncharacterized stress response pathway in Caenorhabditis elegans that acts through an evolutionarily conserved motif, termed ESRE, for ethanol and stress response element. We characterize here the requirements for ESRE gene expression and show that the ESRE network is regulated by a conserved SWI/SNF family nucleosome remodeling complex termed PBAF. Depletion of PBAF subunits SWSN-7/BAF200 and PBRM-1/BAF180 results in decreased expression of ESRE genes and increased sensitivity to thermal stress. When overexpressed, SWSN-7/BAF200 and PBRM-1/BAF180 led to increased ESRE transcription, enhanced thermotolerance, and induction of a nuclear ESRE-binding activity. Our data support a model in which PBAF is recruited by an ESRE-binding protein to genomic ESRE sites. We also show that the closely related SWI/SNF complex, BAF, which regulates stress induction through DAF-16/FOXO, does not contribute to ESRE gene expression or bind directly to ESRE sites. To our knowledge, this is the first report demonstrating direct and specific regulation of a stress response network by the PBAF nucleosome-remodeling complex in vivo in metazoa. In addition, we show that PBAF cooperates with the histone demethylase, JMJC-1/NO66, to promote expression of ESRE genes following stress.