Horizontal transfer of retrotransposons between bivalves and other aquatic species of multiple phyla.

The LTR retrotransposon Steamer is a selfish endogenous element in the soft-shell clam genome that was first detected because of its dramatic amplification in bivalve transmissible neoplasia afflicting the species. We amplified and sequenced related retrotransposons from the genomic DNA of many other bivalve species, finding evidence of horizontal transfer of retrotransposons from the genome of one species to another. First, the phylogenetic tree of the Steamer-like elements from 19 bivalve species is markedly discordant with host phylogeny, suggesting frequent cross-species transfer throughout bivalve evolution. Second, sequences nearly identical to Steamer were identified in the genomes of Atlantic razor clams and Baltic clams, indicating recent transfer. Finally, a search of the National Center for Biotechnology Information sequence database revealed that Steamer-like elements are present in the genomes of completely unrelated organisms, including zebrafish, sea urchin, acorn worms, and coral. Phylogenetic incongruity, a patchy distribution, and a higher similarity than would be expected by vertical inheritance all provide evidence for multiple long-distance cross-phyla horizontal transfer events. These data suggest that over both short- and long-term evolutionary timescales, Steamer-like retrotransposons, much like retroviruses, can move between organisms and integrate new copies into new host genomes.

A phylogenomic framework, evolutionary timeline and genomic resources for comparative studies of decapod crustaceans.

Comprising over 15 000 living species, decapods (crabs, shrimp and lobsters) are the most instantly recognizable crustaceans, representing a considerable global food source. Although decapod systematics have received much study, limitations of morphological and Sanger sequence data have yet to produce a consensus for higher-level relationships. Here, we introduce a new anchored hybrid enrichment kit for decapod phylogenetics designed from genomic and transcriptomic sequences that we used to capture new high-throughput sequence data from 94 species, including 58 of 179 extant decapod families, and 11 of 12 major lineages. The enrichment kit yields 410 loci (greater than 86 000 bp) conserved across all lineages of Decapoda, more clade-specific molecular data than any prior study. Phylogenomic analyses recover a robust decapod tree of life strongly supporting the monophyly of all infraorders, and monophyly of each of the reptant, 'lobster' and 'crab' groups, with some results supporting pleocyemate monophyly. We show that crown decapods diverged in the Late Ordovician and most crown lineages diverged in the Triassic-Jurassic, highlighting a cryptic Palaeozoic history, and post-extinction diversification. New insights into decapod relationships provide a phylogenomic window into morphology and behaviour, and a basis to rapidly and cheaply expand sampling in this economically and ecologically significant invertebrate clade.

Worms that suck: Phylogenetic analysis of Hirudinea solidifies the position of Acanthobdellida and necessitates the dissolution of Rhynchobdellida.

Annelids possessing a posterior sucker and a fixed number of somites - most famously leeches (Hirudinida), but also crayfish worms (Branchiobdellida) and salmonid parasites (Acanthobdellida) - form a clade; however, determining the relationships between these orders has proven challenging. Here, we compile the largest molecular phylogenetic dataset yet analysed for these groups, including new sequences for key taxa. We find robust model-based support for a clade formed by Hirudinida and Acanthobdellida, contrasting the largest prior studies. We determine that conflicting prior studies included contaminant sequences for Acanthobdella peledina. In addition to this broad-scale comparison, the size of our dataset grants us invaluable insight into the internal relationships of leeches and crayfish worms. Of particular importance, a largely marine clade of leeches (Piscicolidae and Ozobranchidae) is recovered as sister to all remaining Hirudinida. This necessitates the dissolution of the paraphyletic suborder Rhynchobdellida into two new suborders (Oceanobdelliformes and Glossiphoniiformes). Likewise, we decompose Arhynchobdellida into its respective suborders: Hirudiniformes, Erpobdelliformes, and the new, monotypic, Americobdelliformes.

Phylogenetic analysis of Placobdella (Hirudinea: Rhynchobdellida: Glossiphoniidae) with consideration of COI variation.

