Host-Endosymbiont Genome Integration in a Deep-Sea Chemosymbiotic Clam

Comparative study of integrated bio-responses in deep-sea and nearshore mussels upon abiotic condition changes: Insight into distinct regulation and adaptation

Deep-sea mussels, one of the dominant species in most deep-sea ecosystems, have long been used as model organisms to investigate the adaptations and symbiotic relationships of deep-sea macrofauna under laboratory conditions due to their ability to survive under atmospheric pressure. However, the impact of additional abiotic conditions beyond pressure, such as temperature and light, on their physiological characteristics remains unknown. In this study, deep-sea mussels (Gigantidas platifrons) from cold seep of the South China Sea, along with nearshore mussels (Mytilus coruscus) from the East China Sea, were reared in unfavorable abiotic conditions for up to 8 days. Integrated biochemical indexes including antioxidant defense, immune ability and energy metabolism were investigated in the gill and digestive gland, while cytotoxicity was determined in hemocytes of both types of mussels. The results revealed mild bio-responses in two types of mussels in the laboratory, represented by the effective antioxidant defense with constant total antioxidant capability level and malondialdehyde content. There were also disparate adaptations in deep-sea and nearshore mussels. In deep-sea mussels, significantly increased immune response and energy reservation were observed in gills, together with the elevated cytotoxicity in hemocytes, implying the more severe biological adaptation was required, mainly due to the symbiotic bacteria loss under laboratory conditions. On the contrary, insignificant biological responses were exhibited in nearshore mussels except for the increased energy consumption, indicating the trade-off strategy to use more energy to deal with potential stress. Overall, this comparative study highlights the basal bio-responses of deep-sea and nearshore mussels out of their native environments, providing evidence that short-term culture of both mussels under easily achievable laboratory conditions would not dramatically alter their biological status. This finding will assist in broadening the application of deep-sea mussels as model organism in future research regardless of the specialized research equipment.

Endosymbioses Have Shaped the Evolution of Biological Diversity and Complexity Time and Time Again

Life on Earth comprises prokaryotes and a broad assemblage of endosymbioses. The pages of Molecular Biology and Evolution and Genome Biology and Evolution have provided an essential window into how these endosymbiotic interactions have evolved and shaped biological diversity. Here, we provide a current perspective on this knowledge by drawing on decades of revelatory research published in Molecular Biology and Evolution and Genome Biology and Evolution, and insights from the field at large. The accumulated work illustrates how endosymbioses provide hosts with novel phenotypes that allow them to transition between adaptive landscapes to access environmental resources. Such endosymbiotic relationships have shaped and reshaped life on Earth. The early serial establishment of mitochondria and chloroplasts through endosymbioses permitted massive upscaling of cellular energetics, multicellularity, and terrestrial planetary greening. These endosymbioses are also the foundation upon which all later ones are built, including everything from land-plant endosymbioses with fungi and bacteria to nutritional endosymbioses found in invertebrate animals. Common evolutionary mechanisms have shaped this broad range of interactions. Endosymbionts generally experience adaptive and stochastic genome streamlining, the extent of which depends on several key factors (e.g. mode of transmission). Hosts, in contrast, adapt complex mechanisms of resource exchange, cellular integration and regulation, and genetic support mechanisms to prop up degraded symbionts. However, there are significant differences between endosymbiotic interactions not only in how partners have evolved with each other but also in the scope of their influence on biological diversity. These differences are important considerations for predicting how endosymbioses will persist and adapt to a changing planet.

Pan-evolutionary and regulatory genome architecture delineated by an integrated macro- and microsynteny approach

