Next Article in Journal
The N-Linked Glycosylation Site N191 Is Necessary for PKA Signal Transduction in Eel Follicle-Stimulating Hormone Receptor
Next Article in Special Issue
Transcriptomic and Metabolomic Analysis of Korean Pine Cell Lines with Different Somatic Embryogenic Potential
Previous Article in Journal
Using Staphylococcus aureus Cas9 to Expand the Scope of Potential Gene Targets for Genome Editing in Soybean
Previous Article in Special Issue
Uncovering Novel Features of the Pc Locus in Horn Development from Gene-Edited Holstein Cattle by RNA-Sequencing Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Gene Expression Detects the Factors Influencing the Reproductive Success and the Survival Rates of Paracentrotus lividus Offspring

1
Stazione Zoologica Anton Dohrn, Department of Ecosustainable Marine Biotechnology, Via Ammiraglio Ferdinando Acton n. 55, 80133 Napoli, Italy
2
Stazione Zoologica Anton Dohrn, Department of Ecosustainable Marine Biotechnology, Ischia Marine Centre, 80077 Napoli, Italy
3
Department of Biology, University of Naples Federico II, Via Cinthia, 80126 Naples, Italy
4
General Fisheries Commission for the Mediterranean (GFCM) Palazzo Blumenstihl Via Vittoria Colonna, 1, 00193 Rome, Italy
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2022, 23(21), 12790; https://doi.org/10.3390/ijms232112790
Submission received: 29 July 2022 / Revised: 11 October 2022 / Accepted: 14 October 2022 / Published: 24 October 2022
(This article belongs to the Special Issue 21st Anniversary of IJMS: Advances in Molecular Genetics and Genomics)

Abstract

:
The increase in the demand for Paracentrotus lividus roe, a food delicacy, causes increased pressure on its wild stocks. In this scenario, aquaculture facilities will mitigate the effects of anthropogenic pressures on the wild stocks of P. lividus. Consequently, experimental studies should be conducted to enhance techniques to improve efficient aquaculture practices for these animals. Here, we for the first time performed molecular investigations on cultured sea urchins. We aimed at understanding if maternal influences may significantly impact the life of future offspring, and how the culture conditions may impact the development and growth of cultured specimens. Our findings demonstrate that the outcomes of in vitro fertilization of P. lividus are influenced by maternal influences, but these effects are largely determined by culture conditions. In fact, twenty-three genes involved in the response to stress and skeletogenesis, whose expressions were measured by Real Time qPCR, were differently expressed in sea urchins cultured in two experimental conditions, and the results were largely modified in offspring deriving from two groups of females. The findings herein reported will be critical to develop protocols for the larval culture of the most common sea urchin, both for research and industrial production purposes for mass production.

1. Introduction

1.1. Sea Urchins and Aquaculture

Sea urchins are a resource for the scientific research, besides their key role in the ecology of Mediterranean shallow ecosystems and their fishery importance. In fact, the Mediterranean Paracentrotus lividus has both an ecological and an economic importance, and an emerging demand for fresh animals characterizes the fish markets, because their gonads (also known as “roe”) are considered a food delicacy worldwide and a substantial source of revenue. The roe consists not only of immature and mature gametes but also of germinal and connective tissues. Consequently, the number of gametes is not the only factor related to gonadal size, because the growth of nutritive phagocytes and accessory tissues also influence the size of gonads. This is important because it is possible to enhance the gonadal growth of farmed individuals by favoring the accumulation of nutrients independently from their actual reproductive power [1]. Aquaculture of echinoderms (sea urchins and sea cucumbers) is defined as “echinoculture”, but this term often refers specifically to sea urchin culture in the laboratory. Wild populations of P. lividus are presently facing changing environmental conditions and overexploitation of their stocks. This issue urgently needs action to ensure correct management of the resource. The implementation of eco-friendly and efficient aquaculture systems could help in overcoming the overexploitation of marine resources, including sea urchins. In addition, sea urchins play a pivotal role in scientific research [2] because they are widely used for ecotoxicological, physiological and embryological studies due to their easy handling in the laboratory and to the transparency of their embryos, ensuring that the first developmental stages can be easily monitored [3,4,5]. Sea urchin larvae are also employed as live feed for other cultured organisms, such as fishes and filter feeders (e.g., shellfish larvae). Consequently, they represent a valuable resource for the culture of other economically-relevant species [6].

1.2. Factors Influencing Larval Development

Several factors, such as diet composition, larval density and quality of the culture environment, may deeply impact the reproductive success of this species [7,8,9,10]. Variable feed abundance, feeding frequencies and species of live microalgae used for larval diets can severely influence larval growth, survival and metamorphic success [8,9,11,12] according to the cultured species of sea urchin. In addition, the nutritional quality of microalgae can vary depending on culture factors influencing their cell size, digestibility and biochemical composition [13]. For example, Chaetoceros spp. and Tisochrysis lutea contain relatively high amounts of long-chain PUFAs, but neither are rich in both EPA and DHA [14], and in their research were found to decreasein the absence of Rhodomonas lens. Castilla-Gavillán et al. [10] demonstrated that P. lividus larvae fed on Rhodomonas sp. contained higher total lipid content than those fed on other microalgae. Similar dynamics could have positive implications for the development of larval echinoderms. Previous studies [15] indicated that Dunaliella sp. may be sufficient to sustain the complete development of P. lividus up to metamorphosis, although small additions of other feeds (e.g., R. lens) may represent a useful improvement towards complete success of the larviculture. In contrast, some microalgae may have deleterious effects on the growth of post-larvae [16].
However, maternal and zygotic factors also influence and control the early development of sea urchins. Maternal factors include messenger RNA and proteins expressed during oogenesis. mRNA is also necessary before activation of the embryonic genome for early development, as they are involved in the regulation of metabolism, cell cycle and development [17]. When genome activation occurs in embryos, zygotic factors begin to be expressed. The first known developmental genes are expressed at the reaching of 16-cell stage; before this stage any developmental process is completely driven by maternal factors [18,19]. Besides their genetic contribution, mother sea urchins control the phenotypic development of their offspring in response to environmental conditions [20]. These maternal influences, which result in the combined effect of maternal phenotype and genotype [21], can have significant effects on the fitness and performance of the offspring [22,23,24]. In addition, they can influence the sea urchin population ecology [21,25]. However, factors experienced by sea urchin females during early development can also affect the phenotypes and fitness of their offspring.
In addition, the size of culture vessels and the density of larvae may represent key factors in determining rates of survival of larvae up to metamorphosis. An experiment conducted in small-scale 1.3-L culture vessels revealed an inverse relationship between larval density of a sea urchin (Diadema antillarum) and growth [26]. Previous authors obtained variable rates of survival by approaching different larval-rearing systems [27]. Different studies were characterized by variable feeding protocols, fertilization techniques and rearing devices [28]. In particular, previous research has indicated culture experiences within tanks of different volumes, besides the influence of additional rearing devices (e.g., aeration systems, palettes, procedures, etc.). Hence, we aimed to investigate if tanks of the same shape but different sizes might influence the survival of larvae. In fact, keeping constant the larval density and the concentration of microalgae feed, tanks of different volumes are characterized by variable hydrodynamic conditions. For example, larger tanks need faster bubbling to ensure sufficient gas exchanges, when compared to smaller tanks, but aeration may interfere with the limited swimming capabilities of planktonic larvae, pushing them into some sectors on the bottom. However, an appropriate determination of the most favorable stocking density within a production-oriented recirculating aquaculture system (RAS) has not been performed. Consequently, the size of culture vessels, the quality of larval nutrition and the stocking density remain extremely relevant for the larviculture of sea urchins, especially when considering the feasibility of production for fisheries, restoration and scientific research. For these reasons, we tested the effect of the tank size in order to detect any influence on survival rates during the early development of these sea urchins.
We had two distinct aims in our manuscript: (i) we aimed at understanding if maternal influences may significantly impact the life of future offspring; (ii) we aimed at detecting how and if slight variations in the culture conditions may impact the development and growth of cultured specimens. In addition, we aimed at evaluating the best production sizes, in terms of tank volumes, to ensure correct growth and higher survival rates of P. lividus embryos. On the whole, we aimed at understanding how genic resources interact with environmental stressors to ensure maximum reproductive performance in a marine invertebrate. Furthermore, and for the first time, we tested the effect of these production units by performing molecular analyses (Real-Time qPCR) to check the expression levels of several genes involved in stress responses and skeletogenesis. In particular, we followed the variation of expression levels of twenty-three genes, which were first isolated in sea urchin embryos in response to different stressors [3,5,29,30,31,32,33]. In this study, gene expression levels are not responses related to a stimulus but must be intended as higher or lower levels of gene expression.

