Environmental and molecular regulation of asexual reproduction in the sea anemone Nematostella vectensis

Cnidarians exhibit incredible reproductive diversity, with most capable of sexual and asexual reproduction. Here, we investigate factors that influence asexual reproduction in the burrowing sea anemone Nematostella vectensis, which can propagate asexually by transverse fission of the body column. By altering culture conditions, we demonstrate that the presence of a burrowing substrate strongly promotes transverse fission. In addition, we show that animal size does not affect fission rates, and that the plane of fission is fixed along the oral–aboral axis of the polyp. Homeobox transcription factors and components of the TGFβ, Notch, and FGF signalling pathways are differentially expressed in polyps undergoing physal pinching suggesting they are important regulators of transverse fission. Gene ontology analyses further suggest that during transverse fission the cell cycle is suppressed, and that cell adhesion and patterning mechanisms are downregulated to promote separation of the body column. Finally, we demonstrate that the rate of asexual reproduction is sensitive to population density. Collectively, these experiments provide a foundation for mechanistic studies of asexual reproduction in Nematostella, with implications for understanding the reproductive and regenerative biology of other cnidarian species.

arguably makes investigation of their reproductive strategies important for our understanding of ocean ecology and conservation. We are interested in understanding the mechanisms that regulate asexual reproduction in this diverse group. In different organisms, asexual reproduction can be obligate or facultative. Although cnidarians exhibit diverse reproductive strategies, most species are capable of facultative asexual reproduction [1]. Determining how sexual and asexual reproduction are regulated is an important first step to begin understanding the mechanisms that dictate an organism's reproductive health and how it allocates resources between alternative reproductive strategies. While sexual reproduction has been relatively well studied in cnidarians, particularly coral spawning, causal factors that affect rates of asexual reproduction are not yet well characterized in any species [2], and most studies are descriptive-focusing on how the process of asexual reproduction occurs (e.g. budding, fissioning) [3][4][5].
Environmental factors have been found to influence asexual reproduction in non-cnidarian species including stochastic environmental conditions [6], mate availability [7], seasonality [8] and population density [7,9]. Individual attributes such as physical condition [10] and body size [11,12] have also been shown to have an influence on asexual reproduction. Within cnidarians, seasonality can affect asexual reproduction in Aurelia [13], population density and body size have been shown to regulate the rates of bud formation in Hydra [14,15], and feeding has been demonstrated to increase asexual reproduction in Hydra, Aurelia and Nematostella [4,13,14]. Efforts to identify molecular mechanisms that regulate asexual reproduction are limited to the window during which a clone is being generated. Within Bilateria, gene expression has been characterized during asexual fissioning in the annelid Pristina leidyi [16], and functional experiments argue TGFβ, Wnt signalling, and homeobox transcription factors regulate fissioning in Planaria [12,17,18]. Within Cnidaria, FGF signalling has been implicated in extracellular matrix degradation leading to clone detachment during asexual reproduction in the coral Pocillopora acuta [19], and both FGF and Notch regulate bud site position, outgrowth and detachment in Hydra [20][21][22]. Despite the current efforts, there is still a lack of data about the regulators of asexual reproduction, especially within the cnidarians, arguing additional efforts in currently studied species and additional species will facilitate uncovering broad phylogenetic patterns and individual specific regulators of asexual reproduction [2].
Here, we provide an initial characterization of factors that influence asexual reproduction at both the individual and population levels in the cnidarian sea anemone Nematostella vectensis. Nematostella undergoes asexual reproduction by transverse fission, where physal pinching leads to the separation of the lower body column ( physa) which regenerates into a clonal individual [23,24] (figure 1a). Additionally, Nematostella readily spawn in laboratory culture, making it well suited for future studies aimed at investigating condition-dependent aspects of cnidarian facultative asexual reproduction. To the best of our knowledge, this is the first study to characterize factors that affect facultative asexual reproduction in an anthozoan cnidarian. We investigated how environmental, physical and social conditions of Nematostella influenced transverse fission and identified genes differentially expressed in steady-state and actively fissioning animals. Lastly, we compared gene expression in high-and lowdensity populations which we demonstrate to have distinct rates of asexual reproduction. We found no changes in gene expression based on population density, arguing that differential gene expression does not play a large role in regulating the rate of asexual reproduction by transverse fission.