The forthcoming massive genome data generated by the Earth BioGenome Project will open up a new era of comparative genomics, for which genome synteny analysis provides an important framework. Profiling genome synteny represents an essential step in elucidating genome architecture, regulatory blocks/elements and their evolutionary history. Here we describe PanSyn, ( https://github.com/yhw320/PanSyn ), the most comprehensive and up-to-date genome synteny pipeline, providing step-by-step instructions and application examples to demonstrate its usage. PanSyn inherits both basic and advanced functions from existing popular tools, offering a user-friendly, highly customized approach for genome macrosynteny analysis and integrated pan-evolutionary and regulatory analysis of genome architecture, which are not yet available in public synteny software or tools. The advantages of PanSyn include: (i) advanced microsynteny analysis by functional profiling of microsynteny genes and associated regulatory elements; (ii) comprehensive macrosynteny analysis, including the inference of karyotype evolution from ancestors to extant species; and (iii) functional integration of microsynteny and macrosynteny for pan-evolutionary profiling of genome architecture and regulatory blocks, as well as integration with external functional genomics datasets from three- or four-dimensional genome and ENCODE projects. PanSyn requires basic knowledge of the Linux environment and Perl programming language and the ability to access a computer cluster, especially for large-scale genomic comparisons. Our protocol can be easily implemented by a competent graduate student or postdoc and takes several days to weeks to execute for dozens to hundreds of genomes. PanSyn provides yet the most comprehensive and powerful tool for integrated evolutionary and functional genomics.

Scallop-bacteria symbiosis from the deep sea reveals strong genomic coupling in the absence of cellular integration

Previous studies have revealed tight metabolic complementarity between bivalves and their endosymbiotic chemosynthetic bacteria, but little is known about their interactions with ectosymbionts. Our analysis of the ectosymbiosis between a deep-sea scallop (Catillopecten margaritatus) and a gammaproteobacterium showed that bivalves could be highly interdependent with their ectosymbionts as well. Our microscopic observation revealed abundant sulfur-oxidizing bacteria (SOB) on the surfaces of the gill epithelial cells. Microbial 16S rRNA gene amplicon sequencing of the gill tissues showed the dominance of the SOB. An analysis of the SOB genome showed that it is substantially smaller than its free-living relatives and has lost cellular components required for free-living. Genomic and transcriptomic analyses showed that this ectosymbiont relies on rhodanese-like proteins and SOX multienzyme complex for energy generation, mainly on the Calvin-Benson-Bassham (CBB) cycle and peripherally on a phosphoenolpyruvate carboxylase for carbon assimilation. Besides, the symbiont encodes an incomplete tricarboxylic acid (TCA) cycle. Observation of the scallop's digestive gland and its nitrogen metabolism pathways indicates it does not fully rely on the ectosymbiont for nutrition. Analysis of the host's gene expression provided evidence that it could offer intermediates for the ectosymbiont to complete its TCA cycle and some amino acid synthesis pathways using exosomes, and its phagosomes, endosomes, and lysosomes might be involved in harvesting nutrients from the symbionts. Overall, our study prompts us to rethink the intimacy between the hosts and ectosymbionts in Bivalvia and the evolution of chemosymbiosis in general.

Hemoglobin wonders: a fascinating gas transporter dive into molluscs

Hemoglobin (Hb) has been identified in at least 14 molluscan taxa so far. Research spanning over 130 years on molluscan Hbs focuses on their genes, protein structures, functions, and evolution. Molluscan Hbs are categorized into single-, two-, and multiple-domain chains, including red blood cell, gill, and extracellular Hbs, based on the number of globin domains and their respective locations. These Hbs exhibit variation in assembly, ranging from monomeric and dimeric to higher-order multimeric forms. Typically, molluscan Hbs display moderately high oxygen affinity, weak cooperativity, and varying pH sensitivity. Hb's potential role in antimicrobial pathways could augment the immune defense of bivalves, which may be a complement to their lack of adaptive immunity. The role of Hb as a respiratory protein in bivalves likely originated from the substitution of hemocyanin. Molluscan Hbs demonstrate adaptive evolution in response to environmental changes via various strategies (e.g. increasing Hb types, multimerization, and amino acid residue substitutions at key sites), enhancing or altering functional properties for habitat adaptation. Concurrently, an increase in Hb assembly diversity, coupled with a downward trend in oxygen affinity, is observed during molluscan differentiation and evolution. Hb in Protobranchia, Heteroconchia, and Pteriomorphia bivalves originated from separate ancestors, with Protobranchia inheriting a relative ancient molluscan Hb gene. In bivalves, extracellular Hbs share a common origin, while gill Hbs likely emerged from convergent evolution. In summary, research on molluscan Hbs offers valuable insights into the origins, biological variations, and adaptive evolution of animal Hbs.

Reduced chemosymbiont genome in the methane seep thyasirid and the cooperated metabolisms in the holobiont under anaerobic sediment

Previous studies have deciphered the genomic basis of host-symbiont metabolic complementarity in vestimentiferans, bathymodioline mussels, vesicomyid clams and Alviniconcha snails, yet little is known about the chemosynthetic symbiosis in Thyasiridae-a family of Bivalvia regarded as an excellent model in chemosymbiosis research due to their wide distribution in both deep-sea and shallow-water habitats. We report the first circular thyasirid symbiont genome, named Candidatus Ruthturnera sp. Tsphm01, with a size of 1.53 Mb, 1521 coding genes and 100% completeness. Compared to its free-living relatives, Ca. Ruthturnera sp. Tsphm01 genome is reduced, lacking components for chemotaxis, citric acid cycle and de novo biosynthesis of small molecules (e.g. amino acids and cofactors), indicating it is likely an obligate intracellular symbiont. Nevertheless, the symbiont retains complete genomic components of sulphur oxidation and assimilation of inorganic carbon, and these systems were highly and actively expressed. Moreover, the symbiont appears well-adapted to anoxic environment, including capable of anaerobic respiration (i.e. reductions of DMSO and nitrate) and possession of a low oxygen-adapted type of cytochrome c oxidase. Analysis of the host transcriptome revealed its metabolic complementarity to the incomplete metabolic pathways of the symbiont and the acquisition of nutrients from the symbiont via phagocytosis and exosome. By providing the first complete genome of reduced size in a thyasirid symbiont, this study enhances our understanding of the diversity of symbiosis that has enabled bivalves to thrive in chemosynthetic habitats. The resources will be widely used in phylogenetic, geographic and evolutionary studies of chemosynthetic bacteria and bivalves.

The genome of the deep-sea anemone Actinernus sp. contains a mega-array of ANTP-class homeobox genes

Members of the phylum Cnidaria include sea anemones, corals and jellyfish, and have successfully colonized both marine and freshwater habitats throughout the world. The understanding of how cnidarians adapt to extreme environments such as the dark, high-pressure deep-sea habitat has been hindered by the lack of genomic information. Here, we report the first chromosome-level deep-sea cnidarian genome, of the anemone Actinernus sp., which was 1.39 Gbp in length and contained 44 970 gene models including 14 806 tRNA genes and 30 164 protein-coding genes. Analyses of homeobox genes revealed the longest chromosome hosts a mega-array of Hox cluster, HoxL, NK cluster and NKL homeobox genes; until now, such an array has only been hypothesized to have existed in ancient ancestral genomes. In addition to this striking arrangement of homeobox genes, analyses of microRNAs revealed cnidarian-specific complements that are distinctive for nested clades of these animals, presumably reflecting the progressive evolution of the gene regulatory networks in which they are embedded. Also, compared with other sea anemones, circadian rhythm genes were lost in Actinernus sp., which likely reflects adaptation to living in the dark. This high-quality genome of a deep-sea cnidarian thus reveals some of the likely molecular adaptations of this ecologically important group of metazoans to the extreme deep-sea environment. It also deepens our understanding of the evolution of genome content and organization of animals in general and cnidarians in particular, specifically from the viewpoint of key developmental control genes like the homeobox-encoding genes, where we find an array of genes that until now has only been hypothesized to have existed in the ancient ancestor that pre-dated both the cnidarians and bilaterians.

The First Genome Survey of the Snail Provanna glabra Inhabiting Deep-Sea Hydrothermal Vents

The snail P. glabra is an endemic species in deep-sea chemosynthetic ecosystems of the Northwest Pacific Ocean. To obtain more genetic information on this species and provide the basis for subsequent whole-genome map construction, a genome survey was performed on this snail from the hydrothermal vent of Okinawa Trough. The genomic size of P. glabra was estimated to be 1.44 Gb, with a heterozygosity of 1.91% and a repeated sequence content of 69.80%. Based on the sequencing data, a draft genome of 1.32 Gb was assembled. Transposal elements (TEs) accounted for 40.17% of the entire genome, with DNA transposons taking the highest proportion. It was found that most TEs were inserted in the genome recently. In the simple sequence repeats, the dinucleotide motif was the most enriched microsatellite type, accounting for 53% of microsatellites. A complete mitochondrial genome of P. glabra with a total length of 16,268 bp was assembled from the sequencing data. After comparison with the published mitochondrial genome of Provanna sp. from a methane seep, 331 potential single nucleotide polymorphism (SNP) sites were identified in protein-coding genes (PCGs). Except for the cox1 gene, nad2, nad4, nad5, and cob genes are expected to be candidate markers for population genetic and phylogenetic studies of P. glabra and other deep-sea snails. Compared with shallow-water species, three mitochondrial genes of deep-sea gastropods exhibited a higher evolutionary rate, indicating strong selection operating on mitochondria of deep-sea species. This study provides insights into the genome characteristics of P. glabra and supplies genomic resources for further studies on the adaptive evolution of the snail in extreme deep-sea chemosynthetic environments.

Bivalves as Emerging Model Systems to Study the Mechanisms and Evolution of Sex Determination: A Genomic Point of View

Bivalves are a diverse group of molluscs that have recently attained a central role in plenty of biological research fields, thanks to their peculiar life history traits. Here, we propose that bivalves should be considered as emerging model systems also in sex-determination (SD) studies, since they would allow to investigate: 1) the transition between environmental and genetic SD, with respect to different reproductive backgrounds and sexual systems (from species with strict gonochorism to species with various forms of hermaphroditism); 2) the genomic evolution of sex chromosomes (SCs), considering that no heteromorphic SCs are currently known and that homomorphic SCs have been identified only in a few species of scallops; 3) the putative role of mitochondria at some level of the SD signaling pathway, in a mechanism that may resemble the cytoplasmatic male sterility of plants; 4) the evolutionary history of SD-related gene (SRG) families with respect to other animal groups. In particular, we think that this last topic may lay the foundations for expanding our understanding of bivalve SD, as our current knowledge is quite fragmented and limited to a few species. As a matter of fact, tracing the phylogenetic history and diversity of SRG families (such as the Dmrt, Sox, and Fox genes) would allow not only to perform more targeted functional experiments and genomic analyses, but also to foster the possibility of establishing a solid comparative framework.

