Characterization of large numbers of single nucleotide polymorphisms (SNPs) throughout a genome has the power to refine the understanding of population demographic history and to identify genomic regions under selection in natural populations. To this end, population genomic approaches that harness the power of next-generation sequencing to understand the ecology and evolution of marine invertebrates represent a boon to test long-standing questions in marine biology and conservation. We employed restriction-site-associated DNA sequencing (RAD-seq) to identify SNPs in natural populations of the sea anemone Nematostella vectensis, an emerging cnidarian model with a broad geographic range in estuarine habitats in North and South America, and portions of England. We identified hundreds of SNP-containing tags in thousands of RAD loci from 30 barcoded individuals inhabiting four locations from Nova Scotia to South Carolina. Population genomic analyses using high-confidence SNPs resulted in a highly-resolved phylogeography, a result not achieved in previous studies using traditional markers. Plots of locus-specific FST against heterozygosity suggest that a majority of polymorphic sites are neutral, with a smaller proportion suggesting evidence for balancing selection. Loci inferred to be under balancing selection were mapped to the genome, where 90% were located in gene bodies, indicating potential targets of selection. Results from analyses with and without a reference genome supported similar conclusions, further supporting RAD-seq as a method that can be efficiently applied to species lacking existing genomic resources. We discuss the utility of RAD-seq approaches in burgeoning Nematostella research as well as in other cnidarian species, particularly corals, to determine phylogeographic relationships of populations and identify regions of the genome undergoing selection.

BackgroundCircadian rhythms in behavior and physiology are the observable phenotypes from cycles in expression of, interactions between, and degradation of the underlying molecular components. In bilaterian animals, the core molecular components include Timeless-Timeout, photoreceptive cryptochromes, and several members of the basic-loop-helix-Per-ARNT-Sim (bHLH-PAS) family. While many of core circadian genes are conserved throughout the Bilateria, their specific roles vary among species. Here, we identify and experimentally study the rhythmic gene expression of conserved circadian clock members in a sea anemone in order to characterize this gene network in a member of the phylum Cnidaria and to infer critical components of the clockwork used in the last common ancestor of cnidarians and bilaterians.Methodology/Principal FindingsWe identified homologs of circadian regulatory genes in the sea anemone Nematostella vectensis, including a gene most similar to Timeout, three cryptochromes, and several key bHLH-PAS transcription factors. We then maintained N. vectensis either in complete darkness or in a 12 hour light: 12 hour dark cycle in three different light treatments (blue only, full spectrum, blue-depleted). Gene expression varied in response to light cycle and light treatment, with a particularly strong pattern observed for NvClock. The cryptochromes more closely related to the light-sensitive clade of cryptochromes were upregulated in light treatments that included blue wavelengths. With co-immunoprecipitation, we determined that heterodimerization between CLOCK and CYCLE is conserved within N. vectensis. Additionally, we identified E-box motifs, DNA sequences recognized by the CLOCK:CYCLE heterodimer, upstream of genes showing rhythmic expression.Conclusions/SignificanceThis study reveals conserved molecular and functional components of the circadian clock that were in place at the divergence of the Cnidaria and Bilateria, suggesting the animal circadian clockwork is more ancient than previous data suggest. Characterizing circadian regulation in a cnidarian provides insight into the early origins of animal circadian rhythms and molecular regulation of environmentally cued behaviors.

BackgroundCalcium-activated photoproteins are luciferase variants found in photocyte cells of bioluminescent jellyfish (Phylum Cnidaria) and comb jellies (Phylum Ctenophora). The complete genomic sequence from the ctenophore Mnemiopsis leidyi, a representative of the earliest branch of animals that emit light, provided an opportunity to examine the genome of an organism that uses this class of luciferase for bioluminescence and to look for genes involved in light reception. To determine when photoprotein genes first arose, we examined the genomic sequence from other early-branching taxa. We combined our genomic survey with gene trees, developmental expression patterns, and functional protein assays of photoproteins and opsins to provide a comprehensive view of light production and light reception in Mnemiopsis.ResultsThe Mnemiopsis genome has 10 full-length photoprotein genes situated within two genomic clusters with high sequence conservation that are maintained due to strong purifying selection and concerted evolution. Photoprotein-like genes were also identified in the genomes of the non-luminescent sponge Amphimedon queenslandica and the non-luminescent cnidarian Nematostella vectensis, and phylogenomic analysis demonstrated that photoprotein genes arose at the base of all animals. Photoprotein gene expression in Mnemiopsis embryos begins during gastrulation in migrating precursors to photocytes and persists throughout development in the canals where photocytes reside. We identified three putative opsin genes in the Mnemiopsis genome and show that they do not group with well-known bilaterian opsin subfamilies. Interestingly, photoprotein transcripts are co-expressed with two of the putative opsins in developing photocytes. Opsin expression is also seen in the apical sensory organ. We present evidence that one opsin functions as a photopigment in vitro, absorbing light at wavelengths that overlap with peak photoprotein light emission, raising the hypothesis that light production and light reception may be functionally connected in ctenophore photocytes. We also present genomic evidence of a complete ciliary phototransduction cascade in Mnemiopsis.ConclusionsThis study elucidates the genomic organization, evolutionary history, and developmental expression of photoprotein and opsin genes in the ctenophore Mnemiopsis leidyi, introduces a novel dual role for ctenophore photocytes in both bioluminescence and phototransduction, and raises the possibility that light production and light reception are linked in this early-branching non-bilaterian animal.

