Front. Genet. Frontiers in Genetics Front. Genet. 1664-8021 Frontiers Media S.A. 10.3389/fgene.2021.717538 Genetics Original Research Identification of Sex-Specific Markers Through 2b-RAD Sequencing in the Sea Urchin (Mesocentrotus nudus) Cui Zhouping 1 † Zhang Jian 1 2 † Sun Zhihui 1 * Liu Bingzheng 1 Zhao Chong 1 Chang Yaqing 1 * 1Key Laboratory of Mariculture and Stock Enhancement in North China Sea, Ministry of Agriculture and Rural Affairs, Dalian Ocean University, Dalian, China 2School of Life Science, Liaoning Normal University, Dalian, China

Edited by: Beide Fu, Institute of Hydrobiology (CAS), China

Reviewed by: Juan D. Gaitan-Espitia, The University of Hong Kong, China; Chuangju Li, Chinese Academy of Fishery Sciences (CAFS), China

 *Correspondence: Zhihui Sun, sunzhihui@dlou.edu.cn Yaqing Chang, changlab@hotmail.com

†These authors have contributed equally to this work

This article was submitted to Livestock Genomics, a section of the journal Frontiers in Genetics

 05 08 2021 2021 12 717538 31 05 2021 14 07 2021 Copyright © 2021 Cui, Zhang, Sun, Liu, Zhao and Chang. 2021 Cui, Zhang, Sun, Liu, Zhao and Chang

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.

Sex-specific markers play an important role in revealing sex-determination mechanism. Sea urchin (Mesocentrotus nudus) is an economically important mariculture species in several Asian countries and its gonads are the sole edible parts for people. However, growth rate and immunocompetence differ by sex in this species, sex-specific markers have not been identified, and the sex-determination mechanism of sea urchin remains undetermined. In this study, type IIB endonuclease restriction-site associated DNA sequencing (2b-RAD-seq) and a genome survey of M. nudus were performed, and three female-specific markers and three female heterogametic single nucleotide polymorphism (SNP) loci were identified. We validated these sex-specific markers via PCR amplification in a large number of individuals, including wild and artificially bred populations. Several open reading frames (ORFs) were predicted, although there are no potential genes known for sex determination and sex differentiation within the scaffold in which the sex-specific markers are located. Importantly, the female-specific sequences and female heterozygous SNP loci indicate that a female heterogametic and male homogametic ZW/ZZ sex-determination system should exist in M. nudus. The results provide a solid basis for revealing the sex-determination mechanism of this species, and open up new possibilities for developing sex-control breeding in sea urchin.

 2b-RAD sex-specific marker sea urchin genome survey SNP Introduction

Sexual dimorphism is common in aquaculture animals, and growth rate, immunocompetence, and body size generally differ significantly by sex (Mei and Gui, 2015; Jiang et al., 2017). Sea urchin is widely distributed and commercially harvested all over the world (Inomata et al., 2016), and the gonads are the only edible parts for people. In addition, edible sea urchins are important mariculture species in several Asian countries. This group of animals shows gender-specific differences in immune response, for example, the immunocompetence of female Paracentrotus lividus urchins might be superior to that of males (Arizza et al., 2013). Moreover, the gonadal growth rate in Abatus cavernosus urchins is faster in males than in females (Gil et al., 2020), and the types and the contents of free amino acids, as well as lipid concentrations, are also closely linked to sexual dimorphism in urchins (Osako et al., 2007; Martínez-Pita et al., 2010; Díaz de Vivar et al., 2019). However, most sea urchins lack conspicuous sexually dimorphic phenotypes, making it impossible to distinguish the phenotypic sex of living animals. This dramatically increases the costs and reduces the efficiency of genetic breeding. In addition to commercial value, sea urchins are excellent organisms for the study of early embryonic development (Juliano et al., 2010), interactions among genes (Yuh et al., 1998; Oliveri et al., 2002), evolution (Kober and Bernardi, 2013), fertilization (Chun et al., 2014), and immunity (Branco et al., 2014). Therefore, it is important to reveal the genetic mechanism of sex determination in this species. To do so, it is critical to develop reliable genetic sex-identification markers (Kobayashi et al., 2013; Liu et al., 2013; Mei and Gui, 2015; Chen et al., 2018). A successful system would improve breeding efficiency and make it possible to implement sex-control breeding in these echinoderms (Wang et al., 2009; Martínez et al., 2014).

