Aryl hydrocarbon receptor (AHR) in the cnidarian Nematostella vectensis: comparative expression, protein interactions, and ligand binding

Persistent organic pollutants (POPs) in marine crustaceans: Bioaccumulation, physiological and cellular responses

Persistent organic pollutants (POPs) are ubiquitous in marine ecosystems. These compounds can be accumulated in water, sediments and organisms, persist in time, and have toxic effects in human and wildlife. POPs can be uptaken and bioaccumulated by crustaceans, affecting different physiological processes, including energy metabolism, immunity, osmoregulation, excretion, growth, and reproduction. Nonetheless, animals have evolved sub-cellular mechanisms for detoxification and protection from chemical stress. POPs induce the activity of enzymes involved in xenobiotic metabolism and antioxidant systems, that in vertebrates are importantly regulated at gene expression (transcriptional) level. However, the activation and control of these enzyme systems upon the exposure to POPs have been scarcely studied in invertebrate species, including crustaceans. Herein, we summarize various aspects of the bioaccumulation of POPs in marine crustaceans and their physiological effects. We specially focus on the regulation of xenobiotics metabolism and antioxidant enzymes as key sub-cellular mechanisms for detoxification and protection from chemical stress.

Molecular Evolution of Aryl Hydrocarbon Receptor Signaling Pathway Genes

The Aryl hydrocarbon receptor is an ancient transcriptional factor originally discovered as a sensor of dioxin. In addition to its function as a receptor of environmental toxicants, it plays an important role in development. Although a significant amount of research has been carried out to understand the AHR signal transduction pathway and its involvement in species' susceptibility to environmental toxicants, none of them to date has comprehensively studied its evolutionary origins. Studying the evolutionary origins of molecules can inform ancestral relationships of genes. The vertebrate genome has been shaped by two rounds of whole-genome duplications (WGD) at the base of vertebrate evolution approximately 600 million years ago, followed by lineage-specific gene losses, which often complicate the assignment of orthology. It is crucial to understand the evolutionary origins of this transcription factor and its partners, to distinguish orthologs from ancient non-orthologous homologs. In this study, we have investigated the evolutionary origins of proteins involved in the AHR pathway. Our results provide evidence of gene loss and duplications, crucial for understanding the functional connectivity of humans and model species. Multiple studies have shown that 2R-ohnologs (genes and proteins that have survived from the 2R-WGD) are enriched in signaling components relevant to developmental disorders and cancer. Our findings provide a link between the AHR pathway's evolutionary trajectory and its potential mechanistic involvement in pathogenesis.

Step in Time: Conservation of Circadian Clock Genes in Animal Evolution

Over the past few decades, the molecular mechanisms responsible for circadian phenotypes of animals have been studied in increasing detail in mammals, some insects, and other invertebrates. Particular circadian proteins and their interactions are shared across evolutionary distant animals, resulting in a hypothesis for the canonical circadian clock of animals. As the number of species for which the circadian clockwork has been described increases, the circadian clock in animals driving cyclical phenotypes becomes less similar. Our focus in this review is to develop and synthesize the current literature to better understand the antiquity and evolution of the animal circadian clockwork. Here, we provide an updated understanding of circadian clock evolution in animals, largely through the lens of conserved genes characterized in the circadian clock identified in bilaterian species. These comparisons reveal extensive variation within the likely composition of the core clock mechanism, including losses of many genes, and that the ancestral clock of animals does not equate to the bilaterian clock. Despite the loss of these core genes, these species retain circadian behaviors and physiology, suggesting novel clocks have evolved repeatedly. Additionally, we highlight highly conserved cellular processes (e.g., cell division, nutrition) that intersect with the circadian clock of some animals. The conservation of these processes throughout the animal tree remains essentially unknown, but understanding their role in the evolution and maintenance of the circadian clock will provide important areas for future study.

