Niche-dependent sponge hologenome expression profiles and the host-microbes interplay: a case of the hawaiian demosponge Mycale Grandis.

BACKGROUND: Most researches on sponge holobionts focus primarily on symbiotic microbes, yet data at the level of the sponge hologenome are still relatively scarce. Understanding of the sponge host and its microbial gene expression profiles and the host-microbes interplay in different niches represents a key aspect of sponge hologenome. Using the Hawaiian demosponge Mycale grandis in different niches as a model, i.e. on rocks, on the surface of coral Porites compressa, under alga Gracilaria salicornia, we compared the bacterial and fungal community structure, functional gene diversity, expression pattern and the host transcriptome by integrating open-format (deep sequencing) and closed-format (GeoChip microarray) high-throughput techniques. RESULTS: Little inter-niche variation in bacterial and fungal phylogenetic diversity was detected for M. grandis in different niches, but a clear niche-dependent variability in the functional gene diversity and expression pattern of M. grandis host and its symbiotic microbiota was uncovered by GeoChip microarray and transcriptome analyses. Particularly, sponge host genes related to innate immunity and microbial recognition showed a strong correlation with the microbial symbionts' functional gene diversity and transcriptional richness in different niches. The cross-niche variability with respect to the symbiont functional gene diversity and the transcriptional richness of M. grandis holobiont putatively reflects the interplay of niche-specific selective pressure and the symbiont functional diversity. CONCLUSIONS: Niche-dependent gene expression profiles of M. grandis hologenome and the host-microbes interplay were suggested though little inter-niche variation in bacterial and fungal diversity was detected, particularly the sponge innate immunity was found to be closely related to the symbiotic microbes. Altogether, these findings provide novel insights into the black box of one sponge holobiont in different niches at the hologenome level.

Clownfish larvae exhibit faster growth, higher metabolic rates and altered gene expression under future ocean warming.

Increasing ocean temperatures have been demonstrated to have a range of negative impacts on coral reef fishes. However, despite a wealth of studies of juvenile/adult reef fish, studies of how early developmental stages respond to ocean warming are limited. As overall population persistence is influenced by the development of early life stages, detailed studies of larval responses to ocean warming are essential. Here, in an aquaria-based study we investigate how temperatures associated with future warming and present-day marine heatwaves (+3 °C) impact the growth, metabolic rate, and transcriptome of 6 discrete developmental stages of clownfish larvae (Amphiprion ocellaris). A total of 6 clutches of larvae were assessed, with 897 larvae imaged, 262 larvae undergoing metabolic testing and 108 larvae subject to transcriptome sequencing. Our results show that larvae reared at +3 °C grow and develop significantly faster and exhibit higher metabolic rates than those in control conditions. Finally, we highlight the molecular mechanisms underpinning the response of larvae from different developmental stages to higher temperatures, with genes associated with metabolism, neurotransmission, heat stress and epigenetic reprogramming differentially expressed at +3 °C. Overall, these results indicate that clownfish development could be altered under future warming, with developmental rate, metabolic rate, and gene expression all affected. Such changes may lead to altered larval dispersal, changes in settlement time and increased energetic costs.

Transcriptome assemblies of two deep-sea octocorals Calyptrophora lyra and Chrysogorgia stellata from West Pacific seamount, Godin Guyot.

This is the first report of the transcriptome assemblies of the deep-sea octocorals Calyptrophora lyra and Chrysogorgia stellata, which were collected in a survey of the West Pacific seamounts area. We sequenced the transcriptomes of C. lyra and C. stellata using the Illumina NovaSeq 6000 System. De novo assembly and analysis of the coding regions predicted 193,796 unigenes from the total 116,441,796 reads of C. lyra and 235,513 unigenes from the total 122,031,866 reads of C. stellata. Our data are a valuable resource with which to understand the ecological and biological characteristics of the West Pacific deep-sea corals. The data will also contribute to the study of deep-sea environments as extreme and limited habitats and provide direction for future research and further insight into the organismal responses of deep-sea corals to environmental changes.

