Gene expression kinetics of Exaiptasia pallida innate immune response to Vibrio parahaemolyticus infection

doi:10.21203/rs.3.rs-24961/v2

Download PDF

Research article

Gene expression kinetics of Exaiptasia pallida innate immune response to Vibrio parahaemolyticus infection

https://doi.org/10.21203/rs.3.rs-24961/v2

This work is licensed under a CC BY 4.0 License

Journal Publication

published 09 Nov, 2020

Read the published version in BMC Genomics →

You are reading this older preprint version

Read the latest preprint version →

Background: Recent sequencing projects on early-diverging metazoans such as cnidarians, have unveiled a rich innate immunity gene repertoire; however, little is known about immunity gene regulation in the host’s early response against marine bacterial pathogens over time. Here, we used RNA-seq on the sea anemone Exaiptasia pallida (Ep) strain CC7 as a model to depict the innate immune response during the onset of infection with the marine pathogenic bacteria Vibrio parahaemolyticus (Vp) clinical strain O3:K6, and lipopolysaccharides (LPS) exposure. Pairwise and time series analyses identified the genes responsive to infection as well as the kinetics of innate immune genes over time. Comparisons between the responses to live Vp and purified LPS was then performed.

Results: Gene expression and functional analyses detected hundreds to thousands of genes responsive to the Vp infection after 1, 3, 6 and 12 hours, including a few shared with the response to LPS. Our results bring to light the first indications that non-canonical cytoplasmic pattern recognition receptors (PRRs) such as NOD-like and RIG-I-like receptor homologs take part in the immune response of Ep. Over-expression of several members of the lectin-complement pathways in parallel with novel transmembrane and Ig containing ficolins (CniFLs) suggest an active defense against the pathogen. Although lacking typical Toll-like receptors (TLRs), Ep activates a TLR-like pathway including the up-regulation of MyD88, TRAF6, NF-κB and AP-1 genes, which are not induced under LPS treatment and therefore suggest an alternative ligand-to-PRR trigger. Two cytokine-dependent pathways involving Tumor necrosis factor receptors (TNFRs) and several other potential downstream signaling genes likely lead to inflammation and/or apoptosis. Finally, both the extrinsic and intrinsic apoptotic pathways were strongly supported by over-expression of effector and executioner genes.

Conclusions: To our knowledge, this pioneering study is first to follow the kinetics of the innate immune response in a cnidarian during the onset of infection with a bacterial pathogen. Overall, our findings reveal the involvement of both novel immune gene candidates such as NLRs, RLRs and CniFLs, and previously identified TLR-like and apoptotic pathways in anthozoan innate immunity with a large amount of transcript-level evidence.

Epigenetics & Genomics

innate immunity

bacteria

infection

transcriptomics

Cnidaria

The immune system is an ancient and complex network that works to detect and defend an organism against pathogens, as well as to regulate the interactions between microbes and hosts. In vertebrates, the immune system consists of a two-fold mechanism: 1) the innate immunity which provides early protection to the host, and 2) the adaptive immune system capable of mounting a specific attack against invading pathogens and which displays memory. The current consensus is that invertebrates lack the components of the well characterized vertebrate adaptive immune system. However, with the recent proliferation of sequencing projects on early-diverging metazoans, the ancestral immune repertoire is clearly appearing more complex than previously suspected. As the sister group of bilaterians, cnidarians are the most early-diverging animals with true tissue layers and a gastro-vascular cavity, also known as the coelenteron. This game-changing internal feature is especially relevant to the evolution of innate immunity, as it marks the apparition of an internal location for interactions between host cells and microorganisms. Therefore, questions regarding the evolution of key bilaterian traits, such as pathogen and symbiont recognition, can be effectively explored using cnidarians. Indeed, cnidarians, including corals and anemones, have co-evolved since the Paleozoic era (570-245 Mya) with a multitude of unicellular algae, bacteria and viruses, and possess an innate immune system with an extensive list of genes homologous to those of the innate immunity of vertebrates (Hemmrich et al. 2007; Brennan and Gilmore 2018; Mansfield and Gilmore 2019)[1–3]. However, when looking at the composition of specific immune pathways such as the TLR-to-NF-κB pathway across evolutionary divergent genera of cnidarians, key genes can be missing (Brennan and Gilmore 2018). Therefore, comparative studies of the innate immune response in corals, Hydra, Nematostella, and Exaiptasia have the potential to help us understand how certain aspects of innate immunity evolved and may also point to alternative paths to activation of well conserved key factors of immunity.

Comparisons of the genomes of Acropora corals in symbiosis with algae, the aposymbiotic anemone Nematostella vectensis and the freshwater polyp Hydra vulgaris with those of other metazoans show increased complexity in gene families including Toll/Interleukin-1 Receptor- (TIR), NAIP/C2TA/HET-E/TP1- (NACHT), Tumor Necrosis Factor- (TNF) and Death-domains (DD) containing proteins (Miller et al. 2007; Shinzato et al. 2011; Hamada et al. 2013; Lange et al. 2010)[4–7]. In aiptasia’s (i.e. Exaiptasia pallida as renamed in Grajales et al. 2014)[8] transcriptome/genome (Sunagawa et al. 2009; Baumgarten et al. 2015, respectively)[9, 10], an active panel of immune-related genes is involved in the establishment of host-algae symbiosis (Poole et al. 2016)[11], and the disruption of this symbiosis also known as bleaching (Detournay et al. 2012)[12], which warrants the question whether the same is true against pathogen invasion.

Overall, studies of innate immunity in cnidarians have focused on its evolutionary history (Hemmrich et al. 2007; Miller et al. 2007), its role in bacteria sensing and aging (Franzenburg et al. 2012; Bosch 2013)[13, 14], self-/non-self-recognition (Dishaw and Litman 2009; Puill-Stephan et al. 2012)[15, 16], or within ecological contexts such as host-algae symbiosis (Davy et al. 2012)[17], bleaching (Weis 2019)[18], wound healing (van de Water et al. 2015)[19] and diseases (Mydlarz et al. 2006, 2016; Palmer and Traylor-Knowles 2012; van de Water et al. 2018)[20–22]. Although around twenty studies have investigated the innate immune response of cnidarian species to bacterial infection, only a few used the model Exaiptasia and none investigated the progress of the innate immune response over the onset of infection.

For example, Exaiptasia was used as a surrogate to better understand the diseases’ characteristics and interactions between coral pathogens such as Serratia marcescens (Sm) and Vibrio spp., and their hosts (Zaragoza et al. 2014; Alagely et al. 2011; Krediet et al. 2014)[23–25]. Later, Poole and colleagues (2016) used quantitative real-time PCR to show that genes of the complement pathway positively respond to both colonization with Symbiodinaceae and infection with Sm in aiptasia. More recently, Roesel and Vollmer (2019)[26] used RNA-seq to compare the transcriptomic response of symbiotic and aposymbiotic (menthol-bleached) Ep to 24 hours exposure to the bacterial pathogen Vibrio coralliilyticus (Vc). They found thousands of genes significantly differentially expressed suggesting roles of complement and coagulation cascades, NOD/Toll receptor signaling, and apoptosis in the innate immune response of Ep to Vibrio exposure.

Purified LPS - a Gram-negative specific endotoxin - has been used to induce a more controlled innate immune response in coral species (Fuess et al. 2017, 2016; Kvennefors et al. 2010; Williams et al. 2017)[27–29] and in anemones (Detournay et al. 2012; Perez and Weis 2006; Neubauer et al. 2016; Brennan et al. 2017)[12, 30–32]. For example, in the coral Acropora millepora, the expression of a lectin homolog dubbed Millectin was found to bind to pathogens and be activated under both LPS and peptidoglycan treatments (Kvennefors et al. 2010). Palmer and colleagues (2011)[33] exposed three coral species to the endotoxin LPS and showed that different parameters of innate immunity such as prophenoloxidase activity and melanin concentration could vary greatly between species. Very recently, Connelly and colleagues (2020)[34] used RNA-seq to show that the innate immune responses to LPS treatment in two congeneric coral species were genotype-specific and only involved the up-regulation of 6 transcripts, confirming the difficulty of studying the innate immune response in corals. In Exaiptasia, LPS and peptidoglycan failed to trigger significant changes in gene expression (Poole et al. 2016). Likewise, in 293-human cells, LPS did not activate a reconstituted TLR-to-NF-𝜅B pathway from both Hydra and Nematostella, although bacteria and flagellin did (Brennan et al. 2017; Bosch et al. 2009).