Identification of Modulators of the C. elegans Aryl Hydrocarbon Receptor and Characterization of Transcriptomic and Metabolic AhR-1 Profiles

The Aryl hydrocarbon Receptor (AhR) is a xenobiotic sensor in vertebrates, regulating the metabolism of its own ligands. However, no ligand has been identified to date for any AhR in invertebrates. In C. elegans, the AhR ortholog, AHR-1, displays physiological functions. Therefore, we compared the transcriptomic and metabolic profiles of worms expressing AHR-1 or not and investigated the putative panel of chemical AHR-1 modulators. The metabolomic profiling indicated a role for AHR-1 in amino acids, carbohydrates, and fatty acids metabolism. The transcriptional profiling in neurons expressing AHR-1, identified 95 down-regulated genes and 76 up-regulated genes associated with neuronal and metabolic functions in the nervous system. A gene reporter system allowed us to identify several AHR-1 modulators including bacterial, dietary, or environmental compounds. These results shed new light on the biological functions of AHR-1 in C. elegans and perspectives on the evolution of the AhR functions across species.

The Ah Receptor: Adaptive Metabolism, Ligand Diversity, and the Xenokine Model

The Ah receptor (AHR) has been studied for almost five decades. Yet, we still have many important questions about its role in normal physiology and development. Moreover, we still do not fully understand how this protein mediates the adverse effects of a variety of environmental pollutants, such as the polycyclic aromatic hydrocarbons (PAHs), the chlorinated dibenzo-p-dioxins ("dioxins"), and many polyhalogenated biphenyls. To provide a platform for future research, we provide the historical underpinnings of our current state of knowledge about AHR signal transduction, identify a few areas of needed research, and then develop concepts such as adaptive metabolism, ligand structural diversity, and the importance of proligands in receptor activation. We finish with a discussion of the cognate physiological role of the AHR, our perspective on why this receptor is so highly conserved, and how we might think about its cognate ligands in the future.

The molecular mechanism of AhR-ARNT-XREs signaling pathway in the detoxification response induced by polycyclic aromatic hydrocarbons (PAHs) in clam Ruditapes philippinarum

The aryl hydrocarbon receptor (AhR) has been known primarily for its role in the regulation of several drug and xenobiotic metabolizing enzymes to mitigate environmental stresses. In this study, we interfere the expression of AhR gene to investigate the mechanism of AhR signaling pathway in the detoxification and antioxidation defense system that induced by Polycyclic Aromatic Hydrocarbons (PAHs) exposure by RNA interference (RNAi). The gene expressions of aryl hydrocarbon receptor nuclear translocator (ARNT), heat shock protein 90 (Hsp90) were evaluated after being exposed to benzo(a)pyrene (BaP) (4 μg/L) for 5 days and the positive correlations between AhR, ARNT, HSP90 indirectly indicating that AhR may have the ability to bind to ligands such as PAHs in Ruditapes philippinarum (R. philippinarum). Besides, the activities of detoxification enzymes were determined to investigate the role of AhR signaling pathway played in the metabolic detoxification. What's more, the gene expressions of protein kinase C (PKC) signaling pathway, mitogen-activated protein kinase (MAPKs) signaling pathway, NF-E2-related factor 2 (Nrf2) signaling pathway and antioxidant defense system indicated that AhR may regulate the Nrf2-Keap1 signaling pathway through Kelch-like ECH-associated protein-1 (Keap1) and MAPKs, PKC signaling pathways. In conclusion, adoption of RNA interference technology to explore the role of RpAhR gene played in the detoxification and antioxidation defense system under the PAHs stress at different time points can informe molecular endpoints for application towards ecotoxicology monitoring of bivalves.

Transcriptional remodelling upon light removal in a model cnidarian: Losses and gains in gene expression

