Timing and Scope of Genomic Expansion within Annelida: Evidence from Homeoboxes in the Genome of the Earthworm Eisenia fetida

Annelid Comparative Genomics and the Evolution of Massive Lineage-Specific Genome Rearrangement in Bilaterians

The organization of genomes into chromosomes is critical for processes such as genetic recombination, environmental adaptation, and speciation. All animals with bilateral symmetry inherited a genome structure from their last common ancestor that has been highly conserved in some taxa but seemingly unconstrained in others. However, the evolutionary forces driving these differences and the processes by which they emerge have remained largely uncharacterized. Here, we analyze genome organization across the phylum Annelida using 23 chromosome-level annelid genomes. We find that while many annelid lineages have maintained the conserved bilaterian genome structure, the Clitellata, a group containing leeches and earthworms, possesses completely scrambled genomes. We develop a rearrangement index to quantify the extent of genome structure evolution and show that, compared to the last common ancestor of bilaterians, leeches and earthworms have among the most highly rearranged genomes of any currently sampled species. We further show that bilaterian genomes can be classified into two distinct categories-high and low rearrangement-largely influenced by the presence or absence, respectively, of chromosome fission events. Our findings demonstrate that animal genome structure can be highly variable within a phylum and reveal that genome rearrangement can occur both in a gradual, stepwise fashion, or rapid, all-encompassing changes over short evolutionary timescales.

A genome resource for the marine annelid Platynereis dumerilii

The marine annelid Platynereis dumerilii is a model organism used in many research areas including evolution and development, neurobiology, ecology and regeneration. Here we present the genomes of P. dumerilii and of the closely related P. massiliensis and P. megalops, to facilitate comparative genomic approaches and help explore Platynereis biology. We used long-read sequencing technology and chromosomal-conformation capture along with extensive transcriptomic resources to obtain and annotate a draft genome assembly of ~1.47 Gbp for P. dumerilii, of which more than half represent repeat elements. We predict around 29,000 protein-coding genes, with relatively large intron sizes, over 38,000 non-coding genes, and 580 miRNA loci. We further explore the high genetic variation (~3% heterozygosity) within the Platynereis species complex. Gene ontology reveals the most variable loci to be associated with pigmentation, development and immunity. The current work sets the stage for further development of Platynereis genomic resources.

Spatial Colinear but Broken Temporal Expression of Duplicated ParaHox Genes in Asexually Reproducing Annelids, Nais communis and Pristina longiseta

ParaHox genes are key developmental regulators involved in the patterning of the digestive tract along the anteroposterior axis and the development of the nervous system. Most studies have focused on the function of these genes in embryogenesis, while their expression patterns in postembryonic development often remain unknown. In this study, we identified for the first time all ParaHox orthologs in two naidid oligochaetes, N. communis and P. longiseta, and described their expression patterns during normal growth and fission in these animals. We showed that Gsx and Cdx are presented by two paralogs, while Xlox is a single copy gene in both species. Using whole-mount in situ hybridization, we also found that orthologs, except for the Xlox gene, have similar activity patterns with minor differences in details, while the expression patterns of paralogs can differ significantly. However, all these genes are involved in axial patterning and/or in tissue remodeling during growth and asexual reproduction in naidids. Moreover, during paratomic fission, these genes are expressed with spatial colinearity but temporal colinearity is broken. The results of this study may be evidence of the functional diversification of duplicated genes and suggest involvement of the ParaHox genes in whole-body patterning during growth and asexual reproduction in annelids.

Gene Loss may have Shaped the Cnidarian and Bilaterian Hox and ParaHox Complement

Hox and ParaHox transcription factors are important for specifying cell fates along the primary body axes during the development of most animals. Within Cnidaria, much of the research on Hox/ParaHox genes has focused on Anthozoa (anemones and corals) and Hydrozoa (hydroids) and has concentrated on the evolution and function of cnidarian Hox genes in relation to their bilaterian counterparts. Here we analyze together the full complement of Hox and ParaHox genes from species representing all four medusozoan classes (Staurozoa, Cubozoa, Hydrozoa, and Scyphozoa) and both anthozoan classes (Octocorallia and Hexacorallia). Our results show that Hox genes involved in patterning the directive axes of anthozoan polyps are absent in the stem leading to Medusozoa. For the first time, we show spatial and temporal expression patterns of Hox and ParaHox genes in the upside-down jellyfish Cassiopea xamachana (Scyphozoa), which are consistent with diversification of medusozoan Hox genes both from anthozoans and within medusozoa. Despite unprecedented taxon sampling, our phylogenetic analyses, like previous studies, are characterized by a lack of clear homology between most cnidarian and bilaterian Hox and Hox-related genes. Unlike previous studies, we propose the hypothesis that the cnidarian-bilaterian ancestor possessed a remarkably large Hox complement and that extensive loss of Hox genes was experienced by both cnidarian and bilaterian lineages.

A Chromosome-level Genome Assembly of the Highly Heterozygous Sea Urchin Echinometra sp. EZ Reveals Adaptation in the Regulatory Regions of Stress Response Genes

