Selection Shapes the Genomic Landscape of Introgressed Ancestry in a Pair of Sympatric Sea Urchin Species.

A growing number of recent studies have demonstrated that introgression is common across the tree of life. However, we still have a limited understanding of the fate and fitness consequence of introgressed variation at the whole-genome scale across diverse taxonomic groups. Here, we implemented a phylogenetic hidden Markov model to identify and characterize introgressed genomic regions in a pair of well-diverged, nonsister sea urchin species: Strongylocentrotus pallidus and Strongylocentrotus droebachiensis. Despite the old age of introgression, a sizable fraction of the genome (1% to 5%) exhibited introgressed ancestry, including numerous genes showing signals of historical positive selection that may represent cases of adaptive introgression. One striking result was the overrepresentation of hyalin genes in the identified introgressed regions despite observing considerable overall evidence of selection against introgression. There was a negative correlation between introgression and chromosome gene density, and two chromosomes were observed with considerably reduced introgression. Relative to the nonintrogressed genome-wide background, introgressed regions had significantly reduced nucleotide divergence (dXY) and overlapped fewer protein-coding genes, coding bases, and genes with a history of positive selection. Additionally, genes residing within introgressed regions showed slower rates of evolution (dN, dS, dN/dS) than random samples of genes without introgressed ancestry. Overall, our findings are consistent with widespread selection against introgressed ancestry across the genome and suggest that slowly evolving, low-divergence genomic regions are more likely to move between species and avoid negative selection following hybridization and introgression.

Establishment of Parentage Identification Method for Sea Urchin Strongylocentrotus intermedius Based on SSR-seq Technology.

To establish a parentage identification method for Strongylocentrotus intermedius, 15 microsatellite loci and simple sequence repeat sequencing (SSR-seq) technology were used to perform SSR sequencing and typing of the validation population with known pedigree information and the simulation population. Cervus v3.0 was used for gene frequency statistics, simulated analysis, and parentage identification analysis. The results showed that, in validation population, using 15 microsatellite loci, the highest success rate of parent pairs identification was 86%, the highest success rate of female parent identification was 93%, and the highest success rate of male parent identification was 90%. The simulated population was analyzed using 12-15 loci, and the identification rate was up to 90%. In cases where accurate parentage was not achieved, individuals could exhibit genetic similarities with 1-3 male or female parents. Individuals identified as lacking a genetic relationship can be selected as parents to prevent inbreeding. This study shows that parent pairs or single parents of most offspring can be identified successfully using these 15 selected loci. The results lay a foundation for the establishment of a parentage identification method for S. intermedius.

Trans-Arctic vicariance in Strongylocentrotus sea urchins.

The sea urchins Strongylocentotus pallidus and S. droebachiensis first invaded the Atlantic Ocean from the Pacific following the opening of the Bering seaway in the late Miocene. While trans-Arctic dispersal during the Pleistocene is thought to have maintained species' integrity, a recent genomic analysis identified a reproductively isolated cryptic species within S. droebachiensis. Based on previous studies, the distribution of one of these lineages (S. droebachiensis W) includes the shallow water habitats of the northwest Atlantic and Pacific, while the other (S. droebachiensis E) is found throughout the shallow habitat in the northeast but is mostly restricted to deep habitats (>65 m) in the northwest Atlantic. However, since genetic variation within S. droebachiensis has been largely unstudied in the north Pacific and Arctic oceans, the biogeography of the cryptic species is not well known, and it is difficult to identify the mechanisms driving population subdivision and speciation. Here we use population genetic analyses to characterize the distribution of each species, and to test hypotheses about the role of vicariance in the evolution of systematic and genomic divergence within the genus. We collected individuals of all three Strongylocentrotus species (n = 365) from 10 previously unsampled locations in the northeast Pacific and north Atlantic (Labrador Sea and Norway), and generated mtDNA sequence data for a 418 bp fragment of cytochrome c oxidase subunit I (COI). To assess the biogeography of all three species, we combined our alignment with five previously published data sets (total n = 789) and used statistical parsimony and maximum likelihood to identify species and characterize their distribution within and among oceans. Patterns of haplotype sharing, pairwise F ST , and hierarchical analyses of molecular variance (AMOVA) identified trans-Arctic dispersal in S. pallidus and S. droebachiensis W, but other than 5 previously reported singletons we failed to detect additional mtDNA haplotypes of S. droebachiensis E in the north Pacific. Within the Atlantic, patterns of habitat segregation suggests that temperature may play a role in limiting the distribution of S. droebachiensis E, particularly throughout the warmer coastal waters along the coast of Nova Scotia. Our results are consistent with the cycles of trans-Arctic dispersal and vicariance in S. pallidus and S. droebachiensis W, but we suggest that the evolution of Atlantic populations of S. droebachiensis E has been driven by persistent trans-Arctic vicariance that may date to the initial invasion in the late Pliocene.