Several methods have been applied to identify sex-specific molecular markers, including traditional methods and high-throughput sequencing methods. Amplified fragment length polymorphism (AFLP) is a traditional method to score random SNP markers across an entire genome. However, this approach is expensive and requires complex experiments (Griffiths and Orr, 1999; Flachowsky et al., 2001). With the development of next-generation sequencing, several methods have been developed to obtain thousands of SNPs in a single experiment, such as double-digest restriction-site-associated DNA sequencing (dd-RAD-seq), sequence-based genotyping (SBG) and type IIB endonucleases restriction-site-associated DNA sequencing (2b-RAD-seq). Dd-RAD-seq technology can generate number of SNPs allowing for the construction of phylogenetic trees, detection of genetic variation of a population and accurate genetic maps (Zhou et al., 2014; Puchta et al., 2018). Due to the multiplexing a large number of samples, dd-RAD-seq requires the appropriate tools for analysis (Peterson et al., 2012). Using SBG can achieve genome-wide SNP discovery and genotyping of large populations without the information of a reference genome sequence, and this approach had been successfully used in the plants (Truong et al., 2012). Compared to the other traditional RAD-seq methods, 2b-RAD-seq is cost-effective and simple method, and allows for nearly every restriction site in the genome to be screened in parallel (Wang et al., 2012). Some studies have already successfully applied the 2b-RAD approach to explore sex-specific DNA fragments and SNPs in a great diversity of aquatic organisms (Liu et al., 2018; Zhou et al., 2019; Zhu et al., 2021). In addition to quickly identify genetic sex, sex-specific markers also provide insights into the sex-determination system in various animals (Liu et al., 2018). In several species of Scylla crab, sex-specific SNPs have been identified by 2b-RAD-seq, among which female are heterogamety and male are homogamety, while indicating that a ZW/ZZ sex-determination system may exist in the above organisms (Shi et al., 2018). In Mastacembelus armatus, two male-specific 2b-RAD tags and two SNPs were identified through 2b-RAD-seq, and reported that the Tcl gene was a candidate sex-determination gene (Xue et al., 2020). However, few studies have identified sex-specific markers in sea urchins.

Mesocentrotus nudus belongs to the Strongylocentrotidae family of Echinodermata, and is mainly distributed and farmed in coastal areas of several Asian countries, such as China, South Korea, and Japan (Agatsuma, 2020). To further understand the sex-determination mechanism in this species, we used 2b-RAD-seq to identify sex-specific markers. Subsequently, we verified the candidate markers in a large number of samples, including laboratory breeding population, cultured population, and wild populations, and perfect accordance was observed between sex-specific makers and sexual phenotype. Our results can be used to significantly improve breeding efficiency, and provide useful information for understanding the sex-determination mechanisms in sea urchins.

 Materials and Methods Sample Collection and Preparation

A total of 19 sea urchin individuals were collected from the Key Laboratory of Mariculture and Stock Enhancement in North China Sea, Ministry of Agriculture and Rural Affairs, Dalian Ocean University. In addition, aquaculture populations of M. nudus including 95 individuals were purchased from markets in Dalian Haibao Fishery Limited Company, China. Finally, 19 wild sea urchins were fished near Xixiakou, Shandong Province, China. Tube feet tissues from each individual were sampled and stored in absolute ethanol at −20°C for genomic DNA extraction. To determine the gender of each sea urchin, the corresponding gonad tissues were sampled in 4% paraformaldehyde (PFA) for histological examination.

 DNA Extraction

Genomic DNA was extracted following the protocols of a marine animal tissue genomic DNA extraction kit (Tiangen, DP324). The quality and concentration of DNA was determined by 1% agarose gel electrophoresis and using a SimpliNano spectrophotometer (Biochrom, United Kingdom).

 Gonadal Histology and Identifying of Gender

Samples of gonad tissues were fixed in PFA at 4°C, then washed three times with phosphate buffered solution (PBS) at room temperature. After incubating the samples in 30% sucrose for 3 h, they were embedded in optimal cutting temperature (OCT) compound and sectioned at 3 μm. Then they were stained with hematoxylin/eosin. The gender of each sea urchin was determined by observing the gonadal structure with a microscope (Leica DM4B microscope).

 Library Construction and 2b-RAD Sequencing

The genomic DNA of 20 individuals, including 10 males and 10 females, was used for 2b-RAD-seq. The 2b-RAD libraries were constructed at Qingdao OE Biotech Co., Ltd. (Qingdao, China) following previous methods with slight modification (Wang et al., 2012). Briefly, 200 ng genomic DNA from each individual was digested with 1 U BsaXI restriction enzyme, then 10 μL digestion product was ligated with adaptor. Next, adaptor primers were used to amplify the ligation products by PCR. The quality and concentration of each PCR product was determined on an 8% polyacrylamide gel. Purified PCR products were digested with 2U SapI restriction enzyme at 37°C for 30 min, and the digested products were further purified using streptavidin-coated magnetic beads. The supernatant was incubated with T4 DNA ligase at 16°C for 45 min, followed by purification of the ligation products. Barcodes were introduced by PCR with barcode-bearing primers, the PCR products were purified before sequencing. Finally, each library was pooled for sequencing using the Illumina Nova PE150PE platform.

 Data Filtering, Survey Genome, and Screening of Candidate Sex-Specific Markers

Paired-end reads were merged using Pear software (version 0.9.6) (Zhang et al., 2014), and terminal tag sequences were excluded from each output read. Reads with ambiguous bases (N) exceeding 8%, poor-quality reads, and those without restriction sites were excluded from the raw data. Ultimately, clean reads 33 bp in length were obtained. For the genome survey, DNA was extracted from one female and one male M. nudus. Libraries were constructed and sequenced on the Illumina NovaSeq platform (Illumina Inc., United States). Subsequently, 150 bp paired-end reads were clustered to build reference sequences for further locus genotyping using Ustacks software version 1.34 (Catchen et al., 2011). For genome survey, K-mer distribution was estimated by using jellyfish1 (version: 2.2.6)1 with parameters “-m 17 –C.” The genome size and heterozygosity ratio were estimated by the GenomeScope2. To identify candidate sex-specific markers, high-quality 33 bp reads were aligned against the genome survey using SOAP2 software (Li et al., 2009), and sequences with fewer than three reads were precluded (Li et al., 2009). Finally, candidate sex-specific markers for each individual were generated for further analyses.

 Validation of Sex Specific Sequences