Echinometra is the most widespread genus of sea urchin and has been the focus of a wide range of studies in ecology, speciation, and reproduction. However, available genetic data for this genus are generally limited to a few select loci. Here, we present a chromosome-level genome assembly based on 10x Genomics, PacBio, and Hi-C sequencing for Echinometra sp. EZ from the Persian/Arabian Gulf. The genome is assembled into 210 scaffolds totaling 817.8 Mb with an N50 of 39.5 Mb. From this assembly, we determined that the E. sp. EZ genome consists of 2n = 42 chromosomes. BUSCO analysis showed that 95.3% of BUSCO genes were complete. Ab initio and transcript-informed gene modeling and annotation identified 29,405 genes, including a conserved Hox cluster. E. sp. EZ can be found in high-temperature and high-salinity environments, and we therefore compared E. sp. EZ gene families and transcription factors associated with environmental stress response ("defensome") with other echinoid species with similar high-quality genomic resources. While the number of defensome genes was broadly similar for all species, we identified strong signatures of positive selection in E. sp. EZ noncoding elements near genes involved in environmental response pathways as well as losses of transcription factors important for environmental response. These data provide key insights into the biology of E. sp. EZ as well as the diversification of Echinometra more widely and will serve as a useful tool for the community to explore questions in this taxonomic group and beyond.

Environmental Epigenetics in Soil Ecosystems: Earthworms as Model Organisms

One of the major emerging concerns within ecotoxicology is the effect of environmental pollutants on epigenetic changes, including DNA methylation, histone modifications, and non-coding RNAs. Epigenetic mechanisms regulate gene expression, meaning that the alterations of epigenetic marks can induce long-term physiological effects that can even be inherited across generations. Many invertebrate species have been used as models in environmental epigenetics, with a special focus on DNA methylation changes caused by environmental perturbations (e.g., pollution). Among soil organisms, earthworms are considered the most relevant sentinel organisms for anthropogenic stress assessment and are widely used as standard models in ecotoxicological testing of soil toxicity. In the last decade, several research groups have focused on assessing the impact of environmental stress on earthworm epigenetic mechanisms and tried to link these mechanisms to the physiological effects. The aim of this review is to give an overview and to critically examine the available literature covering this topic. The high level of earthworm genome methylation for an invertebrate species, responsiveness of epigenome to environmental stimuli, availability of molecular resources, and the possibility to study epigenetic inheritance make earthworms adequate models in environmental epigenomics. However, there are still many knowledge gaps that need to be filled in, before we can fully explore earthworms as models in this field. These include detailed characterization of the methylome using next-generation sequencing tools, exploration of multigenerational and transgenerational effects of pollutants, and information about other epigenetic mechanisms apart from DNA methylation. Moreover, the connection between epigenetic effects and phenotype has to be further explored.

Importance of clitellar tissue in the regeneration ability of earthworm Eudrilus eugeniae

Among the annelids, earthworms are renowned for their phenomenal ability to regenerate the lost segments. The adult earthworm Eudrilus eugeniae contains 120 segments and the body segments of the earthworm are divided into pre-clitellar, clitellar and post-clitellar segments. The present study denoted that clitellum plays vital role in the successful regeneration of the species. We have performed histological studies to identify among the three skin layers of the earthworm, which cellular layer supports the blastema formation and regeneration of the species. The histological evidences denoted that the proliferation of the longitudinal cell layer at the amputation site is crucial for the successful regeneration of the earthworm and it takes place only in the presence of an intact clitellum. Besides we have performed clitellar transcriptome analysis of the earthworm Eudrilus eugeniae to monitor the key differentially expressed genes and their associated functions and pathways controlling the clitellar tissue changes during both anterior and posterior regeneration of the earthworm. A total of 4707 differentially expressed genes (DEGs) were identified between the control clitellum and clitellum of anterior regenerated earthworms and 4343 DEGs were detected between the control clitellum and clitellum of posterior regenerated earthworms. The functional enrichment analysis confirmed the genes regulating the muscle mass shape and structure were significantly downregulated and the genes associated with response to starvation and anterior-posterior axis specification were significantly upregulated in the clitellar tissue during both anterior and posterior regeneration of the earthworm. The RNA sequencing data of clitellum and the comparative transcriptomic analysis were helpful to understand the complex regeneration process of the earthworm.

Duplication of spiralian-specific TALE genes and evolution of the blastomere specification mechanism in the bivalve lineage

Background

Despite the conserved pattern of the cell-fate map among spiralians, bivalves display several modified characteristics during their early development, including early specification of the D blastomere by the cytoplasmic content, as well as the distinctive fate of the 2d blastomere. However, it is unclear what changes in gene regulatory mechanisms led to such changes in cell specification patterns. Spiralian-TALE (SPILE) genes are a group of spiralian-specific transcription factors that play a role in specifying blastomere cell fates during early development in limpets. We hypothesised that the expansion of SPILE gene repertoires influenced the evolution of the specification pattern of blastomere cell fates.

Results

We performed a transcriptome analysis of early development in the purplish bifurcate mussel and identified 13 SPILE genes. Phylogenetic analysis of the SPILE gene in molluscs suggested that duplications of SPILE genes occurred in the bivalve lineage. We examined the expression patterns of the SPILE gene in mussels and found that some SPILE genes were expressed in quartet-specific patterns, as observed in limpets. Furthermore, we found that several SPILE genes that had undergone gene duplication were specifically expressed in the D quadrant, C and D quadrants or the 2d blastomere. These expression patterns were distinct from the expression patterns of SPILE in their limpet counterparts.

Conclusions

These results suggest that, in addition to their ancestral role in quartet specification, certain SPILE genes in mussels contribute to the specification of the C and D quadrants. We suggest that the expansion of SPILE genes in the bivalve lineage contributed to the evolution of a unique cell fate specification pattern in bivalves.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13227-021-00181-2.