Isolation and characterization of four pathogenic strains of Vibrio spp. from farmed sea urchins, Strongylocentrotus intermedius, in China.

The sea urchin Strongylocentrotus intermedius, famous for its gonadal quality, is one of the most important farmed species in the sea area of northern China. Since 2020, outbreaks of black peristomial membrane disease (commonly called black mouth disease) have frequently occurred in spring and winter in cultured S. intermedius. In this study, we isolated the predominant bacteria from different tissues of diseased sea urchins from a North China farm in the spring of 2021. Four pathogenic strains (named SIBMPM01, SIBMPM02, SIBMPM03 and SIBMCF01) were obtained and characterized by Gram staining, morphological observation, artificial infection tests, and metabolic characteristics. Our results showed that: 1) all obtained strains belonged to the genus Vibrio and had morphological differences. Phylogenetic analysis indicated that the four obtained strains might be novel Vibrio species. 2) Laboratory-based artificial infection tests showed that sea urchins infected with either SIBMPM01, SIBMPM02, SIBMPM03 or SIBMCF01 exhibited pathological symptoms of a black peristomial membrane in a dosage-dependent and temperature-dependent manner. The virulence of SIBMCF01 was greater than those of the others. 3) Metabolic characterization data showed that SIBMPM01, SIBMPM02, SIBMPM03 and SIBMCF01 shared similar metabolic characteristics. 4) Antimicrobial susceptibility tests demonstrated that the four obtained strains were all sensitive to ampicillin, doxycycline, norfloxacin, ofloxacin, furazolidone and chloramphenicol. SIBMPM01 was specifically sensitive to neomycin, and SIBMCF01 was specifically sensitive to carboxybenzyl penicillin.

Evidence for association of Vibrio echinoideorum with tissue necrosis on test of the green sea urchin Strongylocentrotus droebachiensis.

"Sea urchin lesion syndrome" is known as sea urchin disease with the progressive development of necrotic epidermal tissue and loss of external organs, including appendages on the outer body surface. Recently, a novel strain, Vibrio echinoideorum has been isolated from the lesion of green sea urchin (Strongylocentrotus droebachiensis), an economically important mariculture species in Norway. V. echinoideorum has not been reported elsewhere in association with green sea urchin lesion syndrome. Therefore, in this study, an immersion based bacterial challenge experiment was performed to expose sea urchins (wounded and non-wounded) to V. echinoideorum, thereby mimicking a nearly natural host-pathogen interaction under controlled conditions. This infection experiment demonstrated that only the injured sea urchins developed the lesion to a significant degree when exposed to V. echinoideorum. Pure cultures of the employed bacterial strain were recovered from the infected animals and its identity was confirmed by the MALDI-TOF MS spectra profiling. Additionally, the hemolytic phenotype of V. echinoideorum substantiated its virulence potential towards the host, and this was also supported by the cytolytic effect on red spherule cells of sea urchin. Furthermore, the genome sequence of V. echinoideorum was assumed to encode potential virulence genes and were subjected to in silico comparison with the established virulence factors of Vibrio vulnificus and Vibrio tasmaniensis. This comparative virulence profile provided novel insights about virulence genes and their putative functions related to chemotaxis, adherence, invasion, evasion of the host immune system, and damage of host tissue and cells. Thus, it supports the pathogenicity of V. echinoideorum. In conclusion, the interaction of V. echinoideorum with injured sea urchin facilitates the development of lesion syndrome and therefore, revealing its potentiality as an opportunistic pathogen.

