The Significance of Genotypic Diversity in Coral Competitive Interaction: A Transcriptomic Perspective

Competitive interactions shape coral assemblages and govern the dynamics of coral ecosystems. Although competition is an ecological concept, the outcomes of competitive interactions are ultimately determined by patterns of gene expression. These patterns are subject to genotypic variation on both sides of any interaction. Such variation is typically treated as “noise”, but it is sometimes possible to identify patterns within it that reveal important hidden factors in an experiment. To incorporate genotypic variation into the investigation of coral competitive interactions, we used RNA-sequencing to study changes in gene expression in a hard coral (Porites cylindrica) resulting from non-contact competition experiment with a soft coral (Lobophytum pauciflorum). Hard coral genotype explained the largest proportion of variation between samples; however, it was also possible to detect gene expression changes in 76 transcripts resulting from interaction with the soft coral. In addition, we found a group of 20 short secreted proteins that were expressed as a coordinated unit in three interacting Porites-Lobophytum pairs. The presence of this secretion response was idiosyncratic in that it could not be predicted based on polyp behaviour, or the genotype of hard or soft coral alone. This study illustrates the significance of individual variation as a determinant of competitive behaviour, and also provides some intriguing glimpses into the molecular mechanisms employed by hard corals competing at a distance.


INTRODUCTION
Competition is an important ecological interaction, especially in highly productive tropical systems such as rainforests and coral reefs where it is a driver of ecosystem dynamics (Connell et al., 2004;Álvarez-Noriega et al., 2018). Competition also plays an important role in determining the impacts of climate change and other anthropogenic impacts on these systems. On coral-reefs, where hard corals (Scleractinia) have historically been dominant, many locations have seen shifts in species composition in favour of other reef taxa such as macroalgae (Roff and Mumby, 2012), octocorals (Lenz et al., 2015;Lasker et al., 2020), zoanthids (Cruz et al., 2016), and sponges (Bell et al., 2013). These shifts may themselves result from altered competition between reef taxa, and on non-Scleractinia dominated reefs the frequency of interactions between scleractinians and other major reef taxa is increased (Ladd et al., 2019). Even on reefs where Scleractinia still dominate, the role of competition, and the nature of competitive hierarchies is changing (Horwitz et al., 2017;Johnston et al., 2020), leading to shifts in the structure and function of reef communities.
The effects of competitive interactions between reef taxa are challenging to measure because the outcomes of competition play out slowly, are not strictly hierarchical (Precoda et al., 2017) and can reverse over time (Bak et al., 1982). In addition, most field surveys and experiments to date have relied on visible signs to detect competitive interactions and determine their order of dominance. Visible competitive strategies of corals include: overtopping to starve competitors of light; deployment of mesenteric filaments to externally digest a competitor; and elongation of polyps or development of sweeper tentacles to enable contact followed nematocyst discharge [reviewed by Lang and Chornesky, 1990;Chadwick and Morrow, 2011;Yosef et al., 2020). Although these physical signs are reliable indicators of competition when competitors are in contact, it is now clear that a wide range of reef taxa including scleractinian corals, octocorals, sponges, and algae Sammarco et al., 1983;Fearon and Cameron, 1996;Koh and Sweatman, 2000;Chadwick and Morrow, 2011) all produce toxins that could mediate competitive interactions without close contact. Despite abundant evidence of non-contact competitive capability in a variety of reef taxa, research on competitive strategies in corals has overwhelmingly focussed on interactions that involve contact (Chornesky, 1983;Sebens and Miles, 1988;Tanner, 1995;Fleury et al., 2004;Shearer et al., 2012). This bias could lead to underestimates of key competitive interactions, especially those for which non-contact competition is the primary mode.
Molecular techniques such as transcriptomics and metabolomics have the potential to resolve key gaps in our understanding of competition between reef taxa but have so far seen little use in coral competition research (except: Shearer et al., 2012Shearer et al., , 2014Quinn et al., 2016). Importantly, these techniques can directly measure molecules involved in both defensive and aggressive responses to competition and can therefore be employed to study non-contact interactions or interactions that do not generate clear physical effects. This is underscored by the results of Shearer et al. (2012) who studied molecular responses of Acropora millepora to four species of macroalgae and found the greatest change in gene expression in a competitive regime that showed the least physiological evidence of competitive impact. In addition, molecular analyses may reveal the mechanisms that underpin individual variation in competitive outcome that have been shown to exist between and within species. Interspecific competitive outcomes between reef taxa can be difficult to predict, with highly idiosyncratic dominance relationships between individual species pairs (Precoda et al., 2017). Dominance relationships may also depend on variations in genotype or physiological state of individual competitors. Although this has not been explicitly explored in the context of competition, evidence from molecular studies across a range of other extrinsic factors suggests that such intraspecific variation in response to stressors is likely to be high (Marshall and Baird, 2000;Loya et al., 2001;Obura, 2001;Fitt et al., 2009;Hughes et al., 2017;Sekizawa et al., 2017;Wright et al., 2017).
In this study we explore the transcriptomic response of Porites cylindrica, a hard coral, to competition with Lobophytum pauciflorum, a soft coral. The experiment was designed to investigate non-contact competition, although, as we show below, some limited contact via elongated mesenterial filaments was also observed. Since both coral types can be fragmented, we were able to pair each genotype of Porites with five genotypes of the competing soft coral. This design allows us to describe a core molecular response, which appears to be consistent across competing pairs, as well as a more specialised response involving up-regulation of secreted proteins that is restricted to a subset of competing pairs.

