Rocío Pérez-Portela
Carlos Leiva- 1Departament de Biologia Evolutiva, Ecologia i Ciències Ambientals, Universitat de Barcelona, Barcelona, Spain
- 2Research Institute of Biodiversity (IRBIO), Universitat de Barcelona, Barcelona, Spain
- 3University of Guam Marine Laboratory, Mangilao, GU, United States
The abundance of the black sea urchin, Arbacia lixula, has been increasing during the last decades likely related to global warming. This thermophilous species has a leading role in maintaining marine barrens in the Mediterranean with the consequent negative impact on coastal rocky ecosystems due to its grazing activity. In this study, we used transcriptomic data from coelomocytes (the cell effectors of the immune system) of females and males of this sea urchin to study potential differences in performance between sexes under laboratory conditions. Differential adaptations, responses to environmental stressors, and resistance against pathogens between sexes may lead to different outcomes in the ongoing expansion of this species in the Mediterranean Sea. Differential expression analyses demonstrated the existence of 120 transcripts, corresponding to 119 genes and two isoforms of the same gene, differentially expressed between coelomocytes of females and males, being 73 up-regulated in males and 47 up-regulated in females. The differential expression patterns were retrieved from a diversity of genes that play different roles related to the immune response due to their antibacterial activity, immune cell activation, cell to cell interaction, intracellular signaling, and detoxification functioning, among others. Our results point out a higher energetic demand of male coelomocytes due to a higher immune activity than females, whereas females have more efficient molecular systems to avoid oxidative stress caused by infections. In conclusion, our study provides evidence of sex-based differences in the expression of genes related to the immune and stress responses in coelomocytes of the sea urchin A. lixula.
Introduction
Sea urchins are among the most relevant keystone species in marine ecosystems, with a plethora of examples of the structuring effects of their grazing activity in sublittoral habitats, from tropical seas to kelp forests (e.g., Palacín et al., 1998; Edmunds and Carpenter, 2001; Hernández et al., 2008; Eklöf et al., 2008; Ling et al., 2009). Their overgrazing activity in coastal rocky assemblages causes regime shifts from productive macroalgal beds to impoverished urchin barrens, with anthropogenic stressors such as overfishing and climate change strengthening the resilience of urchin barrens (Hernández et al., 2008; Ling et al., 2009; Ling et al., 2015; Bernal-Ibáñez et al., 2021). Due to the invaluable ecosystem services provided by macroalgal and kelp beds (e.g., Smale et al., 2013; Rocha et al., 2015; Krause-Jensen and Duarte, 2016), it is of utmost importance to understand all details of the biology of the habitat-modifying echinoid species.
In the Mediterranean Sea, the two sea urchin species producing barrens are Paracentrotus lividus and Arbacia lixula (Palacín et al., 1998; Micheli et al., 2005), having A. lixula a leading role in maintaining the barren state once phase shift occurred (Bonaviri et al., 2011). Arbacia lixula is a highly mobile species on open rocky surfaces that “scratches” on the substrate to feed on encrusting coralline algae and benthic invertebrates. This feeding behavior plays an important role in avoiding the macroalgae re-colonization, which finally promotes the persistence of large bare areas (Bonaviri et al., 2011; Agnetta et al., 2013). As a relatively new colonizer in expansion in the Mediterranean Sea (Wangensteen et al., 2012; Pérez-Portela et al., 2019), the thermophilous A. lixula has increased its abundance during the last century in the Mediterranean (Petit et al., 1950; Francour et al., 1994). Following this trend, and due to its superior physiological performance in warmer conditions (Privitera et al., 2011; Wangensteen et al., 2013a), global warming will likely favor A. lixula expansion and dominance in the whole Mediterranean basin (Gianguzza et al., 2011), with its associated negative effects due to its role on the maintaining barrens areas, already evident in the Southern-Central Mediterranean (Bonaviri et al., 2011).
Given this projection, an enhancing number of studies are using A. lixula as a model to investigate its genetic variation and adaptation patterns (Wangensteen et al., 2012; Pérez-Portela et al., 2019; Carreras et al., 2021), trophic ecology (Wangensteen et al., 2011), biological cycles (Wangensteen et al., 2013a) and responses to future climatic conditions and stressors (e.g. Pérez-Portela et al., 2020a). During the last decade, a wide array of techniques and approaches has been implemented, for instance, to measure its biological response to underwater high-frequency noise (Vazzana et al., 2020), neurotoxicity (Maisano et al., 2015), transcriptomic responses to temperature changes (Pérez-Portela et al., 2020a), and the effects of environmental changes on larval development, settlement and juvenile survival (Wangensteen et al., 2013b; Gianguzza et al., 2014). For these studies, understanding all factors that may affect the individual response to different treatments is essential, but in most of them, the variable sex is not considered despite its potential relevance in modulating the individual responses.
Genomic, transcriptomic, and proteomic studies of echinoderm coelomocytes have demonstrated the large complexity and flexibility of their innate immune system (e.g. Dheilly et al., 2013; Smith et al., 2018; Zhang et al., 2019). In this work, we use transcriptomic data from coelomocytes of females and males of A. lixula to study potential differences in the immune system between sexes. The coelomocytes, the cell effectors of the immune system, consist of heterogeneous cell populations contained in the coelomic fluid, which are also involved in other physiological functions such as the stress response (e.g. Matranga et al., 2000; Pérez-Portela et al., 2016; Smith et al., 2018). In turn, these functions are considered fundamental elements of adaptation to new environments and face environmental shifts associated with climate change (e.g., Palumbi et al., 2014; Cunning et al., 2018; Liu et al., 2018). Previous studies in the sea urchin P. lividus, regarding the differential immune response according to sex, demonstrated that the immune system is more active in females than in males, with higher levels of cytotoxicity and hemolysis activity in females, likely related to a different cell composition within the coelomic fluid and a greater number of phagocytes and uncolored spherulocytes in females (Arizza et al., 2013). Additionally, other studies with marine invertebrates have also evidenced stronger immune responses in females than males (Nguyen et al., 2018, and references therein), although exceptions also exist (Jiang et al., 2017). Indeed, focusing on gene expression responses, more recent studies have identified patterns of sex-specific differentially expressed genes in different marine invertebrates (e.g., Li et al., 2019; Artal et al., 2020; Kadiene et al., 2020), which are also present in somatic tissues (e.g., Rotllant et al., 2017). Hence, it is crucial to understand the sex-specific expression of A. lixula coelomocytes, as different adaptation performance, responses to environmental stressors, and resistance against pathogens of males and females may lead to different outcomes in the ongoing expansion of this species in the Mediterranean Sea. Arbacia lixula lacks sexual dimorphism and has populations with a sex ratio close to 1:1, although slightly deviations in favor of males (Elakkermi et al., 2021) or females (Klaoudatos et al., 2022) have been detected in some populations.
Significant differences in transcriptomic patterns and gene expression have been described between gonads of males and females in A. lixula based on differential functioning of these organs, ovaries and testis (Pérez-Portela et al., 2016), but in the present study, we compare differences between sexes using the same somatic tissue type. Our initial hypothesis is that, following that observed in P. lividus, females of A. lixula would immunologically be more active than males, and therefore both sexes would display differential patterns of transcriptomic expression in coelomocytes.
Material and Methods
Raw reads and a reference transcriptome used in this study were obtained from Pérez-Portela et al, (2020a), available at Mendeley Data https://doi.org/10.17632/5673n552yj.1 and NCBI BioProject number PRJNA642021 (BioSamples SAMN15375803- SAMN15375808). We used transcriptomic information of three females and three males maintained under control conditions in laboratory.
Sampling and Experimental Design
Details of the sampling and experimental design are described in Pérez-Portela et al, (2020a). Briefly, adult specimens of A. lixula (~48-52 mm in diameter) were collected from a subtidal population of Blanes, in Cala Sant Francesc (41°40′22.47″N, 2°48′10.81″E, North-Western Mediterranean), to minimize potential variation due to genetic divergence among populations. Once in the laboratory, sea urchins were maintained in an aquarium with airflow and flow-through running seawater at 13˚C. For coelomocytes collection and RNA extraction, sea urchins were pulled out from the aquaria, quickly dissected under RNAase-free conditions, and coelomocyte fluid was collected. Sex determination in sea urchins was done by visual examination of the gonads, and posterior histological analyses. Coelomocytes of each specimen (a total of six: three females and three males) were gathered from the coelomic fluid and RNA was extracted following a protocol optimized elsewhere (Pérez-Portela and Riesgo, 2013; Pérez-Portela et al., 2016). High-quality mRNA was isolated and cDNA libraries were constructed with an average insert size of 300 pb. Libraries were multiplexed and sequenced in an Illumina HiSeq2000 Sequencer (101 base paired-end).
The software FASTQC v. 0.10.0 (www.bioinformatics.babraham.ac.uk) was used to visualize and measure the quality of the raw reads generated. Adapters and bases with low quality (phred scores <33) were trimmed off, and a length filter was applied retaining sequences of >25 bases using TrimGalore v. 0.2.6 (www.bioinformatics.babraham.ac.uk).
Differential Expression Analyses and Annotation
Reads from each sample were aligned against a reference transcriptome for coelomocytes of A. lixula (reference_assembly_CT+T22.fasta, available at https://data.mendeley.com/datasets/5673n552yj/1, Pérez-Portela et al., 2020b). This reference transcriptome was de novo assembled by Pérez-Portela and coauthors (2020b) from 12 different coelomocyte samples of A. lixula, and it was the most complete one among those available for coelomocytes in this species (Pérez-Portela et al. 2020b). Only transcripts over 500 pb from the reference transcriptome were used for mapping our reads, and only paired reads after trimming were mapped using Bowtie2 v. 2.2.1 (Langmead and Salzberg, 2012) as implemented in Trinity (Grabherr et al., 2011). RSEM v. 1.2.11 (Li and Dewey 2011) was then run to generate a table with read counts, and unmapped reads were discarded. We retained information on differential expression of all isoforms detected for a given gene (or component) because different isoforms of the same gene may have different functions.
Differential expression (DE) analyses between females and males (three replicates each) were performed with the package DESeq2 (Love et al., 2014) in R v 4.0 (www.CRAN.R-project.org; R Development Core Team 2008). Before analyses of differential gene expression, read counts were internally normalized by the “DEseq()” function, and then a negative binomial model was fit to accurately estimate differential expression. The significance value for multiple comparisons was adjusted to 0.01 for False Discovery Rate (FDR) with the function “padj” (Benjamini- Hochberg adjustment) as implemented in DESeq2. Visualization of the significant outcomes of isoforms differentially expressed (up- and down-regulated) between females and males was obtained with a heatmap performed with the gplots package of R (Warnes et al., 2009).
Using the Gene Ontology (GO) annotation results for the reference transcriptome, we obtained both the gene annotations and the GO terms associated with the differentially expressed isoforms. GO terms and their associated p-value (adjusted) were then summarized in the REVIGO webserver (Supek et al., 2011), and graphically represented with the treemap R package (Tennekes and Ellis, 2017). The size of the rectangles was adjusted to reflect the adjusted p-value to visually detect the most significant terms. Differentially expressed isoforms without blast hit, unknown function and/or without annotation were isolated, and the InterProScan 5 software (Jones et al., 2014), that predicts protein family membership and the presence of functional domains and sites, was run at the SUPERFAMILY level (De Lima Morais et al., 2011). The InterProScan was run as implemented in the Blast2GO software v 3.31.1 (Conesa et al., 2005).
Results
Differential Expression Between Sexes
A total of 283,880,164 reads from all samples, and 66,985 transcripts from the reference transcriptome, including all genes and isoforms, were used for the differential gene expression analysis in this study. All samples had over 23,000,000 reads (ranging from 23,074,766 to 84,053,170) after filtering and trimming. A total of 120 transcripts, corresponding to 119 genes and two isoforms of the same gene, were differentially expressed (p-adjusted value < 0.01) between coelomocytes of females and males, being 73 up-regulated in males and 47 up-regulated in females. The heatmap that describes the hierarchical clustering of all transcripts related to their expression differences between sexes, including three replicates for each sex, is displayed in Figure 1. The heatmap shows a clear clustering of all replicates per sex.
Figure 1 Heatmap based on the differentially expressed transcripts (DE) between females and males of A. lixula. The clustering on the top of the heatmap represents the similarity among the samples analyzed. F1-F3 and M1-M3 represent the three female and three male replicates, respectively. The term “Value” represents the gene expression value.
Functional Annotation
Out of the 120 differentially expressed transcripts, 40 (~33.3%) had blast hit against known genes. The InterProScan analyses could only predict protein domains of three uncharacterized transcripts. A complete list of differentially expressed annotated transcripts is presented in Table 1, including transcript identification code (id), log2fold change, p- adjusted value (padj), and gene annotation, and GO-terms.
Table 1 List of the DE annotated transcripts between sexes in A. lixula.
Among the transcripts differentially expressed between females and males, we can find some codifying genes described in other studies as directly related to the immune response, such as: the zinc finger ccch domain-containing protein 14-like (log2fold change= -8.53), the ras-related c3 botulinum toxin substrate 1-like (log2fold change= -8.17), the hepatic lectin (log2fold change= -2.68), and the bactericidal permeability-increasing (log2fold change= -3.3), all of them up-regulated in males. For females, we found up-regulated the echinonectin (log2fold change= 9.07) and ferric-chelate reductase 1 like (log2fold change= 6.78).
Other genes codifying for potential proteins involved in the immune activity, due to its role in cell to cell interaction and intracellular signaling, were also DE between sexes, including the activating signal cointegrator 1 complex subunit 3-like (log2fold change= -7.35) and the low-density lipoprotein receptor-related protein 2-like isoform 2 (log2fold change= -1.00) both up-regulated in males. Transcripts related to the transmembrane transport functions as two isoforms of the major facilitator superfamily domain-containing protein 7-a-like (isoforms 1 and 2, please see c275529_g1_i1 and c275529_g1_i2 in Table 1) were also identified as differentially expressed between sexes (log2fold change= 8.19 and -8.05, respectively). Several genes codifying for integral components of the cellular membranes were also DE, including the transmembrane protein 85-like and alphap integrin precursor, both up-regulated in males (log2fold change= -8.22 and -9.03, respectively), and the ADP, ATP carrier protein-like gene up-regulated in females (log2fold change= 7.51), among others. As expected, the GO terms associated with these DE genes were also different between males and females. Out of the 35 GO categories recovered, 22 were derived from genes up-regulated in males, and 13 from genes up-regulated in females. Figure 2 depicts the up-regulated GO categories associated with the DE genes. The GO term categories found can give us a general idea of the biological processes, cellular components, and molecular functions differentially represented between sexes. Among the most significant GO terms up- and/or down-represented in females and males, we find some terms related to intracellular signaling and transport such as “transmembrane transport”, “protein binding”, “integral component of membranes”, “membrane”; those related with cell to cell interaction such as “cell adhesion”; with the stress-responses such as “ferric-chelate reductase activity” and “metal ion binding”; with energy metabolism as “peptide cross-linking”, “NADH dehydrogenase complex” and “protein-glutamine gamma-glutamyl transferase activity”; and other general categories such as the biological processes of “RNA-dependent DNA biosynthetic process”, and their associated cellular components and molecular functions.
Figure 2 Gene ontology (GO) treemaps for DE genes annotated in A. lixula. Those up-regulated in males and females are represented separately for Biological processes, Cellular components, and Molecular functions. The size of the rectangles reflects the adjusted p- values associated with the differentially regulated categories.
Discussion
Our transcriptomic analysis from coelomocytes of the sea urchin A. lixula demonstrated the existence of significant differences in the expression of some genes, and isoforms, between females and males. These patterns of differential expression were retrieved from a diversity of genes that seem to play different roles during the immune and stress responses in sea urchins.
A study in P. lividus identified a different composition of the coelomocyte populations between sexes, with females being immunologically more active (Arizza et al., 2013). In our study, we detect differences in transcriptional levels of expression between females and males of A. lixula, but a higher number of up-regulated genes are associated with the innate immune response in males. Although the relationship between levels of transcription and abundance of proteins, and the functional consequences of differential expression of some particular genes are, on some occasions, still unclear (Feder and Walser, 2005), some studies have shown the correlation between expression levels and protein abundance under experimental conditions (Maier et al., 2009). Therefore, changes in gene expression can be used as an approximation of the physiological performance of a tissue and/or organism under laboratory-controlled conditions (Feder and Walser, 2005; Buckley et al., 2006; Schoville et al., 2012), and so, in this work, the number of immune-related genes up-regulated can be considered as a proxy to infer the immune activity. Hence, A. lixula results suggest that, contrary to that described in P. livudus, males might be immunologically more competent than females, or at least, both sexes use different immune pathways. Immune-competence analyses with sea cucumber Apostichopus japonicus, also evidenced a more powerful immune system in males over females (Jiang et al., 2017), which is in line with our results in A. lixula.
Although many genes are potentially involved in the innate immune response in sea urchins (e.g. Dheilly et al., 2013; Zhang et al., 2019), in this study, we only discuss some of the most relevant ones and/or those playing a well-known role in the immunity, cell adhesion, cell signaling, stress, and detoxification. Regarding the genes directly involved in the immune response and cell adhesion, at least, four of them were up-regulated in males: the bactericidal permeability-increasing protein, the zinc finger ccch domain-containing protein 14-like, ras-related c3 botulinum toxin substrate 1-like, and a hepatic lectin. The bactericidal permeability-increasing gene encodes for a lipopolysaccharide-binding protein with potent cytotoxicity against Gram-negative bacteria, which has been isolated in a diversity of echinoderm tissues (Shao et al., 2015), but a DE between sexes has never been described before. The zinc finger ccch domain-containing protein 14-like has a crucial role in the regulation of cytokines production, immune cell activation, and immune homeostasis (Fu and Blackshear, 2017), and it is associated with the molecular function GO term “protein binding”. The ras-related c3 botulinum toxin substrate 1-like gene showed significant differences in expression between males and females and its foldchange was among the lowest in the whole dataset. This gene codifies for a plasma membrane protein, involved in cell signaling, a member of the Rho family of small GTPases that regulates processes of inflammation and apoptosis (Hirohama et al., 2020), and stimulates the NADPH oxidase activity and phagocytosis in macrophages (Dheilly et al., 2013). The hepatic lectin gene codifies for a molecule secreted by the coelomocytes to the coelomic fluid as a part of the immune humoral system. Lectins are a large group of proteins and glycoproteins, soluble or integral components of membranes, that act as mediators of the immune system in invertebrates, including echinoderms, due to their function in the cell adhesion, and cell to cell and cell-extracellular matrix interactions (Dheilly et al., 2013). They can specifically recognize and bind to the microbial surface and lead to phagocytosis (Dheilly et al., 2013; Zhang et al., 2019).
In females, two genes potentially related to the immune response were detected as “up regulated”: the echinonectin and the ferric-chelate reductase 1-like. The echinonectin gene encodes for an echinoderm adhesion protein, a lectin, that was initially described during the embrionary development, and that also plays an important role in skeleton formation (Mann et al., 2008). Nevertheless, more recent studies might suggest its role in the defense against pathogens as other lectins (Dheilly et al., 2013) due to its important role in cell adhesion. The ferric-chelate reductase 1-like gene encodes a molecule linked to the membrane function and iron regulation that reduces Fe+3 to Fe+2, and is associated to the molecular functioning GO terms “ferric-chelate reductase activity” and “metal ion binding”. In sea urchins, iron plays a critical role in the immune defense against pathogens stimulating cell proliferation, and nitric oxide and reactive oxygen elements (Dheilly et al., 2013). In vertebrates, iron sequestration avoids it to be used by microorganisms (Brock and Mulero, 2000; Ganz, 2009), a function that iron-binding proteins could be also playing in sea urchins (Dheilly et al., 2013; Zhang et al., 2019). However, the specific involvement of the ferric-chelate reductase 1 is not mentioned in previous studies on echinoderms. On the other hand, some studies with sea stars have detected expression of this gene in females and signals of selection on it, promoting adaption to reduce oxidative stress (Guerra et al., 2021). This gene is particularly interesting since in other marine invertebrates, such as the mollusk Ruditapes philippinarum, females showed more efficient pathways against oxidative stress caused by infections than males (Matozzo and Marin, 2010), and so our results from A. lixula could also support this idea.
We also highlight the importance of gene isoforms in the coelomocytes’ functioning. We detected the existence of two isoforms of the major facilitator superfamily domain-containing protein 7-a-like coding gene, involved in the biological processes of “transmembrane transport” and “integral component of membranes”, with an opposite pattern of expression between sexes and almost equivalent foldchange. This gene expression pattern suggests the sex-specificity of these two isoforms. Previous transcriptomic analyses in A. lixula unraveled a rich diversity of gene isoforms (Pérez-Portela et al., 2016; Pérez-Portela et al., 2020a), tissue-specificity (Pérez-Portela et al., 2016), or even opposite patterns of expression under different thermal treatments in coelomocytes (Pérez-Portela et al., 2020a). Hence, it does not seem rare to detect the existence of isoforms with opposite patterns of differential expression between sexes, as in other marine invertebrates, such as abalones, the presence of gene isoforms was also sex-dependent (Kim et al., 2017).
In general, and independently of the particular genes up- or down- regulated per sex, different GO terms inferred from A. lixula were also identified in the immune response against bacteria artificially inoculated in Strongylocentrotus intermedius (Zhang et al., 2019), and coelomocytes’ proteomic analyses of Strongylocentrotus purpuratus (Dheilly et al., 2013), including the “transmembrane transport”, “membrane”, “integral component of membrane”, “cell adhesion”, “protein binding” and “metal ion binding”, among others. These similar results with other sea urchin species support the relevant role of the cell interactions throughout the cellular membranes for the immune response in both sexes of A. lixula. Another interesting molecular pathway potentially linked to immune functions is the lipid metabolism due to its relevance in the development of the immunological response (Miller et al., 2003). Within it, we detected the low-density lipoprotein receptor-related protein 2-like gene differentially expressed between sexes of A. lixula. These cell receptors are involved in capturing low-density lipoproteins for intra- and extra-cellular signaling and are also identified in coelomocytes of S. purpuratus (Dheilly et al., 2013). Finally, we identified the up-regulation of the transglutaminase-like protein and nadh dehydrogenase codifying genes in males, both associated with the energy metabolism and the GO terms “peptide-cross linking”, “NADH dehydrogenase complex” and “protein-glutamine gamma-glutamyl-transferase activity”. This finding, as in other studies on echinoderms (e.g. Dheilly et al., 2013; Zhang et al., 2019), points out a higher energetic demand in coelomocytes of males due to a higher immune activity.
In conclusion, our study provides evidence of sex-based differences in the expression of some genes directly and indirectly related to the immune and stress responses. The potential higher immune activity in males than in females may lead to an advantage against pathogens but under the cost of a higher energetic demand. On the other hand, females seem to be better adapted to the oxidative stress associated with infections. Both sexes did not display significant differences in transcriptomic responses to thermal stress under controlled conditions in an aquarium (Pérez-Portela et al., 2020a) but did it under control conditions. Whether the absence of transcriptomic divergence between sexes under stressful thermal conditions is due to the low replication of one of the sexes (three replicates per sex) or the dilution of the potential of those differences under stress, is a question that we cannot answer with the current data. Future studies should focus on the differential immune response against pathogens between sexes under variable environmental conditions, to provide a complete picture of the relevance of these sex-associated differences in transcriptional patterns between females and males of A. lixula.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI [accession: PRJNA642021].
Author Contributions
RP-P. designed this study and performed the analyses. RP-P. and CL created figures and tables and wrote the manuscript. Both authors agreed with the current version of the manuscript.
Funding
This research was funded by the Spanish Government projects: Adaptive- PGC2018-100735-B-I00 (MCIU/AEI/FEDER, UE) and PopCOmics CTM2017-88080 (MINECO/AEI/FEDER, UE), and a Ramón y Cajal contract to R.P-P. (RYC2018-025070-I).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We thank Dr A. Riesgo for her helpful comments during the analysis of the data. This research was funded by the Spanish Government projects: Adaptive- PGC2018-100735-B-I00 (MCIU/AEI/FEDER, UE) and PopCOmics CTM2017-88080 (MINECO/AEI/FEDER, UE), and a Ramón y Cajal contract to R.P-P. (RYC2018-025070-I).
References
Agnetta D., Bonaviri C., Badalamenti F., Scianna C., Vizzini S., Gianguzza P.. (2013). Functional traits of two co-occurring sea urchins across a barren/forest patch system. J. Sea Res. 76, 170–177. doi: 10.1016/j.seares.2012.08.009
Arizza V., Vazzana M., Schillaci D., Russo D., Giaramita F. T., Parrinello N. (2013). Gender Differences in the Immune System Activities of Sea Urchin Paracentrotus Lividus. Comp. Biochem. Physiol. Part A. Mol. Integr. Physiol. 164, 447–455. doi: 10.1016/j.cbpa.2012.11.021
Artal M. C., Pereira K. D., Luchessi A. D., Okura V. K., Henry T. B., Marques-Souza H., et al. (2020). Transcriptome Analysis in Parhyale Hawaiensis Reveal Sex-Specific Responses to AgNP and AgCl Exposure. Environ. pollut. 260, 113963. doi: 10.1016/j.envpol.2020.113963
Bernal-Ibáñez A., Cacabelos E., Melo R., Gestoso I. (2021). The Role of Sea-Urchins in Marine Forests From Azores, Webbnesia and Cabo Verde: Human Pressures, Climate-Change Effects and Restoration Opportunities. Front. Mar. Sci. 8. doi: 10.3389/fmars.2021.649873
Bonaviri C., Fernández T. V., Fanelli G., Badalamenti F., Gianguzza P. (2011). Leading Role of the Sea Urchin Arbacia Lixula in Maintaining the Barren State in Southwestern Mediterranean. Mar. Biol. 158, 2505–2513. doi: 10.1007/s00227-011-1751-2
Brock J. H., Mulero V. (2000). Cellular and Molecular Aspects of Iron and Immune Function. Proc. Nutr. Soc 59, 537–540. doi: 10.1017/S002966510000077X
Buckley B. A., Gracey A. Y., Somero G. N. (2006). The Cellular Response to Heat Stress in the Goby Gillichthys Mirabilis: A cDNA Microarray and Protein-Level Analysis. J. Exp. Biol. 209, 2660–2677. doi: 10.1242/jeb.02292
Carreras C., Ordóñez V., García-Cisneros À., Fuentes W., Palacín C., Pascual M., et al. (2021). The Two Sides of the Mediterranean: Population Genomics of the Black Sea Urchin Arbacia Lixula (Linnaeus 1758) in a Warming Sea. Front. Mar. Sci. 8. doi: 10.3389/fmars.2021.739008
Conesa A., Götz S., Garcia-Gomez J. M., Terol J., Talon M., Robles M. (2005). Blast2GO: A Universal Tool for Annotation, Visualization and Analysis in Functional Genomics Research. Bioinformatics 21, 3674–3676. doi: 10.1093/bioinformatics/bti610
Cunning R., Bay R. A., Gillette P., Baker A. C., Traylor-Knowles N. (2018). Comparative Analysis of the Pocillopora Damicornis Genome Highlights Role of Immune System in Coral Evolution. Sci. Rep. 8, 1–10. doi: 10.1038/s41598-018-34459-8
De Lima Morais D. A., Fang H., Rackham O. J. L., Wilson D., Pethica R., Chothia C., et al. (2011). SUPERFAMILY 1.75 Including a Domain Centric Gene Ontology Method. Nucleic Acids Res. 39, D427–D434. doi: 10.1093/nar/gkq1130
Dheilly N. M., Raftos D. A., Haynes P. A., Smith L. C., Nair S. V. (2013). Shotgun Proteomics of Coelomic Fluid From the Purple Sea Urchin, Strongylocentrotus Purpuratus. Dev. Comp. Immunol. 40, 35–50. doi: 10.1016/j.dci.2013.01.007
Edmunds P. J., Carpenter R. C. (2001). Recovery of Diadema Antillarum Reduces Macroalgal Cover and Increases Abundance of Juvenile Corals on a Caribbean Reef. Proc. Natl. Acad. Sci. U.S.A. 98, 5067–5071. doi: 10.1073/pnas.071524598
Eklöf J. S., de la Torre-Castro M., Gullström M., Uku J., Muthiga N., Lyimo T., et al. (2008). Sea Urchin Overgrazing of Seagrasses: A Review of Current Knowledge on Causes, Consequences, and Management. Estuar. Coast. Shelf. Sci. 79, 569–580. doi: 10.1016/j.ecss.2008.05.005
Elakkermi M., Mezali K., Soualili D. L. (2021). Interpopulation Variability of the Reproductive Cycle of Arbacia Lixula (Echinodermata: Echinoidea) in the Mostaganem Shallow-Water Area (South-Western Mediterranean). Reg. Stud. Mar. 45, 101810. doi: 10.1016/j.rsma.2021.101810
Feder M. E., Walser J. C. (2005). The Biological Limitations of Transcriptomics in Elucidating Stress and Stress Responses. J. Evol. Biol. 18, 901–910. doi: 10.1111/j.1420-9101.2005.00921.x
Francour P., Boudouresque C. F., Harmelin J. G., Harmelin-Vivien M. L., Quignard J. P. (1994). Are the Mediterranean Waters Becoming Warmer? Information From Biological Indicators. Mar. pollut. Bull. 28, 523–526. doi: 10.1016/0025-326X(94)90071-X
Fu M., Blackshear P. J. (2017). RNA-Binding Proteins in Immune Regulation: A Focus on CCCH Zinc Finger Proteins. Nat. Rev. Immunol. 17, 130–143. doi: 10.1038/nri.2016.129
Ganz T. (2009). Iron in Innate Immunity: Starve the Invaders. Curr. Opin. Immunol. 2, 63–67. doi: 10.1016/j.coi.2009.01.011
Gianguzza P., Agnetta D., Bonaviri C., Di Trapani F., Visconti G., Gianguzza F., et al. (2011). The Rise of Thermophilic Sea Urchins and the Expansion of Barren Grounds in the Mediterranean Sea. Chem. Ecol. 27, 129–134. doi: 10.1080/02757540.2010.547484
Gianguzza P., Visconti G., Gianguzza F., Vizzini S., Sarà G., Dupont S. (2014). Temperature Modulates the Response of the Thermophilous Sea Urchin Arbacia Lixula Early Life Stages to CO2-Driven Acidification. Mar. Environ. Res. 93, 70–77. doi: 10.1016/j.marenvres.2013.07.008
Grabherr M. G., Haas B. J., Yassour M., Levin J. Z., Thompson D. A., Amit I., et al. (2011). Trinity: Reconstructing a Full-Length Transcriptome Without a Genome From RNA-Seq Data. Nat. Biotechnol. 29, 644–652. doi: 10.1038/nbt.1883
Guerra V. I., Haynes G., Byrne M., Hart M. W. (2021). Selection on Genes Associated With the Evolution of Divergent Life Histories: Gamete Recognition or Something Else? Evol. Dev. 23, 423–438. doi: 10.1111/ede.12392
Hernández J. C., Clemente S., Sangil C., Brito A. (2008). The key role of the sea urchin Diadema aff. antillarum in controlling macroalgae assemblages throughout the Canary Islands (eastern subtropical Atlantic): an spatio-temporal approach. Mar. Environ. Res. 66, 259−270. doi: 10.1016/j.marenvres.2008.03.002
Hirohama D., Kawarazaki W., Nishimoto M., Ayuzawa N., Marumo T., Shibata S., et al. (2020). PGI2 Analog Attenuates Salt-Induced Renal Injury Through the Inhibition of Inflammation and Rac1-MR Activation. Int. J. Mol. Sci. 21, 4433. doi: 10.3390/ijms21124433
Jiang J., Zhou Z., Dong Y., Gao S., Sun H., Chen Z., et al. (2017). Comparative Analysis of Immunocompetence Between Females and Males in the Sea Cucumber Apostichopus Japonicus. Fish. Shellfish. Immunol. 63, 438–443. doi: 10.1016/j.fsi.2017.02.038
Jones P., Binns D., Chang H.-Y., Fraser M., Li W., McAnulla C., et al. (2014). InterProScan 5: Genome-Scale Protein Function Classification. Bioinformatics 30, 1236–1240. doi: 10.1093/bioinformatics/btu031
Kadiene E. U., Ouddane B., Gong H. Y., Kim M. S., Lee J. S., Pan Y. J., et al. (2020). Differential Gene Expression Profile of Male and Female Copepods in Response to Cadmium Exposure. Ecotoxicol. Environ. Saf. 204, 111048. doi: 10.1016/j.ecoenv.2020.111048
Kim M., Rhee J. S., Kim T. H., Lee J. S., Choi A. Y., Choi B. S., et al. (2017). Alternative Splicing Profile and Sex-Preferential Gene Expression in the Female and Male Pacific Abalone Haliotis Discus Hannai. Genes 8, 99. doi: 10.3390/genes8030099
Klaoudatos D., Tziantziou L., Lolas A., Neofitou N., Vafidis D. (2022). Population Characteristics of the Upper Infralittoral Sea Urchin Arbacia Lixula (Linnaeus 1758) in Eastern Mediterranean (Central Greece): An Indicator Species for Coastal Water Quality. J. Mar. Sci. Eng. 10, 395. doi: 10.3390/jmse10030395
Krause-Jensen D., Duarte C. M. (2016). Substantial Role of Macroalgae in Marine Carbon Sequestration. Nat. Geosci. 9, 737–742. doi: 10.1038/ngeo2790
Langmead B., Salzberg S. L. (2012). Fast Gapped-Read Alignment With Bowtie 2. Nat. Methods 9, 357–359. doi: 10.1038/nmeth.1923
Li B., Dewey C. N. (2011). RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics 12, 1−16
Li N., Arief N., Edmands S. (2019). Effects of Oxidative Stress on Sex-Specific Gene Expression in the Copepod Tigriopus Californicus Revealed by Single Individual RNA-Seq. Comp. Biochem. Physiol. Part D. Genomics Proteomics 31, 100608. doi: 10.1016/j.cbd.2019.100608
Ling S. D., Johnson C. R., Frusher S. D., Ridgway K. R. (2009). Overfishing Reduces Resilience of Kelp Beds to Climate-Driven Catastrophic Phase Shift. Proc. Natl. Acad. Sci. U.S.A. 106, 22341–22345. doi: 10.1073/pnas.0907529106
Ling S. D., Scheibling R. E., Rassweiler A., Johnson C. R., Shears N., Connell S. D., et al. (2015). Global Regime Shift Dynamics of Catastrophic Sea Urchin Overgrazing. Philos. Trans. R. Soc Lond. B. Biol. Sci. 370, 20130269. doi: 10.1098/rstb.2013.0269
Liu Z., Li M., Yi Q., Wang L., Song L. (2018). The Neuroendocrine-Immune Regulation in Response to Environmental Stress in Marine Bivalves. Front. Physiol. 9. doi: 10.3389/fphys.2018.01456
Love M. I., Huber W., Anders S. (2014). Moderated Estimation of Fold Change and Dispersion for RNA-Seq Data With Deseq2. Genome Biol. 15, 550. doi: 10.1186/s13059-014-0550-8
Maier T., Güell M., Serrano L. (2009). Correlation of mRNA and Protein in Complex Biological Samples. FEBS Lett. 583, 3966–3973. doi: 10.1016/j.febslet.2009.10.036
Maisano M., Cappello T., Catanese E., Vitale V., Natalotto A., Giannetto A., et al. (2015). Developmental Abnormalities and Neurotoxicological Effects of CuO NPs on the Black Sea Urchin Arbacia Lixula by Embryotoxicity Assay. Mar. Environ. Res. 111, 121–127. doi: 10.1016/j.marenvres.2015.05.010
Mann K., Poustka A. J., Mann M. (2008). The Sea Urchin (Strongylocentrotus Purpuratus) Test and Spine Proteomes. Proteome Sci. 6, 1–10. doi: 10.1186/1477-5956-6-22
Matozzo V., Marin M. G. (2010). First Evidence of Gender-Related Differences in Immune Parameters of the Clam Ruditapes Philippinarum (Mollusca, Bivalvia). Mar. Biol. 157, 1181–1189. doi: 10.1007/s00227-010-1398-4
Matranga V., Toia G., Bonaventura R., Müller W. E. (2000). Cellular and Biochemical Responses to Environmental and Experimentally Induced Stress in Sea Urchin Coelomocytes. Cell Stress Chaperones. 5, 113. doi: 10.1379/1466-1268(2000)005<0113:cabrte>2.0.co;2
Micheli F., Benedetti-Cecchi L., Gambaccini S., Bertocci I., Borsini C., Osio G. C., et al. (2005). Cascading Human Impacts, Marine Protected Areas, and the Structure of Mediterranean Reef Assemblages. Ecol. Monogr. 75, 81–102. doi: 10.1890/03-4058
Miller Y. I., Chang M. K., Binder C. J., Shaw P. X., Witztum J. L. (2003). Oxidized Low Density Lipoprotein and Innate Immune Receptors. Curr. Opin. Lipidol. 14, 437–445. doi: 10.1097/00041433-200310000-00004
Nguyen T. V., Alfaro A. C., Merien F., Young T., Grandiosa R. (2018). Metabolic and Immunological Responses of Male and Female New Zealand Greenshell™ Mussels (Perna Canaliculus) Infected With Vibrio Sp. J. Invertebr. Pathol. 157, 80–89. doi: 10.1016/j.jip.2018.08.008
Palacín C., Giribet G., Carner S., Dantart L., Turon X. (1998). Low Densities of Sea Urchins Influence the Structure of Algal Assemblages in the Western Mediterranean. J. Sea. Res. 39, 281–290. doi: 10.1016/S1385-1101(97)00061-0
Palumbi S. R., Barshis D. J., Traylor-Knowles N., Bay R. A. (2014). Mechanisms of Reef Coral Resistance to Future Climate Change. Science 344, 895–898. doi: 10.1126/science.1251336
Pérez-Portela R., Riesgo A. (2013). Optimizing Preservation Protocols to Extract High-Quality RNA From Different Tissues of Echinoderms for Next-Generation Sequencing. Mol. Ecol. Resour. 13, 884–889. doi: 10.1111/1755-0998.12122
Pérez-Portela R., Riesgo A., Wangensteen O. S., Palacín C., Turon X. (2020a). Enjoying the Warming Mediterranean: Transcriptomic Responses to Temperature Changes of a Thermophilous Keystone Species in Benthic Communities. Mol. Ecol. 29, 3299–3315. doi: 10.1111/mec.15564
Pérez-Portela R., Riesgo A., Wangensteen O. S., Palacín C., Turon X. (2020b). “Reference Transcriptomes of Coelomocytes”, Mendeley Data, V1. doi: 10.17632/5673n552yj.1
Pérez-Portela R., Turon X., Riesgo A. (2016). Characterization of the Transcriptome and Gene Expression of Four Different Tissues in the Ecologically Relevant Sea Urchin Arbacia Lixula Using RNAseq. Mol. Ecol. Resour. 16, 794–808. doi: 10.1111/1755-0998.12500
Pérez-Portela R., Wangensteen O. S., Garcia-Cisneros A., Valero-Jiménez C., Palacín C., Turon X. (2019). Spatio-Temporal Patterns of Genetic Variation in Arbacia Lixula, a Thermophilous Sea Urchin in Expansion in the Mediterranean. Heredity 122, 244–259. doi: 10.1038/s41437-018-0098-6
Petit G., Deboutteville C., Bougis P. (1950). Le Fichier Faunistique Du Laboratoire Arago. Vie. Milieu. 3, 356–360.
Privitera D., Noli M., Falugi C., Chiantore M. (2011). Benthic Assemblages and Temperature Effects on Paracentrotus Lividus and Arbacia Lixula Larvae and Settlement. J. Exp. Mar. Biol. Ecol. 407, 6–11. doi: 10.1016/j.jembe.2011.06.030
Rocha J., Yletyinen J., Biggs R., Blenckner T., Peterson G. (2015). Marine Regime Shifts: Drivers and Impacts on Ecosystems Services. Philos. Trans. R. Soc Lond. B. Biol. Sci. 370, 20130273. doi: 10.1098/rstb.2013.0273
Rotllant G., Nguyen T. V., Sbragaglia V., Rahi L., Dudley K. J., Hurwood D., et al. (2017). Sex and Tissue Specific Gene Expression Patterns Identified Following De Novo Transcriptomic Analysis of the Norway Lobster, Nephrops Norvegicus. BMC Genom. 18, 1–14. doi: 10.1186/s12864-017-3981-2
Schoville S. D., Barreto F. S., Moy G. W., Wolff A., Burton R. S. (2012). Investigating the Molecular Basis of Local Adaptation to Thermal Stress: Population Differences in Gene Expression Across the Transcriptome of the Copepod Tigriopus Californicus. BMC Evol. Biol. 12, 170. doi: 10.1186/1471-2148-12-170
Shao Y., Li C., Che Z., Zhang P., Zhang W., Duan X., et al. (2015). Cloning and Characterization of Two Lipopolysaccharide-Binding Protein/Bactericidal Permeability–Increasing Protein (LBP/BPI) Genes From the Sea Cucumber Apostichopus Japonicus With Diversified Function in Modulating ROS Production. Dev. Comp. Immunol. 52, 88–97. doi: 10.1016/j.dci.2015.04.015
Smale D. A., Burrows M. T., Moore P., O’Connor N., Hawkins S. J. (2013). Threats and Knowledge Gaps for Ecosystem Services Provided by Kelp Forests: A Northeast Atlantic Perspective. Ecol. Evol. 3, 4016–4038. doi: 10.1002/ece3.774
Smith L. C., Arizza V., Hudgell M. A. B., Barone G., Bodnar A. G., Buckley K. M., et al. (2018). “Echinodermata: The Complex Immune System in Echinoderms.”, in Advances in Comparative Immunology (Cham: Springer), 409–501. doi: 10.1007/978-3-319-76768-0_13
Supek F., Bošnjak M., Škunca N., Šmuc T. (2011). REVIGO Summarizes and Visualizes Long Lists of Gene Ontology Terms. PLoS One 6, e21800. doi: 10.1371/journal.pone.0021800
Vazzana M., Mauro M., Ceraulo M., Dioguardi M., Papale E., Mazzola S., et al. (2020). Underwater High Frequency Noise: Biological Responses in Sea Urchin Arbacia Lixula (Linnaeus 1758). Comp. Biochem. Physiol. Part A. Mol. Integr. Physiol. 242, 110650. doi: 10.1016/j.cbpa.2020.110650
Wangensteen O. S., Dupont S., Casties I., Turon X., Palacín C. (2013b). Some Like it Hot: Temperature and pH Modulate Larval Development and Settlement of the Sea Urchin Arbacia Lixula. J. Exp. Mar. Biol. Ecol. 449, 304–311. doi: 10.1016/j.jembe.2013.10.007
Wangensteen O. S., Turon X., Casso M., Palacín C. (2013a). The Reproductive Cycle of the Sea Urchin Arbacia Lixula in Northwest Mediterranean: Potential Influence of Temperature and Photoperiod. Mar. Biol. 160, 3157–3168. doi: 10.1007/s00227-013-2303-8
Wangensteen O. S., Turon X., García-Cisneros A., Recasens M., Romero J., Palacín C. (2011). A Wolf in Sheep’s Clothing: Carnivory in Dominant Sea Urchins in the Mediterranean. Mar. Ecol. Prog. Ser. 441, 117–128. doi: 10.3354/meps09359
Wangensteen O. S., Turon X., Pérez-Portela R., Palacín C. (2012). Natural or Naturalized? Phylogeography Suggests That the Abundant Sea Urchin Arbacia Lixula is a Recent Colonizer of the Mediterranean. PloS One 7, e45067. doi: 10.1371/journal.pone.0045067
Warnes G. R., Bolker B., Bonebakker L., Gentleman R., Huber W., Liaw A., et al. (2009). “Gplots: Various R Programming Tools for Plotting Data,” in R Package Version. 2:1.
Keywords: Sea urchin, transcriptomics, immune cells, coelomocytes, gene expression, sex-specificity,Arbacia lixula
Citation: Pérez-Portela R and Leiva C (2022) Sex-Specific Transcriptomic Differences in the Immune Cells of a Key Atlantic-Mediterranean Sea Urchin. Front. Mar. Sci. 9:908387. doi: 10.3389/fmars.2022.908387
Received: 30 March 2022; Accepted: 07 June 2022;
Published: 18 July 2022.
Edited by:
Andrew Stanley Mount, Clemson University, United StatesReviewed by:
Owen S. Wangensteen, UiT The Arctic University of Norway, NorwayJosé Carlos Hernández, University of La Laguna, Spain
Copyright © 2022 Pérez-Portela and Leiva. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Rocío Pérez-Portela, rocio_perez@ub.edu