Placobdella is a genus of blood-feeding leeches in the family Glossiphoniidae. Historically, species of Placobdella have posed difficulty for systematists owing to a lack of informative morphological characters and the preponderance of inadequate or incomplete species descriptions. Here, we conduct a phylogenetic analysis of 55 individuals representing 20 of the 24 currently recognized nominal taxa using COI, ND1, 12S rDNA and ITS sequences under parsimony, maximum likelihood and Bayesian inference. We also examine the isolated COI phylogeny for the genus using an expanded dataset encompassing three additional species not included in the concatenated dataset. Finally, we assess genetic variation at the COI locus to validate initial specimen identifications and estimate how COI variation may reflect species boundaries. We conclude that Placobdella is a monophyletic group that places as the sister group to a clade formed by the genera Haementeria and Helobdella. We discuss the evolutionary implications of several internal relationships that are robustly resolved by all three optimality criteria, paying particular attention to the apparent fluidity of morphological characters exhibited by members of Placobdella. We also find preliminary evidence for the presence of cryptic and undescribed diversity within the genus.

Evidence of horizontal transmission of the cancer-associated Steamer retrotransposon among ecological cohort bivalve species.

Bivalve specimens from legacy frozen tissue collections, and others freshly obtained, were surveyed for the presence of the Steamer long terminal repeat (LTR)-retrotransposon associated with disseminated hemic neoplasia of the soft-shelled clam Mya areneria. Of 22 species investigated using primers for the pol region, only Atlantic M. arenaria, Atlantic and North Sea razor clams Ensis directus, and Baltic clams Macoma balthica from the North Sea were found to possess copies of Steamer in their genomes. Notably, close relatives like Mya truncata and Siliqua patula did not exhibit evidence of Steamer. Amplified Steamer sequences were uniformly identical in all M. areneria specimens, and were highly variable across specimens of E. directus. Variation in the latter included nucleotide polymorphisms among and within individuals as well as length variation in 2 specimens corresponding to the deletion of a predicted stable hairpin structure. Results implicate Atlantic razor clams as the proximal source for horizontal transmission of Steamer among ecologically similar yet markedly distantly related bivalves. The consequences of cross-species transmission of the Steamer retrotransposon are unknown, and the finding of Steamer in 3 bivalve species suggests that further spread is possible.

In silico hybridization enables transcriptomic illumination of the nature and evolution of Myxozoa.

BACKGROUND: The Myxozoa, a group of oligocellular, obligate endoparasites, has long been poorly understood in an evolutionary context. Recent genome-level sequencing techniques such as RNA-seq have generated large amounts of myxozoan sequence data, providing valuable insight into their evolutionary history. However, sequences from host tissue contamination are present in next-generation sequencing reactions of myxozoan tissue, and differentiating between the two has been inadequately addressed. In order to shed light on the genetic underpinnings of myxozoan biology, assembled contigs generated from these studies that derived from the myxozoan must be decoupled from transcripts derived from host tissue and other contamination. This study describes a pipeline for categorization of transcripts asmyxozoan based on similarity searching with known host and parasite sequences, explores the extent to which host contamination is present in previously existing myxozoan datasets, and implements this pipeline on a newly sequenced transcriptome of Myxobolus pendula, a parasite of the common creek chub gill arch. METHODS: The insilico hybridization pipeline uses iterative BLAST searching and database-driven e-value comparison to categorize transcripts as deriving from host, parasite, or other contamination. Functional genetic analysis of M. pendula was conducted using further BLAST searching, Hidden Markov Modeling, and sequence alignment and phylogenetic reconstruction. RESULTS: Three RNA libraries of encysted M. pendula plasmodia were sequenced and subjected to the method. Nearly half of the final set of contiguous assembly sequences (47.3 %) was identified as putative myxozoan transcripts. Putative contamination was also identified in at least 1/3(rd) of previously published myxozoan transcripts. The set of M. pendula transcripts was mined for a range of biologically insightful genes, including taxonomically restricted nematocyst structural proteins and nematocyst proteins identified through mass tandem spectrometry of other cnidarians. Several novel findings emerged, including a fourth myxozoan minicollagen gene, putative myxozoan toxin proteins,and extracellular matrix glycoproteins. CONCLUSIONS: This study serves as a model for the handling of next-generation myxozoan sequence. The need for careful categorization was demonstrated in both previous and new sets of myxozoan sequences. The final set of confidently assigned myxozoan transcripts can be mined for any biologically relevant gene or gene family without spurious misidentification of host contamination as a myxozoan homolog. As exemplified by M. pendula, the repertoire of myxozoan polar capsules may be more complex than previously thought, with an additional minicollagen homolog and putative expression of toxin proteins.

iDNA from terrestrial haematophagous leeches as a wildlife surveying and monitoring tool - prospects, pitfalls and avenues to be developed.