2. Results

2.1. Survival Rates

Even if the survival rates of larvae deriving from females A-B in the first three weeks of growth were not significantly different from those deriving from females C–D, when analyzed by means of paired t-tests (Table 1):
The performances of the two sets of larvae exhibited different patterns of development (Figure 1).
In particular, pool A-B exhibited higher survival rates in the smaller tanks (Figure 1A), completing the settlement after 17 days. In contrast, pool C-D exhibited lower survival rates during the first two weeks and the number of swimming larvae was strongly reduced after two weeks, due to higher mortality. In fact, in the last two weeks, larvae were absent in both groups of females considered. These differences were less evident in larger tanks (Figure 1B), where the number of larvae continuously decreased in the first 20 days, exhibiting a low number of settled individuals in both A-B and C-D groups (Figure 1B). The larval stage progression was initially similar in the two pools and the embryonic stages developed synchronously through the pluteus stage. However, from the onset of feeding through the larval stages, the development became increasingly less synchronous between pools A-B vs. C-D, and even among tanks of different sizes. After the first week, larvae from the pool A-B exhibited faster development, showed dark-green guts and after 10–12 days reached a six-arm stage, while larvae from the pool C-D exhibited slower development, about half of them arrested their development at a four-arm stage, and were characterized by partially-empty guts. In particular, larvae cultured in smaller tanks reached the advanced rudiment stage after about two weeks, while larvae in larger tanks showed lower mobility, sank to the bottom more often and their guts exhibited a pale color.
To evaluate the effect of maternal influences and the effect of tank size on the survival of larvae (removing the influence of settlers on the density of swimming larvae) only the first week of culture was analyzed and the slopes of linear regressions were compared (Figure 2).
This representation permits the evaluation that in smaller tanks, during the first week, an 100% survival was exhibited by larvae deriving from females A and B, whilst, in the same period, females C and D produced a decline of larvae, down to 2.6% of the initial pool. When slopes were compared, the differences in the survival patterns of the two larval pools were significant at p < 0.05. In larger tanks a continuous decrease of the swimming larvae was observed during the first week for both larval pools (Figure 2B) and the differences in the survival patterns of the two larval pools were not significant (p > 0.05).

2.2. Molecular Analyses of Embryos Deriving from Various Females

The variation of expression levels of twenty-three genes, involved in the stress response and skeletogenesis and previously analyzed [3,5,34,35,36,37,38] (see Table S3) was followed by Real Time qPCR (see Table S4).

2.2.1. Genes Involved in Stress Response

Concerning the 18 genes analyzed, the results show that the plutei deriving from female A did not have variations in expression for any of the genes analyzed. The same results were found for the plutei deriving from female B, with the only exception being CASP8 gene, which decreased its expression level. In the case of plutei deriving from female C, fifteen genes decreased their expression levels: ARF1 (−3.5), CASP8 (−5.5), caspase 3/7 (−3.5), ERCC3 (−3.9), GRHPR (−2.4), hsp60 (−2.4), hsp70 (−6.4), HIF1A (−2.0), PARP-1 (−3.0), p53 (−4.8) and 14-3-3 ε (−5.5). Only cytb gene showed an increase (2.5) of its expression level. Moreover, in the plutei deriving from female D, nine genes decreased their expression levels (CASP8, −2.7; ERCC3, −3.2; hsp60, −2.8; hsp70, −2.9; p38MAPK, −2.3; SDH, −2.7; p53, −4.4; and 14-3-3 ε, −3.4).

2.2.2. Skeletogenic Genes

None of the five genes analyzed was switched on in the plutei deriving from females A and B. Three genes were down-regulated in the plutei deriving from female C (SM30, −3.9; SM50, −4.0; uni, 4.0), whereas plutei deriving from female D showed a down-regulation of SM50 (−7.3) and an up-regulation of univin (4.0).

2.3. Network Analysis

Interactomic analysis showed that eighteen genes of the twenty-three analyzed from Real Time qPCR are connected by functional point of view: ADP-ribosylation factor 1 (ARF1), caspase-8 (CASP8), Sp-Cspe3/7L (CASPASE 3/7), cytochrome b (CYTB), DNA-methyltransferase 1 (MTase), ERCC excision repair 3 (ERCC3), glutamine synthetase (GS), glyoxylate reductase/hydroxypyruvate reductase (GRHPR), heat shock protein 56 (HSP56), heat shock protein 60 (HSP60), heat shock protein 70 (HSP70), hypoxia inducible factor 1-alpha (HIF1A), nuclear factor kappa light-chain-enhancer of activated B cells (NF-kB), Poly(ADP-ribose) polymerase 1 (PARP-1), p38 mitogen-activated protein kinase (p38 MAPK), succinate dehydrogenase (SDH), tumor protein p53 (p53). Only the skeletogenic genes 14-3-3 epsilon protein (14-3-3 ε), bone morphogenetic protein 5-7 (BMP5-7), and nectin (Nec) weren’t included in network analysis because they have no orthologous genes in humans. Moreover, GRHPR is not linked to any gene under analysis.
As shown in the network reported in Figure 3, the genes are correlated as follows: ARF1 with hsp70, p53 and PARP1; Caspase 3/7 with CASP8, hsp70, p38 mapk, PARP1; CASP8 with Caspase 3/7, p38 mapk, hsp70, p53, NF-kB, PARP1; p38 mapk with CASP8, Caspase 3/7, hsp60, hsp70, p53 and PARP1; Hsp70 with CASP8, p53, NF-kB, MTase, HIF1A, hsp60, SDH; p53 with p38 mapk, CASP8, hsp70, HIF1A, ERCC3, NF-kB, MTase, PARP1, hsp56; PARP1 with ARF1, Caspase 3/7, p38 mapk, CASP8, p53, HIF1A, NF-kB, ERCC3 and MTase; Hsp56 with p53; NF-kB with PARP1, p38 mapk, p53, hsp70, HIF1A, ERCC3, MTase; MTase with PARP1, NF-kB, p53, hsp70, ERCC3; ERCC3 with MTase, PARP1 and p53; HIF1A with NF-kB, PARP1, MTase, p53, hsp70 and SDH; SDH with hsp70, HIF1A, hsp60 and CytB; CytB with SDH and hsp60; Hsp60 with CytB, SDH, p38 mapk, hsp70, p53 and GS; GS with hsp60 and hsp70. The corresponding human orthologous genes are also reported.