Competition Experiment
Molecular and behavioural responses to non-contact competition were investigated using an experiment conducted at Orpheus (Goolboddi) Island Research Station, in the central Great Barrier Reef,Australia (18'34 'S;146'29'E). Five colonies of the soft coral Lobophytum and 18 nubbins (∼3 cm) from each of three colonies of the hard coral Porites were collected with a bonecutter from reefs around Orpheus Island (GBRMPA Permit No. G16/38499.1). The Porites nubbins were fixed onto ceramic tiles. Each soft coral colony was cut into 12 pieces containing one or two lobes/fingers (∼5 cm) and these were placed on top of the tiles but not attached. The hard and soft coral fragments were then allowed to recover for three weeks prior to the start of the experiment.
After the acclimatization period, corals were placed in tanks (1,300 ml, open system, 400 ml/min of 10 µm filtered sea water) for 60-days. In each tank, a soft and a hard coral piece were placed purposely <3 cm apart to prevent corals touching each other, and simulate a non-contact interaction from the start of the experiment, while isolated hard and soft corals were used as controls ( Figure 1A). This pair-wise design was built with five biological replicates of the soft coral (colonies: La, Lb, Lc, Ld, and Le) and three biological replicates of Porites (colonies: Pd, Pe, and Pf). Combinations will be referred to by listing the hard coral followed by the soft coral (e.g., Pd-La) while controls are denoted with a C (e.g., Pd-C). In total, the experiment was composed of 15 biological combinations of interacting corals (e.g., Pd-La and Pe-La), and 3 hard coral controls (Pd-C, Pe-C, and Pf-C). For each combination and control there were two technical replicates/clones, resulting in a total of 36 experimental tanks.