Annelid genomes: Enchytraeus crypticus, a soil model for the innate (and primed) immune system

Enchytraeids (Annelida) are soil invertebrates with worldwide distribution that have served as ecotoxicology models for over 20 years. We present the first high-quality reference genome of Enchytraeus crypticus, assembled from a combination of Pacific Bioscience single-molecule real-time and Illumina sequencing platforms as a 525.2 Mbp genome (910 gapless scaffolds and 18,452 genes). We highlight isopenicillin, acquired by horizontal gene transfer and conferring antibiotic function. Significant gene family expansions associated with regeneration (long interspersed nuclear elements), the innate immune system (tripartite motif-containing protein) and response to stress (cytochrome P450) were identified. The ACE (Angiotensin-converting enzyme) - a homolog of ACE2, which is involved in the coronavirus SARS-CoV-2 cell entry - is also present in E. crypticus. There is an obvious potential of using E. crypticus as a model to study interactions between regeneration, the innate immune system and aging-dependent decline.

Study on the regulatory mechanism of the earthworm microbial community in vitro and in vivo under cadmium stress

In this paper, cadmium (Cd) stress tests were performed on Eisenia fetida in sterile artificial soil, and its regulatory mechanism between microbial communities in vivo and in vitro after Cd stress was explored. In the test, 0, 50, 100, 125, 250 and 500 mg kg^-1 Cd stress concentrations were implemented. After long-term and short-term stress, the microbes in the earthworms and the soil were cultured with ECO plates. The data statistics of carbon source utilization intensity were carried out using the method developed by our team. CCA was scientifically integrated into TOPSIS to establish a new data analysis model to find the regulatory nodes after stress (Ning et al., 2020). Macro gene sequencing technology revealed that the species with the highest absolute abundance in the microbial communities in vivo and in vitro were all unnamed new species. It was confirmed that the HBA gene, NEUROD1 gene and ABCA3 gene were the regulatory genes of the microbial community in the earthworms under Cd stress, while the TC.FEV.OM gene and cheBR gene were the main regulatory genes of the microbial community in the soil. These results provide a scientific and theoretical reference and model basis for the bioremediation of Cd-contaminated soil and the detoxification mechanism of earthworms.

Amynthas corticis genome reveals molecular mechanisms behind global distribution

Earthworms (Annelida: Crassiclitellata) are widely distributed around the world due to their ancient origination as well as adaptation and invasion after introduction into new habitats over the past few centuries. Herein, we report a 1.2 Gb complete genome assembly of the earthworm Amynthas corticis based on a strategy combining third-generation long-read sequencing and Hi-C mapping. A total of 29,256 protein-coding genes are annotated in this genome. Analysis of resequencing data indicates that this earthworm is a triploid species. Furthermore, gene family evolution analysis shows that comprehensive expansion of gene families in the Amynthas corticis genome has produced more defensive functions compared with other species in Annelida. Quantitative proteomic iTRAQ analysis shows that expression of 147 proteins changed in the body of Amynthas corticis and 16 S rDNA sequencing shows that abundance of 28 microorganisms changed in the gut of Amynthas corticis when the earthworm was incubated with pathogenic Escherichia coli O157:H7. Our genome assembly provides abundant and valuable resources for the earthworm research community, serving as a first step toward uncovering the mysteries of this species, and may provide molecular level indicators of its powerful defensive functions, adaptation to complex environments and invasion ability.

Cryptic Clitellata: Molecular Species Delimitation of Clitellate Worms (Annelida): An Overview

Methods for species delimitation using molecular data have developed greatly and have become a staple in systematic studies of clitellate worms. Here we give a historical overview of the data and methods used to delimit clitellates from the mid-1970s to today. We also discuss the taxonomical treatment of the cryptic species, including the recommendation that cryptic species, as far as possible, should be described and named. Finally, we discuss the prospects and further development of the field.

High-quality genome assembly of Metaphire vulgaris

Earthworms enrich the soil and protect the health of their ecological environment. Previous studies on these invertebrates determined their protein content, hormone secretions, medicinal value, and ecological habits, but their whole genomic sequence remains incomplete. We performed whole genome sequencing of Metaphire vulgaris (Chen, 1930), which belongs to the genus Metaphire of the family Megascolecidae. The genome assembly was 729 Mb, with a N50 contig size of 4.2 Mb. In total, 559 contigs were anchored to 41 chromosomes according to the results of Hi-C (High-throughput Chromosome Conformation Capture) technology, which was confirmed by karyological analysis. A comparison of the genomic sequences and genes indicated that there was a whole-genome duplication in M. vulgaris followed by several chromosome fusion events. Hox genes and lumbrokinase genes were identified as partial clusters surrounding the genome. Our high-quality genome assembly of M. vulgaris will provide valuable information for gene function and evolutionary studies in earthworms.

A tale of two leeches: Toward the understanding of the evolution and development of behavioral neural circuits

In the animal kingdom, behavioral traits encompass a broad spectrum of biological phenotypes that have critical roles in adaptive evolution, but an EvoDevo approach has not been broadly used to study behavior evolution. Here, we propose that, by integrating two leech model systems, each of which has already attained some success in its respective field, it is possible to take on behavioral traits with an EvoDevo approach. We first identify the developmental changes that may theoretically lead to behavioral evolution and explain why an EvoDevo study of behavior is challenging. Next, we discuss the pros and cons of the two leech model species, Hirudo, a classic model for invertebrate neurobiology, and Helobdella, an emerging model for clitellate developmental biology, as models for behavioral EvoDevo research. Given the limitations of each leech system, neither is particularly strong for behavioral EvoDevo. However, the two leech systems are complementary in their technical accessibilities, and they do exhibit some behavioral similarities and differences. By studying them in parallel and together with additional leech species such as Haementeria, it is possible to explore the different levels of behavioral development and evolution.