3. Discussion

Investigations of other echinoderms reported that the performance of one life history stage led to (positive or negative) carryover effects on the subsequent life stages [39]. In fact, carryover effects arising within a generation can influence next generations [40,41,42]. Carryover effects can also arise across a generation, an effect known as transgenerational plasticity (TGP) [43]. Due to environmental stresses [44] experienced by one or both parents during the development of their gametes [23,41], a transgenerational response can lead to phenotypic changes in the offspring [43,45]. Various evidence indicates that long-term exposure of adult sea urchins to thermal shocks impacts the embryo development of their offspring [39,46,47]. Here we aimed at checking and eventually confirming this hypothesis regarding Paracentrotus lividus embryos and larvae.
During the whole experiment the main abiotic descriptors were kept stable and did not achieve significant variations among replicates and treatments. Maternal effects on the survival rates of larvae were definitely observed in smaller culture vessels. Larger vessels prompt a negative influence on the survival, which is likely to be superimposed onto the material influence itself. This evidence indicates that both the maternal effects and the size of the culture tanks are critical to determine the reproductive success of this species, but the influence of the latter overwhelms that of the former.
Evidently the size of the tank matters, because the two experimental conditions produced significantly different results in terms of survival rates of larvae. The effect might be explained by taking into account that circulation of water and oxygen stratification change drastically according to the shape and the size of the culture vessels. In particular, deeper tanks exhibit a lower circulation of water and a higher stratification of gases, because the rate of aeration (air bubbling) cannot be proportionally changed to avoid damage to the delicate larvae swimming closer to the bubbling tubes [48,49,50]. We might also consider the need for larvae to remain in the upper levels and to avoid strong contact with the bottom of the tank, since it may contain algal and bacterial films which produce mucilage that promotes clogging and aggregation of larvae (personal observation; [51]). Since larger tanks require higher levels of aeration to guarantee sufficient gas exchanges, the fluxes produced by the air stone also promote the creation of sink areas where larvae aggregate close to the bottom. In contrast, the aeration may be very weak in the smaller tanks, permitting the larvae to swim towards the surface without interference. In fact, in smaller containers most larvae were continuously visible under the surface of the water, while in larger tanks they were actively transported by the currents and could accumulate in areas with lower perturbations, close to the bottom (personal observation). Attempts to reduce this effect by decreasing the flux of the air produced lower gas exchanges which were insufficient to guarantee consistency of the abiotic factors (including the saturation of oxygen in the seawater) that represented a pre-condition to evaluate other sources of stress [49,52].
As reported, the density of larvae deriving from females C-D exhibited a decrease in smaller tanks, as opposed to larvae deriving from females A-B, which performed best. In contrast, both groups of larvae showed mortality rates higher than 60% in the first week when cultured in larger tanks, and the differences were not significant in this condition. This may lead to the conclusion that the first factor determining survival rate is represented by the environmental conditions and that, in an optimal environment, maternal influences become determinant. The overlap of the growth conditions to the maternal influence may also explain the interannual differences in larval recruitment of this species observed in the field [53].
Molecular analyses were performed to understand if there were differences in the levels of gene expression associated with offspring deriving from different P. lividus females (see Figure 4). Molecular data strongly supported the low success of embryonic development for the plutei deriving from females 3 and 4. In fact, plutei deriving from females A and B did not show significant variations in their gene expressions.
In contrast, plutei deriving from females C and D exhibited a remarkable switch on in almost all the genes analyzed that are involved in stress response and in skeletogenic process, indicating negative maternal influences, as showed by the results of Real Time qPCR. In fact, the down-regulation of genes in the plutei deriving from P. lividus females C and D was evident compared to those deriving from females A and B. These effects in marine organisms are not new [23], because they have been widely investigated in various sea urchin species [54,55,56], but to the best of our knowledge such relationships have never been investigated in P. lividus [17].
Interactomic analyses showed that stress genes were functionally correlated among themselves. GRHPR was the only gene that appeared to have no correlation with others. Previous investigations indicated that this gene is correlated to HSPA4 and GLUL through CAT [38], which was also analyzed in this research.
Changes of gene expression induced by environmental stressors on somatic cells can affect the physiology of exposed individuals. As a consequence, some alterations can be propagated to subsequent generations through the germline, occurring during the gametogenesis [57,58,59]. In particular, environmental stresses experienced by one or both parents may lead to aberrant patterns in F1, which can become evident from the early developmental stages [56,60]. It is known that a network of regulatory genes, the so called “defensome” [3,30,61], regulate the mechanisms of specification and differentiation during embryonic development.
Our findings show how changes in the levels of gene expression may be used as an early indicator of stressful conditions for sea urchins, and more generally for marine invertebrates and the marine environment itself. Furthermore, these findings demonstrate that changes in gene expression can help understand the health state of adult sea urchins. The identification of key genes and the molecular pathways in which they are involved represent fundamental tools for understanding how marine organisms attempt to protect themselves against toxicants in order to avoid deleterious consequences and irreversible damage. Furthermore, changes in the response to environmental stressors may contribute to an adaptive survival advantage of organisms through beneficial gene expression [39].
The findings herein demonstrated will be useful for the culture of the most common Mediterranean sea urchin, both for research purposes and for the needs of aquaculture. In fact, despite several established methods of commercially producing echinoderms from gametes [62,63], a scalable hatchery process for P. lividus has not yet been established. Attempts to culture this species over the last decades have been met with varying degrees of success, and speculative factors including environmental toxins, water quality and general culture methods precluded reliable development through metamorphosis [26]. This study was intended to develop P. lividus larval culture protocols within a laboratory or an industrial production plant capable of scaled production for aquaculture, restoration or scientific research. The culture conditions demonstrated are adequate for the larval development of P. lividus, but variables of interest, i.e., maternal influences and size of the culture vessels, are demonstrated to be critical factors to reach sufficient survival rates and complete development of this species. However, further investigations aimed at improving yields are necessary, as well as to understand how larval nutrition affects settlement and post-settlement success in order to further improve commercial production viability.

