Comparative transcriptome analysis reveals insights into the streamlined genomes of haplosclerid demosponges

A novel target-enriched multilocus assay for sponges (Porifera): Red Sea Haplosclerida (Demospongiae) as a test case

With declining biodiversity worldwide, a better understanding of species diversity and their relationships is imperative for conservation and management efforts. Marine sponges are species-rich ecological key players on coral reefs, but their species diversity is still poorly understood. This is particularly true for the demosponge order Haplosclerida, whose systematic relationships are contentious due to the incongruencies between morphological and molecular phylogenetic hypotheses. The single gene markers applied in previous studies did not resolve these discrepancies. Hence, there is a high need for a genome-wide approach to derive a phylogenetically robust classification and understand this group's evolutionary relationships. To this end, we developed a target enrichment-based multilocus probe assay for the order Haplosclerida using transcriptomic data. This probe assay consists of 20,000 enrichment probes targeting 2956 ultraconserved elements in coding (i.e. exon) regions across the genome and was tested on 26 haplosclerid specimens from the Red Sea. Our target-enrichment approach correctly placed our samples in a well-supported phylogeny, in agreement with previous haplosclerid molecular phylogenies. Our results demonstrate the applicability of high-resolution genomic methods in a systematically complex marine invertebrate group and provide a promising approach for robust phylogenies of Haplosclerida. Subsequently, this will lead to biologically unambiguous taxonomic revisions, better interpretations of biological and ecological observations and new avenues for applied research, conservation and managing declining marine diversity.

Phylomitogenomics bolsters the high-level classification of Demospongiae (phylum Porifera)

Class Demospongiae is the largest in the phylum Porifera (Sponges) and encompasses nearly 8,000 accepted species in three subclasses: Keratosa, Verongimorpha, and Heteroscleromorpha. Subclass Heteroscleromorpha contains ∼90% of demosponge species and is subdivided into 17 orders. The higher level classification of demosponges underwent major revision as the result of nearly three decades of molecular studies. However, because most of the previous molecular work only utilized partial data from a small number of nuclear and mitochondrial (mt) genes, this classification scheme needs to be tested by larger datasets. Here we compiled a mt dataset for 136 demosponge species-including 64 complete or nearly complete and six partial mt-genome sequences determined or assembled for this study-and used it to test phylogenetic relationships among Demospongiae in general and Heteroscleromorpha in particular. We also investigated the phylogenetic position of Myceliospongia araneosa, a highly unusual demosponge without spicules and spongin fibers, currently classified as Demospongiae incertae sedis, for which molecular data were not available. Our results support the previously inferred sister-group relationship between Heteroscleromorpha and Keratosa + Verongimorpha and suggest five main clades within Heteroscleromorpha: Clade C0 composed of order Haplosclerida; Clade C1 composed of Scopalinida, Sphaerocladina, and Spongillida; Clade C2 composed of Axinellida, Biemnida, Bubarida; Clade C3 composed of Tetractinellida; and Clade C4 composed of Agelasida, Clionaida, Desmacellida, Merliida, Suberitida, Poecilosclerida, Polymastiida, and Tethyida. The inferred relationships among these clades were (C0(C1(C2(C3+C4)))). Analysis of molecular data from M. araneosa placed it in the C3 clade as a sister taxon to the highly skeletonized tetractinellids Microscleroderma sp. and Leiodermatium sp. Molecular clock analysis dated divergences among the major clades in Heteroscleromorpha from the Cambrian to the Early Silurian, the origins of most heteroscleromorph orders in the middle Paleozoic, and the most basal splits within these orders around the Paleozoic to Mesozoic transition. Overall, the results of this study are mostly congruent with the accepted classification of Heteroscleromorpha, but add temporal perspective and new resolution to phylogenetic relationships within this subclass.

An ependymin-related blue carotenoprotein decorates marine blue sponge

