Organ transcriptomes of the lucinid clam Loripes orbiculatus (Poli, 1791) provide insights into their specialised roles in the biology of a chemosymbiotic bivalve

Comparative transcriptomics revealed parallel evolution and innovation of photosymbiosis molecular mechanisms in a marine bivalve

Photosymbioses between heterotrophic hosts and autotrophic symbionts are evolutionarily prevalent and ecologically significant. However, the molecular mechanisms behind such symbioses remain less elucidated, which hinders our understanding of their origin and adaptive evolution. This study compared gene expression patterns in a photosymbiotic bivalve (Fragum sueziense) and a closely related non-symbiotic species (Trigoniocardia granifera) under different light conditions to detect potential molecular pathways involved in mollusc photosymbiosis. We discovered that the presence of algal symbionts greatly impacted host gene expression in symbiont-containing tissues. We found that the host immune functions were suppressed under normal light compared with those in the dark. In addition, we found that cilia in the symbiont-containing tissues play important roles in symbiont regulation or photoreception. Interestingly, many potential photosymbiosis genes could not be annotated or do not exhibit orthologues in T. granifera transcriptomes, indicating unique molecular functions in photosymbiotic bivalves. Overall, we found both novel and known molecular mechanisms involved in animal-algal photosymbiosis within bivalves. Given that many of the molecular pathways are shared among distantly related host lineages, such as molluscs and cnidarians, it indicates that parallel and/or convergent evolution is instrumental in shaping host-symbiont interactions and responses in these organisms.

Metabolically-versatile Ca. Thiodiazotropha symbionts of the deep-sea lucinid clam Lucinoma kazani have the genetic potential to fix nitrogen

Lucinid clams are one of the most diverse and widespread symbiont-bearing animal groups in both shallow and deep-sea chemosynthetic habitats. Lucinids harbor Ca. Thiodiazotropha symbionts that can oxidize inorganic and organic substrates such as hydrogen sulfide and formate to gain energy. The interplay between these key metabolic functions, nutrient uptake and biotic interactions in Ca. Thiodiazotropha is not fully understood. We collected Lucinoma kazani individuals from next to a deep-sea brine pool in the eastern Mediterranean Sea, at a depth of 1150 m and used Oxford Nanopore and Illumina sequencing to obtain high-quality genomes of their Ca. Thiodiazotropha gloverae symbiont. The genomes served as the basis for transcriptomic and proteomic analyses to characterize the in situ gene expression, metabolism and physiology of the symbionts. We found genes needed for N2 fixation in the deep-sea symbiont's genome, which, to date, were only found in shallow-water Ca. Thiodiazotropha. However, we did not detect the expression of these genes and thus the potential role of nitrogen fixation in this symbiosis remains to be determined. We also found the high expression of carbon fixation and sulfur oxidation genes, which indicate chemolithoautotrophy as the key physiology of Ca. Thiodiazotropha. However, we also detected the expression of pathways for using methanol and formate as energy sources. Our findings highlight the key traits these microbes maintain to support the nutrition of their hosts and interact with them.

Molecular characterization of MyD88 as a potential biomarker for pesticide-induced stress in Bombyx mori

The widespread use of pesticides hampers the immune system of non-target organisms, however, there is a lack of common biomarkers to detect such effects. Myeloid differentiation primary response factor 88 (MyD88) is a crucial junction protein in the Toll-like receptor signaling pathway, which plays an important role in the inflammatory response. In this study, we investigated MyD88 as a potential biomarker for pesticide-induced stress. Phylogenetic analysis revealed that MyD88 was a conserved protein in the evolution of vertebrates and invertebrates. MyD88s usually have death domain (DD) and Toll/interleukin-1 receptor (TIR) domain. Bombyx mori (B. mori) is an important economic insect that is sensitive to toxic substances. We found microbial pesticides enhanced the expression level of MyD88 in B. mori. Transcriptome analysis demonstrated that MyD88 expression level was increased in the fatbody after dinotefuran exposure, a third-generation neonicotinoid pesticide. Moreover, the expression of MyD88 was upregulated in fatbody and midgut by imidacloprid, a first-generation neonicotinoid pesticide. Additionally, insect growth regulator (IGR) pesticides, such as methoprene and fenoxycarb, could induce MyD88 expression in the fatbody of B. mori. These results indicated that MyD88 is a potential biomarker for pesticide-induced stress in B. mori. This study provides novel insights into screening common biomarkers for multiple pesticide stresses and important implications for the development of more sustainable pest management strategies.

Specification of hemocyte subpopulations based on immune-related activities and the production of the agglutinin MkC1qDC in the bivalve Modiolus kurilensis