Complete genome analysis reveals environmental adaptability of sulfur-oxidizing bacterium Thioclava nitratireducens M1-LQ-LJL-11 and symbiotic relationship with deep-sea hydrothermal vent Chrysomallon squamiferum

One sulfur-oxidizing bacterium Thioclava sp. M1-LQ-LJL-11 was isolated from the gill of Chrysomallon squamiferum collected from 2700 m deep hydrothermal named Longqi on the southwest Indian Ocean ridge. In order to understand its survival mechanism in hydrothermal extreme environment and symbiotic relationship with its host, the complete genome of strain M1-LQ-LJL-11 was sequenced and analyzed. A total of 6117 Mb of valid data was obtained, including 4096 coding genes, 61 non coding genes, including 9 rRNAs (among them, there are 3 in 23S rRNA, 3 in 5S rRNA, and 3 in 16S rRNA.), 52 tRNAs and 35 genomic islands. Strain M1-LQ-LJL-11 contains one chromosome and two plasmids. In the genome annotation information of the strain, we found 28 genes including cys sox, sor, sqr, tst related to sulfur metabolism and 17 metal resistance genes. Interestingly, a pair of quorum sensing system which probably regulating biofilm formation located in chromosome was found. These genes are critical for self-adaptation against severe environment as well as host survival. This study provides a basis understanding for the adaptive strategies of deep-sea hydrothermal bacteria and symbiotic relationship with its host in extreme environments through gene level.

Early genome erosion and internal phage-symbiont-host interaction in the endosymbionts of a cold-seep tubeworm

Endosymbiosis with chemosynthetic Gammaproteobacteria is widely recognized as an adaptive mechanism of siboglinid tubeworms, yet evolution of these endosymbionts and their driving forces remain elusive. Here, we report a finished endosymbiont genome (HMS1) of the cold-seep tubeworm Sclerolinum annulatum. The HMS1 genome is small in size, with abundant prophages and transposable elements but lacking gene sets coding for denitrification, hydrogen oxidization, oxidative phosphorylation, vitamin biosynthesis, cell pH and/or sodium homeostasis, environmental sensing, and motility, indicative of early genome erosion and adaptive evolution toward obligate endosymbiosis. Unexpectedly, a prophage embedded in the HMS1 genome undergoes lytic cycle. Highly expressed ROS scavenger and LexA repressor genes indicate that the tubeworm host likely activates the lysogenic phage into lytic cycle through the SOS response to regulate endosymbiont population and harvest nutrients. Our findings indicate progressive evolution of Sclerolinum endosymbionts toward obligate endosymbiosis and expand the knowledge about phage-symbiont-host interaction in deep-sea tubeworms.

Multiple and diversified transposon lineages contribute to early and recent bivalve genome evolution

BACKGROUND

Transposable elements (TEs) can represent one of the major sources of genomic variation across eukaryotes, providing novel raw materials for species diversification and innovation. While considerable effort has been made to study their evolutionary dynamics across multiple animal clades, molluscs represent a substantially understudied phylum. Here, we take advantage of the recent increase in mollusc genomic resources and adopt an automated TE annotation pipeline combined with a phylogenetic tree-based classification, as well as extensive manual curation efforts, to characterize TE repertories across 27 bivalve genomes with a particular emphasis on DDE/D class II elements, long interspersed nuclear elements (LINEs), and their evolutionary dynamics.

RESULTS

We found class I elements as highly dominant in bivalve genomes, with LINE elements, despite less represented in terms of copy number per genome, being the most common retroposon group covering up to 10% of their genome. We mined 86,488 reverse transcriptases (RVT) containing LINE coming from 12 clades distributed across all known superfamilies and 14,275 class II DDE/D-containing transposons coming from 16 distinct superfamilies. We uncovered a previously underestimated rich and diverse bivalve ancestral transposon complement that could be traced back to their most recent common ancestor that lived ~ 500 Mya. Moreover, we identified multiple instances of lineage-specific emergence and loss of different LINEs and DDE/D lineages with the interesting cases of CR1- Zenon, Proto2, RTE-X, and Academ elements that underwent a bivalve-specific amplification likely associated with their diversification. Finally, we found that this LINE diversity is maintained in extant species by an equally diverse set of long-living and potentially active elements, as suggested by their evolutionary history and transcription profiles in both male and female gonads.

CONCLUSIONS