Here, RNA-seq sequencing on clones of the symbiotic Ep strain CC7 was used to identify the innate immune response of the host to bacterial infection. To challenge the innate immune system of Exaiptasia, the Gram-negative halophilic Vp was chosen because marine zoonotic pathogens such as Vp strain O3:K6 can have disastrous consequences for human populations. Indeed, VP is the leading cause of marine-borne illnesses worldwide (Livny et al. 2014)[35] causing gastroenteritis via the consumption of contaminated seafood (Czachor, 1992; Lipp and Rose, 1997)[36, 37]. On the other hand, environmental strains naturally occur in the mucus of coral species (Chimetto et al. 2008; Nithyanand et al. 2009)[38, 39], which make Vp a member of the rich bacterial community associated with cnidarians. This study focuses on the gene regulation at each time point and the gene expression kinetics of immune-related genes in response to 1, 3, 6 and 12-hour exposures to the pathogen, and contact with LPS for 3, 6 and 12 hours. Our results improve our understanding of innate immunity through the measurement of gene activity in response to infection in one of the most basally ancestral tissue-forming animals and highlights several fascinating gene targets for future functional studies.

Physiological response to infection

In order to study Ep’s innate immune response to Vp infection over time, we first visually monitored the response of anemones over 48 hours. Animals removed from the treatment medium after 12 or 24 hours, and subsequently left in filtered seawater (FSW) for another 36 or 24 hours, recovered their original appearances after the temporary assault (Supplementary Figure 1). Physiologically, Ep’s response to infection consists of condensed tentacles (Figure 1), increased mucus production, shriveling and cell detachment. Under the microscope, Vp in contact with the host can be seen swarming all over the outward surface of the anemones, which shows signs of tissue damage and invasion at 12 and 24 hours post-inoculation. After 24 hours in the infectious medium some mortality among Ep specimens was observed. To investigate the early innate immune response that allows Ep to defend itself against Vp’s attack, we focused our gene expression analyses on the first 12 hours of infection, when all anemone specimens remained alive.

Transcriptome assembly and annotation

Sample libraries - from 8 to 28 million raw reads across 23 samples (one sample failed quality control) - included both control and infected anemones at all four time points, and were assembled into 321,413 contiguous sequences or contigs, resulting in the holobiont de novo transcriptome (i.e. genome-free assembly). This transcriptome assembly was then aligned to the peptide sequences predicted from the aposymbiotic (deprived of algae) Ep strain CC7 genome from the Pringle laboratory at Stanford University (Baumgarten et al. 2015) to generate our Exaiptasia-specific reference transcriptome. This alignment identified 132,179 contigs (41% of the holobiont transcriptome) specific to Ep using BLASTX, which were then exclusively used in the rest of the analyses. Our reference transcriptome includes contigs corresponding to non-overlapping fragments of the same transcript. Thus, contigs represent both complete and partial gene sequences with 1 to 47 ‘isoforms’ per gene (as determined by the TRINITY assembler). The length of contigs ranges from 201 to 49,073 bp with an average length of 950 bp. On average 74% of reads per sample align to the reference transcriptome (number of reads per sample mapping to the reference transcriptome are given in Supplementary Table 1).

Among the reference transcriptome contigs, 130,037 (98%) retrieve a match from the NCBI non-redundant database (e-value < 10^-3 and % similarity > 35), 84,369 (64%) were annotated with Gene Ontology (GO) terms, 17,192 with KEGG enzyme codes, and 61,980 with InterPro IDs. GO annotations for the whole reference transcriptome resulted in 9,240 unique GO terms. Among these, 2,761 belong to 82 different “immune_class” categories of the CateGOrizer tool (Supplementary Figure 2).

Transcriptome-wide trends

RNA-seq data was produced for 32 anemones corresponding to three biological replicates at each time point and for each experimental group: controls and infected anemones at the four time points (minus one control anemone which did not pass quality control), and anemones exposed to LPS only for 3, 6, and 12 hours. The data was analyzed using two methods: 1) a pairwise comparison of gene expression levels between treated and control anemones at each time point, and 2) a time series analysis of the changes in gene expression over time for each experimental group.

The general trends primarily revealed from the pairwise gene expression analyses between control and infected samples at each of the four time points show that: 1) more genes are up-regulated than down-regulated throughout the host response to infection; 2) the number of regulated genes increases with time; and 3) the regulatory activity in response to the infection over 12 hours ranges from ≤ 1% to 3% of the whole transcriptome.

Table 1. Proportions of DECs of interest across the response to Vp infection - Number of DECs (robust Exact test at FDR<0.05) detected at each time point for the holobiont transcriptome, the transcriptome of aiptasia alone, and genes belonging to the four GO classes of interest in this study. The percentage of DECs belonging to the aiptasia-specific response is given in parentheses next to each DECs number.

The proportion of the aiptasia-specific DECs within the holobiont (i.e. the host and its symbionts) response ranges from 70 to 83% across up-/down-regulated DECs at the four time points of the experiment (Table 1). The numbers of up- and down-regulated DECs are similar after 1-hour post-infection, but the percentage of up-regulated DECs dramatically increased to 81% of the Ep gene activity after 3 hours. This relationship weakens with time but remains high (73 and 63% of DECs up-regulated) over the next two time points. Among the aiptasia-specific DECs, 2% of up- and < 1% of down-regulated DECs represented the immune related GO terms. To complement the list of DECs potentially involved in the innate immune response of Ep, more permissive lists were created to include genes associated with the recognition of the pathogen, the transduction of a defensive signal and the degradation of the invader. Thus, 2-8% of the aiptasia-specific up-regulated DECs represented the receptor activity, signal transduction and peptidase activity GO classes, while the same GO classes were represented by 2-6% of the down-regulated DECs. The up- and down-regulated genesets were then compared to each other across time points to identify DECs unique or shared between the different time points via Venn diagrams (Supplementary Figure 3). These Venn diagrams show that: 1) there are more genes unique to each time point than genes detected twice or several times; 2) there are very few genes continuously regulated during 12 hours of infection; and 3) later time points share more genes. These trends suggest a progression in the host response even though anemones were kept in the same infectious medium over 12 hours. As Vp cells settle down and anemones concentrate them from the medium, host-pathogen physical interaction likely decreases with time, while signal transduction from first pathogen contact to effector activation progresses.

Progression of innate immunity over 12 hours

The Fisher’s Exact Test (FET) performed on the up- and down-regulated lists of DECs detected at each time point highlight many significantly enriched GO classes unrelated to the immune response of Ep, but rather its maintenance of homeostasis and the energetic cost of its defense against pathogenic invasion. Although those aspects of Ep’s response to bacterial infection are of interest as well, we chose to focus on innate immunity. Therefore, Figure 2 only shows the results of enrichment analyses performed on up-regulated DECs and GO categories related to different aspects of the immune response. The enrichment analyses’ results on down-regulated DECs are shown in Supplementary Figure 4.

The results of the enrichment analyses strongly support an innate immune response activity already detectable after 1 hour and involving up-regulated genes linked to the ‘Toll signaling pathway’ biological process (GO_BP) and ‘TNFR-binding’ molecular function (GO_MF) represented throughout the experiment (Figure 2; FET: FDR < 0.05; Supplementary Table 2 for statistical values). As the immune response progresses, the ‘regulation of apoptotic process’ becomes evident after 3 hours, and remains apparent through the rest of the experiment. Among the genes indicating apoptotic activity, some belong to the ‘regulation of: NF-κB transcription factor, interferon-beta production, and MDA-5 signaling pathway’ GO_BP - as well as the GO_MF ‘caspase binding’ highlighted at 3, 6, and 12 hours post-infection. The role of the TNF pathway in the immune response of Ep gathers even more support at 3 and 6 hours, with the enrichment of ‘TNF-mediated signaling’. The ‘binding to AP-1 adaptor complex’ is highlighted at these two time points as well. At 6 hours, ‘activity of endopeptidase’ and ‘endopeptidase inhibitor’ was detected for the first time and was evident again at 12 hours. Finally, after 12 hours of infection, while activity within TNF receptor, apoptosis, and proteolysis pathways is still strong, additional peptidase activities join in the host response. In summary, these results reveal a noticeable chronology in the immune activity of Ep against the pathogen, starting with the involvement of TNFRs, then the activation of apoptosis followed by proteolysis. The genes supporting these functional categories, as well as DECs without GO annotations but containing immunity-related conserved domains such as DD, Ig, lectin, leucine-rich repeat (LRR), NACHT, scavenger receptor (SR), TIR, and TNF represent an interesting source of novel innate immunity candidates.

Profiling of seven gene expression clusters

