Mutations in Growth-Related Genes Induced by EMS Treatment in Scallops

Background: The goal of genetic breeding is to select variants with mutations that are related to expected traits, such as fast growth. Artificial induction has been widely used to obtain strains with more mutations for further selection. Ethylmethylsulfone (EMS) is one of the most commonly used chemical mutagens in plant and microorganism breeding. However, the application of EMS mutagenesis in shellfish has not been reported. The aim of this study is to evaluate the potential use of EMS as a mutagen in scallop breeding, especially in characterization of mutations in growth-related genes.

Results: Our results indicated that hatching of about 50% of fertilized eggs was blocked by treatment with 20 mM EMS for 3 h and the resulted larvae developed normally into adult stages. We then evaluated the mutagenic effects of EMS by sequencing the genomes of 4 adult scallops from the control group and 12 from the treatment group at 8 months after fertilization. On average, after removing shared types of mutations, there were 1,151,380 ± 258,188 SNPs (Single Nucleotide Polymorphisms) and 229,256 ± 51,714 InDels (insertion-deletion) in each animal in the EMS treatment group, while there were only134841 ± 10,115 SNPs and 42,605 ± 5,136 InDels in the control group. The average mutation rate in the genome of the EMS treatment group (0.0137 ± 0.0013%) was about 9 times that of the control group (0.0015 ± 0.0002%). GO (Gene Ontology) annotation and KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analyses revealed that mutations induced by EMS occurred evenly in most biological processes, cellular components and functions, as well in most pathways. However, significant lower percentage of mutations were found in the exonic region, in non-synonymous or Stopgain/Stoploss SNPs and in coding domains, suggesting apparent DNA repair or selection during grow-out stage. Analyses of the growth-related genes with mutations indicated that mutations in MFS (Major Facilitator Superfamily) and Tubulin were only found in the large-sized group (Five largest scallops: Treated-1, Treated-2, Treated-3, Treated-4, and Treated-5) and Homeobox and Socs (Suppressor of cytokine signaling) only in the small group (Two smallest scallops: Treated-11 and Treated-12). These results suggested that these genes may be involved in the regulation of growth in these animals, although further verification is certainly warranted.

Conclusion: Treatment of fertilized eggs with 20 mM EMS for 3 h induced 9 times more mutations in scallop genomes. We found that mutations in MFS and Tubulin may be related to fast growth in the large-sized group and those mutations in Homeobox and SOCs may be involved in the slow growth in the small-sized scallops. EMS can be used to accelerate selection of economically important traits in molluscs.

Background

The purpose of genetic breeding in aquaculture is to select variants with superior economically important traits, such as rapid growth, disease resistance and high survival. No matter which breeding method is used, successful selection depends on the availability of mutations in the base population. Mutations are the primary agents of long-term evolution. Genetic variation and naturally occurring mutations play an important role in evolution (Hou et al., 2019). Although naturally occurring mutations represent a major force driving evolution, the occurrence rate of spontaneous mutation is extremely low (Lian et al., 2020). Studies have shown that the natural mutation rate of specific loci in fish is generally lower than 1.0 × 10^–6 (Kuroyanagi et al., 2013). The effective way to obtain a high mutation rate is by artificial induction. Artificially induced genetic variations represent important supplementary variation in plant breeding programs complementary to sources from natural origins (Lu et al., 2016). In recent years, artificial mutagenesis has been widely used in animal and plant breeding as well as in microbe breeding. For example, application of physical and chemical mutagenesis has led to great progress in the selection of algal strains with increased lipid production (Ornuma et al., 2018).

Physical and chemical approaches are the most commonly used artificial induction methods in the construction of plant mutant library (Chen et al., 2018). Compared with physical mutagenesis, chemical mutagenesis is more cost-effective and does not require complicated equipment. Ethylmethylsulfone (EMS) is one of the most commonly used chemical mutagens in breeding of plants such as wheat and Maize (Tran et al., 2020; Xiong et al., 2020). By alkylation of guanine, EMS tends to induce more point mutations from G:C to T:A (Till et al., 2007; Sikora et al., 2011). EMS has been used in the selection of new strains with increased lipid productivity in Chlorella sp. (Ornuma et al., 2018) and also new strain with elevated tolerance to salt in wheats (Johanna et al., 2020). In recent years, EMS chemical mutagenesis has also begun to be applied in animals such as Drosophila (Lee et al., 2016) and nematodes (Song et al., 2014). However, EMS has not been used in genetic breeding of molluscs.

The bay scallop Argopecten irradians is a commercially important bivalve widely cultured in northern China (Ning et al., 2018). The cultured population of the bay scallops originated from only a small number of brood stocks survived from the last attempt of 3 consecutive introductions in early 1980s (Zhang et al., 1986). Efforts have been devoted to improve the stocks through re-introduction of more brood stocks, hybridization and selective breeding. Selective breeding has been proven to be an effective approach for selection of variants with superior traits, mainly from populations with naturally occurring mutations (Wang et al., 2018). However, selective breeding based on naturally occurred mutations is limited by decreasing genetic diversity caused by continuous inbreeding and thus lack of selectable mutants (Zhao et al., 2019). Therefore, it is necessary to broaden the genetic diversity of the bay scallops by artificial induction. In the present study, we aimed to evaluate the feasibility of using EMS as a mutagen in the genetic breeding of scallops and characterize the mutations in growth-related genes.

Materials and Methods

Animals