Organismal responses to light:dark cycles can result from two general processes: (a) direct response to light or (b) a free-running rhythm (i.e., a circadian clock). Previous research in cnidarians has shown that candidate circadian clock genes have rhythmic expression in the presence of diel lighting, but these oscillations appear to be lost quickly after removal of the light cue. Here, we measure whole-organism gene expression changes in 136 transcriptomes of the sea anemone Nematostella vectensis, entrained to a light:dark environment and immediately following light cue removal to distinguish two broadly defined responses in cnidarians: light entrainment and circadian regulation. Direct light exposure resulted in significant differences in expression for hundreds of genes, including more than 200 genes with rhythmic, 24-hr periodicity. Removal of the lighting cue resulted in the loss of significant expression for 80% of these genes after 1 day, including most of the hypothesized cnidarian circadian genes. Further, 70% of these candidate genes were phase-shifted. Most surprisingly, thousands of genes, some of which are involved in oxidative stress, DNA damage response and chromatin modification, had significant differences in expression in the 24 hr following light removal, suggesting that loss of the entraining cue may induce a cellular stress response. Together, our findings suggest that a majority of genes with significant differences in expression for anemones cultured under diel lighting are largely driven by the primary photoresponse rather than a circadian clock when measured at the whole animal level. These results provide context for the evolution of cnidarian circadian biology and help to disassociate two commonly confounded factors driving oscillating phenotypes.

Dynamic remodeling of the interactomes of Nematostella vectensis Hsp70 isoforms under heat shock

Heat shock protein 70s (Hsp70s) are a highly conserved class of molecular chaperones that fold a large proportion of the proteome. Nematostella vectensis (Nv) is an estuarine sea anemone that has emerged as a model species to characterize molecular responses to physiological stressors due to its exposure to diverse, extreme abiotic conditions. Previous transcriptional data has shown dramatic differences among expression profiles of three NvHsp70 isoforms (NvHsp70A, B and D) under stress but it is unknown if, and to what extent, the client proteins for these chaperones differ. In order to determine client specificity, NvHsp70A, B and D were expressed in Saccharomyces cerevisiae budding yeast lacking native Hsp70 and interacting proteins for each Hsp70 were determined with mass spectrometry in yeast ambient and heat shock conditions. Our analyses showed <50% of identified interacting proteins were common to all three anemone Hsp70s and 3-18% were unique to an individual Hsp70. Mapping of temperature induced interactions suggest that under stress a proportion of clients are transferred from NvHsp70A and NvHsp70D to NvHsp70B. Together, these data suggest a diverse set of interacting proteins for Hsp70 isoforms that likely determines the precise functions for Hsp70s in organismal acclimation and potentially adaptation. BIOLOGICAL SIGNIFICANCE: Although the Hsp70 family of molecular chaperones has been studied for >50 years, it is still not fully understood why organisms encode and express many highly-similar Hsp70 isoforms. The prevailing theory is that these isoforms have identical function, but are expressed under unique cellular conditions that include heat shock to cope with increased number of unfolded/misfolded proteins. The sea anemone Nematostella vectensis encodes three Hsp70 isoforms A, B and D that when expressed in yeast demonstrate unique functionalities. This study provides the interactome of NvHsp70s A, B and D and demonstrates that Hsp70 isoforms, while highly similar in sequence, have unique co-chaperone and client interactors.

Ultraviolet radiation significantly enhances the molecular response to dispersant and sweet crude oil exposure in Nematostella vectensis