Application of SNP in Genetic Sex Identification and Effect of Estradiol on Gene Expression of Sex-Related Genes in Strongylocentrotus intermedius.

Sea urchin (Strongylocentrotus intermedius) is an economically important mariculture species in Asia, and its gonads are the only edible part. The efficiency of genetic breeding in sea urchins is hampered due to the inability to distinguish gender by appearance. In this study, we first identified a sex-associated single nucleotide polymorphism (SNP) by combining type IIB endonuclease restriction site-associated DNA sequencing (2b-RAD-seq) and genome survey. Importantly, this SNP is located within spata4, a gene specifically expressed in male. Knocking down of spata4 by RNA interference (RNAi) in male individuals led to the downregulation of other conserved testis differentiation-related genes and germ cell marker genes. We also revealed that sex ratio in this validated culture population of S. intermedius is not 1:1. Moreover, after a 58-day feeding experiment with estradiol, the expression levels of several conserved genes that are related to testis differentiation, ovary differentiation, and estrogen metabolism were dynamically changed. Taken together, our results will contribute toward improving breeding efficiency, developing sex-controlled breeding, and providing a solid base for understanding sex determination mechanisms in sea urchins.

ATAC-Seq for Assaying Chromatin Accessibility Protocol Using Echinoderm Embryos.

Cis-regulatory elements (CREs) and transcription factors (TFs) associated with them determine temporal and spatial domains of gene expression. Therefore, identification of these CREs and TFs is crucial to elucidating transcriptional programs across taxa. With chromatin accessibility facilitating transcription factor access to DNA, the identification of regions of open chromatin sheds light both on the function of the regulatory elements and their evolution, thus allowing the recognition of potential CREs. Buenrostro and colleagues have developed a novel method for exploring chromatin accessibility: assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq), which can be used for the purpose of identifying putative CREs. This method was shown to have considerable advantages when compared to traditional methods such as sequence conservation analyses or functional assays. Here we present the adaptation of the ATAC-seq method to echinoderm species and discuss how it can be used for CRE discovery.

Distant hybrids of Heliocidaris crassispina (♀) and Strongylocentrotus intermedius (♂): identification and mtDNA heteroplasmy analysis.

BACKGROUND: Distant hybridization between the sea urchin Heliocidaris crassispina (♀) and the sea urchin Strongylocentrotus intermedius (♂) was successfully performed under laboratory conditions. A new variety of hybrid sea urchin (HS hybrid) was obtained. However, the early-development success rates for the HS hybrids were significantly lower than those of purebred H. crassispina or S. intermedius offspring. In addition, it was difficult to distinguish the HS-hybrid adults from the pure H. crassispina adults, which might lead to confusion in subsequent breeding attempts. In this study, we attempted to develop a method to quickly and effectively identify HS hybrids, and to preliminarily investigate the molecular mechanisms underlying the poor early-development success rates in the HS hybrids. RESULTS: The hybrid sea urchins (HS hybrids) were identified both morphologically and molecularly. There were no significant differences in the test height to test diameter ratios between the HS hybrids and the parents. The number and arrangement of ambulacral pore pairs in the HS hybrids differed from those of the parental lines, which might serve as a useful morphological character for the identification of the HS hybrids. A primer pair that identified the HS hybrids was screened by comparing the mitochondrial genomes of the parental lines. Moreover, paternal leakage induced mitochondrial DNA heteroplasmy in the HS hybrids, which might explain the low rates of early development success in these hybrids. CONCLUSIONS: The distant-hybrid sea urchins were accurately identified using comparative morphological and molecular genetic methods. The first evidence of mtDNA heteroplasmy after the distant hybridization of an echinoderm was also provided.