We found that bivalves host an exceptional diversity of transposons compared to other molluscs. Their LINE complement could mainly follow a "stealth drivers" model of evolution where multiple and diversified families are able to survive and co-exist for a long period of time in the host genome, potentially shaping both recent and early phases of bivalve genome evolution and diversification. Overall, we provide not only the first comparative study of TE evolutionary dynamics in a large but understudied phylum such as Mollusca, but also a reference library for ORF-containing class II DDE/D and LINE elements, which represents an important genomic resource for their identification and characterization in novel genomes.

A draft genome assembly of reef-building octocoral Heliopora coerulea

Coral reefs are under existential threat from climate change and anthropogenic impacts. Genomic studies have enhanced our knowledge of resilience and responses of some coral species to environmental stress, but reference genomes are lacking for many coral species. The blue coral Heliopora is the only reef-building octocoral genus and exhibits optimal growth at a temperature close to the bleaching threshold of scleractinian corals. Local and high-latitude expansions of Heliopora coerulea were reported in the last decade, but little is known about the molecular mechanisms underlying its thermal resistance. We generated a draft genome of H. coerulea with an assembled size of 429.9 Mb, scaffold N50 of 1.42 Mb and BUSCO completeness of 94.9%. The genome contains 239.1 Mb repetitive sequences, 27,108 protein coding genes, 6,225 lncRNAs, and 79 miRNAs. This reference genome provides a valuable resource for in-depth studies on the adaptive mechanisms of corals under climate change and the evolution of skeleton in cnidarian.

A bacterial symbiont in the gill of the marine scallop Argopecten irradians irradians metabolizes dimethylsulfoniopropionate

Microbial lysis of dimethylsulfoniopropionate (DMSP) is a key step in marine organic sulfur cycling and has been recently demonstrated to play an important role in mediating interactions between bacteria, algae, and zooplankton. To date, microbes that have been found to lyse DMSP are largely confined to free-living and surface-attached bacteria. In this study, we report for the first time that a symbiont (termed "Rhodobiaceae bacterium HWgs001") in the gill of the marine scallop Argopecten irradians irradians can lyse and metabolize DMSP. Analysis of 16S rRNA gene sequences suggested that HWgs001 accounted for up to 93% of the gill microbiota. Microscopic observations suggested that HWgs001 lived within the gill tissue. Unlike symbionts of other bivalves, HWgs001 belongs to Alphaproteobacteria rather than Gammaproteobacteria, and no genes for carbon fixation were identified in its small genome. Moreover, HWgs001 was found to possess a dddP gene, responsible for the lysis of DMSP to acrylate. The enzymatic activity of dddP was confirmed using the heterologous expression, and in situ transcription of the gene in scallop gill tissues was demonstrated using reverse-transcription PCR. Together, these results revealed a taxonomically and functionally unique symbiont, which represents the first-documented DMSP-metabolizing symbiont likely to play significant roles in coastal marine ecosystems.

Hologenome analysis reveals independent evolution to chemosymbiosis by deep-sea bivalves

BACKGROUND

Bivalves have independently evolved a variety of symbiotic relationships with chemosynthetic bacteria. These relationships range from endo- to extracellular interactions, making them ideal for studies on symbiosis-related evolution. It is still unclear whether there are universal patterns to symbiosis across bivalves. Here, we investigate the hologenome of an extracellular symbiotic thyasirid clam that represents the early stages of symbiosis evolution.

RESULTS

We present a hologenome of Conchocele bisecta (Bivalvia: Thyasiridae) collected from deep-sea hydrothermal vents with extracellular symbionts, along with related ultrastructural evidence and expression data. Based on ultrastructural and sequencing evidence, only one dominant Thioglobaceae bacteria was densely aggregated in the large bacterial chambers of C. bisecta, and the bacterial genome shows nutritional complementarity and immune interactions with the host. Overall, gene family expansions may contribute to the symbiosis-related phenotypic variations in different bivalves. For instance, convergent expansions of gaseous substrate transport families in the endosymbiotic bivalves are absent in C. bisecta. Compared to endosymbiotic relatives, the thyasirid genome exhibits large-scale expansion in phagocytosis, which may facilitate symbiont digestion and account for extracellular symbiotic phenotypes. We also reveal that distinct immune system evolution, including expansion in lipopolysaccharide scavenging and contraction of IAP (inhibitor of apoptosis protein), may contribute to the different manners of bacterial virulence resistance in C. bisecta.

CONCLUSIONS

Thus, bivalves employ different pathways to adapt to the long-term co-existence with their bacterial symbionts, further highlighting the contribution of stochastic evolution to the independent gain of a symbiotic lifestyle in the lineage.

Interactions among deep-sea mussels and their epibiotic and endosymbiotic chemoautotrophic bacteria: Insights from multi-omics analysis

Endosymbiosis with Gammaproteobacteria is fundamental for the success of bathymodioline mussels in deep-sea chemosynthesis-based ecosystems. However, the recent discovery of Campylobacteria on the gill surfaces of these mussels suggests that these host-bacterial relationships may be more complex than previously thought. Using the cold-seep mussel ( Gigantidas haimaensis) as a model, we explored this host-bacterial system by assembling the host transcriptome and genomes of its epibiotic Campylobacteria and endosymbiotic Gammaproteobacteria and quantifying their gene and protein expression levels. We found that the epibiont applies a sulfur oxidizing (SOX) multienzyme complex with the acquisition of soxB from Gammaproteobacteria for energy production and switched from a reductive tricarboxylic acid (rTCA) cycle to a Calvin-Benson-Bassham (CBB) cycle for carbon assimilation. The host provides metabolic intermediates, inorganic carbon, and thiosulfate to satisfy the materials and energy requirements of the epibiont, but whether the epibiont benefits the host is unclear. The endosymbiont adopts methane oxidation and the ribulose monophosphate pathway (RuMP) for energy production, providing the major source of energy for itself and the host. The host obtains most of its nutrients, such as lysine, glutamine, valine, isoleucine, leucine, histidine, and folate, from the endosymbiont. In addition, host pattern recognition receptors, including toll-like receptors, peptidoglycan recognition proteins, and C-type lectins, may participate in bacterial infection, maintenance, and population regulation. Overall, this study provides insights into the complex host-bacterial relationships that have enabled mussels and bacteria to thrive in deep-sea chemosynthetic ecosystems.

Metal bioaccumulation and transfer in benthic species based on isotopic signatures of individual amino acids in South China Sea cold seep environments