According to the flank sequence of sex-specific markers, specific primers were designed using online software2 according to the flank sequences of the sex-specific marker. Then PCR reactions were carried out in a total volume of 20 μL, containing 10 μL 2 X Taq buffer (KangWei), 1 μL forward and reverse primer (10 mM), and 50–100 ng template DNA, adding sterile water to reach the final volume. The PCR amplification program was as follows: pre-denaturation at 94°C for 2 min, then 30 cycles of 94°C for 30 s, 58°C for 30 s, and a final extension of 2 min at 72°C. The PCR products were determined via 1% agarose gel electrophoresis and sequenced by Tsingke Biotechnology (Tianjin, China).

 Results 2b-RAD-Seq

In total, the 2b-RAD-seq produced 137,684,505 raw reads, 76,661,950 for females and 61,022,555 for males. After filtering low-quality reads, 72,125,643 clean reads for females and 56,991,085 clean reads for males were obtained. Based on the read quality and quantity, two female and two male samples were selected as references, clean reads were clustered and assembled. Subsequently, the reads of 20 individuals were aligned with reference samples, with an average alignment rate of 72.04%. After alignment analysis, 2,384,989 tags with an average depth of 36.78 × were obtained from the 20 individuals, and 43,922 polymorphic SNP markers were also screened (Table 1).

Summary of 2b-RAD data obtained from 20 M. nudus individuals.

Name Gender Raw reads Raw bases Clean bases Tags Sequencing depth
F1 Female 7,416,917 7,211,370 6,997,835 119,303 39.37
F2 Female 7,416,917 7,211,370 6,977,732 131,127 36.83
F3 Female 7,416,917 7,211,370 6,964,296 115,465 40.94
F4 Female 7,416,917 7,211,370 6,962,948 116,298 40.65
F5 Female 7,416,917 7,211,370 6,973,281 121,724 37.85
F6 Female 7,915,473 7,696,999 7,465,187 119,963 41.41
F7 Female 7,915,473 7,696,999 7,475,314 131,311 38.88
F8 Female 7,915,473 7,696,999 7,432,458 116,989 42.4
F9 Female 7,915,473 7,696,999 7,440,918 118,280 42.67
F10 Female 7,915,473 7,696,999 7,435,674 123,705 40.27
M1 Male 5,672,850 5,486,113 5,343,733 116,336 31.43
M2 Male 5,672,850 5,486,113 5,341,739 126,894 29.32
M3 Male 5,672,850 5,486,113 5,318,442 111,114 33
M4 Male 5,672,850 5,486,113 5,307,239 111,615 33.18
M5 Male 5,672,850 5,486,113 5,282,468 115,196 30.92
M6 Male 6,531,661 6,290,517 6,072,984 118,336 34.98
M7 Male 6,531,661 6,290,517 6,095,777 127,515 33.51
M8 Male 6,531,661 6,290,517 6,105,512 111,501 37.54
M9 Male 6,531,661 6,290,517 6,077,488 112,770 36.65
M10 Male 6,531,661 6,290,517 6,045,703 119,547 34.15
Average / 6,884,225.3 6,671,249.8 6,455,836.4 119,249.5 36.79
Total / 137,684,505 133,424,995 129116728 2,384,989 /
Genome Survey

Because the 2b-RAD tags were only 27 bp in length, too short to design PCR primers, and the genomic information for M. nudus is unavailable, a male and a female individual were subjected to short-read de novo sequencing. This generated 648,230,824 raw reads from the female; after filtering low-quality reads, 620,114,248 clean reads were acquired with a GC content of 37.06%. Additionally, a total of 669,375,890 raw reads were obtained from male individual, and 639,002,820 clean reads with a GC content of 37.23% were produced after filtering. The estimated genome size using SOAP de novo K-mer module (K = 17) was 626,279,328 bp with high heterozygosity (1.3%). The SRA raw reads were deposited in the GenBank public database with the accession number of PRJNA741812. Using these clean reads, a draft genome of each sex was assembled (1,582,413 scaffolds in female and 5,953,609 scaffolds in male). The longest scaffold was 52,656 bp with N50 lengths of 1,579 bp in the female genome, and 10,758 bp with N50 lengths of 150 bp in the male genome. Moreover, there were 78,550 (4.96%) and 943 (0.01%) scaffolds with lengths of more than 2 kb in the female and male draft genome, respectively (Table 2).

Summery of genome survey dates.