Inferring Tunicate Relationships and the Evolution of the Tunicate Hox Cluster with the Genome of Corella inflata

Tunicates, the closest living relatives of vertebrates, have served as a foundational model of early embryonic development for decades. Comparative studies of tunicate phylogeny and genome evolution provide a critical framework for analyzing chordate diversification and the emergence of vertebrates. Toward this goal, we sequenced the genome of Corella inflata (Ascidiacea, Phlebobranchia), so named for the capacity to brood self-fertilized embryos in a modified, "inflated" atrial chamber. Combining the new genome sequence for Co. inflata with publicly available tunicate data, we estimated a tunicate species phylogeny, reconstructed the ancestral Hox gene cluster at important nodes in the tunicate tree, and compared patterns of gene loss between Co. inflata and Ciona robusta, the prevailing tunicate model species. Our maximum-likelihood and Bayesian trees estimated from a concatenated 210-gene matrix were largely concordant and showed that Aplousobranchia was nested within a paraphyletic Phlebobranchia. We demonstrated that this relationship is not an artifact due to compositional heterogeneity, as had been suggested by previous studies. In addition, within Thaliacea, we recovered Doliolida as sister to the clade containing Salpida and Pyrosomatida. The Co. inflata genome provides increased resolution of the ancestral Hox clusters of key tunicate nodes, therefore expanding our understanding of the evolution of this cluster and its potential impact on tunicate morphological diversity. Our analyses of other gene families revealed that several cardiovascular associated genes (e.g., BMP10, SCL2A12, and PDE2a) absent from Ci. robusta, are present in Co. inflata. Taken together, our results help clarify tunicate relationships and the genomic content of key ancestral nodes within this phylogeny, providing critical insights into tunicate evolution.

Genome and single-cell RNA-sequencing of the earthworm Eisenia andrei identifies cellular mechanisms underlying regeneration

The earthworm is particularly fascinating to biologists because of its strong regenerative capacity. However, many aspects of its regeneration in nature remain elusive. Here we report chromosome-level genome, large-scale transcriptome and single-cell RNA-sequencing data during earthworm (Eisenia andrei) regeneration. We observe expansion of LINE2 transposable elements and gene families functionally related to regeneration (for example, EGFR, epidermal growth factor receptor) particularly for genes exhibiting differential expression during earthworm regeneration. Temporal gene expression trajectories identify transcriptional regulatory factors that are potentially crucial for initiating cell proliferation and differentiation during regeneration. Furthermore, early growth response genes related to regeneration are transcriptionally activated in both the earthworm and planarian. Meanwhile, single-cell RNA-sequencing provides insight into the regenerative process at a cellular level and finds that the largest proportion of cells present during regeneration are stem cells.

The mechanisms regulating regeneration of the earthworm are unclear. Here, the authors use genomic and transcriptomic analysis of the earthworm Eisenia andrei together with Hi-C analysis to identify genes involved and show activation of LINE2 transposable elements on regeneration.

The draft genome of a new Verminephrobacter eiseniae strain: a nephridial symbiont of earthworms

Purpose

Verminephrobacter is a genus of symbiotic bacteria that live in the nephridia of earthworms. The bacteria are recruited during the embryonic stage of the worm and transferred from generation to generation in the same manner. The worm provides shelter and food for the bacteria. The bacteria deliver micronutrients to the worm. The present study reports the genome sequence assembly and annotation of a new strain of Verminephrobacter called Verminephrobacter eiseniae msu.

Methods

We separated the sequences of a new Verminephrobacter strain from the whole genome of Eisenia fetida using the sequence of V. eiseniae EF01-2, and the bacterial genome was assembled using the CLC Workbench. The de novo-assembled genome was annotated and analyzed for the protein domains, functions, and metabolic pathways. Besides, the multigenome comparison was performed to interpret the phylogenomic relationship of the strain with other proteobacteria.

Result

The FastqSifter sifted a total of 593,130 Verminephrobacter genomic reads. The de novo assembly of the reads generated 1832 contigs with a total genome size of 4.4 Mb. The Average Nucleotide Identity denoted the bacterium belongs to the species V. eiseniae, and the 16S rRNA analysis confirmed it as a new strain of V. eiseniae. The AUGUSTUS genome annotation predicted a total of 3809 protein-coding genes; of them, 3805 genes were identified from the homology search.

Conclusion

The bioinformatics analysis confirmed the bacterium is an isolate of V. eiseniae, and it was named Verminephrobacter eiseniae msu. The whole genome of the bacteria can be utilized as a useful resource to explore the area of symbiosis further.

Exploration of dipeptidyl-peptidase IV (DPP IV) inhibitors in a low-molecular mass extract of the earthworm Eisenia fetida and identification of the inhibitors as amino acids like methionine, leucine, histidine, and isoleucine