4. Materials and Method

4.1. Gamete Collection and Embryo Culture

Adult sea urchins were collected by scuba divers around the island of Procida (Italy) in February 2021, at 10 m depth. Collected specimens were transported in an insulated box to the laboratory and maintained in tanks with circulating sea water for ten days until testing. Sea urchins were injected with 1 mL of 2 M KCl through the peribuccal membrane to stimulate the emission of gametes. After the production of gametes, the specimens were restored in aerated recirculated tanks and released after a few days into the same area of collection. Eggs were collected in a glass dish, washed with filtered sea water (FSW) and kept in FSW until use. Concentrated sperm was collected from males by means of a plastic pipette and kept undiluted at 4 °C until use.
Eggs were fertilized and were kept at 20 °C in a controlled temperature chamber on a 12 h–12 h light–dark cycle. Forty-eight hours post-fertilization the embryos reached the pluteus stage, and these larvae were distributed in equal concentrations (1.5 plutei mL−1) in triplicate tanks of two different volumes: (i) small scale, in 2 L beakers; (ii) larger scale, in 10 L tanks. Embryos deriving from individual females were paired two by two and kept separate, aiming at detecting any variation due to maternal influences and if the initial quality of the gametes influenced the larval development. In particular, the eggs deriving from four females were randomly pooled into two groups (females A-B and females C-D) prior to being fertilized with a pool of sperms deriving from five males. Four females were randomly chosen because TR-qPCR technique is useful to detect individual variations in the gene expression, avoiding the flattening of results due to variable gene expressions detected in larger populations.
Larvae were fed on Dunaliella tertiolecta (5000 cells ml−1; [15]) along with a composed algal food (SHG Snow Reef, Super High Group, Ovada, Italy; www.superhigroup.com accessed on 28 July 2022) containing yeast, fish oils, Chlorella and various HUFAs reproducing the composition of natural marine snow in a suspension of particles sizing 5–10 μm [15]. Larvae were filtered every two days onto a 60 μm mesh sieve and transferred to a renewed culture medium using a Pasteur pipette. At the same time intervals, two 50 mL samples were collected to evaluate and record the survival rates, check the larval stage progression, the overall health conditions (mobility, malformations) and the gut fullness. During the experiment, seawater samples were collected biweekly and nitrites and phosphates were measured (using a photometer AL 450, Aqualytic, Alkimia Srl), as well as salinity and temperature (see Supplementary Table S1 in the Supplementary Materials).

4.2. Statistical Analyses

The survival rates of larvae deriving from two sets of females (pooling of individuals A-B vs. C-D) were calculated using GraphPad Prism 8 software (GraphPad Software, San Diego, CA, USA, www.graphpad.com, accessed 1 September 2021) using the following equation for one phase decay:
Y = ( Y 0 Plateau ) exp   ( K X ) + Plateau
where: Y0 is the Y value when X (time) is zero; Plateau is the Y value at infinite times; K is the rate constant, expressed in reciprocal of the X-axis time units. The survival rates obtained with A-B vs. C-D larvae, as well as those obtained in smaller vs. larger tanks, were compared and two-sided p-values were obtained by paired t-tests.
The slopes obtained from different individuals (A-B vs. C-D) in different types of culture tanks (smaller vs. larger) were compared by linear regression on the data limited to the first week, to avoid any influence of settlement (older larvae settle and the number of swimming individuals is reduced) and to determine the slope of a best-fit line. A value of p < 0.05 was chosen as a threshold level for significance. All data have been analyzed with the aid of GraphPad Prism 8.0 software.

4.3. RNA Extraction and cDNA Synthesis

About 5000 fertilized eggs from each of the four sea urchin females were collected at the pluteus stage, corresponding to 48 h post-fertilization (hpf). These samples were centrifuged at 3500 relative centrifugal force for 15 min in a swing-out rotor at 4 °C. Embryos were then placed in 60 microliters of the RNA-later (Qiagen, Hilden, Germany) and then frozen in liquid nitrogen. Samples were kept at −80 °C until use. Total RNA was extracted using Aurum Total RNA Mini Kit (Bio–Rad, Hercules, CA, USA), according to the manufacturer’s instructions. The amount of total RNA extracted was estimated by the absorbance at 260 nm and the purity by 260/280 and 260/230 nm ratios, using a NanoDrop spectrophotometer (ND–1000 UV–VIS Spectrophotometer; NanoDrop Technologies, Wilmington, DE, USA). The integrity of RNA was evaluated by observing the rRNA subunits (28S and 18S) on agarose gel electrophoresis. For each sample, 1000 ng of total RNA was retrotranscribed with an iScript cDNA Synthesis kit (Bio–Rad, Milan, Italy), following the manufacturer’s instructions.

4.4. Real Time qPCR Experiments

The variations in the expression of twenty-three genes were followed by Real Time qPCR (see Supplementary Table S2 in the Supplementary Materials for their biological functions and Supplementary Table S3 for the sequence of the primers) in sea urchin embryos deriving from each of the four females collected at the pluteus stage, reached at 48 hpf. cDNA (1 μL) was used as a template in a reaction containing a final concentration of 0.3 mM for each primer and 1× FastStart SYBR Green master mix (total volume of 10 μL) (Applied Biosystems, Monza, Italy). PCR amplifications were performed in a ViiATM7 Real Time PCR System (Applied Biosystems, Monza, Italy) thermal cycler using the following thermal profile: 95 °C for 10 min, one cycle for cDNA denaturation; 95 °C for 15 s and 60 °C for 1 min, 40 cycles for amplification; one cycle for final elongation at 72 °C for 5 min; one cycle for melting curve analysis (from 60 °C to 95 °C) to verify the presence of a single product. Each assay included a no-template control for each primer pair. To capture intra-assay variability, all Real-Time qPCR reactions were carried out in triplicate. Fluorescence was measured using ViiATM7 software (Applied Biosystems, Monza, Italy). For all Real-Time qPCR experiments the data of each cDNA sample were normalized with the mRNA level of the ubiquitin [5] and 18S rRNA [64] genes, used as housekeeping genes. The expression of each gene was analyzed and normalized against the internal control by REST program (Relative Expression Software Tool), based on the Pfaffl method, reporting the expression values of the genes of interest with respect to the housekeeping genes, using each time one of four different females as the reference group [65,66]. Relative expression ratios above two cycles were considered significant. Experiments are means between duplicates. Statistical analysis was performed using GraphPad Prism version 8.00 for Windows (GraphPad Software, San Diego, CA, USA).
The network analysis on genes under analysis was performed by STRING [67], which aims to identify relationships on the basis of associated functions and data mining from experimental studies reported in the literature. Since sea urchin genes are not annotated in the STRING database, human orthologues were used to search for P. lividus genes.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms232112790/s1.

Author Contributions