Invertebrate-derived DNA (iDNA) from terrestrial haematophagous leeches has recently been proposed as a powerful non-invasive tool with which to detect vertebrate species and thus to survey their populations. However, to date little attention has been given to whether and how this, or indeed any other iDNA-derived data, can be combined with state-of-the-art analytical tools to estimate wildlife abundances, population dynamics and distributions. In this review, we discuss the challenges that face the application of existing analytical methods such as site-occupancy and spatial capture-recapture (SCR) models to terrestrial leech iDNA, in particular, possible violations of key assumptions arising from factors intrinsic to invertebrate parasite biology. Specifically, we review the advantages and disadvantages of terrestrial leeches as a source of iDNA and summarize the utility of leeches for presence, occupancy, and spatial capture-recapture models. The main source of uncertainty that attends species detections derived from leech gut contents is attributable to uncertainty about the spatio-temporal sampling frame, since leeches retain host-blood for months and can move after feeding. Subsequently, we briefly address how the analytical challenges associated with leeches may apply to other sources of iDNA. Our review highlights that despite the considerable potential of leech (and indeed any) iDNA as a new survey tool, further pilot studies are needed to assess how analytical methods can overcome or not the potential biases and assumption violations of the new field of iDNA. Specifically we argue that studies to compare iDNA sampling with standard survey methods such as camera trapping, and those to improve our knowledge on leech (and other invertebrate parasite) physiology, taxonomy, and ecology will be of immense future value.

Diversity of features of the female reproductive system and other morphological characters in leeches (Citellata, Hirudinida) in phylogenetic conception.

An epistemological-evolutionary conception of leeches (Hirudinida) based on features of the female reproductive system in combination with other morphological characters is presented in the spirit of the cladistic school of taxonomy. Characters relating to the structure of the ovary and the course of oogenesis in leeches were interpreted in this manner, for the first time. Each study was conducted on type species of higher taxonomic groups of true leeches. Results of analyses using features of the reproductive system only as well as in combination with other morphological characters show Piscicolidae and Glossiphoniidae as sister clades making Rhynchobdellida a monophyletic group. Also, Hirudiniformes and Erpobdelliformes appeared to be sister clades within Arhynchobdellida. The relationship between the outgroup specimens and leeches remained unresolved, because both Acanthobdella peledina and branchiobdellidans appeared to be in an equivocal relationship to hirudinidans. Characters concerning the structure of the female reproductive system and course of oogenesis thus appeared to be useful, although conservative, for reconstruction of leech phylogeny, and they well reflect phylogenetic relationships of Hirudinida at the family level.

Phylogenomics of Annelida revisited: a cladistic approach using genome-wide expressed sequence tag data mining and examining the effects of missing data.

We present phylogenomic analyses of the most comprehensive molecular character set compiled for Annelida and its constituent taxa, including over 347 000 aligned nucleotide sites for 39 taxa. The nucleotide data set was recovered using a pre-existing amino acid data set of almost 48 000 aligned sites as a backbone for tBLASTn searches against NCBI. In addition, orthology determinations of the loci in the original amino acid data set were scrutinized using an All vs All Reciprocal Best Hit approach, employing BLASTp, and examining for statistical interdependency among the loci. This approach revealed considerable sequence redundancy among the loci in the original data set and a new data set was compiled, with the redundancy removed. The newly compiled nucleotide data set, the original amino acid data set, and the new reduced amino acid data set were subjected to parsimony analyses and two forms of bootstrap resampling. The last-named data set also was analysed using a maximum-likelihood approach. There were two main objectives to these analyses: (i) to examine the general topology, including support, resulting from the analyses of the new data sets and (ii) to assess the consistency of the branching patterns across optimality criteria by comparison with previous probabilistic approaches. The phylogenetic hypotheses resulting from analyses of the three data sets are largely unsupported, reflecting the continued difficulty of finding numerous, reliable, and suitable loci for a group as ancient as Annelida. Resulting parsimonious hypotheses disagree, in some respects, with the previous probabilistic approaches; Sedentaria and, in most cases, Errantia are not supported as monophyletic groups but Pleistoannelida is recovered as a (unsupported) monophyletic group in one of the three parsimony analyses as well as the likelihood analysis. In addition, we performed missing data titration studies to estimate the impact of missing data on overall support and support for specific clades.

Systematics and evolution of syllids (Annelida, Syllidae).