We have reported previously that the water extract of the earthworm Eisenia fetida has inhibitory effect on human dipeptidyl-peptidase IV (DPP IV) in vitro. Here we studied to identify DPP IV inhibitors in a low-molecular mass extract (designated U3EE) under 3 kDa prepared from the water extract. U3EE showed 50 % inhibition (IC₅₀) at the concentration of 5.3 ± 0.3 mg/mL. An inhibitory active fraction obtained by solid-phase extraction of U3EE was separated into three parts by reversed-phase HPLC. These parts were shown by GC/MS to be composed of ten (Ala, Gly, Thr, Ser, Asn, Asp, Lys, His, Orn, and cystine), two (Leu and Ile), and one (Met) amino acids, respectively. Among them, Met, Leu, and His showed strong inhibition with IC₅₀ values of 3.4 ± 0.3, 6.1 ± 0.3 and 14.7 ± 1.2 mM, respectively; Ala, Lys, Orn, and Ile showed rather weaker inhibition than those, while the others showed no inhibition. Met, Leu, and Ile were competitive inhibitors and His was a mixed-type one. DPP IV inhibition by U3EE might be due to additive and/or synergistic effects of the inhibitory amino acids, suggesting that it could be useful as pharmaceutical and supplement for diabetes prevention.

Hox gene expression during development of the phoronid Phoronopsis harmeri

Background

Phoronida is a small group of marine worm-like suspension feeders, which together with brachiopods and bryozoans form the clade Lophophorata. Although their development is well studied on the morphological level, data regarding gene expression during this process are scarce and restricted to the analysis of relatively few transcription factors. Here, we present a description of the expression patterns of Hox genes during the embryonic and larval development of the phoronid Phoronopsis harmeri.

Results

We identified sequences of eight Hox genes in the transcriptome of Ph. harmeri and determined their expression pattern during embryonic and larval development using whole mount in situ hybridization. We found that none of the Hox genes is expressed during embryonic development. Instead their expression is initiated in the later developmental stages, when the larval body is already formed. In the investigated initial larval stages the Hox genes are expressed in the non-collinear manner in the posterior body of the larvae: in the telotroch and the structures that represent rudiments of the adult worm. Additionally, we found that certain head-specific transcription factors are expressed in the oral hood, apical organ, preoral coelom, digestive system and developing larval tentacles, anterior to the Hox-expressing territories.

Conclusions

The lack of Hox gene expression during early development of Ph. harmeri indicates that the larval body develops without positional information from the Hox patterning system. Such phenomenon might be a consequence of the evolutionary intercalation of the larval form into an ancestral life cycle of phoronids. The observed Hox gene expression can also be a consequence of the actinotrocha representing a “head larva”, which is composed of the most anterior body region that is devoid of Hox gene expression. Such interpretation is further supported by the expression of head-specific transcription factors. This implies that the Hox patterning system is used for the positional information of the trunk rudiments and is, therefore, delayed to the later larval stages. We propose that a new body form was intercalated to the phoronid life cycle by precocious development of the anterior structures or by delayed development of the trunk rudiment in the ancestral phoronid larva.

Draft Genome of the Asian Buffalo Leech Hirudinaria manillensis

The Asian Buffalo leech, Hirudinaria manillensis, is an aquatic sanguivorous species distributed widely in Southeast Asia. H. manillensis has long been used clinically for bloodletting and other medical purposes. Recent studies have focused on artificial culturing, strain optimization, and the identification and development new drugs based on the anticoagulant effects of H. manillensis bites; however, data regarding its genome remain unclear. This study aimed to determine the genome sequence of an adult Asian Buffalo leech. We generated a draft assembly of 151.8 Mb and a N50 scaffold of 2.28 Mb. Predictions indicated that the assembled genome contained 21,005 protein-coding genes. Up to 17,865 genes were annotated in multiple databases including Gene Ontology. Sixteen anticoagulant proteins with a Hirudin or Antistasin domain were identified. This study is the first to report the whole-genome sequence of the Asian Buffalo leech, an important sanguivorous leech of clinical significance. The quality of the assembly is comparable to those of other annelids. These data will help further the current understanding of the biological mechanisms and genetic characteristics of leeches and serve as a valuable resource for future studies.

Oppia nitens C.L. Koch, 1836 (Acari: Oribatida): Current Status of Its Bionomics and Relevance as a Model Invertebrate in Soil Ecotoxicology

The oribatid soil mite Oppia nitens C.L. Koch, 1836, is a model microarthropod in soil ecotoxicity testing. This species has a significant role in supporting soil functions and as a suitable indicator of soil contamination. Despite its significance to the environment and to ecotoxicology, however, very little is known of its biology, ecology, and suborganismal responses to contaminants in the soil. In the present review, we present detailed and critical insights into the biology and ecology of O. nitens in relation to traits that are crucial to its adaptive responses to contaminants in soil. We used a species sensitivity distribution model to rank the species sensitivity to heavy metals (cadmium and zinc) and neonicotinoids (imidacloprid and thiacloprid) compared with other standardized soil invertebrates. Although the International Organization for Standardization and Environment and Climate Change Canada are currently standardizing a protocol for the use of O. nitens in soil toxicity testing, we believe that O. nitens is limited as a model soil invertebrate until the molecular pathways associated with its response to contaminants are better understood. These pathways can only be elucidated with information from the mites' genome or transcriptome, which is currently lacking. Despite this limitation, we propose a possible molecular pathway to metal tolerance and a putative adverse outcome pathway to heavy metal toxicity in O. nitens. Environ Toxicol Chem 2019;38:2593-2613. © 2019 SETAC.

Functional role of pax6 during eye and nervous system development in the annelid Capitella teleta