Little is known about venom in young developmental stages of animals. The appearance of toxins and stinging cells during early embryonic stages in the sea anemone Nematostella vectensis suggests that venom is already expressed in eggs and larvae of this species. Here, we harness transcriptomic, biochemical and transgenic tools to study venom production dynamics in Nematostella. We find that venom composition and arsenal of toxin-producing cells change dramatically between developmental stages of this species. These findings can be explained by the vastly different interspecific interactions of each life stage, as individuals develop from a miniature non-feeding mobile planula to a larger sessile polyp that predates on other animals and interact differently with predators. Indeed, behavioral assays involving prey, predators and Nematostella are consistent with this hypothesis. Further, the results of this work suggest a much wider and dynamic venom landscape than initially appreciated in animals with a complex life cycle.

SummaryIon selectivity of metazoan voltage-gated Na+ channels is critical for neuronal signaling and has long been attributed to a ring of four conserved amino acids that constitute the ion selectivity filter (SF) at the channel pore. Yet, in addition to channels with a preference for Ca2+ ions, the expression and characterization of Na+ channel homologs from the sea anemone Nematostella vectensis, a member of the early-branching metazoan phylum Cnidaria, revealed a sodium-selective channel bearing a noncanonical SF. Mutagenesis and physiological assays suggest that pore elements additional to the SF determine the preference for Na+ in this channel. Phylogenetic analysis assigns the Nematostella Na+-selective channel to a channel group unique to Cnidaria, which diverged >540 million years ago from Ca2+-conducting Na+ channel homologs. The identification of Cnidarian Na+-selective ion channels distinct from the channels of bilaterian animals indicates that selectivity for Na+ in neuronal signaling emerged independently in these two animal lineages.

Gene families, which encode toxins, are found in many poisonous animals, yet there is limited understanding of their evolution at the nucleotide level. The release of the genome draft sequence for the sea anemone Nematostella vectensis enabled a comprehensive study of a gene family whose neurotoxin products affect voltage-gated sodium channels. All gene family members are clustered in a highly repetitive approximately 30-kb genomic region and encode a single toxin, Nv1. These genes exhibit extreme conservation at the nucleotide level which cannot be explained by purifying selection. This conservation greatly differs from the toxin gene families of other animals (e.g., snakes, scorpions, and cone snails), whose evolution was driven by diversifying selection, thereby generating a high degree of genetic diversity. The low nucleotide diversity at the Nv1 genes is reminiscent of that reported for DNA encoding ribosomal RNA (rDNA) and 2 hsp70 genes from Drosophila, which have evolved via concerted evolution. This evolutionary pattern was experimentally demonstrated in yeast rDNA and was shown to involve unequal crossing-over. Through sequence analysis of toxin genes from multiple N. vectensis populations and 2 other anemone species, Anemonia viridis and Actinia equina, we observed that the toxin genes for each sea anemone species are more similar to one another than to those of other species, suggesting they evolved by manner of concerted evolution. Furthermore, in 2 of the species (A. viridis and A. equina) we found genes that evolved under diversifying selection, suggesting that concerted evolution and accelerated evolution may occur simultaneously.

Distinguishing natural versus anthropogenic dispersal of organisms is essential for determining the native range of a species and implementing an effective conservation strategy. For cryptogenic species with limited historical records, molecular data can help to identify introductions.

The establishment of host-bacterial colonization during development is a fundamental process influencing the fitness of many organisms, but the factors controlling community membership and influencing the establishment of the microbial ecosystem during development are poorly understood. The starlet sea anemone Nematostella vectensis serves as a cnidarian model organism due to the availability of laboratory cultures and its high tolerance for broad ranges of salinity and temperature. Here, we show that the anemone's epithelia are colonized by diverse bacterial communities and that the composition of its microbiota is tightly coupled to host development. Environmental variations led to robust adjustments in the microbial composition while still maintaining the ontogenetic core signature. In addition, analysis of bacterial communities of Nematostella polyps from five different populations revealed a strong correlation between host biogeography and bacterial diversity despite years of laboratory culturing. These observed variations in fine-scale community composition following environmental change and for individuals from different geographic origins could represent the microbiome's contribution to host acclimation and potentially adaptation, respectively, and thereby contribute to the maintenance of homeostasis due to environmental changes.