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.

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.

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.

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.

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.

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.

Debugging diversity - a pan-continental exploration of the potential of terrestrial blood-feeding leeches as a vertebrate monitoring tool.

The use of environmental DNA (eDNA) has become an applicable noninvasive tool with which to obtain information about biodiversity. A subdiscipline of eDNA is iDNA (invertebrate-derived DNA), where genetic material ingested by invertebrates is used to characterize the biodiversity of the species that served as hosts. While promising, these techniques are still in their infancy, as they have only been explored on limited numbers of samples from only a single or a few different locations. In this study, we investigate the suitability of iDNA extracted from more than 3,000 haematophagous terrestrial leeches as a tool for detecting a wide range of terrestrial vertebrates across five different geographical regions on three different continents. These regions cover almost the full geographical range of haematophagous terrestrial leeches, thus representing all parts of the world where this method might apply. We identify host taxa through metabarcoding coupled with high-throughput sequencing on Illumina and IonTorrent sequencing platforms to decrease economic costs and workload and thereby make the approach attractive for practitioners in conservation management. We identified hosts in four different taxonomic vertebrate classes: mammals, birds, reptiles and amphibians, belonging to at least 42 different taxonomic families. We find that vertebrate blood ingested by haematophagous terrestrial leeches throughout their distribution is a viable source of DNA with which to examine a wide range of vertebrates. Thus, this study provides encouraging support for the potential of haematophagous terrestrial leeches as a tool for detecting and monitoring terrestrial vertebrate biodiversity.

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.

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.

Marine Leech Anticoagulant Diversity and Evolution.

Leeches (Annelida: Hirudinea) possess powerful salivary anticoagulants and, accordingly, are frequently employed in modern, authoritative medicine. Members of the almost exclusively marine family Piscicolidae account for 20% of leech species diversity, and they feed on host groups (e.g., sharks) not encountered by their freshwater and terrestrial counterparts. Moreover, some species of Ozobranchidae feed on endangered marine turtles and have been implicated as potential vectors for the tumor-associated turtle herpesvirus. In spite of their ecological importance and unique host associations, there is a distinct paucity of data regarding the salivary transcriptomes of either of these families. Using next-generation sequencing, we profiled transcribed, putative anticoagulants and other salivary bioactive compounds that have previously been linked to blood feeding from 7 piscicolid species (3 elasmobranch feeders; 4 non-cartilaginous fish feeders) and 1 ozobranchid species (2 samples). In total, 149 putative anticoagulants and bioactive loci were discovered in varying constellations throughout the different samples. The putative anticoagulants showed a broad spectrum of described antagonistic pathways, such as inhibition of factor Xa and platelet aggregation, which likely have similar bioactive roles in marine fish and turtles. A transcript with homology to ohanin, originally isolated from king cobras, was found in Cystobranchus vividus but is otherwise unknown from leeches. Estimation of selection pressures for the putative anticoagulants recovered evidence for both positive and purifying selection along several isolated branches in the gene trees, and positive selection was also estimated for a few select codons in a variety of marine species. Similarly, phylogenetic analyses of the amino acid sequences for several anticoagulants indicated divergent evolution.

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.

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.

The mitogenome of the bed bug Cimex lectularius (Hemiptera: Cimicidae).

We report the extraction of a bed bug mitogenome from high-throughput sequencing projects originally focused on the nuclear genome of Cimex lectularius. The assembled mitogenome has a similar AT nucleotide composition bias found in other insects. Phylogenetic analysis of all protein-coding genes indicates that C. lectularius is clearly a member of a paraphyletic Cimicomorpha clade within the Order Hemiptera.

Comparative Mitogenomics of Leeches (Annelida: Clitellata): Genome Conservation and Placobdella-Specific trnD Gene Duplication.