Marine animals display diverse vibrant colors, but the mechanisms underlying their specific coloration remain to be clarified. Blue coloration is known to be achieved through a bathochromic shift of the orange carotenoid astaxanthin (AXT) by the crustacean protein crustacyanin, but other examples have not yet been well investigated. Here, we identified an ependymin (EPD)-related water-soluble blue carotenoprotein responsible for the specific coloration of the marine blue sponge Haliclona sp. EPD was originally identified in the fish brain as a protein involved in memory consolidation and neuronal regeneration. The purified blue protein, designated as EPD-related blue carotenoprotein-1, was identified as a secreted glycoprotein. We show that it consists of a heterodimer that binds orange AXT and mytiloxanthin and exhibits a bathochromic shift. Our crystal structure analysis of the natively purified EPD-related blue carotenoprotein-1 revealed that these two carotenoids are specifically bound to the heterodimer interface, where the polyene chains are aligned in parallel to each other like in β-crustacyanin, although the two proteins are evolutionary and structurally unrelated. Furthermore, using reconstitution assays, we found that incomplete bathochromic shifts occurred when the protein bound to only AXT or mytiloxanthin. Taken together, we identified an EPD in a basal metazoan as a blue protein that decorates the sponge body by binding specific structurally unrelated carotenoids.

Divergent Ca2+/calmodulin feedback regulation of CaV1 and CaV2 voltage-gated calcium channels evolved in the common ancestor of Placozoa and Bilateria

Ca_V1 and Ca_V2 voltage-gated calcium channels evolved from an ancestral Ca_V1/2 channel via gene duplication somewhere near the stem animal lineage. The divergence of these channel types led to distinguishing functional properties that are conserved among vertebrates and bilaterian invertebrates and contribute to their unique cellular roles. One key difference pertains to their regulation by calmodulin (CaM), wherein bilaterian Ca_V1 channels are uniquely subject to pronounced, buffer-resistant Ca²⁺/CaM-dependent inactivation, permitting negative feedback regulation of calcium influx in response to local cytoplasmic Ca²⁺ rises. Early diverging, nonbilaterian invertebrates also possess Ca_V1 and Ca_V2 channels, but it is unclear whether they share these conserved functional features. The most divergent animals to possess both Ca_V1 and Ca_V2 channels are placozoans such as Trichoplax adhaerens, which separated from other animals over 600 million years ago shortly after their emergence. Hence, placozoans can provide important insights into the early evolution of Ca_V1 and Ca_V2 channels. Here, we build upon previous characterization of Trichoplax Ca_V channels by determining the cellular expression and ion-conducting properties of the Ca_V1 channel orthologue, TCa_V1. We show that TCa_V1 is expressed in neuroendocrine-like gland cells and contractile dorsal epithelial cells. In vitro, this channel conducts dihydropyridine-insensitive, high-voltage–activated Ca²⁺ currents with kinetics resembling those of rat Ca_V1.2 but with left-shifted voltage sensitivity for activation and inactivation. Interestingly, TCa_V1, but not TCa_V2, exhibits buffer-resistant Ca²⁺/CaM-dependent inactivation, indicating that this functional divergence evolved prior to the emergence of bilaterian animals and may have contributed to their unique adaptation for cytoplasmic Ca²⁺ signaling within various cellular contexts.

Actinoporin-like Proteins Are Widely Distributed in the Phylum Porifera

Actinoporins are proteinaceous toxins known for their ability to bind to and create pores in cellular membranes. This quality has generated interest in their potential use as new tools, such as therapeutic immunotoxins. Isolated historically from sea anemones, genes encoding for similar actinoporin-like proteins have since been found in a small number of other animal phyla. Sequencing and de novo assembly of Irish Haliclona transcriptomes indicated that sponges also possess similar genes. An exhaustive analysis of publicly available sequencing data from other sponges showed that this is a potentially widespread feature of the Porifera. While many sponge proteins possess a sequence similarity of 27.70–59.06% to actinoporins, they show consistency in predicted structure. One gene copy from H. indistincta has significant sequence similarity to sea anemone actinoporins and possesses conserved residues associated with the fundamental roles of sphingomyelin recognition, membrane attachment, oligomerization, and pore formation, indicating that it may be an actinoporin. Phylogenetic analyses indicate frequent gene duplication, no distinct clade for sponge-derived proteins, and a stronger signal towards actinoporins than similar proteins from other phyla. Overall, this study provides evidence that a diverse array of Porifera represents a novel source of actinoporin-like proteins which may have biotechnological and pharmaceutical applications.

Microbiome diversity and host immune functions influence survivorship of sponge holobionts under future ocean conditions

