Conserved Motifs and Prediction of Regulatory Modules in Caenorhabditis elegans

Towards a comprehensive regulatory map of Mammalian Genomes

Genome mapping studies have generated a nearly complete collection of genes for the human genome, but we still lack an equivalently vetted inventory of human regulatory sequences. Cis-regulatory modules (CRMs) play important roles in controlling when, where, and how much a gene is expressed. We developed a training data-free CRM-prediction algorithm, the Mammalian Regulatory MOdule Detector (MrMOD) for accurate CRM prediction in mammalian genomes. MrMOD provides genome position-fixed CRM models similar to the fixed gene models for the mouse and human genomes using only genomic sequences as the inputs with one adjustable parameter - the significance p-value. Importantly, MrMOD predicts a comprehensive set of high-resolution CRMs in the mouse and human genomes including all types of regulatory modules not limited to any tissue, cell type, developmental stage, or condition. We computationally validated MrMOD predictions used a compendium of 21 orthogonal experimental data sets including thousands of experimentally defined CRMs and millions of putative regulatory elements derived from hundreds of different tissues, cell types, and stimulus conditions obtained from multiple databases. In ovo transgenic reporter assay demonstrates the power of our prediction in guiding experimental design. We analyzed CRMs located in the chromosome 17 using unsupervised machine learning and identified groups of CRMs with multiple lines of evidence supporting their functionality, linking CRMs with upstream binding transcription factors and downstream target genes. Our work provides a comprehensive base pair resolution annotation of the functional regulatory elements and non-functional regions in the mammalian genomes.

Males induce premature demise of the opposite sex by multifaceted strategies

Interactions between the sexes negatively impact health in many species. In Caenorhabditis, males shorten the lifespan of the opposite sex-hermaphrodites or females. Here we use transcriptomic profiling and targeted screens to systematically uncover conserved genes involved in male-induced demise in C. elegans. Some genes (for example, delm-2, acbp-3), when knocked down, are specifically protective against male-induced demise. Others (for example, sri-40), when knocked down, extend lifespan with and without males, suggesting general mechanisms of protection. In contrast, many classical long-lived mutants are impacted more negatively than wild type by the presence of males, highlighting the importance of sexual environment for longevity. Interestingly, genes induced by males are triggered by specific male components (seminal fluid, sperm and pheromone), and manipulating these genes in combination in hermaphrodites induces stronger protection. One of these genes, the conserved ion channel delm-2, acts in the nervous system and intestine to regulate lipid metabolism. Our analysis reveals striking differences in longevity in single sex versus mixed sex environments and uncovers elaborate strategies elicited by sexual interactions that could extend to other species.

Evolutionary Dynamics of the SKN-1 → MED → END-1,3 Regulatory Gene Cascade in Caenorhabditis Endoderm Specification

Gene regulatory networks and their evolution are important in the study of animal development. In the nematode, Caenorhabditis elegans, the endoderm (gut) is generated from a single embryonic precursor, E. Gut is specified by the maternal factor SKN-1, which activates the MED → END-1,3 → ELT-2,7 cascade of GATA transcription factors. In this work, genome sequences from over two dozen species within the Caenorhabditis genus are used to identify MED and END-1,3 orthologs. Predictions are validated by comparison of gene structure, protein conservation, and putative cis-regulatory sites. All three factors occur together, but only within the Elegans supergroup, suggesting they originated at its base. The MED factors are the most diverse and exhibit an unexpectedly extensive gene amplification. In contrast, the highly conserved END-1 orthologs are unique in nearly all species and share extended regions of conservation. The END-1,3 proteins share a region upstream of their zinc finger and an unusual amino-terminal poly-serine domain exhibiting high codon bias. Compared with END-1, the END-3 proteins are otherwise less conserved as a group and are typically found as paralogous duplicates. Hence, all three factors are under different evolutionary constraints. Promoter comparisons identify motifs that suggest the SKN-1, MED, and END factors function in a similar gut specification network across the Elegans supergroup that has been conserved for tens of millions of years. A model is proposed to account for the rapid origin of this essential kernel in the gut specification network, by the upstream intercalation of duplicate genes into a simpler ancestral network.