Cold seeps are deep-sea 'oases' with dense and dominant coexisting populations of large mussels and tubeworms under extreme environments. Under such natural source of high metal concentrations, the present study investigated the metal bioaccumulation and transfer with trophic positions in six benthic species by the isotopic δ¹⁵N and δ¹³C signatures in the active Haima cold seep, South China Sea. Comparing the isotopic signatures of bulk-tissue and amino acids by compound-specific isotopic analysis (CSIA-AA), we found that the bulk trophic (TP_B) values in the benthos except mussels were significantly higher than those of CSIA-based TP_Glu-Phe values. The estimated CSIA-based TP_Glu-Phe values showed a relatively compressed food chain with much changeable and unique amino acid isotopic heterogeneity, followed slim tubeworms (1.20)<mussels (1.38)<clams (1.52)<brittle stars (1.82)<giant tubeworms (2.16)<shrimps (2.31). All species accumulated relatively high concentrations of Fe, Zn, Cu, and Cr, especially for Zn in clams. Pearson correlation analysis showed that most metals had no significant relationship between their bioaccumulation and trophic positions, whereas Hg showed a significantly positive bioaccumulation through trophic transfer in such a compressed food chain. Water exposure was a major metal source rather than bacterial assimilation for most metals in the cold seep higher consumers. Hyperaccumulation of specific metals in some tissues of different benthos indicated different metal overflows in the Haima cold seep (As and Ni for tubeworms, Zn and Cd for clam gills, Ag and Cu for mussel gills). This study demonstrated high metal adaptations in different species and stable isotopic characteristics of amino acid metabolism in a natural high metal source of an active deep-sea cold seep, which is important for deep-sea development and environmental protection.

Immune diversity in lophotrochozoans, with a focus on recognition and effector systems

Lophotrochozoa is one of the most species-rich but immunologically poorly explored phyla. Although lack of acquired response in a narrow sense, lophotrochozoans possess various genetic mechanisms that enhance the diversity and specificity of innate immune system. Here, we review the recent advances of comparative immunology studies in lophotrochozoans with focus on immune recognition and effector systems. Haemocytes and coelomocytes are general important yet understudied player. Comparative genomics studies suggest expansion and functional divergence of lophotrochozoan immune reorganization systems is not as "homogeneous and simple" as we thought including the large-scale expansion and molecular divergence of pattern recognition receptors (PRRs) (TLRs, RLRs, lectins, etc.) and signaling adapters (MyD88s etc.), significant domain recombination of immune receptors (RLR, NLRs, lectins, etc.), extensive somatic recombination of fibrinogenrelated proteins (FREPs) in snails. Furthermore, there are repeatedly identified molecular mechanisms that generate immune effector diversity, including high polymorphism of antimicrobial peptides and proteins (AMPs), reactive oxygen and nitrogen species (RONS) and cytokines. Finally, we argue that the next generation omics tools and the recently emerged genome editing technicism will revolutionize our understanding of innate immune system in a comparative immunology perspective.

Ancient homomorphy of molluscan sex chromosomes sustained by reversible sex-biased genes and sex determiner translocation

Contrary to classic theory prediction, sex-chromosome homomorphy is prevalent in the animal kingdom but it is unclear how ancient homomorphic sex chromosomes avoid chromosome-scale degeneration. Molluscs constitute the second largest, Precambrian-originated animal phylum and have ancient, uncharacterized homomorphic sex chromosomes. Here, we profile eight genomes of the bivalve mollusc family of Pectinidae in a phylogenetic context and show 350 million years sex-chromosome homomorphy, which is the oldest known sex-chromosome homomorphy in the animal kingdom, far exceeding the ages of well-known heteromorphic sex chromosomes such as 130-200 million years in mammals, birds and flies. The long-term undifferentiation of molluscan sex chromosomes is potentially sustained by the unexpected intertwined regulation of reversible sex-biased genes, together with the lack of sexual dimorphism and occasional sex chromosome turnover. The pleiotropic constraint of regulation of reversible sex-biased genes is widely present in ancient homomorphic sex chromosomes and might be resolved in heteromorphic sex chromosomes through gene duplication followed by subfunctionalization. The evolutionary dynamics of sex chromosomes suggest a mechanism for 'inheritance' turnover of sex-determining genes that is mediated by translocation of a sex-determining enhancer. On the basis of these findings, we propose an evolutionary model for the long-term preservation of homomorphic sex chromosomes.

A high-quality, haplotype-phased genome reconstruction reveals unexpected haplotype diversity in a pearl oyster

Homologous chromosomes in the diploid genome are thought to contain equivalent genetic information, but this common concept has not been fully verified in animal genomes with high heterozygosity. Here we report a near-complete, haplotype-phased, genome assembly of the pearl oyster, Pinctada fucata, using hi-fidelity (HiFi) long reads and chromosome conformation capture data. This assembly includes 14 pairs of long scaffolds (>38 Mb) corresponding to chromosomes (2n = 28). The accuracy of the assembly, as measured by an analysis of k-mers, is estimated to be 99.99997%. Moreover, the haplotypes contain 95.2% and 95.9%, respectively, complete and single-copy BUSCO genes, demonstrating the high quality of the assembly. Transposons comprise 53.3% of the assembly and are a major contributor to structural variations. Despite overall collinearity between haplotypes, one of the chromosomal scaffolds contains megabase-scale non-syntenic regions, which necessarily have never been detected and resolved in conventional haplotype-merged assemblies. These regions encode expanded gene families of NACHT, DZIP3/hRUL138-like HEPN, and immunoglobulin domains, multiplying the immunity gene repertoire, which we hypothesize is important for the innate immune capability of pearl oysters. The pearl oyster genome provides insight into remarkable haplotype diversity in animals.

Effects of scale worm parasitism on interactions between the symbiotic gill microbiome and gene regulation in deep sea mussel hosts