Female Male
Scaffold (length) Scaffold (number) Scaffold (lentgh) Scaffold (lentgh)
Number (>2,000 bp) – 78,550 – 943
Max 52,656 – 10,758 –
N50 1,579 101,006 150 1,882,406
Total 817,576,203 1,582,413 1,076,599,617 5,953,609
Identification of Candidate Sex-Specific Tags and SNPs

To identify putative sex-specific makers, the 2b-RAD tags and SNPs were aligned against the genome survey. After data filtering, we obtained 13 putative female-specific tags, which were present in 10 female samples (Table 3). In addition, 20,201 putative SNP loci were also obtained for further analysis. A smaller p-value of SNP means a higher correlation with sex; hence, we screened the top 10 sex-specific SNP locus with the lowest p-values for further analysis. Notably, the females appeared heterozygous whereas the males appeared homozygous for all 10-candidate sex-specific SNPs (Table 4), these data suggest that a ZW/ZZ sex-determination system exists in M. nudus.

Information of candidate sex-specific tags.

Tag name Ref ID Sequence
Female-tag 1 ref7866-1 AACTCTGGAACCTTCCCTCCATCTTCC
Female-tag 2 ref62548-1 CCAAAAGCTACTTCATCTCCTGCAAGG
Female-tag 3 ref72559-1 CCTTACATTACGCTTTCTCCGACTGGC
Female-tag 4 ref75945-1 CGGCGGAAGACCGGAGCTCCATGTCAT
Female-tag 5 ref78854-1 CTAGATCTAACCTAGTCTCCTATGTTT
Female-tag 6 ref80771-1 CTCCACAGTACTTGTACTCCTTGACTC
Female-tag 7 ref102102-1 GCCTGGAGCACTGAATCTCCGTGCAAA
Female-tag 8 ref106304-1 GGACCTTTCACCTCTTCTCCTGTATGA
Female-tag 9 ref110721-1 GGTAGGAGCACTTTACCTCCCTCCTAA
Female-tag 10 ref129620-1 TATAGTTTTACCTTTTCTCCTTGTTCT
Female-tag 11 ref155669-1 TGTTTGCTCACTGGAGCTCCTTATGGA
Female-tag 12 ref159947-1 TTCAGAAAAACACTAGCTCCCTAAGAC
Female-tag 13 ref161941-1 TTCTAGAAAACTCTCACTCCACGAAAA

Information of top 10 candidate sex-specific SNPs.

SNP name Ref ID SNP position Reference base Alternative base Female Male
1 ref32333 24 A G A/G A/A
2 ref32333 26 T A T/A T/T
3 ref68728 2 C T C/T C/C
4 ref53064 3 C T C/T C/C
5 ref107153 1 G A G/A G/G
6 ref42855 9 A C C/C A/A
7 ref125173 6 T A T/T T/A
8 ref161496 8 G T G/G G/T
9 ref86238 22 T C T/C T/T
10 ref59409 3 T C T/C T/T
Validation of Candidate Sex-Specific Tags and SNPs

To confirm the authenticity of the candidate sex-specific tags and SNPs, primers were designed according to the flanking sequences of the 2b-RAD tags and SNPs from the genome survey (Table 5). Then, PCR reactions were performed using specific primers for 10 female and 9 male laboratory breeding sea urchins. As shown in Figure 1, the PCR products from 3 of the 13 female tags (female tags 3, 7, and 8) amplified one expected sizes band in females but not in males. Next, we purified these PCR products and sequenced them via Sanger sequencing. The nucleotide sequences were consistent with the female genomic sequences, further confirming the female-specific segments. Considering that the other 10 tags amplified products from both sexes (data not shown), these 10 tags were excluded from further study. In addition, the 10 candidate SNP loci were validated via PCR assays using specific primers from 20 individuals, and all PCR products were sequenced. Three of the candidate sex-specific SNPs were confirmed (Figure 2), and the others were excluded.

primers used to verify candidate sex-specific tags and SNPs.

Sex-specific makers Primer ID Primer sequences (5′–3′)
ref72559-1 Female-tag 3-F Female-tag 3-R GGAGTAGTGTCCCAATATCCG ACTACGACCCTGTGTTTGTTT
ref102102-1 Female-tag 7-F Female-tag 7-R TGGCACTTTGGTGACAATACA ACAACCGCGAGTCTTGAAAT
ref106304-1 Female-tag 8-F Female-tag 8-R ACTACCCTACCACAAAAGCA TTCATGCCTGATTCCGGTTG
ref32333 SNP1-F SNP1-R CAAGGTGCTCTTGCGTGTA GGCTAAGCACGGTTATGAGA
ref32333 SNP2-F SNP2-R CAAGGTGCTCTTGCGTGTA GGCTAAGCACGGTTATGAGA
ref68728 SNP3-F SNP3-R TGAACAGAGCAATCCAGCCAT TTTGTGAAGCTGCGTGCAAAT

Validation of three candidate female-specific 2b-RAD-tags in 10 females and 9 males. Label M represents the DL2,000 DNA marker.

Validation of candidate sex-specific SNP loci in 10 females and 9 males. (A) Three female-specific SNP markers in the sequencing chromatograms of the sex-related sequence of M. nudus. (B) Statistics of sexual phenotype and SNP genotype of the 10 females and 9 males.

 Genetic Sex Identification in <italic>M. nudus</italic> Using Sex-Specific Tags