A large, combined phylogenetic analysis (including morphological and molecular data from 18S rDNA, 16S rDNA and cytochrome c oxidase subunit I), with the highest number of species and genera of Syllidae studied to date (213 terminals), is examined. The data were explored with different parameters and optimality criteria (parsimony, likelihood, and bayesian inference). The monophyly of Syllidae and most of the traditional subfamilies is supported. The subfamily Eusyllinae is polyphyletic, as currently delineated, but it is herein reorganized and its diagnosis modified to be a valid group. Additional well supported clades arise. The phylogenetic relationships of the well known and established genera, as well as several enigmatic genera (e.g. Anguillosyllis, Paraopisthosyllis and Parahaplosyllis), the position of which in syllid taxonomy was uncertain or dubious to date, are clarified. The results corroborate previous hypotheses about the evolution of the reproductive and brooding modes. Within Syllinae, the nature of the stolon is phylogenetically informative. The classification of the whole family is revised and discussed on the basis of this phylogenetic hypothesis. © The Willi Hennig Society 2011.

Bacterial symbiont and salivary peptide evolution in the context of leech phylogeny.

The evolutionary history of leeches is employed as a general framework for understanding more than merely the systematics of this charismatic group of annelid worms, and serves as a basis for understanding blood-feeding related correlates ranging from the specifics of gut-associated bacterial symbionts to salivary anticoagulant peptides. A variety of medicinal leech families were examined for intraluminal crop bacterial symbionts. Species of Aeromonas and Bacteroidetes were characterized with DNA gyrase B and 16S rDNA. Bacteroidetes isolates were found to be much more phylogenetically diverse and suggested stronger evidence of phylogenetic correlation than the gammaproteobacteria. Patterns that look like co-speciation with limited taxon sampling do not in the full context of phylogeny. Bioactive compounds that are expressed as gene products, like those in leech salivary glands, have 'passed the test' of evolutionary selection. We produced and bioinformatically mined salivary gland EST libraries across medicinal leech lineages to experimentally and statistically evaluate whether evolutionary selection on peptides can identify structure-function activities of known therapeutically relevant bioactive compounds like antithrombin, hirudin and antistasin. The combined information content of a well corroborated leech phylogeny and broad taxonomic coverage of expressed proteins leads to a rich understanding of evolution and function in leech history.

The Origin and Evolution of Antistasin-like Proteins in Leeches (Hirudinida, Clitellata).

Bloodfeeding is employed by many parasitic animals and requires specific innovations for efficient feeding. Some of these innovations are molecular features that are related to the inhibition of hemostasis. For example, bloodfeeding insects, bats, and leeches release proteins with anticoagulatory activity through their salivary secretions. The antistasin-like protein family, composed of serine protease inhibitors with one or more antistasin-like domains, is tightly linked to inhibition of hemostasis in leeches. However, this protein family has been recorded also in non-bloodfeeding invertebrates, such as cnidarians, mollusks, polychaetes, and oligochaetes. The present study aims to 1) root the antistasin-like gene tree and delimit the major orthologous groups, 2) identify potential independent origins of salivary proteins secreted by leeches, and 3) identify major changes in domain and/or motif structure within each orthologous group. Five clades containing leech antistasin-like proteins are distinguishable through rigorous phylogenetic analyses based on nine new transcriptomes and a diverse set of comparative data: the trypsin + leukocyte elastase inhibitors clade, the antistasin clade, the therostasin clade, and two additional, unnamed clades. The antistasin-like gene tree supports multiple origins of leech antistasin-like proteins due to the presence of both leech and non-leech sequences in one of the unnamed clades, but a single origin of factor Xa and trypsin + leukocyte elastase inhibitors. This is further supported by three sequence motifs that are exclusive to antistasins, the trypsin + leukocyte elastase inhibitor clade, and the therostasin clade, respectively. We discuss the implications of our findings for the evolution of this diverse family of leech anticoagulants.

Unringing a bell: metazoan phylogenomics and the partition bootstrap.

Several large phylogenomic analyses have recently cast doubt on long-held beliefs about early metazoan phylogenetic patterns. Those data sets, and the relative bootstrap support for various controversial clades, are reanalysed in the context of parsimony, yielding results that are at considerable odds with the original likelihood or Bayesian findings. Discrepancies are considered in light of the tendency of RAxML to overestimate support values by virtue (sic) of its lazy search algorithm and its autocorrelated pseudoreplication as well as the extraordinary ability for Bayesian analyses to be led astray by missing data. In addition to standard nonparametric bootstrapping as a measure of support, a new strategy involving resampling loci as units, partition bootstrap support, is introduced as a more defensible alternative to resampling nonindependent sites. © The Willi Hennig Society 2009.