Diverse adaptations to the challenging deep sea environment are expected to be found across all deep sea organisms. Scale worms Branchipolynoe pettiboneae are believed to adapt to the deep sea environment by parasitizing deep sea mussels; this biotic interaction is one of most known in the deep sea chemosynthetic ecosystem. However, the mechanisms underlying the effects of scale worm parasitism on hosts are unclear. Previous studies have revealed that the microbiota plays an important role in host adaptability. Here, we compared gill-microbiota, gene expression and host-microorganism interactions in a group of deep sea mussels (Gigantidas haimaensis) parasitized by scale worm (PA group) and a no parasitic control group (NPA group). The symbiotic microorganism diversity of the PA group significantly decreased than NPA group, while the relative abundance of chemoautotrophic symbiotic bacteria that provide the host with organic carbon compounds significantly increased in PA. Interestingly, RNA-seq revealed that G. haimaensis hosts responded to B. pettiboneaei parasitism through significant upregulation of protein and lipid anabolism related genes, and that this parasitism may enhance host mussel nutrient anabolism but inhibit the host’s ability to absorb nutrients, thus potentially helping the parasite obtain nutrients from the host. In an integrated analysis of the interactions between changes in the microbiota and host gene dysregulation, we found an agreement between the microbiota and transcriptomic responses to B. pettiboneaei parasitism. Together, our findings provide new insights into the effects of parasite scale worms on changes in symbiotic bacteria and gene expression in deep sea mussel hosts. We explored the potential role of host-microorganism interactions between scale worms and deep sea mussels, and revealed the mechanisms through which scale worm parasitism affects hosts in deep sea chemosynthetic ecosystem.

Evolution of EGF-like and Zona pellucida domains containing shell matrix proteins in mollusks

Several types of shell matrix proteins (SMPs) have been identified in molluskan shells. Their diversity is the consequence of various molecular processes, including domain shuffling and gene duplication. However, the evolutionary origin of most SMPs remains unclear. In this study, we investigated the evolutionary process EGF-like and zona pellucida (ZP) domains containing SMPs. Two types of the proteins (EGF-like protein (EGFL) and EGF-like and ZP domains containing protein (EGFZP)) were found in the pearl oyster, Pinctada fucata. In contrast, only EGFZP was identified in the gastropods. Phylogenetic analysis and genomic arrangement studies showed that EGFL and EGFZP formed a clade in bivalves, and their encoding genes were localized in tandem repeats on the same scaffold. In P. fucata, EGFL genes were expressed in the outer part of mantle epithelial cells are related to the calcitic shell formation. However, in both P. fucata and the limpet Nipponacmea fuscoviridis, EGFZP genes were expressed in the inner part of the mantle epithelial cells are related to aragonitic shell formation. Furthermore, our analysis showed that in P. fucata, the ZP domain interacts with eight SMPs that have various functions in the nacreous shell mineralization. The data suggest that the ZP domain can interact with other SMPs, and EGFL evolution in pterimorph bivalves represents an example of neo-functionalization that involves the acquisition of a novel protein through gene duplication.

Genomic Survey and Resources for the Boring Giant Clam Tridacna crocea

The boring giant clam Tridacna crocea is an evolutionary, ecologically, economically, and culturally important reef-dwelling bivalve targeted by a profitable ornamental fishery in the Indo-Pacific Ocean. In this study, we developed genomic resources for T. crocea. Using low-pass (=low-coverage, ~6×) short read sequencing, this study, for the first time, estimated the genome size, unique genome content, and nuclear repetitive elements, including the 45S rRNA DNA operon, in T. crocea. Furthermore, we tested if the mitochondrial genome can be assembled from RNA sequencing data. The haploid genome size estimated using a k-mer strategy was 1.31–1.39 Gbp, which is well within the range reported before for other members of the family Cardiidae. Unique genome content estimates using different k-mers indicated that nearly a third and probably at least 50% of the genome of T. crocea was composed of repetitive elements. A large portion of repetitive sequences could not be assigned to known repeat element families. Taking into consideration only annotated repetitive elements, the most common were classified as Satellite DNA which were more common than Class I-LINE and Class I-LTR Ty3-gypsy retrotransposon elements. The nuclear ribosomal operon in T. crocea was partially assembled into two contigs, one encoding the complete ssrDNA and 5.8S rDNA unit and a second comprising a partial lsrDNA. A nearly complete mitochondrial genome (92%) was assembled from RNA-seq. These newly developed genomic resources are highly relevant for improving our understanding of the biology of T. crocea and for the development of conservation plans and the fisheries management of this iconic reef-dwelling invertebrate.

Genome of a giant isopod, Bathynomus jamesi, provides insights into body size evolution and adaptation to deep-sea environment

Background

The deep-sea may be regarded as a hostile living environment, due to low temperature, high hydrostatic pressure, and limited food and light. Isopods, a species-rich group of crustaceans, are widely distributed across different environments including the deep sea and as such are a useful model for studying adaptation, migration, and speciation. Similar to other deep-sea organisms, giant isopods have larger body size than their shallow water relatives and have large stomachs and fat bodies presumably to store organic reserves. In order to shed light on the genetic basis of these large crustaceans adapting to the oligotrophic environment of deep-sea, the high-quality genome of a deep-sea giant isopod Bathynomus jamesi was sequenced and assembled.

Results

B. jamesi has a large genome of 5.89 Gb, representing the largest sequenced crustacean genome to date. Its large genome size is mainly attributable to the remarkable proliferation of transposable elements (84%), which may enable high genome plasticity for adaptive evolution. Unlike its relatives with small body size, B. jamesi has expanded gene families related to pathways of thyroid and insulin hormone signaling that potentially contribute to its large body size. Transcriptomic analysis showed that some expanded gene families related to glycolysis and vesicular transport were specifically expressed in its digestive organs. In addition, comparative genomics and gene expression analyses in six tissues suggested that B. jamesi has inefficient lipid degradation, low basal metabolic rate, and bulk food storage, suggesting giant isopods adopt a more efficient mechanism of nutrient absorption, storage, and utilization to provide sustained energy supply for their large body size.

Conclusions

Taken together, the giant isopod genome may provide a valuable resource for understanding body size evolution and adaptation mechanisms of macrobenthic organisms to deep-sea environments.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01302-6.

Comparative transcriptomic analysis of in situ and onboard fixed deep-sea limpets reveals sample preparation-related differences

Precise gene expression reflects the molecular response of deep-sea organisms to their harsh living environments. However, changes in environmental factors during lifting samples from the deep sea to a research vessel can also affect gene expression. By using the transcriptomic approach, we compared the gene expression profiles of the onboard fixed with the in situ fixed samples of the deep-sea limpet Bathyacmaea lactea. Our results revealed that the concomitant stress during conventional deep-sea sampling without RNA in situ fixation greatly influenced the gene expression. Various biological activities, such as cell and tissue structure, lysosomal activity, fluid balance, and unsaturated fatty acid metabolism, were perturbed, suggesting that the sampling stress has exerted systemic impacts on the life of the limpets. These findings clearly illustrate that deep-sea samples without RNA in situ fixation can easily lead to biased results in gene expression analysis, which requires to be appropriately addressed in future studies.