The time series analysis comparing the changes in expression of DECs responding to Vp infection over time identified seven profiles (Figure 3). These profiles consist of plots showing the median expression levels for groups or clusters of DECs ranging from 13 to 75 contigs. Expression profiles diverge based on the changes with time but also the overall levels of expression of the DEC clusters. For example, cluster 1, 2 and 3 show very similar trends over time but display widely different median levels of expression. The overall trends across the seven profiles seem to mainly diverge at the 3 and 12-hour time points, while median expression levels are specific to each cluster. When looking at the genes represented in each cluster, annotations with relevance to innate immunity were noticed in each cluster. For example, Cluster 1 includes DECs homologous to transcription factor MafB-like, NF-κB, and ETS-related transcription factor Elf-3; Cluster 3: TNFR superfamily member 27, TRAF 2-like, NF-κB p100 subunit-like, complement C2-like, and complement factor B precursor; and Cluster 6: bcl-2-like protein 1, TNFR superfamily member 1A-like, TRAF 3, Cell death protein 3, WD repeat domain phosphoinositide-interacting protein 2-like, Protein NLRC5, and Macrophage scavenger receptor types I and II. A full list of the DECs detected by the time series analysis and organized by cluster is given in Supplementary Table 3.

Identification of innate immunity genes

The pairwise gene expression analyses detected hundreds to thousands of DECs at each time point, mounting to a total of 6,842 unique DECs responding to the infection (Supplementary Table 4 and 5). These strong effects of the infection on the transcriptome of Exaiptasia are illustrated in Figure 4. Firstly, the effects of infection are clearly separating on dimension 1 with controls and infected samples distinctly separated at each time point. Secondly, this separation between treatments becomes greater and greater over time also on dimension 1. Bacterial infection is progressive, therefore as the pathogens invade the anemones’ tissue, the transcriptomic activity in the host, reflected here by the number of DECs detected, increases over time. In addition, a subtler influence of time is illustrated by the separation between the treatments on dimension 2. This separation is initially detectable at the 3-hour time point, and increases obviously at the 12-hour time point (Figure 4). The fold changes illustrated in Figure 4 are provided in Supplementary Table 6, which follows the organization chosen by Zhang and colleagues (2018)[40] in a relevant study that analyzed the transcriptomic response of the amphioxus - Branchiostoma belcheri - to Vp.

Among the most regulated and well supported (i.e. showing high similarity) DECs (Supplementary Table 6), several gene families including coagulation factor, lectin, NLR, SR, interferon regulatory factor (IRF), and TNF were highly represented. For example, the most highly regulated immune related DECs at 1-hour post-infection are the IRF1-like sequences (FC = 795.7, Robust Exact Test - FDR = 6.52E^-3; and FC = -98.1, Robust Exact Test - FDR = 4.09E^-2), which are accompanied by isoforms up-regulated at all three other time points (Supplementary Table 4). At 3 hours post-infection, the TNF receptor-associated factor 5 (TRAF5) -like DEC is down-regulated 230.6 times (Robust Exact Test, FDR = 1.79E^-2), but isoforms are up-regulated at 3, 6 and 12 hours (Supplementary Table 4). From the same gene family, but after 6 hours post-infection, the TRAF4 DEC shows the highest up-regulation by a fold change of 488.2 (Robust Exact Test, FDR = 1.04E^-3). In the absence of TLRs in the genome of aiptasia, genes containing SR, Lectin, or LRR domains belonging to the macrophage scavenger receptor (MSR), lectin or NLR families potentially represent the pathogen recognition capacity in Ep, discussed further below. Furthermore, ligands and receptors of the TNF family suggest the induction of potential cytokine dependent pathways possibly connecting to an apoptosis-like pathway reported below.

Apoptotic gene activity

Similar to the innate immunity genes illustrated in Supplementary Table 6, highly regulated genes associated with apoptosis are shown in Supplementary Table 7. DECs likely involved in the apoptotic process were already detected after 1-hour post-infection and further progression of the apoptotic signal was observed at each following time point. For example, at the 3-hour time point, a homolog of ‘v-ets erythroblastosis virus E26 oncogene’ is up-regulated 249.4 times (FDR = 3.68E^-5) and was still up-regulated at 6 and 12 hours (Supplementary Table 4). A homolog of the key apoptotic executor – CASP3 – is detected as the overall most up-regulated DEC of this analysis, with a fold change of 1250.5 after 6 hours post-infection (Robust Exact Test - FDR = 3.06E^-17), in addition to isoforms also detected at 3 and 12 hours (Supplementary Table 4). Briefly, the genes in Supplementary Table 7 suggest that both the extrinsic and intrinsic apoptotic pathways are playing a role in the anemone’s response. The former is likely represented by caspase 8 (CASP8), Ced3 (CASP9) and CASP3, and the latter by mitochondria-associated genes such as Bcl-2, Bcl-2-associated X protein (Bax), and Bcl-2 antagonist/killer (Bak), among others. DECs representing the Ced3 gene take part in the expression profiles of clusters 6 and 7, while DECs for Bcl-2 are represented by clusters 1, 3 and 6 (Figure 3). In addition, the apoptotic initiator Apaf-1 is also represented in the expression profiles of clusters 2 and 7 (Figure 3).

Lack of TLR-like pathway activation under LPS

In order to identify the genes potentially responding to the pathogen-associated molecular pattern molecule (PAMP) LPS, lists of DECs produced as a result of 3, 6 and 12 hours of LPS treatment were compared to the DECs responding to the Vp infection at the same time points. Only a few DECs were found in common to both treatments (the intersections in Figure 5), and neither TLR homologs nor orthologs of the Hydra’s HyLRR-2 and HyTRR-1 genes (Augustin et al. 2010)[41] were identified in the Ep’s responses to LPS or Vp (Supplementary Table 8 versus Table 4, respectively). Based on these comparisons (Figure 5), a TLR-independent TLR-like pathway possibly represented by homologs of MyD88, TRAF, TBK1, NIK, IRF, NF-κB and AP-1, is activated under Vp infection but not LPS exposure. On the other hand, exposure to LPS induces TNF family homologs such as TNFR superfamily member 16 (TNFRSF-mb16) and LPS-induced TNFα factor up-regulation after 6 and 12 hours, respectively. Moreover, searching for potential receptors of PAMP among the DECs triggered by LPS reveals an active regulation of many genes potentially involved in sensorial pathways including G-protein coupled- and various neuronal receptors (Supplementary Table 8).

This time series RNA-seq experiment was conducted to investigate the molecular innate immune response of the symbiotic Exaiptasia pallida strain CC7 to Vibrio parahaemolyticus strain O3:K6 infection over 12 hours. Analyses focus on three main aspects of the innate immune response of the host: 1) the genes differentially expressed in infected anemones, 2) the gene expression changes over the onset of infection, and 3) the comparison between the responses to Vp and exposure to the endotoxins LPS. Punctual pairwise gene expression and time series analyses of Exaiptasia’s transcriptomic response to the virulent pathogen Vp identified 6,842 differentially expressed contigs responsive to bacterial infection across four time points. Below we discuss pathway level kinetics that contribute pertinent information to previously described innate immune gene families as well as novel innate immune gene candidates.

Physiological response to bacterial challenge

Ep exposed to the clinical O3:K6 Vp strain at 10⁸ cfu/ml starts to show mortality after 24 hours. Similar to this cross-pathogenicity, another study using the gram-negative human pathogen - Serratia marcescens - at the same concentration showed that Ep has a survival rate of 70% after 24 hours (Krediet et al. 2014). Our observations of Ep’s physiological response to Vp also concur with the descriptions of the same species infected with the coral pathogens Vc and V. shiloi (Zaragoza et al. 2014). Although darkening of the tissues during those infections was linked to melanin production, our search for evidence of melanin biosynthesis at the transcript level failed to identify homologs in this dataset. Anemones exposed to Vp show darkening of their tissues, retraction of tentacles, and mortality after 24 hours post-inoculation. Thus, our results suggest that the clinical Vp strain O3:K6 responsible for gastroenteritis in humans remains pathogenic in saltwater at 27 degrees celsius, and can successfully infect and start killing aiptasia anemones after 24 hours. Brown and Rodriguez-Lanetty (2015)[42] determined that an exposure to 10⁸ cfu/ml of Vc at 30 °C for three days was a sub-lethal treatment from which 100% of Ep exposed could recover. A similar result was observed here but after only 12 hours of exposure to 10⁸ cfu/ml of Vp at 27 °C. Based on these results, Ep likely possesses an innate immune system capable of efficiently fighting back seriously challenging and varied pathogenic assaults.

Functional support for non-canonical cytoplasmic PRRs