The animals used in this study were ‘Bohai Red’ scallops, a new bay scallop strain selected from the hybrids between the bay scallop and the Peruvian scallop, which is widely cultured in northern China (Wang et al., 2017). The scallops were cultured in a scallop farm located in Yangma Island, Mouping, Shandong. They were conditioned to mature in the scallop hatchery of Yantai Spring-Sea AquaSeed Co., Ltd. located in Laizhou, Shandong.

Determination of 50% Lethal Concentration of EMS Treatment

When the scallops became mature, they were induced to spawn individually in beakers of 1 L by exposing them to air for 30 min followed by a temperature shock from 18 to 23°C. The spawning scallops were watched carefully to collect eggs and sperm separately. Eggs or sperm from different scallops were pooled together. Then sperm were mixed with eggs to obtain fertilized eggs for subsequent experiments.

After fertilization, fertilized eggs were divided into six groups and a certain amount of EMS stock solution was added to the fertilized egg suspension so that the concentrations of EMS in the containers were 0, 10, 20, 30, 40, 50, and 60 mM. Three replicates were set for each group. EMS stock solution was purchased from Sigma (Ronkonkoma, United States). The fertilized eggs were exposed to EMS solutions of different concentrations for 3 h before they were filtered and washed with fresh seawater and kept in water bath for hatching at 23°C. The fertilized eggs were stirred every half hour.

Twenty-eight hours after fertilization, three 1 ml samples were taken from each container and observed under a microscope. Hatching rates, the percentages of D-formed larvae in total fertilized eggs, were determined for each container. The 50% lethal concentration of EMS was determined based on the hatching rates of fertilized eggs exposed to different concentrations of EMS.

Animal Rearing and Sampling

The shell height and shell length of the larvae in each EMS concentration group were measured on Day 3, 6, 9, and 10 after fertilization to examine the effects of EMS on larval growth. Five 1 ml samples were randomly taken from each container for measurements. Ten days after fertilization, when over 50% of the larvae developed eyespots, the larvae were set on plastic collectors and allowed for metamorphosis. After metamorphosis, the set juveniles were taken first to the ponds and then to the open sea for nursery. When the juveniles reached 2 cm in shell height in late Aug., they were dispersed into lantern nets for grow-out.

About 8 months after fertilization, about 20,000 scallops were harvested and screened for individuals with apparent phenotypic mutation. As large phenotypic variations existed in shell height, shell length and shell width, we used the sum of shell height, shell length and shell width to select the largest and smallest individuals. As a result, 5 largest (Treated-1, Treated-2, Treated-3, Treated-4, and Treated-5), 2 smallest scallops (Treated-11 and Treated-12) and 5 normal-sized (Treated-6, Treated-7, Treated-8, Treated-9, and Treated-10) were selected from the two ends and middle of the normal distribution of individual shell lengths for whole genomic resequencing. Four scallops from the untreated control group were also randomly selected for whole genome sequencing. The sizes of the sequenced scallops were given in Table 1.

TABLE 1

TABLE 1. Sizes of the sequenced individuals in the control and treatment group.

Whole Genome Sequencing

After harvest, the scallops were dissected to collect their adductor muscles. DNA was extracted from each adductor muscle sample. And the qualified DNA was randomly broken into 350bp fragments by Covaris crusher, and NDM607-01 was used to build the database. The whole library was prepared by terminal repair, ployA tail addition, sequencing connector addition, purification and PCR amplification. After the library construction was completed, Qubit2.0 was used for preliminary quantification, and the library was diluted to 1 ng/μl. Agilent 2100 was then used to test the insert size of the library, and QPCR was used to accurately quantify the effective concentration of the library (effective concentration of the library >2 nM). And sequenced using paired-end Illumina platform.

Evaluation of EMS-Induced Mutations