ChiIV3 Acts as a Novel Target of WRKY40 to Mediate Pepper Immunity Against Ralstonia solanacearum Infection

ChiIV3, a chitinase of pepper (Capsicum annuum), stimulates cell death in pepper plants. However, there are only scarce reports on its role in resistance against bacterial wilt disease such as that caused by Ralstonia solanacearum and their transcriptional regulation. In this study, the silencing of ChiIV3 in pepper plants significantly reduced the resistance to R. solanacearum. The transcript of ChiIV3 was induced by R. solanacearum inoculation (RSI) as well as exogenous application of methyl jasmonate and abscisic acid. The bioinformatics analysis revealed that the ChiIV3 promoter consists of multiple stress-related cis elements, including six W-boxes and one MYB1AT. With the 5' deletion assay in the ChiIV3 promoter, the W4-box located from -640 to -635 bp was identified as the cis element that is required for the response to RSI. In addition, the W4-box element was shown to be essential for the binding of the ChiIV3 promoter by the WRKY40 transcription factor, which is known to positively regulate the defense response to R. solanacearum. Site-directed mutagenesis in the W4-box sequence impaired the binding of WRKY40 to the ChiIV3 promoter. Subsequently, the transcription of ChiIV3 decreased in WRKY40-silenced pepper plants. These results demonstrated that the expression of the defense gene ChiIV3 is controlled through multiple modes of regulation, and WRKY40 directly binds to the W4-box element of the ChiIV3 promoter region for its transcriptional regulation.

A two-level model for the role of complex and young genes in the formation of organism complexity and new insights into the relationship between evolution and development

Background

How genome complexity affects organismal phenotypic complexity is a fundamental question in evolutionary developmental biology. Previous studies proposed various contributing factors of genome complexity and tried to find the connection between genomic complexity and organism complexity. However, a general model to answer this question is lacking. Here, we introduce a ‘two-level’ model for the realization of genome complexity at phenotypic level.

Results

Five representative species across Protostomia and Deuterostomia were involved in this study. The intrinsic gene properties contributing to genome complexity were classified into two generalized groups: the complexity and age degree of both protein-coding and noncoding genes. We found that young genes tend to be simpler; however, the mid-age genes, rather than the oldest genes, show the highest proportion of high complexity. Complex genes tend to be utilized preferentially in each stage of embryonic development, with maximum representation during the late stage of organogenesis. This trend is mainly attributed to mid-age complex genes. In contrast, young genes tend to be expressed in specific spatiotemporal states. An obvious correlation between the time point of the change in over- and under-representation and the order of gene age was observed, which supports the funnel-like model of the conservation pattern of development. In addition, we found some probable causes for the seemingly contradictory ‘funnel-like’ or ‘hourglass’ model.

Conclusions

These results indicate that complex and young genes contribute to organismal complexity at two different levels: Complex genes contribute to the complexity of individual proteomes in certain states, whereas young genes contribute to the diversity of proteomes in different spatiotemporal states. This conclusion is valid across the five species investigated, indicating it is a conserved model across Protostomia and Deuterostomia. The results in this study also support ‘funnel-like model’ from a new viewpoint and explain why there are different evo–devo relation models.

Electronic supplementary material

The online version of this article (10.1186/s13227-018-0111-4) contains supplementary material, which is available to authorized users.

The combinatorial control of alternative splicing in C. elegans