The chromosome-scale genome assembly of the yellowtail clownfish Amphiprion clarkii provides insights into the melanic pigmentation of anemonefish.

Anemonefish are an emerging group of model organisms for studying genetic, ecological, evolutionary, and developmental traits of coral reef fish. The yellowtail clownfish Amphiprion clarkii possesses species-specific characteristics such as inter-species co-habitation, high intra-species color variation, no anemone specificity, and a broad geographic distribution, that can increase our understanding of anemonefish evolutionary history, behavioral strategies, fish-anemone symbiosis, and color pattern evolution. Despite its position as an emerging model species, the genome of A. clarkii is yet to be published. Using PacBio long-read sequencing and Hi-C chromatin capture technology, we generated a high-quality chromosome-scale genome assembly initially comprised of 1,840 contigs with an N50 of 1,203,211 bp. These contigs were successfully anchored into 24 chromosomes of 843,582,782 bp and annotated with 25,050 protein-coding genes encompassing 97.0% of conserved actinopterygian genes, making the quality and completeness of this genome the highest among all published anemonefish genomes to date. Transcriptomic analysis identified tissue-specific gene expression patterns, with the brain and optic lobe having the largest number of expressed genes. Further analyses revealed higher copy numbers of erbb3b (a gene involved in melanocyte development) in A. clarkii compared with other anemonefish, thus suggesting a possible link between erbb3b and the natural melanism polymorphism observed in A. clarkii. The publication of this high-quality genome, along with A. clarkii's many unique traits, position this species as an ideal model organism for addressing scientific questions across a range of disciplines.

A chromosome-scale genome assembly of the false clownfish, Amphiprion ocellaris.

The false clownfish Amphiprion ocellaris is a popular fish species and an emerging model organism for studying the ecology, evolution, adaptation, and developmental biology of reef fishes. Despite this, high-quality genomic resources for this species are scarce, hindering advanced genomic analyses. Leveraging the power of PacBio long-read sequencing and Hi-C chromosome conformation capture techniques, we constructed a high-quality chromosome-scale genome assembly for the clownfish A. ocellaris. The initial genome assembly comprised of 1,551 contigs of 861.42 Mb, with an N50 of 863.85 kb. Hi-C scaffolding of the genome resulted in 24 chromosomes containing 856.61 Mb. The genome was annotated with 26,797 protein-coding genes and had 96.62% completeness of conserved actinopterygian genes, making this genome the most complete and high quality among published anemonefish genomes. Transcriptomic analysis identified tissue-specific gene expression patterns, with the brain and optic lobe having the largest number of expressed genes. Further, comparative genomic analysis revealed 91 genome elements conserved only in A. ocellaris and its sister species Amphiprion percula, and not in other anemonefish species. These elements are close to genes that are involved in various nervous system functions and exhibited distinct expression patterns in brain tissue, potentially highlighting the genetic toolkits involved in lineage-specific divergence and behaviors of the clownfish branch. Overall, our study provides the highest quality A. ocellaris genome assembly and annotation to date, whilst also providing a valuable resource for understanding the ecology and evolution of reef fishes.

First transcriptome assembly of a newly discovered vent mussel, Gigantidas vrijenhoeki, at Onnuri Vent Field on the northern Central Indian Ridge.

This is the first report of a transcriptome assembly of a newly discovered hydrothermal vent mussel, Gigantidas vrijenhoeki (Bivalvia: Mytilidae), on the Central Indian Ridge. Gigantidas vrijenhoeki was identified from material collected at the newly discovered Onnuri Vent Field (OVF) on the Central Indian Ridge in 2018, and was reported as a new species, distinct from another dominant hydrothermal vent mussel, Bathymodiolus marisindicus, in 2020. We sequenced the transcriptome of G. vrijenhoeki using the Illumina HiSeq X System. De novo assembly and analysis of the coding regions predicted 25,405 genes, 84.76% of which was annotated by public databases. The transcriptome of G. vrijenhoeki will be a valuable resource in studying the ecological and biological characteristics of this new species, which is distinct from other deep-sea mussels. These data should also support the investigation of the relationship between the environmental conditions of hydrothermal vents and the unique distribution of G. vrijenhoeki in the OVF of the Central Indian Ridge.