Divergent paths in the evolutionary history of maternally transmitted clam symbionts

Vertical transmission of bacterial endosymbionts is accompanied by virtually irreversible gene loss that results in a progressive reduction in genome size. While the evolutionary processes of genome reduction have been well described in some terrestrial symbioses, they are less understood in marine systems where vertical transmission is rarely observed. The association between deep-sea vesicomyid clams and chemosynthetic Gammaproteobacteria is one example of maternally inherited symbioses in the ocean. Here, we assessed the contributions of drift, recombination and selection to genome evolution in two extant vesicomyid symbiont clades by comparing 15 representative symbiont genomes (1.017-1.586 Mb) to those of closely related bacteria and the hosts' mitochondria. Our analyses suggest that drift is a significant force driving genome evolution in vesicomyid symbionts, though selection and interspecific recombination appear to be critical for maintaining symbiont functional integrity and creating divergent patterns of gene conservation. Notably, the two symbiont clades possess putative functional differences in sulfide physiology, anaerobic respiration and dependency on environmental vitamin B12, which probably reflect adaptations to different ecological habitats available to each symbiont group. Overall, these results contribute to our understanding of the eco-evolutionary processes shaping reductive genome evolution in vertically transmitted symbioses.

Comparative analysis of the Mercenaria mercenaria genome provides insights into the diversity of transposable elements and immune molecules in bivalve mollusks

BACKGROUND

The hard clam Mercenaria mercenaria is a major marine resource along the Atlantic coasts of North America and has been introduced to other continents for resource restoration or aquaculture activities. Significant mortality events have been reported in the species throughout its native range as a result of diseases (microbial infections, leukemia) and acute environmental stress. In this context, the characterization of the hard clam genome can provide highly needed resources to enable basic (e.g., oncogenesis and cancer transmission, adaptation biology) and applied (clam stock enhancement, genomic selection) sciences.

RESULTS

Using a combination of long and short-read sequencing technologies, a 1.86 Gb chromosome-level assembly of the clam genome was generated. The assembly was scaffolded into 19 chromosomes, with an N50 of 83 Mb. Genome annotation yielded 34,728 predicted protein-coding genes, markedly more than the few other members of the Venerida sequenced so far, with coding regions representing only 2% of the assembly. Indeed, more than half of the genome is composed of repeated elements, including transposable elements. Major chromosome rearrangements were detected between this assembly and another recent assembly derived from a genetically segregated clam stock. Comparative analysis of the clam genome allowed the identification of a marked diversification in immune-related proteins, particularly extensive tandem duplications and expansions in tumor necrosis factors (TNFs) and C1q domain-containing proteins, some of which were previously shown to play a role in clam interactions with infectious microbes. The study also generated a comparative repertoire highlighting the diversity and, in some instances, the specificity of LTR-retrotransposons elements, particularly Steamer elements in bivalves.

CONCLUSIONS

The diversity of immune molecules in M. mercenaria may allow this species to cope with varying and complex microbial and environmental landscapes. The repertoire of transposable elements identified in this study, particularly Steamer elements, should be a prime target for the investigation of cancer cell development and transmission among bivalve mollusks.

The genome of the zebra mussel, Dreissena polymorpha: a resource for comparative genomics, invasion genetics, and biocontrol

The zebra mussel, Dreissena polymorpha, continues to spread from its native range in Eurasia to Europe and North America, causing billions of dollars in damage and dramatically altering invaded aquatic ecosystems. Despite these impacts, there are few genomic resources for Dreissena or related bivalves. Although the D. polymorpha genome is highly repetitive, we have used a combination of long-read sequencing and Hi-C-based scaffolding to generate a high-quality chromosome-scale genome assembly. Through comparative analysis and transcriptomics experiments, we have gained insights into processes that likely control the invasive success of zebra mussels, including shell formation, synthesis of byssal threads, and thermal tolerance. We identified multiple intact steamer-like elements, a retrotransposon that has been linked to transmissible cancer in marine clams. We also found that D. polymorpha have an unusual 67 kb mitochondrial genome containing numerous tandem repeats, making it the largest observed in Eumetazoa. Together these findings create a rich resource for invasive species research and control efforts.

Genome-wide identification and characteristic analysis of PGRP gene family in Tegillarca granosa reveals distinct immune response of the invasive pathogen

The peptidoglycan recognition proteins (PGRPs) are conserved innate immune molecular in invertebrates and vertebrates, which play important roles in immune system by recognize the peptidoglycans of bacterial cell walls. Although PGRPs have been extensively characterized in insects, a systematic analysis of PGRPs in bivalves is lacking. In the present study, the phylogenic relationships, gene structures and expression profiles of PGRPs in marine bivalves were analyzed. The results indicated that the most PGRPs of bivalves were predicted to degrade the peptidoglycans and prevent excessive immunostimulation of bacteria. In addition, the results of the present study showed that the protein diversity of PGRPs in most marine bivalves was mainly generated by the alternative splicing of genes, however the alternative splicing of PGRP gene family was absent in Tegillarca granosa. The differences of PGRPs might be related to the genetic and environmental differences of marine bivalves. Spatiotemporal expression profiling in T. granosa suggested that PGRPs play important roles in the immune response of invasive pathogens. The present study describes a comprehensive view of PGRPs in the blood clam T. granosa and provides a foundation for functional characterization of this gene family in innate immune of marine bivalves.

Novel insights on obligate symbiont lifestyle and adaptation to chemosynthetic environment as revealed by the giant tubeworm genome