Draft genome of the European medicinal leech Hirudo medicinalis (Annelida, Clitellata, Hirudiniformes) with emphasis on anticoagulants.

The European medicinal leech has been used for medicinal purposes for millennia, and continues to be used today in modern hospital settings. Its utility is granted by the extremely potent anticoagulation factors that the leech secretes into the incision wound during feeding and, although a handful of studies have targeted certain anticoagulants, the full range of anticoagulation factors expressed by this species remains unknown. Here, we present the first draft genome of the European medicinal leech, Hirudo medicinalis, and estimate that we have sequenced between 79-94% of the full genome. Leveraging these data, we searched for anticoagulation factors across the genome of H. medicinalis. Following orthology determination through a series of BLAST searches, as well as phylogenetic analyses, we estimate that fully 15 different known anticoagulation factors are utilized by the species, and that 17 other proteins that have been linked to antihemostasis are also present in the genome. We underscore the utility of the draft genome for comparative studies of leeches and discuss our results in an evolutionary context.

Multilocus Metabarcoding of Terrestrial Leech Bloodmeal iDNA Increases Species Richness Uncovered in Surveys of Vertebrate Host Biodiversity.

Leech-derived invertebrate DNA (iDNA) has been successfully leveraged to conduct surveys of vertebrate host biodiversity across the Indo Pacific. However, this technique has been limited methodologically, typically only targeting mammalian 16S rDNA, or both 16S and vertebrate 12S rDNA for leech host determination. To improve the taxonomic richness of vertebrate host species in iDNA surveys, we re-analyze datasets from Bangladesh, Cambodia, China, and Madagascar through metabarcoding via next generation sequencing (NGS) of 12S, 16S (2 types, one designed to target mammals and the other, residual eDNA), nicotinamide adenine dinucleotide hydride dehydrogenase 2 (ND2), and cytochrome c oxidase subunit 1 (COI). With our 5 primer sets, we identify 41 unique vertebrate hosts to the species level, among 1,200 leeches analyzed, along with an additional 13 taxa to the family rank. Within our 41 taxa, we note that adding ND2 and COI loci increased species richness detection by 25%. NGS has emerged as more efficient than Sanger sequencing for large scale metabarcoding applications and, with the decline in cost of NGS, our pooled sample multilocus protocol is an attractive option for iDNA biodiversity surveys.

Poly-paraphyly of Hirudinidae: many lineages of medicinal leeches.

BACKGROUND: Medicinal leeches became infamous for their utility in bloodletting popularized in the 19th century, and have seen a recent resurgence in post-operative treatments for flap and replantation surgeries, and in terms of characterization of salivary anticoagulants. Notorious throughout the world, the quintessential leech family Hirudinidae has been taken for granted to be monophyletic, as has the non-bloodfeeding family Haemopidae. RESULTS: This study is the first to evaluate molecular evidence from hirudinid and haemopid leeches in a manner that encompasses the global scope of their taxonomic distributions. We evaluated the presumed monophyly of the Hirudinidae and assessed previous well-accepted classification schemes. The Hirudinidae were found not to be monophyletic, falling instead into two distinct and unrelated clades. Members of the non-bloodfeeding family Haemopidae were scattered throughout the tree and among traditional hirudinid genera. A combination of nuclear 18S rDNA and 28S rDNA with mitochondrial 12S rDNA and cytochrome c oxidase I were analyzed with Parsimony and with Bayesian methods. CONCLUSION: The family Hirudinidae must be refined to include only the clade containing Hirudo medicinalis (European medicinal leech) and related leeches irrespective of bloodfeeding behavior. A second clade containing Macrobdella decora (North American medicinal leech) and its relatives may yet be recognized in Semiscolecidae in order to avoid paraphyly. The African distribution of species from each of the divergent hirudinid clades suggests that a deep divergence took place in the history of the medicinal leeches hundreds of millions of years ago.

Ideating iDNA: Lessons and limitations from leeches in legacy collections.

Indirect methods for conducting faunal inventories present great promise, and genomic inventories derived from environmental sources (eDNA) are improving. Invertebrate ingested DNA (iDNA) from terrestrial leeches in the family Haemadipsidae has shown potential for surveying vertebrates and biodiversity monitoring in protected areas. Here we present an initial, and critical, evaluation of the limitations and biases of current iDNA protocols for biodiversity monitoring using both standard and NGS barcoding approaches. Key findings include the need for taxon relevant multi-locus markers and reference databases. In particular, the limitations of available reference databases have profound potential to mislead and bias eDNA and iDNA results if not critically interpreted. Nevertheless, there is great potential for recovery of amplifiable DNA from gut contents of invertebrate museum specimens which may reveal both temporal patterns and cryptic diversity in protected areas with increased efficiency. Our analyses of ingested DNA (iDNA) from both freshly stored and previously collected (legacy) samples of terrestrial leeches successfully identified vertebrates from Myanmar, Australia and Madagascar and indicate the potential to characterize microbial communities, pathogen diversity and interactions at low cost.