Transcriptomes shed light on transgenerational and developmental effects of ocean warming on embryos of the sea urchin Strongylocentrotus intermedius.

Ocean warming increasingly endangers the fitness of marine invertebrates. Transgenerational effects (TE) potentially mitigate the impacts of environmental stress on the embryos of marine invertebrates. The molecular mechanisms, however, remain largely unknown. Using high-throughput RNA sequencing technology, we investigated the gene expression patterns of embryos (the gastrula stage) of the sea urchin Strongylocentrotus intermedius at different developmental temperatures, whose parents were exposed to long-term (15 months) elevated temperature (A) or not (B). The temperatures at which adults were held for ~4 weeks prior to the start of the experiment (21 °C for A and 18 °C for B) were also used for the development of offspring (high: 21 °C and ambient (laboratory): 18 °C) resulting in four experimental groups (HA and HB at 21 °C, and LA and LB at 18 °C). The embryos were sampled ~24 h after fertilization. All samples were in the gastrula stage. Twelve mRNA libraries (groups HA, HB, LA, LB, 3 replicates for each group) were established for the following sequencing. Embryos whose parents were exposed to elevated temperatures or not showed 1891 significantly different DEGs (differentially expressed genes) at the ambient developmental temperature (LB vs LA, LB as control) and 2203 significantly different DEGs at the high developmental temperature (HB vs HA, HB as control), respectively. This result indicates complex molecular mechanisms of transgenerational effects of ocean warming, in which a large number of genes are involved. With the TE, we found 904 shared DEGs in both LB vs LA (LB as control) and HB vs HA (HB as control) changed in the same direction of expression (i.e., up- or down-regulated), indicating that parental exposed temperatures affect the expression of these genes in the same manner regardless of the development temperature. With developmental exposure, we found 198 shared DEGs in both HB vs LB (HB as control) and HA vs LA (HA as control) changed in the same direction of expression. Of the 198 DEGs, more genes were up-regulated at high developmental temperature. Interestingly, embryos whose parents were exposed to high temperature showed fewer differently expressed DEGs between high and low developmental temperatures than the individuals whose parents were exposed to ambient temperature. The results indicate that gene expressions are probably depressed by the transgenerational effect of ocean warming. The roles of hsp70 and hnf6 in thermal acclimation are highlighted for future studies. The present study provides new insights into the molecular mechanisms of the transgenerational and developmental effects of ocean warming on the embryos of sea urchins.

Transcriptome sequencing reveals phagocytosis as the main immune response in the pathogen-challenged sea urchin Strongylocentrotus intermedius.

The clarification of host immune responses to causative bacteria of spotting disease in the sea urchin Strongylocentrotus intermedius is vital to preventing and controlling this disease, especially to selective breeding for disease resistance. For this purpose, sea urchins were challenged with the causative bacterium Vibrio sp. to obtain spotting diseased and undiseased samples. We conducted next-generation sequencing to assess the key genes/pathways in control (CG), diseased (DG), and undiseased (UG) groups. A total of 454.1 million clean reads were obtained and assembled into 23,899 UniGenes with an N50 of 1359 bp, with 86.11% of them matching the genome sequence of the sea urchin S. purpuratus. A total of 8415 UniGenes were mapped to the non-redundant database. Salmon expression analysis revealed 725 significantly differentially expressed genes (DEGs) among CG, DG, and UG. These DEGs were enriched into 72 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, including a core set of immune correlated pathways notably in the phagosome, vitamin digestion and absorption, Wnt signaling, and Notch signaling pathways. DG was evidently upregulated in these immune pathways and could enhance phagocytosis directly or indirectly. Thus, phagocytosis was the main coelomic cellular immune response in S. intermedius challenged by spotting disease causative bacterium. The expression patterns of 10 DEGs were confirmed via RT-qPCR, and the expression levels were consistent with the results of RNA-seq. Furthermore, 9899 SSRs were identified, and 123,692, 151,827, and 114,368 candidate SNPs were identified from CG, DG, and UG, respectively. These results provide basic information for our understanding of sea urchin antibacterial immunity.