Mitochondrial DNA sequences, often in combination with nuclear markers and morphological data, are frequently used to unravel the phylogenetic relationships, population dynamics and biogeographic histories of a plethora of organisms. The information provided by examining complete mitochondrial genomes also enables investigation of other evolutionary events such as gene rearrangements, gene duplication and gene loss. Despite efforts to generate information to represent most of the currently recognized groups, some taxa are underrepresented in mitochondrial genomic databases. One such group is leeches (Annelida: Hirudinea: Clitellata). Herein, we expand our knowledge concerning leech mitochondrial makeup including gene arrangement, gene duplication and the evolution of mitochondrial genomes by adding newly sequenced mitochondrial genomes for three bloodfeeding species: Haementeria officinalis, Placobdella lamothei and Placobdella parasitica. With the inclusion of three new mitochondrial genomes of leeches, a better understanding of evolution for this organelle within the group is emerging. We found that gene order and genomic arrangement in the three new mitochondrial genomes is identical to previously sequenced members of Clitellata. Interestingly, within Placobdella, we recovered a genus-specific duplication of the trnD gene located between cox2 and atp8. We performed phylogenetic analyses using 12 protein-coding genes and expanded our taxon sampling by including GenBank sequences for 39 taxa; the analyses confirm the monophyletic status of Clitellata, yet disagree in several respects with other phylogenetic hypotheses based on morphology and analyses of non-mitochondrial data.

Genome assembly and geospatial phylogenomics of the bed bug Cimex lectularius.

The common bed bug (Cimex lectularius) has been a persistent pest of humans for thousands of years, yet the genetic basis of the bed bug's basic biology and adaptation to dense human environments is largely unknown. Here we report the assembly, annotation and phylogenetic mapping of the 697.9-Mb Cimex lectularius genome, with an N50 of 971 kb, using both long and short read technologies. A RNA-seq time course across all five developmental stages and male and female adults generated 36,985 coding and noncoding gene models. The most pronounced change in gene expression during the life cycle occurs after feeding on human blood and included genes from the Wolbachia endosymbiont, which shows a simultaneous and coordinated host/commensal response to haematophagous activity. These data provide a rich genetic resource for mapping activity and density of C. lectularius across human hosts and cities, which can help track, manage and control bed bug infestations.

Comparative Transcriptomic Analyses of Three Species of Placobdella (Rhynchobdellida: Glossiphoniidae) Confirms a Single Origin of Blood Feeding in Leeches.

One of the recalcitrant questions regarding the evolutionary history of clitellate annelids involves the feeding preference of the common ancestor of extant rhynchobdellid (proboscis bearing) and arhynchobdellid (jaw bearing) leeches. Whereas early evidence, based on morphological data, pointed towards independent acquisitions of blood feeding in the 2 orders, molecular-based phylogenetic data suggest that the ancestor of modern leeches was a sanguivore. Here, we use a comparative transcriptomic approach in order to increase our understanding of the diversity of anticoagulation factors for 3 species of the genus Placobdella, for which comparative data have been lacking, and inspect these in light of archetypal anticoagulant data for both arhynchobdellid and other rhynchobdellid species. Notwithstanding the varying levels of host specificity displayed by the 3 different species of Placobdella, transcriptomic profiles with respect to anticoagulation factors were largely similar -this despite the fact that Placobdella kwetlumye only retains a single pair of salivary glands, as opposed to the 2 pairs more common in the genus. Results show that 9 different anticoagulant proteins and an additional 5 putative antihemostasis proteins are expressed in salivary secretions of the 3 species. In particular, an ortholog of the archetypal, single-copy, anticoagulant hirudin (not previously available as comparative data for rhynchobdellids) is present in at least 2 of 3 species examined, corroborating the notion of a single origin of blood feeding in the ancestral leech.

Multiplexed pyrosequencing of nine sea anemone (Cnidaria: Anthozoa: Hexacorallia: Actiniaria) mitochondrial genomes.

Six complete and three partial actiniarian mitochondrial genomes were amplified in two semi-circles using long-range PCR and pyrosequenced in a single run on a 454 GS Junior, doubling the number of complete mitogenomes available within the order. Typical metazoan mtDNA features included circularity, 13 protein-coding genes, 2 ribosomal RNA genes, and length ranging from 17,498 to 19,727 bp. Several typical anthozoan mitochondrial genome features were also observed including the presence of only two transfer RNA genes, elevated A + T richness ranging from 54.9 to 62.4%, large intergenic regions, and group 1 introns interrupting NADH dehydrogenase subunit 5 and cytochrome c oxidase subunit I, the latter of which possesses a homing endonuclease gene. Within the sea anemone Alicia sansibarensis, we report the first mitochondrial gene order rearrangement within the Actiniaria, as well as putative novel non-canonical protein-coding genes. Phylogenetic analyses of all 13 protein-coding and 2 ribosomal genes largely corroborated current hypotheses of sea anemone interrelatedness, with a few lower-level differences.