Normal development requires the right splice variants to be made in the right tissues at the right time. The core splicing machinery is engaged in all splicing events, but which precise splice variant is made requires the choice between alternative splice sites—for this to occur, a set of splicing factors (SFs) must recognize and bind to short RNA motifs in the pre-mRNA. In C. elegans, there is known to be extensive variation in splicing patterns across development, but little is known about the targets of each SF or how multiple SFs combine to regulate splicing. Here we combine RNA-seq with in vitro binding assays to study how 4 different C. elegans SFs, ASD-1, FOX-1, MEC-8, and EXC-7, regulate splicing. The 4 SFs chosen all have well-characterised biology and well-studied loss-of-function genetic alleles, and all contain RRM domains. Intriguingly, while the SFs we examined have varied roles in C. elegans development, they show an unexpectedly high overlap in their targets. We also find that binding sites for these SFs occur on the same pre-mRNAs more frequently than expected suggesting extensive combinatorial control of splicing. We confirm that regulation of splicing by multiple SFs is often combinatorial and show that this is functionally significant. We also find that SFs appear to combine to affect splicing in two modes—they either bind in close proximity within the same intron or they appear to bind to separate regions of the intron in a conserved order. Finally, we find that the genes whose splicing are regulated by multiple SFs are highly enriched for genes involved in the cytoskeleton and in ion channels that are key for neurotransmission. Together, this shows that specific classes of genes have complex combinatorial regulation of splicing and that this combinatorial regulation is critical for normal development to occur.

Alternative splicing (AS) is a highly regulated process that is crucial for normal development. It requires the core splicing machinery, but the specific choice of splice site during AS is controlled by splicing factors (SFs) such as ELAV or RBFOX proteins that bind to specific sequences in pre-mRNAs to regulate usage of different splice sites. AS varies across the C. elegans life cycle and here we study how diverse SFs combine to regulate AS during C. elegans development. We selected 4 RRM-containing SFs that are all well studied and that have well-characterised loss-of-function genetic alleles. We find that these SFs regulate many of the same targets, and that combinatorial interactions between these SFs affect both individual splicing events and organism-level phenotypes including specific effects on the neuromuscular system. We further show that SFs combine to regulate splicing of an individual pre-mRNA in two distinct modes—either by binding in close proximity or by binding in a defined order on the pre-mRNA. Finally, we find that the genes whose splicing are most likely to be regulated by multiple SFs are genes that are required for the proper function of the neuromuscular system. These genes are also most likely to have changing AS patterns across development, suggesting that their splicing regulation is highly complex and developmentally regulated. Taken together, our data show that the precise splice variant expressed at any point in development is often the outcome of regulation by multiple SFs.

Understanding Haemonchus contortus Better Through Genomics and Transcriptomics

Parasitic roundworms (nematodes) cause substantial mortality and morbidity in animals globally. The barber's pole worm, Haemonchus contortus, is one of the most economically significant parasitic nematodes of small ruminants worldwide. Although this and related nematodes can be controlled relatively well using anthelmintics, resistance against most drugs in common use has become a major problem. Until recently, almost nothing was known about the molecular biology of H. contortus on a global scale. This chapter gives a brief background on H. contortus and haemonchosis, immune responses, vaccine research, chemotherapeutics and current problems associated with drug resistance. It also describes progress in transcriptomics before the availability of H. contortus genomes and the challenges associated with such work. It then reviews major progress on the two draft genomes and developmental transcriptomes of H. contortus, and summarizes their implications for the molecular biology of this worm in both the free-living and the parasitic stages of its life cycle. The chapter concludes by considering how genomics and transcriptomics can accelerate research on Haemonchus and related parasites, and can enable the development of new interventions against haemonchosis.

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

A modular system of DNA enhancer elements mediates tissue-specific activation of transcription by high dietary zinc in C. elegans

Zinc is essential for biological systems, and aberrant zinc metabolism is implicated in a broad range of human diseases. To maintain homeostasis in response to fluctuating levels of dietary zinc, animals regulate gene expression; however, mechanisms that mediate the transcriptional response to fluctuating levels of zinc have not been fully defined. Here, we identified DNA enhancer elements that mediate intestine-specific transcriptional activation in response to high levels of dietary zinc in C. elegans. Using bioinformatics, we characterized an evolutionarily conserved enhancer element present in multiple zinc-inducible genes, the high zinc activation (HZA) element. The HZA was consistently adjacent to a GATA element that mediates expression in intestinal cells. Functional studies using transgenic animals demonstrated that this modular system of DNA enhancers mediates tissue-specific transcriptional activation in response to high levels of dietary zinc. We used this information to search the genome and successfully identified novel zinc-inducible genes. To characterize the mechanism of enhancer function, we demonstrated that the GATA transcription factor ELT-2 and the mediator subunit MDT-15 are necessary for zinc-responsive transcriptional activation. These findings define new mechanisms of zinc homeostasis and tissue-specific regulation of transcription.

Collaborative regulation of development but independent control of metabolism by two epidermis-specific transcription factors in Caenorhabditis elegans

Cell fate specification is typically initiated by a master regulator, which is relayed by tissue-specific regulatory proteins (usually transcription factors) for further enforcement of cell identities, but how the factors are coordinated among each other to "finish up" the specification remains poorly understood. Caenorhabditis elegans epidermis specification is initiated by a master regulator, ELT-1, that activates its targets, NHR-25 and ELT-3, two epidermis-specific transcription factors that are important for development but not for initial specification of epidermis, thus providing a unique paradigm for illustrating how the tissue-specific regulatory proteins work together to enforce cell fate specification. Here we addressed the question through contrasting genome-wide in vivo binding targets between NHR-25 and ELT-3. We demonstrate that the two factors bind discrete but conserved DNA motifs, most of which remain in proximity, suggesting formation of a complex between the two. In agreement with this, gene ontology analysis of putative target genes suggested differential regulation of metabolism but coordinated control of epidermal development between the two factors, which is supported by quantitative analysis of expression of their specific or common targets in the presence or absence of either protein. Functional validation of a subset of the target genes showed both activating and inhibitory roles of NHR-25 and ELT-3 in regulating their targets. We further demonstrated differential control of specification of AB and C lineage-derived epidermis. The results allow us to assemble a comprehensive gene network underlying C. elegans epidermis development that is likely to be widely used across species and provides insights into how tissue-specific transcription factors coordinate with one another to enforce cell fate specification initiated by its master regulator.

Transcription factor binding to Caenorhabditis elegans first introns reveals lack of redundancy with gene promoters

Gene expression is controlled through the binding of transcription factors (TFs) to regulatory genomic regions. First introns are longer than other introns in multiple eukaryotic species and are under selective constraint. Here we explore the importance of first introns in TF binding in the nematode Caenorhabditis elegans by combining computational predictions and experimentally derived TF–DNA interaction data. We found that first introns of C. elegans genes, particularly those for families enriched in long first introns, are more conserved in length, have more conserved predicted TF interactions and are bound by more TFs than other introns. We detected a significant positive correlation between first intron size and the number of TF interactions obtained from chromatin immunoprecipitation assays or determined by yeast one-hybrid assays. TFs that bind first introns are largely different from those binding promoters, suggesting that the different interactions are complementary rather than redundant. By combining first intron and promoter interactions, we found that genes that share a large fraction of TF interactions are more likely to be co-expressed than when only TF interactions with promoters are considered. Altogether, our data suggest that C. elegans gene regulation may be additive through the combined effects of multiple regulatory regions.

A negative-feedback loop between the detoxification/antioxidant response factor SKN-1 and its repressor WDR-23 matches organism needs with environmental conditions

Negative-feedback loops between transcription factors and repressors in responses to xenobiotics, oxidants, heat, hypoxia, DNA damage, and infection have been described. Although common, the function of feedback is largely unstudied. Here, we define a negative-feedback loop between the Caenorhabditis elegans detoxification/antioxidant response factor SKN-1/Nrf and its repressor wdr-23 and investigate its function in vivo. Although SKN-1 promotes stress resistance and longevity, we find that tight regulation by WDR-23 is essential for growth and reproduction. By disabling SKN-1 transactivation of wdr-23, we reveal that feedback is required to set the balance between growth/reproduction and stress resistance/longevity. We also find that feedback is required to set the sensitivity of a core SKN-1 target gene to an electrophile. Interestingly, the effect of feedback on target gene induction is greatly reduced when the stress response is strongly activated, presumably to ensure maximum activation of cytoprotective genes during potentially fatal conditions. Our work provides a framework for understanding the function of negative feedback in inducible stress responses and demonstrates that manipulation of feedback alone can shift the balance of competing animal processes toward cell protection, health, and longevity.

Genome-wide prediction and functional validation of promoter motifs regulating gene expression in spore and infection stages of Phytophthora infestans

Most eukaryotic pathogens have complex life cycles in which gene expression networks orchestrate the formation of cells specialized for dissemination or host colonization. In the oomycete Phytophthora infestans, the potato late blight pathogen, major shifts in mRNA profiles during developmental transitions were identified using microarrays. We used those data with search algorithms to discover about 100 motifs that are over-represented in promoters of genes up-regulated in hyphae, sporangia, sporangia undergoing zoosporogenesis, swimming zoospores, or germinated cysts forming appressoria (infection structures). Most of the putative stage-specific transcription factor binding sites (TFBSs) thus identified had features typical of TFBSs such as position or orientation bias, palindromy, and conservation in related species. Each of six motifs tested in P. infestans transformants using the GUS reporter gene conferred the expected stage-specific expression pattern, and several were shown to bind nuclear proteins in gel-shift assays. Motifs linked to the appressoria-forming stage, including a functionally validated TFBS, were over-represented in promoters of genes encoding effectors and other pathogenesis-related proteins. To understand how promoter and genome architecture influence expression, we also mapped transcription patterns to the P. infestans genome assembly. Adjacent genes were not typically induced in the same stage, including genes transcribed in opposite directions from small intergenic regions, but co-regulated gene pairs occurred more than expected by random chance. These data help illuminate the processes regulating development and pathogenesis, and will enable future attempts to purify the cognate transcription factors.

The genus Phytophthora includes over one hundred species of plant pathogens that have devastating effects worldwide in agriculture and natural environments. Its most notorious member is P. infestans, which causes the late blight diseases of potato and tomato. Their success as pathogens is dependent on the formation of specialized cells for plant-to-plant transmission and host infection, but little is known about how this is regulated. Recognizing that changes in gene expression drive the formation of these cell types, we used a computational approach to predict the sequences of about one hundred transcription factor binding sites associated with expression in either of five life stages, including several types of spores and infection structures. We then used a functional testing strategy to prove their biological activity by showing that the DNA motifs enabled the stage-specific expression of a transgene. Our work lays the groundwork for dissecting the molecular mechanisms that regulate life-stage transitions and pathogenesis in Phytophthora. A similar approach should be useful for other plant and animal pathogens.

Computational prediction of the polyQ and CAG repeat spinocerebellar ataxia network based on sequence identity to untranslated regions

Computational prediction of biological networks would be a tremendous asset to systems biology and personalized medicine. In this paper, we use a moving window bioinformatic screen to identify transcripts with partial identity to the 5' and 3'UTRs of the polyQ spinocerebellar ataxia (SCA) genes ATXN1, ATXN2, ATXN3, ATXN7, TBP and CACNA1A and the CAG repeat expansion gene PPP2R2B. We find that the bioinformatic screen enriches for transcripts that encode proteins that interact and that have functions relevant to polyQ SCA. Transcription control and RNA binding are the primary functional groups represented in the proteins from the combined screens. The insulin growth factor pathway, the WNT pathway, long term potentiation, melanogenesis and ATM mediated DNA repair pathways were identified as important pathways. UGUUU repeats were identified as an abundant motif in the SCA network and PAXIP1, CELF2, CREBBP, EBF1, PLEKHG4, SRSF4, C5orf42, NFIA, STK24, and YWHAG were identified as statistically significant proteins in the polyQ and PPP2R2B network.