To further validate the accuracy of the three female-specific tags, we performed PCR reactions on 95 sea urchins from breeding population outside the laboratory and wild population. Initially, the gender of each sea urchin was determined by gonadal histology, and 38 females (Figure 3) and 57 males (Figure 4) were identified in breeding population, 10 females and 9 males were identified in wild population (Figure 5). Then the three female-specific tags, including female-tag 3, female-tag 7, female-tag 8, were used to identify the genetic sex of these individuals. As expected, all three female-specific tags were only detected in the females, but not in males (Figures 6–8). These data suggest that the three female-specific markers could be used to identify the genetic sex of M. nudus with a very high rate of accuracy.

Identification of 38 females by histological detection of gonads in the breeding population. NP, nutritive phagocytes; Oo, oogonia.

Identification of 57 males by histological detection of gonads in the breeding population. NP, nutritive phagocytes; SPC, spermatocyte.

Identification of 10 females and 9 males by histological detection of gonads in the wild population. NP, nutritive phagocytes; Oo, oogonia; SPC, spermatocyte.

Electrophoretogram of the amplified products using female-tag 3 for genotypic sex identification in M. nudus. Label M represents the DL2,000 DNA marker (500 and 250 bp DNA ladder are shown).

Electrophoretogram of the amplified products using female-tag 7 for genotypic sex identification in M. nudus. Label M represents the DL2,000 DNA marker (500 and 250 bp DNA ladder are shown).

Electrophoretogram of the amplified products using female-tag 8 for genotypic sex identification in M. nudus. Label M represents the DL2,000 DNA marker (500 and 250 bp DNA ladder are shown).

 Sex-Specific Sequence Annotation

The three female specific tags (female-tag 3, female-tag 7, female-tag 8) locate at scaffold 2838205, 655827, 439353, respectively. The scaffold 2838205 is 717 bp in length, and 655827, 439353 are 4,146 and 5,496 bp, respectively. To better understand these scaffolds, we predicated homologous genes by BLAST program. The databases of Softberry and NCBI, genomes of Strongylocentrotus purpuratus were used for annotation of the scaffolds. Although several open reading frames (ORFs) were predicted, there are no potential genes known for sex determination or differentiation within any scaffold (Figure 9). In addition, potential ORFs fail to be identified in the scaffold 2838205 due to a short sequence.

Annotation of sex-specific sequences obtained in M. nudus. Poly A, poly adenylic acid; CDS, coding sequence; TSS, transcriptional start site. (A) Annotation ofscaffold 655827. (B) Annotation of scaffold 439353.

 Discussion

Mesocentrotus nudus is an economically important species in several Asian countries, but its benefits differ by sex and the gender of animals cannot be easily ascertained. In addition, sex-specific markers have not been identified, and the sex-determination mechanism of sea urchin remains undetermined. According to previous reports, the chromosomes of sea urchins are too small to be analyzed accurately, and heteromorphic sex chromosomes cannot be observed in most species (Saotome, 1987; Saotome et al., 2002; Eno et al., 2009). As an exception, a heteromorphic chromosome has been observed in Paracentrotus lividus sea urchins, for which an XX/XY sex-determination system has also been deduced (Lipani et al., 1996). In many animals, sex determination system has been deduced using sex-specific markers which developed by 2b-RAD-seq. In redtail catfish (Mystus wyckioides), some Y-specific fragments and X-specific fragments have been identified by 2b-RAD-seq, and XX/XY sex determination system was deduced because of the X-specific fragment showing dosage effect between females and males (Zhou et al., 2019). In this study, we revealed that all of the identified sex-specific markers were female-specific, and the identified SNP markers were heterozygous in females but homozygous in males (Figures 1, 2), our data suggest that a ZW/ZZ sex-determination system exists in M. nudus. However, it has been reported that the diploid chromosome number of M. nudus is 42 (2n = 42), and no heteromorphic sex chromosomes had not been observed (Saotome, 1987). In the future, we will use these female-specific sequences for chromosome fluorescence in situ hybridization (FISH) and attempt to locate sex chromosomes. Our results may help elucidate the origin of sex-determination systems.

Although it is difficult to observe heteromorphic sex chromosomes, sex-determining regions have been reported in several aquaculture animals (Wen et al., 2020; Xue et al., 2020). Liu et al. (2018) obtained a male-specific sequence 8,661 bp in length from both bighead carp (Hypophthalmichthys nobilis) and silver carp (Hypophthalmichthys molitrix) via 2b-RAD-seq, and suggested that it was located on the Y chromosome or in a nascent Y-linked region. Similarly, using male-specific 2b-RAD tags, Kirubakaran et al. (2019) identified a potential sex-determining region of 9 kb in Atlantic cod, and suggested that the zkY gene located in this region was a candidate sex-determining gene. In the present study, we identified three female-specific regions, providing a solid basis to reveal the sex-determination mechanisms in sea urchin. After homologous analysis by BLAST, several ORFs were predicted in these female-specific scaffolds, but no potential genes known for sex determination or differentiation within any scaffold. This phenome is similar to most previously reports that even though the sex-specific markers were extended into long segments within hundreds of genes, no sex related genes have been obtained (Pan et al., 2015; Shi et al., 2018; Zhou et al., 2019; Fang et al., 2020). In the future, we will use genome walking to analyze the differences in sequences between male and female individuals. The complete genome should also be sequenced to further elucidate the linkages among the three sex-specific sequences, as well as to reveal the sex-determination mechanism in sea urchin.