Within cells, intracellular PRRs detect pathogens and are critical in the response to microbial components in the cytosol. Here, homologs of NLRs, NLRC-like genes were significantly regulated at each time point (Supplementary Table 4) and include up-regulated sequences with NACHT-like domains and LRR as well as a down-regulated sequence with trans-membrane domains (TMD). The NLRs are cytoplasmic proteins, which recognize bacterial peptidoglycans and trigger proinflammatory and antimicrobial immune response. Among the NLR family, the NLR family CARD domain containing protein 4 (NLRC4) is the best-known mammalian member capable of triggering the inflammasome upon detection of microbial ligands (Duncan and Canna 2018)[43]. Interestingly, a large and complex NLR repertoire is present in early-diverging metazoans, including cnidarians, which exhibit an especially rich family of NACHT/NB-ARC containing genes in comparison with many later-diverging metazoans (Lange et al. 2011; Hamada et al. 2013). In actinarians, the diverse repertoire of NLR-like genes and their novel domain architecture has been hypothesized to perform novel roles in immune signaling and/or pathogen recognition, which are yet to be described (Baumgarten et al. 2015; van der Burg et al. 2016)[44]. Among the 78 NLRC contigs regulated during Ep’s infection, some contain LRR motifs (Supplementary Table 4), suggesting a capacity to bind to PAMP. These genes warrant further functional investigations to determine whether they play a role in the recognition of pathogens within the cells.

Additional genes were identified as potential cytoplasmic PRRs, including interferon-induced helicase C domain-containing protein 1 (IFIH1) and RIG-I-like receptor 3 (RLR-3), which both contain the RIG domain of the RLRs and play a role in the intracellular recognition of viruses. In total eight IFIH1 sequences and three RLR-3 sequences were detected as up-regulated at 3, 6 and 12 hours (Supplementary Table 4). As the pathogen Vp acquired several new genomic islands through horizontal transfer, which contain phage-encoded virulence factors (Ceccarelli et al. 2013)[45], it is possible that the Ep genes mentioned above are induced by viral-like genetic material brought into the host cells by invading Vp. To our knowledge, this is the first report of infection-induced RLR in a cnidarian.

Involvement of the lectin-complement pathway

In the response to Vp infection, we found that anemones up-regulated Ficolin-1 homologs containing TMD and Ig (unlike bilaterian ficolins; van der Burg et al. 2016) at 3, 6, and 12 hours (Supplementary Table 4), as well as several components of the lectin-complement pathway. For example, although no regulation of MASP was detected, complement component C2, C3 and C4 were up-regulated over time. Moreover, aiptasia’s factor B (Ep_Bf-1) of the alternative pathway was up-regulated in the infected anemones at 3, 6 and 12 hours (Supplementary Table 4). Poole et al. (2016) used quantitative real-time PCR to show up-regulation of both Ep’s Bf-1 and MASP genes in response to 24-hour exposure to low and high concentrations of the pathogen Serratia marcescens, and concluded that concentration of microbes, and symbiotic state influence complement gene expression. The co-expression of TMD/Ig-containing Ficolin-1 and genes of the lectin-complement pathways detected here supports a potential role of CniFLs in the initial recognition of pathogens (Baumgarten et al. 2015), and highlights both the lectin and alternative pathways active roles in the defense of anthozoans against pathogenic bacteria. Nevertheless, it remains to be demonstrated whether the CniFLs bind to certain PAMP and can trigger the lectin-complement pathway, potentially stimulating phagocytes in aiptasia.

Several other lectins were detected during the response to Vp infection. For example, L-rhamnose-binding lectins (Ep_RBLs) contigs, homologs of collectin-12 (colec12), and a single C-type lectin LRR-containing sequence were up-regulated (Supplementary Table 4). RBL interact with various types of bacteria in fish and invertebrates to induce proinflammatory cytokines such as IL-1β, TNFα and IL-8, and enhance macrophage phagocytosis (Watanabe et al. 2009)[46]. The scavenger receptor lectin-like family is rich and diverse in anthozoans (Neubauer et al. 2016)[31], plays an important role in host-symbiont interactions during onset of symbiosis (Davy et al. 2012), and can bind to pathogens (Kvennefors et al. 2008)[47]. Our results highlight the diversity of lectins potentially involved in the recognition, potential activation of phagocytosis, and proinflammatory signaling upon contact with a pathogen in aiptasia.

TLR-like pathway without TLR

Here, neither homologs of canonical vertebrate TLR or IL-1R, nor orthologs of the Hydra’s HyTRR, HyLRR genes (Augustin et al. 2010) were identified in the Ep transcriptomic responses. TLR and IL-1R-like homologs have been identified in several other anthozoans (Poole and Weis 2014)[48], however many of those genes lack the ectodomain LRR and only possess TMD- and TIR-domains (Baumgarten et al. 2015). Aiptasia’s genome possesses homologs of the major PRR types including NLRs, RLRs, and C-type lectin receptors (CLRs) but lack TLRs. However, TLR-to-NF-κB pathway components are definitely present in Ep, as well as in A. digitifera and O. faveolata corals, Hydra and Nematostella (reviewed in Brennan and Gilmore 2018). Several of these genes include homologs with well conserved domains such as MyD88, TRAF6, IRFs, AP-1, and NF-κB, and were positively responsive to Vp infection (Supplementary Table 4). Functional conservation of the TLR-to-NF-κB was confirmed in Nematostella (Brennan et al. 2017) and O. faveolata (Williams et al. 2018)[29]. Furthermore, in Hydra, MyD88-deficient polyps were more susceptible to infection with the human pathogen Pseudomonas aeruginosa (Franzenburg et al. 2012). Taken together, these results suggest that in Exaiptasia, the recruitment of MyD88 and downstream signaling to activation of NF-κB is induced by Vp, and like in Hydra must involve different PRRs than prototypical TLR homologs. Based on these findings, we hypothesize that an ancestral TLR-to-NF-κB signaling pathway responsive to bacterial stimuli evolved before transmembrane TLRs in early metazoans.

Potential TNFR signaling to inflammation and apoptosis

Another important aspect revealed by this study is the activation of cytokine-dependent signaling pathways in the innate immune response to Vp, two of which could lead to apoptosis and an inflammatory response. Indeed, the regulation of the ‘TNF-mediated signaling pathway’ was enriched at 3, 6 and 12 h post infection. One of these pathways likely starts with the activation of the TNFR1, which is the receptor for the proinflammatory cytokine TNFα and the cytotoxic protein lymphotoxin-alpha in humans. The Ep_TNFR1 DEC (GenBank: KXJ17455) is up-regulated 223.9 times after 1 hour, then again at a lower level after 3, 6 and 12 hours. Several TRAFs including TRAF2, -5 and the MAP3K14 among other MAPK3 homologs, are potentially acting downstream of Ep_TNFR1 and are up-regulated at 6 and 12 hours. Notably, more support for the TNFR1 pathway is shown with the TNFAIP3 and the TNIP1 homologs that are both up-regulated at 3, 6, 12 hours. TNFAIP3 is an essential component of the ubiquitin-editing protein complex that ensures the transient nature of inflammatory signaling pathways, while TNIP1 inhibits the NF-κB activation and TNF-induced NF-κB-dependent gene expression by regulating the activity of TNFAIP3. Finally, the transcription factor AP-1 is up-regulated at each time point and may influence the immune response via activation of Ep_TNFR1.

The other cytokine-dependent pathway may involve the TNFRSF-mb16 (GenBank: KXJ29589) homologs, which contain the Death domain and are up-regulated at 3, 6 and 12 hours. The Ep_TNFRSF-mb16 sequence resembles the Fas receptor (FasR or TNFRSF6 gene) which plays a central role upon Fas ligand (FasL) binding in the physiological regulation of programmed cell death in humans via its interaction with FADD and CASP8. In the coral A. digitifera, thirteen TNFR-like sequences containing the DD were identified and proposed as mediators of apoptosis through caspase activation (Quistad and Traylor-Knowles 2016)[49]. Moreover, the conservation of the TNF-induced apoptotic response was demonstrated in the coral Acropora digitifera (Quistad et al. 2014)[50]. We also identified candidates potentially acting up- and down-stream of FasR, such as matrix metalloproteinase and FasLG homologs.

Apoptosis as defense mechanism