Transcriptome analysis to characterize the genes related to gonad growth and fatty acid metabolism in the sea urchin Strongylocentrotus intermedius.

BACKGROUND: Sea urchin gonads of both sexes, commonly termed "roe", are highly valued seafood delicacies, and Strongylocentrotus intermedius is considered one of the tastiest sea urchins. In order to produce high-quality gonads for consumption and clarify the mechanism of gonad growth and development of the sea urchin, more genetic information, especially at the transcriptome level, is needed. OBJECTIVE: A more thorough understanding of sea urchin gonad growth and development in both sexes could enable regulation of these processes at several stages with the aim of suppressing gametogenesis in order to produce high-quality gonads for consumption. METHODS: The adult sea urchins S. intermedius were cultured for 3 months, and were sampled for the gonadal transcriptome analysis which has been performed on the RNAs of three male and female adults of S. intermedius in each gonad development stage. RESULTS: Illumina sequencing raw sequence data was deposited in the NCBI Sequence Read Archive (SRA) database (PRJNA532998). It generated 560,196,356 raw reads and 548,956,944 clean reads were acquired, which were assembled into 107,850 transcripts with 44,124 genes. Comparative analysis showed the differentially expressed genes (DEGs) from 114 to 2566. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were used to determine the functional significance of these DEGs. We have selected 9 genes related to growth and 12 genes related to fatty acid biosynthesis and metabolism in sea urchin gonads. CONCLUSION: These data for sea urchins were intended to provide markers for gonad growth and development that can be accumulated for use in aquaculture applications.

Trans-Species Polymorphism in Mitochondrial Genome of Camarodont Sea Urchins.

Mitochondrial (mt) genomes of the sea urchins Strongylocentrotus intermedius and Mesocentrotus nudus demonstrate the identical patterns of intraspecific length variability of the ND6 gene, consisting of 489 bp (S variant) and 498 bp (L variant), respectively. For both species, the ND6 length difference is due to the 488A>G substitution, which changes the stop codon TAG in S variant for a tryptophan codon TGG in L variant and elongates the corresponding ND6 protein by three additional amino acids, Trp-Leu-Trp. The phylogenetic analysis based on mt genomes of sea urchins and related echinoderm groups from GenBank has shown the S and L ND6 variants as shared among the camarodont sea urchins; the rest of the echinoderms demonstrate the S variant only. The data suggest that the ND6 488A>G substitution can be the first example of the trans-species polymorphism in sea urchins, persisting at least since the time of the Odontophora diversification at the Eocene/Oligocene boundary (approximately 34 million years ago), which was characterized by an abrupt climate change and significant global ocean cooling. Alternative hypotheses, including the convergent RNA editing and/or codon reassignment, are not supported by direct comparisons of the ND6 gene sequences with the corresponding transcripts using the basic local alignment search tool (BLAST) of full sea urchin transcriptomes.

Molecular characterization and expression of SiFad1 in the sea urchin (Strongylocentrotus intermedius).