The transcription factor Pax6 is an important regulator of early animal development. Loss of function mutations of pax6 in a range of animals result in a reduction or complete loss of the eye, a reduction of a subset of neurons, and defects in axon growth. There are no studies focusing on the role of pax6 during development of any lophotrochozoan representative, however, expression of pax6 in the developing eye and nervous system in a number of species suggest that pax6 plays a highly conserved role in eye and nervous system formation. We investigated the functional role of pax6 during development of the marine annelid Capitella teleta. Expression of pax6 transcripts in C. teleta larvae is similar to patterns found in other animals, with distinct subdomains in the brain and ventral nerve cord as well as in the larval and juvenile eye. To perturb pax6 function, two different splice-blocking morpholinos and a translation-blocking morpholino were used. Larvae resulting from microinjections with either splice-blocking morpholino show a reduction of the pax6 transcript. Development of both the larval eyes and the central nervous system architecture are highly disrupted following microinjection of each of the three morpholinos. The less severe phenotype observed when only the homeodomain is disrupted suggests that presence of the paired domain is sufficient for partial function of the Pax6 protein. Preliminary downstream target analysis confirms disruption in expression of some components of the retinal gene regulatory network, as well as disruption of genes involved in nervous system development. Results from this study, taken together with studies from other species, reveal an evolutionarily conserved role for pax6 in eye and neural specification and development.

Sponges Lack ParaHox Genes

Addressing the origin of axial-patterning machinery is essential for understanding the evolution of animal form. Historically, sponges, a lineage that branched off early in animal evolution, were thought to lack Hox and ParaHox genes, suggesting that these critical axial-patterning genes arose after sponges diverged. However, a recent study has challenged this long-held doctrine by claiming to identify ParaHox genes (Cdx family) in two calcareous sponge species, Sycon ciliatum and Leucosolenia complicata. We reanalyzed the main data sets in this paper and analyzed an additional data set that expanded the number of bilaterians represented and removed outgroup homeodomains. As in the previous study, our Neighbor-Joining analyses of the original data sets recovered a clade that included sponge and Cdx genes, whereas Bayesian analyses placed these sponge genes within the NKL subclass of homeodomains. Unlike the original study, only one of our two maximum-likelihood analyses was congruent with Cdx genes in sponges. Our analyses of our additional data set led to the sponge genes consistently being placed within the NKL subclass of homeodomains regardless of method or model. Our results show more support for these sponge genes belonging to the NKL subclass, and therefore imply that Hox and ParaHox genes arose after Porifera diverged from the rest of animals.

Data on genome annotation and analysis of earthworm Eisenia fetida

The present article reports the complete draft genome annotation of earthworm Eisenia fetida, obtained from the manuscript entitled "Timing and Scope of Genomic Expansion within Annelida: Evidence from Homeoboxes in the Genome of the Earthworm E. fetida" (Zwarycz et al., 2015) and provides the data on the repetitive elements, protein coding genes and noncoding RNAs present in the genome dataset of the species. The E. fetida protein coding genes were predicted from AUGUSTUS gene prediction and subsequently annotated based on their sequence similarity, Gene Ontology (GO) functional terms, InterPro domains, Clusters of Orthologous Groups (COGs) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways information. The genome wide comparison of orthologous clusters and phylogenomic analysis of the core genes were performed to understand the events of genome evolution and genomic diversity between E. fetida and its related metazoans. In addition, the genome dataset was screened to identify the crucial stem cell markers, regeneration specific genes and immune-related genes and their functionally enriched GO terms were predicted from Fisher׳s enrichment analysis. The E. fetida genome annotation data containing the GFF (general feature format) annotation file, predicted coding gene sequences and translated protein sequences were deposited to the figshare repository under the DOI: https://doi.org/10.6084/m9.figshare.6142322.v1.

Large scale changes in the transcriptome of Eisenia fetida during regeneration

Earthworms show a wide spectrum of regenerative potential with certain species like Eisenia fetida capable of regenerating more than two-thirds of their body while other closely related species, such as Paranais litoralis seem to have lost this ability. Earthworms belong to the phylum Annelida, in which the genomes of the marine oligochaete Capitella telata and the freshwater leech Helobdella robusta have been sequenced and studied. Herein, we report the transcriptomic changes in Eisenia fetida (Indian isolate) during regeneration. Following injury, E. fetida regenerates the posterior segments in a time spanning several weeks. We analyzed gene expression changes both in the newly regenerating cells and in the adjacent tissue, at early (15days post amputation), intermediate (20days post amputation) and late (30 days post amputation) by RNAseq based de novo assembly and comparison of transcriptomes. We also generated a draft genome sequence of this terrestrial red worm using short reads and mate-pair reads. An in-depth analysis of the miRNome of the worm showed that many miRNA gene families have undergone extensive duplications. Sox4, a master regulator of TGF-beta mediated epithelial-mesenchymal transition was induced in the newly regenerated tissue. Genes for several proteins such as sialidases and neurotrophins were identified amongst the differentially expressed transcripts. The regeneration of the ventral nerve cord was also accompanied by the induction of nerve growth factor and neurofilament genes. We identified 315 novel differentially expressed transcripts in the transcriptome, that have no homolog in any other species. Surprisingly, 82% of these novel differentially expressed transcripts showed poor potential for coding proteins, suggesting that novel ncRNAs may play a critical role in regeneration of earthworm.