The sex-specific markers identified in our study were used to identify gender of more than one hundred of sea urchins with 100% accuracy (Figures 6–8). However, the need for DNA sample preparation limits the efficiency of this method. A direct PCR method would significantly reduce the identification time (Wang et al., 1993; Satya et al., 2013; Gao et al., 2020). This method uses NaOH-incubated products rather than purified DNA as template, and DNA can be extracted in 30 min from the tube feet of sea urchins. Furthermore, we compared the effectiveness of PCR using two different DNA extraction methods (using a DNA extraction kit and the alkaline extraction DNA procedure) and perfect consistency was observed between the two (data not shown). However, most aquaculture farms do not have PCR equipment and other required devices, which limits the application of this method. The loop-mediated isothermal amplification (LAMP) method can quickly amplify DNA with high specificity and efficiency under isothermal conditions, and LAMP equipment is not expensive (Notomi et al., 2000; Nagamine et al., 2002). Therefore, it would be of great value to develop a rapid sex-identification system that combines LAMP technology with female-specific markers.

 Conclusion

In conclusion, we identified sex-specific markers and SNPs via 2b-RAD-seq and a genome survey in M. nudus. These were found to belong to three female-specific scaffolds and were successfully used to accurately identify genetic sex of M. nudus. In addition, a ZW/ZZ sex-determination system was deduced for M. nudus, based on the presence of female-specific markers and the sex-specific SNPs in females appeared heterozygous but homozygous in males. Taken together, our results will help improve breeding efficiency, and provide a powerful tool for uncovering sex-determination mechanisms in sea urchins.

 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: https://www.ncbi.nlm.nih.gov/genbank/, PRJNA741812.

 Ethics Statement

Ethical review and approval was not required for the animal study because sea urchin belongs to invertebrates, all animal experiments in this study are carried out with the permission of the local government.

 Author Contributions

ZS and YC: conceptualization and resources. ZC and JZ: methodology. ZC, JZ, and BL: investigation. ZC and CZ: data curation. ZS, ZC, and JZ: writing—original draft preparation. ZS: writing—review and editing, funding acquisition. All authors have read and agreed to the published version of the manuscript.

 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.