Estuarine organisms are subjected to combinations of anthropogenic and natural stressors, which together can reduce an organisms' ability to respond to either stress or can potentiate or synergize the cellular impacts for individual stressors. Nematostella vectensis (starlet sea anemone) is a useful model for investigating novel and evolutionarily conserved cellular and molecular responses to environmental stress. Using RNA-seq, we assessed global changes in gene expression in Nematostella in response to dispersant and/or sweet crude oil exposure alone or combined with ultraviolet radiation (UV). A total of 110 transcripts were differentially expressed by dispersant and/or crude oil exposure, primarily dominated by the down-regulation of 74 unique transcripts in the dispersant treatment. In contrast, UV exposure alone or combined with dispersant and/or oil resulted in the differential expression of 1133 transcripts, of which 436 were shared between all four treatment combinations. Most significant was the differential expression of 531 transcripts unique to one or more of the combined UV/chemical exposures. Main categories of genes affected by one or more of the treatments included enzymes involved in xenobiotic metabolism and transport, DNA repair enzymes, and general stress response genes conserved among vertebrates and invertebrates. However, the most interesting observation was the induction of several transcripts indicating de novo synthesis of mycosporine-like amino acids and other novel cellular antioxidants. Together, our data suggest that the toxicity of oil and/or dispersant and the complexity of the molecular response are significantly enhanced by UV exposure, which may co-occur for shallow water species like Nematostella.

Aryl hydrocarbon receptor (AHR): "pioneer member" of the basic-helix/loop/helix per-Arnt-sim (bHLH/PAS) family of "sensors" of foreign and endogenous signals

The basic-helix/loop/helix per-Arnt-sim (bHLH/PAS) family comprises many transcription factors, found throughout all three kingdoms of life; bHLH/PAS members "sense" innumerable intracellular and extracellular "signals" - including endogenous compounds, foreign chemicals, gas molecules, redox potential, photons (light), gravity, heat, and osmotic pressure. These signals then initiate downstream signaling pathways involved in responding to that signal. The term "PAS", abbreviation for "per-Arnt-sim" was first coined in 1991. Although the mouse Arnt gene was not identified until 1991, evidence of its co-transcriptional binding partner, aryl hydrocarbon receptor (AHR), was first reported in 1974 as a "sensor" of foreign chemicals, up-regulating cytochrome P450 family 1 (CYP1) and other enzyme activities that usually metabolize the signaling chemical. Within a few years, AHR was proposed also to participate in inflammation. The mouse [Ah] locus was shown (1973-1989) to be relevant to chemical carcinogenesis, mutagenesis, toxicity and teratogenesis, the mouse Ahr gene was cloned in 1992, and the first Ahr(-/-) knockout mouse line was reported in 1995. After thousands of studies from the early 1970s to present day, we now realize that AHR participates in dozens of signaling pathways involved in critical-life processes, affecting virtually every organ and cell-type in the animal, including many invertebrates.

Applying evolutionary genetics to developmental toxicology and risk assessment

Evolutionary thinking continues to challenge our views on health and disease. Yet, there is a communication gap between evolutionary biologists and toxicologists in recognizing the connections among developmental pathways, high-throughput screening, and birth defects in humans. To increase our capability in identifying potential developmental toxicants in humans, we propose to apply evolutionary genetics to improve the experimental design and data interpretation with various in vitro and whole-organism models. We review five molecular systems of stress response and update 18 consensual cell-cell signaling pathways that are the hallmark for early development, organogenesis, and differentiation; and revisit the principles of teratology in light of recent advances in high-throughput screening, big data techniques, and systems toxicology. Multiscale systems modeling plays an integral role in the evolutionary approach to cross-species extrapolation. Phylogenetic analysis and comparative bioinformatics are both valuable tools in identifying and validating the molecular initiating events that account for adverse developmental outcomes in humans. The discordance of susceptibility between test species and humans (ontogeny) reflects their differences in evolutionary history (phylogeny). This synthesis not only can lead to novel applications in developmental toxicity and risk assessment, but also can pave the way for applying an evo-devo perspective to the study of developmental origins of health and disease.

Human and rodent aryl hydrocarbon receptor (AHR): from mediator of dioxin toxicity to physiologic AHR functions and therapeutic options