Substrate promotes asexual reproduction in Nematostella
Asexual reproduction is rarely observed in laboratory populations of Nematostella, but it is well documented in the wild [25]. Nematostella is an infaunal species, but laboratory populations are typically maintained without substrate to burrow into, suggesting that substrate may promote transverse fissioning in Nematostella. To test this, we compared the number of clonal progeny produced in cultures of 10 females with and without gravel substrate (n = 5 replicates). After two months, only two clonal progeny were observed within the five populations without substrate. By contrast, all five replicates with substrate increased by 25 or more clonal individuals (figure 1b; t-test: t 8 = −12.08, p < 0.0001). We conclude that substrate strongly promotes fission behaviour and we therefore included substrate in all subsequent asexual reproduction experiments.
To gain preliminary insight into whether there may be differences in transverse fission with and without substrate, we attempted to capture and image fissioning animals. Two time-lapse videos were created to show royalsocietypublishing.org/journal/rsos R. Soc. Open Sci. 10: 230152 examples of physal pinching via transverse fission in an individual in substrate and an individual removed from substrate once the process began (electronic supplementary material, videos S1 and S2). The animal kept in substrate appears to rotate the body column in a twisting motion at the site of pinching until the physal fragment detaches (electronic supplementary material, video S1). The animal removed from substrate did not exhibit the same twisting motion (electronic supplementary material, video S2), indicating that the presence of substrate may provide a mechanical advantage for fissioning animals.

Animal size does not affect fission rates
One previous study observed a highly variable latency to fission for Nematostella [4], but this account did not consider animal size as a potential explanatory variable. To address this possibility, a NvLWamide-like:: mCherry neural reporter line was used to isolate animals of different sizes based on the number of longitudinal neurons each animal possessed [26,27]. Twenty-nine per cent of animals observed during a three-week window underwent transverse fission. Binary logistic regression showed that starting animal size did not predict whether an animal asexually reproduced (figure 2a; r 2 = 0.005, animal size b = −0.003 ± 0.008, p = 0.71). The latency of animals that did fission ranged from 1 to 18 days, and was not related to animal size (figure 2b; Pearson correlation: r (20) = −0.06, p = 0.80). Some individuals fissioned multiple times over the three-week period. The decrease in size after fissioning, and repeated observations of the same animal showed no consistent pattern linking animal size and latency to fission (figure 2b, compare points of the same colour). Lastly, animals fissioned one to four times over a threeweek period, and there was no correlation between animal size and rate of asexual reproduction (figure 2c; Pearson correlation: r (11) = 0.26, p = 0.40). Taken together, these results demonstrate that size does not determine whether an animal will fission, how fast, or how often in Nematostella.

Physal pinching occurs at a set point along the body column
To determine if the physal pinching that precedes transverse fission occurs in a fixed location along the oral-aboral axis, we exploited the even distribution along the oral-aboral axis of longitudinal neurons expressing the NvLWamide-like::mcherry transgene [28]. The even distribution of neurons allowed us to royalsocietypublishing.org/journal/rsos R. Soc. Open Sci. 10: 230152 quantify the number of neurons present within the oral ( parent) and aboral (clone) remnants within less than 24 h of a fission event to determine the pre-fission size of the parent animals, and to express the position of the fission site as a percentage of the total body size. The difference between parent size before and after transverse fission increased with increasing adult size, indicating that clone size scaled with parent size (figure 3a). However, regardless of starting size, individuals lost 26.7% of their longitudinal neurons to their clones and there was relatively low variability between individuals in this value (±1.3% s.e.m.) (figure 3b). Individuals that underwent consecutive fissioning events showed the same degree of consistency in terms of the position of the fission event along the oral-aboral axis (figure 3a; e.g. compare aqua lines from animal 2). These data suggest that the position of physal pinching during transverse fission is fixed along the oral-aboral axis of the polyp.