The sponge-associated microbial community contributes to the overall health and adaptive capacity of the sponge holobiont. This community is regulated by the environment and the immune system of the host. However, little is known about the effect of environmental stress on the regulation of host immune functions and how this may, in turn, affect sponge-microbe interactions. In this study, we compared the bacterial diversity and immune repertoire of the demosponge, Neopetrosia compacta, and the calcareous sponge, Leucetta chagosensis, under varying levels of acidification and warming stress based on climate scenarios predicted for 2100. Neopetrosia compacta harbors a diverse microbial community and possesses a rich repertoire of scavenger receptors while L. chagosensis has a less diverse microbiome and an expanded range of pattern recognition receptors and immune response-related genes. Upon exposure to RCP 8.5 conditions, the microbiome composition and host transcriptome of N. compacta remained stable, which correlated with high survival (75%). In contrast, tissue necrosis and low survival (25%) of L. chagosensis was accompanied by microbial community shifts and downregulation of host immune-related pathways. Meta-analysis of microbiome diversity and immunological repertoire across poriferan classes further highlights the importance of host-microbe interactions in predicting the fate of sponges under future ocean conditions.

Mitochondrial evolution in the Demospongiae (Porifera): Phylogeny, divergence time, and genome biology

The sponge class Demospongiae is the most speciose and morphologically diverse in the phylum Porifera, and the species within it are vital components of a range of ecosystems worldwide. Despite their ubiquity, a number of recalcitrant problems still remain to be solved regarding their phylogenetic inter-relationships, the timing of their appearance, and their mitochondrial biology, the latter of which is only beginning to be investigated. Here we generated 14 new demosponge mitochondrial genomes which, alongside previously published mitochondrial resources, were used to address these issues. In addition to phylogenomic analysis, we have used syntenic data and analysis of coding regions to forge a framework for understanding the inter-relationships between Demospongiae sub-classes and orders. We have also leveraged our new resources to study the mitochondrial biology of these clades in terms of codon usage, optimisation and gene expression, to understand how these vital cellular components may have contributed to the success of the Porifera. Our results strongly support a sister relationship between Keratosa and (Verongimorpha + Heteroscleromorpha), contradicting previous studies using nuclear markers. Our study includes one species of Clionaida, and show for the first time support for a grouping of Suberitida+(Clionaida+(Tethyida + Poecilosclerida). The findings of our phylogenetic analyses are supported by in-depth examination of structural and coding-level evidence from our mitochondrial data. A time-calibrated phylogeny estimated the origin of Demospongiae in the Cambrian (~529 Mya), and suggests that most demosponge order crown-groups emerged in the Mesozoic. This work therefore provides a robust basis for considering demosponge phylogenetic relationships, as well as essential mitochondrial data for understanding the biological basis for their success and diversity.

Receptors Mediating Host-Microbiota Communication in the Metaorganism: The Invertebrate Perspective

Multicellular organisms live in close association with a plethora of microorganism, which have a profound effect on multiple host functions. As such, the microbiota and its host form an intimate functional entity, termed the metaorganism or holobiont. But how does the metaorganism communicate? Which receptors recognize microbial signals, mediate the effect of the microbiota on host physiology or regulate microbiota composition and homeostasis? In this review we provide an overview on the function of different receptor classes in animal host-microbiota communication. We put a special focus on invertebrate hosts, including both traditional invertebrate models such as Drosophila melanogaster and Caenorhabditis elegans and “non-model” invertebrates in microbiota research. Finally, we highlight the potential of invertebrate systems in studying mechanism of host-microbiota interactions.

Sponge Behavior and the Chemical Basis of Responses: A Post-Genomic View

Sponges perceive and respond to a range of stimuli. How they do this is still difficult to pin down despite now having transcriptomes and genomes of an array of species. Here we evaluate the current understanding of sponge behavior and present new observations on sponge activity in situ. We also explore biosynthesis pathways available to sponges from data in genomes/transcriptomes of sponges and other non-bilaterians with a focus on exploring the role of chemical signaling pathways mediating sponge behavior and how such chemical signal pathways may have evolved. Sponge larvae respond to light but opsins are not used, nor is there a common photoreceptor molecule or mechanism used across sponge groups. Other cues are gravity and chemicals. In situ recordings of behavior show that both shallow and deep-water sponges move a lot over minutes and hours, and correlation of behavior with temperature, pressure, oxygen, and water movement suggests that at least one sponge responds to changes in atmospheric pressure. The sensors for these cues as far as we know are individual cells and, except in the case of electrical signaling in Hexactinellida, these most likely act as independent effectors, generating a whole-body reaction by the global reach of the stimulus to all parts of the animal. We found no evidence for use of conventional neurotransmitters such as serotonin and dopamine. Intriguingly, some chemicals synthesized by symbiont microbes could mean other more complex signaling occurs, but how that interplay might happen is not understood. Our review suggests chemical signaling pathways found in sponges do not reflect loss of a more complex set.

Evolution of the main skeleton-forming genes in sponges (phylum Porifera) with special focus on the marine Haplosclerida (class Demospongiae)

The skeletons of sponges (Phylum Porifera) are comprised of collagen, often embedded with small siliceous structures (spicules) arranged in various forms to provide strength and flexibility. The main proteins responsible for the formation of the spicules in demosponges are the silicateins, which are related to the cathepsins L of other animals. While the silicatein active site, necessary for the formation of biosilica crystals, is characterized by the amino acids SHN, different variants of the silicatein genes have been found, some that retain SHN at the active site and some that don't. As part of an effort to further understand skeleton formation in marine sponges of the order Haplosclerida, a search for all silicatein variants were made in Irish species representing the main clades of this large sponge group. For this task, transcriptomes were sequenced and de novo assembled from Haliclona oculata, H. simulans and H. indistincta. Silicatein genes were identified from these and all available genomes and transcriptomes from Porifera. These were analysed along with all complete silicateins from GenBank. Silicateins were only found in species belonging to the class Demospongiae but excluding Keratosa and Verongimorpha and there was significant duplication and diversity of these genes. Silicateins showing SHN at the active site were polyphyletic. Indeed silicatein sequences were divided into six major clades (CHNI, CHNII, CHNIII, SHNI, SHNII and C/SQN). In those clades where haplosclerids were well represented the silicatein phylogeny reflected previous ribosomal and mitochondrial topologies. The most basal silicatein clade (CHNI) contained sequences only from marine haplosclerids and freshwater sponges while one silicatein from H. indistincta was more related to cathepsins L (outgroup) than to the overall silicatein clade indicating the presence of an old silicatein or an intermediary form. This data could suggest that marine haplosclerids were one of the first groups of extant demosponges to acquire silicatein genes. Furthermore, we suggest that the paucity of spicule types in this group may be due to their single copy of SHNI variants, and the lack of a silintaphin gene.

Characterization of mitochondrial proteomes of nonbilaterian animals

Mitochondria require ~1,500 proteins for their maintenance and proper functionality, which constitute the mitochondrial proteome (mt-proteome). Although a few of these proteins, mostly subunits of the electron transport chain complexes, are encoded in mitochondrial DNA (mtDNA), the vast majority are encoded in the nuclear genome and imported to the organelle. Previous studies have shown a continuous and complex evolution of mt-proteome among eukaryotes. However, there was less attention paid to mt-proteome evolution within Metazoa, presumably because animal mtDNA and, by extension, animal mitochondria are often considered to be uniform. In this analysis, two bioinformatic approaches (Orthologue-detection and Mitochondrial Targeting Sequence prediction) were used to identify mt-proteins in 23 species from four nonbilaterian phyla: Cnidaria, Ctenophora, Placozoa, and Porifera, as well as two choanoflagellates, the closest animal relatives. Our results revealed a large variation in mt-proteome in nonbilaterian animals in size and composition. Myxozoans, highly reduced cnidarian parasites, possessed the smallest inferred mitochondrial proteomes, while calcareous sponges possessed the largest. About 513 mitochondrial orthologous groups were present in all nonbilaterian phyla and human. Interestingly, 42 human mitochondrial proteins were not identified in any nonbilaterian species studied and represent putative innovations along the bilaterian branch. Several of these proteins were involved in apoptosis and innate immunity, two processes known to evolve within Metazoa. Conversely, several proteins identified as mitochondrial in nonbilaterian phyla and animal outgroups were absent in human, representing cases of possible loss. Finally, a few human cytosolic proteins, such as histones and cytosolic ribosomal proteins, were predicted to be targeted to mitochondria in nonbilaterian animals. Overall, our analysis provides the first step in characterization of mt-proteomes in nonbilaterian animals and understanding evolution of animal mt-proteome. © 2018 IUBMB Life, 70(12):1289-1301, 2018.

Differential expression of immune receptors in two marine sponges upon exposure to microbial-associated molecular patterns

The innate immune system helps animals to navigate the microbial world. The response to microbes relies on the specific recognition of microbial-associated molecular patterns (MAMPs) by immune receptors. Sponges (phylum Porifera), as early-diverging animals, provide insights into conserved mechanisms for animal-microbe crosstalk. However, experimental data is limited. We adopted an experimental approach followed by RNA-Seq and differential gene expression analysis in order to characterise the sponge immune response. Two Mediterranean species, Aplysina aerophoba and Dysidea avara, were exposed to a “cocktail” of MAMPs (lipopolysaccharide and peptidoglycan) or to sterile artificial seawater (control) and sampled 1 h, 3 h, and 5 h post-treatment for RNA-Seq. The response involved, first and foremost, a higher number of differentially-expressed genes in A. aerophoba than D. avara. Secondly, while both species constitutively express a diverse repertoire of immune receptors, they differed in their expression profiles upon MAMP challenge. The response in D. avara was mediated by increased expression of two NLR genes, whereas the response in A. aerophoba involved SRCR and GPCR genes. From the set of annotated genes we infer that both species activated apoptosis in response to MAMPs while in A. aerophoba phagocytosis was additionally stimulated. Our study assessed for the first time the transcriptomic responses of sponges to MAMPs and revealed conserved and species-specific features of poriferan immunity as well as genes potentially relevant to animal-microbe interactions.

A de novo approach to disentangle partner identity and function in holobiont systems

BACKGROUND

Study of meta-transcriptomic datasets involving non-model organisms represents bioinformatic challenges. The production of chimeric sequences and our inability to distinguish the taxonomic origins of the sequences produced are inherent and recurrent difficulties in de novo assembly analyses. As the study of holobiont meta-transcriptomes is affected by challenges invoked above, we propose an innovative bioinformatic approach to tackle such difficulties and tested it on marine models as a proof of concept.

RESULTS

We considered three holobiont models, of which two transcriptomes were previously published and a yet unpublished transcriptome, to analyze and sort their raw reads using Short Read Connector, a k-mer based similarity method. Before assembly, we thus defined four distinct categories for each holobiont meta-transcriptome: host reads, symbiont reads, shared reads, and unassigned reads. Afterwards, we observed that independent de novo assemblies for each category led to a diminution of the number of chimeras compared to classical assembly methods. Moreover, the separation of each partner's transcriptome offered the independent and comparative exploration of their functional diversity in the holobiont. Finally, our strategy allowed to propose new functional annotations for two well-studied holobionts (a Cnidaria-Dinophyta, a Porifera-Bacteria) and a first meta-transcriptome from a planktonic Radiolaria-Dinophyta system forming widespread symbiotic association for which our knowledge is considerably limited.

CONCLUSIONS

In contrast to classical assembly approaches, our bioinformatic strategy generates less de novo assembled chimera and allows biologists to study separately host and symbiont data from a holobiont mixture. The pre-assembly separation of reads using an efficient tool as Short Read Connector is an effective way to tackle meta-transcriptomic challenges and offers bright perpectives to study holobiont systems composed of either well-studied or poorly characterized symbiotic lineages and ultimately expand our knowledge about these associations.

Think like a sponge: The genetic signal of sensory cells in sponges

A complex genetic repertoire underlies the apparently simple body plan of sponges. Among the genes present in poriferans are those fundamental to the sensory and nervous systems of other animals. Sponges are dynamic and sensitive animals and it is intuitive to link these genes to behaviour. The proposal that ctenophores are the earliest diverging metazoan has led to the question of whether sponges possess a 'pre-nervous' system or have undergone nervous system loss. Both lines of thought generally assume that the last common ancestor of sponges and eumetazoans possessed the genetic modules that underlie sensory abilities. By corollary extant sponges may possess a sensory cell homologous to one present in the last common ancestor, a hypothesis that has been studied by gene expression. We have performed a meta-analysis of all gene expression studies published to date to explore whether gene expression is indicative of a feature's sensory function. In sponges we find that eumetazoan sensory-neural markers are not particularly expressed in structures with known sensory functions. Instead it is common for these genes to be expressed in cells with no known or uncharacterized sensory function. Indeed, many sensory-neural markers so far studied are expressed during development, perhaps because many are transcription factors. This suggests that the genetic signal of a sponge sensory cell is dissimilar enough to be unrecognizable when compared to a bilaterian sensory or neural cell. It is possible that sensory-neural markers have as yet unknown functions in sponge cells, such as assembling an immunological synapse in the larval globular cell. Furthermore, the expression of sensory-neural markers in non-sensory cells, such as adult and larval epithelial cells, suggest that these cells may have uncharacterized sensory functions. While this does not rule out the co-option of ancestral sensory modules in later evolving groups, a distinct genetic foundation may underlie the sponge sensory system.

Physiology and Evolution of Voltage-Gated Calcium Channels in Early Diverging Animal Phyla: Cnidaria, Placozoa, Porifera and Ctenophora

Voltage-gated calcium (Ca_v) channels serve dual roles in the cell, where they can both depolarize the membrane potential for electrical excitability, and activate transient cytoplasmic Ca²⁺ signals. In animals, Ca_v channels play crucial roles including driving muscle contraction (excitation-contraction coupling), gene expression (excitation-transcription coupling), pre-synaptic and neuroendocrine exocytosis (excitation-secretion coupling), regulation of flagellar/ciliary beating, and regulation of cellular excitability, either directly or through modulation of other Ca²⁺-sensitive ion channels. In recent years, genome sequencing has provided significant insights into the molecular evolution of Ca_v channels. Furthermore, expanded gene datasets have permitted improved inference of the species phylogeny at the base of Metazoa, providing clearer insights into the evolution of complex animal traits which involve Ca_v channels, including the nervous system. For the various types of metazoan Ca_v channels, key properties that determine their cellular contribution include: Ion selectivity, pore gating, and, importantly, cytoplasmic protein-protein interactions that direct sub-cellular localization and functional complexing. It is unclear when these defining features, many of which are essential for nervous system function, evolved. In this review, we highlight some experimental observations that implicate Ca_v channels in the physiology and behavior of the most early-diverging animals from the phyla Cnidaria, Placozoa, Porifera, and Ctenophora. Given our limited understanding of the molecular biology of Ca_v channels in these basal animal lineages, we infer insights from better-studied vertebrate and invertebrate animals. We also highlight some apparently conserved cellular functions of Ca_v channels, which might have emerged very early on during metazoan evolution, or perhaps predated it.

Gene Expression Dynamics Accompanying the Sponge Thermal Stress Response

Marine sponges are important members of coral reef ecosystems. Thus, their responses to changes in ocean chemistry and environmental conditions, particularly to higher seawater temperatures, will have potential impacts on the future of these reefs. To better understand the sponge thermal stress response, we investigated gene expression dynamics in the shallow water sponge, Haliclona tubifera (order Haplosclerida, class Demospongiae), subjected to elevated temperature. Using high-throughput transcriptome sequencing, we show that these conditions result in the activation of various processes that interact to maintain cellular homeostasis. Short-term thermal stress resulted in the induction of heat shock proteins, antioxidants, and genes involved in signal transduction and innate immunity pathways. Prolonged exposure to thermal stress affected the expression of genes involved in cellular damage repair, apoptosis, signaling and transcription. Interestingly, exposure to sublethal temperatures may improve the ability of the sponge to mitigate cellular damage under more extreme stress conditions. These insights into the potential mechanisms of adaptation and resilience of sponges contribute to a better understanding of sponge conservation status and the prediction of ecosystem trajectories under future climate conditions.

Detection of Prokaryotic Genes in the Amphimedon queenslandica Genome

Horizontal gene transfer (HGT) is common between prokaryotes and phagotrophic eukaryotes. In metazoans, the scale and significance of HGT remains largely unexplored but is usually linked to a close association with parasites and endosymbionts. Marine sponges (Porifera), which host many microorganisms in their tissues and lack an isolated germ line, are potential carriers of genes transferred from prokaryotes. In this study, we identified a number of potential horizontally transferred genes within the genome of the sponge, Amphimedon queenslandica. We further identified homologs of some of these genes in other sponges. The transferred genes, most of which possess catalytic activity for carbohydrate or protein metabolism, have assimilated host genome characteristics and are actively expressed. The diversity of functions contributed by the horizontally transferred genes is likely an important factor in the adaptation and evolution of A. queenslandica. These findings highlight the potential importance of HGT on the success of sponges in diverse ecological niches.