In aiptasia's response to infection, a functional enrichment for ‘regulation of apoptotic process’ could be detected at 3 hours. The major cell death pathways, including apoptotic cell death, are a crucial barrier against microbial infection (Bergsbaken et al. 2009; Lamkanfi and Dixit 2010; Zitvogel et al. 2010)[51–53]. For instance, in response to bacterial infection, apoptosis or programmed cell death is used in the host innate immune response to: 1) eliminate pathogens at the early stage of infection without emitting alarm signals, and 2) induce dendritic cells (DCs) to engulf apoptotic bodies containing infected microbes (Elliott and Ravichandran 2010)[54]. Among the 216 DECs supporting the enrichment, many key players of apoptosis as well as genes with potential connections to the apoptotic pathway were identified at each time point. For example, the Ep_CASP8 is up-regulated at 3 and 6 hours and likely functions within the cells. If interacting with the Ep_TNFRSF proteins, Ep_CASP8 would represent the extrinsic (receptor-mediated) pathway to apoptosis. On the other hand, the intrinsic (mitochondria-mediated) pathway is strongly supported by the up-regulation of ATF4, Bcl-2, Bcl-W, Bax, Bak, Apaf-1 as well as CASP9 homologs in Ep (Supplementary Table 4). Both pathways could engage the apoptotic executor - CASP3, which is up-regulated for the duration of the infection with a peak of 1,250 in fold change at 6 hours. These results strongly suggest that apoptosis plays an important role in the innate immune response of Ep against bacterial pathogens. As Ep is capable of recovering from a punctual exposure to Vp, it is probable that apoptosis plays a part in the removal of invading bacteria and the survival of the host. Further work will determine if and how Ep is capable of containing a bacterial invasion via programmed cell death.

Comparison between responses to Vp infection and LPS

The genes of the TLR-to-NF-κB like pathway identified in the Ep’s response to Vp infection were not activated by LPS. Hydra polyps exposed to LPS treatments activate the expression of the antimicrobial peptides in a dose-dependent manner, suggesting that LPS is an inducer of immune defenses in cnidarians (Bosch et al. 2009). Indeed, studies on corals have shown that LPS can trigger an immune response in several species (Fuess et al. 2016, 2017; Kvennefors et al. 2010)[28, 55], and more specifically genes of the TLR-to-NF-κB pathway (Williams et al. 2018; Connelly et al. 2020). However, LPS failed to stimulate an engineered Hydra LRR protein to activate NF-κB signaling in human cells, whereas flagellin did (Bosch et al. 2009)[56]. Similarly, LPS did not stimulate the sole Nematostella TLR to activate a reconstituted Nv-TLR-to-NF-κB pathway in 293T cells, while flagellin and heat-inactivated Vc did (Brennan et al. 2017). Finally, in an attempt to induce innate immune response in Exaiptasia, Poole and colleagues (2016) exposed anemones to LPS and peptidoglycan, but did not detect any significant changes in gene expression, and decided to use live Serratia marcescens instead. Interestingly, in this study the LPS-induced TNFα factor was up-regulated in response to both LPS and Vp. Overall, these results suggest that in cnidarians lacking TLR homologs, LPS can still trigger the production of cytokines and antimicrobial peptides independently of the TLR-to-NF-κB pathway.

Novel genes responsive to infection

Several of the genes similar to the potential novel actinarian immune genes (NG1, NG2, and NG3) proposed by van der Burg et al. (2016) based on domain architectures, are here annotated F-box only protein 11, 4 and F-box/LRR-repeat protein 14, and found up-regulated at all four time points. This study also reveals for the first time in a cnidarian that diverse G-protein coupled receptor (GPCR) homologs are highly regulated all along the infection and in response to LPS exposure (Supplementary Table 8). A variety of GPCRs are expressed in T cells and were recently linked to an important role in the mediation of immunity in humans (Mashaghi et al. 2016; Wang 2018; Lu and Cyster 2019)[57–59]. Thus Ep_GPCRs constitute fascinating molecular targets towards the identification of ancestral immune cells in early-diverging metazoans and/or the activation of innate immunity via sensory pathways.

Anthozoan innate immune response to pathogen infection

In the present study, we detected many more up-regulated DECs than down-regulated ones, but also the correlation between the stress intensity and the number of DECs was: the greater the stress over time, the greater the up-regulation (Table 1). Vidal-Dupiol and colleagues (2011a, b; 2014)[60–62] have looked at the physiological and gene expression response of the coral Pocillopora damicornis (Pd) under temperature-dependent infection with Vc. Among the immune genes detected in Pd exposed to non-virulent Vc many also took part in this anemone’s response to Vp including NF-κB, AP-1, MyD88, and C3. However, these genes show the opposite regulation under the virulent Vc treatment, leading the authors to conclude the exact opposite to what was found in our study, “the greater the stress, the greater the down-regulation”. Therefore, besides the biological differences between anemones and corals, the confounding effects of temperature-induced bleaching and Vc’s virulence, the durational difference (i.e. days versus hours), as well as the various virulence mechanisms used by the different pathogens (Vc versus Vp), all contribute to the contradictory relationships found in Vidal-Dupiol et al. (2014) and our experiment.

Heat shock proteins (HSPs) are also among the actively regulated genes during a significant portion of the onset of Vp infection (Supplementary Table 4), which suggests that they play a role in defending against pathogenic invasion. This supports the findings in an experiment to determine the capacity of Ep for immunological priming and memory, in which Brown and Rodriguez-Lanetty (2015) used 2D fluorescent gel electrophoresis and mass spectrometry to detect differentially expressed HSPs in anemones exposed to sub-lethal Vc treatment for four weeks. Here, HSP70 and HSP90 DECs were up-regulated at 3, 6 and 12 hours in the response to Vp (Supplementary Table 4). Higher HSP expression levels in primed anemones were suggested to be an evolutionary trait that helps anemones to resist successive encounters with pathogens as long as they occur within a six week period (Brown and Rodriguez-Lanetty 2015)[42] and were thought to be more likely relevant to immunological priming conditions than pathogen clearance. However, the transcriptional changes of HSP transcripts observed here also suggest a role in the early defense against pathogens.

Our data provide strong support for the involvement of complement and coagulation cascades, NOD-like receptor signaling, lysosome, and apoptosis throughout the 12 hours of infection (Figure 2; Supplementary Table 6 and 7). These were also the pathways most affected in the responses of Ep to Vc detected using RNA-seq to explore the genetic links between anthozoan-algal symbiosis and innate immunity (Roesel and Vollmer 2019). Similar to our results, more DECs were found up-regulated than down-regulated in response to 24 hours exposure to Vc (i.e. 2,155 versus 1,552 DECs, respectively). Using KEGG pathway enrichment analysis, one pathway was significantly over-represented: the complement pathway under Vc treatment only. Our results build on Roesel and Vollmer’s (2019) results by pointing to the functional equivalence to canonical PRRs in cnidarians as well as many potential downstream immune signaling genes.

To our knowledge, this is the first time series study of the transcriptomic response to bacterial infection in a cnidarian model. Thousands of differentially expressed sequences were responsive to the infection over 12 hours. Among the most interesting discoveries were indications that several innovations, including possible neofunctionalization, conservation of signaling pathways over conservation of particular key genes, and involvement of lineage-specific immunity genes, are represented in the innate immune response of aiptasia. Recognition of bacterial pathogens is likely fulfilled by putative novel immune receptors. Transduction of immune signals through a TLR-to-NF-κB like pathway is supported. Two TNFR pathways could possibly lead to an inflammatory response and apoptosis. And uncharacterized actinarian-specific genes with novel immune-domain architectures are responsive to bacterial infection. Altogether these results reveal the genes responsive to bacterial infection among innate immunity candidates in anemones, as well as potential novel players offering many new exciting research subjects to advance knowledge of innate immunity in early-diverging metazoans.

Experimental design

Individual clones of the sea anemone Ep strain CC7 were acquired from the Pringle laboratory at Stanford University and kept at the Scientific Center of Monaco under controlled conditions in the Ecosystems and Immunity Laboratory of the Biomedical Department. Anemones are maintained in three liter tanks inside a culture chamber (Percival Scientific, USA) at constant 27°C with 20 μmol.m^-2.s^-1 of light (12:12h). FSW at 0.45 μm is renewed every week and anemones are fed twice a week with artemias, except on the week of the experiment during which they fast.

For the infection experiment, a single anemone or biological replicate of approximately half a centimeter in diameter was placed in each well of eight 6-well plates filled with 10 ml of 0.22 μm FSW. Anemones were kept still overnight until fixed to the bottom.