S.F. performed the experimental work and the statistical analyses and produced the first draft of the manuscript. F.G. performed and supervised the experimental work in Procida. R.E. contributed to molecular experiments and aided the production of the first draft and its revision and contributed statistical analyses and molecular investigations. B.P. and M.G. contributed to the experimental work in Procida. A.D.C. planned the research and supervised the production of the manuscript. M.C. planned the research, performed molecular tests, supervised the students and contributed to the production of the first draft and its final revision. V.Z. planned the research, revised the manuscript and supervised the experimental work in Procida. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We acknowledge the private firm Echinoidea in Procida Island (Gulf of Naples), managed by Chiara and Michele Trapanese, where in vitro fertilization of the sea urchin Paracentrotus lividus experiments were performed, including larval culture at various production scales. Thanks to the technical help of Mario Loffredo e Domenico Mattera. Francesca Glaviano was supported by a PhD (PhD in Biology, University of Naples Federico II) fellowship at the Stazione Zoologica Anton Dohrn. Roberta Esposito was supported by a PhD (PhD in Biology, University of Naples Federico II) fellowship funded by the Photosynthesis 2.0 project of the Stazione Zoologica Anton Dohrn.

Conflicts 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.

References

  1. Böttger, S.A.; Devin, M.G.; Walker, C.W. Suspension of Annual Gametogenesis in North American Green Sea Urchins (Strongylocentrotus Droebachiensis) Experiencing Invariant Photoperiod—Applications for Land-Based Aquaculture. Aquaculture 2006, 261, 1422–1431. [Google Scholar] [CrossRef]
  2. Adonin, L.; Drozdov, A.; Barlev, N.A. Sea Urchin as a Universal Model for Studies of Gene Networks. Front. Genet. 2021, 11, 627259. [Google Scholar] [CrossRef] [PubMed]
  3. Marrone, V.; Piscopo, M.; Romano, G.; Ianora, A.; Palumbo, A.; Costantini, M. Defensome against Toxic Diatom Aldehydes in the Sea Urchin Paracentrotus Lividus. PLoS ONE 2012, 7, e31750. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Ruocco, N.; Costantini, M.; Santella, L. New Insights into Negative Effects of Lithium on Sea Urchin Paracentrotus Lividus Embryos. Sci. Rep. 2016, 6, 32157. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Romano, G.; Costantini, M.; Buttino, I.; Ianora, A.; Palumbo, A. Nitric Oxide Mediates the Stress Response Induced by Diatom Aldehydes in the Sea Urchin Paracentrotus Lividus. PLoS ONE 2011, 6, e25980. [Google Scholar] [CrossRef] [Green Version]
  6. Luis, O.; Delgado, F.; Gago, J. Year-Round Captive Spawning Performance of the Sea Urchin Paracentrotus Lividus: Relevance for the Use of Its Larvae as Live Feed. Aquat. Living Resour. 2005, 18, 45–54. [Google Scholar] [CrossRef]
  7. Hammer, H.; Hammer, B.; Watts, S.; Lawrence, A.; Lawrence, J. The Effect of Dietary Protein and Carbohydrate Concentration on the Biochemical Composition and Gametogenic Condition of the Sea Urchin Lytechinus Variegatus. J. Exp. Mar. Biol. Ecol. 2006, 334, 109–121. [Google Scholar] [CrossRef]
  8. Carboni, S.; Vignier, J.; Chiantore, M.; Tocher, D.R.; Migaud, H. Effects of Dietary Microalgae on Growth, Survival and Fatty Acid Composition of Sea Urchin Paracentrotus Lividus throughout Larval Development. Aquaculture 2012, 324, 250–258. [Google Scholar] [CrossRef]
  9. Cárcamo, P.F.; Candia, A.I.; Chaparro, O.R. Larval Development and Metamorphosis in the Sea Urchin Loxechinus Albus (Echinodermata: Echinoidea): Effects of Diet Type and Feeding Frequency. Aquaculture 2005, 249, 375–386. [Google Scholar] [CrossRef]
  10. Castilla-Gavilán, M.; Buzin, F.; Cognie, B.; Dumay, J.; Turpin, V.; Decottignies, P. Optimising Microalgae Diets in Sea Urchin Paracentrotus Lividus Larviculture to Promote Aquaculture Diversification. Aquaculture 2018, 490, 251–259. [Google Scholar] [CrossRef]
  11. Kelly, M.S.; Hunter, A.J.; Scholfield, C.L.; McKenzie, J.D. Morphology and Survivorship of Larval Psammechinus Miliaris (Gmelin) (Echinodermata: Echinoidea) in Response to Varying Food Quantity and Quality. Aquaculture 2000, 183, 223–240. [Google Scholar] [CrossRef]
  12. Liu, H.; Kelly, M.S.; Cook, E.J.; Black, K.; Orr, H.; Zhu, J.X.; Dong, S.L. The Effect of Diet Type on Growth and Fatty Acid Composition of the Sea Urchin Larvae, II. Psammechinus Miliaris (Gmelin). Aquaculture 2007, 264, 263–278. [Google Scholar] [CrossRef]
  13. Guedes, A.C.; Malcata, F.X. Nutritional Value and Uses of Microalgae in Aquaculture. Aquaculture 2012, 10, 59–78. [Google Scholar]
  14. Brown, M.R.; Jeffrey, S.W.; Volkman, J.K.; Dunstan, G.A. Nutritional Properties of Microalgae for Mariculture. Aquaculture 1997, 151, 315–331. [Google Scholar] [CrossRef]
  15. Zupo, V.; Glaviano, F.; Caramiello, D.; Mutalipassi, M. Effect of Five Benthic Diatoms on the Survival and Development of Paracentrotus Lividus Post-Larvae in the Laboratory. Aquaculture 2018, 495, 13–20. [Google Scholar] [CrossRef]
  16. Ruocco, N.; Costantini, S.; Zupo, V.; Lauritano, C.; Caramiello, D.; Ianora, A.; Budillon, A.; Romano, G.; Nuzzo, G.; D’Ippolito, G. Toxigenic Effects of Two Benthic Diatoms upon Grazing Activity of the Sea Urchin: Morphological, Metabolomic and de Novo Transcriptomic Analysis. Sci. Rep. 2018, 8, 5622. [Google Scholar] [CrossRef] [Green Version]
  17. Kipryushina, Y.O.; Yakovlev, K.V. Maternal Control of Early Patterning in Sea Urchin Embryos. Differentiation 2020, 113, 28–37. [Google Scholar] [CrossRef]
  18. Tu, Q.; Brown, C.T.; Davidson, E.H.; Oliveri, P. Sea Urchin Forkhead Gene Family: Phylogeny and Embryonic Expression. Dev. Biol. 2006, 300, 49–62. [Google Scholar] [CrossRef] [Green Version]
  19. Yaguchi, S.; Yaguchi, J.; Angerer, R.C.; Angerer, L.M. A Wnt-FoxQ2-Nodal Pathway Links Primary and Secondary Axis Specification in Sea Urchin Embryos. Dev. Cell 2008, 14, 97–107. [Google Scholar] [CrossRef] [Green Version]
  20. Burton, T.; McKelvey, S.; Stewart, D.C.; Armstrong, J.D.; Metcalfe, N.B. Early Maternal Experience Shapes Offspring Performance in the Wild. Ecology 2013, 94, 618–626. [Google Scholar] [CrossRef] [Green Version]
  21. Venturelli, P.A.; Murphy, C.A.; Shuter, B.J.; Johnston, T.A.; van Coeverden de Groot, P.J.; Boag, P.T.; Casselman, J.M.; Montgomerie, R.; Wiegand, M.D.; Leggett, W.C. Maternal Influences on Population Dynamics: Evidence from an Exploited Freshwater Fish. Ecology 2010, 91, 2003–2012. [Google Scholar] [CrossRef] [PubMed]
  22. Green, B.S. Maternal Effects in Fish Populations. Adv. Mar. Biol. 2008, 54, 1–105. [Google Scholar] [PubMed]
  23. Marshall, D.J. Transgenerational Plasticity in the Sea: Context-dependent Maternal Effects across the Life History. Ecology 2008, 89, 418–427. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Uller, T. Developmental Plasticity and the Evolution of Parental Effects. Trends Ecol. Evol. 2008, 23, 432–438. [Google Scholar] [CrossRef] [PubMed]
  25. Benton, T.G.; St. Clair, J.J.H.; Plaistow, S.J. Maternal Effects Mediated by Maternal Age: From Life Histories to Population Dynamics. J. Anim. Ecol. 2008, 77, 1038–1046. [Google Scholar] [CrossRef]
  26. Leber, K.M.; Lorenzen, K.; Main, K.L.; Moe, M., Jr.; Vaughan, D.; Capo, T.; Bardales, A.; Gillette, P. Developing Restoration Methods to Aid in Recovery of a Key Herbivore, Diadema Antillarum, on Florida Coral Reefs-2008/2009 Final Report (April 1, 2008 to March 31, 2009); Mote Marine Laboratory: Sarasota, FL, USA; University of Miami: Coral Gables, FL, USA, 2009. [Google Scholar]
  27. Grosjean, P. Growth Model of the Reared Sea Urchin Paracentrotus Lividus (Lamarck, 1816). Ph.D. Thesis, Université de Mons, Mons, Belgium, 2001. [Google Scholar]
  28. Grosjean, P.; Spirlet, C.; Jangoux, M. Experimental Study of Growth in the Echinoid Paracentrotus Lividus (Lamarck, 1816) (Echinodermata). J. Exp. Mar. Biol. Ecol. 1996, 201, 173–184. [Google Scholar] [CrossRef]
  29. Varrella, S.; Romano, G.; Ianora, A.; Bentley, M.G.; Ruocco, N.; Costantini, M. Molecular Response to Toxic Diatom-Derived Aldehydes in the Sea Urchin Paracentrotus Lividus. Mar. Drugs 2014, 12, 2089–2113. [Google Scholar] [CrossRef] [Green Version]
  30. Varrella, S.; Romano, G.; Ruocco, N.; Ianora, A.; Bentley, M.G.; Costantini, M. First Morphological and Molecular Evidence of the Negative Impact of Diatom-Derived Hydroxyacids on the Sea Urchin Paracentrotus Lividus. Toxicol. Sci. 2016, 151, 419–433. [Google Scholar] [CrossRef] [Green Version]
  31. Ruocco, N.; Costantini, S.; Zupo, V.; Romano, G.; Ianora, A.; Fontana, A.; Costantini, M. High-Quality RNA Extraction from the Sea Urchin Paracentrotus Lividus Embryos. PLoS ONE 2017, 12, e0172171. [Google Scholar] [CrossRef]
  32. Morroni, L.; Sartori, D.; Costantini, M.; Genovesi, L.; Magliocco, T.; Ruocco, N.; Buttino, I. First Molecular Evidence of the Toxicogenetic Effects of Copper on Sea Urchin Paracentrotus Lividus Embryo Development. Water Res. 2019, 160, 415–423. [Google Scholar] [CrossRef]
  33. Albarano, L.; Zupo, V.; Guida, M.; Libralato, G.; Caramiello, D.; Ruocco, N.; Costantini, M. PAHs and PCBs Affect Functionally Intercorrelated Genes in the Sea Urchin Paracentrotus Lividus Embryos. Int. J. Mol. Sci. 2021, 22, 12498. [Google Scholar] [CrossRef] [PubMed]
  34. Russo, R.; Pinsino, A.; Costa, C.; Bonaventura, R.; Matranga, V.; Zito, F. The Newly Characterized Pl-jun Is Specifically Expressed in Skeletogenic Cells of the Paracentrotus Lividus Sea Urchin Embryo. FEBS J. 2014, 281, 3828–3843. [Google Scholar] [CrossRef] [PubMed]
  35. Varrella, S.; Romano, G.; Costantini, S.; Ruocco, N.; Ianora, A.; Bentley, M.G.; Costantini, M. Toxic Diatom Aldehydes Affect Defence Gene Networks in Sea Urchins. PLoS ONE 2016, 11, e0149734. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Ruocco, N.; Varrella, S.; Romano, G.; Ianora, A.; Bentley, M.G.; Somma, D.; Leonardi, A.; Mellone, S.; Zuppa, A.; Costantini, M. Diatom-Derived Oxylipins Induce Cell Death in Sea Urchin Embryos Activating Caspase-8 and Caspase 3/7. Aquat. Toxicol. 2016, 176, 128–140. [Google Scholar] [CrossRef]
  37. Ruocco, N.; Annunziata, C.; Ianora, A.; Libralato, G.; Manfra, L.; Costantini, S.; Costantini, M. Toxicity of Diatom-Derived Polyunsaturated Aldehyde Mixtures on Sea Urchin Paracentrotus Lividus Development. Sci. Rep. 2019, 9, 517. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Esposito, R.; Ruocco, N.; Albarano, L.; Ianora, A.; Manfra, L.; Libralato, G.; Costantini, M. Combined Effects of Diatom-Derived Oxylipins on the Sea Urchin Paracentrotus Lividus. Int. J. Mol. Sci. 2020, 21, 719. [Google Scholar] [CrossRef] [Green Version]
  39. Ross, P.M.; Parker, L.; Byrne, M. Transgenerational Responses of Molluscs and Echinoderms to Changing Ocean Conditions. ICES J. Mar. Sci. 2016, 73, 537–549. [Google Scholar] [CrossRef] [Green Version]
  40. Byrne, M.; Prowse, T.A.A.; Sewell, M.A.; Dworjanyn, S.; Williamson, J.E.; Vaïtilingon, D. Maternal Provisioning for Larvae and Larval Provisioning for Juveniles in the Toxopneustid Sea Urchin Tripneustes Gratilla. Mar. Biol. 2008, 155, 473–482. [Google Scholar] [CrossRef]
  41. Kovalchuk, I. Transgenerational Epigenetic Inheritance in Animals. Front. Genet. 2012, 3, 76. [Google Scholar] [CrossRef] [Green Version]
  42. Munday, P.L.; Warner, R.R.; Monro, K.; Pandolfi, J.M.; Marshall, D.J. Predicting Evolutionary Responses to Climate Change in the Sea. Ecol. Lett. 2013, 16, 1488–1500. [Google Scholar] [CrossRef]
  43. Byrne, M.; Selvakumaraswamy, P.; Ho, M.A.; Woolsey, E.; Nguyen, H.D. Sea Urchin Development in a Global Change Hotspot, Potential for Southerly Migration of Thermotolerant Propagules. Deep. Sea Res. Part II Top. Stud. Oceanogr. 2011, 58, 712–719. [Google Scholar]
  44. Shama, L.N.S.; Strobel, A.; Mark, F.C.; Wegner, K.M. Transgenerational Plasticity in Marine Sticklebacks: Maternal Effects Mediate Impacts of a Warming Ocean. Funct. Ecol. 2014, 28, 1482–1493. [Google Scholar] [CrossRef]
  45. Hamdoun, A.; Epel, D. Embryo Stability and Vulnerability in an Always Changing World. Proc. Natl. Acad. Sci. USA 2007, 104, 1745–1750. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Suckling, C.C.; Clark, M.S.; Richard, J.; Morley, S.A.; Thorne, M.A.S.; Harper, E.M.; Peck, L.S. Adult Acclimation to Combined Temperature and pH Stressors Significantly Enhances Reproductive Outcomes Compared to Short-term Exposures. J. Anim. Ecol. 2015, 84, 773–784. [Google Scholar] [CrossRef]
  47. Zhao, C.; Zhang, L.; Shi, D.; Ding, J.; Yin, D.; Sun, J.; Zhang, B.; Zhang, L.; Chang, Y. Transgenerational Effects of Ocean Warming on the Sea Urchin Strongylocentrotus Intermedius. Ecotoxicol. Environ. Saf. 2018, 151, 212–219. [Google Scholar] [CrossRef]
  48. Schlegel, P.; Havenhand, J.N.; Gillings, M.R.; Williamson, J.E. Individual Variability in Reproductive Success Determines Winners and Losers under Ocean Acidification: A Case Study with Sea Urchins. PLoS ONE 2012, 7, e53118. [Google Scholar] [CrossRef] [Green Version]
  49. Sato, K.N.; Andersson, A.J.; Day, J.M.D.; Taylor, J.R.A.; Frank, M.B.; Jung, J.-Y.; McKittrick, J.; Levin, L.A. Response of Sea Urchin Fitness Traits to Environmental Gradients across the Southern California Oxygen Minimum Zone. Front. Mar. Sci. 2018, 5, 258. [Google Scholar] [CrossRef]
  50. Vafidis, D.; Antoniadou, C.; Kyriakouli, K. Reproductive Cycle of the Edible Sea Urchin Paracentrotus Lividus (Echinodermata: Echinoidae) in the Aegean Sea. Water 2019, 11, 1029. [Google Scholar] [CrossRef] [Green Version]
  51. Zhang, W.; Chang, Y.; Luo, S.; Zhou, H.; Tian, X.; Ding, J.; Chen, X. Effects of Biofilms as the Main and as a Supplementary Food on the Survival, Somatic Growth and Gonad Enhancement of Sea Urchin Strongylocentrotus Intermedius. Aquac. Int. 2014, 22, 925–936. [Google Scholar] [CrossRef]
  52. Rinde, E.; Christie, H.; Fagerli, C.W.; Bekkby, T.; Gundersen, H.; Norderhaug, K.M.; Hjermann, D.Ø. The Influence of Physical Factors on Kelp and Sea Urchin Distribution in Previously and Still Grazed Areas in the NE Atlantic. PLoS ONE 2014, 9, e100222. [Google Scholar] [CrossRef] [Green Version]
  53. López, S.; Turon, X.; Montero, E.; Palacín, C.; Duarte, C.M.; Tarjuelo, I. Larval Abundance, Recruitment and Early Mortality in Paracentrotus Lividus (Echinoidea). Interannual Variability and Plankton-Benthos Coupling. Mar. Ecol. Prog. Ser. 1998, 172, 239–251. [Google Scholar] [CrossRef] [Green Version]
  54. Wang, H.; Ding, J.; Ding, S.; Chang, Y. Integrated Metabolomic and Transcriptomic Analyses Identify Critical Genes in Eicosapentaenoic Acid Biosynthesis and Metabolism in the Sea Urchin Strongylocentrotus Intermedius. Sci. Rep. 2020, 10, 1697. [Google Scholar] [CrossRef] [PubMed]
  55. Picard, V.; Mulner-Lorillon, O.; Bourdon, J.; Morales, J.; Cormier, P.; Siegel, A.; Bellé, R. Model of the Delayed Translation of Cyclin B Maternal MRNA after Sea Urchin Fertilization. Mol. Reprod. Dev. 2016, 83, 1070–1082. [Google Scholar] [CrossRef] [Green Version]
  56. Shi, D.; Zhao, C.; Chen, Y.; Ding, J.; Zhang, L.; Chang, Y. Transcriptomes Shed Light on Transgenerational and Developmental Effects of Ocean Warming on Embryos of the Sea Urchin Strongylocentrotus Intermedius. Sci. Rep. 2020, 10, 7931. [Google Scholar] [CrossRef] [PubMed]
  57. Klosin, A.; Lehner, B. Mechanisms, Timescales and Principles of Trans-Generational Epigenetic Inheritance in Animals. Curr. Opin. Genet. Dev. 2016, 36, 41–49. [Google Scholar] [CrossRef] [PubMed]
  58. Bonduriansky, R.; Crean, A.J.; Day, T. The Implications of Nongenetic Inheritance for Evolution in Changing Environments. Evol. Appl. 2012, 5, 192–201. [Google Scholar] [CrossRef]
  59. Bertoldo, M.J.; Locatelli, Y.; O’Neill, C.; Mermillod, P.; Bertoldo, M.J.; Locatelli, Y.; O’Neill, C.; Mermillod, P. Impacts of and Interactions between Environmental Stress and Epigenetic Programming during Early Embryo Development. Reprod. Fertil. Dev. 2015, 27, 1125–1136. [Google Scholar] [CrossRef]
  60. Masullo, T.; Biondo, G.; Natale, M.D.; Tagliavia, M.; Bennici, C.D.; Musco, M.; Ragusa, M.A.; Costa, S.; Cuttitta, A.; Nicosia, A. Gene Expression Changes after Parental Exposure to Metals in the Sea Urchin Affect Timing of Genetic Programme of Embryo Development. Biology 2021, 10, 103. [Google Scholar] [CrossRef]
  61. Goldstone, J.V.; Hamdoun, A.; Cole, B.J.; Howard-Ashby, M.; Nebert, D.W.; Scally, M.; Dean, M.; Epel, D.; Hahn, M.E.; Stegeman, J.J. The Chemical Defensome: Environmental Sensing and Response Genes in the Strongylocentrotus Purpuratus Genome. Dev. Biol. 2006, 300, 366–384. [Google Scholar]
  62. Harris, L.G.; Eddy, S.D. Sea Urchin Ecology and Biology. In Echinoderm Aquaculture; Wiley: Hoboken, NJ, USA, 2015; Part I, Chapter 1; pp. 1–24. [Google Scholar]
  63. McBride, S.C. Sea Urchin Aquaculture. In American Fisheries Society Symposium; American Fisheries Society: Bethesda, MD, USA, 2005; Volume 46, p. 179. [Google Scholar]
  64. Ragusa, M.A.; Costa, S.; Gianguzza, M.; Roccheri, M.C.; Gianguzza, F. Effects of Cadmium Exposure on Sea Urchin Development Assessed by SSH and RT-QPCR: Metallothionein Genes and Their Differential Induction. Mol. Biol. Rep. 2013, 40, 2157–2167. [Google Scholar]
  65. Pfaffl, M.W. A New Mathematical Model for Relative Quantification in Real-Time RT–PCR. Nucleic Acids Res. 2001, 29, e45. [Google Scholar] [CrossRef] [PubMed]
  66. Pfaffl, M.W.; Horgan, G.W.; Dempfle, L. Relative Expression Software Tool (REST©) for Group-Wise Comparison and Statistical Analysis of Relative Expression Results in Real-Time PCR. Nucleic Acids Res. 2002, 30, e36. [Google Scholar] [CrossRef] [PubMed]
  67. Szklarczyk, D.; Gable, A.L.; Lyon, D.; Junge, A.; Wyder, S.; Huerta-Cepas, J.; Simonovic, M.; Doncheva, N.T.; Morris, J.H.; Bork, P. STRING V11: Protein–Protein Association Networks with Increased Coverage, Supporting Functional Discovery in Genome-Wide Experimental Datasets. Nucleic Acids Res. 2019, 47, D607–D613. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Average survival rates of larvae obtained from the two pools of individuals (A+B vs. C+D) in smaller (A) and larger (B) tanks.