Broad geographic sampling and DNA barcoding do not support the presence of Helobdella stagnalis (Linnaeus, 1758) (Clitellata: Glossiphoniidae) in North America.

The description of Helobdella stagnalis (Linnaeus, 1758) has emphasized the presence of a nuchal, chitinous scute located on the dorsal surface in the first third of the body as the diagnostic character for the species. Historically, identifications of species of Helobdella have relied heavily on this character and, as a result, Helobdella stagnalis has been reported from an inordinately broad geographic range, including Europe, Asia, Africa, North America, and South America. In addition to a few earlier investigations, a recent analysis showed that great genetic distances (orders of magnitude greater than previous estimations of intraspecific divergence in leeches) are present between scute-bearing specimens identified as H. stagnalis from Europe and North America, implying that H. stagnalis does not occur in North America. The present study expands the geographic boundaries of taxon sampling for both European and North American taxa, and re-examines the phylogenetic relationships and cytochrome c oxidase subunit I (COI) variation within scute-bearing species of the genus Helobdella. Our analyses include specimens putatively identified as "Helobdella stagnalis" from Sweden, Norway, Iceland, England, France, Italy, Slovenia, Turkey, Russia, and Iran, as well as numerous localities covering Canada and the USA. Our results corroborate previous studies in that European and west Asian specimens form a clade, including the neotype, which is separate from North American taxa. To alleviate future taxonomic confusion, we redescribe H. stagnalis and designate a neotype from the inferred type locality. The designation of a neotype stabilizes the taxonomy of scute-bearing leeches of the genus Helobdella and enables us to definitively correct erroneous identifications reported in previous studies. We also note that at least four lineages of scute-bearing, North American species of Helobdella lack formal descriptions.

Characterization of the digestive-tract microbiota of Hirudo orientalis, a european medicinal leech.

FDA-approved, postoperative use of leeches can lead to bacterial infections. In this study, we used culture-dependent and culture-independent approaches to characterize the digestive-tract microbiota of Hirudo orientalis. Surprisingly, two Aeromonas species, A. veronii and A. jandaei, were cultured. Uncultured Rikenella-like bacteria were most similar to isolates from Hirudo verbana.

The transcriptome of the Bermuda fireworm Odontosyllis enopla (Annelida: Syllidae): A unique luciferase gene family and putative epitoky-related genes.

The Bermuda fireworm Odontosyllis enopla exhibits an extremely tight circalunar circadian behavior that results in an impressive bioluminescent mating swarm, thought to be due to a conventional luciferase-mediated oxidation of a light-emitting luciferin. In addition, the four eyes become hypertrophied and heavily pigmented, and the nephridial system is modified to store and release gametes and associated secretions. In an effort to elucidate transcripts related to bioluminescence, circadian or circalunar periodicity, as well as epitoky-related changes of the eyes and nephridial system, we examined the transcriptomic profile of three female O. enopla during a bioluminescent swarm in Ferry Reach, Bermuda. Using the well-characterized luciferase gene of the Japanese syllid Odontosyllis undecimdonta as a reference, a complete best-matching luciferase open reading frame (329 amino acids in length) was found in all three individuals analyzed in addition to numerous other paralogous sequences in this new gene family. No photoproteins were detected. We also recovered a predicted homolog of 4-coumarate-CoA ligase (268 amino acids in length) that best matched luciferase of the firefly Luciola with the best predicted template being the crystal structure of luciferase for Photinus pyralis, the common eastern firefly. A wide variety of genes associated with periodicity were recovered including predicted homologs of clock, bmal1, period, and timeless. Several genes corresponding to putative epitoky-related changes of the eyes were recovered including predicted homologs of a phototransduction gene, a retinol dehydrogenase and carotenoid isomerooxygenase as well as a visual perception related retinal rod rhodopsin-sensitive cGMP 3',5'-cyclic phosphodiesterase. Genes correlating to putative epitoky-related changes of the nephridia included predicted homologs of nephrocystin-3 and an egg-release sex peptide receptor.