Vp strain RIMD 2210633 serotype O3:K6 came from the laboratory of Prof. Kodama (Osaka University, Japan). Routinely, frozen Vp are cultured at 37°C in Luria-Bertani liquid medium (LB) supplemented up to a final concentration of 1.5 g/L NaCl. For the infection treatment, Vp was then grown again at 27°C in LB supplemented with NaCl to a final of 3.5 g/L and without agitation for 24 hours. The Vp cultures were centrifuged at 2,500 rpm for 10 minutes, the bacterial pellet resuspended in 0.2 μm FSW to be concentrated 10 times, and 500 μL were dispensed into the wells containing the anemones in 9.5 ml of FSW. The final concentration of Vp was 10⁸ cfu/ml. Preliminary kinetics experiments using this concentration of pathogen were performed to choose relevant time points at which samples would be fixed for further molecular analyses. For the LPS treatment, ultrapure LPS from E. coli K12 (Invivogen cat. tlrl-peklps) was added to reach a final concentration of 10 ug/ml in the medium containing the anemones.

After Vp or LPS added to treatment wells, anemones were placed in the Percival chamber at 27 °C under gentle orbital agitation of 500 rpm for 1, 3, 6 and 12 hours. At each time point, whole anemones were detached, transferred in 1.5 ml tubes and frozen in liquid nitrogen for further molecular analyses. The LPS treated samples were collected only after 3, 6 and 12 hours, whereas the Vp treatments also were sampled at 1 hour in addition to the other three time points. Three independent anemones (n=3) were collected at each time point for each condition to result in 33 samples total (Vp exposure: 3 biological replicates x 2 conditions x 4 time points and LSP: 3 biological replicates x 1 condition x 3 time points).

Monitoring of physiological appearance

Photographs were taken using a stereomicroscope Stemi 305 trino (Zeiss) mounted with a digital camera 5MP Toupcam (Touptek). Photos were taken at each time point under bright light and the 0.8X magnification to monitor the overall appearance of the animals. Further observations were recorded at 5X magnification to assess the progression of the infection and the integrity of the animal tissue. Animals were considered dead when not contracting after being touched.

RNA isolations

Frozen anemone tissues were disrupted and cells lysed in Trizol reagent (Ambion life technologies) to isolate RNA. To do so, approximately 100 uL of 1mm diameter glass beads in 1.5 ml tubes filled with 1 ml of lysis buffer were used for anemones weighing approximately 100 mg. Tissue was completely disrupted using a Precellys 24 tissue lyser (Bertin Technologies) at 500 bpm for 2x 5 seconds. Samples were then immediately placed on ice and total RNA extraction protocol was executed as follows:

Two chloroform phase separations were performed on ice in order to precipitate the polysaccharides from the aqueous phase as much as possible. 300 uL of chloroform was added to the tissue lysate the first time, hand shaken for 15 seconds and settled at room temperature for 5 minutes, then centrifuged at 13 000 g for 15 minutes at 4°C. 200 uL of chloroform was added to the aqueous phase the second time and the previous steps repeated. The second aqueous phase was then added to an equal volume of chilled 100% ethanol and transferred onto a column for purification. The purification, including a DNase treatment step, followed the Direct-zol RNA MiniPrep Plus protocol (cat# R207 from Zymoresearch). RNA samples were then sent to Eurofins (Europe) for quality control on the Bioanalyzer (Agilent), cDNA library preparation using polyA selection (96 RNA-seq library from total RNA protocol, Illumina) and sequencing. cDNA libraries were sequenced as 2x 100 bp paired-end reads on the Illumina HiSeq2500 with chemistry v4 and high-output run mode.

Sequencing data processing

Raw sequences in fastq files were processed to remove Illumina adapters, low quality bases and short reads using the Trimmomatic-0.36 tool (Bolger et al. 2014)[63]. New reference transcriptomes were assembled de novo using all of the Exaitasia pallida libraries produced here and either guided by the Ep genome produced by the Pringle Laboratory (version 1.1) or genome-free using Trinity v2.6.5 release (Grabherr et al. 2011)[64]. Contiguous sequences or contigs were used to retrieve gene descriptions using blatsx against the non-redundant NCBI database with the following filter: E-value < 10^-3, and only BLAST description annotation excluding the terms ‘unknown’, ‘hypothetical’, or ‘uncharacterized’. In addition, coding regions were predicted using TransDecoder and the predicted peptide sequences were used to retrieve further gene annotations from several different databases using BLASTP against Uniprot, Pfam, Swissprot and UniRef90. HMMER and the Pfam database were also used to identify conserved protein domains present in our predicted peptide database. Depictions of domain architecture were based on conserved Domain Architecture Retrieval Tool (NCBI) results. Transmembrane domains were identified via the TMHMM (v.2.0) of the DTU Bioinformatics server. Contigs belonging to the anemone were isolated from the rest of sequences using BLASTX against the Ep predicted protein database produced by the Pringle Laboratory (Baumgarten et al. 2015). Transcript-level quantification was estimated for gene and isoform expression by allocating multi-mapping reads among all transcripts of the reference transcriptome using the RSEM (Li and Dewey 2011)[65] and Bowtie-2 (Langmead and Salzberg 2012)[66] packages.

Two different gene expression analyses were performed on the count data: 1) a pairwise differential expression analysis (PDEA) comparing infected or LPS treated samples to controls using the ‘edgeR’ Bioconductor package (Robinson et al. 2010)[67], and 2) a time course expression analysis (TCEA) detecting genes with significant expression profile differences using the ‘maSigPro’ Bioconductor package (Conesa et al. 2006)[68], both implemented in the OmicsBox (Version 1.3.11) platform. Genes with low read counts across libraries (i.e. samples) were filtered out on a count-per-million basis equal to 1. The raw library sizes were normalized using the Trimmed mean of M values method from edgeR. For the PDEA, a simple design was used to compare treated versus control conditions at each time point. The Exact Test (in edgeR) with the strengthened estimation against outlier features option was used to detect differential expression. For the TCEA, a multiple series time course design was used to detect differentially expressed genes using a false discovery rate correction < 0.05 for variable selection in the stepwise regression. A cutoff value for the R-square of the regression model was set to 0.9. The optimal number of clusters (i.e. similar gene expression profiles) was determined using the model-based clustering method (in maSigPro 1.58.0). Results according to the selection procedure consisted in: 1) genes with differential expression at the first time point, 2) genes with a linear differential expression over time, 3) genes with a transitory differential expression, 4) genes with differential expression over time in controls, and 5) genes with differential expression profiles between treated and control samples.

The pairwise differential expression analysis comparing the transcript abundance in control anemones to anemones in contact of the pathogen was performed using edgeR (Robinson et al. 2010) on the Blast2GO platform (Götz et al. 2008)[69]. The exact test based on the quantile-adjusted conditional maximum likelihood method for single-factor experimental design was applied here to detect DECs between control and infected anemones at different time points. Command lines used for the transcriptome assembly are provided in Supplementary File 1. The tables of DECs in Supplementary Table 2 and 3 show DECs with fold changes > 5 or < -5, homolog peptide sequence similarity values > 70% and represent subsets of Supplementary Tables 3 and 4.

Functional analyses

Immune class GO categorization

To visualize the progress of aiptasia’s immune response to the infection, GO terms for differentially expressed contigs at each time point were categorized according to the ‘immune-class’ categorization method by Zhi-Liang Hu (2008)[70]. Briefly, the CateGOrizer tool allows simplifying the visualization of functional annotations from gene lists by counting all paths between child and parent GO terms. Here, all paths between child and parent GO terms (using the cumulative count method: multiple paths possible between a single child and a single parent) were counted for the whole transcriptome as well as for the lists of DECs corresponding to each time point. The percentage of ‘immune_class’ GO term paths was then calculated as follows: GO path % = nbr. of GO paths in DECs list for GO:000xyz / nbr. of GO paths in transcriptome for GO:000xyz.

Targeted Gene Ontology