A Revised Spiralian Homeobox Gene Classification Incorporating New Polychaete Transcriptomes Reveals a Diverse TALE Class and a Divergent Hox Gene

The diversity of mechanisms and capacity for regeneration across the Metazoa present an intriguing challenge in evolutionary biology, impacting on the burgeoning field of regenerative medicine. Broad taxonomic sampling is essential to improve our understanding of regeneration, and studies outside of the traditional model organisms have proved extremely informative. Within the historically understudied Spiralia, the Annelida have an impressive variety of tractable regenerative systems. The biomeralizing, blastema-less regeneration of the head appendage (operculum) of the serpulid polychaete keelworm Spirobranchus (formerly Pomatoceros) lamarcki is one such system. To profile potential regulatory mechanisms, we classified the homeobox gene content of opercular regeneration transcriptomes. As a result of retrieving several difficult-to-classify homeobox sequences, we performed an extensive search and phylogenetic analysis of the TALE and PRD-class homeobox gene content of a broad selection of lophotrochozoan genomes. These analyses contribute to our increasing understanding of the diversity, taxonomic extent, rapid evolution, and radical flexibility of these recently discovered homeobox gene radiations. Our expansion and integration of previous nomenclature systems helps to clarify their cryptic orthology. We also describe an unusual divergent S. lamarcki Antp gene, a previously unclassified lophotrochozoan orphan gene family (Lopx), and a number of novel Nk class orphan genes. The expression and potential involvement of many of these lineage- and clade-restricted homeobox genes in S. lamarcki operculum regeneration provides an example of diversity in regenerative mechanisms, as well as significantly improving our understanding of homeobox gene evolution.

Annotation of nerve cord transcriptome in earthworm Eisenia fetida

In annelid worms, the nerve cord serves as a crucial organ to control the sensory and behavioral physiology. The inadequate genome resource of earthworms has prioritized the comprehensive analysis of their transcriptome dataset to monitor the genes express in the nerve cord and predict their role in the neurotransmission and sensory perception of the species. The present study focuses on identifying the potential transcripts and predicting their functional features by annotating the transcriptome dataset of nerve cord tissues prepared by Gong et al., 2010 from the earthworm Eisenia fetida. Totally 9762 transcripts were successfully annotated against the NCBI nr database using the BLASTX algorithm and among them 7680 transcripts were assigned to a total of 44,354 GO terms. The conserve domain analysis indicated the over representation of P-loop NTPase domain and calcium binding EF-hand domain. The COG functional annotation classified 5860 transcript sequences into 25 functional categories. Further, 4502 contig sequences were found to map with 124 KEGG pathways. The annotated contig dataset exhibited 22 crucial neuropeptides having considerable matches to the marine annelid Platynereis dumerilii, suggesting their possible role in neurotransmission and neuromodulation. In addition, 108 human stem cell marker homologs were identified including the crucial epigenetic regulators, transcriptional repressors and cell cycle regulators, which may contribute to the neuronal and segmental regeneration. The complete functional annotation of this nerve cord transcriptome can be further utilized to interpret genetic and molecular mechanisms associated with neuronal development, nervous system regeneration and nerve cord function.

Expansion of TALE homeobox genes and the evolution of spiralian development

Spiralians, including molluscs, annelids and platyhelminths, share a unique development process that includes the typical geometry of early cleavage and early segregation of cell fate in blastomeres along the animal-vegetal axis. However, the molecular mechanisms underlying this early cell fate segregation are largely unknown. Here, we report spiralian-specific expansion of the three-amino-acid loop extension (TALE) class of homeobox genes. During early development, some of these TALE genes are expressed in staggered domains along the animal-vegetal axis in the limpet Nipponacmea fuscoviridis and the polychaete Spirobranchus kraussii. Inhibition or overexpression of these genes alters the developmental fate of blastomeres, as predicted by the gene expression patterns. These results suggest that the expansion of novel TALE genes plays a critical role in the establishment of a novel cell fate segregation mechanism in spiralians.

Four paralog gelsolin genes are differentially expressed in the earthworm Lumbricus terrestris

We have identified and characterized four distinct variants of the gelsolin-related protein (EWAM P1-P4) in the earthworm L. terrestris. All of these proteins biochemically qualify as gelsolins since they sever actin filaments in a calcium dependent manner. P1, P2 and P3 are present in the Lumbricus body wall muscle whereas in the gizzard muscle P3 and P4 were found. P1-P4 are encoded by four paralog genes and are differentially expressed in various muscle cell tissues. While the genes for P1 and P2 contain one intron, there was no intron in both P3 and P4 genes. The coding sequences consist of 1104bp (368 amino acids) for P1/P4 and 1101bp (367 amino acids) for P2/P3. Corresponding genes were confirmed by northern blot analysis which revealed three (calculated lengths: 3100, 2300 and 2100 nucleotides) and two (calculated lengths: 2300 and 1700 nucleotides) mRNA transcripts in the body wall and the gizzard, respectively. EWAM mRNA was localized by fluorescence in situ hybridization in the body wall and the gizzard muscle. P1 mRNA was detected in the inner proximal layers of both the circular and longitudinal muscle of the body wall whereas in the gizzard no significant staining was observed for P1. P2-P4 mRNAs were abundant in the outer distal layers of both the circular and the longitudinal muscles of both body wall and gizzard. The differential expression of four paralog gelsolin genes suggests a functional adaptation of different muscle cells with respect to actin filament turnover and modulation of its polymer state.

Clustered brachiopod Hox genes are not expressed collinearly and are associated with lophotrochozoan novelties