Differential gene expression provides insight about mechanisms of transverse fission
Bulk RNA sequencing on fissioning ( pinching) and steady-state (not pinching) animals was used to identify genes differentially expressed during transverse fission. An average of 31 million reads were generated per sample which were mapped to the NVEC200 Nematostella genome [29]. A principal component analysis (PCA) revealed that samples segregated by pinching status along PC1, which explained 53% of the variance (figure 4a). A total of 1577 differentially expressed genes (DEGs) were identified with 517 genes upregulated and 1060 downregulated in pinching animals (figure 4b).
Within the list of DEGs, we identified some transcripts known to regulate fissioning or asexual reproduction in other species (see Introduction). For example, hox genes have been implicated in transverse fission in planarians and we identified five predicted homeobox transcripts that were upregulated in pinching animals. Upregulation of TGFβ pathway components is needed for fission initiation in Planaria [12], but we identified a putative TGFβ target (TGFβ-induced protein ig-h3-like) that was downregulated in our pinching animals. Similarly, elements of the notch pathway are upregulated during Hydra budding [21] but the five notch-like protein transcripts we found were all downregulated. Some FGF pathway components are upregulated during asexual reproduction in other cnidarians [19,20,22], and we identified two FGF-like ligands that were upregulated in pinching animals. However, an FGF-like ligand and multiple FGF-like receptors were downregulated, making the role of FGF in Nematostella fissioning less clear. Additionally, because transverse fission leads to the detachment of the physa from the rest of the body column, we predicted that transcripts related to cell adhesion and extracellular matrix organization would be differentially expressed. As predicted, transcripts that may promote cell adhesion were downregulated, and transcripts that likely regulate tissue organization and cytoskeletal architecture were differentially expressed. See electronic supplementary material, table S1, royalsocietypublishing.org/journal/rsos R. Soc. Open Sci. 10: 230152 for transcript IDs and annotations. Together these data imply that downregulation of transcripts related to maintaining body integrity likely aid in the separation of the physal fragment and suggest roles for major developmental signalling pathways in the Nematostella transverse fission process.

Gene ontology analysis uncovered functional characteristics of transverse fission behaviour
To provide a broader understanding of how fissioning is regulated, a rank-based gene ontology (GO) enrichment analysis of the DEGs found between pinching and non-pinching animals was conducted.   royalsocietypublishing.org/journal/rsos R. Soc. Open Sci. 10: 230152 the cell cycle becomes tightly regulated, most likely suppressed, while cell adhesion and patterning mechanisms are downregulated to promote effective separation of the physa from the body column. Additional studies are needed to fully understand how the genes identified in this study regulate transverse fissioning.

