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.

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.

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 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.

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.

Specificity and transcriptional activity of microbiota associated with low and high microbial abundance sponges from the Red Sea.

Marine sponges are generally classified as high microbial abundance (HMA) and low microbial abundance (LMA) species. Here, 16S rRNA amplicon sequencing was applied to investigate the diversity, specificity and transcriptional activity of microbes associated with an LMA sponge (Stylissa carteri), an HMA sponge (Xestospongia testudinaria) and sea water collected from the central Saudi Arabia coast of the Red Sea. Altogether, 887 068 denoised sequences were obtained, of which 806 661 sequences remained after quality control. This resulted in 1477 operational taxonomic units (OTUs) that were assigned to 27 microbial phyla. The microbial composition of S. carteri was more similar to that of sea water than to that of X. testudinaria, which is consistent with the observation that the sequence data set of S. carteri contained many more possibly sea water sequences (~24%) than the X. testudinaria data set (~6%). The most abundant OTUs were shared between all three sources (S. carteri, X. testudinaria, sea water), while rare OTUs were unique to any given source. Despite this high degree of overlap, each sponge species contained its own specific microbiota. The X. testudinaria-specific bacterial taxa were similar to those already described for this species. A set of S. carteri-specific bacterial taxa related to Proteobacteria and Nitrospira was identified, which are likely permanently associated with S. carteri. The transcriptional activity of sponge-associated microorganisms correlated well with their abundance. Quantitative PCR revealed the presence of Poribacteria, representing typical sponge symbionts, in both sponge species and in sea water; however, low transcriptional activity in sea water suggested that Poribacteria are not active outside the host context.

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.

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.

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.

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.

The evolution of ultraconserved elements with different phylogenetic origins.

BACKGROUND: Ultraconserved elements of DNA have been identified in vertebrate and invertebrate genomes. These elements have been found to have diverse functions, including enhancer activities in developmental processes. The evolutionary origins and functional roles of these elements in cellular systems, however, have not yet been determined. RESULTS: Here, we identified a wide range of ultraconserved elements common to distant species, from primitive aquatic organisms to terrestrial species with complicated body systems, including some novel elements conserved in fruit fly and human. In addition to a well-known association with developmental genes, these DNA elements have a strong association with genes implicated in essential cell functions, such as epigenetic regulation, apoptosis, detoxification, innate immunity, and sensory reception. Interestingly, we observed that ultraconserved elements clustered by sequence similarity. Furthermore, species composition and flanking genes of clusters showed lineage-specific patterns. Ultraconserved elements are highly enriched with binding sites to developmental transcription factors regardless of how they cluster. CONCLUSION: We identified large numbers of ultraconserved elements across distant species. Specific classes of these conserved elements seem to have been generated before the divergence of taxa and fixed during the process of evolution. Our findings indicate that these ultraconserved elements are not the exclusive property of higher modern eukaryotes, but rather transmitted from their metazoan ancestors.

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.

A genetic screen for modifiers of Drosophila caspase Dcp-1 reveals caspase involvement in autophagy and novel caspase-related genes.

BACKGROUND: Caspases are cysteine proteases with essential functions in the apoptotic pathway; their proteolytic activity toward various substrates is associated with the morphological changes of cells. Recent reports have described non-apoptotic functions of caspases, including autophagy. In this report, we searched for novel modifiers of the phenotype of Dcp-1 gain-of-function (GF) animals by screening promoter element- inserted Drosophila melanogaster lines (EP lines). RESULTS: We screened approximately 15,000 EP lines and identified 72 Dcp-1-interacting genes that were classified into 10 groups based on their functions and pathways: 4 apoptosis signaling genes, 10 autophagy genes, 5 insulin/IGF and TOR signaling pathway genes, 6 MAP kinase and JNK signaling pathway genes, 4 ecdysone signaling genes, 6 ubiquitination genes, 11 various developmental signaling genes, 12 transcription factors, 3 translation factors, and 11 other unclassified genes including 5 functionally undefined genes. Among them, insulin/IGF and TOR signaling pathway, MAP kinase and JNK signaling pathway, and ecdysone signaling are known to be involved in autophagy. Together with the identification of autophagy genes, the results of our screen suggest that autophagy counteracts Dcp-1-induced apoptosis. Consistent with this idea, we show that expression of eGFP-Atg5 rescued the eye phenotype caused by Dcp-1 GF. Paradoxically, we found that over-expression of full-length Dcp-1 induced autophagy, as Atg8b-GFP, an indicator of autophagy, was increased in the eye imaginal discs and in the S2 cell line. Taken together, these data suggest that autophagy suppresses Dcp-1-mediated apoptotic cell death, whereas Dcp-1 positively regulates autophagy, possibly through feedback regulation. CONCLUSIONS: We identified a number of Dcp-1 modifiers that genetically interact with Dcp-1-induced cell death. Our results showing that Dcp-1 and autophagy-related genes influence each other will aid future investigations of the complicated relationships between apoptosis and autophagy.

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.

bZIPDB: a database of regulatory information for human bZIP transcription factors.

BACKGROUND: Basic region-leucine zipper (bZIP) proteins are a class of transcription factors (TFs) that play diverse roles in eukaryotes. Malfunctions in these proteins lead to cancer and various other diseases. For detailed characterization of these TFs, further public resources are required. DESCRIPTION: We constructed a database, designated bZIPDB, containing information on 49 human bZIP TFs, by means of automated literature collection and manual curation. bZIPDB aims to provide public data required for deciphering the gene regulatory network of the human bZIP family, e.g., evaluation or reference information for the identification of regulatory modules. The resources provided by bZIPDB include (1) protein interaction data including direct binding, phosphorylation and functional associations between bZIP TFs and other cellular proteins, along with other types of interactions, (2) bZIP TF-target gene relationships, (3) the cellular network of bZIP TFs in particular cell lines, and (4) gene information and ontology. In the current version of the database, 721 protein interactions and 560 TF-target gene relationships are recorded. bZIPDB is annually updated for the newly discovered information. CONCLUSION: bZIPDB is a repository of detailed regulatory information for human bZIP TFs that is collected and processed from the literature, designed to facilitate analysis of this protein family. bZIPDB is available for public use at http://biosoft.kaist.ac.kr/bzipdb.

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.

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.

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.

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.