The Draft Genome of an Octocoral, Dendronephthya gigantea.

Coral reefs composed of stony corals are threatened by global marine environmental changes. However, soft coral communities of octocorallian species, appear more resilient. The genomes of several cnidarians species have been published, including from stony corals, sea anemones, and hydra. To fill the phylogenetic gap for octocoral species of cnidarians, we sequenced the octocoral, Dendronephthya gigantea, a nonsymbiotic soft coral, commonly known as the carnation coral. The D. gigantea genome size is ∼276 Mb. A high-quality genome assembly was constructed from PacBio long reads (29.85 Gb with 108× coverage) and Illumina short paired-end reads (35.54 Gb with 128× coverage) resulting in the highest N50 value (1.4 Mb) reported thus far among cnidarian genomes. About 12% of the genome is repetitive elements and contained 28,879 predicted protein-coding genes. This gene set is composed of 94% complete BUSCO ortholog benchmark genes, which is the second highest value among the cnidarians, indicating high quality. Based on molecular phylogenetic analysis, octocoral and hexacoral divergence times were estimated at 544 MYA. There is a clear difference in Hox gene composition between these species: unlike hexacorals, the Antp superclass Evx gene was absent in D. gigantea. Here, we present the first genome assembly of a nonsymbiotic octocoral, D. gigantea to aid in the comparative genomic analysis of cnidarians, including stony and soft corals, both symbiotic and nonsymbiotic. The D. gigantea genome may also provide clues to mechanisms of differential coping between the soft and stony corals in response to scenarios of global warming.

Phenotypic and molecular consequences of stepwise temperature increase across generations in a coral reef fish.

Global warming will have far-reaching consequences for marine species over coming decades, yet the magnitude of these effects may depend on the rate of warming across generations. Recent experiments show coral reef fishes can compensate the metabolic challenges of elevated temperature when warm conditions are maintained across generations. However, the effects of a gradual temperature increase across generations remain unknown. In the present study, we analysed metabolic and molecular traits in the damselfish Acanthochromis polyacanthus that were exposed to +1.5°C in the first generation and +3.0°C in the second (Step +3.0°C). This treatment of stepwise warming was compared to fish reared at current-day temperatures (Control), second-generation fish of control parents reared at +3.0°C (Developmental +3.0°C) and fish exposed to elevated temperatures for two generations (Transgenerational +1.5°C and Transgenerational +3.0°C). Hepatosomatic index, oxygen consumption and liver gene expression were compared in second-generation fish of the multiple treatments. Hepatosomatic index increased in fish that developed at +3.0°C, regardless of the parental temperature. Routine oxygen consumption of Step +3.0°C fish was significantly higher than Control; however, their aerobic scope recovered to the same level as Control fish. Step +3.0°C fish exhibited significant upregulation of genes related to mitochondrial activity and energy production, which could be associated with their increased metabolic rates. These results indicate that restoration of aerobic scope is possible when fish experience gradual thermal increase across multiple generations, but the metabolic and molecular responses are different from fish reared at the same elevated thermal conditions in successive generations.

Genome-Based Analyses of Six Hexacorallian Species Reject the "Naked Coral" Hypothesis.