Metabolism of aryl hydrocarbons and toxicity of dioxins led to the discovery of the aryl hydrocarbon receptor (AHR). Tremendous advances have been made on multiplicity of AHR signaling and identification of endogenous ligands including the tryptophan metabolites FICZ and kynurenine. However, human AHR functions are still poorly understood due to marked species differences as well as cell-type- and cell context-dependent AHR functions. Observations in dioxin-poisoned individuals may provide hints to physiologic AHR functions in humans. Based on these observations three human AHR functions are discussed: (1) Chemical defence and homeostasis of endobiotics. The AHR variant Val381 in modern humans leads to reduced AHR affinity to aryl hydrocarbons in comparison with Neanderthals and primates expressing the Ala381 variant while affinity to indoles remains unimpaired. (2) Homeostasis of stem/progenitor cells. Dioxins dysregulate homeostasis in sebocyte stem cells. (3) Modulation of immunity. In addition to microbial defence, AHR may be involved in a ‘disease tolerance defence pathway’. Further characterization of physiologic AHR functions may lead to therapeutic options.

Assessing the identity and expression level of the cytochrome P450 20A1 (CYP20A1) gene in the BPA-, BDE-47, and WAF-exposed copepods Tigriopus japonicus and Paracyclopina nana

CYP20A1 is a member of the cytochrome P450 (CYP) superfamily, identified as an orphan P450 without any assigned biological function; hence, its continued status as an "orphan" gene. In order to address this shortcoming in our understanding of this superfamily, we sought to characterize the CYP20A1 gene in the copepods Tigriopus japonicus (Tj-CYP20A1) and Paracyclopina nana (Pn-CYP20A1) at their mRNA transcriptional level. We assessed the response of this gene's expression in various developmental stages and in response to treatment with bisphenol A (BPA), 2, 2', 4, 4'-tetrabromodiphenyl ether (BDE-47), and water accommodated fractions (WAFs) of crude oil. As shown in the vertebrate CYP20A1, both Tj-CYP20A1 and Pn-CYP20A1 contained characteristic conserved motifs and domain regions (I helix, K helix and heme-binding motifs) with unusual amino acid sequences apparent in their gene structure. Also molecular characterization of the putative responsive elements in the promoter regions was performed. We observed transcriptional up-regulation of these genes during post-embryonic developmental stages including sex-specific up-regulation in adults. In addition, concentration- and time-dependent mRNA transcripts in response to xenobiotics (BPA, BDE-47, and WAFs) were seen. This study focuses on the molecular elucidation of CYP20A1 genes and their interactions with xenobiotics in the copepods T. japonicus and P. nana that provides important insight into the biological importance of CYP20A1 in invertebrates.

Diversity as Opportunity: Insights from 600 Million Years of AHR Evolution

The aryl hydrocarbon receptor (AHR) was for many years of interest only to pharmacologists and toxicologists. However, this protein has fundamental roles in biology that are being revealed through studies in diverse animal species. The AHR is an ancient protein. AHR homologs exist in most major groups of modern bilaterian animals, including deuterostomes (chordates, hemichordates, echinoderms) and the two major clades of protostome invertebrates [ecdysozoans (e.g. arthropods and nematodes) and lophotrochozoans (e.g. molluscs and annelids)]. AHR homologs also have been identified in cnidarians such as the sea anemone Nematostella and in the genome of Trichoplax, a placozoan. Bilaterians, cnidarians, and placozoans form the clade Eumetazoa, whose last common ancestor lived approximately 600 million years ago (MYA). The presence of AHR homologs in modern representatives of all these groups indicates that the original eumetazoan animal possessed an AHR homolog. Studies in invertebrates and vertebrates reveal parallel functions of AHR in the development and function of sensory neural systems, suggesting that these may be ancestral roles. Vertebrate animals are characterized by the expansion and diversification of AHRs, via gene and genome duplications, from the ancestral protoAHR into at least five classes of AHR-like proteins: AHR, AHR1, AHR2, AHR3, and AHRR. The evolution of multiple AHRs in vertebrates coincided with the acquisition of high-affinity binding of halogenated and polynuclear aromatic hydrocarbons and the emergence of adaptive functions involving regulation of xenobiotic-metabolizing enzymes and roles in adaptive immunity. The existence of multiple AHRs may have facilitated subfunction partitioning and specialization of specific AHR types in some taxa. Additional research in diverse model and non-model species will continue to enrich our understanding of AHR and its pleiotropic roles in biology and toxicology.