Figure 1. Average survival rates of larvae obtained from the two pools of individuals (A+B vs. C+D) in smaller (A) and larger (B) tanks.
Ijms 23 12790 g001
Figure 2. Average larval survival rates evaluated during the first week in (A) smaller and (B) largertanks, referring to the recruits of females A+B vs. those of females C+D. The linear regressions are superimposed.
Figure 2. Average larval survival rates evaluated during the first week in (A) smaller and (B) largertanks, referring to the recruits of females A+B vs. those of females C+D. The linear regressions are superimposed.
Ijms 23 12790 g002
Figure 3. (A) Interactomic analysis by STRING (https://string–db.org/; accessed on 30 April 2022). The network graphically displays the relationship between genes. The biological relationships between genes are indicated by different colors. Known interactions: reported by database = light blue and determined experimentally = pink. Expected interactions: gene proximity = green; gene fusion = red; and genes with similar pattern = light blue. (B) Homo sapiens gene names and the corresponding P. lividus orthologous genes. The most significant relations among genes (confidence score cut-off = 900) displaying experimental evidence are highlighted.
Figure 3. (A) Interactomic analysis by STRING (https://string–db.org/; accessed on 30 April 2022). The network graphically displays the relationship between genes. The biological relationships between genes are indicated by different colors. Known interactions: reported by database = light blue and determined experimentally = pink. Expected interactions: gene proximity = green; gene fusion = red; and genes with similar pattern = light blue. (B) Homo sapiens gene names and the corresponding P. lividus orthologous genes. The most significant relations among genes (confidence score cut-off = 900) displaying experimental evidence are highlighted.
Ijms 23 12790 g003
Figure 4. Heatmaps (Heatmapper available at www.heatmapper.ca, accessed on 28 July 2022) showing the expression profiles and hierarchical clustering of twenty-three genes analyzed by Real Time qPCR in P. lividus embryos deriving from females A, B, C and D. Color code: red = up-regulated genes with respect to the control; green = down-regulated genes.
Figure 4. Heatmaps (Heatmapper available at www.heatmapper.ca, accessed on 28 July 2022) showing the expression profiles and hierarchical clustering of twenty-three genes analyzed by Real Time qPCR in P. lividus embryos deriving from females A, B, C and D. Color code: red = up-regulated genes with respect to the control; green = down-regulated genes.
Ijms 23 12790 g004
Table 1. Paired t-test (two-tailed p value) to compare the larval survivorships (A-B vs. C-D) in smaller and larger tanks. ** indicates p < 0.01.
Table 1. Paired t-test (two-tailed p value) to compare the larval survivorships (A-B vs. C-D) in smaller and larger tanks. ** indicates p < 0.01.
Paired t TestSmallerLarger
pvalue0.08300.1169
pvalue summaryn.s.n.s.
t1.9801.758
df88
Mean of differences (B-A)−21.93−8.778
SD of differences33.2214.98
SEM of differences11.074.994
95% confidence interval−47.46 to 3.610−20.29 to 2.738
R squared (partial eta squared)0.32890.2786
Correlation coefficient (r)0.76080.8595
p value (one tailed)0.00860.0015
p value summary****
Significantly effective pairingYesYes
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Federico, S.; Glaviano, F.; Esposito, R.; Pinto, B.; Gharbi, M.; Di Cosmo, A.; Costantini, M.; Zupo, V. Gene Expression Detects the Factors Influencing the Reproductive Success and the Survival Rates of Paracentrotus lividus Offspring. Int. J. Mol. Sci. 2022, 23, 12790. https://doi.org/10.3390/ijms232112790

AMA Style

Federico S, Glaviano F, Esposito R, Pinto B, Gharbi M, Di Cosmo A, Costantini M, Zupo V. Gene Expression Detects the Factors Influencing the Reproductive Success and the Survival Rates of Paracentrotus lividus Offspring. International Journal of Molecular Sciences. 2022; 23(21):12790. https://doi.org/10.3390/ijms232112790

Chicago/Turabian Style

Federico, Serena, Francesca Glaviano, Roberta Esposito, Bruno Pinto, Maissa Gharbi, Anna Di Cosmo, Maria Costantini, and Valerio Zupo. 2022. "Gene Expression Detects the Factors Influencing the Reproductive Success and the Survival Rates of Paracentrotus lividus Offspring" International Journal of Molecular Sciences 23, no. 21: 12790. https://doi.org/10.3390/ijms232112790

APA Style

Federico, S., Glaviano, F., Esposito, R., Pinto, B., Gharbi, M., Di Cosmo, A., Costantini, M., & Zupo, V. (2022). Gene Expression Detects the Factors Influencing the Reproductive Success and the Survival Rates of Paracentrotus lividus Offspring. International Journal of Molecular Sciences, 23(21), 12790. https://doi.org/10.3390/ijms232112790

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop