Insights into bilaterian evolution from three spiralian genomes

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.

Draft genome of Semisulcospira libertina, a species of freshwater snail

Semisulcospira libertina, a species of freshwater snail, is widespread in East Asia. It is important as a food source. Additionally, it is a vector of clonorchiasis, paragonimiasis, metagonimiasis, and other parasites. Although S. libertina has ecological, commercial, and clinical importance, its whole-genome has not been reported yet. Here, we revealed the genome of S. libertina through de novo assembly. We assembled the whole-genome of S. libertina and determined its transcriptome for the first time using Illumina NovaSeq 6000 platform. According to the k-mer analysis, the genome size of S. libertina was estimated to be 3.04 Gb. Using RepeatMasker, a total of 53.68% of repeats were identified in the genome assembly. Genome data of S. libertina reported in this study will be useful for identification and conservation of S. libertina in East Asia.

Violet bioluminescent Polycirrus sp. (Annelida: Terebelliformia) discovered in the shallow coastal waters of the Noto Peninsula in Japan

Terebellidae worms have large numbers of tentacles responsible for various biological functions. Some Terebellidae worms whose tentacles emit light are found around the world, including exceptional violet-light-emitting Polycirrus spp. found in Europe and North America. However, there is no video-recorded observation of the luminous behavior of such unique species in nature, and the genetic information related to their ecology are lacking. Here, for the first time, we video-recorded the violet-light-emitting behavior of an undescribed Japanese worm in its natural habitat. The worm was designated as Polycirrus sp. ISK based on morphological observations, and the luminescence spectrum showed a peak at 444 nm, which is an exceptionally short wavelength for bioluminescence in a shallow coastal water environment. An analysis of differentially expressing genes based on separate RNA-Seq analysis for the tentacles and the rest of body revealed the specific expression of genes that are probably involved in innate immunity in the tentacles exposed to predators. We also found a Renilla luciferase homologous gene, but coelenterazine was not detected in the worm extract by analyses using a liquid chromatography and a recombinant Renilla luciferase. These results will promote an understanding of the ecology and luminescence mechanisms of luminous Polycirrus spp.

The Nereid on the rise: Platynereis as a model system