Transcriptomic evaluation of the American oyster, Crassostrea virginica, deployed during the Deepwater Horizon oil spill: Evidence of an active hydrocarbon response pathway

Estuarine organisms were impacted by the Deepwater Horizon oil spill which released ∼5 million barrels of crude oil into the Gulf of Mexico in the spring and summer of 2010. Crassostrea virginica, the American oyster, is a keystone species in these coastal estuaries and is routinely used for environmental monitoring purposes. However, very little is known about their cellular and molecular responses to hydrocarbon exposure. In response to the spill, a monitoring program was initiated by deploying hatchery-reared oysters at three sites along the Alabama and Mississippi coast (Grand Bay, MS, Fort Morgan, AL, and Orange Beach, AL). Oysters were deployed for 2-month periods at five different time points from May 2010 to May 2011. Gill and digestive gland tissues were harvested for gene expression analysis and determination of aliphatic and polycyclic aromatic hydrocarbon (PAH) concentrations. To facilitate identification of stress response genes that may be involved in the hydrocarbon response, a nearly complete transcriptome was assembled using Roche 454 and Illumina high-throughput sequencing from RNA samples obtained from the gill and digestive gland tissues of deployed oysters. This effort resulted in the assembly and annotation of 27,227 transcripts comprised of a large assortment of stress response genes, including members of the aryl hydrocarbon receptor (AHR) pathway, Phase I and II biotransformation enzymes, antioxidant enzymes and xenobiotic transporters. From this assembly several potential biomarkers of hydrocarbon exposure were chosen for expression profiling, including the AHR, two cytochrome P450 1A genes (CYP1A-like 1 and CYP1A-like 2), Cu/Zn superoxide dismutase (CuZnSOD), glutathione S-transferase theta (GST theta) and multidrug resistance protein 3 (MRP3). Higher expression levels of GST theta and MRP3 were observed in gill tissues from all three sites during the summer to early fall 2010 deployments. Linear regression analysis indicated a statistically significant relationship between total PAH levels in digestive gland tissue samples with CYP1A-like 2, CuZnSOD, GST theta and MRP3 induction. These observations provide evidence of a potentially conserved AHR pathway in invertebrates and yield new insight into the development of novel biomarkers for use in environmental monitoring activities.

Dose- and time-dependent expression of aryl hydrocarbon receptor (AhR) and aryl hydrocarbon receptor nuclear translocator (ARNT) in PCB-, B[a]P-, and TBT-exposed intertidal copepod Tigriopus japonicus

The aryl hydrocarbon receptor (AhR) and aryl hydrocarbon nuclear translocator (ARNT) genes from the copepod Tigriopus japonicus (Tj) were cloned to examine their potential functions in the invertebrate putative AhR-CYP signaling pathway. The amino acid sequences encoded by the Tj-AhR and Tj-ARNT genes showed high similarity to homologs of Daphnia and Drosophila, ranging from 68% and 70% similarity for the AhR genes to 56% for the ARNT genes. To determine whether Tj-AhR and Tj-ARNT are modulated by environmental pollutants, transcriptional expression of Tj-AhR and Tj-ARNT was analyzed in response to exposure to five concentrations of polychlorinated biphenyl (PCB 126) (control, 10, 50, 100, 500 μg L(-1)), benzo[a]pyrene (B[a]P) (control, 5, 10, 50, 100 μg L(-1)), and tributyltin (TBT) (control, 1, 5, 10, 20 μg L(-1)) 24h after exposure. A time-course experiment (0, 3, 6, 12, 24h) was performed to analyze mRNA expression patterns after exposure to PCB, B[a]P, and TBT. T. japonicus exhibited dose-dependent and time-dependent upregulation of Tj-AhR and Tj-ARNT in response to pollutant exposure, and the degree of expression was dependent on the pollutant, suggesting that pollutants such as PCB, B[a]P, and TBT modulate expression of Tj-AhR and Tj-ARNT genes in the putative AhR-CYP signaling pathway.