Collection and Analysis of Porites Behavioural Data
Behavioural observations were recorded throughout the competition experiment to determine if Porites interacting with Lobophytum were showing signs of competitive behaviour or if their polyp activity was affected by the interaction. Porites polyp activity and competitive behaviour were observed three times per day to avoid bias due to highly variable diel polyp activity patterns (Levy et al., 2006). Observation times were between: 8 am-11 am, 12 pm-4 pm, and 6 pm-9 pm.  Supplementary Table 1. Polyp activity measurements were taken starting from day eight of the experiment and continuing until day 60.
In addition to basic polyp activity, competitive behaviour of Porites towards Lobophytum, such as elongated polyps (Sammarco et al., 1985;Rinkevich and Sakamaki, 2001) and/or mesenteric filaments, was also recorded.
These data were analysed using a cumulative link mixed effect model (clmm) with the package "ordinal" (Christensen, 2015) in the statistics program R (R Core Team, 2016), to determine if competition affected Porites polyp activity. A range of models were explored within this framework and the most parsimonious was selected on the basis of AIC 1 . The final model (Eq. 1) included polyp activity (Activity) as an ordered factor (Closed < Partially-open < Open) dependent on the following fixed effects; time categorized in eight groups of ∼5 days each (Time), the Porites colony the nubbin came from (Colony) and the nubbin's treatment (competition or control) (Treatment). In addition, tank was modelled as a random effect (Tank).