Fatty acid desaturases (Fads) are a key enzyme in the process of biosynthesis of highly unsaturated fatty acids (HUFAs). In this study, we cloned the full-length sequence of the SiFad1 gene (SiFad1) and analyzed its expression profiles during different developmental stages and in different tissues of Strongylocentrotus intermedius. The full-length cDNA of SiFad1 is composed of 1086 bp, with a putative open reading frame of 885 bp encoding a polypeptide of 294 amino acid (AA) residues. The predicted molecular mass of SiFad1 is 34.67 kDa and its theoretical pI is 8.41. The presence of conserved motifs including three histidine boxes (HXXXH, HXXHH, XXXHH), a FA_desaturases domain and three transmembrane domains suggests that SiFad1 belongs to the microsomal fatty acid desaturases family. Its tissue distribution showed that the highest expression of SiFad1 is in the intestine and the weakest expression is in Aristotle's lantern of S. intermedius. Time-course expression measurements in different developmental stages showed the highest expression of SiFad1 occurs in the gastrula and the weakest expression in the juvenile sea urchin. Knock-down of SiFad1 by specific siRNA revealed that the significantly depressed expression of Elovl5 had decreased in the coelomocytes, intestines and gonads at 24 h post transfection, indicating that the downstream target gene of SiFad1 is Elovl5 and SiFad1 and Elovl5 have positive regulatory effects. When we examined the changes in fatty acids in the gonads before and after interference, the results showed that after 24 h of interference, the content of C20:4n-6 produced by SiFad1 had decreased. Taken together, these results will enable us to understand the role of SiFad1 in fatty acid anabolism, which will help us to understand the fatty acid synthesis pathways and regulatory mechanisms of Strongylocentrotus intermedius and provide a theoretical experimental basis for improving the ability of sea urchins to synthesize fatty acids and cultivating sea urchins of higher quality and nutritional value.

De novo transcriptomic analysis of gonad of Strongylocentrotus nudus and gene discovery for biosynthesis of polyunsaturated fatty acids.

BACKGROUND: Strongylocentrotus nudus is an important cultured sea urchin species in north China, because its gonad is rich in unsaturated fatty acids, particularly long polyunsaturated fatty acids (LC-PUFAs). These PUFAs play pleiotropic and crucial roles in a wide range of biological process. OBJECTIVE: However, the genes contributing to biosynthesis PUFAs have not been elucidated yet, and the molecular mechanism relative to the difference in PUFA composition between male and female gonad as been revealed but the corresponding has not been understood. METHODS: In this paper, solexa sequencing based transcriptomic approach was used to identify and characterize the key genes relative to PUFA synthesis and further conducted different expressed genes between male and female gonad. RESULTS: A total of 130,124 transcripts and 189330 unigenes were de novo assembled from 64.32 Gb data. Next, these unigenes were subjected to functional annotation by mapping to six public databases, and this process revealed a lot of genes involving in lipid metabolism. In addition, three types of fatty acids front-end desaturase and three species of very long fatty acids elongase were identified and the pathway for PUFA biosynthesis was hypothesized. Last, comparative analysis revealed the higher expression level of Δ5 desaturase, Δ6 desaturase, ELOVL-4, -6 and -7 in male gonad compared with female. CONCLUSION: This results could plausible explain the differ in composition of PUFAs between male and female gonad of sea urchin.

Transcriptome profiling reveals key roles of phagosome and NOD-like receptor pathway in spotting diseased Strongylocentrotus intermedius.

Spotting disease is a common disease in the process of aquaculture and restocking of the sea urchin Strongylocentrotus intermedius and leads to mass mortality. To characterize the molecular processes and candidate genes related to spotting disease in S. intermedius, we conducted next-generation sequencing to assess the key genes/pathways in spotting diseased sea urchin (DUG) compared to healthy ones (HUG). A total of 321.1 million clean reads were obtained and assembled into 93,877 Unigenes with an N50 of 1185 bp, in which 86.48% of them matched to the genome sequence of the sea urchin S. purpuratus and 27,456 Unigenes mapped to Nr database. Salmon expression analysis revealed 1557 significantly differently expressed genes (DEGs) between DUG and HUG. These DEGs were enriched into 151 KEGG pathways including a core set of immune correlated pathways notably in phagosome and NOD-like receptor signaling. DUG displayed an obvious downregulation in these immune pathways. The expression patterns of six DEGs were confirmed by RT-qPCR, and the expressions were consistent with the results of RNA-seq. Furthermore, 15,990 SSRs were identified and a total of 235,249 and 295,567 candidate SNPs were identified from DUG and HUG, respectively. All these results provided basic information for our understanding of spotting disease outbreak in sea urchin.