Scleractinian corals are the foundation species of the coral-reef ecosystem. Their calcium carbonate skeletons form extensive structures that are home to millions of species, making coral reefs one of the most diverse ecosystems of our planet. However, our understanding of how reef-building corals have evolved the ability to calcify and become the ecosystem builders they are today is hampered by uncertain relationships within their subclass Hexacorallia. Corallimorpharians have been proposed to originate from a complex scleractinian ancestor that lost the ability to calcify in response to increasing ocean acidification, suggesting the possibility for corals to lose and gain the ability to calcify in response to increasing ocean acidification. Here, we employed a phylogenomic approach using whole-genome data from six hexacorallian species to resolve the evolutionary relationship between reef-building corals and their noncalcifying relatives. Phylogenetic analysis based on 1,421 single-copy orthologs, as well as gene presence/absence and synteny information, converged on the same topologies, showing strong support for scleractinian monophyly and a corallimorpharian sister clade. Our broad phylogenomic approach using sequence-based and sequence-independent analyses provides unambiguous evidence for the monophyly of scleractinian corals and the rejection of corallimorpharians as descendants of a complex coral ancestor.

Shedding light on cell compartmentation in the candidate phylum Poribacteria by high resolution visualisation and transcriptional profiling.

Assigning functions to uncultivated environmental microorganisms continues to be a challenging endeavour. Here, we present a new microscopy protocol for fluorescence in situ hybridisation-correlative light and electron microscopy (FISH-CLEM) that enabled, to our knowledge for the first time, the identification of single cells within their complex microenvironment at electron microscopy resolution. Members of the candidate phylum Poribacteria, common and uncultivated symbionts of marine sponges, were used towards this goal. Cellular 3D reconstructions revealed bipolar, spherical granules of low electron density, which likely represent carbon reserves. Poribacterial activity profiles were retrieved from prokaryotic enriched sponge metatranscriptomes using simulation-based optimised mapping. We observed high transcriptional activity for proteins related to bacterial microcompartments (BMC) and we resolved their subcellular localisation by combining FISH-CLEM with immunohistochemistry (IHC) on ultra-thin sponge tissue sections. In terms of functional relevance, we propose that the BMC-A region may be involved in 1,2-propanediol degradation. The FISH-IHC-CLEM approach was proven an effective toolkit to combine -omics approaches with functional studies and it should be widely applicable in environmental microbiology.

Sex Change in Clownfish: Molecular Insights from Transcriptome Analysis.

Sequential hermaphroditism is a unique reproductive strategy among teleosts that is displayed mainly in fish species living in the coral reef environment. The reproductive biology of hermaphrodites has long been intriguing; however, very little is known about the molecular pathways underlying their sex change. Here, we provide the first de novo transcriptome analyses of a hermaphrodite teleost´s undergoing sex change in its natural environment. Our study has examined relative gene expression across multiple groups-rather than just two contrasting conditions- and has allowed us to explore the differential expression patterns throughout the whole process. Our analysis has highlighted the rapid and complex genomic response of the brain associated with sex change, which is subsequently transmitted to the gonads, identifying a large number of candidate genes, some well-known and some novel, involved in the process. The present study provides strong evidence of the importance of the sex steroidogenic machinery during sex change in clownfish, with the aromatase gene playing a central role, both in the brain and the gonad. This work constitutes the first genome-wide study in a social sex-changing species and provides insights into the genetic mechanism governing social sex change and gonadal restructuring in protandrous hermaphrodites.

Characterization of piRNAs across postnatal development in mouse brain.

PIWI-interacting RNAs (piRNAs) are responsible for maintaining the genome stability by silencing retrotransposons in germline tissues- where piRNAs were first discovered and thought to be restricted. Recently, novel functions were reported for piRNAs in germline and somatic cells. Using deep sequencing of small RNAs and CAGE of postnatal development of mouse brain, we identified piRNAs only in adult mouse brain. These piRNAs have similar sequence length as those of MILI-bound piRNAs. In addition, we predicted novel candidate regulators and putative targets of adult brain piRNAs.

Hologenome analysis of two marine sponges with different microbiomes.

BACKGROUND: Sponges (Porifera) harbor distinct microbial consortia within their mesohyl interior. We herein analysed the hologenomes of Stylissa carteri and Xestospongia testudinaria, which notably differ in their microbiome content. RESULTS: Our analysis revealed that S. carteri has an expanded repertoire of immunological domains, specifically Scavenger Receptor Cysteine-Rich (SRCR)-like domains, compared to X. testudinaria. On the microbial side, metatranscriptome analyses revealed an overrepresentation of potential symbiosis-related domains in X. testudinaria. CONCLUSIONS: Our findings provide genomic insights into the molecular mechanisms underlying host-symbiont coevolution and may serve as a roadmap for future hologenome analyses.