RNA Sequencing, Assembly and Transcript Quantification
Porites nubbins from one of the two replicates in each experimental condition were randomly sampled for RNA sequencing after 30 days of interaction with Lobophytum to determine the effects of competition on Porites gene expression. Samples were taken by quickly crushing the nubbin with a hammer and immediately, snap-freezing it in liquid nitrogen. Samples were stored at −80 • C until required. RNA was extracted with TRIzol Reagent (Ambion, catalogue#15596-026). RNA quality checks and library preparation were performed as described in supplementary methods. High-quality RNA extractions were obtained for nubbins from colonies of Porites Pd and Pf. It was not possible to extract RNA from nubbins of colony Pe, therefore 12 samples (10 nubbins in competition and two nubbins in control) from colonies Pd and Pf were sequenced. Samples were sequenced by the Australian Genome Research Facility (AGRF: Melbourne, Australia) using an Illumina HiSeq2500, to obtain approximately 14.5 million reads (100 bp paired-end) per sample. Reads were checked for quality, adapter content and other sequencing artefacts using FastQC (version 0.11.9). All sequencing data have been deposited with Genbank under bioproject (PRJNA706467).
A de novo transcriptome assembly for Porites was constructed by adapting evidence-based best practices (MacManes, 2016) to deal with data from two distinct genotypes (Pd and Pf). Data for each genotype were processed separately for read error correction with Rcorrector (Song and Florea, 2015), followed by initial assembly with Trinity (version2.4.0; Grabherr et al., 2011) using options to enable read trimming and normalization. Independent assemblies produced by this process were then merged together using the software TransFuse (version0.5.0 2 ) with a 95% identity threshold for merging clusters. The merged transcriptome was analysed with the software TransRate (version1.0.3; Smith-Unna et al., 2016), which scores contigs based on agreement with mapped raw reads. All high quality (called "good" by TransRate) contigs retained from this process were subjected to analysis with software Psytrans 3 to remove those likely to originate from algal symbionts (family Symbiodiniaceae). The completeness of the clean assembly was assessed with the software BUSCO version4 (Simão et al., 2015).
To assess the effectiveness of Psytrans, and to check for additional contaminating organisms, we performed two additional analyses on the transcriptome data. Firstly, diamond blastx (version 2.0.7.145; Buchfink et al., 2015) was used to identify the best match and its corresponding phylum of origin in the NCBI nr database (E-value < 0.01 in very-sensitive mode) for all transcripts remaining after processing with Psytrans. Secondly, the lowest common ancestor was inferred for all reads using kraken (version 1; Wood and Salzberg, 2014) for each sample. See supplementary methods for detailed information on construction of the kraken database.
Corrected reads (from Rcorrector) were first trimmed using Trimmomatic (version 0.36) and then mapped to the final transcriptome assembly using Bowtie2 with recommended settings (end to end alignments, report all alignments, min alignment score 0.3) to suit downstream quantification. Corset (version1.05; Davidson and Oshlack, 2014) was then used to cluster transcripts likely to have originated from the same gene and count reads assigned to clusters.
Results from Corset were used to identify a set of 17,093 transcripts suitable for genetic analysis. These transcripts were expressed with at least three reads in all samples and belonged to a singleton cluster (likely a single gene copy with no alternative splicing variants). Analysis of mapped reads for these contigs with ANGSD (Korneliussen et al., 2014) revealed 30,289 SNPs which were then used to calculate the relatedness statistic, theta between each pair of samples with ngsRelate (Hanghøj et al., 2019). This analysis confirmed that the two colonies used for sequencing (Pd and Pf) were unrelated to each other (distinct genotypes; theta = 0) while individual fragments from a single colony were clones (theta = 0.5).
Preliminary differential expression analysis with DESeq2 (Love et al., 2014) revealed a small set of clusters that were exclusively expressed in three samples and that appeared to be of barnacle origin (>90% sequence similarity to barnacle transcripts, via blastx to the NCBI non redundant protein database). In order to identify and remove all barnacle clusters, minimap2 (Li, 2018) was used to map all transcripts to a genomic database consisting of the Porites lutea (Robbins et al., 2019), Cladocopium goreaui (Liu et al., 2018), and Amphibalanus amphitrite genomes (Kim et al., 2019). Mapping was performed using the xsplice option to allow gapped alignment and up to 10% sequence divergence. Any transcript that produced a valid alignment to the barnacle (Amphibalanus amphitrite) genome was deemed to be of barnacle origin. This resulted in identification of 2,297 transcripts and 492 clusters likely to be of barnacle origin, the expression of which was almost entirely restricted to three samples (Pf-Lc, Pd-Le, and Pd-La) (Supplementary Figure 2). Although all Porites nubbins were subject to frequent visual inspection throughout the experiment, it is likely that these barnacle reads originated from commensal, coral associated barnacles (Tsang et al., 2014), which tolerate overgrowth by the coral skeleton and coexist with the coral (Liu et al., 2016). Our analysis assumes that this relationship (likely commensal) had minimal effect on gene expression in the affected nubbins.
After removal of all barnacle clusters the gene expression profiles of barnacle-affected samples clustered together with other samples suggesting that the effect of barnacles on host gene expression was minimal. We therefore retained barnacle-affected samples for further analysis.
Statistical analysis of cleaned gene expression data was performed using DESeq2 with the goal of detecting transcript clusters consistently differentially expressed between control and treatment nubbins. To do this, a single factor capturing all experimental sample groupings (Porites genotype, Lobophytum genotype, and Control) was fitted and hypothesis testing was performed based on a contrast between all conditions involving competitive interactions (non-control samples) and controls. Clusters were deemed to be differentially expressed under competition if they had an adjusted p-value (padj) < 0.1 under this contrast. Complete details and code used to perform this analysis are provided as an online repository 4 .
Additional manual annotation was performed to supplement results for 76 transcript clusters found to be differentially expressed under competition. For each of these transcripts InterProScan (version5.48-83.0; Jones et al., 2014) was first used to identify conserved domains. If the transcript had a blast hit to a protein from SwissProt, the domain structure of this hit was compared with the domain structure of the transcript to determine whether genuine functional homology could be inferred. If domain structure was not conserved, an attempt was made to infer function based on conserved domains alone. Cnidarian-specific functional information was then identified by using google scholar to search for papers containing the combination of conserved domain names, gene names and the words coral or cnidarian. Manual curation of these search results yielded a small number of gene expression studies with relevant functional information. Finally, an attempt was made to assign the transcript to one or more of the following functional categories: immune response, stress response, secreted proteins and toxin. Supplementary Table 2 lists evidence of homology and inferred functions for these 76 differentially expressed transcripts.

Polyp Activity and Behaviour
Analysis of the polyp activity data (Figure 1C and Supplementary Table 3) showed that polyps were more likely to be closed or partially closed in colonies under competition compared with controls (p < 0.01; based on a clmm) and differed significantly between Porites genotypes (p < 0.001). Polyp activity also changed over time but this seemed to vary between Porites genotypes, with Pd and Pf showing a gradual increase in open polyps during the first 30 days of the experiment, whereas genotype Pe showed little change.
In addition to basic polyp activity, competitive behaviour in the form of mesenteric filament formation was observed for six of the ten competing Porites nubbins, including representatives of genotypes Pd (four nubbins) and Pf (two nubbins) but not for Pe (Supplementary Table 4). Although filaments were clearly visible when present (Figure 1B) they were short-lived and it is likely that some instances of filament formation were not observed.

Transcriptome Assembly
Assembly of the Porites cylindrica transcriptome with Trinity resulted in 532,484 and 502,263 raw transcripts for Pf and Pd genotypes, respectively. Merging these assemblies with transfuse resulted in a total of 709,417 contigs and of these 422,222 were found to be good contigs by transrate (see section "Methods"). Splitting this assembly using Psytrans resulted in a coral fraction (340,399 contigs) and a Symbiodiniaceae fraction (81,823 contigs). The average mapping rate of the raw corrected reads to the combined coral-Symbiodiniaceae transcriptome was 83.6% while the mapping rate to the coral-only fraction was 53.4% due to a high proportion of reads being of Symbiodiniaceae origin (Supplementary Table 5). The coral-only assembly had a longest contig of 43 kb and average contig length of 992 bp. It contained complete copies of 95.7% of metazoan BUSCO genes (35.2% single copy; 60.5% duplicated) as well as 2.2% fragmented BUSCOs. This percentage of completeness is similar to that observed in other de novo coral transcriptome assemblies from adult tissue, such as Acropora gemmifera (94.4%; Oldach and Vize, 2018).
Assessment of the taxonomic composition of assembled transcripts with blastx was challenging because only 155,814 of the 340,399 transcripts (46%) in our transcriptome could be classified via homology, probably due to poor representation of cnidarian sequences in the nr database. Of the transcripts classified to phylum level by blastx, 86% were of cnidarian origin and tended to have better matches (higher bitscores) than the remaining transcripts matching other phyla (Supplementary  Figure 4). The two most abundant phyla in this longtail of non-cnidarian sequences were Arthropoda (6%) and Chlorophyta (4%) possibly reflecting the presence of barnacle and Symbiodiniaceae transcripts. Analysis with kraken recovered the known pattern of barnacle contamination in three samples, but otherwise suggested that taxonomic composition was consistent (Supplementary Figure 3) and not confounded with key groupings identified in our differential expression analysis (see below).

Effect of Genotype and Competition on Gene Expression
Gene expression data obtained at 30 days after the onset of competition were dominated by differences between the two sequenced Porites genotypes (Pd and Pf). Nubbins from these two colonies separated clearly into two groups along the first principal component of a PCA, and accounted for 81% of total variance in expression for the top 500 most variable genes (Figure 2A). Variation in gene expression due to competition with Lobophytum appeared to be associated with the second principal component of this PCA but could not simply be attributed to Porites or Lobophytum genotype alone. Sample clustering based on 76 genes differentially expressed between control nubbins and those in competition ( Figure 2B) revealed four groups of samples, the composition of which was linked to position along PC2. In particular, sample group three (dark-blue) (Pf-Ld, Pf-La, and Pd-Lb) included all those at one extreme of PC2, while control samples (sample group one; grey) occupied a central position and the remaining samples (groups 4 and 2; red and light-blue, respectively) occupied the other extreme of PC2 ( Figure 2B).
Row clustering also revealed key groups of genes with expression profiles that differ between these sample groups ( Figure 2B). The most striking example of this is demonstrated by row cluster four, which included genes that were exclusively expressed under competition in all samples from column group three (dark-blue) and not at all in other samples. Genes in row clusters one and two had less consistent expression within sample groups but generally partitioned genes into those that were overexpressed (row cluster one) or under expressed (row cluster two) in competition versus controls. Finally, row cluster five captured six genes which differed strongly between samples from genotype Pd and Pf, and where the response to competition was in the opposite direction depending on Porites genotype ( Figure 2B).

Genes Differentially Expressed in Response to Non-contact Competition
A total of 76 transcripts were found to be differentially expressed in response to non-contact competition (DESeq2 adj p-value < 0.1). Of these, it was possible to obtain functional annotations for 30 by combining information from homology to proteins with known function, the presence of conserved domains, and gene expression or population genetic studies in cnidarians (see Supplementary Table 2 for a full list including a summary of evidence). Very few could be assigned meaningful gene ontology terms because although some (24) had highly significant blast matches to SwissProt proteins, analysis with InterProScan revealed differences in domain structure that cast doubt on conservation of function. Nevertheless, it was possible to manually assign 35 functionally annotated genes as having putative roles in immunity (10) or cellular stress (9) including response to ROS and unfolded proteins. In addition, 20 genes were classified as secreted proteins and four of these had hallmarks of toxins (short, secreted and with ShKT, CRISP or protease domains).
Examination of correspondence between row groups in Figure 2B and functional categories revealed a striking enrichment [14/23 (60%) transcripts] in numbers of secreted proteins in row group four compared with differentially expressed genes as a whole [20/76 (20%) transcripts]. None of the proteins in this group were close homologs to each other (no clusters at 70% similarity with cd-hit) but many had the hallmarks of toxins, including two with ShKT domains, one secreted peptidase and five others that were secreted and short (<200 AA).
Transcripts related to stress response and immunity had more complex and variable expression profiles than the secreted group (row cluster 4) and were therefore not associated with specific row clusters in Figure 2B. When considering the distribution of transcripts according to their effect size rank (Figure 3), transcripts putatively assigned to immunity were more often up-regulated (8/10 genes) in response to competition while those associated with the core cellular stress response included a more equal mix of genes up and down regulated (three up-regulated; five down-regulated). The strong tendency for secreted proteins to be up-regulated under competition reflects enrichment for those in row group 4 ( Figure 2B).

DISCUSSION
The dominant effect of colony as a driver of gene expression profiles in this study agrees with, and adds to, a growing body of evidence that this can be a major source of transcriptional variability in corals (Seneca et al., 2010;Granados-Cifuentes et al., 2013;Parkinson et al., 2018). While this overall inter-colony variability can be broadly viewed as a source of noise, it also has important ecological implications, especially in the context of competition where it could underpin the apparently idiosyncratic responses seen in past experiments (Precoda et al., 2017). Here we show that, at the transcriptional level, this variability is not random, but governed by latent factors that divide samples into groups with distinct expression profiles. Even with the relatively small number of samples analysed here it was possible to identify a distinct subset defined by exclusive expression of a suite of secreted proteins (Figure 2).
This association between a distinct pattern of expression and functional category (secretion) suggests that the genes encoding these proteins are part of a coordinated response, perhaps via a common transcriptional regulatory mechanism. Although this response was not correlated with any behavioural factors such as mesenterial filaments or polyp activity, the sequences of genes involved provide some indication of its functional significance. One possibility, supported by the presence of putative toxins among these secreted proteins, is that they are involved in aggression towards the competitor. If this is the case, the mechanism of target delivery for putative toxins remains unclear as the only means of contact between competitors was via mesenterial filaments. These filaments were observed sporadically, and their presence was not correlated with samples showing a secretion response, which suggests that this mode of delivery is unlikely, although it cannot be ruled out. Another working hypothesis is that these candidate toxins could reach the competitor through tentacle contact. Tentacle attack is usually seen through the development of sweeper tentacles where the tip is enriched in nematocyst and other toxins, as documented for Galaxea (Yosef et al., 2020) but this strategy has neither been registered for Porites nor was it observed in this experiment. Another possibility is that these proteins could be secreted into the surrounded Positive ranks represent transcripts up-regulated in competition compared to controls. Transcripts belonging to multiple categories are shown as multi-coloured bars. Blue points represent putative toxins. Category "Stress" includes transcripts associated with response to reactive oxygen species as well as those involved in the unfolded protein response. Labelled points correspond to transcripts with strong evidence of homology to named proteins. CAT, catalase; AOSL, allene oxide synthase-lipoxygenase; SYT, synaptotagmin-like; TNFR, tumor necrosis factor receptor; CALR, calreticulin. water and reach a competitor via diffusion through the water column. Although there is evidence to suggest that a wide range of cnidarians, including Porites cylindrica, may use allelopathy a form of non-contact competition-to inhibit growth  and alter settlement of other corals (Maida et al., 1995;Da-Anoy et al., 2017), the toxic molecules involved have either not been identified, or are small organic compounds Aceret et al., 1995) rather than proteins or peptides. Although peptides are metabolically expensive to produce (and therefore less likely to be released into the surrounding water) they can be highly specific and potent, and are known to be used by cyanobacteria as allelochemicals (Gross, 2003). Our results indicate these short-secreted proteins could be part of a peptide-mediated allelopathic response in corals and warrant further investigation to clarify their roles in competition and determine their bioactivity and mode of delivery.
In addition to the secretion response observed for a subset of Porites/Lobophytum pairs, it was also possible to identify genes with consistent responses to non-contact competition across all Porites samples. This "core" response included genes with putative roles in immunity and with roles in responding to cellular stress. The closest comparable experiments are those of Shearer et al. (2012Shearer et al. ( , 2014 who explored gene expression responses in acroporid corals (Acropora millepora and Montipora digitata) to acute (6 and 48 h) and longer-term (20 days) competition with a range of macroalgal species. Short-term exposure experiments highlighted differential expression of genes across a wide range of putative roles, including immunity and cellular stress, but found that responses were highly species specific. Long-term exposure resulted in differential expression of a suite of genes involved in cellular stress, however, the specific genes involved were almost entirely different from those found differentially expressed by Porites nubbins in this experiment. Specifically, they included up-regulation of six heat shock proteins which we did not observe, and a change in expression of calreticulin in the opposite direction to that seen here.
Traditionally Porites spp. have been regarded as weak or non-aggressive competitors (Sheppard, 1979;Dizon and Yap, 2005), but both the morphological and molecular data presented here challenge this idea. Mesenteric attack ( Figure 1B) is a form of aggression employed by some corals in competitive interactions, but not previously documented for any Porites spp. The molecular data implying secretion of candidate toxins by Porites under challenge are also inconsistent with the response being passive, and investigation of the properties of these novel proteins should be a priority. Although genotype is typically the strongest influence in transcriptomic experiments on corals, the work presented here illustrates the complexity and heterogeneity of responses in competitive interactions, and the extent to which these cannot be accounted for by genotype alone.

DATA AVAILABILITY STATEMENT
The transcriptomic sequencing data presented in this study are deposited in Genbank under bioproject PRJNA706467.

AUTHOR CONTRIBUTIONS
NA, AM, and DM designed the study. NA and AM performed sampling and process samples. IC and NA performed transcriptomic data analysis and the manuscript writing. RJ, IC, and NA performed the statistical analysis. DM, IC, NA, and RJ edited the draft manuscript. All authors read and approved the final version of the manuscript.

FUNDING
This study was funded through the Ecuador National Secretariat of Higher Education, Science, Technology, and Innovation (SENESCYT) doctoral scholarship; College of Public Health, Medical and Veterinary Sciences-SSA funds; and ARC Centre of Excellence for Coral Reef Studies.