We aligned the raw data to the Bay scallop reference genome (BioSample accessions: SAMN08322131; Liu et al., 2020 using the mem algorithm in BWA (Burrows-Wheeler Alignment) (Li and Durbin, 2009). We then used Samtools to convert the file format and sort the subsequent alignment sequences (Li et al., 2009). Subsequently, GATK (The Genome Analysis Toolkit) was used to remove repeat sequences and to detect SNPs and InDels. SNPs and InDels specific to each group were obtained by removing shared mutations identified by comparing their mutation information (including CHROM, POS, REF, ALT, and GENO) in both the control and treatment group. Then we annotated each variation with annovar to obtain information for further identification of synonymous and nonsynonymous mutations. Finally, SNPs/InDels were annotated to the reference genome using the annovar package (Alexander et al., 2009; Vilella et al., 2009; Wang et al., 2010; Danecek et al., 2011; Yang et al., 2011; Schiffels and Durbin, 2014; Zhang et al., 2018). For GO and KEGG enrichment analysis, the final results were obtained by performing enrichment analysis using hypergeometric test, followed by FDR (False Discovery Rate) correction of the p-values.

Screening of Mutations in Growth-Related Genes

We compared the genes with nonsynonymous SNP sites in the large- and small-sized individuals against controls and normal-sized individuals to find the mutated genes that appeared in only large- or small-sized individuals. From these mutated genes, we then searched for genes that are potentially related to growth.

Statistics and Data Analyses

Statistical analyses of the data were performed using the SPSS (Statistical Product and Service Solutions) statistical package (SPSS Science, Chicago, Illinois). Briefly, the raw data on size parameters in different treated group and hatching rate were first tested for normality of distribution using the Kolmogorov-Smirnov and Shapiro-Wilk tests, and for homogeneity of the variance using the Levene tests. If the raw data failed these tests, square root or natural logarithm transformations were then performed. When both the normality and homogeneity were met, the differences in growth and hatching rate were analyzed with one-way analysis of variance (ANOVA) followed by multiple comparison tests (Tukey’s test) (Wang et al., 2010).

Results

Determination of the 50% Lethal Concentration of EMS Treatment

To figure out the appropriate dose for EMS treatment, we first determined the 50% lethal concentration of EMS by examining the hatching rates of fertilized eggs treated with different concentration of EMS. Our results showed that EMS inhibited hatching rate in a dose-dependent manner. In the treatment groups, some fertilized eggs either developed into abnormal embryos or ceased to develop at early stages (Figure 1). Compared with the control group, about 50% of the fertilized eggs was blocked from hatching by 20 mM EMS (Figure 2). Therefore, we chose 20 mM EMS in subsequent experiments.

FIGURE 1

FIGURE 1. Hatching rate of fertilized eggs treated with different concentrations of EMS.

FIGURE 2

FIGURE 2. D-formed larvae at 28 h after fertilization in the control (A) and the EMS treatment group (B).

Effects of EMS Treatment on Larval Growth

EMS exhibited slightly inhibitory effects on larval growth. As can be seen in Figure 3, the growth curve of shell height and shell length of the treatment group (20 mM EMS) were significantly different from that of the control group during larval stage (p < 0.05, ANOVA). Shell height on day 6 and day 10 and shell length on Day 6, 9, and 10 in the control group were significantly larger than those of the EMS treatment group (p < 0.05, t-test).

FIGURE 3

FIGURE 3. Effects of EMS on growth in shell height (A) and shell width (B) of scallop larvae.

Genomic Analyses of EMS-Induced Mutations

To identify the mutations induced by EMS that may have contributed to important traits such as growth, we re-sequenced the genomes of the largest and smallest individuals from the treatment group, as well as those of the control group. In total, we obtained 156.9 Gb clean data from these 16 animals. The Q30 values of all samples were higher than 80%, indicating that the data were of low error rate and high quality (Table 2). The average mapping rate to the reference genome sequence was 99.54% with an average sequencing depth of 11.19× (Table 3).

TABLE 2

TABLE 2. Genome sequencing data.

TABLE 3

TABLE 3. Summary of assembly results.

SNPs

The average number of SNPs in the control group (5,866,344 ± 95,080) is not significantly different from that of the EMS treatment group (5,626,416 ± 297,062) (Table 4). However, after removing the shared SNPs with the control group, significantly more SNPs were found in the EMS treatment group than the control group (p < 0.05, t-test). As can be seen from Table 4, the number of SNPs in the control group was 134,841 ± 10,115 while that in the EMS treatment group was 1,151,380 ± 258,188. The SNPs are mainly located in intergenic regions, followed by intronic regions, and the rest were located in upstream and exon region. The Ts/Tv (Transitions/Transversions) value in the EMS treatment group was significantly higher than that in the control group (p < 0.05) (Table 5). (Transitions means a transformation between nucleotides of the same type; Transversions means a transformation between nucleotides of the different type).

TABLE 4

TABLE 4. SNP information in the control and EMS treatment group after calling.

TABLE 5

TABLE 5. SNP number in the control and EMS treated group after removing shared types with controls.

The results also showed that the average number of non-synonymous SNP in the treatment group (17,020 ± 3,449) were much higher than those in the controls (2471 ± 275). The average number of Stopgain/Stoploss SNPs in the treatment group (358 ± 60) was also higher than that of the control group (64 ± 9) (Table 6).

TABLE 6

TABLE 6. SNP mutation types in the control and EMS treatment groups.

GO enrichment analysis for the genes with SNPs located in the upstream and exon regions revealed that genes were significantly enriched in the GO terms of protein binding, cell-cell adhesion via plasma-membrane adhesion molecules, homophilic-cell adhesion via plasma-membrane adhesion molecules, binding and cell-cell adhesion in the EMS treatment group, but the genes were only enriched in GO term of protein binding in the control group (Figure 4).

FIGURE 4

FIGURE 4. (A) GO annotation analysis of SNPs located in upstream and exon region in the control group. (B) GO annotation analysis of SNPs located in upstream and exon region in the EMS treatment group.

GO and KEGG analyses of the genes with mutations showed that there were no distinct distribution for the GO terms that were classified into the biological processes, cellular components and molecular functions between the EMS-treated group and control group, as well as for the KEGG pathways (Figure 5).

FIGURE 5

FIGURE 5. (A) GO annotation analysis of nonsynonymous SNPs in the control group. (B) GO annotation analysis of nonsynonymous SNPs in the treatment group. (C) KEGG enrichment analysis of nonsynonymous SNPs in the control group. (D) KEGG enrichment analysis of nonsynonymous SNPs in the treatment group.

InDels

The average numbers of InDels in the genome and in the CDs of the control group (1,548,571 ± 28,720 in the genome and 9,016 ± 310 in the CDs) was also not significantly different from that of the EMS treatment group (1,462,511 ± 83,433 in the genome and 8,518 ± 552 in the CDs) (Table 7). After removing shared InDels in the controls, the average numbers of InDels in the genome and in the CDs of the EMS treatment group (229,256 ± 51,714 in the genome and 1798 ± 302 in the CDs) were much higher than those in the control group (42,605 ± 5,136 in the genome and 640 ± 15 in the CDs) (Table 8).

TABLE 7

TABLE 7. InDels in the control and EMS treatment group after calling.

TABLE 8

TABLE 8. InDels in the control and EMS treated group after removing shared types with controls.

Go and KEGG analysis of InDels showed that the number of genes enriched in the treated group was almost twice that of the control group. And mutations were evenly distributed across nearly all biological processes, cellular components and molecular functions, and nearly all pathways (Figure 6).

FIGURE 6

FIGURE 6. (A) GO annotation analysis of InDels in the control group. (B) GO annotation analysis of Indels in the treatment group. (C) KEGG enrichment analysis of Indels in the control group. (D) KEGG enrichment analysis of Indels in the treatment group.

Mutation Rates Induced by EMS Treatment

In this study, statistics were performed on data after removing mutations shared between control and treated groups. we have found that the average mutation rate as revealed from SNPs of the control group was 0.00115 ± 0.00015%, while that of the EMS treatment group was 0.01144 ± 0.00222%. And the average mutation rate as revealed from InDels of the control group was 0.000365 ± 0.000064%, while that of the EMS treatment group was 0.002281 ± 0.000473%. Taken together, the average mutation rate in the genome of the EMS treatment group (0.013721 ± 0.001347%) was much higher than that in the control group (0.001518 ± 0.000215%) (p < 0.05). Our results showed that the average mutation rate of the treatment group was about 9 times that of the control group.

We also found that the mutations in the treatment group tend to occur in non-coding or intronic regions, resulting in synonymous mutations or mutations with fewer Stopgain/Stoploss. The percentage of exonic SNPs in the total SNPs in the EMS treatment group (3.47 ± 0.12%) was significantly lower than that of the control group (3.63 ± 0.13%; p < 0.05) while the percentage of intronic SNPs in the total SNPs was not different between them (p > 0.05) (Table 5). The percentage of non-synonymous SNPs in the EMS treatment group (1.49 ± 0.06%) was significantly lower than that of the control group (1.83 ± 0.10%) while the percentage of synonymous SNPs in the EMS treatment group (1.93 ± 0.06%) was significantly higher than that of the control group (1.73 ± 0.03%) (p < 0.001 for both comparisons) (Table 6). Similarly, the percentage of Stopgain/Stoploss SNPs in the treatment group (0.032 ± 0.002%) was also significantly lower than that of the control group (0.048 ± 0.006%) (p < 0.001; Table 7). The percentage of InDels in the coding domains out of the total InDels in the whole genome in the treatment group (0.81 ± 0.16%) was much lower than that of the control group (1.52 ± 0.15%) (p < 0.001; Table 8).

Mutations in Growth-Related Genes

To find the potential mutations that may be involved in growth regulation, we attempted to screen the mutations in growth-related genes that existed in only the large-sized animals or small-sized animals but not in the control group that were not exposed to EMS or normal-sized animals that were treated with EMS. We made a matrix using the information of nonsynonymous SNPs by first filtering for genes that contained only nonsynonymous SNP sites in the treated group and then further filtering for genes that contained only nonsynonymous SNP sites in the large individual group as well as only in the small individual group. Our results showed that mutations in major facilitator superfamily (MFS) and Tubulin were found only in the large-sized group while mutations in Homeobox and Socs were found only in the small-sized group.

Discussion

Effects of EMS Treatment on Hatching and Larval Growth

In this study, we found that EMS had inhibitory effects on hatching and larval growth. Similar results have also been observed in previous studies (Hou et al., 2019). We speculate that the mutation induced by EMS may affect the transcription of mRNA and disorder the expression of protein, consistent with the principle of EMS mutagenesis.

Feasibility of Using EMS as a Mutagen in Scallops

As naturally occurring mutations are often limited for genetic breeding programs, artificial mutation induction has become a potent method for the development of novel germplasm (Penna et al., 2012). In this study, we examined the feasibility of using EMS as a mutagen in scallops breeding. We demonstrated that by treating fertilized eggs at 20 mM EMS for 3 h, we were able to induce 9 times more mutations than naturally occurring mutations. Of these mutations, about 3.47% occurred in the coding region and 1.48% were nonsynonymous mutations. These mutations may cause alteration in amino acid sequences in genes that are related to important traits and thus provide rich variations for subsequent selective breeding. Since only the mutations in the gonads may be inherited to offspring, screening of mutated individuals with remarkable phenotypes (or expected traits) in the next generation is thus needed for breeding. The mutated individuals may also provide new opportunities for the exploration of the genes related to the mutated traits. This is demonstrated by the subsequent analyses of growth-related genes utilizing the largest and smallest individuals in the treatment group, as discussed later.

Features of EMS-Induced Mutations in Scallops

In this study, EMS induced extensive variations in the genome, with a high mutation rate of 0.014 ± 0.001%. Mutations in the format of both SNPs and InDels were induced in the mutated individuals. SNPs accounted for the majority of mutations, indicating that point mutations were the main type of variation induced by EMS treatment. Similar results have been found in Oryzias latipes induced by ENU (N-ethyl-N-nitrosourea), another commonly used chemical mutagen (Yoshihito et al., 2006).

After removing the shared SNPs of the controls, we found that the Ts/Tv values of both the treatment group and the control group were slightly larger than 1.0, indicating more transitions than transversions. The cause for more transitions than transversions may lie in that fact that 5-methylcytosine residues are the most prone mutation sites in the genome and can be easily deaminated to form thymine. The Ts/Tv values in both the control and the treatment group were significantly higher than the theoretical value of Ts/Tv, which is 0.5, a phenomenon named conversion deviation. It was believed that the formation of conversion deviation may be a way to reduce harmful mutations due to the selection in the long-term evolutionary process (Wakeley, 1996). Therefore, we speculate that the individuals with increased Ts/Tv in the treatment group may had a better chance of survival than those with lower Ts/Tv.

Our results also showed that there were more nonsynonymous than synonymous SNPs in the control group. In the contrary, more synonymous than nonsynonymous SNPs were found in the EMS treatment group. Nonsynonymous mutations are thought to be largely deleterious due to their property of changing amino acids. For example, it has been found that the ratio of nonsynonymous mutation/synonymous mutation involved in cancer related genes was significantly lower than that of normal genes (Chu and Wei, 2019). It is thus possible that individuals with more nonsynonymous mutations eventually did not survive, resulting in a low nonsynonymous mutation/synonymous in the EMS treatment group.

Previous studies have shown that Stopgain/Stoploss SNPs may alter protein functions (Boschiero et al., 2018). Our results also showed that the numbers of Stopgain/Stoploss SNPs in the individuals treated with EMS were significantly more than those in the control group, indicating that EMS treatment may also affect protein functions by these Stopgain/Stoploss SNPs.

Besides the SNPs in coding regions, the importance of those SNPs located in the non-coding regions cannot be overlooked. In human, almost 90% of disease-related gene variants are found to be located in the non-coding region of the genome. The variations in non-coding regions can cause diseases through alterations in gene expression by changing important regulatory elements (Rojano et al., 2016). In fact, the abundance of SNPs in noncoding regions may serve as an important source of variations for breeding programs (Arabnejad et al., 2018).

Mutation Sites in the EMS Treatment Group

Our results showed that the mutations induced by EMS treatment distributed evenly within the genome, among all the biological processes, cellular components, molecular functions as well as pathways. While EMS seemed to induce mutations randomly in the genome, the mutations in the survivors of the treatment group tend to be less lethal to the animals. More mutations occurred in the intronic or non-coding domains, resulting in fewer non-synonymous mutations that may alter the amino acid sequences of functional proteins. The reasons for this observation may be caused by DNA repair after EMS induction and higher mortality in the individuals with lethal mutations in functional proteins. Recently, Monroe et al. (2022) indicated that mutations tend to occur at sites that are less critical to the survival of plants. The organisms seemed to have developed certain mechanisms to ensure that their important domains not affected by mutations. Our results apparently also support this hypothesis.

Mutations in the Growth-Related Genes

In this study, we found that mutations in MFS and Tubulin were only found in large-sized individuals, suggesting that these mutations in MFS and Tubulin may be the major cause for the rapid growth in these large-sized individual. Similarly, the mutations in Homeobox and Socs were only found in small-sized individuals, suggesting that these mutations may be related to the slow growth in these small-sized individuals.

The major facilitator superfamily (MFS) is a ubiquitous group of proteins involved in the transport of a wide range of compounds (Elena et al., 2006). MFS gene is a meristem regulating gene in rice development and mutation in MFS may result in additional organ growth (Ren et al., 2017). Tubulin has been reported to be involved in cell elongation in flax (Gavazzi et al., 2017). The homeobox family is a large and diverse superclass of genes, many of which may act as transcription factors that play important roles in embryogenesis, tissue differentiation and in animals (Pearson et al., 2005; Di et al., 2015; Hench et al., 2015; Fu et al., 2021) or as master regulators for developmental genes (Ramanathan et al., 2018). Homeobox gene may also be involved in DNA repair during growth (Feltes, 2019). The SOCs are key negative regulators of cytokine and growth factor signaling (Wang et al., 2019) which may exert negative effects on cytokine signaling pathways involved in immunity, growth and development (Liu et al., 2016; Tian et al., 2020). Studies have shown that Socs gene are required for cytokine signal attenuation in mammary epithelial cells and act to limit proliferation through a negative feedback mechanism (Arun et al., 2015). In human cancer cells, activation of SOCs by epigenetic modulation through histone acetylation may induce apoptosis (Sanaei et al., 2020). Despite these lines of indirect evidence, the functions of these growth-related genes have not been well studied. Therefore, further studies are certainly needed to understand how mutations in these genes affect growth in the scallops.

Data Availability Statement

The original contributions presented in the study are publicly available. This data can be found here: (https://submit.ncbi.nlm.nih.gov/subs/bioproject/ BioProject ID: PRJNA810888).

Author Contributions

CHW and BL were responsible for implementation of experiments and wrote the articles. MC and JHN were responsible for providing language help. XL was responsible for writing assistance. CDW designed the experiments and revised the manuscript.

Funding

This work was funded by Focally Deployed Program, Center for Ocean Mega-Science, Chinese Academy of Sciences (NO.COMS2019R03), National Natural Science Foundation (Grant No. 31972791), the Earmarked Fund for Agriculture Seed Improvement Project of Shandong Province (No. 2020LZGC016), the Earmarked Fund for Shandong Modern Agro-Industry Technology Research System (No. SDAIT-14).

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.

References

Alexander, D. H., Novembre, J., and Lange, K. (2009). Fast Model-Based Estimation of Ancestry in Unrelated Individuals. Genome Res. 19 (9), 1655–1664. doi:10.1101/gr.094052.109

PubMed Abstract | CrossRef Full Text | Google Scholar

Arabnejad, M., Dawkins, B. A., Bush, W. S., White, B. C., Harkness, A. R., and McKinney, B. A. (2018). Transition-Transversion Encoding and Genetic Relationship Metric in ReliefF Feature Selection Improves Pathway Enrichment in GWAS. BioData Mining 11, 23. doi:10.1186/s13040-018-0186-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Arun, S. J., Thomson, P. C., Sheehy, P. A., Khatkar, M. S., Raadsma, H. W., and Williamson, P. (2015). Targeted Analysis Reveals an Important Role of JAK-STAT-SOCS Genes for Milk Production Traits in Australian Dairy Cattle. Front. Genet. 6, 342. doi:10.3389/fgene.2015.00342

PubMed Abstract | CrossRef Full Text | Google Scholar

Boschiero, C., Moreira, G. C. M., Gheyas, A. A., Godoy, T. F., Gasparin, G., Mariani, P. D. S. C., et al. (2018). Genome-Wide Characterization of Genetic Variants and Putative Regions under Selection in Meat and Egg-Type Chicken Lines. BMC Genomics 19, 83. doi:10.1186/s12864-018-4444-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, C., Cui, Q. Z., Huang, S. W., Wang, S. H., Liu, X. H., Lu, X. Y., et al. (2018). An EMS Mutant Library for Cucumber. J. Integr. Agric. 17 (07), 1612–1619. doi:10.1016/S2095-3119(17)61765-9

CrossRef Full Text | Google Scholar

Chu, D., and Wei, L. (2019). Nonsynonymous, Synonymous and Nonsense Mutations in Human Cancer-Related Genes Undergo Stronger Purifying Selections Than Expectation. BMC Cancer 19, 359. doi:10.1186/s12885-019-5572-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Danecek, P., Auton, A., Abecasis, G., Albers, C. A., Banks, E., DePristo, M. A., et al. (2011). The Variant Call Format and VCFtools. Bioinformatics 27 (15), 2156–2158. doi:10.1093/bioinformatics/btr330

PubMed Abstract | CrossRef Full Text | Google Scholar

Di, Z., Yu, Y., Wu, Y., Hao, P., He, Y., Zhao, H., et al. (2015). Genome-Wide Analysis of Homeobox Genes from Mesobuthus Martensii Reveals Hox Gene Duplication in Scorpions. Insect Biochem. Mol. Biol. 61, 25–33. doi:10.1016/j.ibmb.2015.04.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Elena, L. E., Teresa, G. J., Carmen, C., and Covadonga, V. (2006). A Novel MFS Transporter Encoding Gene in Fusarium Verticillioides Probably Involved in Iron-Siderophore Transport. Mycol. Res. 110, 1102–1110. doi:10.1016/j.mycres.2006.06.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Feltes, B. C. (2019). Architects Meets Repairers: The Interplay between Homeobox Genes and DNA Repair. DNA Repair 73, 34–48. doi:10.1016/j.dnarep.2018.10.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, S. J., Zhang, J. L., and Xu, H. J. (2021). A Genome-Wide Identification and Analysis of the Homeobox Genes in the Brown Planthopper, Nilaparvata Lugens (Hemiptera: Delphacidae). Arch. Insect Biochem. Physiol. 108 (1), e21833. doi:10.1002/arch.21833

PubMed Abstract | CrossRef Full Text | Google Scholar

Gavazzi, F., Pigna, G., Braglia, L., Gianì, S., Breviario, D., and Morello, L. (2017). Evolutionary Characterization and Transcript Profiling of β-tubulin Genes in Flax (Linum usitatissimum L.) during Plant Development. BMC Plant Biol. 17, 237. doi:10.1186/s12870-017-1186-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Hench, J., Henriksson, J., Abou-Zied, A. M., Lüppert, M., Dethlefsen, J., Mukherjee, K., et al. (2015). The Homeobox Genes of Caenorhabditis E and Insights into Their Spatio-Temporal Expression Dynamics during Embryogenesis. PLoS One 10, e0126947. doi:10.1371/journal.pone.0126947

PubMed Abstract | CrossRef Full Text | Google Scholar

Hou, J. L., Zhang, X. Y., Wang, G. X., Sun, Z. H., Du, W., Zhao, Y. X., et al. (2019). Novel Breeding Approach for Japanese Flounder Using Atmosphere and Room Temperature Plasma Mutagenesis Tool. BMC Genomics 20, 323. doi:10.1186/s12864-019-5681-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Johanna, L., Shahriar, S. S., Sameer, H., Nick, S., Mats, T., Olof, O., et al. (2020). Development and Characterization of an EMS-Mutagenized Wheat Population and Identification of Salt-Tolerant Wheat Lines. BMC Plant Biol. 20, 18. doi:10.1186/s12870-019-2137-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Kuroyanagi, M., Katayama, T., Imai, T., Yamamoto, Y., Chisada, S.-i., Yoshiura, Y., et al. (2013). New Approach for Fish Breeding by Chemical Mutagenesis: Establishment of TILLING Method in Fugu (Takifugu Rubripes) with ENU Mutagenesis. BMC Genomics 14, 786. doi:10.1186/1471-2164-14-786

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, C.-H., Rimesso, G., Reynolds, D. M., Cai, J., and Baker, N. E. (2016). Whole-Genome Sequencing and iPLEX MassARRAY Genotyping Map an EMS-Induced Mutation Affecting Cell Competition in Drosophila Melanogaster. G3 Genes, Genomes, Genet. 6, 3207–3217. doi:10.1534/g3.116.029421

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., and Durbin, R. (2009). Fast and Accurate Short Read Alignment with Burrows-Wheeler Transform. Bioinformatics 25, 1754–1760. doi:10.1093/bioinformatics/btp324

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., et al. (2009). The Sequence Alignment/Map Format and SAMtools. Bioinformatics 25, 2078–2079. doi:10.1093/bioinformatics/btp352

PubMed Abstract | CrossRef Full Text | Google Scholar

Lian, X., Liu, Y., Guo, H., Fan, Y., Wu, J., Guo, H., et al. (2020). Ethyl Methanesulfonate Mutant Library Construction in Gossypium Hirsutum L. For Allotetraploid Functional Genomics and Germplasm Innovation. Plant J. 103, 858–868. doi:10.1111/tpj.14755

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C.-Z., He, A.-Y., Chen, L.-Q., Limbu, S. M., Wang, Y.-W., Zhang, M.-L., et al. (2016). Molecular Characterization and Immune Response to Lipopolysaccharide (LPS) of the Suppressor of Cytokine Signaling (SOCS)-1, 2 and 3 Genes in Nile tilapia (Oreochromis Niloticus). Fish Shellfish Immunol. 50, 160–167. doi:10.1016/j.fsi.2016.01.027

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, X., Li, C., Chen, M., Liu, B., Yan, X. J., Ning, J. H., et al. (2020). Draft Genomes of Two Atlantic Bay Scallop Subspecies Argopecten Irradians Irradians and A. i. Concentricus. Sci. Data 7 (1), 99. doi:10.1038/s41597-020-0441-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, Y., Dai, S., Gu, A., Liu, M., Wang, Y., Luo, S., et al. (2016). Microspore Induced Doubled Haploids Production from Ethyl Methanesulfonate (EMS) Soaked Flower Buds Is an Efficient Strategy for Mutagenesis in Chinese Cabbage. Front. Plant Sci. 7, 1780. doi:10.3389/fpls.2016.01780

PubMed Abstract | CrossRef Full Text | Google Scholar

Monroe, J. G., Srikant, T., Carbonell-Bejerano, P., Becker, C., Lensink, C., Exposito-Alonso, M., et al. (2022). Mutation Bias Reflects Natural Selection in Arabidopsis T. Nature 602, 101–105. doi:10.1038/s41586-021-04269-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Ning, X., Feng, L. Y., Li, X., Wang, S., Zhang, M., Wang, S., et al. (2018). The Scallop IGF2 mRNA-Binding Protein Gene PyIMP and Association of a Synonymous Mutation with Growth Traits. Genes Genet. Syst. 93, 91–100. doi:10.1266/ggs.17-00028

PubMed Abstract | CrossRef Full Text | Google Scholar

Ornuma, T., Wilawan, N., Patiruj, J., and Sontichai, C. (2018). EMS-Induced Mutation Followed by Quizalofop-Screening Increased Lipid Productivity in Chlorella Sp. Bioproc. Biosyst. Eng. 41, 613–619. doi:10.1007/s00449-018-1896-1

CrossRef Full Text | Google Scholar

Pearson, J. C., Lemons, D., and McGinnis, W. (2005). Modulating Hox Gene Functions during Animal Body Patterning. Nat. Rev. Genet. 6, 893–904. doi:10.1038/nrg1726

PubMed Abstract | CrossRef Full Text | Google Scholar

Penna, S., Vitthal, S. B., and Yadav, P. V. (2012). In Vitro Mutagenesis and Selection in Plant Tissue Cultures and Their Prospects for Crop Improvement. Bioremediation, Biodiversity, Bioavailability 6, 6–14. doi:10.1371/journal.pgen.1008243

CrossRef Full Text | Google Scholar

Ramanathan, A., Srijaya, T. C., Sukumaran, P., ZainRosnah, R. B., and Abu Kasim, N. H. (2018). Homeobox Genes and Tooth Development: Understanding the Biological Pathways and Applications in Regenerative Dental Science. Arch. Oral Biol. 85, 23–39. doi:10.1016/j.archoralbio.2017.09.033

PubMed Abstract | CrossRef Full Text | Google Scholar

Ren, D. Y., Li, Y. F., Wang, Z., Xu, F. F., Guo, S., Zhao, F. M., et al. (2017). Identification and Gene Mapping of a Multi-Floret Spikelet 1 (Mfs1) Mutant Associated with Spikelet Development in rice. J. Integr. Agric. 11, 1574–1579. doi:10.1016/S2095-3119(12)60160-9

CrossRef Full Text | Google Scholar

Rojano, E., Ranea, J. A., and Perkins, J. R. (2016). Characterisation of Non-Coding Genetic Variation in Histamine Receptors Using AnNCR-SNP. Amino Acids 48, 2433–2442. doi:10.1007/s00726-016-2265-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Sanaei, M., Kavoosi, F., and Behjoo, H. (2020). Effect of Valproic Acid and Zebularine on SOCS-1 and SOCS-3 Gene Expression in colon Carcinoma SW48 Cell Line. Exp. Oncol. 42, 183–187. doi:10.32471/exp-oncology.2312-8852.vol-42-no-3.15113

PubMed Abstract | CrossRef Full Text | Google Scholar

Schiffels, S., and Durbin, R. (2014). Inferring Human Population Size and Separation History from Multiple Genome Sequences. Nat. Genet. 46 (8), 919–925. doi:10.1038/ng.3015

PubMed Abstract | CrossRef Full Text | Google Scholar

Sikora, P., Chawade, A., Larsson, M., Olsson, J., and Olsson, O. (2011). Mutagenesis as a Tool in Plant Genetics, Functional Genomics, and Breeding. Int. J. Plant Genomics 2011, 314829. doi:10.1155/2011/314829

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, S., Guo, Y., Zhang, X., Zhang, J., and Ma, E. (2014). Screening and SNP Mapping of Copper-Resistant Mutations in C. Elegans. Yi Chuan 36 (12), 1261–1268. doi:10.3724/SP.J.1005.2014.1261

PubMed Abstract | CrossRef Full Text | Google Scholar

Tian, B., Tang, D., Wu, J., Liang, M., Hao, D., and wei, Q. (2020). Molecular Characterization, Expression Pattern and Evolution of Nine Suppressors of Cytokine Signaling (SOCS) Gene in the Swamp Eel (Monopterus Albus). Fish Shellfish Immunol. 96, 177–189. doi:10.1016/j.fsi.2019.11.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Till, B. J., Cooper, J., Tai, T. H., Colowit, P., Greene, E. A., Henikoff, S., et al. (2007). Discovery of Chemically Induced Mutations in Rice by TILLING. BMC Plant Biol. 7, 19. doi:10.1186/1471-2229-7-19

PubMed Abstract | CrossRef Full Text | Google Scholar

Tran, Q. H., Bui, N. H., Kappel, C., Dau, N. T. N., Nguyen, L. T., Tran, T. T., et al. (2020). Mapping-by-Sequencing via MutMap Identifies a Mutation in ZmCLE7 Underlying Fasciation in a Newly Developed EMS Mutant Population in an Elite Tropical Maize Inbred. Genes 11 (3), 281. doi:10.3390/genes11030281

PubMed Abstract | CrossRef Full Text | Google Scholar

Vilella, A. J., Severin, J., Ureta-Vidal, A., Heng, L., Durbin, R., and Birney, E. (2009). EnsemblCompara GeneTrees: Complete, Duplication-Aware Phylogenetic Trees in Vertebrates. Genome Res. 19, 327–335. doi:10.1101/gr.073585.107

PubMed Abstract | CrossRef Full Text | Google Scholar

Wakeley, J. (1996). The Excess of Transitions Among Nucleotide Substitutions: New Methods of Estimating Transition Bias Underscore its Significance. Trends Ecol. Evol. 11, 158–162. doi:10.1016/0169-5347(96)10009-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, B., Wangkahart, E., Secombes, C. J., and Wang, T. (2019). Insights into the Evolution of the Suppressors of Cytokine Signaling (SOCs) Gene Family in Vertebrates. Mol. Biol. Evol. 36, 393–411. doi:10.1093/molbev/msy230

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, C. D., Liu, B., Liu, X., Ma, B., Zhao, Y. M., Zhao, X., et al. (2017). Selection of a New Scallop Strain, the Bohai Red, from the Hybrid between the bay Scallop and the Peruvian Scallop. Aquaculture 479, 250–255. doi:10.1016/j.aquaculture.2017.05.045

CrossRef Full Text | Google Scholar

Wang, K., Li, M., and Hakonarson, H. (2010). ANNOVAR: Functional Annotation of Genetic Variants from High-Throughput Sequencing Data. Nucleic Acids Res. 38, e164. doi:10.1093/nar/gkq603

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, Y., Sun, G., Zeng, Q., Chen, Z., Hu, X., Li, H., et al. (2018). Predicting Growth Traits with Genomic Selection Methods in Zhikong Scallop (Chlamys Farreri). Mar. Biotechnol. 20, 769–779. doi:10.1007/s10126-018-9847-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Xiong, H., Zhou, C., Guo, H., Xie, Y., Zhao, L., Gu, J., et al. (2020). Transcriptome Sequencing Reveals Hotspot Mutation Regions and Dwarfing Mechanisms in Wheat Mutants Induced by γ-ray Irradiation and EMS. J. Radiat. Res. 61 (1), 44–57. doi:10.1093/jrr/rrz075

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J., Lee, S. H., Goddard, M. E., and Visscher, P. M. (2011). GCTA: A Tool for Genome-Wide Complex Trait Analysis. Am. J. Hum. Genet. 88, 76–82. doi:10.1016/j.ajhg.2010.11.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoshihito, T., Shunichi, T., Makoto, F. S., Yasuhiro, K., Takeshi, T., Takao, S., et al. (2006). Generation of Medaka Gene Knockout Models by Target-Selected Mutagenesis. Genome Biol. 7, 53–64. doi:10.1186/gb-2006-7-12-r116

CrossRef Full Text | Google Scholar

Zhang, C., Dong, S.-S., Xu, J.-Y., He, W.-M., and Yang, T.-L. (2018). PopLDdecay: a Fast and Effective Tool for Linkage Disequilibrium Decay Analysis Based on Variant Call Format Files. Bioinformatics 35, 1786–1788. doi:10.1093/bioinformatics/bty875

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, F. S., He, Y. C., Liu, X. S., and Ma, J. G. (1986). On the Introduction of Spat-Rearing and Experimental Culture of bay Scallop, Argopecten Irradians (Lamarck). Aquaculture 57, 1–4. doi:10.1016/0044-8486(86)90244-9

CrossRef Full Text | Google Scholar

Zhao, L., Li, Y., Lou, J., Yang, Z., Liao, H., Fu, Q., et al. (2019). Transcriptomic Profiling Provides Insights into Inbreeding Depression in Yesso Scallop Patinopecten Yessoensis. Mar. Biotechnol. 21, 623–633. doi:10.1007/s10126-019-09907-9

PubMed Abstract | CrossRef Full Text | Google Scholar