Asexual reproduction increases at low population density
An overarching hypothesis in species with facultative asexual reproduction is that there should exist specific environmental conditions that act to promote asexual reproduction, otherwise it is difficult to conceive of how multiple modes of reproduction could be evolutionarily maintained [6,30]. In the wild, tidal events can collect and trap Nematostella in isolated salt marsh pools at varying population densities. We hypothesized that population density would affect rates of asexual reproduction and predicted that asexual reproduction would increase in animals kept at low density due to reduced mating opportunity. Animals were placed into either a high-or low-density population (40 or 4 animals, respectively). After three weeks, the number of clones produced per starting individuals was determined for each density treatment. Although there was variation in the overall number of clones produced, with the least asexual reproduction occurring April-May, individuals in low-density treatments always exhibited greater asexual reproduction than those in high-density treatments regardless of time of year (figure 6a). On average, individuals from lowdensity populations generated significantly more clones than those from high-density populations (figure 6a'; t-test: t 32 = 3.2, p = 0.003).
We questioned if differences in metabolic waste build-up might explain the density-based differences in asexual reproduction. Hydrogen peroxide (H 2 O 2 ) is a natural by-product of oxidative metabolism in aerobic organisms which can act as a signalling molecule or cause oxidative damage to biomolecules depending on several factors including local concentration [31,32]. Nematostella are exposed to varying levels of H 2 O 2 from cellular respiration and from the tendency of peroxides to accumulate in their natural salt marsh habitats [33]. To test if observed differences in fission rates were due to differences in H 2 O 2 , we repeated the above experiment but added additional low-and high-density treatments that received H 2 O 2 in their water (final concentration 0.00025% or 82 µM). This concentration is not lethal and does not inhibit regeneration in Nematostella [33]. Like before, there was variability in the amount of asexual reproduction that occurred, but in nearly all replicates low-density animals produced more clones than high-density ones, regardless of whether H 2 O 2 was present (figure 6b). On average, there was an overall effect of density. Both low-density treatments produced significantly more clones than both high-density treatments (figure 6b'; mixed ANOVA-main effect density: F 1,24 = 4.1, p = 0.05) and this was not influenced by the addition of H 2 O 2 (mixed ANOVA-H 2 O 2 X density: F 1,24 = 0.007, p = 0.94). The amount of ammonia present across all treatments was also measured to determine if differences in this nitrogenous waste might explain density-based differences. Ammonia levels were low and did not differ between any treatments (figure 6c; mixed ANOVA-H 2 O 2 X density: F 1,24 = 0.06, p = 0.81). We conclude that Nematostella can sense their population density and then adjust their asexual reproduction accordingly. The mechanism by which this occurs is unknown, but differences in the metabolic waste products ammonia and H 2 O 2 are not responsible. royalsocietypublishing.org/journal/rsos R. Soc. Open Sci. 10: 230152 2.7. Gene expression did not vary based on population density To better understand the effect that density has on Nematostella, we performed bulk mRNA sequencing on whole individuals collected from low-and high-density treatments at the end of the experiment period. None of the sampled animals were undergoing physal pinching at the time of collection. We predicted that there would be differential gene expression based on density treatment alone, but that DEGs might not be related to the process of asexual reproduction since this process is unpredictable and asynchronous across individuals. Sequencing yielded an average of 50 million raw reads per sample which were mapped to the Nematostella NVE transcriptome v. 2.0 [34]. Using transcripts per million normalized read counts, the average distance between samples was estimated using a PCA that showed that samples did not separate along PC1 (60% variance) or PC2 (11% variance) based on population density treatment ( figure 7a). Further, a differential gene expression analysis showed that no transcripts were differentially expressed between samples from high-and low-density treatments (figure 7b), indicating that within our experimental context changes in gene expression do not explain density-based differences in asexual reproduction.

Discussion
One of our primary goals was to expand our knowledge of cnidarian reproduction by characterizing affecters of asexual reproduction in the cnidarian Nematostella. It was not surprising that substrate royalsocietypublishing.org/journal/rsos R. Soc. Open Sci. 10: 230152 enhanced and is effectively required for transverse fission to occur since they are naturally burrowing sea anemones. Our preliminary observations that animals that fission in substrate appear to twist, make it enticing to speculate that substrate offers a mechanical advantage to promote fissioning success. However, it is not clear if the attempts at fissioning are similar in substrate compared to no substrate cultures and substrate could also play a larger role by increasing the rate of physal pinching that leads to transverse fission. Regardless of the mechanism, it is important to note that our findings that substrate and low population density positively influence fission rates are immediately applicable for the quick expansion of laboratory cultures of Nematostella. Regarding animal size, we did not detect any size requirement related to transverse fission behaviour. Since when tracking individual animals there was variability in whether transverse fission occurred, there must be other factors that influence whether or not animals in isolation fission. We speculate that body condition is a likely factor because increases and decreases in feeding are known to influence fission rates in Nematostella [4]. Future experiments might consider controlling body condition, in part by isolating individuals and monitoring their food intake in advance of starting experiments. The rate of Nematostella asexual reproduction by transverse fission increases in populations at lower density. This observation argues that Nematostella are capable of sensing the number of nearby conspecifics and increase their rates of asexual reproduction to increase fitness when competition from other individuals is low or mating potential is decreased. This finding presents a number of obvious questions. For example, how do Nematostella sense population density and subsequently alter fission rates, or do sex ratios influence asexual reproduction rates at different population densities? Some insights come from our observations. There was no difference in gene expression between animals at high and low density which argues that density does not alter transcription rates, suggesting an alternative mechanism for density-dependent asexual reproduction rates. However, since none of the animals we collected were in the process of fissioning, we cannot exclude the possibility that densitydependent physal pinching is initiated locally by a minor change in gene expression that would be difficult to detect. Although, given that we detected 1577 DEGs between fissioning and steady-state individuals we believe that this explanation is less likely. Another possibility is that Nematostella use sensory information to interpret their social environment and in turn regulate their fission behaviour. Nematostella lack vision, and it is unlikely that density is determined by touch since burrowed adults rarely move around and individuals were placed equidistant from one another in our experiments. Our favoured hypothesis is that detection of a chemical signal in the environment initiates a cascade of events that mediate physal pinching and subsequently transverse fission. This hypothesis is supported by the fact that chemical communication is a common sensory modality in aquatic organisms [35,36] and it has been implicated in the initiation and synchronization of spawning in corals [37][38][39][40]. The chemical signal is likely not differences in metabolic waste build-up since there was no effect of H 2 O 2 or ammonia on density-based differences in fission rates, but we acknowledge that other waste products may be co-opted for discerning population density.
Regardless of what regulates asexual reproduction, we know that it must have a downstream effect on initiating physal pinching-which we have identified to occur at a set point along the oral-aboral axis (approx. bottom 27% of parent is lost to clone). Set fission locations have also been identified in Bilateria [12] and Hydra, where budding occurs in the bottom one-third of the body column [41]. This similarity royalsocietypublishing.org/journal/rsos R. Soc. Open Sci. 10: 230152 between Hydra budding location and Nematostella physal pinching location is intriguing because it suggests that there may be a conservation of mechanisms that designate the asexual field in cnidaria.
It is unknown what causes physal pinching at this set location. In Planaria, the pinching that precedes fissioning is neurally regulated [11,12], and it is enticing to speculate that in Nematostella neural regulation controls localized constriction by causing columnar ring muscle to constrict [42] at the pinching site.
Another goal was to establish a framework for exploration of the molecular mechanisms governing asexual reproduction in Nematostella. Understanding the dynamics of transcription during transverse fission is a first step to characterizing this behaviour. Once Nematostella are actively fissioning, we observe changes in gene expression that differ from non-fissioning animals, which we hypothesize are important for complete detachment of the physa from the body column. We identified 1577 transcripts that were differentially expressed between pinching and steady-state animals, including predicted homeobox and notch-like protein transcripts as well as TGFβ and FGF pathway components-all of which have been implicated in asexual reproduction in other bilaterians and cnidarians (see Introduction). Functional experiments and further characterization of differential gene expression throughout the fission process are needed to improve the molecular model of physal pinching in Nematostella, and to determine whether mechanisms of asexual reproduction are homologous despite obvious differences in cloning mechanics between many taxa (e.g. budding in Hydra versus transverse fission in Nematostella). For example, FGF and notch signalling help define the boundary between bud and parent in Hydra [20,21] but it is unknown if these pathways serve a similar role in Nematostella. Additionally, DEGs related to the biological processes of negative cell cycle regulation and cell adhesion should be further explored. The downregulation of the latter probably aids in detachment of the physa, especially given that FGF signalling has been linked to bud detachment in Hydra [22] as well as extracellular matrix degradation and subsequent asexual polyp detachment in the coral Pocillopora acuta [19]. These novel experiments measuring changes in gene expression at the site of constriction and detachment provide us with the first description of transcriptional changes during the process of asexual reproduction in the specific tissue undergoing transition. One limitation of this analysis is isolating only the tissue in question. Future experiments might consider a similar analysis of neighbouring tissue during transverse fissioning to determine if there are differences in transcript expression in the portion of the animal that will develop a head versus the portion of the animal that will develop a physa.
Our findings can lend insight into whether the factors and mechanisms that govern asexual reproduction are conserved across the Cnidaria, but continued work in Nematostella will allow us to draw broader conclusions. The potential to use Nematostella as a model to identify the conditions that promote asexual versus sexual reproduction to better our understanding of cnidarian facultative asexual reproduction should also be explored. It is our opinion that any improvement on our knowledge of the proximate and ultimate factors that govern cnidarian reproduction are of value because it will allow us to better understand a phylum that exhibits incredible reproductive diversity, plays essential ecological roles, and can be applied to management plans or research into strategies aimed to bolster declining wild populations.

Animal husbandry
Nematostella laboratory stocks were maintained in mixed sex groups kept in glass bowls filled with 1/3X artificial salt water (1/3X ASW; pH 8.1-8.2) and kept in a dark 17°C incubator. Water changes were performed weekly, and stocks were fed four times a week with freshly hatched artemia (24-28 h old cultures). To culture artemia for feeding, 1 tsp of brine shrimp eggs (Brine Shrimp Direct-Premium grade brine shrimp eggs) were added to 500 ml of 1/3X ASW and aerated.

Substrate and asexual reproduction
For this experiment, a colony of adult female Nematostella were clonally propagated in glass dishes over several generations in the laboratory as described by Hand & Uhlinger [43] and were formally naive to substrata. One hundred individuals of similar body size were subset evenly into 10 groups and placed in two-cup Pyrex dishes containing either an aquarium gravel substrate (Onyx Sand by Seachem) (n = 5) or no substrate (n = 5). The dishes and gravel were washed and autoclaved prior to usage. After placement royalsocietypublishing.org/journal/rsos R. Soc. Open Sci. 10: 230152 of animals into their respective conditions, the experiment was left to run continuously for two months. During this time, animals were hand fed green mussel (approx. 2 mm 3 in size) three times weekly with sterile forceps followed by weekly water changes where old water was decanted and replaced with fresh 1/3X ASW. The substrate was never removed or washed during the experiment. At the conclusion of the study, animals were re-counted using a mesh window screen to sift through the substrate. Videos of fission behaviour were generated from time lapse images taken at 2 min intervals of an individual kept in substrate and an individual removed from substrate, after pinching had begun, over the course of 11 h and 4 h, respectively. Imaging was done with either an sCMOS camera with an attached lens (Thor Labs Quantalux CS2100M-USB and MVL5M23) or a Nikon dissection microscope (Nikon SMZ1270) with a USB camera attachment (Levenhuk M300). To help visualize fission behaviour of the individual kept in substrate, a clear polydimethylsiloxane (PDMS) based substrate was made. The PDMS (Sylgard184 Ellsworth Adhesives 184 SIL ELAST KIT 0.5KG) was mixed at a ratio of 10:1 ( prepolymer to curing agent), degassed until clear, and cured at 80°C for 20 min before cutting the solid slab into small pieces.

Animal size and asexual reproduction
Longitudinal neurons were counted in NvLWamide-like::mCherry expressing animals at the start of the threeweek long experiment period for the latency to fission, rate of fission, and fission location experiments. Neurons were counted by hand on live animals using a dissection scope equipped with fluorescence (Nikon SMZ 1270) [28]. Animals were relaxed with MgCl 2 (7.14% w/v) in 1/3X ASW before counts were taken. After counting, Nematostella were individually placed in Petri dishes (100 × 15 mm) filled with 38 g of gravel substrate (Super Naturals Aquarium Gravel by CaribSea) and 38 ml of 1/3X ASW and kept in a 22°C incubator with a 12 : 12 light-dark cycle. To control the amount of food animals received based on their size, replicates were fed 10% of their starting longitudinal neuron count in artemia twice a week (e.g. 100 starting longitudinal neurons = 10 artemia twice a week). Each replicate received a water change weekly. Individuals were checked daily for clones by gently looking through the substrate under a dissection scope. If asexual reproduction had occurred, the adult and clone fragments were relaxed with MgCl 2 and longitudinal neurons were counted.

Transcriptomic analysis of animals during physal pinching
Animals were placed in gravel substrate and maintained in the same manner as those from the 'substrate and asexual reproduction' experiment described above. Individuals were either sampled as physal pinching occurred ( pinching) or in steady-state animals (not pinching) (n = 3 pooled individuals per treatment). Since there is currently no established staging protocol for physal pinching in Nematostella, we do not know at what point during the fission process each individual was collected. Animals were put into RNAlater and stored at −20°C until RNA isolation. Total RNA was isolated from all samples using an RNAqueous kit (Ambion) followed by DNAse treatment with a DNA-free kit (Invitrogen). RNA quantity was assessed using both a Qubit Fluorometer with associated RNA Broad Range kit (Thermo Fisher Scientific) and an Agilent 2100 Bioanalyzer using a 'pico' kit, each according to their respective manufacturer's instructions. Libraries were prepared from 14.6 to 100 ng of high-quality total RNA by the Sequencing and Discovery Genomics core at the Stowers Institute for Medical Research following manufacturer protocols for the Illumina TruSeq Stranded mRNA Library Prep kit and Illumina TruSeq RNA Single Indexes Sets A and B. Resulting short fragment libraries were checked for quality and quantity using an Agilent Bioanalyzer and a Qubit Fluorometer. Libraries were pooled at equal molar concentrations, quantified and sequenced as 75 bp single reads on a high-output flow cell using an Illumina NextSeq 500 instrument and NextSeq Control Software 2.2.0.4. Following sequencing, Illumina Primary Analysis v. RTA 2.4.11 and bcl2fastq2 v. 2.20 were run to demultiplex reads for all libraries and generate FASTQ files. Sequencing produced an average of 31 062 832 reads per sample, ranging from 26 312 701 to 33 855 476. Sequencing results were deposited in the NCBI GEO database (GSE223794). The filtered reads were aligned to the Nvec200 Nematostella genome [29] with STAR aligner (v. 2.7.3a) using Ensembl 99 gene models. Reads mapping to the reference at the same start position with 100% alignment identity of the transcript were regarded as PCR duplicates and discarded. TPM values were generated using RSEM (v. 1.3.0).
All statistical analyses were completed in the R environment (R3.5.0). The raw read-counts-per-gene table was filtered to contain only genes with at least three counts and used to perform outlier detection with the arrayQualityMetrics R package [44] with three different outlier detection methods: distance royalsocietypublishing.org/journal/rsos R. Soc. Open Sci. 10: 230152 between arrays, boxplots and MA plots. No samples were excluded from downstream analysis. The filtered count data were analysed using the R package DESeq2 [45], which produced a normalized counts-pergene table used in subsequent analyses. Normalized count data were transformed using the rlog function. Corrections were made for multiple testing at a cut-off of 0.1. Using a generalized linear model with Wald pairwise comparisons, DEGs were determined between samples that were pinching or not pinching. The VEGAN package was used to visualize and calculate significance and generate a PCA of variancestabilized gene expression [46]. Functional summaries of DEGs were determined by rank-based GO enrichment analysis, using signed, unadjusted log-transformed p-values ( positive if upregulated, negative if downregulated) with the GO_MWU package (https://github.com/z0on/GO_MWU). This method uses a Mann-Whitney U test and measures whether each GO term is significantly enriched in upregulated or downregulated genes based on their delta rank (quantitative shift in rank) rather than looking for GO terms among significant genes only.

Population density and asexual reproduction
Laboratory cultured wild-type Nematostella from mixed sex populations were used for all population density experiments. Prior to being assigned to a high-or low-density treatment (n = 40 or 4 animals, respectively), all animals were relaxed in MgCl 2 and their body length was measured using a clear plastic ruler placed under their Petri dish. Within each replicate, average animal length and length range was controlled across population density treatments. Animals were placed approximately equidistant in their assigned treatment Petri dish (145 × 25 mm) filled with 100 g of gravel substrate (Super Naturals Aquarium Gravel by CaribSea) and 200 ml of 1/3X ASW and kept in a 22°C incubator with a 12 : 12 light-dark cycle for the three-week long experiment period. High-and lowdensity treatments were fed 400 and 40 µl of artemia culture twice a week, respectively. Water changes were performed weekly and a 20 ml water top-off occurred mid-week. H 2 O 2 density treatments had H 2 O 2 added to their 1/3X ASW every time water was changed or topped off (final concentration 0.00025% or 82 µM). Replicates were checked daily for egg clutches which were removed if found. Ammonia levels were tested on water removed from each treatment dish during weekly water changes using an API brand NH 3 /NH 4 + test kit. At the end of the experiment period, all animals were removed from each dish and counted to determine the number of clones that were generated per starting individual. Dishes were checked multiple times under a dissection scope to ensure no animals were overlooked.

Transcriptomic analysis of population density animals
Total RNA was extracted from animals at the end of the population density experiment period. In total, 10 high-and nine low-density animals were collected from four different population density replicates. None of the sampled individuals were undergoing physal pinching at the time of collection. Each animal was flash frozen in liquid nitrogen, processed with a liquid-nitrogen-cooled mortar and pestle, stabilized in 1 ml of TRI Reagent (Sigma-Aldrich), and then stored at −80°C until total RNA was isolated as previously published [47]. RNA concentration and purity were measured using a NanoDrop spectrophotometer (all samples OD:260/280 and OD:260/230 ≥ 2) and sample integrity was confirmed using an Agilent 2100 Bioanalyzer (all samples RIN ≥ 8.0) following manufacturer's instructions. RNA was stored at −80°C until needed. Sample cDNA libraries were prepared from an average of 21 to 274 ng of high-quality RNA. Library preparation and sequencing were done by Novogene Inc. Sample libraries were constructed following manufacturer protocols using an Illumina Stranded mRNA Prep kit, which uses poly A enrichment to purify mRNA from total RNA, and an Illumina NovaSeq 6000 was used for sequencing ( paired end, 150 bp reads). An average of 50 324 344 reads were produced per sample, ranging from 38 628 894 to 70 726 914. Reads were quality filtered and adapter trimmed using Trimmomatic [48]. A read was discarded if there was adapter contamination, if it contained more than 10% uncertain nucleotides, or if there were more than 50% low-quality nucleotides (Qscore ≤ 5). Sequencing data can be retrieved from the NCBI SRA database (PRJNA925453). Sequencing data were uploaded to the Galaxy web platform, and the public server at usegalaxy.org was used to map and analyse the sequencing data [49]. Filtered reads were mapped and quantified using the Salmon quant tool with default parameters (v. 1.5.1) [50] using the Nematostella NVE transcriptome v. 2.0 [34]. Transcripts per million were used to quantify transcript abundance. The count tables generated with Salmon quant were used for differential transcript analysis between high-and low-royalsocietypublishing.org/journal/rsos R. Soc. Open Sci. 10: 230152 density samples using the Bioconductor DESeq2 tool using default settings (v. 2.11.40.7) [45]. A false discovery rate of p ≤ 0.05 was used to determine differential expression (adjusted for multiple testing using the Benjamini-Hochberg procedure).
Ethics. This work did not require ethical approval or permits. Animals were maintained, handled and studied in accordance with the advice provided by the Association for the Study of Animal Behaviour/Animal Behaviour Society Guidelines for the Use of Animals in Research (Animal Behaviour, 2020, 159, I-XI).
Data accessibility. All data were generated and analysed by the authors of this study. Individual replicate data are provided in the figures of the paper for all behaviour experiments as well as in electronic supplementary material, data file 1, and on Figshare [51]. RNA sequencing results for pinching versus steady-state animals were deposited in the NCBI GEO database (GSE223794). Population density RNA sequencing data have been uploaded to the NCBI SRA database (PRJNA925453).