Funding. This research was funded by the National Natural Science Foundation of China (Grant No. 31802276).

 References Agatsuma Y. (2020). “Mesocentrotus nudus,” in Developments in Aquaculture and Fisheries Science, ed. Lawrence J. M. (Amsterdam: Elsevier Press), 627–641. 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. A Mol. Integr. Physiol. 164 447–455. 10.1016/j.cbpa.2012.11.021 23220062 Branco P. C. Figueiredo D. A. L. da Silva J. R. M. C. (2014). New insights into innate immune system of sea urchin: coelomocytes as biosensors for environmental stress. OA Biology 18:2. Catchen J. M. Amores A. Hohenlohe P. Cresko W. Postlethwait J. H. (2011). Stacks: building and genotyping loci de novo from short-read sequences. G3 1 171–182. 10.1534/g3.111.000240 22384329 Chen J. Fan Z. Tan D. Jiang D. Wang D. (2018). A review of genetic advances related to sex control and manipulation in tilapia. J. World Aquac. Soc. 49 277–291. 10.1111/jwas.12479 Chun J. T. Limatola N. Vasilev F. Santella L. (2014). Early events of fertilization in sea urchin eggs are sensitive to actin-binding organic molecules. Biochem. Biophys. Res. Commun. 450 1166–1174. 10.1016/j.bbrc.2014.06.057 24960199 Díaz de Vivar M. E. Zárate E. V. Rubilar T. Epherra L. Avaro M. G. Sewell M. A. (2019). Lipid and fatty acid profiles of gametes and spawned gonads of Arbacia dufresnii (Echinodermata: Echinoidea): sexual differences in the lipids of nutritive phagocytes. Mar. Biol. 166:96. 10.1007/s00227-019-3544-y Eno C. C. Bottger S. A. Walker C. W. (2009). Methods for karyotyping and for localization of developmentally relevant genes on the chromosomes of the purple sea urchin, Strongylocentrotus purpuratus. Biol. Bull. 217 306–312. 10.1086/BBLv217n3p306 20040754 Fang S. Zhang Y. Shi X. Zheng H. Li S. Zhang Y. (2020). Identification of male-specific SNP markers and development of PCR-based genetic sex identification technique in crucifix crab (Charybdis feriatus) with implication of an XX/XY sex determination system. Genomics 112 404–411. 10.1016/j.ygeno.2019.03.003 30851358 Flachowsky H. Schumann E. Weber W. E. Peil A. (2001). Application of AFLP for the detection of sex-specific markers in hemp. Plant Breed. 120 305–309. 10.1046/j.1439-0523.2001.00620.x Gao F. X. Shi Y. Duan W. Lu W. J. Huang W. Zhang X. J. (2020). A rapid and reliable method for identifying genetic sex in obscure pufferfish (Takifugu obscurus). Aquaculture 519:734749. 10.1016/j.aquaculture.2019.734749 Gil D. G. Zaixso H. E. Tolosano J. A. (2020). Sex-specific differences in gonopore and gonadal growth trajectories in the brooding sea urchin, Abatus cavernosus (Spatangoida). Invertebr. Biol. 139:e12278. 10.1111/ivb.12278 Griffiths R. Orr K. (1999). The use of amplified fragment length polymorphism (AFLP) in the isolation of sex-specific markers. Mol. Ecol. 8 671–674. 10.1046/j.1365-294x.1999.00578.x 10327660 Inomata E. Murata Y. Matsui T. Agatsuma Y. (2016). Gonadal production and quality in the sea urchin Mesocentrotus nudus fed a high-protein concentrated red alga Pyropia yezoensis. Aquaculture 454 184–191. 10.1016/j.aquaculture.2015.12.003 Jiang J. Zhou Z. Dong Y. Gao S. Sun H. Chen Z. (2017). Comparative analysis of immunocompetence between females and males in the sea cucumber Apostichopus japonicus. Fish Shellfish Immunol. 63 438–443. 10.1016/j.fsi.2017.02.038 28238861 Juliano C. E. Yajima M. Wessel G. M. (2010). Nanos functions to maintain the fate of the small micromere lineage in the sea urchin embryo. Dev. Biol. 337 220–232. 10.1016/j.ydbio.2009.10.030 19878662 Kirubakaran T. G. Andersen Ø De Rosa M. C. Andersstuen T. Hallan K. Kent M. P. (2019). Characterization of a male specific region containing a candidate sex determining gene in Atlantic cod. Sci. Rep. 9 1–9. 10.1038/s41598-018-36748-8 30644412 Kobayashi Y. Nagahama Y. Nakamura M. (2013). Diversity and plasticity of sex determination and differentiation in fishes. Sex Dev. 7 115–125. 10.1159/000342009 22948719 Kober K. M. Bernardi G. (2013). Phylogenomics of strongylocentrotid sea urchins. BMC Evol. Biol. 13:88. 10.1186/1471-2148-13-88 23617542 Li R. Yu C. Li Y. Lam T. Yiu S. Kristiansen K. (2009). SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics 25 1966–1967. 10.1093/bioinformatics/btp336 19497933 Lipani C. Vitturi R. Sconzo G. Barbata G. (1996). Karyotype analysis of the sea urchin Paracentrotus lividus (Echinodermata): evidence for a heteromorphic chromosome sex mechanism. Mar Biol. 127 67–72. 10.1007/BF00993645 Liu H. Guan B. Xu J. Hou C. Tian H. Chen H. (2013). Genetic manipulation of sex ratio for the large-scale breeding of YY super-male and XY all-male yellow catfish (Pelteobagrus fulvidraco (Richardson)). Mar. Biotechnol. 15 321–328. 10.1007/s10126-012-9487-7 23053056 Liu H. Pang M. Yu X. Zhou Y. Tong J. Fu B. (2018). Sex-specific markers developed by next-generation sequencing confirmed an XX/XY sex determination system in bighead carp (Hypophthalmichehys nobilis) and silvercarp (Hypophthalmichthys molitrix). DNA Res. 25 257–264. 10.1093/dnares/dsx054 29315393 Martínez P. Viñas A. M. Sánchez L. Díaz N. Ribas L. Piferrer F. (2014). Genetic architecture of sex determination in fish: applications to sex ratio control in aquaculture. Front. Genet. 5:340. 10.3389/fgene.2014.00340 25324858 Martínez-Pita I. García F. J. Pita M. L. (2010). Males and females gonad fatty acids of the sea urchins Paracentrotus lividus and Arbacia lixula (Echinodermata). Helgol. Mar. Res. 64 135–142. 10.1007/s10152-009-0174-7 Mei J. Gui J. F. (2015). Genetic basis and biotechnological manipulation of sexual dimorphism and sex determination in fish. Sci. China Life Sci. 58 124–136. 10.1007/s11427-014-4797-9 25563981 Nagamine K. Hase T. Notomi T. (2002). Accelerated reaction by loop-mediated isothermal amplification using loop primers. Mol. Cell. Probes 16 223–229. 10.1006/mcpr.2002.0415 12144774 Notomi T. Okayama H. Masubuchi H. Yonekawa T. Watanabe K. Amino N. (2000). Loop-mediated isothermal amplification of DNA. Nucleic Acids Res. 28:e63. 10.1093/nar/28.12.e63 10871386 Oliveri P. Carrick D. M. Davidson E. H. (2002). A regulatory gene network that directs micromere specification in the sea urchin embryo. Dev. Biol. 246 209–228. 10.1006/dbio.2002.0627 12027443 Osako K. Fujii A. Ruttanapornvareesakul Y. Nagano N. Kuwahara K. Okamoto A. (2007). Differences in free amino acid composition between testis and ovary of sea urchin Anthocidaris crassispina during gonadal development. Fish. Sci. 73 660–667. 10.1111/j.1444-2906.2007.01379.x Pan Z. J. Li X. Y. Zhou F. J. Qiang X. G. Gui J. F. (2015). Identification of sex-specific markers reveals male heterogametic sex determination in Pseudobagrus ussuriensis. Mar. Biotechnol. 17 441–451. 10.1007/s10126-015-9631-2 25981673 Peterson B. K. Weber J. N. Kay E. H. Fisher H. S. Hoekstra H. E. (2012). Double digest RADseq: an inexpensive method for de novo SNP discovery and genotyping in model and non-model species. PLoS One 7:e37135. 10.1371/journal.pone.0037135 22675423 Puchta M. Bolc P. Piechota U. (2018). Review of genome sampling methods in sequencing libraries preparation protocols. Agron. Sci. 73 93–105. 10.24326/asx.2018.4.9 Saotome K. (1987). Chromosome numbers in 8 Japanese species of sea urchins. Zoo Sci. 4 483–487. Saotome K. Kamimura R. Kurokura H. Hirano R. (2002). Male chromosomes of sea urchin hybrid andromerogones created with cryopreserved sperm. Zoo Sci. 19 185–189. 10.2108/zsj.19.185 12012781 Satya P. Mitra S. Ray D. P. Mahapatra B. S. Karan M. Jana S. (2013). Rapid and inexpensive NaOH based direct PCR for amplification of nuclear and organelle DNA from ramie (Boehmeria nivea), a bast fibre crop containing complex polysaccharides. Ind. Crops Prod. 50 532–536. 10.1016/j.indcrop.2013.07.049 Shi X. Waiho K. Li X. Ikhwanuddin M. Miao G. Lin F. (2018). Female-specific SNP markers provide insights into a WZ/ZZ sex determination system for mud crabs Scylla paramamosain, S. tranquebarica and S. serrata with a rapid method for genetic sex identification. BMC Genom. 19:981. 10.1186/s12864-018-5380-8 30594128 Truong H. T. Ramos A. M. Yalcin F. de Ruiter M. van der Poel H. J. Huvenaars K. H. (2012). Sequence-based genotyping for marker discovery and co-dominant scoring in germplasm and populations. PLoS One 7:e37565. 10.1371/journal.pone.0037565 22662172 Wang D. Mao H. L. Chen H. X. Liu H. Q. Gui J. F. (2009). Isolation of Y- and X-linked SCAR markers in yellow catfish and application in the production of all-male populations. Anim. Genet. 40 978–981. 10.1111/j.1365-2052.2009.01941.x 19694649 Wang H. Qi M. Cutler A. J. (1993). A simple method of preparing plant samples for PCR. Nucleic Acids Res. 21:4153. 10.1093/nar/21.17.4153 8371994 Wang S. Meyer E. McKay J. K. Matz M. V. (2012). 2b-RAD: a simple and flexible method for genome-wide genotyping. Nat. Methods 9 808–810. 10.1038/nmeth.2023 22609625 Wen M. Feron R. Pan Q. Guguin J. Jouanno E. Herpin A. (2020). Sex chromosome and sex locus characterization in goldfish, Carassius auratus (Linnaeus, 1758). BMC Genom. 21:552. 10.1186/s12864-020-06959-3 32781981 Xue L. Z. Guo X. F. Zhou Y. L. Wang Z. W. Fan H. P. Li D. P. (2020). Screening and characterization of sex-specific markers by 2b-RAD sequencing in zig-zag eel (Mastacembelus armatus) with implication of XY sex determination system. Aquaculture 528:735550. 10.1016/j.aquaculture.2020.735550 Yuh C. H. Bolouri H. Davidson E. H. (1998). Genomic cis-regulatory logic: experimental and computational analysis of a sea urchin gene. Science 279 1896–1902. 10.1126/science.279.5358.1896 9506933 Zhang J. Kobert K. Flouri T. Stamatakis A. (2014). PEAR: a fast and accurate Illumina Paired-End reAd mergeR. Bioinformatics 30 614–620. 10.1093/bioinformatics/btt593 24142950 Zhou X. Xia Y. Ren X. Chen Y. Huang L. Huang S. (2014). Construction of a SNP-based genetic linkage map in cultivated peanut based on large scale marker development using next-generation double-digest restriction-site-associated DNA sequencing (ddRADseq). BMC Genomics 15:351. 10.1186/1471-2164-15-351 24885639 Zhou Y. Wu J. Wang Z. Li G. Mei J. Zhou L. (2019). Identification of sex-specific markers and heterogametic XX/XY sex determination system by 2b-RAD sequencing in redtail catfish (Mystus wyckioides). Aquac. Res. 50 2251–2266. 10.1111/are.14106 Zhu C. Liu H. Cheng L. Pan Z. Chang G. Wu N. (2021). Identification of sex-specific sequences through 2b-RAD sequencing in Pseudobagrus ussuriensis. Aquaculture 539:736639. 10.1016/j.aquaculture.2021.736639

http://www.genome.umd.edu/jellyfish.html

http://biotools.nubic.northwestern.edu/OligoCalc.html