Bivalves, such as Modiolus are used as indicator organisms to monitor the state of the marine environment. Even though hemocytes are known to play a key role in the adaptive and protective mechanisms of bivalves, these cells are poorly studied in horse-mussel Modiolus kurilensis. In this paper, we present classification of horse-mussel hemocytes based on their immune functions, including the production of specific immune-related molecules, as well as their morphological composition after isolation by density gradient centrifugation. An effective fractionation protocol was adapted to separate four hemocyte subpopulations with distinct morphofunctional profiles. First subpopulation consisted of small under-differentiated hemoblasts (2.20 ± 0.85%) with a bromodeoxyuridine positive nucleus, and did not show any immune reactivity. Second was represented by agranulocytes (24.11 ± 2.40%), with evenly filled cytoplasm containing a well-developed protein-synthesizing apparatus, polysomes, smooth endoplasmic reticulum and mitochondria, and positively stained for myeloperoxidase, acidic proteins, glycogen and neutral polysaccharides. Third subpopulation consisted of eosinophilic granulocytes (62.64 ± 9.32%) that contained the largest number of lysosomes, peroxisomes and vesicles with contents of different density, and showed the highest phosphatase, reactive oxygen species (ROS) and phagocytic activities. Lastly, fourth group, basophilic granulocytes (14.21 ± 0.34%), are main producers of lectin-like protein MkC1qDC, recently discovered in M. kurilensis and characterized by pronounced antibacterial and anticancer activity. These cells characterized by intracytoplasmic of the MkC1qDC localization, forming granule-like bodies visualized with specific antibody. Both granulocytes and agranulocytes showed phagocytic activity and ROS production, and these reactions were more pronounced for eosinophilic granulocytes, suggesting that this group is the key element of the cell-mediated immune response of M. kurilensis. Our results support a concept of bivalve's hemocyte specification with distinct phenotypes.

Impact of adaptive filtering on power and false discovery rate in RNA-seq experiments

Background

In RNA-sequencing studies a large number of hypothesis tests are performed to compare the differential expression of genes between several conditions. Filtering has been proposed to remove candidate genes with a low expression level which may not be relevant and have little or no chance of showing a difference between conditions. This step may reduce the multiple testing burden and increase power.

Results

We show in a simulation study that filtering can lead to some increase in power for RNA-sequencing data, too aggressive filtering, however, can lead to a decline. No uniformly optimal filter in terms of power exists. Depending on the scenario different filters may be optimal. We propose an adaptive filtering strategy which selects one of several filters to maximise the number of rejections. No additional adjustment for multiplicity has to be included, but a rule has to be considered if the number of rejections is too small.

Conclusions

For a large range of simulation scenarios, the adaptive filter maximises the power while the simulated False Discovery Rate is bounded by the pre-defined significance level. Using the adaptive filter, it is not necessary to pre-specify a single individual filtering method optimised for a specific scenario.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04928-z.

Differential regulation of degradation and immune pathways underlies adaptation of the ectosymbiotic nematode Laxus oneistus to oxic-anoxic interfaces

Eukaryotes may experience oxygen deprivation under both physiological and pathological conditions. Because oxygen shortage leads to a reduction in cellular energy production, all eukaryotes studied so far conserve energy by suppressing their metabolism. However, the molecular physiology of animals that naturally and repeatedly experience anoxia is underexplored. One such animal is the marine nematode Laxus oneistus. It thrives, invariably coated by its sulfur-oxidizing symbiont Candidatus Thiosymbion oneisti, in anoxic sulfidic or hypoxic sand. Here, transcriptomics and proteomics showed that, whether in anoxia or not, L. oneistus mostly expressed genes involved in ubiquitination, energy generation, oxidative stress response, immune response, development, and translation. Importantly, ubiquitination genes were also highly expressed when the nematode was subjected to anoxic sulfidic conditions, together with genes involved in autophagy, detoxification and ribosome biogenesis. We hypothesize that these degradation pathways were induced to recycle damaged cellular components (mitochondria) and misfolded proteins into nutrients. Remarkably, when L. oneistus was subjected to anoxic sulfidic conditions, lectin and mucin genes were also upregulated, potentially to promote the attachment of its thiotrophic symbiont. Furthermore, the nematode appeared to survive oxygen deprivation by using an alternative electron carrier (rhodoquinone) and acceptor (fumarate), to rewire the electron transfer chain. On the other hand, under hypoxia, genes involved in costly processes (e.g., amino acid biosynthesis, development, feeding, mating) were upregulated, together with the worm's Toll-like innate immunity pathway and several immune effectors (e.g., bactericidal/permeability-increasing proteins, fungicides). In conclusion, we hypothesize that, in anoxic sulfidic sand, L. oneistus upregulates degradation processes, rewires the oxidative phosphorylation and reinforces its coat of bacterial sulfur-oxidizers. In upper sand layers, instead, it appears to produce broad-range antimicrobials and to exploit oxygen for biosynthesis and development.