miRNA Repertoires of Demosponges Stylissa carteri and Xestospongia testudinaria.

MicroRNAs (miRNAs) are small regulatory RNAs that are involved in many biological process in eukaryotes. They play a crucial role in modulating genetic expression of their targets, which makes them integral components of transcriptional regulatory networks. As sponges (phylum Porifera) are commonly considered the most basal metazoan, the in-depth capture of miRNAs from these organisms provides additional clues to the evolution of miRNA families in metazoans. Here, we identified the core proteins involved in the biogenesis of miRNAs, and obtained evidence for bona fide miRNA sequences for two marine sponges Stylissa carteri and Xestospongia testudinaria (11 and 19 respectively). Our analysis identified several miRNAs that are conserved amongst demosponges, and revealed that all of the novel miRNAs identified in these two species are specific to the class Demospongiae.

Transcriptome and proteome dynamics in larvae of the barnacle Balanus Amphitrite from the Red Sea.

BACKGROUND: The barnacle Balanus amphitrite is widely distributed in marine shallow and tidal waters, and has significant economic and ecological importance. Nauplii, the first larval stage of most crustaceans, are extremely abundant in the marine zooplankton. However, a lack of genome information has hindered elucidation of the molecular mechanisms of development, settlement and survival strategies in extreme marine environments. We sequenced and constructed the genome dataset for nauplii to obtain comprehensive larval genetic information. We also investigated iTRAQ-based protein expression patterns to reveal the molecular basis of nauplii development, and to gain information on larval survival strategies in the Red Sea marine environment. RESULTS: A nauplii larval transcript dataset, containing 92,117 predicted open reading frames (ORFs), was constructed and used as a reference for the proteome analysis. Genes related to translation, oxidative phosphorylation and cytoskeletal development were highly abundant. We observed remarkable plasticity in the proteome of Red Sea larvae. The proteins associated with development, stress responses and osmoregulation showed the most significant differences between the two larval populations studied. The synergistic overexpression of heat shock and osmoregulatory proteins may facilitate larval survival in intertidal habitats or in extreme environments. CONCLUSIONS: We presented, for the first time, comprehensive transcriptome and proteome datasets for Red Sea nauplii. The datasets provide a foundation for future investigations focused on the survival mechanisms of other crustaceans in extreme marine environments.

The co-transcriptome of uropathogenic Escherichia coli-infected mouse macrophages reveals new insights into host-pathogen interactions.

Urinary tract infections (UTI) are among the most common infections in humans. Uropathogenic Escherichia coli (UPEC) can invade and replicate within bladder epithelial cells, and some UPEC strains can also survive within macrophages. To understand the UPEC transcriptional programme associated with intramacrophage survival, we performed host-pathogen co-transcriptome analyses using RNA sequencing. Mouse bone marrow-derived macrophages (BMMs) were challenged over a 24 h time course with two UPEC reference strains that possess contrasting intramacrophage phenotypes: UTI89, which survives in BMMs, and 83972, which is killed by BMMs. Neither of these strains caused significant BMM cell death at the low multiplicity of infection that was used in this study. We developed an effective computational framework that simultaneously separated, annotated and quantified the mammalian and bacterial transcriptomes. Bone marrow-derived macrophages responded to the two UPEC strains with a broadly similar gene expression programme. In contrast, the transcriptional responses of the UPEC strains diverged markedly from each other. We identified UTI89 genes up-regulated at 24 h post-infection, and hypothesized that some may contribute to intramacrophage survival. Indeed, we showed that deletion of one such gene (pspA) significantly reduced UTI89 survival within BMMs. Our study provides a technological framework for simultaneously capturing global changes at the transcriptional level in co-cultures, and has generated new insights into the mechanisms that UPEC use to persist within the intramacrophage environment.