The Nereid Platynereis dumerilii (Audouin and Milne Edwards (Annales des Sciences Naturelles 1:195-269, 1833) is a marine annelid that belongs to the Nereididae, a family of errant polychaete worms. The Nereid shows a pelago-benthic life cycle: as a general characteristic for the superphylum of Lophotrochozoa/Spiralia, it has spirally cleaving embryos developing into swimming trochophore larvae. The larvae then metamorphose into benthic worms living in self-spun tubes on macroalgae. Platynereis is used as a model for genetics, regeneration, reproduction biology, development, evolution, chronobiology, neurobiology, ecology, ecotoxicology, and most recently also for connectomics and single-cell genomics. Research on the Nereid started with studies on eye development and spiralian embryogenesis in the nineteenth and early twentieth centuries. Transitioning into the molecular era, Platynereis research focused on posterior growth and regeneration, neuroendocrinology, circadian and lunar cycles, fertilization, and oocyte maturation. Other work covered segmentation, photoreceptors and other sensory cells, nephridia, and population dynamics. Most recently, the unique advantages of the Nereid young worm for whole-body volume electron microscopy and single-cell sequencing became apparent, enabling the tracing of all neurons in its rope-ladder-like central nervous system, and the construction of multimodal cellular atlases. Here, we provide an overview of current topics and methodologies for P. dumerilii, with the aim of stimulating further interest into our unique model and expanding the active and vibrant Platynereis community.

The potential for using shell proteins in gastropod systematics, assessed in patellogastropod limpets

This investigation of the application of shell protein information to gastropod systematics initially utilized available Lottia gigantea sequences and a transcriptome of Patelloida mimula developed here. Levels of differentiation between predicted sequences of reciprocal best-hit potential homologues in P. mimula and L. gigantea suggested that they could be useful within families, and possibly in higher taxa using some shell-associated proteins, particularly the peroxidases. Subsequently, proteomic analyses of the acid-soluble fraction of extractions from 17 shells and five tissue samples were conducted by combined liquid chromatography/mass spectrometry with nano-electrospray ionization. All proteins with abundance more than 1.2% in the L. gigantea shell proteome were identified with 100% confidence in most extractions by SearchGui/PeptideShaker analyses. In total, 259 of 379 peptides predicted from in silico digestion of L. gigantea shell proteins were represented by validated peptide spectrum matches in one or more specimens. Systematics applications were investigated by analysing metrics such as protein coverage by peptides and phylogenetic analyses of peptide presence/absence. The investigation suggested that diagnostic profiles based on fixed presence/absence differences can be used to separate species pairs. However, further development of analytical techniques and accumulation of reference databases is required for realising fully the systematics potential of the shell proteome.

A highly-contiguous genome assembly of the Eurasian spruce bark beetle, Ips typographus, provides insight into a major forest pest

Conifer-feeding bark beetles are important herbivores and decomposers in forest ecosystems. These species complete their life cycle in nutritionally poor substrates and some can kill enormous numbers of trees during population outbreaks. The Eurasian spruce bark beetle (Ips typographus) can destroy >100 million m3 of spruce in a single year. We report a 236.8 Mb I. typographus genome assembly using PacBio long-read sequencing. The final phased assembly has a contig N50 of 6.65 Mb in 272 contigs and is predicted to contain 23,923 protein-coding genes. We reveal expanded gene families associated with plant cell wall degradation, including pectinases, aspartyl proteases, and glycosyl hydrolases. This genome sequence from the genus Ips provides timely resources to address questions about the evolutionary biology of the true weevils (Curculionidae), one of the most species-rich animal families. In forests of today, increasingly stressed by global warming, this draft genome may assist in developing pest control strategies to mitigate outbreaks.

TIAMMAt: Leveraging biodiversity to revise protein domain models, evidence from innate immunity

Sequence annotation is fundamental for studying the evolution of protein families, particularly when working with nonmodel species. Given the rapid, ever-increasing number of species receiving high-quality genome sequencing, accurate domain modeling that is representative of species diversity is crucial for understanding protein family sequence evolution and their inferred function(s). Here, we describe a bioinformatic tool called Taxon-Informed Adjustment of Markov Model Attributes (TIAMMAt) which revises domain profile hidden Markov models (HMMs) by incorporating homologous domain sequences from underrepresented and nonmodel species. Using innate immunity pathways as a case study, we show that revising profile HMM parameters to directly account for variation in homologs among underrepresented species provides valuable insight into the evolution of protein families. Following adjustment by TIAMMAt, domain profile HMMs exhibit changes in their per-site amino acid state emission probabilities and insertion/deletion probabilities while maintaining the overall structure of the consensus sequence. Our results show that domain revision can heavily impact evolutionary interpretations for some families (i.e., NLR’s NACHT domain), whereas impact on other domains (e.g., rel homology domain and interferon regulatory factor domains) is minimal due to high levels of sequence conservation across the sampled phylogenetic depth (i.e., Metazoa). Importantly, TIAMMAt revises target domain models to reflect homologous sequence variation using the taxonomic distribution under consideration by the user. TIAMMAt’s flexibility to revise any subset of the Pfam database using a user-defined taxonomic pool will make it a valuable tool for future protein evolution studies, particularly when incorporating (or focusing) on nonmodel species.

The Evolution and Diversity of Interleukin-17 Highlight an Expansion in Marine Invertebrates and Its Conserved Role in Mucosal Immunity

The interleukin-17 (IL-17) family consists of proinflammatory cytokines conserved during evolution. A comparative genomics approach was applied to examine IL-17 throughout evolution from poriferans to higher vertebrates. Cnidaria was highlighted as the most ancient diverged phylum, and several evolutionary patterns were revealed. Large expansions of the IL-17 repertoire were observed in marine molluscs and echinoderm species. We further studied this expansion in filter-fed Mytilus galloprovincialis, which is a bivalve with a highly effective innate immune system supported by a variable pangenome. We recovered 379 unique IL-17 sequences and 96 receptors from individual genomes that were classified into 23 and 6 isoforms after phylogenetic analyses. Mussel IL-17 isoforms were conserved among individuals and shared between closely related Mytilidae species. Certain isoforms were specifically implicated in the response to a waterborne infection with Vibrio splendidus in mussel gills. The involvement of IL-17 in mucosal immune responses could be conserved in higher vertebrates from these ancestral lineages.

Reconstruction of proto-vertebrate, proto-cyclostome and proto-gnathostome genomes provides new insights into early vertebrate evolution

Ancient polyploidization events have had a lasting impact on vertebrate genome structure, organization and function. Some key questions regarding the number of ancient polyploidization events and their timing in relation to the cyclostome-gnathostome divergence have remained contentious. Here we generate de novo long-read-based chromosome-scale genome assemblies for the Japanese lamprey and elephant shark. Using these and other representative genomes and developing algorithms for the probabilistic macrosynteny model, we reconstruct high-resolution proto-vertebrate, proto-cyclostome and proto-gnathostome genomes. Our reconstructions resolve key questions regarding the early evolutionary history of vertebrates. First, cyclostomes diverged from the lineage leading to gnathostomes after a shared tetraploidization (1R) but before a gnathostome-specific tetraploidization (2R). Second, the cyclostome lineage experienced an additional hexaploidization. Third, 2R in the gnathostome lineage was an allotetraploidization event, and biased gene loss from one of the subgenomes shaped the gnathostome genome by giving rise to remarkably conserved microchromosomes. Thus, our reconstructions reveal the major evolutionary events and offer new insights into the origin and evolution of vertebrate genomes.

Genomic signatures supporting the symbiosis and formation of chitinous tube in the deep-sea tubeworm Paraescarpia echinospica

Vestimentiferan tubeworms are iconic animals that present as large habitat-forming chitinized tube bushes in deep-sea chemosynthetic ecosystems. They are gutless and depend entirely on their endosymbiotic sulfide-oxidizing chemoautotrophic bacteria for nutrition. Information on the genomes of several siboglinid endosymbionts has improved our understanding of their nutritional supplies. However, the interactions between tubeworms and their endosymbionts remain largely unclear due to a paucity of host genomes. Here, we report the chromosome-level genome of the vestimentiferan tubeworm Paraescarpia echinospica. We found that the genome has been remodeled to facilitate symbiosis through the expansion of gene families related to substrate transfer and innate immunity, suppression of apoptosis, regulation of lysosomal digestion, and protection against oxidative stress. Furthermore, the genome encodes a programmed cell death pathway that potentially controls the endosymbiont population. Our integrated genomic, transcriptomic, and proteomic analyses uncovered matrix proteins required for the formation of the chitinous tube and revealed gene family expansion and co-option as evolutionary mechanisms driving the acquisition of this unique supporting structure for deep-sea tubeworms. Overall, our study provides novel insights into the host’s support system that has enabled tubeworms to establish symbiosis, thrive in deep-sea hot vents and cold seeps, and produce the unique chitinous tubes in the deep sea.

Sawfly Genomes Reveal Evolutionary Acquisitions That Fostered the Mega-Radiation of Parasitoid and Eusocial Hymenoptera

The tremendous diversity of Hymenoptera is commonly attributed to the evolution of parasitoidism in the last common ancestor of parasitoid sawflies (Orussidae) and wasp-waisted Hymenoptera (Apocrita). However, Apocrita and Orussidae differ dramatically in their species richness, indicating that the diversification of Apocrita was promoted by additional traits. These traits have remained elusive due to a paucity of sawfly genome sequences, in particular those of parasitoid sawflies. Here, we present comparative analyses of draft genomes of the primarily phytophagous sawfly Athalia rosae and the parasitoid sawfly Orussus abietinus. Our analyses revealed that the ancestral hymenopteran genome exhibited traits that were previously considered unique to eusocial Apocrita (e.g., low transposable element content and activity) and a wider gene repertoire than previously thought (e.g., genes for CO₂ detection). Moreover, we discovered that Apocrita evolved a significantly larger array of odorant receptors than sawflies, which could be relevant to the remarkable diversification of Apocrita by enabling efficient detection and reliable identification of hosts.

Evolutionary History of DNA Methylation Related Genes in Bivalvia: New Insights From Mytilus galloprovincialis

DNA methylation is an essential epigenetic mechanism influencing gene expression in all organisms. In metazoans, the pattern of DNA methylation changes during embryogenesis and adult life. Consequently, differentiated cells develop a stable and unique DNA methylation pattern that finely regulates mRNA transcription during development and determines tissue-specific gene expression. Currently, DNA methylation remains poorly investigated in mollusks and completely unexplored in Mytilus galloprovincialis. To shed light on this process in this ecologically and economically important bivalve, we screened its genome, detecting sequences homologous to DNA methyltransferases (DNMTs), methyl-CpG-binding domain (MBD) proteins and Ten-eleven translocation methylcytosine dioxygenase (TET) previously described in other organisms. We characterized the gene architecture and protein domains of the mussel sequences and studied their phylogenetic relationships with the ortholog sequences from other bivalve species. We then comparatively investigated their expression levels across different adult tissues in mussel and other bivalves, using previously published transcriptome datasets. This study provides the first insights on DNA methylation regulators in M. galloprovincialis, which may provide fundamental information to better understand the complex role played by this mechanism in regulating genome activity in bivalves.

The genome of Nautilus pompilius illuminates eye evolution and biomineralization

Nautilus is the sole surviving externally shelled cephalopod from the Palaeozoic. It is unique within cephalopod genealogy and critical to understanding the evolutionary novelties of cephalopods. Here, we present a complete Nautilus pompilius genome as a fundamental genomic reference on cephalopod innovations, such as the pinhole eye and biomineralization. Nautilus shows a compact, minimalist genome with few encoding genes and slow evolutionary rates in both non-coding and coding regions among known cephalopods. Importantly, multiple genomic innovations including gene losses, independent contraction and expansion of specific gene families and their associated regulatory networks likely moulded the evolution of the nautilus pinhole eye. The conserved molluscan biomineralization toolkit and lineage-specific repetitive low-complexity domains are essential to the construction of the nautilus shell. The nautilus genome constitutes a valuable resource for reconstructing the evolutionary scenarios and genomic innovations that shape the extant cephalopods.

Whole-genome microsynteny-based phylogeny of angiosperms

Plant genomes vary greatly in size, organization, and architecture. Such structural differences may be highly relevant for inference of genome evolution dynamics and phylogeny. Indeed, microsynteny-the conservation of local gene content and order-is recognized as a valuable source of phylogenetic information, but its use for the inference of large phylogenies has been limited. Here, by combining synteny network analysis, matrix representation, and maximum likelihood phylogenetic inference, we provide a way to reconstruct phylogenies based on microsynteny information. Both simulations and use of empirical data sets show our method to be accurate, consistent, and widely applicable. As an example, we focus on the analysis of a large-scale whole-genome data set for angiosperms, including more than 120 available high-quality genomes, representing more than 50 different plant families and 30 orders. Our 'microsynteny-based' tree is largely congruent with phylogenies proposed based on more traditional sequence alignment-based methods and current phylogenetic classifications but differs for some long-contested and controversial relationships. For instance, our synteny-based tree finds Vitales as early diverging eudicots, Saxifragales within superasterids, and magnoliids as sister to monocots. We discuss how synteny-based phylogenetic inference can complement traditional methods and could provide additional insights into some long-standing controversial phylogenetic relationships.

Genomic insights into the adaptation and evolution of the nautilus, an ancient but evolving "living fossil"

The nautilus, commonly known as a "living fossil," is endangered and may be at risk of extinction. The lack of genomic information hinders a thorough understanding of its biology and evolution, which can shed light on the conservation of this endangered species. Here, we report the first high-quality chromosome-level genome assembly of Nautilus pompilius. The assembled genome size comprised 785.15 Mb. Comparative genomic analyses indicated that transposable elements (TEs) and large-scale genome reorganizations may have driven lineage-specific evolution in the cephalopods. Remarkably, evolving conserved genes and recent TE insertion activities were identified in N. pompilius, and we speculate that these findings reflect the strong adaptability and long-term survival of the nautilus. We also identified gene families that are potentially responsible for specific adaptation and evolution events. Our study provides unprecedented insights into the specialized biology and evolution of N. pompilius, and the results serve as an important resource for future conservation genomics of the nautilus and closely related species.

The genome of the venomous snail Lautoconus ventricosus sheds light on the origin of conotoxin diversity

BACKGROUND

Venoms are deadly weapons to subdue prey or deter predators that have evolved independently in many animal lineages. The genomes of venomous animals are essential to understand the evolutionary mechanisms involved in the origin and diversification of venoms.

RESULTS

Here, we report the chromosome-level genome of the venomous Mediterranean cone snail, Lautoconus ventricosus (Caenogastropoda: Conidae). The total size of the assembly is 3.59 Gb; it has high contiguity (N50 = 93.53 Mb) and 86.6 Mb of the genome assembled into the 35 largest scaffolds or pseudochromosomes. On the basis of venom gland transcriptomes, we annotated 262 complete genes encoding conotoxin precursors, hormones, and other venom-related proteins. These genes were scattered in the different pseudochromosomes and located within repetitive regions. The genes encoding conotoxin precursors were normally structured into 3 exons, which did not necessarily coincide with the 3 structural domains of the corresponding proteins. Additionally, we found evidence in the L. ventricosus genome for a past whole-genome duplication event by means of conserved gene synteny with the Pomacea canaliculata genome, the only one available at the chromosome level within Caenogastropoda. The whole-genome duplication event was further confirmed by the presence of a duplicated hox gene cluster. Key genes for gastropod biology including those encoding proteins related to development, shell formation, and sex were located in the genome.

CONCLUSIONS

The new high-quality L. ventricosus genome should become a reference for assembling and analyzing new gastropod genomes and will contribute to future evolutionary genomic studies among venomous animals.

Mobilizing molluscan models and genomes in biology

Molluscs are among the most ancient, diverse, and important of all animal taxa. Even so, no individual mollusc species has emerged as a broadly applied model system in biology. We here make the case that both perceptual and methodological barriers have played a role in the relative neglect of molluscs as research organisms. We then summarize the current application and potential of molluscs and their genomes to address important questions in animal biology, and the state of the field when it comes to the availability of resources such as genome assemblies, cell lines, and other key elements necessary to mobilising the development of molluscan model systems. We conclude by contending that a cohesive research community that works together to elevate multiple molluscan systems to 'model' status will create new opportunities in addressing basic and applied biological problems, including general features of animal evolution. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Structural and Functional Characterization of the FGF Signaling Pathway in Regeneration of the Polychaete Worm Alitta virens (Annelida, Errantia)

Epimorphic regeneration of lost body segments is a widespread phenomenon across annelids. However, the molecular inducers of the cell sources for this reparative morphogenesis have not been identified. In this study, we focused on the role of fibroblast growth factor (FGF) signaling in the posterior regeneration of Alitta virens. For the first time, we showed an early activation of FGF ligands and receptor expression in an annelid regenerating after amputation. The expression patterns indicate that the entire regenerative bud is competent to FGFs, whose activity precedes the initiation of cell proliferation. The critical requirement of FGF signaling, especially at early stages, is also supported by inhibitor treatments followed by proliferation assay, demonstrating that induction of blastemal cells depends on FGFs. Our results show that FGF signaling pathway is a key player in regenerative response, while the FGF-positive wound epithelium, ventral nerve cord and some mesodermal cells around the gut could be the inducing tissues. This mechanism resembles reparative regeneration of vertebrate appendages suggesting such a response to the injury may be ancestral for all bilaterians.

Genomic Insights into the Origin and Evolution of Molluscan Red-Bloodedness in the Blood Clam Tegillarca granosa

Blood clams differ from their molluscan kins by exhibiting a unique red-blood (RB) phenotype; however, the genetic basis and biochemical machinery subserving this evolutionary innovation remain unclear. As a fundamental step toward resolving this mystery, we presented the first chromosome-level genome and comprehensive transcriptomes of the blood clam Tegillarca granosa for an integrated genomic, evolutionary, and functional analyses of clam RB phenotype. We identified blood clam-specific and expanded gene families, as well as gene pathways that are of RB relevant. Clam-specific RB-related hemoglobins (Hbs) showed close phylogenetic relationships with myoglobins (Mbs) of blood clam and other molluscs without the RB phenotype, indicating that clam-specific Hbs were likely evolutionarily derived from the Mb lineage. Strikingly, similar to vertebrate Hbs, blood clam Hbs were present in a form of gene cluster. Despite the convergent evolution of Hb clusters in blood clam and vertebrates, their Hb clusters may have originated from a single ancestral Mb-like gene as evidenced by gene phylogeny and synteny analysis. A full suite of enzyme-encoding genes for heme synthesis was identified in blood clam, with prominent expression in hemolymph and resembling those in vertebrates, suggesting a convergence of both RB-related Hb and heme functions in vertebrates and blood clam. RNA interference experiments confirmed the functional roles of Hbs and key enzyme of heme synthesis in the maintenance of clam RB phenotype. The high-quality genome assembly and comprehensive transcriptomes presented herein serve new genomic resources for the super-diverse phylum Mollusca, and provide deep insights into the origin and evolution of invertebrate RB.

How to Make a Rodent Giant: Genomic Basis and Tradeoffs of Gigantism in the Capybara, the World's Largest Rodent

Gigantism results when one lineage within a clade evolves extremely large body size relative to its small-bodied ancestors, a common phenomenon in animals. Theory predicts that the evolution of giants should be constrained by two tradeoffs. First, because body size is negatively correlated with population size, purifying selection is expected to be less efficient in species of large body size, leading to increased mutational load. Second, gigantism is achieved through generating a higher number of cells along with higher rates of cell proliferation, thus increasing the likelihood of cancer. To explore the genetic basis of gigantism in rodents and uncover genomic signatures of gigantism-related tradeoffs, we assembled a draft genome of the capybara (Hydrochoerus hydrochaeris), the world's largest living rodent. We found that the genome-wide ratio of nonsynonymous to synonymous mutations (ω) is elevated in the capybara relative to other rodents, likely caused by a generation-time effect and consistent with a nearly neutral model of molecular evolution. A genome-wide scan for adaptive protein evolution in the capybara highlighted several genes controlling postnatal bone growth regulation and musculoskeletal development, which are relevant to anatomical and developmental modifications for an increase in overall body size. Capybara-specific gene-family expansions included a putative novel anticancer adaptation that involves T-cell-mediated tumor suppression, offering a potential resolution to the increased cancer risk in this lineage. Our comparative genomic results uncovered the signature of an intragenomic conflict where the evolution of gigantism in the capybara involved selection on genes and pathways that are directly linked to cancer.

The Crown Pearl: a draft genome assembly of the European freshwater pearl mussel Margaritifera margaritifera (Linnaeus, 1758)

Since historical times, the inherent human fascination with pearls turned the freshwater pearl mussel Margaritifera margaritifera (Linnaeus, 1758) into a highly valuable cultural and economic resource. Although pearl harvesting in M. margaritifera is nowadays residual, other human threats have aggravated the species conservation status, especially in Europe. This mussel presents a myriad of rare biological features, e.g. high longevity coupled with low senescence and Doubly Uniparental Inheritance of mitochondrial DNA, for which the underlying molecular mechanisms are poorly known. Here, the first draft genome assembly of M. margaritifera was produced using a combination of Illumina Paired-end and Mate-pair approaches. The genome assembly was 2.4 Gb long, possessing 105,185 scaffolds and a scaffold N50 length of 288,726 bp. The ab initio gene prediction allowed the identification of 35,119 protein-coding genes. This genome represents an essential resource for studying this species’ unique biological and evolutionary features and ultimately will help to develop new tools to promote its conservation.

Molecular Characterization of Complement Component 3 (C3) in the Pearl Oyster Pinctada fucata Improves Our Understanding of the Primitive Complement System in Bivalve

As the central component in the complement system, complement component 3 (C3) plays essential roles in both the innate and adaptive immune responses. Here, a C3 gene (designated as pf-C3) was obtained from the pearl oyster Pinctada fucata by RT-PCR and rapid amplification of cDNA ends (RACE). The pf-C3 cDNA consists of 5,634 bp with an open reading frame (ORF) of 5,193 bp encoding a protein of 1,730 amino acids with a 19 residue signal peptide. The deduced pf-C3 protein possessed the characteristic structural features present in its homologs and contained the A2M_N_2, ANATO, A2M, A2M_comp, A2M_recep, and C345C domains, as well as the C3 convertase cleavage site, thioester motif, and conserved Cys, His, and Glu residues. Phylogenetic analysis revealed that pf-C3 is closely related to the C3s from other mollusks. Pf-C3 mRNA was expressed in all examined tissues including gill, digestive gland, adductor muscle, mantle and foot, while the highest expression was found in the digestive gland. Following the challenge with Vibrio alginolyticus, pf-C3 expression was significantly induced in hemocytes. Luciferase reporter assays indicated that pf-C3a could activate the NF-κB signal pathway in HEK293T cells. Further knockdown of pf-C3 by specific siRNA could significantly reduce the phagocytosis of V. alginolyticus by hemocytes in vitro. These results would help increase understanding of the function of C3 in the invertebrate immune system and therefore provide new insights into the roles of the primitive complement system in invertebrates.

MolluscDB: a genome and transcriptome database for molluscs

As sequencing becomes more accessible and affordable, the analysis of genomic and transcriptomic data has become a cornerstone of many research initiatives. Communities with a focus on particular taxa or ecosystems need solutions capable of aggregating genomic resources and serving them in a standardized and analysis-friendly manner. Taxon-focussed resources can be more flexible in addressing the needs of a research community than can universal or general databases. Here, we present MolluscDB, a genome and transcriptome database for molluscs. MolluscDB offers a rich ecosystem of tools, including an Ensembl browser, a BLAST server for homology searches and an HTTP server from which any dataset present in the database can be downloaded. To demonstrate the utility of the database and verify the quality of its data, we imported data from assembled genomes and transcriptomes of 22 species, estimated the phylogeny of Mollusca using single-copy orthologues, explored patterns of gene family size change and interrogated the data for biomineralization-associated enzymes and shell matrix proteins. MolluscDB provides an easy-to-use and openly accessible data resource for the research community. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Sticky problems: extraction of nucleic acids from molluscs

Traditional molecular methods and omics-techniques across molluscan taxonomy increasingly inform biology of Mollusca. Recovery of DNA and RNA for such studies is challenged by common biological properties of the highly diverse molluscs. Molluscan biomineralization, adhesive structures and mucus involve polyphenolic proteins and mucopolysaccharides that hinder DNA extraction or copurify to inhibit enzyme-catalysed molecular procedures. DNA extraction methods that employ the detergent hexadecyltrimethylammoniumbromide (CTAB) to remove these contaminants importantly facilitate molecular-level study of molluscs. Molluscan pigments may stain DNA samples and interfere with spectrophotometry, necessitating gel electrophoresis or fluorometry for accurate quantification. RNA can reliably be extracted but the 'hidden break' in 28S rRNA of molluscs (like most protostomes) causes 18S and 28S rRNA fragments to co-migrate electrophoretically. This challenges the standard quality control based on the ratio of 18S and 28S rRNA, developed for deuterostome animals. High-AT content in molluscan rRNA prevents the effective purification of polyadenylated mRNA. Awareness of these matters aids the continuous expansion of molecular malacology, enabling work also with museum specimens and next-generation sequencing, with the latter imposing unprecedented demands on DNA quality. Alternative methods to extract nucleic acids from molluscs are available from literature and, importantly, from communications with others who study the molecular biology of molluscs. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.

Transcriptional profiling of identified neurons in leech

Background

While leeches in the genus Hirudo have long been models for neurobiology, the molecular underpinnings of nervous system structure and function in this group remain largely unknown. To begin to bridge this gap, we performed RNASeq on pools of identified neurons of the central nervous system (CNS): sensory T (touch), P (pressure) and N (nociception) neurons; neurosecretory Retzius cells; and ganglia from which these four cell types had been removed.

Results

Bioinformatic analyses identified 3565 putative genes whose expression differed significantly among the samples. These genes clustered into 9 groups which could be associated with one or more of the identified cell types. We verified predicted expression patterns through in situ hybridization on whole CNS ganglia, and found that orthologous genes were for the most part similarly expressed in a divergent leech genus, suggesting evolutionarily conserved roles for these genes. Transcriptional profiling allowed us to identify candidate phenotype-defining genes from expanded gene families. Thus, we identified one of eight hyperpolarization-activated cyclic-nucleotide gated (HCN) channels as a candidate for mediating the prominent sag current in P neurons, and found that one of five inositol triphosphate receptors (IP3Rs), representing a sub-family of IP3Rs absent from vertebrate genomes, is expressed with high specificity in T cells. We also identified one of two piezo genes, two of ~ 65 deg/enac genes, and one of at least 16 transient receptor potential (trp) genes as prime candidates for involvement in sensory transduction in the three distinct classes of leech mechanosensory neurons.

Conclusions

Our study defines distinct transcriptional profiles for four different neuronal types within the leech CNS, in addition to providing a second ganglionic transcriptome for the species. From these data we identified five gene families that may facilitate the sensory capabilities of these neurons, thus laying the basis for future work leveraging the strengths of the leech system to investigate the molecular processes underlying and linking mechanosensation, cell type specification, and behavior.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07526-0.

Restriction of an intron size en route to endothermy

Ca2+-insensitive and -sensitive E1 subunits of the 2-oxoglutarate dehydrogenase complex (OGDHC) regulate tissue-specific NADH and ATP supply by mutually exclusive OGDH exons 4a and 4b. Here we show that their splicing is enforced by distant lariat branch points (dBPs) located near the 5' splice site of the intervening intron. dBPs restrict the intron length and prevent transposon insertions, which can introduce or eliminate dBP competitors. The size restriction was imposed by a single dominant dBP in anamniotes that expanded into a conserved constellation of four dBP adenines in amniotes. The amniote clusters exhibit taxon-specific usage of individual dBPs, reflecting accessibility of their extended motifs within a stable RNA hairpin rather than U2 snRNA:dBP base-pairing. The dBP expansion took place in early terrestrial species and was followed by a uridine enrichment of large downstream polypyrimidine tracts in mammals. The dBP-protected megatracts permit reciprocal regulation of exon 4a and 4b by uridine-binding proteins, including TIA-1/TIAR and PUF60, which promote U1 and U2 snRNP recruitment to the 5' splice site and BP, respectively, but do not significantly alter the relative dBP usage. We further show that codons for residues critically contributing to protein binding sites for Ca2+ and other divalent metals confer the exon inclusion order that mirrors the Irving-Williams affinity series, linking the evolution of auxiliary splicing motifs in exons to metallome constraints. Finally, we hypothesize that the dBP-driven selection for Ca2+-dependent ATP provision by E1 facilitated evolution of endothermy by optimizing the aerobic scope in target tissues.

Genomics and Transcriptomics of the green mussel explain the durability of its byssus

Mussels, which occupy important positions in marine ecosystems, attach tightly to underwater substrates using a proteinaceous holdfast known as the byssus, which is tough, durable, and resistant to enzymatic degradation. Although various byssal proteins have been identified, the mechanisms by which it achieves such durability are unknown. Here we report comprehensive identification of genes involved in byssus formation through whole-genome and foot-specific transcriptomic analyses of the green mussel, Perna viridis. Interestingly, proteins encoded by highly expressed genes include proteinase inhibitors and defense proteins, including lysozyme and lectins, in addition to structural proteins and protein modification enzymes that probably catalyze polymerization and insolubilization. This assemblage of structural and protective molecules constitutes a multi-pronged strategy to render the byssus highly resistant to environmental insults.

Functions of Forkhead Box O on Glucose Metabolism in Abalone Haliotis discus hannai and Its Responses to High Levels of Dietary Lipid

The forkhead box O (FoxO) subfamily is a member of the forkhead transcription factor family. It has regulation functions in glucose metabolism in mammals and fish. In the present study, a gene of the foxo homolog in abalone Haliotis discus hannai was cloned. A conservative forkhead (FH) domain and a transactivation (FoxO-TAD) domain were identified. Abalone foxo-specific siRNA (small interfering RNA) was injected to investigate the functions of foxo on glucose metabolism. Knockdown of foxo inhibited expression of phosphoenolpyruvate carboxykinase (pepck) and significantly increased expressions of hexokinase (hk) and pyruvate kinase (pk), but it failed to inhibit the relative mRNA level of glucose-6-phosphatase (g6pase). Then, a 100-day feeding trial was conducted to investigate the response of foxo and glucose metabolism in abalone fed with 1.57% (LFD, low-fat diet), 3.82% (MFD, middle-fat diet) and 6.72% (HFD, high-fat diet) of dietary lipid, respectively. The insulin-signaling pathway (AKT) was depressed and FoxO was activated by the HFD, but it did not inhibit glycolysis (hk) or improved gluconeogenesis significantly (pepck and g6pase). At the same time, impaired hepatopancreas glycogen storage raised hemolymph glucose levels. In conclusion, abalone foxo can be regulated by dietary lipid and can regulate gluconeogenesis or glycolysis in response to changes of dietary lipid levels, in which glycogen metabolism plays an important role.

Non-collinear Hox gene expression in bivalves and the evolution of morphological novelties in mollusks

Hox genes are key developmental regulators that are involved in establishing morphological features during animal ontogeny. They are commonly expressed along the anterior-posterior axis in a staggered, or collinear, fashion. In mollusks, the repertoire of body plans is widely diverse and current data suggest their involvement during development of landmark morphological traits in Conchifera, one of the two major lineages that comprises those taxa that originated from a uni-shelled ancestor (Monoplacophora, Gastropoda, Cephalopoda, Scaphopoda, Bivalvia). For most clades, and bivalves in particular, data on Hox gene expression throughout ontogeny are scarce. We thus investigated Hox expression during development of the quagga mussel, Dreissena rostriformis, to elucidate to which degree they might contribute to specific phenotypic traits as in other conchiferans. The Hox/ParaHox complement of Mollusca typically comprises 14 genes, 13 of which are present in bivalve genomes including Dreissena. We describe here expression of 9 Hox genes and the ParaHox gene Xlox during Dreissena development. Hox expression in Dreissena is first detected in the gastrula stage with widely overlapping expression domains of most genes. In the trochophore stage, Hox gene expression shifts towards more compact, largely mesodermal domains. Only few of these domains can be assigned to specific developing morphological structures such as Hox1 in the shell field and Xlox in the hindgut. We did not find traces of spatial or temporal staggered expression of Hox genes in Dreissena. Our data support the notion that Hox gene expression has been coopted independently, and to varying degrees, into lineage-specific structures in the respective conchiferan clades. The non-collinear mode of Hox expression in Dreissena might be a result of the low degree of body plan regionalization along the bivalve anterior-posterior axis as exemplified by the lack of key morphological traits such as a distinct head, cephalic tentacles, radula apparatus, and a simplified central nervous system.

Simple sequence repeats drive genome plasticity and promote adaptive evolution in penaeid shrimp

Simple sequence repeats (SSRs) are rare (approximately 1%) in most genomes and are generally considered to have no function. However, penaeid shrimp genomes have a high proportion of SSRs (>23%), raising the question of whether these SSRs play important functional and evolutionary roles in these SSR-rich species. Here, we show that SSRs drive genome plasticity and adaptive evolution in two penaeid shrimp species, Fenneropenaeus chinensis and Litopenaeus vannamei. Assembly and comparison of genomes of these two shrimp species at the chromosome-level revealed that transposable elements serve as carriers for SSR expansion, which is still occurring. The remarkable genome plasticity identified herein might have been shaped by significant SSR expansions. SSRs were also found to regulate gene expression by multi-omics analyses, and be responsible for driving adaptive evolution, such as the variable osmoregulatory capacities of these shrimp under low-salinity stress. These data provide strong evidence that SSRs are an important driver of the adaptive evolution in penaeid shrimp.

Evolutionary dynamics of transposable elements in bdelloid rotifers

Transposable elements (TEs) are selfish genomic parasites whose ability to spread autonomously is facilitated by sexual reproduction in their hosts. If hosts become obligately asexual, TE frequencies and dynamics are predicted to change dramatically, but the long-term outcome is unclear. Here, we test current theory using whole-genome sequence data from eight species of bdelloid rotifers, a class of invertebrates in which males are thus far unknown. Contrary to expectations, we find a variety of active TEs in bdelloid genomes, at an overall frequency within the range seen in sexual species. We find no evidence that TEs are spread by cryptic recombination or restrained by unusual DNA repair mechanisms. Instead, we find that that TE content evolves relatively slowly in bdelloids and that gene families involved in RNAi-mediated TE suppression have undergone significant expansion, which might mitigate the deleterious effects of active TEs and compensate for the consequences of long-term asexuality.

Analysis of Fox genes in Schmidtea mediterranea reveals new families and a conserved role of Smed-foxO in controlling cell death

The forkhead box (Fox) genes encode transcription factors that control several key aspects of development. Present in the ancestor of all eukaryotes, Fox genes underwent several duplications followed by loss and diversification events that gave rise to the current 25 families. However, few Fox members have been identified from the Lophotrochozoa clade, and specifically from planarians, which are a unique model for understanding development, due to the striking plasticity of the adult. The aim of this study was to identify and perform evolutionary and functional studies of the Fox genes of lophotrochozoan species and, specifically, of the planarian Schmidtea mediterranea. Generating a pipeline for identifying Forkhead domains and using phylogenetics allowed us the phylogenetic reconstruction of Fox genes. We corrected the annotation for misannotated genes and uncovered a new family, the QD, present in all metazoans. According to the new phylogeny, the 27 Fox genes found in Schmidtea mediterranea were classified into 12 families. In Platyhelminthes, family losses were accompanied by extensive gene diversification and the appearance of specific families, the A(P) and N(P). Among the newly identified planarian Fox genes, we found a single copy of foxO, which shows an evolutionary conserved role in controlling cell death.

MolluscDB: an integrated functional and evolutionary genomics database for the hyper-diverse animal phylum Mollusca

Mollusca represents the second largest animal phylum but remains poorly explored from a genomic perspective. While the recent increase in genomic resources holds great promise for a deep understanding of molluscan biology and evolution, access and utilization of these resources still pose a challenge. Here, we present the first comprehensive molluscan genomics database, MolluscDB (http://mgbase.qnlm.ac), which compiles and integrates current molluscan genomic/transcriptomic resources and provides convenient tools for multi-level integrative and comparative genomic analyses. MolluscDB enables a systematic view of genomic information from various aspects, such as genome assembly statistics, genome phylogenies, fossil records, gene information, expression profiles, gene families, transcription factors, transposable elements and mitogenome organization information. Moreover, MolluscDB offers valuable customized datasets or resources, such as gene coexpression networks across various developmental stages and adult tissues/organs, core gene repertoires inferred for major molluscan lineages, and macrosynteny analysis for chromosomal evolution. MolluscDB presents an integrative and comprehensive genomics platform that will allow the molluscan community to cope with ever-growing genomic resources and will expedite new scientific discoveries for understanding molluscan biology and evolution.

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.

A Transcriptomic Approach to the Recruitment of Venom Proteins in a Marine Annelid

The growing number of known venomous marine invertebrates indicates that chemical warfare plays an important role in adapting to diversified ecological niches, even though it remains unclear how toxins fit into the evolutionary history of these animals. Our case study, the Polychaeta Eulalia sp., is an intertidal predator that secretes toxins. Whole-transcriptome sequencing revealed proteinaceous toxins secreted by cells in the proboscis and delivered by mucus. Toxins and accompanying enzymes promote permeabilization, coagulation impairment and the blocking of the neuromuscular activity of prey upon which the worm feeds by sucking pieces of live flesh. The main neurotoxins ("phyllotoxins") were found to be cysteine-rich proteins, a class of substances ubiquitous among venomous animals. Some toxins were phylogenetically related to Polychaeta, Mollusca or more ancient groups, such as Cnidaria. Some toxins may have evolved from non-toxin homologs that were recruited without the reduction in molecular mass and increased specificity of other invertebrate toxins. By analyzing the phylogeny of toxin mixtures, we show that Polychaeta is uniquely positioned in the evolution of animal venoms. Indeed, the phylogenetic models of mixed or individual toxins do not follow the expected eumetazoan tree-of-life and highlight that the recruitment of gene products for a role in venom systems is complex.

Resolving Transcriptional States and Predicting Lineages in the Annelid Capitella teleta Using Single-Cell RNAseq

Evolution and diversification of cell types has contributed to animal evolution. However, gene regulatory mechanisms underlying cell fate acquisition during development remains largely uncharacterized in spiralians. Here we use a whole-organism, single-cell transcriptomic approach to map larval cell types in the annelid Capitella teleta at 24- and 48-h post gastrulation (stages 4 and 5). We identified eight unique cell clusters (undifferentiated precursors, ectoderm, muscle, ciliary-band, gut, neurons, neurosecretory cells, and protonephridia), thus helping to identify uncharacterized molecular signatures such as previously unknown neurosecretory cell markers in C. teleta. Analysis of coregulatory programs in individual clusters revealed gene interactions that can be used for comparisons of cell types across taxa. We examined the neural and neurosecretory clusters more deeply and characterized a differentiation trajectory starting from dividing precursors to neurons using Monocle3 and velocyto. Pseudotime analysis along this trajectory identified temporally-distinct cell states undergoing progressive gene expression changes over time. Our data revealed two potentially distinct neural differentiation trajectories including an early trajectory for brain neurosecretory cells. This work provides a valuable resource for future functional investigations to better understanding neurogenesis and the transitions from neural precursors to neurons in an annelid.

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.

Sequencing of Transcriptome and Small RNA Revealed the Xenobiotic Metabolism-Related Genes and Potential Regulatory miRNA in Asian Tramp Snail

The Asian tramp snail, Bradybaena similaris (Ferusssac), is an invasive land snail species and has been a rising agricultural pest in south of China. As a pest, it also plays a role in transmission of Angiostrongylus cantonensis. However, present studies on this species are rare and the molecular information is limited. For this purpose, we sequenced the transcriptome and small RNA of B. similaris collected from citrus orchards. In total, 89,747 unigenes with an N50 size of 1287 bp and an average length of 817 bp were generated from ∼8.9 Gb transcriptome and 31 Mb clean reads were generated from ∼36 Mb small RNA library. To demonstrate the usefulness of these two datasets, we analyzed a series of genes associated with xenobiotic metabolism and core RNAi machinery. Analysis of the transcripts resulted in annotation of 126 putative genes encoding cytochrome P450 monooxygenases (CYP, 45), carboxyl/cholinesterases (CCE, 13), glutathione-S-transferases (GST, 24), and ATP-binding cassette transporters (ABC, 44). Analysis of the small RNA detected 42 miRNAs. In addition, four genes involved in small RNA pathways (miRNA, piRNA, and siRNA) were identified, and a total of 430 genes that can be targeted by miRNAs were predicted. Moreover, we found that a few miRNAs could target certain genes involved in xenobiotic metabolism. Therefore, we believe that these two datasets and the characterization of the identified/predicted genes will facilitate the molecular study of this species as well as other land snails with agricultural importance.

The Iron-Responsive Genome of the Chiton Acanthopleura granulata

Molluscs biomineralize structures that vary in composition, form, and function, prompting questions about the genetic mechanisms responsible for their production and the evolution of these mechanisms. Chitons (Mollusca, Polyplacophora) are a promising system for studies of biomineralization because they build a range of calcified structures including shell plates and spine- or scale-like sclerites. Chitons also harden the calcified teeth of their rasp-like radula with a coat of iron (as magnetite). Here we present the genome of the West Indian fuzzy chiton Acanthopleura granulata, the first from any aculiferan mollusc. The A. granulata genome contains homologs of many genes associated with biomineralization in conchiferan molluscs. We expected chitons to lack genes previously identified from pathways conchiferans use to make biominerals like calcite and nacre because chitons do not use these materials in their shells. Surprisingly, the A. granulata genome has homologs of many of these genes, suggesting that the ancestral mollusc may have had a more diverse biomineralization toolkit than expected. The A. granulata genome has features that may be specialized for iron biomineralization, including a higher proportion of genes regulated directly by iron than other molluscs. A. granulata also produces two isoforms of soma-like ferritin: one is regulated by iron and similar in sequence to the soma-like ferritins of other molluscs, and the other is constitutively translated and is not found in other molluscs. The A. granulata genome is a resource for future studies of molluscan evolution and biomineralization.

Evolution and Potential Function in Molluscs of Neuropeptide and Receptor Homologues of the Insect Allatostatins

The allatostatins (ASTs), AST-A, AST-B and AST-C, have mainly been investigated in insects. They are a large group of small pleotropic alloregulatory neuropeptides that are unrelated in sequence and activate receptors of the rhodopsin G-protein coupled receptor family (GPCRs). The characteristics and functions of the homologue systems in the molluscs (Buccalin, MIP and AST-C-like), the second most diverse group of protostomes after the arthropods, and of high interest for evolutionary studies due to their less rearranged genomes remains to be explored. In the present study their evolution is deciphered in molluscs and putative functions assigned in bivalves through meta-analysis of transcriptomes and experiments. Homologues of the three arthropod AST-type peptide precursors were identified in molluscs and produce a larger number of mature peptides than in insects. The number of putative receptors were also distinct across mollusc species due to lineage and species-specific duplications. Our evolutionary analysis of the receptors identified for the first time in a mollusc, the cephalopod, GALR-like genes, which challenges the accepted paradigm that AST-AR/buccalin-Rs are the orthologues of vertebrate GALRs in protostomes. Tissue transcriptomes revealed the peptides, and their putative receptors have a widespread distribution in bivalves and in the bivalve Mytilus galloprovincialis, elements of the three peptide-receptor systems are highly abundant in the mantle an innate immune barrier tissue. Exposure of M. galloprovincialis to lipopolysaccharide or a marine pathogenic bacterium, Vibrio harveyi, provoked significant modifications in the expression of genes of the peptide precursor and receptors of the AST-C-like system in the mantle suggesting involvement in the immune response. Overall, our study reveals that homologues of the arthropod AST-systems in molluscs are potentially more complex due to the greater number of putative mature peptides and receptor genes. In bivalves they have a broad and varying tissue distribution and abundance, and the elements of the AST-C-like family may have a putative function in the immune response.

Fish lysozyme gene family evolution and divergent function in early development

Lysozymes are an ancient group of antimicrobial enzymes of the innate immune system. Here we provide a comparative analysis of the evolution and function of lysozymes during early development in fish, the most speciose vertebrate group. In fishes, lineage and species-specific evolution of both C-type (chicken or conventional) and G-type (goose type) genes occurred. Phylogenetic analysis revealed that the teleost lysozyme G-type members group with the tetrapod homologues but the teleost C-type form three different clusters with the tetrapods. Most of the teleost C-type cluster with tetrapod Lyz but there are some that group with the mammalian Lyzl1/2 and LALBA. This suggests that early in gnathostome evolution these genes already existed and that lyzl1/2 and lalba genes are present in fish and tetrapods. Gene synteny analysis to confirm sequence orthologies failed to identify conserved genome regions between teleosts and other vertebrates lysozyme gene regions suggesting that in the ancestral bony fish genome lyz, lyzl1/2, lalba and lyg precursor genes were transposed to different chromosome regions. The homologue of the mammalian lactalbumin (LALBA) gene was identified for the first time in teleosts and was expressed in skin and during egg and larval development. Lysozyme activity was detected in teleost eggs and varied between species and in the gilthead sea bream lyg and lalba transcript abundance differed in eggs and larvae from different brood stock suggesting differences exist in maternal innate immune protection.

The transcriptome of anterior regeneration in earthworm Eudrilus eugeniae

The oligochaete earthworm, Eudrilus eugeniae is capable of regenerating both anterior and posterior segments. The present study focuses on the transcriptome analysis of earthworm E. eugeniae to identify and functionally annotate the key genes supporting the anterior blastema formation and regulating the anterior regeneration of the worm. The Illumina sequencing generated a total of 91,593,182 raw reads which were assembled into 105,193 contigs using CLC genomics workbench. In total, 40,946 contigs were annotated against the NCBI nr and SwissProt database and among them, 15,702 contigs were assigned to 14,575 GO terms. Besides a total of 9389 contigs were mapped to 416 KEGG biological pathways. The RNA-Seq comparison study identified 10,868 differentially expressed genes (DEGs) and of them, 3986 genes were significantly upregulated in the anterior regenerated blastema tissue samples of the worm. The GO enrichment analysis showed angiogenesis and unfolded protein binding as the top enriched functions and the pathway enrichment analysis denoted TCA cycle as the most significantly enriched pathway associated with the upregulated gene dataset of the worm. The identified DEGs and their function and pathway information can be effectively utilized further to interpret the key cellular, genetic and molecular events associated with the regeneration of the worm.

Genomic and Transcriptomic Analyses of Bioluminescence Genes in the Enope Squid Watasenia scintillans

Watasenia scintillans, a sparkling enope squid, has bioluminescence organs to illuminate its body with its own luciferase activity. To clarify the molecular mechanism underlying its scintillation, we analysed high-throughput sequencing data acquired previously and obtained draft genome sequences accomplished with comparative genomic data among the cephalopods. The genome mapped by transcriptome data showed that (1) RNA editing contributed to transcriptome variation of lineage specific genes, such as W. scintillans luciferase, and (2) two types of luciferase enzymes were characterized with reasonable 3D models docked to a luciferin molecule. We report two different types of luciferase in one organism and possibly related to variety of colour types in the W. scintillans fluorescent organs.

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.

Extensive gene rearrangements in the mitogenomes of congeneric annelid species and insights on the evolutionary history of the genus Ophryotrocha

BACKGROUND

Annelids are one the most speciose and ecologically diverse groups of metazoans. Although a significant effort has been recently invested in sequencing genomes of a wide array of metazoans, many orders and families within the phylum Annelida are still represented by a single specimen of a single species. The genus of interstitial annelids Ophryotrocha (Dorvilleidae, Errantia, Annelida) is among these neglected groups, despite its extensive use as model organism in numerous studies on the evolution of life history, physiological and ecological traits. To compensate for the paucity of genomic information in this genus, we here obtained novel complete mitochondrial genomes of six Ophryotrocha species using next generation sequencing. In addition, we investigated the evolution of the reproductive mode in the Ophryotrocha genus using a phylogeny based on two mitochondrial markers (COXI and 16S rDNA) and one nuclear fragment (Histone H3).

RESULTS

Surprisingly, gene order was not conserved among the six Ophryotrocha species investigated, and varied greatly as compared to those found in other annelid species within the class Errantia. The mitogenome phylogeny for the six Ophryotrocha species displayed a separation of gonochoric and hermaphroditic species. However, this separation was not observed in the phylogeny based on the COX1, 16S rDNA, and H3 genes. Parsimony and Bayesian ancestral trait reconstruction indicated that gonochorism was the most parsimonious ancestral reproductive mode in Ophryotrocha spp.

CONCLUSIONS

Our results highlight the remarkably high level of gene order variation among congeneric species, even in annelids. This encourages the need for additional mitogenome sequencing of annelid taxa in order to properly understand its mtDNA evolution, high biodiversity and phylogenetic relationships.

Conservative route to genome compaction in a miniature annelid

The causes and consequences of genome reduction in animals are unclear because our understanding of this process mostly relies on lineages with often exceptionally high rates of evolution. Here, we decode the compact 73.8-megabase genome of Dimorphilus gyrociliatus, a meiobenthic segmented worm. The D. gyrociliatus genome retains traits classically associated with larger and slower-evolving genomes, such as an ordered, intact Hox cluster, a generally conserved developmental toolkit and traces of ancestral bilaterian linkage. Unlike some other animals with small genomes, the analysis of the D. gyrociliatus epigenome revealed canonical features of genome regulation, excluding the presence of operons and trans-splicing. Instead, the gene-dense D. gyrociliatus genome presents a divergent Myc pathway, a key physiological regulator of growth, proliferation and genome stability in animals. Altogether, our results uncover a conservative route to genome compaction in annelids, reminiscent of that observed in the vertebrate Takifugu rubripes.

This study reports the genome of the miniature segmented annelid Dimorphilus gyrociliatus and reveals no drastic changes in genome architecture and regulation, unlike other cases of genome miniaturization.

Genomic Landscape of Mutational Biases in the Pacific Oyster Crassostrea gigas

Mutation is a driving force of evolution that has been shaped by natural selection and is universally biased. Previous studies determined genome-wide mutational patterns for several species and investigated the heterogeneity of mutational patterns at fine-scale levels. However, little evidence of the heterogeneity of mutation rates over large genomic regions was shown. Hence, the mutational patterns of different large-scale genomic regions and their association with selective pressures still need to be explored. As the second most species-rich animal phylum, little is known about the mutational patterns in Mollusca, especially oysters. In this study, the mutational bias patterns are characterized by using whole-genome resequencing data in the Crassostrea gigas genome. I studied the genome-wide relative rates of the pair mutations and found that the predominant mutation is GC -> AT, irrespective of the genomic regions. This analysis reveals that mutational biases were associated with gene expression levels across the C. gigas genome. Genes with higher expression levels and breadth expression patterns, longer coding length, and more exon numbers had relatively higher GC -> AT rates. I also found that genes with larger dN/dS values had relatively higher GC -> AT rates. This work represents the first comprehensive research on the mutational biases in Mollusca species. Here, I comprehensively investigated the relationships between mutational biases with some intrinsic genetic factors and evolutionary indicators and proposed that selective pressures are important forces shaping the mutational biases across the C. gigas genome.

Genome biology of the paleotetraploid perennial biomass crop Miscanthus

Miscanthus is a perennial wild grass that is of global importance for paper production, roofing, horticultural plantings, and an emerging highly productive temperate biomass crop. We report a chromosome-scale assembly of the paleotetraploid M. sinensis genome, providing a resource for Miscanthus that links its chromosomes to the related diploid Sorghum and complex polyploid sugarcanes. The asymmetric distribution of transposons across the two homoeologous subgenomes proves Miscanthus paleo-allotetraploidy and identifies several balanced reciprocal homoeologous exchanges. Analysis of M. sinensis and M. sacchariflorus populations demonstrates extensive interspecific admixture and hybridization, and documents the origin of the highly productive triploid bioenergy crop M. × giganteus. Transcriptional profiling of leaves, stem, and rhizomes over growing seasons provides insight into rhizome development and nutrient recycling, processes critical for sustainable biomass accumulation in a perennial temperate grass. The Miscanthus genome expands the power of comparative genomics to understand traits of importance to Andropogoneae grasses.

The perennial grass Miscanthus is a promising biomass crop. Here, via genomics and transcriptomics, the authors reveal its allotetraploid origin, characterize gene expression associated with rhizome development and nutrient recycling, and describe the hybrid origin of the triploid M. x giganteus.