Global biogeography of chemosynthetic symbionts reveals both localized and globally distributed symbiont groups

In the ocean, most hosts acquire their symbionts from the environment. Due to the immense spatial scales involved, our understanding of the biogeography of hosts and symbionts in marine systems is patchy, although this knowledge is essential for understanding fundamental aspects of symbiosis such as host-symbiont specificity and evolution. Lucinidae is the most species-rich and widely distributed family of marine bivalves hosting autotrophic bacterial endosymbionts. Previous molecular surveys identified location-specific symbiont types that "promiscuously" form associations with multiple divergent cooccurring host species. This flexibility of host-microbe pairings is thought to underpin their global success, as it allows hosts to form associations with locally adapted symbionts. We used metagenomics to investigate the biodiversity, functional variability, and genetic exchange among the endosymbionts of 12 lucinid host species from across the globe. We report a cosmopolitan symbiont species, Candidatus Thiodiazotropha taylori, associated with multiple lucinid host species. Ca. T. taylori has achieved more success at dispersal and establishing symbioses with lucinids than any other symbiont described thus far. This discovery challenges our understanding of symbiont dispersal and location-specific colonization and suggests both symbiont and host flexibility underpin the ecological and evolutionary success of the lucinid symbiosis.

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

Endosymbiosis with chemosynthetic bacteria has enabled many deep-sea invertebrates to thrive at hydrothermal vents and cold seeps, but most previous studies on this mutualism have focused on the bacteria only. Vesicomyid clams dominate global deep-sea chemosynthesis-based ecosystems. They differ from most deep-sea symbiotic animals in passing their symbionts from parent to offspring, enabling intricate coevolution between the host and the symbiont. Here, we sequenced the genomes of the clam Archivesica marissinica (Bivalvia: Vesicomyidae) and its bacterial symbiont to understand the genomic/metabolic integration behind this symbiosis. At 1.52 Gb, the clam genome encodes 28 genes horizontally transferred from bacteria, a large number of pseudogenes and transposable elements whose massive expansion corresponded to the timing of the rise and subsequent divergence of symbiont-bearing vesicomyids. The genome exhibits gene family expansion in cellular processes that likely facilitate chemoautotrophy, including gas delivery to support energy and carbon production, metabolite exchange with the symbiont, and regulation of the bacteriocyte population. Contraction in cellulase genes is likely adaptive to the shift from phytoplankton-derived to bacteria-based food. It also shows contraction in bacterial recognition gene families, indicative of suppressed immune response to the endosymbiont. The gammaproteobacterium endosymbiont has a reduced genome of 1.03 Mb but retains complete pathways for sulfur oxidation, carbon fixation, and biosynthesis of 20 common amino acids, indicating the host’s high dependence on the symbiont for nutrition. Overall, the host–symbiont genomes show not only tight metabolic complementarity but also distinct signatures of coevolution allowing the vesicomyids to thrive in chemosynthesis-based ecosystems.

MicroRNA-Mediated Regulation of Initial Host Responses in a Symbiotic Organ

One of the most important events in an animal's life history is the initial colonization by its microbial symbionts, yet little is known about this event's immediate impacts on the extent of host gene expression or the molecular mechanisms controlling it. MicroRNAs (miRNAs) are short, noncoding RNAs that bind to target mRNAs, rapidly shaping gene expression by posttranscriptional control of mRNA translation and decay. Here, we show that, in the experimentally tractable binary squid-vibrio symbiosis, colonization of the light organ induces extensive changes in the miRNA transcriptome. Examination of the squid genome revealed the presence of evolutionarily conserved genes encoding elements essential for the production and processing of miRNAs. At 24 h postcolonization, 215 host miRNAs were detected in the light organ, 26 of which were differentially expressed in response to the symbionts. A functional enrichment analysis of genes potentially targeted by downregulation of certain miRNAs at the initiation of symbiosis revealed two major gene ontology (GO) term categories, neurodevelopment and tissue remodeling. This symbiont-induced downregulation is predicted to promote these activities in host tissues and is consistent with the well-described tissue remodeling that occurs at the onset of the association. Conversely, predicted targets of upregulated miRNAs, including the production of mucus, are consistent with attenuation of immune responses by symbiosis. Taken together, our data provide evidence that, at the onset of symbiosis, host miRNAs in the light organ drive alterations in gene expression that (i) orchestrate the symbiont-induced development of host tissues, and (ii) facilitate the partnership by dampening the immune response.IMPORTANCE Animals often acquire their microbiome from the environment at each generation, making the initial interaction of the partners a critical event in the establishment and development of a stable, healthy symbiosis. However, the molecular nature of these earliest interactions is generally difficult to study and poorly understood. We report that, during the initial 24 h of the squid-vibrio association, a differential expression of host miRNAs is triggered by the presence of the microbial partner. Predicted mRNA targets of these miRNAs were associated with regulatory networks that drive tissue remodeling and immune suppression, two major symbiosis-induced developmental outcomes in this and many other associations. These results implicate regulation by miRNAs as key to orchestrating the critical transcriptional responses that occur very early during the establishment of a symbiosis. Animals with more complex microbiota may have similar miRNA-driven responses as their association is initiated, supporting an evolutionary conservation of symbiosis-induced developmental mechanisms.