Transcriptome analysis elucidates key developmental components of bryozoan lophophore development.

The most recent phylogenomic study suggested that Bryozoa (Ectoprocta), Brachiopoda, and Phoronida are monophyletic, implying that the lophophore of bryozoans, phoronids and brachiopods is a synapomorphy. Understanding the molecular mechanisms of the lophophore development of the Lophophorata clade can therefore provide us a new insight into the formation of the diverse morphological traits in metazoans. In the present study, we profiled the transcriptome of the Bryozoan (Ectoproct) Bugula neritina during the swimming larval stage (SW) and the early (4 h) and late (24 h) metamorphic stages using the Illumina HiSeq2000 platform. Various genes that function in development, the immune response and neurogenesis showed differential expression levels during metamorphosis. In situ hybridization of 23 genes that participate in the Wnt, BMP, Notch, and Hedgehog signaling pathways revealed their regulatory roles in the development of the lophophore and the ancestrula digestive tract. Our findings support the hypothesis that developmental precursors of the lophophore and the ancestrula digestive tract are pre-patterned by the differential expression of key developmental genes according to their fate. This study provides a foundation to better understand the developmental divergence and/or convergence among developmental precursors of the lophophore of bryozoans, branchiopods and phoronids.

Dynamic epigenetic control of highly conserved noncoding elements.

BACKGROUND: Many noncoding genomic loci have remained constant over long evolutionary periods, suggesting that they are exposed to strong selective pressures. The molecular functions of these elements have been partially elucidated, but the fundamental reason for their extreme conservation is still unknown. RESULTS: To gain new insights into the extreme selection of highly conserved noncoding elements (HCNEs), we used a systematic analysis of multi-omic data to study the epigenetic regulation of such elements during the development of Drosophila melanogaster. At the sequence level, HCNEs are GC-rich and have a characteristic oligomeric composition. They have higher levels of stable nucleosome occupancy than their flanking regions, and lower levels of mononucleosomes and H3.3, suggesting that these regions reside in compact chromatin. Furthermore, these regions showed remarkable modulations in histone modification and the expression levels of adjacent genes during development. Although HCNEs are primarily initiated late in replication, about 10% were related to early replication origins. Finally, HCNEs showed strong enrichment within lamina-associated domains. CONCLUSION: HCNEs have distinct and protective sequence properties, undergo dynamic epigenetic regulation, and appear to be associated with the structural components of the chromatin, replication origins, and nuclear matrix. These observations indicate that such elements are likely to have essential cellular functions, and offer insights into their epigenetic properties.

Revealing microbial functional activities in the Red Sea sponge Stylissa carteri by metatranscriptomics.

Sponges are important components of marine benthic environments and are associated with microbial symbionts that carry out ecologically relevant functions. Stylissa carteri is an abundant, low-microbial abundance species in the Red Sea. We aimed to achieve the functional and taxonomic characterization of the most actively expressed prokaryotic genes in S. carteri. Prokaryotic mRNA was enriched from sponge total RNA, sequenced using Illumina HiSeq technology and annotated using the metagenomics Rapid Annotation using Subsystem Technology (MG-RAST) pipeline. We detected high expression of archaeal ammonia oxidation and photosynthetic carbon fixation by members of the genus Synechococcus. Functions related to stress response and membrane transporters were among the most highly expressed by S. carteri symbionts. Unexpectedly, gene functions related to methylotrophy were highly expressed by gammaproteobacterial symbionts. The presence of seawater-derived microbes is indicated by the phylogenetic proximity of organic carbon transporters to orthologues of members from the SAR11 clade. In summary, we revealed the most expressed functions of the S. carteri-associated microbial community and linked them to the dominant taxonomic members of the microbiome. This work demonstrates the applicability of metatranscriptomics to explore poorly characterized symbiotic consortia and expands our knowledge of the ecologically relevant functions carried out by coral reef sponge symbionts.