Transcriptomes reveal genes involved in covering and sheltering behaviors of the sea urchin Strongylocentrotus intermedius exposed to UV-B radiation.

Although the potential link exists between behavioral responses to UV-B radiation and the maximization of fitness, molecular mechanisms of these UV-B induced behaviors remain poorly understood. For the first time, we investigated the transcriptomes of covered (CB), sheltered (SB) and non-protected (NA) sea urchins Strongylocentrotus intermedius exposed to UV-B radiation. A total of 330 differentially expressed genes were revealed by transcriptome comparisons. By comparing with the group NA, we found 79 up-regulated and 118 down-regulated genes in SB group, as well as 26 up-regulated and 67 down-regulated genes in group CB. There were 34 up-regulated genes and 52 down-regulated genes in group SB, compared with group CB. These differentially expressed genes failed to enrich either Gene Ontology (GO) or Kyoto Encyclopedia of Genes and Genomes (KEGG), only except an enrichment in KEGG. We highlighted TRPA1 and Opsin as key neurobiological genes involved in the molecular mechanisms of covering and sheltering behaviors of sea urchins exposed to UV-B radiation. What's more, other identified genes provide valuable resources for future investigations on the molecular basis of covering and sheltering behaviors of sea urchins.

In silico assessment and structural characterization of antioxidant peptides from major yolk protein of sea urchin Strongylocentrotus nudus.

Sea urchin gonads have been demonstrated to contain major yolk protein (MYP), which can be hydrolyzed by enzymes to release biologically active peptides. The in silico analysis of the MYP sequence in the BIOPEP database showed the presence of fragments with antioxidant activity. The sequence was hydrolyzed by 21 kinds of proteases and 23 antioxidant peptides were obtained. Eight peptides, including Leu-Trp (LW), Arg-Trp (RW), Ala-Trp (AW), Thr-Trp (TW), Ala-Asp-Phe (ADF), Leu-Trp-Lys (LWK), Ser-Asp-Phe (SDF) and Leu-Tyr (LY), were screened and a score over 0.5 was obtained using PeptideRanker. The peptides LW, TW and LWK showed a stronger antioxidant capacity with IC50 values of 8.85, 9.59 and 9.62 mmol L-1, respectively, compared to that of glutathione (10.81 mmol L-1). Furthermore, AW, LW and LY showed Trolox equivalent antioxidant capacity (TEAC) values of 3.07, 1.87 and 1.52 mmol TE per mmol peptide, respectively. These results suggest that the MYP from sea urchin (S. nudus) gonads is a good source of antioxidant peptides with abundant tryptophan.

Transgenerational effects of heavy fuel oil on the sea urchin Strongylocentrotus intermedius considering oxidative stress biomarkers.