Early and late response of Nematostella vectensis transcriptome to heavy metals

Environmental contamination from heavy metals poses a global concern for the marine environment, as heavy metals are passed up the food chain and persist in the environment long after the pollution source is contained. Cnidarians play an important role in shaping marine ecosystems, but environmental pollution profoundly affects their vitality. Among the cnidarians, the sea anemone Nematostella vectensis is an advantageous model for addressing questions in molecular ecology and toxicology as it tolerates extreme environments and its genome has been published. Here, we employed a transcriptome-wide RNA-Seq approach to analyse N. vectensis molecular defence mechanisms against four heavy metals: Hg, Cu, Cd and Zn. Altogether, more than 4800 transcripts showed significant changes in gene expression. Hg had the greatest impact on up-regulating transcripts, followed by Cu, Zn and Cd. We identified, for the first time in Cnidaria, co-up-regulation of immediate-early transcription factors such as Egr1, AP1 and NF-κB. Time-course analysis of these genes revealed their early expression as rapidly as one hour after exposure to heavy metals, suggesting that they may complement or substitute for the roles of the metal-mediating Mtf1 transcription factor. We further characterized the regulation of a large array of stress-response gene families, including Hsp, ABC, CYP members and phytochelatin synthase, that may regulate synthesis of the metal-binding phytochelatins instead of the metallothioneins that are absent from Cnidaria genome. This study provides mechanistic insight into heavy metal toxicity in N. vectensis and sheds light on ancestral stress adaptations.

Developmental and light-entrained expression of melatonin and its relationship to the circadian clock in the sea anemone Nematostella vectensis

Background

The primary hormone of the vertebrate pineal gland, melatonin, has been identified broadly throughout the eukaryotes. While the role for melatonin in cyclic behavior via interactions with the circadian clock has only been reported in vertebrates, comparative research has shown that the transcription-translation loops of the animal circadian clock likely date to the cnidarian-bilaterian ancestor, leaving open significant questions about the evolutionary origin of melatonin signaling in circadian behavior by interacting with the molecular clock.

Results

Expression of melatonin in adult anemones showed peak expression at the end of light period (zeitgeber time (ZT) = 12) when cultured under diel conditions, coinciding with expression of genes and enzyme activity for members of the melatonin synthesis pathway (tryptophan hydroxylase and hydroxyindol-O-methyltransferase), which also showed rhythmic expression. During embryogenesis and juvenile stages, melatonin showed cyclic oscillations in concentration, peaking in midday. Spatial (in situ hybridization) and quantitative (real-time PCR) transcription of clock genes during development of N. vectensis showed these ‘clock’ genes are expressed early in the development, prior to rhythmic oscillations, suggesting functions independent of a function in the circadian clock. Finally, time-course studies revealed that animals transferred from diel conditions to constant darkness lose circadian expression for most of the clock genes within 4 days, which can be reset by melatonin supplementation.

Conclusions

Our results support an ancient role for melatonin in the circadian behavior of animals by showing cyclic expression of this hormone under diel conditions, light-dependent oscillations in genes in the melatonin synthesis pathway, and the function of melatonin in initiating expression of circadian clock genes in the cnidarian N. vectensis. The differences in expression melatonin and the circadian clock gene network in the adult stage when compared with developmental stages of N. vectensis suggests new research directions to characterize stage-specific mechanisms of circadian clock function in animals.