The GO codes corresponding to the parent terms ‘immune response’, ‘receptor activity’, ‘signal transduction’, and ‘peptidase activity’ were used to count all corresponding child terms within the DECs list for each time point. These DECs were then compared between each time point and visualized via a Venn diagram (http://bioinformatics.psb.ugent.be/webtools/Venn/ ; Supplementary Figure 2). The largest fold changes in expression level for the immune response related DECs analyzed in Supplementary Figure 2 are shown in Table 2. The complete list of DECs containing sequence IDs, hits information, folds changes and statistical results is accessible in the Supplementary Table 3.

Enrichment analyses

To detect a significant enrichment for the immunity associated GO terms within the lists of DECs, two tests were applied: one quite permissive – the Fisher’s exact test, and the other more stringent – the gene set enrichment analysis method. Both tests were run using the Blast2GO platform (Götz et al. 2008) and the whole transcriptome annotated with GO terms as the reference gene set. For both tests, only contigs with false discovery rate corrected p-value lower than 0.05 were considered. The enrichment analysis using the Fisher’s exact test was performed on the up- and down-regulated gensets produced via the gene expression analyses for the 1h, 3h, 6h and 12h time points. Results were then reduced to the most specific GO categories applying a lower than 0.05 filter on FDR values to produce Figure 3. Some limitations of the functional analysis (FET) need to be kept in mind and concerns the level of fragmentation of our transcriptome, meaning that several contigs can belong to the same transcript. This can introduce a bias when running enrichment analysis which compares the proportion of genes representing a gene ontology category in a test geneset, with the proportion of genes for that same GO category in a reference geneset. In a case where, a test geneset - made of a large number of DECs representing a ‘i’ number of transcripts is compared to a different geneset representing the same number of transcripts but by fewer DECs: the Fisher’s exact test is more likely to find GO categories significantly represented in the former geneset than the latter. Therefore, the enrichment is somewhat dependent on the level of fragmentation of the genes represented in the geneset. The enrichment for the GO categories discussed here were manually validated and are represented by many relevant genes provided in Supplementary Table 3 and 4.

Exaiptasia pallida (Ep), Vibrio parahaemolyticus (Vp), lipopolysaccharides (LPS), pattern recognition receptors (PRRs), NOD-like receptors (NLRs), RIG-I-like receptors (RLRs), cnidarian ficolin-like (CniFL), Toll-like receptors (TLRs), Serratia marcescens (Sm), Vibrio coralliilyticus (Vc), leucine-rich repeat (LRR), filtered sea water (FSW), colony forming unit (cfu), Gene Ontology (GO), differentially expressed contigs (DECs), false discovery rate (FDR), Fisher’s exact test (FET), biological process (BP), molecular function (MF), cellular component (CC), scavenger receptors (SR), multidimensional scaling (MDS), log Fold Change (logFC), pathogen-associated molecular pattern (PAMP), transmembrane domain (TMD), rhamnose-binding lectins (RBLs), C-type lectins (CLRs), dendritic cells (DCs), novel genes (NGs), G protein-coupled receptors (GPCRs), Pocillopora damicornis (Pd), Luria-Bertani liquid medium (LB), pairwise differential expression analysis (PDEA), time course expression analysis (TCEA)

Ethics approval and consent to participate

Not applicable

Consent to publish

Not applicable

Availability of data and materials

The datasets generated during the current study will be publicly available in the NCBI SRA repository under the submission number: 8148661, on January 1^st 2021.

Competing interests

The authors declare that they have no competing interests

Funding

Funding was entirely provided by the government of Monaco

Authors' contributions

F.S. designed and ran the experiment, performed the bulk of the computational analyses, interpreted the results and wrote the manuscript. D.D. performed data mining for novel domain architectures. D.C. and L.B. supervised all aspects of the project. All authors have read, improved and approved the manuscript.

Acknowledgements

We are grateful to Dr. Emmanuel Lemichez for his participation in early discussions on the experimental design. We are grateful to Dr. John Pringle for sharing the Ep clone CC7 with us. We also thank Mss. Alexia Verduci and Mr. Rodolphe Pontier-Bres for their help in the laboratory, and the Professors P. Rampal, D. Allemand and Dr. P. Warner for reading and commenting on the manuscript. We thank Eurofins Genomics for the sequencing. D.D. was financially supported by a Master II fellowship funded by the CSM as part of a collaboration between the Ecosystems and Immunity team (CSM) and the Université Côte d’Azur.

Hemmrich G, Miller DJ, Bosch TCG. The evolution of immunity: a low-life perspective. Trends in Immunology. 2007;28:449–54.
Brennan JJ, Gilmore TD. Evolutionary Origins of Toll-like Receptor Signaling. Molecular Biology and Evolution. 2018;35:1576–87.
Mansfield KM, Gilmore TD. Innate immunity and cnidarian-Symbiodiniaceae mutualism. Developmental & Comparative Immunology. 2019;90:199–209.
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.
Shinzato C, Shoguchi E, Kawashima T, Hamada M, Hisata K, Tanaka M, et al. Using the Acropora digitifera genome to understand coral responses to environmental change. Nature. 2011;476:320–3.
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. Molecular Biology and Evolution. 2013;30:167–76.
Lange C, Hemmrich G, Klostermeier UC, Lopez-Quintero JA, Miller DJ, Rahn T, et al. Defining the Origins of the NOD-Like Receptor System at the Base of Animal Evolution. Molecular Biology and Evolution. 2011;28:1687–702.
Grajales A, Rodríguez E. Morphological revision of the genus Aiptasia and the family Aiptasiidae (Cnidaria, Actiniaria, Metridioidea). Zootaxa. 2014;3826:55.
Sunagawa S, Wilson EC, Thaler M, Smith ML, Caruso C, Pringle JR, et al. Generation and analysis of transcriptomic resources for a model system on the rise: the sea anemone Aiptasia pallida and its dinoflagellate endosymbiont. BMC Genomics. 2009;10:258.
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 USA. 2015;112:11893–8.
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. doi:10.3389/fmicb.2016.00519.
Detournay O, Schnitzler CE, Poole A, Weis VM. Regulation of cnidarian–dinoflagellate mutualisms: Evidence that activation of a host TGFβ innate immune pathway promotes tolerance of the symbiont. Developmental & Comparative Immunology. 2012;38:525–37.
Franzenburg S, Fraune S, Kunzel S, Baines JF, Domazet-Loso T, Bosch TCG. MyD88-deficient Hydra reveal an ancient function of TLR signaling in sensing bacterial colonizers. Proceedings of the National Academy of Sciences. 2012;109:19374–9.
Bosch TCG. Cnidarian-Microbe Interactions and the Origin of Innate Immunity in Metazoans. Annu Rev Microbiol. 2013;67:499–518.
Dishaw LJ, Litman GW. Invertebrate Allorecognition: The Origins of Histocompatibility. Current Biology. 2009;19:R286–8.
Puill-Stephan E, Willis BL, Abrego D, Raina J-B, van Oppen MJH. Allorecognition maturation in the broadcast-spawning coral Acropora millepora. Coral Reefs. 2012;31:1019–28.
Davy SK, Allemand D, Weis VM. Cell Biology of Cnidarian-Dinoflagellate Symbiosis. Microbiology and Molecular Biology Reviews. 2012;76:229–61.
Weis VM. Cell Biology of Coral Symbiosis: Foundational Study Can Inform Solutions to the Coral Reef Crisis. Integr Comp Biol. 2019;59:845–55.
van de Water JAJM, Ainsworth TD, Leggat W, Bourne DG, Willis BL, van Oppen MJH. The coral immune response facilitates protection against microbes during tissue regeneration. Mol Ecol. 2015;24:3390–404.
Mydlarz LD, Jones LE, Harvell CD. Innate Immunity, Environmental Drivers, and Disease Ecology of Marine and Freshwater Invertebrates. Annu Rev Ecol Evol Syst. 2006;37:251–88.
Palmer CV, Traylor-Knowles N. Towards an integrated network of coral immune mechanisms. Proceedings of the Royal Society B: Biological Sciences. 2012;279:4106–14.
van de Water JAJM, Chaib De Mares M, Dixon GB, Raina J-B, Willis BL, Bourne DG, et al. Antimicrobial and stress responses to increased temperature and bacterial pathogen challenge in the holobiont of a reef-building coral. Mol Ecol. 2018;27:1065–80.
Zaragoza WJ, Krediet CJ, Meyer JL, Canas G, Ritchie KB, Teplitski M. Outcomes of Infections of Sea Anemone Aiptasia pallida with Vibrio spp. Pathogenic to Corals. Microb Ecol. 2014;68:388–96.
Alagely A, Krediet CJ, Ritchie KB, Teplitski M. Signaling-mediated cross-talk modulates swarming and biofilm formation in a coral pathogen Serratia marcescens. ISME J. 2011;5:1609–20.
Krediet CJ, Meyer JL, Gimbrone N, Yanong R, Berzins I, Alagely A, et al. Interactions between the tropical sea anemone Aiptasia pallida and Serratia marcescens, an opportunistic pathogen of corals. Environmental Microbiology Reports. 2014;6:287–92.
Roesel CL, Vollmer SV. Differential gene expression analysis of symbiotic and aposymbiotic Exaiptasia anemones under immune challenge with Vibrio coralliilyticus. Ecol Evol. 2019;9:8279–93.
Fuess LE, Pinzón C JH, Weil E, Grinshpon RD, Mydlarz LD. Life or death: disease-tolerant coral species activate autophagy following immune challenge. Proc R Soc B. 2017;284:20170771.
Kvennefors ECE, Leggat W, Kerr CC, Ainsworth TD, Hoegh-Guldberg O, Barnes AC. Analysis of evolutionarily conserved innate immune components in coral links immunity and symbiosis. Developmental & Comparative Immunology. 2010;34:1219–29.
Williams LM, Fuess LE, Brennan JJ, Mansfield KM, Salas-Rodriguez E, Welsh J, et al. A conserved Toll-like receptor-to-NF-κB signaling pathway in the endangered coral Orbicella faveolata. Developmental & Comparative Immunology. 2018;79:128–36.
Perez S, Weis VM. Nitric oxide and cnidarian bleaching: an eviction notice mediates breakdown of a symbiosis. Journal of Experimental Biology. 2006;209:2804–10.
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.
Brennan JJ, Messerschmidt JL, Williams LM, Matthews BJ, Reynoso M, Gilmore TD. Sea anemone model has a single Toll-like receptor that can function in pathogen detection, NF-κB signal transduction, and development. Proc Natl Acad Sci USA. 2017;114:E10122–31.
Palmer CV, McGinty ES, Cummings DJ, Smith SM, Bartels E, Mydlarz LD. Patterns of coral ecological immunology: variation in the responses of Caribbean corals to elevated temperature and a pathogen elicitor. Journal of Experimental Biology. 2011;214:4240–9.
Connelly MT, McRae CJ, Liu P-J, Traylor-Knowles N. Lipopolysaccharide treatment stimulates Pocillopora coral genotype-specific immune responses but does not alter coral-associated bacteria communities. Developmental & Comparative Immunology. 2020;109:103717.
Livny J, Zhou X, Mandlik A, Hubbard T, Davis BM, Waldor MK. Comparative RNA-Seq based dissection of the regulatory networks and environmental stimuli underlying Vibrio parahaemolyticus gene expression during infection. Nucleic Acids Research. 2014;42:12212–23.
Czachor JS. Unusual aspects of bacterial water-borne illnesses. American family physician. 1992;46:797–804.
Lipp EK, Rose JB. The role of seafood in foodborne diseases in the United States of America. Revue Scientifique et Technique-Office International des Epizooties. 1997;16:620–640.
Chimetto LA, Brocchi M, Thompson CC, Martins RCR, Ramos HR, Thompson FL. Vibrios dominate as culturable nitrogen-fixing bacteria of the Brazilian coral Mussismilia hispida. Systematic and Applied Microbiology. 2008;31:312–9.
Nithyanand P, Pandian SK. Phylogenetic characterization of culturable bacterial diversity associated with the mucus and tissue of the coral Acropora digitifera from the Gulf of Mannar: Culturable bacterial diversity of A. digitifera. FEMS Microbiology Ecology. 2009;69:384–94.
Zhang Q-L, Zhu Q-H, Liang M-Z, Wang F, Guo J, Deng X-Y, et al. Comparative transcriptomic analysis provides insights into antibacterial mechanisms of Branchiostoma belcheri under Vibrio parahaemolyticus infection. Fish & Shellfish Immunology. 2018;76:196–205.
Augustin R, Fraune S, Bosch TCG. How Hydra senses and destroys microbes. Seminars in Immunology. 2010;22:54–8.
Brown T, Rodriguez-Lanetty M. Defending against pathogens – immunological priming and its molecular basis in a sea anemone, cnidarian. Sci Rep. 2015;5:17425.
Duncan JA, Canna SW. The NLRC4 Inflammasome. Immunol Rev. 2018;281:115–23.
van der Burg CA, Prentis PJ, Surm JM, Pavasovic A. Insights into the innate immunome of actiniarians using a comparative genomic approach. BMC Genomics. 2016;17:850.
Ceccarelli D, Hasan NA, Huq A, Colwell RR. Distribution and dynamics of epidemic and pandemic Vibrio parahaemolyticus virulence factors. Front Cell Infect Microbiol. 2013;3. doi:10.3389/fcimb.2013.00097.
Watanabe Y, Tateno H, Nakamura-Tsuruta S, Kominami J, Hirabayashi J, Nakamura O, et al. The function of rhamnose-binding lectin in innate immunity by restricted binding to Gb3. Developmental & Comparative Immunology. 2009;33:187–97.
Kvennefors ECE, Leggat W, Hoegh-Guldberg O, Degnan BM, Barnes AC. An ancient and variable mannose-binding lectin from the coral Acropora millepora binds both pathogens and symbionts. Developmental & Comparative Immunology. 2008;32:1582–92.
Poole AZ, Weis VM. TIR-domain-containing protein repertoire of nine anthozoan species reveals coral–specific expansions and uncharacterized proteins. Developmental & Comparative Immunology. 2014;46:480–8.
Quistad SD, Traylor-Knowles N. Precambrian origins of the TNFR superfamily. Cell Death Discovery. 2016;2:16058.
Quistad SD, Stotland A, Barott KL, Smurthwaite CA, Hilton BJ, Grasis JA, et al. Evolution of TNF-induced apoptosis reveals 550 My of functional conservation. Proc Natl Acad Sci USA. 2014;111:9567–72.
Bergsbaken T, Fink SL, Cookson BT. Pyroptosis: host cell death and inflammation. Nat Rev Microbiol. 2009;7:99–109.
Lamkanfi M, Dixit VM. Manipulation of Host Cell Death Pathways during Microbial Infections. Cell Host & Microbe. 2010;8:44–54.
Zitvogel L, Kepp O, Kroemer G. Decoding Cell Death Signals in Inflammation and Immunity. Cell. 2010;140:798–804.
Elliott MR, Ravichandran KS. Clearance of apoptotic cells: implications in health and disease. The Journal of Cell Biology. 2010;189:1059–70.
Fuess LE, Pinzόn C. JH, Weil E, Mydlarz LD. Associations between transcriptional changes and protein phenotypes provide insights into immune regulation in corals. Developmental & Comparative Immunology. 2016;62:17–28.
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. Developmental & Comparative Immunology. 2009;33:559–69.
Mashaghi A, Marmalidou A, Tehrani M, Grace PM, Pothoulakis C, Dana R. Neuropeptide substance P and the immune response. Cell Mol Life Sci. 2016;73:4249–64.
Wang D. The essential role of G protein-coupled receptor (GPCR) signaling in regulating T cell immunity. Immunopharmacology and Immunotoxicology. 2018;40:187–92.
Lu E, Cyster JG. G‐protein coupled receptors and ligands that organize humoral immune responses. Immunol Rev. 2019;289:158–72.
Vidal-Dupiol J, Ladriere O, Meistertzheim A-L, Foure L, Adjeroud M, Mitta G. Physiological responses of the scleractinian coral Pocillopora damicornis to bacterial stress from Vibrio coralliilyticus. Journal of Experimental Biology. 2011;214:1533–45.
Vidal-Dupiol J, Ladrière O, Destoumieux-Garzón D, Sautière P-E, Meistertzheim A-L, Tambutté E, et al. Innate Immune Responses of a Scleractinian Coral to Vibriosis. J Biol Chem. 2011;286:22688–98.
Vidal-Dupiol J, Dheilly NM, Rondon R, Grunau C, Cosseau C, Smith KM, et al. Thermal Stress Triggers Broad Pocillopora damicornis Transcriptomic Remodeling, while Vibrio coralliilyticus Infection Induces a More Targeted Immuno-Suppression Response. PLoS ONE. 2014;9:e107672.
Bolger AM, Lohse M, Usadel B. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics. 2014;30:2114–20.
Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol. 2011;29:644–52.
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. 2011;:16.
Langmead B, Salzberg SL. Fast gapped-read alignment with Bowtie 2. Nat Methods. 2012;9:357–9.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26:139–40.
Conesa A, Nueda MJ, Ferrer A, Talon M. maSigPro: a method to identify significantly differential expression profiles in time-course microarray experiments. Bioinformatics. 2006;:1096–102.
Gotz S, Garcia-Gomez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, et al. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Research. 2008;36:3420–35.
Hu Z-L, Reecy JM. CateGOrizer: A Web-Based Program to Batch Analyze Gene On- tology Classification Categories. :3.

Download PDF

Journal Publication

published 09 Nov, 2020

Read the published version in BMC Genomics →

Editorial decision: Minor revision
02 Oct, 2020
Editor assigned by journal
16 Sep, 2020
Submission checks completed at journal
15 Sep, 2020
Editor invited by journal
15 Sep, 2020

You are reading this older preprint version

Read the latest preprint version →

Gene expression kinetics of Exaiptasia pallida innate immune response to Vibrio parahaemolyticus infection

Status:

Journal Publication

Version 2

Abstract

Figures

Background

Results

Discussion

Conclusion

Methods

Abbreviations

Declarations

References

Supplementary Files

Status:

Journal Publication

Version 2