Xinghai Zhu1†
Pingping Liu1†
Xiujiang Hou1Junhao Zhang1Jia Lv1Wei Lu1Qifan Zeng1
Xiaoting Huang1,2
Qiang Xing1,2*
Zhenmin Bao1,2- 1MOE Key Laboratory of Marine Genetics and Breeding, College of Marine Life Sciences, Ocean University of China, Qingdao, China
- 2Laboratory for Marine Fisheries Science and Food Production Processes, Qingdao National Laboratory for Marine Science and Technology, Qingdao, China
The increasing sea temperature caused by global warming has resulted in severe mortalities in maricultural scallops. Therefore, improving thermal tolerance has become an active research area in the scallop farming industry. Bay scallop (Argopecten irradians irradians) was introduced into China in 1982 and has developed into a vast aquaculture industry in northern China. To date, genetic studies on thermal tolerance in bay scallops are limited, and no systematic screening of thermal tolerance-related loci or genes has been conducted in this species. In the present study, we conducted a genome-wide association study (GWAS) for thermal tolerance using the Arrhenius break temperature (ABT) indicators of 435 bay scallops and 38,011 single nucleotide polymorphism (SNP) markers. The GWAS identified 1,906 significant thermal tolerance-associated SNPs located in 16 chromosomes of bay scallop. Gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses showed that 638 genes were enriched in 42 GO terms, while 549 annotated genes were enriched in aggregation pathways. Additionally, the SNP (15-5091-20379557-1) with the lowest P value was located in the transcriptional coactivator p15 (PC4) gene, which is involved in regulating DNA damage repair and stabilizing genome functions. Further analysis in another population identified two new thermal tolerance-associated SNPs in the first coding sequence of PC4 in bay scallops (AiPC4). Moreover, AiPC4 expression levels were significantly correlated (r = 0.675–0.962; P < 0.05) with the ABT values of the examined bay scallops. Our data suggest that AiPC4 might be a positive regulator of thermal tolerance and a potential candidate gene for molecular breeding in bay scallop aiming at thermal tolerance improvement.
Introduction
Aquatic products constitute an important source of food consumed all over the world. As of 2019, global aquaculture output was 80.1 million tons (46.8 million tons contributed by China), generating an estimated 238.0 billion dollars (139.1 billion dollars by China) in revenue (FAO, 2019). Many bivalves are economically important aquaculture species that account for approximately 22.8% of the aquatic products market worldwide (Moffitt and Cajas-Cano, 2014). Although a continuous increase in the production of bivalves has been observed in recent years (Steeves et al., 2018), several obstacles (e.g., environmental stress, biotic and abiotic factors) hamper the development of the industry. Many studies revealed that, anoxia tolerance (Zwaan et al., 1995; Chen et al., 2007a), salinity tolerance (Chen and Chen, 2000; Roldan-Carrillo et al., 2005), and tolerance to bacterial pathogens (Richards et al., 2015; Rojas et al., 2016) are closely related to mortality in bivalves, which may have adverse effect on global aquaculture production. Consequently, researchers gave more attention on the tolerance ability of aquatic commercial bivalves for sustainable development of shellfish industry.
Bivalves are poikilotherms, therefore, temperature is one of the most important environmental factors affecting their physiology (Zittier et al., 2015), including growth efficiency (Sicard et al., 1999), reproductive activity (Martínez and Pérez, 2003), immune defense (Chen et al., 2007b), tolerance ability (Chen and Chen, 2000), and even leading to death (Cheng et al., 2020) in bivalves. For instance, the daily shell growth rate of Antarctic scallops (Adamussium colbecki) during summer was severely affected by temperature (Heilmayer et al., 2005). During the reproductive cycle of mussel (Mytilus edulis), variations in temperature significantly affected the development of the gonads (Lubet, 1983). An increase in temperature significantly decreased the activities of its hemocytes, and consequently increased the mortality rate of Pacific oysters (Gagnaire et al., 2006). Most importantly, thermal tolerance in scallops has been a major target trait in breeding because massive mortality usually occurs in summer. High temperature of seawater in summer even caused a 37% decline in the production of Zhikong scallop (MAC, 1999). Consequently, accurate and efficient indicators for thermal tolerance evaluation have been playing a crucial role in accelerating the breeding progress of thermal resistant scallops (Mason, 1958; Granados-Amores et al., 2012; Yang et al., 2013; Xing et al., 2016a; Ning et al., 2019). Specifically, the heart rate (HR)-based Arrhenius break temperature (ABT), the temperature at which the Arrhenius plot exhibits a sharp discontinuity in slope, is regarded as a heterothermic organism’s upper limit temperature (Dahlhoff et al., 1991) and could represent species’ maximal habitat temperature (Weinstein and Somero, 1998). ABT was originally reported in porcelain crabs (Bamber and Depledge, 1997), followed by intertidal animals (Somero, 2002), snails (Stenseng et al., 2005), and fish (Casselman et al., 2012). Xing et al. (2016a) first reported the feasibility of ABT as a rapid, efficient, and non-invasive indicator for bivalve thermal tolerance evaluation (Xing et al., 2016a). Specifically, ABT could accurately indicate thermal limits in interspecific scallops with different thermal adaptation, with different ABT values of 22.03 ± 0.19, 29.10 ± 0.25, 32.20 ± 0.25 and 34.09 ± 0.19°C for Patinopecten yessoensis, Chlamys farreri, Argopecten irradians, and A. ventricosus, respectively. In intraspecific Zhikong scallops, ABT negatively correlated with sizes, weights and ages of individuals, apparently differed between sexual adults, and was significantly affected by spawning behavior and seasonal variations of seawater temperature (Xing et al., 2016a, 2021). As a consequence, ABT has been used as an efficient indicator in predicting the thermal performance of shellfish (Chen et al., 2016; Xing et al., 2019, 2021), and can also be used as an indicator in the breeding of new variants of thermal resistant aquatic species.
Owing to the fast growth rate and high production of the bay scallop (A. irradians irradians), it has become one of the most commercially important shellfish in China since its introduction from the US in the early 1980s (Wang et al., 2007). However, in the past decades, the scallop production industry has suffered high summer mortalities, and the high temperature has been suspected to be one of the main environmental inducers (Xiao et al., 2005; Jonasson et al., 2007). As reported, compared with the heat-resistant scallop A. irradians concentricus distributed in southern China, scallop A. irradians irradians distributed in northern China exhibited a lower thermal tolerance (Liu et al., 2020). Using the shell height as the target trait, genome-wide genotyping method 2b-RAD (Wang et al., 2012), RADtyping software (Fu et al., 2013) and genomic selection (GS) strategies (Sonesson and Meuwissen, 2009; Wang et al., 2018b), researchers have been able to develop the new bay scallop strain named “Haiyifeng 11.” Interestingly, the non-invasive thermal tolerance indicator ABT (Xing et al., 2016a), regarding as thermal limit of the new strain was 34.54 ± 0.53°C, which is approximately 2°C higher than that of the common population. Accordingly, the survival rate of the thermal tolerant strain (81.32%) increased by 13.22% compared to the control group (68.10%). The thermal tolerance of “Haiyifeng 11” improved markedly after four consecutive generations of breeding via GS. However, the major loci or genes and molecular mechanisms underlying thermal tolerance in bay scallop are yet to be identified. In practical application, knowledge of the loci or genes regulating thermal tolerance could be used as potential markers in target genotyping technology, e.g., high-throughput, sequencing-based GoldenGate approach (HD-marker) for scallop selection breeding with better thermal tolerance performance (Lv et al., 2018).
With the advancement in technology, genome-wide association studies (GWAS) have been successfully utilized in identifying loci/genes involved in target polygenic traits; GWAS is a valuable tool for genetic improvement of target phenotypic traits (Korte and Farlow, 2013). In aquatic species, GWAS has been widely used in investigating economically phenotypic traits, especially in economic fishes, including head size (Geng et al., 2016), disease resistance (Zhou et al., 2017), heat endurance (Jin et al., 2017) and hypoxia tolerance (Wang et al., 2017). Besides, GWAS revealed that two SNPs positioned in Salmo salar chromosome 26 showed strong associations to its flesh color (Helgeland et al., 2019), and five SNPs distributed in four chromosomes were strongly associated with acute heat tolerance trait in large yellow croaker (Wu et al., 2021). Except for fish, GWAS has also been proven to be an efficient approach to reveal trait-associated genes and/or markers in shellfish. In Zhikong scallop (C. farreri), GWAS identified two growth-related quantitative trait loci and one potential sex-determination region (Jiao et al., 2014). In Yesso scallop (P. yessoensis), GWAS identified PyBCO-like 1 (beta-carotene oxygenase like 1) associated with carotenoid coloration (Li et al., 2019a), E2F3 (E2F transcription factor 3) which may regulate its growth (Ning et al., 2019), and three candidate genes involved in carotenoid metabolism, which are necessary for the formation of reddish-orange shells (Zhao et al., 2017). Additionally, GWAS analysis of heat tolerance in Pacific abalone revealed that 27 SNPs significantly affected its adhesion ability (Yu et al., 2021). These results provide a list of candidate loci/genes for functional annotation, decipher the genetic architecture of complex traits and help in identifying genetic variants that can be directly used for genetic improvement.
In the present study, we aimed to identify significant SNP markers associated with thermal tolerance in bay scallops by performing GWAS, with the aid of the non-invasive thermal tolerance indicator ABT. GWAS showed that the highest potential locus was located in the gene encoding transcriptional coactivator p15 (PC4), which repairs DNA damages and stabilizes genome function (Das et al., 2006). Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were further carried out. Subsequent sequencing analysis was conducted in another population to verify the SNPs. We further investigated the expression patterns of the PC4 gene (AiPC4) in the tissues of bay scallop and detected a correlation between the expression levels of AiPC4 and ABT values in bay scallops. Our results revealed the importance of AiPC4 in thermal tolerance in bay scallops, and provided valuable loci and candidate genes for the breeding of thermal tolerant scallops.
Materials and Methods
Scallop Management and Miso-RAD Library Construction
“Haiyifeng 11” new strain was artificially bred from 2011 to 2016 by our group and was used for GWAS in the present study. The detailed procedure of artificial selection for “Haiyifeng 11” could be obtained in Supplementary File 1. Briefly, the new strain was initially bred from the founder population of 1,000 bay scallops, which consisted of 500 black-brown shell color individuals from Laizhou, Yantai (37°10′57 N, 119°56′55 E, Shandong Province, China) and another 500 purple-red shell color individuals from Jiaonan, Qingdao (35°52′20 N, 120°2′47 E, Shandong Province, China) in 2011. Subsequently, the new strain was bred via four consecutive generations by GS focusing on traits of shell height from 2012 to 2016.
Healthy 10-month-old “Haiyifeng11” bay scallops (N > 1,000) were contemporaneously breeding populations and were collected from artificial scallop-rearing substrates at Qingdao Jinshatan Fishery Group Co., Qingdao (35°57′23 N, 120°27′28 E, Shandong Province, China) in 2016. After sample collection, encrusted organisms on scallop shells were removed and the scallops were brought to the laboratory following standard procedures (Maeda-Martínez et al., 2000). The scallops were placed in plastic tanks containing filtered and aerated seawater (consistent with natural seawater temperature) for 7 days for acclimation (Bakhmet, 2005). To avoid tank effects, all collected scallops were randomly maintained in several plastic tanks and the water was partially replaced (1/2) daily. The animals were maintained without feed to avoid specific dynamic actions (Secor, 2009). Subsequently, we randomly selected 435 scallops and determined their ABT values using the method described by Xing et al. (2016a). Briefly, the heart rate (HR) of 435 individuals were measured to calculate ABT values. To detect the HR, the non-invasive method was employed, in which CNY-70 optical sensors were fixed with glue (Krazy glue, Westerville, United States) to the shell surface close to cardiocoelom of the scallops. Variations in the infrared signal were amplified, filtered and recorded using AMP03 (Newshift, Lisbon, Portugal) and Powerlab (8/35, ADInstruments, Sydney, Australia). Data were viewed and analyzed using LabChart software (ADInstruments, Sydney, Australia). ABT was calculated using regression analyses to generate the best fit lines on both sides of a putative break point (Dahlhoff et al., 1991; Stillman and Somero, 1996). Then, gills were removed from each scallop following the procedure described by Mao et al. (2013) and then preserved in 100% ethanol at −20°C for genomic analysis. Our experiments were conducted according to the guidelines and regulations established by the Ocean University of China and the local government.
Genomic DNA was extracted from the gills using a standard phenol-trichloromethane method with some modification (Sambrook et al., 1989). Subsequently, libraries were constructed using the Multi-isoRAD method developed by Wang et al. (2016). Briefly, genomic DNA of the 435 selected bay scallops was digested with BsaXI (New England BioLabs, cat. no. R0609) at 37°C for 45 min. Adaptor ligation, PCR amplification, digestion, and ligation were carried out, after which five DNA fragments were set as one group (barcoding and pooling). After amplifying twice, each production of five concatenated tags from five samples was purified using the MinElute PCR Purification Kit (Qiagen, cat. no. 28004) and pooled for pair-end sequencing on the Illumina HiSeq2000 platform.
Sequence Data Processing and Genotyping
To obtain high-quality reads, the raw reads were first preprocessed to remove paired-end (PE) reads with ambiguous base calls (N), long homopolymer regions (>10 bp), and low-quality bases (> 20% of bases with quality score < 10). After matching the forward and reverse reads, PEAR software was used to assemble the PE reads (Zhang et al., 2014). The reads containing five tags, which were chosen by the assembly step, were divided into single-tag datasets. The single tag containing the restriction enzyme site, extracted from single-tag data, was aligned to the chromosome-level reference genome constructed by our group (unpublished) using SOAP (Li et al., 2009). Genotyping analysis of sequencing data was performed using the RADtyping software developed by Fu et al. (2013). To generate genotypes sets for GWAS, only SNPs without missing values (SNPs with coverage < 4 were considered missing), call rate > 60%, and minor allele frequency (MAF) > 0.01 were used for downstream analysis.
Genome-Wide Association Study and Identification of Candidate Genes
Association analysis between genotypes and ABT values was conducted using Plink-1.0.7 software (Purcell et al., 2007) by performing a Wald test using a general linear model (GLM). The Wald statistic P-value for each SNP was summarized in ascending order. The sequence tag of SNPs with P < 0.001 was extracted for BLASTN (Altschul et al., 1997) search against the draft genome of A. irradians irradians (unpublished data) that were annotated by the protein database NR using the code: blastall −p blastn −i input −d genome −o result −a 5 −e 1e-4 −m 8. Additionally, Manhattan plot was generated using “qqman” function on R package (version 3.6.3). Potential SNPs with P < 0.05 were captured for bioinformatics analysis, including GO function annotation and KEGG pathway enrichment analyses1.
SNP Verification
To verify the association between the ABT values and the gene harboring the most significant SNP, another separate cultured population established by artificial fertilization with over 1,000 sexually mature bay scallops at Marine Seed Co., Ltd., Yenta (37°18′18 N, 119°52′29 E, Shandong Province, China) in 2018 was used for screening thermal tolerance associated mutations in AiPC4. From 150 randomly collected individuals in this population, the top 15 (10%) and the bottom 15 (10%) scallops were selected based on their ABT values in decreasing order. Then, the genomic DNA of the 30 scallops was extracted, and the three-prime UTR (chr15-5091-20379557-1), coding sequence (CDS) 1 and 2 (AiPC contains two CDS) of AiPC4 were amplified using the primer pairs listed in Table 1. The PCR reaction of the above three target DNA fragments was carried out as follows: 2 μL of template DNA (50 ng/μL), 15 μL of Q5 High-Fidelity 2 × Master Mix (New England BioLabs, cat. no. M0492S), 0.4 μL (10 pmol/μL) of each forward and reverse primer, and sterile water up to 30 μL volume was added for amplifying. The related parameters in the PCR reaction condition was as follows: pre-denaturing at 98°C for 30 s, followed by denaturing at 98°C for 10 s, annealing at 60/56/56°C for 15/15/20 s, and extension at 72°C for 10/10/20 s, respectively. After 25 cycles, the sample was extended to 72°C for 2 min. The PCR products were sequenced using the Sanger method by Sangon Biotech (Shanghai, China). Sequences amplified by the same primer pair were compared among the 30 scallops using the ClustalW2 multiple alignment program2. When a mutation was detected, comparisons of genotype frequencies between the top 15 and the bottom 15 scallops were performed using Fisher’s exact test to identify the thermal tolerance-associated mutation. P values < 0.05 were considered statistically significant.
Table 1. Sequences of primers used in this study.
Candidate Gene Evaluation
To investigate the expression pattern of AiPC4 and further analyze the correlation between expression levels of AiPC4 and thermal tolerance in bay scallops, 18 randomly selected bay scallops from another separate population (N > 500) at Marine Seed Co., Ltd., Yenta (37°18′18 N, 119°52′29 E, Shandong Province, China) in 2020 were recruited and their tissues (mantle, gill, gonad, kidney, hepatopancreas, striated muscle, and heart) were sampled. Total RNA of the sampled tissues was isolated following the method described by Hu et al. (2006) and then digested with DNase I (TaKaRa, Shiga, Japan). The concentration and purity of the RNA were determined using a Nanovue Plus spectrophotometer (GE Healthcare, NJ, United States), and RNA integrity was assessed by agarose gel electrophoresis. First-strand cDNA was synthesized according to the manufacturer’s protocol using Moloney murine leukemia virus (MMLV) reverse transcriptase (Thermo, United States) in a 20 μL volume with 2 μg of each total RNA sample as the template and 0.5 μg of oligo (dT)18 (TaKaRa Biotechnology, Liaoning, China) as the primer. The mixture was denatured at 65°C for 5 min and then chilled immediately on ice. After adding the reverse transcriptase, reaction buffer, and dNTPs, cDNA was amplified under the following conditions: 42°C for 90 min and 72°C for 10 min. The cDNA was stored at −20°C and diluted to 5 ng/μL for use as the template in quantitative reverse-transcription PCR (qRT-PCR).
The expression levels of AiPC4 in the tissues of the 18 bay scallops were analyzed by qRT-PCR. The first strand cDNA from the tissues was used as the template and specific primers were designed to amplify the AiPC4 fragments. Elongation factor 1 alpha (EF1α) was selected as a reference for tissue samples (Li et al., 2019b). The primers used for the qRT-PCR are listed in Table 1. Three technical repeats were performed for each reaction. All qRT-PCR reactions were performed using a LightCycler 480 system (Roche Diagnostics, Mannheim, Germany) with a total volume of 20 μL containing 1 × Real-time PCR Master Mix containing SYBR Green dye (TOYOBO, Osaka, Japan), 0.2 μM of each primer, and 2 μL cDNA template. The qRT-PCR reaction was carried out as follows: denaturation at 50°C for 2 min, 94°C for 10 min, followed by 40 cycles of 94°C for 15 s and 62°C for 1 min. Dissociation (from 94°C to 62°C) analysis of the amplification products was performed at the end of each PCR reaction to confirm that only one PCR product was amplified and detected.
Data from the real-time PCR were obtained, and the gene EF1A (elongation factor 1α) was used as an internal control. Fold-change of AiPC4 expression levels among tissues was conducted with the gene expression level in gill as the reference. Statistical analysis of the data was performed with SPSS software 21.0 (IBM Co., Armonk, United States) via one-way ANOVA followed by a post hoc test. We used bivariate correlation analysis to determine the relationships between AiPC4 expression levels in scallop tissues and the ABT values. Differences were considered significant if P < 0.05.
Results
ABT Values and Sequencing Data
The ABT values of the 435 bay scallops varied from 33.69 to 35.73°C, with a mean value of 34.65 ± 0.40°C. The ABT values basically follow the law of normal distribution, which is necessary for the subsequent GWAS analysis for quantitative traits (Supplementary Figure 1).
In total, 350.54 Gigabyte (G) clean data of 87 Multi-isoRAD libraries generated 1,949,337,356 reads after sequencing (GenBank: PRJNA689862). The average sequencing depth was 17.2×. After quality control, we obtained 1,908,407,252 high-quality reads, and an average of 45.15% high-quality reads of each individual were uniquely mapped to the bay scallop genome (unpublished data) by SOAP. We genotyped 489,750 SNPs in 435 individuals using RADtyping. Specifically, population genetic analysis via PCA analysis showed that the 435 individuals have no population stratification which is critical for subsequent GWAS.
Genomic Regions Associated With ABT
A total of 489,750 SNP markers obtained from 435 bay scallops were subjected to quality control using PLINK software. After screening for call rate (> 60%) and MAF (> 1%), a set of 38,079 high-quality SNPs were used for GWAS. We successfully mapped 38,011 of these SNPs in the 16 chromosomes of bay scallop. The results of the GWAS showed that 1,906 SNPs were significantly (P < 0.05) correlated with the ABT and 1,165 of these SNPs in the 27 bp tags were directly annotated in bay scallop genome. The top 14 significant (P < 0.001) SNPs are shown in Table 2. Among the top 14 SNPs, the most significantly associated SNP (chr15-5091-20379557-1), which was located on chromosome 15 (Figure 1), was selected for verification by sequencing the PCR product using Sanger method (Sangon Biotech, Shanghai, China).
Table 2. Information of the top fourteen SNPs associated with ABT.
Figure 1. Genome-wide association study (GWAS) of thermal tolerance traits with 38,079 SNPs in bay scallop. The red line indicates the genome-wide threshold as P-value of 0.001.
SNP Verification
The most significantly associated SNP, chr15-5091-20379557-1, was assigned to a gene encoding PC4, which is involved in repairing DNA damages and stabilizing genome function (Das et al., 2006), which may play an important role in thermal tolerance. And this SNP on AiPC4 gene explained 5.41% variance of the total phenotypic variation according to the GWAS result. To determine the relationship between AiPC4 and thermal tolerance in scallops, mutations spanning the transcribed sequence of AiPC4 were screened in the top 15 (ABT = 33.19 ± 1.57°C) and bottom 15 (ABT = 30.47 ± 1.75°C) individuals (based on ABT values) from 150 randomly collected scallops (ABT values ranged from 28.32 to 35.22°C, with a mean value of 31.43 ± 1.67°C). We identified two SNPs, including c.234T < C and c.369A < G in the CDS 1 of AiPC4 (Figure 2). Comparison between the top and bottom groups showed that the genotype frequencies at the two SNPs were significantly different (P = 0.008 and 0.015, respectively) (Table 3). These results support the findings of GWAS that the PC4 locus might be involved in regulating thermal tolerance in bay scallops.
Figure 2. A few neighboring SNP loci on chromosome 15 are independently associated with ABT in bay scallop (A). The genic structure of AiPC4 (B). The purple box indicate the 5′UTRs and 3′UTRs (untranslated region); the white boxes represent CDS (coding sequence); and the red arrows indicate the two verified SNPs. The CDS are shown relative to their lengths and the two indicated SNPs are marked in the accurate positions in the genic sequences.
Table 3. Comparison of genotype frequencies of SNPs in AiPC4 CDS 1 between the top (N = 15) and the bottom (N = 15) scallop groups classified by the ABT values of thermal tolerance.
Identification of Candidate Genes
Among the top 14 SNPs, nine were identified in the following gene regions; PC4, MYHC (myosin heavy chains), CTLM (ste18p Gγ subunit), ALDH (aldehyde dehydrogenases), PDE11 (phosphodiesterase 11), BCL7 (B cell leukemia/lymphoma 7), SCAD (short-chain acyl-CoA dehydrogenase), RALDH1 (retinaldehyde dehydrogenase family 1), and one uncharacterized gene (LOC110459767). At a less stringent significance level (P < 0.05), there were 1,906 SNPs annotated in 950 genes. To investigate the function of these genes, we conducted GO function annotation and KEGG pathway enrichment analyses. The results showed that 638 genes were enriched in 42 GO terms (level 2) (Figure 3), and the dominant terms were cellular process (293 genes) in biological process, membrane (152 genes) in cellular component, and binding (434 genes) in molecular function. The results of KEGG pathway enrichment analysis showed that 549 genes were enriched in metabolism, genetic information processing, environmental information processing, organismal systems, and human disease (Figure 4). It should be noted that both the KEGG pathway annotation and the GO terms showed a similar trend, with environmental information processing to be referred as the most significantly enriched, including “signal transduction,” “signal molecular and interaction,” and “membrane transport.”
Figure 3. The summary of GO function annotation analysis for candidate selective genes that comprised significant SNPs (P < 0.05). Each term described the function of gene cluster and the length of colored bars represent the numbers difference of genes.
Figure 4. The summary of KEGG pathway enrichment analysis performed for biological process, cellular component and molecular function. Functional annotation of candidate selective genes that comprised significant SNPs (P < 0.05).
AiPC4 Expression Patterns and Its Involvement in Thermal Tolerance
We further investigated the expression patterns of AiPC4 in the tissues of adult bay scallops (mantle, gill, kidney, testis, ovary, striated muscle, hepatopancreas, and heart) using qRT-PCR data generated from 18 individuals. As shown in Figure 5, AiPC4 gene was ubiquitously expressed in all the examined tissues; however, the level of expression was predominantly detected in the mantle, ovary, striated muscle, and heart. The expression of AiPC4 was highest in the heart (9.36-fold compared with the expression in gill), followed by the ovary (4.59-fold), striated muscle (2.55-fold), mantle (2.27-fold), hepatopancreas (1.86-fold), testis (1.61-fold), kidney (1.57-fold), and gill (1.00-fold). The expression levels of AiPC4 was significantly (P < 0.05) higher in the heart and ovary than in the other tissues.
Figure 5. Relative expression levels of AiPC4 in tissues of bay scallop. Three replicates were performed for each adult tissue and three technical replicates were conducted for each PCR. The comparison of the expression levels of AiPC4 among different tissues as performed using one way ANOVA followed by with a post hoc test. Bars with different superscripts indicate significant differences (P < 0.05).
To evaluate the possible role of AiPC4 in thermal tolerance, we determined the correlation between the ABT values of 18 bay scallops and the AiPC4 expression levels in their tissues. As shown in Figure 6, the ABT values were significantly (P < 0.05) positively correlated with AiPC4 expression levels in all the tissues examined, with Pearson’s correlation coefficients ranging from 0.675 to 0.962. Moreover, relatively higher correlation coefficients (r > 0.850) were observed in the gonads (ovary, 0.880; testis, 0.962) and in temperature sensitive tissues (mantle, 0.882; heart, 0.950). The correlation results indicated that AiPC4 might be responsible for the regulation of thermal tolerance in bay scallops.
Figure 6. Correlation analysis between AiPC4 expression levels in tissues and ABT values of bay scallops [(A) mantle; (B) gill; (C) testis; (D) ovary; (E) kidney; (F) striated muscle; (G) hepatopancreas; (H) heart].
Discussion
Genome-wide association study has been proven to be a useful approach in elucidating the genetic variation associated with complex quantitative traits and in identifying the genes regulating these traits, thus, providing markers or genes for selective breeding projects (Huang et al., 2010; Korte and Farlow, 2013). In the present study, we employed GLM (Dobson and Barnett, 2012; Hastie, 2017) in PLINK software to evaluate the correlation between SNPs and the ABT of the experimental bay scallops. Previous researches demonstrated that, a quantitative trait affected by multiple genetic loci and non-genetic factors, may be controlled both by major variants with large genetic effects and by minor variants with small effects (Frankham, 1996; Li et al., 2017). Corresponding to this common assumption, the Manhattan plot (Figure 1) in the present study assigned the top 14 SNPs significantly associated with thermal tolerance of bay scallop into several genes at the genome-wide level, indicating that these genes may dominantly contribute together to the thermal tolerance trait with large effects. Some candidate gene families that are closely associated with thermal tolerance in scallops have also been identified by homology alignment in other studies, such as GRP94 (94 kDa glucose-regulated protein) (Wang et al., 2018a), TNFR (tumor necrosis factor receptors) (Xing et al., 2016b), HSP70 (heat shock protein 70) (Yang et al., 2014), and metallothionein 1 (Yang et al., 2013). Specifically, four and six SNPs in the promoter region of HSP70 and HSP90 were identified to be associated with heat tolerance of two bay scallop populations, via elevating transcriptional activities of target genes (Yang et al., 2014, 2015). Similar result in Zhikong scallop has also concluded that the polymorphism at locus + 94 of CfHSP22 was associated with heat tolerance of scallop (Yang et al., 2011). Above studies including the present research revealed that high expressions of these candidate genes were positively correlated with thermal tolerance, and the mutation at specific loci of the candidate genes (e.g., the SNPs in the promoter region or SNPs in the CDS) might also be associated with thermal tolerance in scallops. In the present study, we conducted a GWAS analysis to screen the possible loci associated with thermal tolerance at the whole genome-wide scale, rather than applied a homology alignment or PCR-based method to detect one or several candidate genes in other studies. In addition, previous studies reported that most of the SNPs in the promoter region affected the putative transcription factor binding sites and this might be the reason of the emergence of different expression patterns of candidate genes under heat stress. However, the molecular mechanism behind thermal tolerance in scallops remains unclear and needs to be further investigated. Therefore, it can be inferred that thermal tolerance in aquatic bivalves is a complex quantitative trait that is regulated by multiple major genes through distinct molecular mechanisms. Among the top 14 significant SNPs identified in the present study, chr15-5091-20379557-1 had the highest correlation with thermal tolerance in bay scallop and was located in AiPC4. AiPC4 has been shown to be involved in RNA polymerase II transcription, repairing DNA damages and stabilizing genome functions (Wang et al., 2004; Mortusewicz et al., 2016). The GWAS result provided an SNP marker that could potentially be used as a thermal tolerance marker in selective breeding for thermal tolerance improvement in bay scallop.
The results of the GWAS showed that thermal tolerance in bay scallop was not dominantly controlled by only a few major genes/loci, but was also regulated by plenty of minor genes/loci, which agreed with the common assumption that quantitative traits are controlled by a small number of genes with large effects and a large number of genes with small effects (Lu et al., 2013). At a less stringent significance level (P < 0.05), we obtained 1,906 SNPs in 950 unique genes associated with ABT. These genes were used for GO function annotation and KEGG pathway enrichment analyses. GO and KEGG analyses have been widely utilized in omics data to predict potential gene phenotypes (Zhang et al., 2013) and elucidate important molecular mechanisms (Wang et al., 2019; Feng et al., 2020). In the present study, GO function annotation showed that these genes were primarily annotated in cellular processes, membranes, and binding functions, which were similar to the findings of previous studies on thermal tolerance in other species (Caldwell et al., 2016; Mondal et al., 2019). A reduction in the upper thermal limits of Nile tilapia (Oreochromis niloticus) resulted in abnormalities in the cellular processes of erythrocytes (Islam et al., 2020). The reorganization of membrane glycerophospholipids was reported to be crucial for maintaining membrane function in fruit flies (Drosophila melanogaster) exposed to thermal stress (Overgaard et al., 2008). Additionally, molecular binding activity is sensitive to temperature alteration (Alison et al., 1996; Schmidt et al., 2000). Furthermore, KEGG pathway enrichment analysis identified the processes influenced by thermal stress. Thermal tolerance research on the land snail, Theba pisana, revealed that heat shock protein 90 was more likely to be implicated in signal transduction processes that are activated by heat stress (Mizrahi et al., 2015). Additionally, the metabolic rate of organisms represents a typical response to progressive warming (Abele et al., 2001), and pathogen invasion can have an important effect on their host’s thermal tolerance (Bruneaux et al., 2017). We speculated that, under thermal conditions, bay scallops adapt by altering cellular processes, membrane and binding functions in accordance with the thermal levels, resulting in changes in molecular functions, including binding activity, signal transduction processes, and metabolic rate.
Among the top 14 SNPs identified in this study, chr15-5091-20379557-1, which had the highest correlation with thermal tolerance, was mapped to AiPC4. Subsequent investigation of AiPC4 mutants was conducted by comparing the ABT scallops from the two extreme ends (N = 15, 10%) in a cultured population (N > 500). Although polymorphism of the chr15-5091-20379557-1 was not verified, multiple sequence alignment analysis identified two additional significant thermal tolerance-associated SNPs (c.234T < C and c.369A < G) located in CDS 1 of AiPC4, which was close to the location of chr15-5091-20379557-1. Furthermore, we examined the expression levels of AiPC4 in different tissues of bay scallop. High expression levels of the gene were detected in multiple thermal-related tissues (Garofalo et al., 2019; Zhou et al., 2019), thus, indicating the important role of AiPC4 in the hosts thermal tolerance. Temperature is one of the major environmental factors affecting the physiology of aquatic poikilotherms, such as scallops (Zittier et al., 2015). In the present study, AiPC4 was detected in most of the observed tissues, with high levels of expression in tissues involved in thermal response, such as the heart, striated muscle, and mantle; additionally, the level of expression was high in the ovary, which is necessary for gametogenesis. Scallops possess a typical bivalve heart composed of two auricles and one ventricle, and the contraction of the heart drives its circulatory system, which is necessary for survival under normal and stressful conditions (Shumway and Parsons, 2011). The mantle is considered to be the first barrier against pathogens (Xing et al., 2016b) and is exposed to extrinsic factors, making it more sensitive to temperature variations (Su et al., 2011). Furthermore, swimming performance is mainly controlled by contractile actions (maximal velocity and maximal power output) of the striated adductor muscle. However, the functions of the striated adductor muscle in the bay scallop are strongly influenced by temperature (Olson and Marsh, 1993). Extrinsic factors, such as high temperature, accelerate the oxygen consumption rate of organisms, which could lead to the formation of reactive oxygen species (ROS) and, consequently, result in oxidative damages (Guerra et al., 2012). As reported, PC4 initiates transcriptional activation during TFIIA–TFIID-promoter complex formation and is important for efficient transcription elongation (Kaiser et al., 1995; Calvo and Manley, 2014). In mammals, PC4 acts as a multifaceted factor in oxidative DNA repair and metabolism maintenance (Mortusewicz et al., 2008; Garavís and Calvo, 2017). In particular, changes in environmental factors such as thermal stress can lead to oxidative stress, resulting in cell injury caused by increased ROS. Similar results have also been reported in sea and estuarine animals (Diplodus vulgaris, D. sargus, Dicentrarchus labrax, Gobius niger, and Liza ramada) (Madeira et al., 2013), Antarctic bivalve Yoldia eightsii (Abele et al., 2001), and sea anemone Anemonia viridis (Richier et al., 2006). Therefore, we speculated that the relatively high expression level of AiPC4 in scallops under high temperature stress is necessary for the repair of DNA damages resulted from oxidative stress and the maintenance of genome stability and metabolic functions. Moreover, we observed relatively (P < 0.01) higher correlation between ABT values and the expression levels of AiPC4 in the testis (r = 0.962), heart (r = 0.950), mantle (r = 0.882), ovary (r = 0.880), hepatopancreas (r = 0.880) and striated muscle (r = 0.810) of scallops.
The above observations indicated that the expression levels of AiPC4 might be responsible for the regulation of thermal tolerance in bay scallop. Although the molecular mechanism of thermal tolerance in mammals has been extensively studied, the regulatory effect of PC4 in bivalves is poorly understood, thus, making it an area worthy of further research. The findings of this study indicated that PC4 was involved in the regulation of thermal tolerance in bay scallops.
In conclusion, through GWAS analysis, we identified the most significant SNP associated with thermal tolerance in the bay scallop. The identified SNP was directly assigned to the coding region of PC4, which was reported to be an important regulator of DNA damage repair and genome function stabilization. GO function annotation and KEGG pathway enrichment analyses revealed the possible molecular mechanism of thermal tolerance in bay scallops. Expression pattern examination of the target gene, AiPC4, and the correlation analysis of gene expression and ABT values jointly demonstrated the importance of AiPC4 in thermal tolerance in bay scallops. Taken together, these evidences collectively suggest that PC4 plays a pivotal role in the thermal tolerance in bay scallops.
Data Availability Statement
The genetic sequencing data involved in our manuscript has been uploaded and deposited in the GenBank (accession number: PRJNA689862).
Author Contributions
QX, XHu, and ZB conceived and designed the experiments. WL, PL, and JL collected the samples. QX, XZ, JZ, and XHo performed the experiments. XZ, PL, JZ, QZ, and XHo analyzed the data. QX and XZ wrote the manuscript. QX, XZ, and PL revised the manuscript. All authors have read and approved the final manuscript.
Funding
This work was supported by the National Natural Science Foundation of China (31702341), China Agriculture Research System of MOF and MARA, and the National Infrastructure of Fishery Germplasm Resources.
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.650045/full#supplementary-material
Footnotes
References
Abele, D., Tesch, C., Wencke, P., and Pörtner, H. O. (2001). How does oxidative stress relate to thermal tolerance in the Antarctic bivalve Yoldia eightsi? Antarct. Sci. 13, 111–118. doi: 10.1017/s0954102001000189
Alison, D. M., Brown, K., Lightowlers, R., and Hughes, M. A. (1996). A low-temperature-responsive gene from barley encodes a protein with single-stranded nucleic acid-binding activity which is phosphorylated in vitro. Plant Mol. Biol. 30, 947–959. doi: 10.1007/bf00020806
Altschul, S. F., Madden, T. L., Schäffer, A. A., Zhang, J., Zhang, Z., Miller, W., et al. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389–3402. doi: 10.1093/nar/25.17.3389
Bakhmet, I. N. (2005). Electrophysiological study of effects of environmental factors on respiration and cardiac activity in the fry of the Atlantic salmon Salmo salar. J. Evol. Biochem. Physiol. 41, 176–184. doi: 10.1007/s10893-005-0052-0
Bamber, S. D., and Depledge, M. H. (1997). Responses of shore crabs to physiological challenges following exposure to selected environmental contaminants. Aquat. Toxicol. 40, 79–92. doi: 10.1016/s0166-445x(97)00040-4
Bruneaux, M., Visse, M., Gross, R., Pukk, L., Saks, L., and Vasemägi, A. (2017). Parasite infection and decreased thermal tolerance: impact of proliferative kidney disease on a wild salmonid fish in the context of climate change. Funct. Ecol. 31, 216–226. doi: 10.1111/1365-2435.12701
Caldwell, R. B., Braselmann, H., Schoetz, U., Heuer, S., Scherthan, H., and Zitzelsberger, H. (2016). Positive cofactor 4 (PC4) is critical for DNA repair pathway re-routing in DT40 cells. Sci. Rep. 6:28890.
Calvo, O., and Manley, J. L. (2014). The transcriptional coactivator PC4/Sub1 has multiple functions in RNA polymerase II transcription. EMBO J. 24, 1009–1020. doi: 10.1038/sj.emboj.7600575
Casselman, M. T., Anttila, K., and Farrell, A. P. (2012). Using maximum heart rate as a rapid screening tool to determine optimum temperature for aerobic scope in Pacific salmon Oncorhynchus spp. J. Fish Biol. 80, 358–377. doi: 10.1111/j.1095-8649.2011.03182.x
Chen, J., Mai, K., Ma, H., Wang, X., Deng, D., Liu, X., et al. (2007a). Effects of dissolved oxygen on survival and immune responses of scallop (Chlamys farreri Jones et Preston). Fish Shellfish Immun. 22, 272–281. doi: 10.1016/j.fsi.2006.06.003
Chen, J., and Chen, W. (2000). Salinity tolerance of Haliotis diversicolor supertexta at different salinity and temperature levels. Aquaculture 181, 191–203. doi: 10.1016/s0044-8486(99)00226-4
Chen, M., Yang, H., Delaporte, M., Zhao, S., and Xing, K. (2007b). Immune responses of the scallop Chlamys farreri after air exposure to different temperatures. J. Exp. Mar. Biol. Ecol. 345, 52–60. doi: 10.1016/j.jembe.2007.01.007
Chen, N., Luo, X., Gu, Y., Han, G., Dong, Y., You, W., et al. (2016). Assessment of the thermal tolerance of abalone based on cardiac performance in Haliotis discus hannai, H. gigantea and their interspecific hybrid. Aquaculture 465, 258–264. doi: 10.1016/j.aquaculture.2016.09.004
Cheng, D., Liu, H., Zhang, H., Tan, K., Ye, T., Ma, H., et al. (2020). Effects of thermal stress on mortality and HSP90 expression levels in the noble scallops Chlamys nobilis with different total carotenoid content. Cell Stress Chaperones 25, 105–117. doi: 10.1007/s12192-019-01052-5
Dahlhoff, E., O’Brien, J., Somero, G. N., and Vetter, R. D. (1991). Temperature effects on mitochondria from hydrothermal vent invertebrates: evidence for adaptation to elevated and variable habitat temperatures. Physiol. Zool. 64, 1490–1508. doi: 10.1086/physzool.64.6.30158226
Das, C., Hizume, K., Batta, K., Kumar, B. R. P., Gadad, S. S., Ganguly, S., et al. (2006). Transcriptional coactivator PC4, a chromatin-associated protein, induces chromatin condensation. Mol. Cell. Biol. 26, 8303–8315. doi: 10.1128/mcb.00887-06
Dobson, A. J., and Barnett, A. G. (2012). An introduction to generalized linear models, third edition. Pharm. Stat. 11:272. doi: 10.1002/pst.1502
FAO (2019). Fisheries and Aquaculture Stastics 2017. Rome: Food and Agriculture Organization of the United Nations.
Feng, X., Zhu, L., Qin, Z., Mo, X., and Li, S. (2020). Temporal transcriptome change of Oncomelania hupensis revealed by Schistosoma japonicum invasion. Cell Biosci. 10:58.
Frankham, R. (1996). Introduction to quantitative genetics (4th edition). Am. J. Hum. Genet. 46:1231.
Fu, X., Dou, J., Mao, J., Su, H., Jiao, W., Zhang, L., et al. (2013). RADtyping: an integrated package for accurate de novo codominant and dominant RAD genotyping in mapping populations. PLoS One 8:e79960. doi: 10.1371/journal.pone.0079960I
Gagnaire, B., Frouin, H., Moreau, K., Thomas-Guyon, H., and Renault, T. (2006). Effects of temperature and salinity on haemocyte activities of the Pacific oyster, Crassostrea gigas (Thunberg). Fish Shellfish Immun. 20, 536–547. doi: 10.1016/j.fsi.2005.07.003
Garavís, M., and Calvo, O. (2017). Sub1/PC4, a multifaceted factor: from transcription to genome stability. Curr. Genet. 63, 1023–1035. doi: 10.1007/s00294-017-0715-6
Garofalo, F., Santovito, G., and Amelio, D. (2019). Morpho-functional effects of heat stress on the gills of Antarctic T. bernacchii and C. hamatus. Mar. Pollut. Bull. 141, 194–204. doi: 10.1016/j.marpolbul.2019.02.048
Geng, X., Liu, S., Yao, J., Bao, L., Zhang, J., Li, C., et al. (2016). A genome-wide association study identifies multiple regions associated with head size in catfish. G3 6, 3389–3398. doi: 10.1534/g3.116.032201
Granados-Amores, A., Campa-Córdova, ÁI., Araya, R., Mazón-Suástegui, J. M., and Saucedo, P. E. (2012). Growth, survival and enzyme activity of lions-paw scallop (Nodipecten subnodosus) spat treated with probiotics at the hatchery. Aquac. Res. 43, 1335–1343. doi: 10.1111/j.1365-2109.2011.02937.x
Guerra, C., Zenteno-Savín, T., Maeda-Martínez, A. N., Philipp, E. E. R., and Abele, D. (2012). Changes in oxidative stress parameters in relation to age, growth and reproduction in the short-lived catarina scallop Argopecten ventricosus reared in its natural environment. Comp. Biochem. Physiol. A Mol. Integr. Physiol. 162, 421–430. doi: 10.1016/j.cbpa.2012.04.018
Heilmayer, O., Honnen, C., Jacob, U., Chiantore, M., Cattaneo-Vietti, R., and Brey, T. (2005). Temperature effects on summer growth rates in the Antarctic scallop, Adamussium colbecki. Polar Biol. 28, 523–527. doi: 10.1007/s00300-005-0716-7
Helgeland, H., Sodeland, M., Zoric, N., Torgersen, J. S., Grammes, F., Lintig, J. V., et al. (2019). Genomic and functional gene studies suggest a key role of beta-carotene oxygenase 1 like (bco1l) gene in salmon flesh color. Sci. Rep. 9:20061.
Hu, X., Bao, Z., Hu, J., Shao, M., Zhang, L., Bi, K., et al. (2006). Cloning and characterization of tryptophan 2, 3-dioxygenase gene of Zhikong scallop Chlamys farreri (Jones and Preston 1904). Aquac. Res. 37, 1187–1194. doi: 10.1111/j.1365-2109.2006.01546.x
Huang, X., Wei, X., Sang, T., Zhao, Q., Feng, Q., Zhao, Y., et al. (2010). Genome-wide association studies of 14 agronomic traits in rice landraces. Nat. Genet. 42:961.
Islam, S. M., Zahangir, M. M., Jannat, R., Hasan, M. N., Suchana, S. A., Rohani, M. F., et al. (2020). Hypoxia reduced upper thermal limits causing cellular and nuclear abnormalities of erythrocytes in Nile tilapia, Oreochromis niloticus. J. Therm. Biol. 90:102604. doi: 10.1016/j.jtherbio.2020.102604
Jiao, W., Fu, X., Dou, J., Li, H., Su, H., Mao, J., et al. (2014). High-resolution linkage and quantitative trait locus mapping aided by genome survey sequencing: building up an integrative genomic framework for a bivalve mollusc. DNA Res. 21, 85–101. doi: 10.1093/dnares/dst043
Jin, Y., Zhou, T., Geng, X., Liu, S., Chen, A., Yao, J., et al. (2017). A genome-wide association study of heat stress-associated SNPs in catfish. Anim. Genet. 48, 233–236. doi: 10.1111/age.12482
Jonasson, J. P., Thorarinsdottir, G., Eiriksson, H., Solmundsson, J., and Marteinsdottir, G. (2007). Collapse of the fishery for Iceland scallop (Chlamys islandica) in Breidafjordur, West Iceland. ICES J. Mar. Sci. 64, 298–308. doi: 10.1093/icesjms/fsl028
Kaiser, K., Stelzer, G., and Meisterernst, M. (1995). The coactivator p15 (PC4) initiates transcriptional activation during TFIIA-TFIID-promoter complex formation. EMBO J. 14, 3520–3527. doi: 10.1002/j.1460-2075.1995.tb07358.x
Korte, A., and Farlow, A. (2013). The advantages and limitations of trait analysis with GWAS: a review. Plant Methods 9:29. doi: 10.1186/1746-4811-9-29
Li, H., Su, G., Li, J., and Bao, Z. (2017). An efficient unified model for genome-wide association studies and genomic selection. Genet. Sel. Evol. 49:64.
Li, R., Yu, C., Li, Y., Lam, T. W., Yiu, S. M., Kristiansen, K., et al. (2009). SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics 25, 1966–1967. doi: 10.1093/bioinformatics/btp336
Li, X., Wang, S., Xun, X., Zhang, M., Wang, S., Li, H., et al. (2019a). A carotenoid oxygenase is responsible for muscle coloration in scallop. Biochim. Biophys. Acta Mol. Cell Biol. Lipids 1864, 966–975. doi: 10.1016/j.bbalip.2019.03.003
Li, Y., Zhang, L., Li, R., Zhang, M., Li, Y., Wang, H., et al. (2019b). Systematic identification and validation of the reference genes from 60 RNA-Seq libraries in the scallop Mizuhopecten yessoensis. BMC Genomics. 20:288. doi: 10.1186/s12864-019-5661-x
Liu, X., Li, C., Chen, M., Liu, B., Yan, X., Ning, J., et al. (2020). Draft genomes of two Atlantic bay scallop subspecies Argopecten irradians irradians and A. i. concentricus. Sci. Data 7:99.
Lu, D., Miller, S., Sargolzaei, M., Kelly, M., Vander, V. G., Caldwell, T., et al. (2013). Genome-wide association analyses for growth and feed efficiency traits in beef cattle. J. Anim. Sci. 91, 3612–3633. doi: 10.2527/jas.2012-5716
Lubet, P. E. (1983). Experimental studies on the action of temperature on the reproductive activity of the mussel (Mytilus edulis L., Mollusca, Lamellibranchia). J. Molluscan Stud. 12, 100–105.
Lv, J., Jiao, W., Guo, H., Liu, P., Wang, R., Zhang, L., et al. (2018). HD-Marker: a highly multiplexed and flexible approach for targeted genotyping of more than 10,000 genes in a single-tube assay. Genome Res. 28, 1919–1930. doi: 10.1101/gr.235820.118
Madeira, D., Narciso, L., Cabral, H. N., Vinagre, C., and Diniz, M. S. (2013). Influence of temperature in thermal and oxidative stress responses in estuarine fish. Comp. Biochem. Physiol. A Mol. Integr. Physiol. 166, 237–243. doi: 10.1016/j.cbpa.2013.06.008
Maeda-Martínez, A. M., Obregón-Barboza, V., Navarrete-Del Toro, M. A., Obregón-Barboza, H., and García-Carreño, F. L. (2000). Trypsin-like enzymes from two morphotypes of the ‘living fossil’ Triops (Crustacea: Branchiopoda: Notostraca). Comp. Biochem. Physiol. B Biochem. Mol. Biol. 126, 317–323. doi: 10.1016/s0305-0491(00)00183-8
Martínez, G., and Pérez, H. (2003). Effect of different temperature regimes on reproductive conditioning in the scallop Argopecten purpuratus. Aquaculture 228, 153–167. doi: 10.1016/s0044-8486(03)00321-1
Mao, J., Lv, J., Miao, Y., Sun, C., Hu, L., Zhang, R., et al. (2013). Development of a rapid and efficient method for non-lethal DNA sampling and genotyping in scallops. PLoS One 8:e68096. doi: 10.1371/journal.pone.0068096
Mason, J. (1958). The breeding of the scallop, Pecten maximus (L.), in Manx waters. J. Mar. Biol. Assoc. 37, 653–671. doi: 10.1017/s0025315400005701
Mizrahi, T., Goldenberg, S., Heller, J., and Arad, Z. (2015). Geographic variation in thermal tolerance and strategies of heat shock protein expression in the land snail Theba pisana in relation to genetic structure. Cell Stress Chaperones 21, 219–238. doi: 10.1007/s12192-015-0652-6
Moffitt, C. M., and Cajas-Cano, L. (2014). Blue growth: the 2014 FAO state of world fisheries and aquaculture. Fisheries 39, 552–553. doi: 10.1080/03632415.2014.966265
Mondal, P., Saleem, S., Sikder, S., Kundu, T. K., Biswas, S. C., and Roy, S. (2019). Multifunctional transcriptional coactivator PC4 is a global co-regulator of p53-dependent stress response and gene regulation. J. Biochem. 166, 403–413. doi: 10.1093/jb/mvz050
Mortusewicz, O., Roth, W., Li, N., Cardoso, M. C., Meisterernst, M., and Leonhardt, H. (2008). Recruitment of RNA polymerase II cofactor PC4 to DNA damage sites. J. Cell Biol. 183, 769–776. doi: 10.1083/jcb.200808097
Mortusewicz, O., Evers, B., and Helleday, T. (2016). PC4 promotes genome stability and DNA repair through binding of ssDNA at DNA damage sites. Oncogene 35, 761–770. doi: 10.1038/onc.2015.135
Ning, X., Li, X., Wang, J., Zhang, X., Kong, L., Meng, D., et al. (2019). Genome-wide association study reveals E2F3 as the candidate gene for scallop growth. Aquaculture 511:734216. doi: 10.1016/j.aquaculture.2019.734216
Olson, J. M., and Marsh, R. L. (1993). Contractile properties of the striated adductor muscle in the bay scallop Argopecten irradians at several temperatures. J. Exp. Biol. 176, 175–193. doi: 10.1242/jeb.176.1.175
Overgaard, J., Tomčala, A., Sørensen, J. G., Holmstrup, M., Krogh, P. H., Šimek, P., et al. (2008). Effects of acclimation temperature on thermal tolerance and membrane phospholipid composition in the fruit fly Drosophila melanogaster. J. Insect Physiol. 54, 619–629. doi: 10.1016/j.jinsphys.2007.12.011
Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A. R., Bender, D., et al. (2007). PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559–575. doi: 10.1086/519795
Richards, G. P., Watson, M. A., Needleman, D. S., Church, K. M., and Häse, C. C. (2015). Mortalities of eastern and Pacific oyster larvae caused by the pathogens Vibrio coralliilyticus and Vibrio tubiashii. Appl. Environ. Microbiol. 81, 292–297. doi: 10.1128/aem.02930-14
Rojas, R., Miranda, C. D., Santander, J., and Romero, J. (2016). First report of Vibrio tubiashii associated with a massive larval mortality event in a commercial hatchery of scallop Argopecten purpuratus in Chile. Front. Microbiol. 7:1473. doi: 10.3389/fmicb.2016.01473
Roldan-Carrillo, L. M., Maeda-Martinez, A. N., Masso-Rojas, A., and Sicard-Gonzalez, M. T. (2005). Salinity tolerance and resistance of the Pacific lion’s paw scallop (Nodipecten subnodosus) and the relationships with species distribution and density in a coastal lagoon. J. Shellfish Res. 24, 353–362. doi: 10.2983/0730-8000(2005)24[353:atarot]2.0.co;2
Richier, S., Sabourault, C., Courtiade, J., Zucchini, N., Allemand, D., and Furla, P. (2006). Oxidative stress and apoptotic events during thermal stress in the symbiotic sea anemone, Anemonia viridis. FEBS J. 273, 4186–4198. doi: 10.1111/j.1742-4658.2006.05414.x
Sambrook, J., Fritsch, E. F., and Maniatis, T. (1989). Molecular Cloning a Laboratory Manual, 2nd Edn. New York, NY: Cold Spring Harbor Laboratory Press.
Schmidt, M., Ding, H., Ramamurthy, V., Mukerji, I., and Oliver, D. (2000). Nucleotide binding activity of SecA homodimer is conformationally regulated by temperature and altered by prlD and azi mutations. J. Biol. Chem. 272:15440. doi: 10.1074/jbc.m000605200
Secor, S. M. (2009). Specific dynamic action: a review of the postprandial metabolic response. J. Comp. Physiol. B Biochem. Syst. Environ. Physiol. 179, 1–56. doi: 10.1007/s00360-008-0283-7
Shumway, S. E., and Parsons, J. (2011). Scallops: Biology, Ecology and Aquaculture. Amsterdam: Elsevier.
Sicard, M. T., Maeda-Martinez, A. N., Ormart, P., Reynoso-Granados, T., and Carvalho, L. (1999). Optimum temperature for growth in the catarina scallop (Argopecten ventricosus-circularis, Sowerby II, 1842). J. Shellfish Res. 18, 385–392.
Somero, G. N. (2002). Thermal physiology and vertical zonation of intertidal animals: optima, limits, and costs of living. Integr. Comp. Biol. 42, 780–789. doi: 10.1093/icb/42.4.780
Sonesson, A. K., and Meuwissen, T. H. (2009). Testing strategies for genomic selection in aquaculture breeding programs. Genet. Sel. Evol. 41:37.
Steeves, L. E., Filgueira, R., Guyondet, T., Chassé, J., and Comeau, L. (2018). Past, present, and future: performance of two bivalve species under changing environmental conditions. Front. Mar. Sci. 5:184. doi: 10.3389/fmars.2018.00184
Stenseng, E., Braby, C. E., and Somero, G. N. (2005). Evolutionary and acclimation-induced variation in the thermal limits of heart function in congeneric marine snails (Genus Tegula): implications for vertical zonation. Biol. Bull. 208, 138–144. doi: 10.2307/3593122
Stillman, J., and Somero, G. (1996). Adaptation to temperature stress and aerial exposure in congeneric species of intertidal porcelain crabs (genus Petrolisthes): correlation of physiology, biochemistry and morphology with vertical distribution. J. Exp. Biol. 199, 1845–1855. doi: 10.1242/jeb.199.8.1845
Su, J., Qiu, L., Li, L., Liu, L., Wang, L., and Siva, V. S. (2011). cDNA cloning and characterization of a new member of the tumor necrosis factor receptor family gene from scallop, Chlamys farreri. Mol. Biol. Rep. 38, 4483–4490. doi: 10.1007/s11033-010-0578-0
Wang, J. Y., Sarker, A. H., Cooper, P. K., and Volkert, M. R. (2004). The single-strand DNA binding activity of human PC4 prevents mutagenesis and killing by oxidative DNA damage. Mol. Cell. Biol. 24, 6084–6093. doi: 10.1128/mcb.24.13.6084-6093.2004
Wang, L., Zhang, H., Song, L., and Guo, X. (2007). Loss of allele diversity in introduced populations of the hermaphroditic bay scallop Argopecten irradians. Aquaculture 271, 252–259. doi: 10.1016/j.aquaculture.2007.06.020
Wang, S., Liu, P., Lv, J., Li, Y., Cheng, T., Zhang, L., et al. (2016). Serial sequencing of isolength RAD tags for cost-efficient genome-wide profiling of genetic and epigenetic variations. Nat. Protoc. 11:2189. doi: 10.1038/nprot.2016.133
Wang, S., Li, X., Li, T., Wang, H., Zhang, X., Lou, J., et al. (2018a). The GRP94 gene of Yesso scallop (Patinopecten yessoensis): characterization and expression regulation in response to thermal and bacterial stresses. Fish Shellfish Immun. 80, 443–451. doi: 10.1016/j.fsi.2018.06.018
Wang, S., Meyer, E., Mckay, J. K., and Matz, M. V. (2012). 2b-RAD: a simple and flexible method for genome-wide genotyping. Nat. Methods 9, 808–810. doi: 10.1038/nmeth.2023
Wang, X., Liu, S., Jiang, C., Geng, X., Zhou, T., Li, N., et al. (2017). Multiple across-strain and within-strain QTLs suggest highly complex genetic architecture for hypoxia tolerance in channel catfish. Mol. Genet. Genomics 292, 63–76. doi: 10.1007/s00438-016-1256-2
Wang, Y., Liu, L., and Chen, Z. (2019). Transcriptome profiling of cervical cancer cells acquired resistance to cisplatin by deep sequencing. Artif. Cells Nanomed. Biotechnol. 47, 2820–2829. doi: 10.1080/21691401.2019.1637882
Wang, Y., Sun, G., Zeng, Q., Chen, Z., Hu, X., Li, H., et al. (2018b). Predicting growth traits with genomic selection methods in Zhikong scallop (Chlamys farreri). Mar. Biotechnol. 20, 769–779. doi: 10.1007/s10126-018-9847-z
Weinstein, R. B., and Somero, G. N. (1998). Effects of temperature on mitochondrial function in the Antarctic fish Trematomus bernacchii. J. Comp. Physiol. B Biochem. Syst. Environ. Physiol. 168, 190–196. doi: 10.1007/s003600050136
Wu, Y., Zhou, Z., Pan, Y., Zhao, J., Bai, H., Chen, B., et al. (2021). GWAS identified candidate variants and genes associated with acute heat tolerance of large yellow croaker. Aquaculture 540:736696. doi: 10.1016/j.aquaculture.2021.736696
Xing, Q., Li, Y., Guo, H., Yu, Q., Huang, X., Wang, S., et al. (2016a). Cardiac performance: a thermal tolerance indicator in scallops. Mar. Biol. 163:244.
Xing, Q., Yu, Q., Dou, H., Wang, J., Li, R., Ning, X., et al. (2016b). Genome-wide identification, characterization and expression analyses of two TNFRs in Yesso scallop (Patinopecten yessoensis) provide insight into the disparity of responses to bacterial infections and heat stress in bivalves. Fish Shellfish Immun. 52, 44–56. doi: 10.1016/j.fsi.2016.03.010
Xing, Q., Zhang, L., Li, Y., Zhu, X., Li, Y., Guo, H., et al. (2019). Development of novel cardiac indices and assessment of factors affecting cardiac activity in a Bivalve Mollusc Chlamys farreri. Front. Physiol. 10:293. doi: 10.3389/fphys.2019.00293
Xing, Q., Wang, J., Hu, L., Sun, Y., Huang, X., Zhang, L., et al. (2021). Seasonal variation of the thermal tolerance indicator ABT and the development of a rapid detection method in scallop Chlamys farreri. Aquaculture 531:735960. doi: 10.1016/j.aquaculture.2020.735960
Xiao, J., Ford, S. E., Yang, H., Zhang, G., Zhang, F., and Guo, X. (2005). Studies on mass summer mortality of cultured Zhikong scallops (Chlamys farreri Jones et Preston) in China. Aquaculture 250, 602–615. doi: 10.1016/j.aquaculture.2005.05.002
Yang, C., Wang, L., Jiang, Q., Wang, J., Yue, F., Zhang, H., et al. (2013). The polymorphism in the promoter region of metallothionein 1 is associated with heat tolerance of scallop Argopecten irradians. Gene 526, 429–436. doi: 10.1016/j.gene.2013.05.033
Yang, C., Wang, L., Liu, C., Zhou, Z., Xin, Z., and Song, L. (2015). The polymorphisms in the promoter of hsp90 gene and their association with heat tolerance of bay scallop. Cell Stress Chaperones 20, 297–308. doi: 10.1007/s12192-014-0546-z
Yang, C., Wang, L., Wang, J., Jiang, Q., Qiu, L., Zhang, H., et al. (2014). The Polymorphism in the promoter of HSP70 gene is associated with heat tolerance of two congener endemic bay scallops (Argopecten irradians irradians and A. i. concentricus). PLoS One 9:e102332. doi: 10.1371/journal.pone.0102332
Yang, C., Zhang, L., Wang, L., Zhang, H., Qiu, L., Siva, V. S., et al. (2011). The gln32lys polymorphism in hsp22 of Zhikong scallop Chlamys farreri is associated with heat tolerance. PLoS One 6:e28564. doi: 10.1371/journal.pone.0028564
Yu, F., Peng, W., Tang, B., Zhang, Y., Wang, Y., Gan, Y., et al. (2021). A genome-wide association study of heat tolerance in Pacific abalone based on genome resequencing. Aquaculture 536:736436. doi: 10.1016/j.aquaculture.2021.736436
Zhang, J., Kassian, K., Tomáš, F., and Alexandros, S. (2014). PEAR: a fast and accurate illumina paired-end reAd mergeR. Bioinformatics 30, 614–620. doi: 10.1093/bioinformatics/btt593
Zhang, T., Jiang, M., Chen, L., Niu, B., and Cai, Y. (2013). Prediction of gene phenotypes based on GO and KEGG pathway enrichment scores. Biomed. Res. Int. 2013:870795.
Zhao, L., Li, Y., Li, Y., Yu, J., Liao, H., Wang, S., et al. (2017). A genome-wide association study identifies the genomic region associated with shell color in Yesso scallop, Patinopecten yessoensis. Mar. Biotechnol. 19, 301–309. doi: 10.1007/s10126-017-9751-y
Zhou, C., Zhou, P., Ren, Y., Cao, L., and Wang, J. (2019). Physiological response and miRNA-mRNA interaction analysis in the head kidney of rainbow trout exposed to acute heat stress. J. Therm. Biol. 83, 134–141. doi: 10.1016/j.jtherbio.2019.05.014
Zhou, T., Liu, S., Geng, X., Jin, Y., Jiang, C., Bao, L., et al. (2017). GWAS analysis of QTL for enteric septicemia of catfish and their involved genes suggest evolutionary conservation of a molecular mechanism of disease resistance. Mol. Genet. Genomics 292, 231–242. doi: 10.1007/s00438-016-1269-x
Zittier, Z. M. C., Bock, C., Lannig, G., and Pörtner, H. O. (2015). Impact of ocean acidification on thermal tolerance and acid-base regulation of Mytilus edulis (L.) from the North Sea. J. Exp. Mar. Bio. Ecol. 473, 16–25. doi: 10.1016/j.jembe.2015.08.001
Keywords: genome-wide association study, Argopecten irradians irradians, thermal tolerance, single nucleotide polymorphism, transcriptional coactivator p15
Citation: Zhu X, Liu P, Hou X, Zhang J, Lv J, Lu W, Zeng Q, Huang X, Xing Q and Bao Z (2021) Genome-Wide Association Study Reveals PC4 as the Candidate Gene for Thermal Tolerance in Bay Scallop (Argopecten irradians irradians). Front. Genet. 12:650045. doi: 10.3389/fgene.2021.650045
Received: 06 January 2021; Accepted: 28 June 2021;
Published: 19 July 2021.
Edited by:
Iker Irisarri, University of Göttingen, GermanyReviewed by:
Changwei Shao, Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences (CAFS), ChinaTimothy D. Leeds, United States Department of Agriculture (USDA), United States
Copyright © 2021 Zhu, Liu, Hou, Zhang, Lv, Lu, Zeng, Huang, Xing and Bao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qiang Xing, qiangxing@ouc.edu.cn
†These authors have contributed equally to this work