Stranding of oil onto a coastline after an oil spill threatens the health of marine benthic organisms. Here, the transgenerational effects of exposure to stranded heavy fuel oil (HFO) on the sea urchin Strongylocentrotus intermedius were assessed. The column containing gravel coated with HFO was prepared in the laboratory to simulate HFO-contaminated gravel shorelines. Adult sea urchins were exposed for 21 days to either a HFO-oiled gravel column at the oil loading of 3000 µg oil/g gravel or a non-HFO-oiled gravel column (as the control treatment) and then offspring were either exposed to HFO or ambient seawater conditions. The sublethal exposure to HFO for 21 days induced polycyclic aromatic hydrocarbons (PAHs) accumulation in gonads, accompanied by increased levels of oxidative lipid, protein and DNA damage and a reduction in total antioxidant capacity. Analysis of gametes indicated that both maternal and paternal exposure could result in the transfer of PAHs and DNA damage to their offspring. Parental (maternal, paternal or both) exposure to HFO caused increases in malformation rates of offspring compared to those from control parents under ambient seawater condition. Continued HFO exposure in the offspring resulted in further increased malformation rates compared with those reared in ambient seawater, as well as oxidative lipid, protein and DNA damage. Furthermore, mother exposure history reduced the total antioxidant capacity of larvae to response to continued HFO exposure. Overall, the results suggest an increased sensitivity to toxic effects of HFO in larvae from exposed both parents compared with those from control parents in S. intermedius, which may consequently affect the recruitment and population maintenance.

SLAF-based high-density genetic map construction and QTL mapping for major economic traits in sea urchin Strongylocentrotus intermedius.

Sea urchin (Strongylocentrotus intermedius) has long been a model species for developmental and evolutionary research, but only a few studies have focused on gene mapping. Here, we reported a high-density genetic map containing 4,387 polymorphism specific-length amplified fragment (SLAF) markers spanning 21 linkage groups (LG) for sea urchin. Based on this genetic map and phenotyping data for eight economic traits, 33 potentially significant QTLs were detected on ten different LGs with explanations ranging from 9.90% to 46.30%, partly including 10 QTLs for test diameter, six QTLs for body weight and eight QTLs for Aristotle's lantern weight. Moreover, we found a QTL enrichment LG, LG15, gathering QTLs for test diameter, body weight, gonad weight, light orange-yellow color difference (≥E1) and light yellow color difference (≥E2). Among all QTLs, we genotyped four QTLs for test diameter, Aristotle's lantern weight and body weight using High Resolution Melting (HRM) technology. Finally, we used the verified SNP marker (detected using SLAF sequencing) to explore their marker-assisted selection (MAS) breeding application potential and found that SNP-29 associated tightly with body weight and that heterozygous genotype was a dominant genotype, indicating that SNP-29 was a promising marker for MAS.

Identification and expression analysis of a TLR11 family gene in the sea urchin Strongylocentrotus intermedius.

In this study, a homolog of the TLR11 family gene from the sea urchin Strongylocentrotus intermedius (denoted as SiTLR11) was cloned and characterized. The full-length cDNA of SiTLR11 was 2096-bp long, which included 43 bp of 5' untranslated region (UTR), 238 bp of 3' UTR, and a putative open reading frame of 1815 bp encoding a polypeptide of 604 amino acid residues. Representative domains such as leucine-rich repeat (LRR) (residues 108-249) and a cytoplasmic Toll-interleukin-1 receptor (TIR) (residues 429-571) domains were detected in the predicted amino acid sequence of SiTLR11. SiTLR11 transcript was widely distributed in all the tested tissues, including intestine, tube feet, gonad, coelomocytes, and peristomial membrane, with the highest expression level in the coelomocytes and peristomial membrane. After the sea urchin was injected with polyinosinic:polycytidylic acid (PolyI:C), the expression level of SiTLR11 in the coelomocytes increased significantly, reaching 1.96-fold the level of the control at 12 h, but decreased to level below that of control at 24 and 48 h. Injection of peptidoglycan (PGN) also led to increased expression of SiTLR11, which peaked at 12 h, yielding an increase of 2.19-fold compared to the control group, and continued to increase at 24 and 48 h. However, almost no differences in immunological activity were found in the groups challenged with lipopolysaccharides (LPS), Zymosan A (ZOA), or Vibrio fortis compared to the control. Taken together, the results strongly suggested that SiTLR11 was functionally involved in the immune response triggered by double-stranded RNA (dsRNA) viruses and Gram-positive bacteria.