The mutualism between the giant tubeworm Riftia pachyptila and its endosymbiont Candidatus Endoriftia persephone has been extensively researched over the past 40 years. However, the lack of the host whole genome information has impeded the full comprehension of the genotype/phenotype interface in Riftia. Here we described the high-quality draft genome of Riftia, its complete mitogenome, and tissue-specific transcriptomic data. The Riftia genome presents signs of reductive evolution, with gene family contractions exceeding expansions. Expanded gene families are related to sulphur metabolism, detoxification, anti-oxidative stress, oxygen transport, immune system, and lysosomal digestion, reflecting evolutionary adaptations to the vent environment and endosymbiosis. Despite the derived body plan, the developmental gene repertoire in the gutless tubeworm is extremely conserved with the presence of a near intact and complete Hox cluster. Gene expression analyses establishes that the trophosome is a multi-functional organ marked by intracellular digestion of endosymbionts, storage of excretory products and haematopoietic functions. Overall, the plume and gonad tissues both in contact to the environment harbour highly expressed genes involved with cell cycle, programmed cell death, and immunity indicating a high cell turnover and defence mechanisms against pathogens. We posit that the innate immune system plays a more prominent role into the establishment of the symbiosis during the infection in the larval stage, rather than maintaining the symbiostasis in the trophosome. This genome bridges four decades of physiological research in Riftia, whilst simultaneously provides new insights into the development, whole organism functions and evolution in the giant tubeworm.

Microbiomes and Obligate Symbiosis of Deep-Sea Animals

Microbial communities associated with deep-sea animals are critical to the establishment of novel biological communities in unusual environments. Over the past few decades, rapid exploration of the deep sea has enabled the discovery of novel microbial communities, some of which form symbiotic relationships with animal hosts. Symbiosis in the deep sea changes host physiology, behavior, ecology, and evolution over time and space. Symbiont diversity within a host is often aligned with diverse metabolic pathways that broaden the environmental niche for the animal host. In this review, we focus on microbiomes and obligate symbionts found in different deep-sea habitats and how they facilitate survival of the organisms that live in these environments. In addition, we discuss factors that govern microbiome diversity, host specificity, and biogeography in the deep sea. Finally, we highlight the current limitations of microbiome research and draw a road map for future directions to advance our knowledge of microbiomes in the deep sea. Expected final online publication date for the Annual Review of Animal Biosciences, Volume 10 is February 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Life in the Dark: Phylogenetic and Physiological Diversity of Chemosynthetic Symbioses

Possibly the last discovery of a previously unknown major ecosystem on Earth was made just over half a century ago, when researchers found teaming communities of animals flourishing two and a half kilometers below the ocean surface at hydrothermal vents. We now know that these highly productive ecosystems are based on nutritional symbioses between chemosynthetic bacteria and eukaryotes and that these chemosymbioses are ubiquitous in both deep-sea and shallow-water environments. The symbionts are primary producers that gain energy from the oxidation of reduced compounds, such as sulfide and methane, to fix carbon dioxide or methane into biomass to feed their hosts. This review outlines how the symbiotic partners have adapted to living together. We first focus on the phylogenetic and metabolic diversity of these symbioses and then highlight selected research directions that could advance our understanding of the processes that shaped the evolutionary and ecological success of these associations. Expected final online publication date for the Annual Review of Microbiology, Volume 75 is October 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

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.

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.

Environment and Diet Influence the Bacterial Microbiome of Ambigolimax valentianus, an Invasive Slug in California

Ambigolimax valentianus is an invasive European terrestrial gastropod distributed throughout California. It is a serious pest of gardens, plant nurseries, and greenhouses. We evaluated the bacterial microbiome of whole slugs to capture a more detailed picture of bacterial diversity and composition in this host. We concentrated on the influences of diet and environment on the Ambigolimax valentianus core bacterial microbiome as a starting point for obtaining valuable information to aid in future slug microbiome studies. Ambigolimax valentianus were collected from two environments (gardens or reared from eggs in a laboratory). DNA from whole slugs were extracted and next-generation 16S rRNA gene sequencing was performed. Slug microbiomes differed between environmental sources (garden- vs. lab-reared) and were influenced by a sterile diet. Lab-reared slugs fed an unsterile diet harbored greater bacterial species than garden-reared slugs. A small core microbiome was present that was shared across all slug treatments. This is consistent with our hypothesis that a core microbiome is present and will not change due to these treatments. Findings from this study will help elucidate the impacts of slug-assisted bacterial dispersal on soils and plants, while providing valuable information about the slug microbiome for potential integrated pest research applications.

Genomic insights into the sessile life and biofouling of barnacles (Crustacea: Cirripedia)

Members of the infraclass Cirripedia, commonly called barnacles, are unique among the subphylum Crustacea in that they exhibit a biphasic life cycle with a planktonic larval stage and a sessile adult stage. Understanding their unique sessile life and mechanisms of attachment are hampered by the lack of genomic resources. Here, we present a 746 Mb genome assembly of Lepas anserifera – the first sequenced stalked barnacle genome. We estimate that Cirripedia first arose ~495 million years ago (MYA) and further diversified since Mesozoic. A demographic analysis revealed remarkable population changes of the barnacle in relation to sea-level fluctuations in the last 2 MYA. Comparative genomic analyses revealed the expansion of a number of developmental related genes families in barnacle genomes, such as Br–C, PCP20 and Lola, which are potentially important for the evolution of metamorphosis, cuticle development and central nervous system. Phylogenetic analysis and tissue expression profiling showed the possible roles of gene duplication, functional diversification and co-option in shaping the genomic evolution of barnacles. Overall, our study provides not only a valuable draft genome for comparative genomic analysis of crustacean evolution, but also facilitates studies of biofouling control.

A High-Quality Reference Genome for a Parasitic Bivalve with Doubly Uniparental Inheritance (Bivalvia: Unionida)

From a genomics perspective, bivalves (Mollusca: Bivalvia) have been poorly explored with the exception for those of high economic value. The bivalve order Unionida, or freshwater mussels, has been of interest in recent genomic studies due to their unique mitochondrial biology and peculiar life cycle. However, genomic studies have been hindered by the lack of a high-quality reference genome. Here, I present a genome assembly of Potamilus streckersoni using Pacific Bioscience single-molecule real-time long reads and 10X Genomics-linked read sequencing. Further, I use RNA sequencing from multiple tissue types and life stages to annotate the reference genome. The final assembly was far superior to any previously published freshwater mussel genome and was represented by 2,368 scaffolds (2,472 contigs) and 1,776,755,624 bp, with a scaffold N50 of 2,051,244 bp. A high proportion of the assembly was comprised of repetitive elements (51.03%), aligning with genomic characteristics of other bivalves. The functional annotation returned 52,407 gene models (41,065 protein, 11,342 tRNAs), which was concordant with the estimated number of genes in other freshwater mussel species. This genetic resource, along with future studies developing high-quality genome assemblies and annotations, will be integral toward unraveling the genomic bases of ecologically and evolutionarily important traits in this hyper-diverse group.