Temporal collinearity is often considered the main force preserving Hox gene clusters in animal genomes. Studies that combine genomic and gene expression data are scarce, however, particularly in invertebrates like the Lophotrochozoa. As a result, the temporal collinearity hypothesis is currently built on poorly supported foundations. Here we characterize the complement, cluster, and expression of Hox genes in two brachiopod species, Terebratalia transversa and Novocrania anomalaT. transversa has a split cluster with 10 genes (lab, pb, Hox3, Dfd, Scr, Lox5, Antp, Lox4, Post2, and Post1), whereas N. anomala has 9 genes (apparently missing Post1). Our in situ hybridization, real-time quantitative PCR, and stage-specific transcriptomic analyses show that brachiopod Hox genes are neither strictly temporally nor spatially collinear; only pb (in T. transversa), Hox3 (in both brachiopods), and Dfd (in both brachiopods) show staggered mesodermal expression. Thus, our findings support the idea that temporal collinearity might contribute to keeping Hox genes clustered. Remarkably, expression of the Hox genes in both brachiopod species demonstrates cooption of Hox genes in the chaetae and shell fields, two major lophotrochozoan morphological novelties. The shared and specific expression of Hox genes, together with Arx, Zic, and Notch pathway components in chaetae and shell fields in brachiopods, mollusks, and annelids provide molecular evidence supporting the conservation of the molecular basis for these lophotrochozoan hallmarks.

Embryonic expression patterns of Hox genes in the oligochaete annelid Tubifex tubifex

We have cloned and characterized the expression of seven Hox genes (designated Ttu-lab, Ttu-Dfd, Ttu-Scr1, Ttu-Scr2, Ttu-Lox5, Ttu-Lox4 and Ttu-Lox2) from the oligochaete annelid Tubifex tubifex. RT-PCR analyses show that except for Ttu-Lox4 and Ttu-Lox2 which begin expression as early as cleavage stages, Tubifex Hox genes are expressed during stages 13-18 when embryos undergo germ band formation, segmentation and body elongation. In terms of combination of tissues (or organs) exhibiting positive cells, the Tubifex Hox genes examined in this study are classified into three groups. Ttu-lab, Ttu-Scr1 and Ttu-Lox5 are expressed only in the ventral nerve cord; Ttu-Scr2 and Ttu-Lox4 are expressed not only in the ventral nerve cord but also in distinct lateral segmental tissues; and Ttu-Dfd and Ttu-Lox2 are expressed not only in the segmental ectoderm along the length of the AP body axis but also in the prostomium. Anterior expression boundaries of Ttu-lab, Ttu-Scr1, Ttu-Lox5 and Ttu-Lox4 are at segments 3, 4, 5, and 9, respectively. Anterior expression boundary of Ttu-Scr2 is at segment 2, and Ttu-Dfd and Ttu-Lox2 are expressed even at the anteriormost portion, the prostomium. These observations suggest that as in other annelids, so-called "spatial colinearity" of anterior expression boundaries of Hox genes has been conserved in the oligochaetes. It is also evident that there are some oligochaete Hox genes which violate the spatial colinearity rule.

Low coverage sequencing of three echinoderm genomes: the brittle star Ophionereis fasciata, the sea star Patiriella regularis, and the sea cucumber Australostichopus mollis

BACKGROUND

There are five major extant groups of Echinodermata: Crinoidea (feather stars and sea lillies), Ophiuroidea (brittle stars and basket stars), Asteroidea (sea stars), Echinoidea (sea urchins, sea biscuits, and sand dollars), and Holothuroidea (sea cucumbers). These animals are known for their pentaradial symmetry as adults, unique water vascular system, mutable collagenous tissues, and endoskeletons of high magnesium calcite. To our knowledge, the only echinoderm species with a genome sequence available to date is Strongylocentrotus pupuratus (Echinoidea). The availability of additional echinoderm genome sequences is crucial for understanding the biology of these animals.

FINDINGS

Here we present assembled draft genomes of the brittle star Ophionereis fasciata, the sea star Patiriella regularis, and the sea cucumber Australostichopus mollis from Illumina sequence data with coverages of 12.5x, 22.5x, and 21.4x, respectively.

CONCLUSIONS

These data provide a resource for mining gene superfamilies, identifying non-coding RNAs, confirming gene losses, and designing experimental constructs. They will be important comparative resources for future genomic studies in echinoderms.

An Overview of Hox Genes in Lophotrochozoa: Evolution and Functionality

Hox genes are regulators of animal embryonic development. Changes in the number and sequence of Hox genes as well as in their expression patterns have been related to the evolution of the body plan. Lophotrochozoa is a clade of Protostomia characterized by several phyla which show a wide morphological diversity. Despite that the works summarized in this review emphasize the fragmentary nature of the data available regarding the presence and expression of Hox genes, they also offer interesting insight into the evolution of the Hox cluster and the role played by Hox genes in several phyla. However, the number of genes involved in the cluster of the lophotrochozoan ancestor is still a question of debate. The data presented here suggest that at least nine genes were present while two other genes, Lox4 and Post-2, may either have been present in the ancestor or may have arisen as a result of duplication in the Brachiopoda-Mollusca-Annelida lineage. Spatial and temporal collinearity is a feature of Hox gene expression which was probably present in the ancestor of deuterostomes and protostomes. However, in Lophotrochozoa, it has been detected in only a few species belonging to Annelida and Mollusca.