Gill microbiome structure and function in the chemosymbiotic coastal lucinid Stewartia floridana

Lucinid bivalves harbor environmentally acquired, chemosynthetic, gammaproteobacterial gill endosymbionts. Lucinid gill microbiomes, which may contain other gammaproteobacterial and/or spirochete taxa, remain under-sampled. To understand inter-host variability of the lucinid gill microbiome, specifically in the bacterial communities, we analyzed the microbiome content of Stewartia floridana collected from Florida. Sampled gills contained a monospecific gammaproteobacterial endosymbiont expressing lithoautotrophic, mixotrophic, diazotrophic and C1 compound oxidation-related functions previously characterized in similar lucinid species. Another low-abundance Spirochaeta-like species in ∼72% of the sampled gills was most closely related to Spirochaeta-like species in another lucinid Phacoides pectinatus and formed a clade with known marine Spirochaeta symbionts. The spirochete expressed genes were involved in heterotrophy and the transport of sugars, amino acids, peptides and other substrates. Few muscular and neurofilament genes from the host and none from the gammaproteobacterial and spirochete symbionts were differentially expressed among quadrats predominantly covered with seagrass species or 80% bare sand. Our results suggest that spirochetes are facultatively associated with S. floridana, with potential scavenging and nutrient cycling roles. Expressed stress- and defense-related functions in the host and symbionts also suggest species-species communications, which highlight the need for further study of the interactions among lucinid hosts, their microbiomes and their environment.

The symbiotic 'all-rounders': Partnerships between marine animals and chemosynthetic nitrogen-fixing bacteria

Nitrogen fixation is a widespread metabolic trait in certain types of microorganisms called diazotrophs. Bioavailable nitrogen is limited in various habitats on land and in the sea, and accordingly, a range of plant, animal, and single-celled eukaryotes have evolved symbioses with diverse diazotrophic bacteria, with enormous economic and ecological benefits. Until recently, all known nitrogen-fixing symbionts were heterotrophs such as nodulating rhizobia, or photoautotrophs such as cyanobacteria. In 2016, the first chemoautotrophic nitrogen-fixing symbionts were discovered in a common family of marine clams, the Lucinidae. Chemosynthetic nitrogen-fixing symbionts use the chemical energy stored in reduced sulfur compounds to power carbon and nitrogen fixation, making them metabolic 'all-rounders' with multiple functions in the symbiosis. This distinguishes them from heterotrophic symbionts that require a source of carbon from their host, and their chemosynthetic metabolism distinguishes them from photoautotrophic symbionts that produce oxygen, a potent inhibitor of nitrogenase. In this review, we consider evolutionary aspects of this discovery, by comparing strategies that have evolved for hosting intracellular nitrogen-fixing symbionts in plants and animals. The symbiosis between lucinid clams and chemosynthetic nitrogen-fixing bacteria also has important ecological impacts, as they form a nested symbiosis with endangered marine seagrasses. Notably, nitrogen fixation by lucinid symbionts may help support seagrass health by providing a source of nitrogen in seagrass habitats. These discoveries were enabled by new techniques for understanding the activity of microbial populations in natural environments. However, an animal (or plant) host represents a diverse landscape of microbial niches due to its structural, chemical, immune and behavioural properties. In future, methods that resolve microbial activity at the single cell level will provide radical new insights into the regulation of nitrogen fixation in chemosynthetic symbionts, shedding new light on the evolution of nitrogen-fixing symbioses in contrasting hosts and environments.