Innate immunity in the simplest animals – placozoans
BMC Genomics volume 20, Article number: 5 (2019)
Innate immunity provides the core recognition system in animals for preventing infection, but also plays an important role in managing the relationship between an animal host and its symbiont. Most of our knowledge about innate immunity stems from a few animal model systems, but substantial variation between metazoan phyla has been revealed by comparative genomic studies. The exploration of more taxa is still needed to better understand the evolution of immunity related mechanisms. Placozoans are morphologically the simplest organized metazoans and the association between these enigmatic animals and their rickettsial endosymbionts has recently been elucidated. Our analyses of the novel placozoan nuclear genome of Trichoplax sp. H2 and its associated rickettsial endosymbiont genome clearly pointed to a mutualistic and co-evolutionary relationship. This discovery raises the question of how the placozoan holobiont manages symbiosis and, conversely, how it defends against harmful microorganisms. In this study, we examined the annotated genome of Trichoplax sp. H2 for the presence of genes involved in innate immune recognition and downstream signaling.
A rich repertoire of genes belonging to the Toll-like and NOD-like receptor pathways, to scavenger receptors and to secreted fibrinogen-related domain genes was identified in the genome of Trichoplax sp. H2. Nevertheless, the innate immunity related pathways in placozoans deviate in several instances from well investigated vertebrates and invertebrates. While true Toll- and NOD-like receptors are absent, the presence of many genes of the downstream signaling cascade suggests at least primordial Toll-like receptor signaling in Placozoa. An abundance of scavenger receptors, fibrinogen-related domain genes and Apaf-1 genes clearly constitutes an expansion of the immunity related gene repertoire specific to Placozoa.
The found wealth of immunity related genes present in Placozoa is surprising and quite striking in light of the extremely simple placozoan body plan and their sparse cell type makeup. Research is warranted to reveal how Placozoa utilize this immune repertoire to manage and maintain their associated microbiota as well as to fend-off pathogens.
The discrimination between self and non-self is crucial for defense against infection and thus for maintaining an organism’s genetic and somatic integrity. In animals this defense is accomplished by an innate immune system comprised of pattern recognition receptors (PRRs) specifically recognizing pathogen-associated molecular patterns (PAMPs) . While vertebrates have evolved a highly sophisticated adaptive immune system, the evolutionary older innate immune system constitutes the main defense against novel invading pathogens in vertebrates and invertebrates and has also been described in lower metazoans like Cnidaria (e.g. ). In recent years, the exploration of genomic data from other non-model animals, such as lophotrochozoans or basal chordates, has further extended our understanding about the diversity and evolution of innate immunity components present in animals (e.g. [3,4,5,6]). In addition, it has been recognized that innate immunity plays a fundamental role in managing symbiotic interactions between microorganisms and animal hosts. Accordingly, it has been hypothesized that it has, at least in part, evolved to mediate a bidirectional cross-talk between both partners of a mutualistic relationship . In this view a host’s innate immunity is able to discern, tolerate and promote beneficial microorganisms while simultaneously fending off pathogens. Well investigated examples for such cross-talks among metazoans are the freshwater polyp Hydra, which seems to “actively select and shape its bacterial community” , and the symbiosis between the bobtail squid Euprymna and the bioluminescent bacterium Vibrio fisheri, which is selected and maintained by the host through PRR signaling (reviewed in ).
Amongst the basal animal phyla (Placozoa, Porifera, Cnidaria and Ctenophora), placozoans (for review see [9, 10]) clearly exhibit the simplest organized body plan, comprising only six somatic cell types identified so far [11,12,13]. Potential intracellular bacterial symbionts have been posited in these animals since ultrastructural analyses of Trichoplax adhaerens revealed bacterial enclosures in a subpopulation of fiber cells and within developing oocytes [11, 14, 15]. These findings suggested that placozoans are able to control the density as well as the distribution of their endosymbionts and that these are propagated to the next generation via vegetative and sexual reproduction. In 2013 Driscoll et al.  were able to identify the placozoan endosymbiont. Previously overlooked bacterial sequences in the Trichoplax adhaerens reference genome  were identified as fragments of a rickettsial endosymbiont genome. Rickettsiales belong to the α-Proteobacteria and constitute an order of obligate intracellular bacteria that are incapable to grow and divide outside the host because they have lost many important metabolic pathways [18, 19].
Recently we sequenced and investigated the genome of the placozoan 16S haplotype H2 (Trichoplax sp. H2; ), a close relative of Trichoplax adhaerens, and the genome assembly also yielded the complete genome of a placozoan rickettsial endosymbiont, which is described in Kamm et al., (in prep.). The results of the bacterial genome analyses, in terms of genome size, gene content and composition of metabolic pathways, clearly point to a mutualistic relationship which raises questions about how it has been established and how it is maintained. Relevant to understanding the placozoan holobiont is, if Trichoplax possesses an innate immunity system at all, and how it is organized to recognize harmful microorganisms on the one hand and to interact with beneficial microorganisms on the other. The small number of cell types and lack of specialized cells in placozoans further raises interesting questions about how immunity works in these most simple animals.
Recognition of non-self molecular patterns by the innate immune system takes place either on the extracellular side, via secreted or membrane bound PRRs, or inside a cell via intracellular receptors. A well investigated component of extracellular recognition in vertebrates and invertebrates includes the Toll-like receptors (TLR) , while intracellular recognition seems mainly mediated by the nucleotide binding and oligomerization domain like receptors (NOD-like), the retinoic-acid inducible gene I like receptors (RIG-I-like) or the stimulator of interferon genes receptor (STING) in conjunction with the cytosolic DNA sensor cyclic GMP-AMP synthase (cGAS) [22,23,24]. Other immunity related factors involved in extracellular recognition include so called scavenger receptors , which recognize a variety of ligands, and secreted lectins, that specifically bind to the carbohydrate moieties of bacterial cell walls. Particular classes of the latter are also able to initiate the complement system via the lectin-pathway . None of these components have been described for placozoans. This represents a significant gap since deviations from the well investigated immunity mechanisms of vertebrates appear to be quite common among Metazoa. Deviations may include the absence of particular immunity genes or pathways due to the emergence of a lineage before these had evolved, but also lineage specific gene losses or expansion of gene families. The latter mirrors the various phyla’s different evolutionary routes to accomplish recognition and management of harmful or beneficial microorganisms. We therefore examined the Trichoplax sp. H2 genome for genes involved in innate immune recognition and downstream signaling. The Trichoplax adhaerens reference genome  has not been specifically included in our results because both genomes are highly related  and the published gene models of the latter were found to be less complete and also inferior with respect to multi-domain genes. Overall, the reference genome did not contribute more information to the presence or absence of immunity related genes.
Results and discussion
True toll-like receptors are absent in Placozoa but TLR-like recognition may be organized differently
TLRs are transmembrane (TM) receptors that comprise extracellular leucine-rich-repeats (LRRs), capable of selectively binding foreign ligands like lipoproteins, lipopolysaccharides, peptidoglycans, ssRNA, dsRNA, (unmethylated) CpG-DNA or flagellin, depending on the receptor subtype . Upon binding of the ligand, the cytoplasmic Toll/interleukin-1 receptor (TIR) homology domain associates with the TIR-domain of the adapter protein MyD88, which itself recruits an Interleukin-1 receptor-associated kinase (IRAK) via the association of both of the protein’s Death domains (DDs) . Branches of downstream effectors terminate with the activation of the transcription factors NF-κB, AP-1 or Interferon regulatory factor (IRF), leading to the expression of inflammatory cytokines and other immune genes .
To determine if TLR signaling is present in Placozoa, we screened the predicted proteins of the Trichoplax sp. H2 genome for TLR pathway orthologs via BLAST and KEGG-, EggNOG- and OrthoMCL-mapping, and further validated putative orthologs by phylogenetic analyses and the presence of the necessary domains using the HMM based domain classification of InterProScan 5 . The latter was also used to identify and characterize TIR-domain containing predictions. In addition, modeling with HHpred  was used for certain candidate genes as a more sensitive approach to verify the presence or absence of domains that are otherwise difficult to detect (e.g. because of a higher divergence from the consensus in placozoans).
We found 16 proteins containing a TIR-domain (Additional file 1: Figure S1). Eight of these belong to a class of “evolutionary conserved TIR-domain containing” (ecTIR-DC) proteins that appear to be absent in vertebrates but are widespread among invertebrates and for which a function as adaptors or regulators in immunity-related signaling has been hypothesized . Among these are the only two placozoan TIR-domain proteins that also contain LRRs. Both gene products (classified as ecTIR-DC 14 II in ) are 3000 amino acids (AAs) large, lack a TM domain and contain several other domains, making any functional prediction speculative.
While these analyses detected no canonical TLRs in Trichoplax sp. H2, we further examined if recognition via TLRs may be similarly organized as in Hydra, where two separate membrane proteins appear to serve the same function: one containing the extracellular LRRs, the other containing the intracellular MyD88 recruiting TIR-domain [2, 31]. We found nine proteins with extracellular LRRs and a predicted single TM, another potential candidate additionally contained immunoglobulin-like domains (Additional file 1: Table S1). Modeling with HHpred revealed no further intra- or extracellular folds or domains, potentially pointing to a different function, and in two proteins the predicted secondary structure was found to be most similar to that of bilaterian TLRs. For the role of a cytoplasmic transducer we found no TIR-only prediction with a TM. Instead, we found a possible candidate in a SEFIR-domain containing transmembrane protein (Additional file 1: Figure S1). SEFIR-domains are related to TIR-domains and are also thought to be involved in homotypic interactions with other TIR/SEFIR-containing proteins . However, SEFIR-domains are also typical for interleukin-17 (IL-17) receptors which are further characterized by extracellular fibronectin-III-like domains. In the case of the Trichoplax SEFIR-domain containing transmembrane protein, neither fibronectin-domains nor other immunoglobulin-like folds were detected by InterProScan and Blast searches did not reveal the slightest resemblance to IL-17 receptors. On the contrary, HHpred detected similarity to IL-17 receptors (and weaker similarity to a TLR) in the region of the SEFIR domain. A reported single hit of a fibronectin-domain to the extracellular part was far beyond significance (see Methods) but homology may nevertheless be considered suggestive, given the presence of the SEFIR domain. The relation of this placozoan gene to other known genes, as well as its function, thus remains difficult to deduce. For the potential role of MyD88 we identified two MyD88-like proteins containing a TIR-domain and the DD-related caspase-recruitment domain (CARD) and another MyD88-like protein containing a SEFIR-domain and a DD (Additional file 1: Figure S1). The latter domain composition resembles both, that of CIKS/CIKSL adaptor proteins involved in IL-17 signaling, and that of MyD88 . HHpred modeling found significant similarity to MyD88 in the region of the DD and to interleukin-17 receptors in the region of the SEFIR domain, which, again, complicates speculations about the protein’s function. However, it is also possible that some of the placozoan proteins that have been found to contain only TIR-domains, or those that belong to the widespread ecTIR-DC proteins, serve as alternative signaling adaptors to MyD88 [21, 30, 34].
We can thus summarize that placozoans possess a variety of TIR/SEFIR-domain containing proteins and putative LRR-receptors that may act in the context of TLR recognition and signaling. A situation resembling that of Hydra’s bipartite TLRs. However, although Hydra’s atypical TLRs are also supported by functional assays , we have to note that their proposed TIR-domain containing transmembrane proteins (HyTRR-1, HyTRR-2) have been found to likely adopt an extracellular immunoglobulin-like fold and thus seem rather related to interleukin-1 receptors . If this contradicts a role in TLR recognition or simply emphasizes a derived status of hydrozoans within Cnidaria (anthozoans possess true TLRs ), has yet to be resolved. In the case of placozoans, the recognition arm of TLR signaling remains speculative since true TLR orthologs are absent and functional data are missing. It is also possible that other, yet unknown, receptors provide the input for TLR signaling.
The presence of TLR pathway orthologs indicates primordial TLR signaling
We identified many orthologs of the downstream effectors of the canonical TLR pathway (Fig. 1; Additional file 1: Table S2; Additional file 2: Dataset S1; Additional file 3: Dataset S2), indicating the presence of components for primordial TLR signaling in Placozoa. There are, however, important exceptions: While the pathway leading to AP-1 via mitogen-activated protein kinases (MAPK) seems complete, the path to NF-κB lacks the NF-kappa-B inhibitor alpha (IκBα), which, after phosphorylation by inhibitor of nuclear factor kappa-B kinases (IKKs), enables translocation of NF-κB to the nucleus (reviewed in ). However, the only conserved domains present in IκBα and related proteins are ankyrin-repeats (ANK) and the detection of putative placozoan homologs may thus be impeded if they deviated substantially from the bilaterian consensus. It is also possible that an unrelated protein among the roughly 200 ANK containing proteins in Trichoplax sp. H2 fulfills this function. Furthermore, the here identified placozoan NF-κB-like protein (p100-subunit-like) deserves some attention. While InterProScan identified only ANK repeats, the analysis with HHpred supports complete Rel homology and Rel homology dimerisation domains (RHD, IPT). Interestingly, two adjacent predictions in the genome contain well recognizable RHDs only, which adds up to three clustered members of the NF-κB/Rel family in placozoans, but otherwise these two appear to have no relation to known members of this gene family. The lack of a DD in the NF-κB- p100-like protein further indicates that the function and regulation of NF-κB complexes in Placozoa deviates from canonical pathways.
Since sponges and Cnidaria do possess complete NF-κB (in current NCBI databases, cf. [36, 37]), the question if placozoans have a disrupted NF-κB or represent a stage before acquisition of the DD depends on the actual phylogeny at the base of the Metazoa, which at this point in time is controversial (e.g. see ). However, if Placozoa are basal in the phylogeny, then the evolutionary scenario would be that they represent the stage before acquisition of the DD. If on the other hand, Placozoa were intermediate to sponges and Cnidaria, then the inference would be a disrupted NF-κB.
The absence of TRIF and TBK1 orthologs further indicates that Placozoa lack the MyD88-independent pathway, consistent with the hypothesis of its emergence in chordates , but the most important deviation from the reference pathway is the unclear orthology of complete IRAK family genes. While we found several IRAK-related kinases via BLAST, KEGG- or EggNOG-mapping, all of them lack a DD which is thought to be responsible for association with MyD88 . Many of these IRAK related placozoan protein kinases lie in a tight genomic cluster and thus seem to be placozoan specific paralogs produced by gene duplication. However, while Porifera possess reasonable IRAK orthologs with a DD (e.g. A. queenslandica: XP_011406220.2) this seems not to be the case for Cnidaria (c.f. [37, 40]) where all IRAK-like genes lack a DD, at least the closest cnidarian hits for bilaterian IRAKs in current databases such as UniProt, NCBI or Compagen. Nevertheless, it is generally assumed that TLR signaling is present in Cnidaria and functional analyses in Hydra have shown that several downstream targets of TLR signaling are affected after silencing of MyD88 , indicating that the pathway is functional even in the absence of a clear IRAK ortholog and the same could be true for placozoans. In the absence of IRAK, other kinases containing a DD could serve similar functions. In the predicted gene models of Trichoplax sp. H2 we were, however, only able to identify one kinase containing a Death-like domain (DLD).
Trichoplax sp. H2 contains a rich and unique scavenger receptor repertoire
Scavenger receptors (SR) are another class of cell surface receptors whose functional diversity has only recently been recognized. They are able to bind a large variety of ligands and their primary role seems to be the clearance from unwanted compounds like cellular debris, pathogens or other foreign particles by endocytosis . SRs constitute a structurally heterogeneous group, but most receptors associated with scavenger receptor activity are characterized by the presence of either C-type lectin domains (CTLD), scavenger receptor cysteine-rich domains (SRCR) or the CD36 domain. Generally, the SR repertoire seems to be greater in invertebrates, likely because of the absence of an adaptive immune system [3, 42, 43]. Their capacity to manage internalization of extracellular ligands by endocytosis makes them suitable candidates not only for detecting and fending off invasive pathogens but also to mediate symbiosis with intracellular bacteria.
We found a large repertoire of transmembrane proteins containing SRCR domains or CTLDs: Together with a single CD36 ortholog, the Trichoplax sp. H2 genome harbors 82 putative scavenger receptors, 45 containing SRCR domains, 36 CTLDs (Fig. 2; Additional file 4: Dataset S3A). Like in Cnidaria , no overlap between SRCR domain and CTLD domain containing proteins was observed. As expected for invertebrates, the number of scavenger receptors containing the above domains exceeds by far that of mammals (e.g. 11 in human) and also that of most Cnidaria so far investigated (12–38), with the exception of Aiptasia pallida (72) where it is more or less the same . However, the diversity of domain composition in placozoan SRs is greater than in the latter two animal groups. While the receptors containing either SRCR, CTLD or CD36 are composed of 6 different domains in mammals and 13 in Cnidaria [25, 43], more than 20 different domains are present in placozoan SRs and most combinations appear restricted to the phylum. Accordingly, most scavenger receptors found in Trichoplax sp. H2 cannot be assigned to well recognized mammalian classes, with the exception of Class B, E and I.
The majority of CTLD containing scavenger receptors in placozoans are G protein–coupled receptors
The SR complement of placozoans also differs in another unique aspect as 28 of 36 CTLD containing receptors bear seven transmembrane helices instead of a single transmembrane domain and are G protein–coupled receptors (GPCRs) of the Secretin/Adhesion family. To our knowledge these are the only described GPCRs whose extracellular domains consist of CTLDs. The one exception is Branchiostoma which harbors one adhesion GPCR containing an N-terminal CTLD, albeit together with other extracellular domains . This suggests that CTLD mediated recognition of molecular patterns in Placozoa is linked to GPCR signaling, trafficking and endocytosis . Furthermore, β-arrestin mediated endocytosis also couples GPCRs to the MAPK pathway. Hence, binding and internalization of targets could trigger a variety of cellular events, including defense or management of symbiotic microorganisms. In this context it is noteworthy that we found a large repertoire of arrestins (defined by the presence of N- and C-terminal Arrestin-like domains) in the Trichoplax sp. H2 genome (Additional file 4: Dataset S3B): in addition to one ortholog of β-arrestin, twenty three α-arrestins, related to the mammalian genes ARRDC2–4, are present, seventeen of which reside in a tight cluster on scaffold 76. Compared to mammals (two β-arrestins, two visual arrestins, six α-arrestins ), Nematostella (UniProt: one β-arrestin, three α-arrestins) and Amphimedon (UniProt: one β-arrestin, six α-arrestins) this gene number constitutes a unique expansion of the arrestin family in Placozoa. Even though the role of α-arrestins in GPCR signaling is less well understood, they are also known to interact with GPCRs in concert with β-arrestin . This observation indicates that GPCR mediated endocytosis in placozoans plays an important role for defense in particular and GPCR signaling in general.
Secreted proteins containing fibrinogen-related domains appear to play an important role in extracellular defense
Before pathogens come into contact with cell-surface receptors, the line of first defense consists of secreted lectin-like proteins that are able to recognize carbohydrates associated with bacterial cell walls and either agglutinate pathogens, or coat and opsonise them to enhance phagocytosis. Possible candidates, among others, are secreted proteins that contain CTLDs  or fibrinogen-related domains (FReDs) . Twelve predictions without a TM were identified harboring 1–3 CTLDs, eight of which contain a signal peptide (SP) and are thus likely to be secreted factors (Additional file 4: Dataset S3C). Fibrinogen-related domains on the other hand were previously claimed to be absent in placozoans . However, using the InterProScan annotation we were able to detect 18 FReD containing gene models with SUPERFAMILY member SSF56496 (Fibrinogen, C-terminal globular domain; Additional file 4: Dataset S3C). While 2 of these are large proteins related to fibrillin, the remainder are around 250 AA in length, bear a single FReD and an SP (except one) and are mostly organized in three clusters in the genome. The PANTHER classification system  classifies them as intelectins, a novel type of soluble lectins (X-type lectin), which seem remotely related to ficolins and have been shown to selectively recognize bacterial glycans . Moreover, an additional 15 gene models (11 with SP) are classified as intelectins by PANTHER and all of these, except 2, are genomic neighbors of the above, producing 5 clusters of 2–9 genes each. In 28 of these 31 genes, modeling with HHpred further supports the presence of a FReD (the remaining 3 predictions are probably truncated). This observation indicates that most of the putative intelectins belong to the same gene family but contain FReDs that are difficult to detect as a result of their phylogenetic distance to related sequences, which are the basis of canonical domain models. This could also be the reason why the identified FReDs in placozoan intelectins are comparatively short compared to vertebrate intelectins (up to 70 AAs versus 200 AAs).
Vertebrate and placozoan intelectins thus only align well in the region of the shorter placozoan FReD (restricted to the N-terminal portion of placozoan intelectins) while the latter also show high sequence similarity in the C-terminal portion (Additional file 5: Figure S2; Additional file 6: Dataset S4). Most placozoan sequences share 10 conserved cysteine residues (4–5 with vertebrate intelectins; complete alignment in Additional file 6: Dataset S4), suggesting that they are essential to form either intrachain or interchain disulfide bonds. The latter could serve to form intelectin oligomers, as is the case in human Intelectin-1 , which probably increases agglutination of trapped bacteria. While mammals possess up to 6 intelectins, teleosts have 9 and tunicates 22 intelectins , there are up to 31 intelectins in Trichoplax sp. H2. The apparent expansion of this gene family in placozoans suggests that intelectins play an important role in extracellular defense of the animal.
Some lectins, like the FReD containing ficolins and the CTLD containing mannose-binding protein (MBP), are known to initiate the lectin-complement-pathway  which represents an immediate response system for recognition and defense. In the cnidarian Aiptasia it has also been shown that the complement system is involved in the management of the cnidarian-dinoflagellate symbiosis . It is thus tempting to speculate if the placozoan intelectins or secreted CTLD proteins could serve a similar function. However, intelectins, as well as the placozoan secreted CTLD proteins, lack collagen repeats or coiled-coils that could mediate oligomerization similar to ficolins or MBP, which is needed to trigger the lectin-pathway. Moreover, placozoans apparently lack a complement system at all since most of its factors are absent. For example, neither KEGG pathway mapping nor screening for the necessary domains revealed genes for the key factors Mannan-binding lectin serine protease (MASP) or the complement factors B, C2, C3 or C4. With respect to the factors C3/C4, the only members of the Alpha-2-macroglobulin family present in Trichoplax sp. H2 are two CD109 genes which are not involved in the complement system. We only identified a single C1q-like factor (Additional file 4: Dataset S3C) bearing a clear C1q domain, and collagen repeats that are only recognizable by structural homology prediction with HHpred. In vertebrates, C1q plays an important function in the classical complement pathway . In invertebrates, which lack adaptive immunity, C1q factors are assumed to rather act as lectins in immune recognition and have undergone expansion in several taxa . The latter, apparently, does not apply to Placozoa.
Intracellular defense in Placozoa does not conform to the classical NOD-like receptor pathway and RIG-I-like receptors or cGAS-STING signaling are absent
NLRs mediate recognition once pathogens or their components have entered the cytoplasm [22, 54]. These proteins contain a central NOD domain (NACHT or NB-ARC) for nucleotide binding and self-oligomerization and C-terminal LRRs for molecular pattern recognition. An N-terminal effector domain mediates protein-protein interactions with downstream targets. Effector domains are members of the Death clan like PYD, CARD or the Death Effector Domain (DED). In the canonical NOD pathway, the receptors oligomerize after binding of PAMPs and expose their CARDs, enabling interaction with the CARD containing serine-threonine kinase RIPK2 which activates NF-κB and MAPK signaling . Other NLRs contain tetratricopeptide (TPR), WD40 or ANK repeats instead of LRRs and seem not involved in PAMP recognition (e.g. [54, 55]). RIG-I-like receptors detect viral RNA and are composed of a central DEAD/DEAH box helicase domain, an N-terminal CARD and a C-terminal regulatory domain . A further mechanism for detecting foreign cytosolic DNA is the cGAS-STING pathway .
While genes for cGAS, STING or RIG-I-like receptors were not found in Trichoplax sp. H2, 42 predicted genes containing a central NACHT/NB-ARC-domain (NACHT 3, NB-ARC 39) were identified (Fig. 3; Additional file 4: Dataset S3D). None of these contain LRRs and are thus unlikely to be involved in recognition of invading pathogens. In addition, no RIPK1/2 ortholog could be found in the Trichoplax sp. H2 genome, which is similar to the situation in cnidarians [55, 56]. The C-terminal repeats of placozoan NLRs are always either WD40 or TPR type repeats. Effector domains are predominantly CARDs or Death-like domains (DLDs) and in two cases DDs. These domain combinations of placozoan NLRs mostly resembles that of Hydra, which also does not contain any NOD-like receptors with LRRs, in contrast to anthozoans [54, 55].
Apaf-1 and controlled suicide as a means of defense
The NOD-like repertoire of placozoans also differs remarkably in another aspect to bilaterians and cnidarians: while most animals possess only one homolog of the Apoptotic protease activating factor-1 (Apaf-1), placozoans harbor several paralogs in their genome. If we consider DLDs as difficult to detect putative CARDs, Trichoplax sp. H2 harbors twelve Apaf-1 homologs, possibly even more because several Apaf-1 genes reside in clusters of NB-ARC containing genes but some gene models are missing either WD40 repeats or a CARD/DLD (Fig. 3; Additional file 4: Dataset S3D).
Apaf-1 plays a key role for the initiation of the mitochondrial pathway of apoptosis . It forms an oligomeric apoptosome and ignites the action of downstream caspases. The trigger is the released danger signal cytochrome c resulting from mitochondrial damage, a possible result of an infection. Apaf-1 thus enables controlled suicide of impaired cells. The death signal is passed from Apaf-1 to Caspase-9 via the association of both of the protein’s CARDs. Although H2 seems not to possess a clear Caspase-9 ortholog, it harbors two neighboring genes containing a caspase-like domain and an N-terminal CARD (Additional file 4: Dataset S3D), suggesting a similar function in apoptosis.
The presence of many Apaf-1 paralogs in placozoans suggests that controlled cell death after damage (e.g. by infection) plays a fundamental role in maintaining somatic and genetic integrity of the animal. Rather than triggering inflammatory responses, the cell is sacrificed once pathogens have escaped defense mechanisms and captured the cell. There are two reasons this strategy would make sense. First, the remarkable regenerative capacity of placozoans would be a contributing factor to how this hypothetical mechanism works and second, placozoans do not possess any organs whose accurate function has to be secured at all costs. Even the sacrifice of the entire animal seems plausible, if we assume that a local population mainly reproduces asexually and consists of genetically identical individuals. The selective advantage for a local genotype could thus lie in a rapid suicide of single infected animals to prevent spreading a disease to the entire population.
The immunity related pathways and genes found in the placozoan 16S haplotype H2 show several deviations from the common features of innate immune systems that have been characterized from phylogenetic branches of Bilateria to date. Due to the inconsistencies in related gene repertoires among other basal Metazoa (Porifera, Cnidaria) it is not always clear if these deviations (e.g. unclear orthology of IRAK and incomplete NF-κB) are the result of gene or domain loss, or because Placozoa simply predate the origin of such genes. Most remarkable is the presence of many genes for downstream TLR signaling while true Toll-like receptors are absent. Although we found candidate genes that may play a role in TLR recognition, this has to be proven experimentally. It is also possible that yet unknown receptors provide the input for TLR signaling in Placozoa. Other identified genes clearly constitute placozoan specific expansions of their immune repertoire (e.g. scavenger receptors, intelectins, Apaf-1). The distinct gene repertoire of placozoans further indicates that different metazoan phyla have embarked on different evolutionary routes to accomplish recognition of non-self in particular fields, but in general they use similar or the same core mechanisms. In any case, the complexity and amount of the found immunity genes suggests that Placozoa have sophisticated means to protect themselves against microbial invasion and to manage and maintain their associated microbiota - though experimental data is still needed for verification and to understand how they utilize this repertoire.
Functional annotation of Trichoplax sp. H2 predicted genes
The sequencing and annotation of the Trichoplax sp. H2 genome has been described in Kamm et al. (2018) . Briefly, gene prediction was carried out using the evidence based Maker pipeline  and functional annotation of the predicted proteins was carried out using InterProScan 5 (5.19–58.0)  with the following analyses: CDD-3.14, SignalP_EUK-4.1, PIRSF-3.01, Pfam-29.0, SignalP_GRAM_POSITIVE-4.1, TMHMM-2.0c, PRINTS-42.0, ProSiteProfiles-20.119, PANTHER-10.0, Coils-2.2.1, Hamap-201,605.11, ProSitePatterns-20.119, SUPERFAMILY-1.75, ProDom-2006.1, SMART-7.1, SignalP_GRAM_NEGATIVE-4.1, Gene3D-3.5.0 and TIGRFAM-15.0. Signal peptide (SP) and transmembrane (TM) predictions for innate immunity related genes were also carried out using the online server of Phobius . Annotation of predicted proteins also included BLASTP  searches against Swiss-Prot (cutoff e-value 1e-5) and KEGG pathway mapping using KAAS . Assignment to orthologous groups was done using the online servers of OrthoMCL (http://orthomcl.org/orthomcl/proteomeUpload.do ) and eggNOG-mapper (http://eggnogdb.embl.de/#/app/emapper ).
Identification of innate immunity related genes in Trichoplax sp. H2
For the identification and description of innate immunity genes in placozoans we preferentially used the predictions of the Trichoplax sp. H2 genome . Gene models of the Trichoplax adhaerens reference genome  were not specifically included because 1. both genomes are highly related and are expected to contain more or less the same genes, 2. the reference genome is less complete and 3. the corresponding predictions are less complete, especially regarding multi-domain genes. If orthologs of important immunity related genes could not be identified in Trichoplax sp. H2, the reference genome and the respective transcriptomes (see Kamm et al. 2018 ) of the two species were screened for their presence. Overall, neither the Trichoplax adhaerens reference genome nor the transcriptomes provided additional information.
Candidate orthologous genes of the Toll-like receptor (TLR) pathway were identified via KEGG mapping (see Kamm et al. 2018 ), best BLASTP hits of the predicted proteins against Swiss-Prot and by blasting respective metazoan orthologs against the predicted proteins. The candidate genes were further verified by mapping the predicted proteins to orthologous groups via OrthoMCL and EggNOG and by the presence of all necessary domains using the output of the InterProScan 5 functional annotation. The InterProScan results were also used to identify immunity and molecular pattern recognition related gene products based on their domain composition. Thus, in the functional annotation we looked for gene products containing (I) leucine-rich repeats (LRRs: IPR032675), the Toll-Interleukin receptor homology (TIR: IPR000157) or SEFIR (IPR013568) domain and could fulfill Toll-like receptor function; (II) the C-type lectin domain (CTLD: IPR001304), the scavenger receptor cysteine-rich (SRCR: IPR001190) domain or the CD36 (IPR002159) domain of putative scavenger receptors; (III) fibrinogen-related domains (FReDs: IPR002181); (IV) NACHT (IPR007111) or NB-ARC (IPR002182) domains of NOD-like receptors; (V) the C-terminal domain of retinoic acid inducible gene I (RIG-I: IPR021673) and (VI) the Mab-21 domain (IPR024810) of cyclic GMP-AMP synthase or IPR029158/IPR033952 of the stimulator of interferon genes protein for the presence of cGAS-STING signaling. Key complement factors were searched for via KEGG mapping and by screening the H2 proteins’ InterProScan annotation for the necessary domains (as defined by the InterProScan annotation of Homo sapiens orthologs in Swiss-Prot). For certain candidate immunity genes, modeling with HHpred , implemented in the MPI Bioinformatics Toolkit , was used as a more sensitive approach to ensure the presence or absence of domains that are otherwise difficult to detect because of a potential higher divergence from the consensus in placozoans. HHpred was run with default parameters and hits with a probability ≥90% and an E-value ≤0.1 were considered significant. Signal peptide and transmembrane region predictions were taken from the InterProScan 5 annotation (SignalP_EUK-4.1 , TMHMM-2.0c ) and also inferred using the Phobius web server .
Orthology validation of TLR pathway genes using phylogenetic analyses
The putative placozoan orthologs of the TLR pathway (Table S2) were subjected to phylogenetic analysis with potential orthologs from their respective gene family members. The potential gene family member sequences were obtained using a modification of the method outlined in Rosenfeld et al., 2016 . Briefly, a cDNA collection of complete and annotated genomes from NCBI was compiled using the taxa (and links) in Additional file 7: Dataset S5. These cDNA sequences were then processed with Transdecoder  to produce the best ORFs for each gene and collected into a single BLAST database. Because the human cDNAs contain the most reliable annotation, these ORFs were then queried against the full database of animal ORFs using TBLASTX to identify matches (cutoff value 1e-10). The most representative ORF for a gene in a particular taxon was determined by taking the BLAST hit with the lowest eValue. The matches for each of the genes per taxon were collected so that we had individual files with all of the orthologs of a particular human gene. Each of the TLR pathway genes of H2 was then used as a query for this database to give us a collection of putative orthologs for the H2 gene. This procedure resulted in orthologous groups with unique entries for a maximum number of animal taxa.
Phylogenetic matrices were built by combining the orthologous groups into a single matrix and by aligning the amino acid sequences with MAFFT (default settings). Poorly aligned regions in these alignments were then masked using Homo sapiens sequences as a guide. Phylogenetic analysis was accomplished using parsimony in TNT . We chose this method because it inherently relies on diagnostic changes to place terminals in a phylogenetic context . The TNT analyses used 1000 ratchet replicates saving 1000 trees at each replicate and consensus trees were used to summarize the phylogenetic searches. For a group of orthologs that make up a gene family (e.g. such as MEP1, TRAF1, TRAF2, TRAF3, TRAF4, TRAF6 and TRAF7) we combined the genes for these gene family members into a single matrix and analyzed them as a gene family unit.
Apoptotic protease activating factor-1
Cyclic GMP-AMP synthase
C-type lectin domain
Death Effector Domain
Evolutionary conserved TIR-domain containing
G protein-coupled receptor
Inhibitor of nuclear factor kappa-B kinase
Rel homology dimerisation domain
Interleukin-1 receptor-associated kinase
Interferon regulatory factor
NF-kappa-B inhibitor alpha
Mitogen-activated protein kinase
Mannan-binding lectin serine protease
Myeloid differentiation primary response 88
Nuclear factor kappa-B
Nucleotide binding and oligomerization domain
Pathogen-associated molecular pattern
Pattern recognition receptor
Rel homology domain
Retinoic-acid inducible gene I like receptor
Scavenger receptor cysteine-rich domain
Stimulator of interferon genes receptor
Toll/interleukin-1 receptor homology domain
Medzhitov R, Janeway CA. Decoding the patterns of self and nonself by the innate immune system. Science. 2002;296:298–300. https://doi.org/10.1126/science.1068883.
Miller DJ, Hemmrich G, Ball EE, Hayward DC, Khalturin K, Funayama N, et al. The innate immune repertoire in cnidaria--ancestral complexity and stochastic gene loss. Genome Biol. 2007;8:R59. https://doi.org/10.1186/gb-2007-8-4-r59.
Huang S, Yuan S, Guo L, Yu Y, Li J, Wu T, et al. Genomic analysis of the immune gene repertoire of amphioxus reveals extraordinary innate complexity and diversity. Genome Res. 2008;18:1112–26. https://doi.org/10.1101/gr.069674.107.
Huang S, Wang X, Yan Q, Guo L, Yuan S, Huang G, et al. The evolution and regulation of the mucosal immune complexity in the basal chordate amphioxus. J Immunol. 2011;186:2042–55. https://doi.org/10.4049/jimmunol.1001824.
Halanych KM, Kocot KM. Repurposed transcriptomic data facilitate discovery of innate immunity toll-like receptor (TLR) genes across Lophotrochozoa. Biol Bull. 2014;227:201–9. https://doi.org/10.1086/BBLv227n2p201.
Gerdol M, Luo Y-J, Satoh N, Pallavicini A. Genetic and molecular basis of the immune system in the brachiopod Lingula anatina. Dev Comp Immunol. 2018;82:7–30. https://doi.org/10.1016/j.dci.2017.12.021.
Chu H, Mazmanian SK. Innate immune recognition of the microbiota promotes host-microbial symbiosis. Nat Immunol. 2013;14:668–75. https://doi.org/10.1038/ni.2635.
Bosch TCG. What hydra has to say about the role and origin of symbiotic interactions. Biol Bull. 2012;223:78–84. https://doi.org/10.1086/BBLv223n1p78.
Schierwater B. My favorite animal,Trichoplax adhaerens. BioEssays. 2005;27:1294–302. https://doi.org/10.1002/bies.20320.
Schierwater B, DeSalle R. Placozoa. Curr Biol. 2018;28:R97–8. https://doi.org/10.1016/j.cub.2017.11.042.
Grell KG, Benwitz G. Die Ultrastruktur von Trichoplax adhaerens F. E Schulze Cytobiologie. 1971;4:216–40.
Jakob W, Sagasser S, Dellaporta S, Holland P, Kuhn K, Schierwater B. The Trox-2 Hox/ParaHox gene of Trichoplax (Placozoa) marks an epithelial boundary. Dev Genes Evol. 2004;214:170–5. https://doi.org/10.1007/s00427-004-0390-8.
Smith CL, Varoqueaux F, Kittelmann M, Azzam RN, Cooper B, Winters CA, et al. Novel cell types, neurosecretory cells, and body plan of the early-diverging metazoan Trichoplax adhaerens. Curr Biol. 2014;24:1565–72. https://doi.org/10.1016/j.cub.2014.05.046.
Grell KG. Eibildung und Furchung von Trichoplax adhaerens F.E. Schulze (Placozoa). Z Morph Tiere. 1972;73:297–314.
Grell KG, Benwitz G. Ergänzende Untersuchungen zur Ultrastruktur von Trichoplax adhaerens F.E. Schulze (Placozoa). Zoomorphology. 1981;98:47–67.
Driscoll T, Gillespie JJ, Nordberg EK, Azad AF, Sobral BW. Bacterial DNA sifted from the Trichoplax adhaerens (Animalia: Placozoa) genome project reveals a putative rickettsial endosymbiont. Genome Biol Evol. 2013;5:621–45. https://doi.org/10.1093/gbe/evt036.
Srivastava M, Begovic E, Chapman J, Putnam NH, Hellsten U, Kawashima T, et al. The Trichoplax genome and the nature of placozoans. Nature. 2008;454:955–60. https://doi.org/10.1038/nature07191.
Renvoisé A, Merhej V, Georgiades K, Raoult D. Intracellular Rickettsiales: insights into manipulators of eukaryotic cells. Trends Mol Med. 2011;17:573–83. https://doi.org/10.1016/j.molmed.2011.05.009.
Ettema TJG, Andersson SGE. The alpha-proteobacteria: the Darwin finches of the bacterial world. Biol Lett 2009;5 March:429–432.
Kamm K, Osigus H-J, Stadler PF, DeSalle R, Schierwater B. Trichoplax genomes reveal profound admixture and suggest stable wild populations without bisexual reproduction. Sci Rep. 2018;8:11168. https://doi.org/10.1038/s41598-018-29400-y.
Takeda K, Akira S. Toll-like receptors in innate immunity. Int Immunol. 2005;17:1–14. https://doi.org/10.1093/intimm/dxh186.
Kanneganti T-D, Lamkanfi M, Núñez G. Intracellular NOD-like receptors in host defense and disease. Immunity. 2007;27:549–59. https://doi.org/10.1016/j.immuni.2007.10.002.
Kawai T, Akira S. Toll-like receptor and RIG-I-like receptor signaling. Ann N Y Acad Sci. 2008;1143:1–20. https://doi.org/10.1196/annals.1443.020.
Margolis SR, Wilson SC, Vance RE. Evolutionary origins of cGAS-STING signaling. Trends Immunol. 2017;38:733–43. https://doi.org/10.1016/j.it.2017.03.004.
Canton J, Neculai D, Grinstein S. Scavenger receptors in homeostasis and immunity. Nat Rev Immunol. 2013;13:621–34. https://doi.org/10.1038/nri3515.
Fujita T. Evolution of the lectin-complement pathway and its role in innate immunity. Nat Rev Immunol. 2002;2:346–53. https://doi.org/10.1038/nri800.
Sun H, Bristow BN, Qu G, Wasserman SA. A heterotrimeric death domain complex in toll signaling. Proc Natl Acad Sci U S A. 2002;99:12871–6. https://doi.org/10.1073/pnas.202396399.
Jones P, Binns D, Chang H-Y, Fraser M, Li W, McAnulla C, et al. InterProScan 5: genome-scale protein function classification. Bioinformatics. 2014;30:1236–40. https://doi.org/10.1093/bioinformatics/btu031.
Söding J. Protein homology detection by HMM-HMM comparison. Bioinformatics. 2005;21:951–60. https://doi.org/10.1093/bioinformatics/bti125.
Gerdol M, Venier P, Edomi P, Pallavicini A. Diversity and evolution of TIR-domain-containing proteins in bivalves and Metazoa: new insights from comparative genomics. Dev Comp Immunol. 2017;70:145–64. https://doi.org/10.1016/j.dci.2017.01.014.
Bosch TCG, Augustin R, Anton-Erxleben F, Fraune S, Hemmrich G, Zill H, et al. Uncovering the evolutionary history of innate immunity: the simple metazoan Hydra uses epithelial cells for host defence. Dev Comp Immunol. 2009;33:559–69. https://doi.org/10.1016/j.dci.2008.10.004.
Novatchkova M, Leibbrandt A, Werzowa J, Neubüser A, Eisenhaber F. The STIR-domain superfamily in signal transduction, development and immunity. Trends Biochem Sci. 2003;28:226–9. https://doi.org/10.1016/S0968-0004(03)00067-7.
Rosani U, Varotto L, Gerdol M, Pallavicini A, Venier P. IL-17 signaling components in bivalves: comparative sequence analysis and involvement in the immune responses. Dev Comp Immunol. 2015;52:255–68. https://doi.org/10.1016/j.dci.2015.05.001.
Xu F, Zhang Y, Li J, Zhang Y, Xiang Z, Yu Z. Expression and function analysis of two naturally truncated MyD88 variants in the Pacific oyster Crassostrea gigas. Fish Shellfish Immunol. 2015;45:510–6. https://doi.org/10.1016/j.fsi.2015.04.034.
Hayden MS, Ghosh S. NF-κB, the first quarter-century: remarkable progress and outstanding questions. Genes Dev. 2012;26:203–34. https://doi.org/10.1101/gad.183434.111.
Gauthier M, Degnan BM. The transcription factor NF-kappaB in the demosponge Amphimedon queenslandica: insights on the evolutionary origin of the Rel homology domain. Dev Genes Evol. 2008;218:23–32. https://doi.org/10.1007/s00427-007-0197-5.
Song X, Jin P, Qin S, Chen L, Ma F. The evolution and origin of animal toll-like receptor signaling pathway revealed by network-level molecular evolutionary analyses. PLoS One. 2012;7:e51657. https://doi.org/10.1371/journal.pone.0051657.
Schierwater B, Holland PWH, Miller DJ, Stadler PF, Wiegmann BM, Wörheide G, et al. Never ending analysis of a century old evolutionary debate: “Unringing” the Urmetazoon bell. Front Ecol Evol. 2016;4:5. https://doi.org/10.3389/fevo.2016.00005.
Yang M, Yuan S, Huang S, Li J, Xu L, Huang H, et al. Characterization of bbtTICAM from amphioxus suggests the emergence of a MyD88-independent pathway in basal chordates. Cell Res. 2011;21:1410–23. https://doi.org/10.1038/cr.2011.156.
Hemmrich G, Miller DJ, Bosch TCG. The evolution of immunity: a low-life perspective. Trends Immunol. 2007;28:449–54. https://doi.org/10.1016/j.it.2007.08.003.
Franzenburg S, Fraune S, Künzel S, Baines JF, Domazet-Loso T, Bosch TCG. MyD88-deficient Hydra reveal an ancient function of TLR signaling in sensing bacterial colonizers. Proc Natl Acad Sci U S A. 2012;109:19374–9. https://doi.org/10.1073/pnas.1213110109.
Stuart LM, Ezekowitz RA. Phagocytosis and comparative innate immunity: learning on the fly. Nat Rev Immunol. 2008;8:131–41. https://doi.org/10.1038/nri2240.
Neubauer EF, Poole AZ, Weis VM, Davy SK. The scavenger receptor repertoire in six cnidarian species and its putative role in cnidarian-dinoflagellate symbiosis. PeerJ. 2016;4:e2692. https://doi.org/10.7717/peerj.2692.
Nordström KJV, Fredriksson R, Schiöth HB. The amphioxus (Branchiostoma floridae) genome contains a highly diversified set of G protein-coupled receptors. BMC Evol Biol. 2008;8:9. https://doi.org/10.1186/1471-2148-8-9.
Magalhaes AC, Dunn H, Ferguson SSG. Regulation of GPCR activity, trafficking and localization by GPCR-interacting proteins. Br J Pharmacol. 2012;165:1717–36. https://doi.org/10.1111/j.1476-5381.2011.01552.x.
Puca L, Brou C. α-arrestins - new players in notch and GPCR signaling pathways in mammals. J Cell Sci 2014;127 Pt 7:1359–1367. https://doi.org/10.1242/jcs.142539.
Kerrigan AM, Brown GD. C-type lectins and phagocytosis. Immunobiology. 2009;214:562–75. https://doi.org/10.1016/j.imbio.2008.11.003.
Hanington PC, Zhang S-M. The primary role of fibrinogen-related proteins in invertebrates is defense, not coagulation. J Innate Immun. 2011;3:17–27. https://doi.org/10.1159/000321882.
Thomas PD, Campbell MJ, Kejariwal A, Mi H, Karlak B, Daverman R, et al. PANTHER: a library of protein families and subfamilies indexed by function. Genome Res. 2003;13:2129–41. https://doi.org/10.1101/gr.772403.
Wesener DA, Wangkanont K, McBride R, Song X, Kraft MB, Hodges HL, et al. Recognition of microbial glycans by human intelectin-1. Nat Struct Mol Biol. 2015;22:603–10. https://doi.org/10.1038/nsmb.3053.
Yan J, Xu L, Zhang Y, Zhang C, Zhang C, Zhao F, et al. Comparative genomic and phylogenetic analyses of the intelectin gene family: implications for their origin and evolution. Dev Comp Immunol. 2013;41:189–99. https://doi.org/10.1016/j.dci.2013.04.016.
Poole AZ, Kitchen SA, Weis VM. The role of complement in cnidarian-dinoflagellate Symbiosis and immune challenge in the sea Anemone Aiptasia pallida. Front Microbiol. 2016;7:519. https://doi.org/10.3389/fmicb.2016.00519.
Gerdol M, Manfrin C, De Moro G, Figueras A, Novoa B, Venier P, et al. The C1q domain containing proteins of the Mediterranean mussel Mytilus galloprovincialis: a widespread and diverse family of immune-related molecules. Dev Comp Immunol. 2011;35:635–43. https://doi.org/10.1016/j.dci.2011.01.018.
Hamada M, Shoguchi E, Shinzato C, Kawashima T, Miller DJ, Satoh N. The complex NOD-like receptor repertoire of the coral Acropora digitifera includes novel domain combinations. Mol Biol Evol. 2013;30:167–76. https://doi.org/10.1093/molbev/mss213.
Lange C, Hemmrich G, Klostermeier UC, López-Quintero JA, Miller DJ, Rahn T, et al. Defining the origins of the NOD-like receptor system at the base of animal evolution. Mol Biol Evol. 2011;28:1687–702. https://doi.org/10.1093/molbev/msq349.
Baumgarten S, Simakov O, Esherick LY, Liew YJ, Lehnert EM, Michell CT, et al. The genome of Aiptasia, a sea anemone model for coral symbiosis. Proc Natl Acad Sci U S A. 2015;112:11893–8. https://doi.org/10.1073/pnas.1513318112.
Reubold TF, Wohlgemuth S, Eschenburg S. Crystal structure of full-length Apaf-1: how the death signal is relayed in the mitochondrial pathway of apoptosis. Structure. 2011;19:1074–83. https://doi.org/10.1016/j.str.2011.05.013.
Campbell MS, Holt C, Moore B, Yandell M. Genome annotation and curation using MAKER and MAKER-P. Curr Protoc Bioinformatics. 2014;48:4.11.1–39. https://doi.org/10.1002/0471250953.bi0411s48.
Käll L, Krogh A, Sonnhammer ELL. Advantages of combined transmembrane topology and signal peptide prediction--the Phobius web server. Nucleic Acids Res. 2007;35(Web Server):W429–32. https://doi.org/10.1093/nar/gkm256.
Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, et al. BLAST+: architecture and applications. BMC Bioinformatics. 2009;10:421. https://doi.org/10.1186/1471-2105-10-421.
Moriya Y, Itoh M, Okuda S, Yoshizawa AC, Kanehisa M. KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 2007;35(Web Server):W182–5. https://doi.org/10.1093/nar/gkm321.
Li L, Stoeckert CJ, Roos DS. OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003;13:2178–89. https://doi.org/10.1101/gr.1224503.
Huerta-Cepas J, Forslund K, Pedro Coelho L, Szklarczyk D, Juhl Jensen L, von Mering C, et al. Fast genome-wide functional annotation through orthology assignment by eggNOG-mapper. Mol Biol Evol. 2017. https://doi.org/10.1093/molbev/msx148.
Zimmermann L, Stephens A, Nam S-Z, Rau D, Kübler J, Lozajic M, et al. A completely Reimplemented MPI bioinformatics toolkit with a new HHpred server at its Core. J Mol Biol. 2018;430:2237–43. https://doi.org/10.1016/j.jmb.2017.12.007.
Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8:785–6. https://doi.org/10.1038/nmeth.1701.
Krogh A, Larsson B, von Heijne G, Sonnhammer EL. Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes. J Mol Biol. 2001;305:567–80. https://doi.org/10.1006/jmbi.2000.4315.
Rosenfeld J, Foox J, DeSalle R. Data for constructing insect genome content matrices for phylogenetic analysis and functional annotation. Data Br. 2016;6:279–81. https://doi.org/10.1016/j.dib.2015.12.015.
TransDecoder (Find Coding Regions Within Transcripts). https://transdecoder.github.io/.
Goloboff PA, Catalano SA. TNT version 1.5, including a full implementation of phylogenetic morphometrics. Cladistics. 2016;32:221–38. https://doi.org/10.1111/cla.12160.
Sarkar IN, Thornton JW, Planet PJ, Figurski DH, Schierwater B, DeSalle R, An automated phylogenetic key for classifying homeoboxes. Mol Phylogenet Evol. 2002;24:388–99. http://www.ncbi.nlm.nih.gov/pubmed/12220982.
Genome sequencing, assembly, annotation and comparative genomic analyses of Trichoplax sp. H2 were supported by grants from the German Science Foundation to B.S. (DFG Schi-277/26, 548 Schi-277/27, Schi-277/29).
Availability of data and materials
All data are either included in this published article and its supplementary information files or available in the genome publication of Trichoplax sp. H2 (Kamm et al. 2018 ). The annotated genome of Trichoplax sp. H2 has been deposited at DDBJ/ENA/GenBank under the accession NOWV00000000. The version described in this paper is version NOWV01000000. Individual genes or products described in this paper are indicated by their locus_tag.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. TIR/SEFIR-domain containing proteins in Trichoplax sp. H2. Table S1. Possible candidates for the extracellular part of a hypothetical bipartite TLR-like receptor in Trichoplax sp. H2. Table S2. The TLR pathway related gene repertoire in Trichoplax sp. H2 apart from TIR-domain containing genes. (PDF 83 kb)
Dataset S1. Trees of the phylogenetic validation of the TLR pathway genes. (PDF 1122 kb)
Dataset S2. Conserved domain alignments of the TLR pathway genes. (PDF 2446 kb)
Dataset S3. Spreadsheet containing the protein domain summary and the locus_tags of: A The Trichoplax sp. H2 scavenger receptors B The Trichoplax sp. H2 arrestins C The Trichoplax sp. H2 lectin-like secreted proteins (CTLD, FReD or C1q containing) D The Trichoplax sp. H2 NOD-like receptors. (XLSX 32 kb)
Figure S2. Vertebrate and placozoan intelectins aligned in the region of the fibrinogen-related domain (FReD) of placozoan intelectins. (PDF 117 kb)
Dataset S4. Complete alignment of the vertebrate and placozoan intelectins shown in Additional file 5: Figure S2 (Fasta format). (FA 15 kb)
Dataset S5. Spreadsheet with the list of taxa used for phylogenetic validation of the TLR pathway genes. Links are given for the NCBI download source of the respective cDNA databases. (XLSX 18 kb)
About this article
Cite this article
Kamm, K., Schierwater, B. & DeSalle, R. Innate immunity in the simplest animals – placozoans. BMC Genomics 20, 5 (2019). https://doi.org/10.1186/s12864-018-5377-3