Eun Soo Noh1
Sathiyamoorty Subramaniyam2
Sunghyun Cho2
Young-Ok Kim1Choul-Ji Park3Jeong-Ho Lee4
Bo-Hye Nam1
Younhee Shin2*- 1Biotechnology Research Division, National Institute of Fisheries Science, Geoje, Republic of Korea
- 2Research and Development Center, Yongin-si, Gyeonggi-do, Republic of Korea
- 3Fisheries Seed and Breeding Research Institute, National Institute of Fisheries Science, Busan, Republic of Korea
- 4Fish Genetics and Breeding Research Center, National Institute of Fisheries Science, Geoje, Republic of Korea
Introduction
Abalone (Haliotis discus hannai) is an expensive seafood in Asian countries and an important species in fishery and mariculture industries, generating marginal revenue for the Chinese and South Korean economies. South Korea is the second largest producer of abalones, followed by China (Nam et al., 2016). Thus far, approximately 70 species of Haliotis have been discovered worldwide; among them, 7 have commercial importance and 6 are naturally distributed in South Korea (i.e., Haliotis discus hannai, Haliotis discus, Haliotis madaka, Haliotis gigantea, Haliotis diversicolor, and Haliotis diversicolor supertexta) (Adachi and Okumura, 2012). In particular, the species Haliotis discus hannai is farmed widely in coastal regions of South Korea (Im and Kim, 2020). Two major factors affect abalone production in the sea: overfishing, due to its high market value, and increased atmospheric CO2, resulting in rising sea temperatures (Peter, 2016). While overfishing can be addressed by imposing strict laws and establishing a mariculture system, mitigating the changes in sea temperature is more difficult. Moreover, temperature fluctuations in sea water cause high mortality in marine cage-based abalone cultivation, particularly in the summer in coastal regions of South Korea (Kang et al., 2019). Thus, genetic/genome-assisted breeding could be a reasonable solution to increase abalone production in natural sea and mariculture systems. A draft reference genome of abalone is available, which could help determine genotypes for specific phenotypic traits (Nam et al., 2017). Many scientific reports have addressed the establishment of molecular datasets related to the physiological process of heat resistant traits, i.e., transcriptome (Tripp-Valdez et al., 2019; Kyeong et al., 2020; Kim et al., 2021), proteome (Kang et al., 2019), and metabolome (Xu et al., 2020) analyses. These datasets can be used to elucidate preliminary gene markers such as heat shock proteins (HSPs), which function in transcription and translation in abalone (Kyeong et al., 2020). Heat stress also alters energy metabolism and increases susceptibility to various pathogens, such as Vibrio parahemolyticus (Nam et al., 2016; Crosson et al., 2020), affecting the reproduction and growth of abalones (Swezey et al., 2020) as well as the metabolic rate in the digestive tract (Frederick et al., 2022). To our knowledge, no genotyping studies or datasets are available for Haliotis discus hannai other than a population assessment (Nam et al., 2021). In this study, we used genotyping-by-sequencing (GBS) to observe the genotypes associated with heat stress, and establish genotypic chip and machine learning (ML)-based prediction models to predict heat-sensitive abalones for breeding purposes with 96 single-nucleotide polymorphisms (SNPs).
Significance of the data
Genotype data was generated for heat-resistant and -sensitive abalones from three different populations. Initially, 96 SNPs selected from the GBS dataset and those used to build the targeted Fluidigm chip were used for validation in two additional populations. These datasets were subjected to ML methods to develop a predictive model for heat-resistant and -sensitive abalone phenotypes. This dataset could be used to generate a genetic library for genome-assisted breeding for abalone.
Materials and methods
Experimental design for GBS and the fluidigm
A total of 400 abalones were selected from the Genetic and Breeding Research Center, NIFS, in Geoje, South Korea, and from a commercial farm in Wando and Heanam, South Korea (Supplementary Figure S1). The average shell length was 54.69 ± 3.78 mm, the shell width was 37.09 ± 2.63 mm, and the body weight was 16.44 ± 3.87 g. The experiment lasted 20 days (23 September 2016, to 12 October 2016). The abalones were maintained in a tank (1.2 × 3 × 0.8 m) with a constant flow of seawater at the Genetic and Breeding Research Center. The temperature was maintained at ≤ 24°C for the first 7 days to acclimate. The temperature was increased by 1°C per day for the next 7 days, until reaching a maximum temperature of 31°C, using an Aquatron system (Yoowon Electronics, Seoul, South Korea). During this time, the dissolved oxygen level was maintained at 7.8 ± 0.5 mg/L. The temperature was maintained at 31°C for an additional 6 days, during which the abalone mortality rate of the general population was 50%. Dead abalones were collected immediately for sampling. Foot muscles were collected and stored in ethanol. The heat resistant of the abalones was measured by survival time. The animals were categorized into two groups: heat-resistant and heat-sensitive (Figure 1A). Finally, 156 heat-resistant and 107 heat-senstive abalones were randomly selected from F4 population for sequencing. A similar experimental procedure was carried out for another two populations in 2019 and 2020 as shown in Table 1. In total, 2,282 samples were included from all three studies (Table 1; Supplementary Table S3).
FIGURE 1. Sampling of machine learning predicted summaries. (A). Heat resistance experiment for sample selection and defining the class as case and control. (B). Sampling locations in South Korea. (C). Sample diversity observed with 96 SNPs based on principal component analysis (PCA). (D). Machine learning models predicted the accuracy from the pooled dataset.
TABLE 1. Selected SNP-assisted machine learning based prediction summary from three different combinations of datasets.
GBS library preparation and sequencing
Total genomic DNA from the 263 samples was extracted from muscle tissue, quantified, and normalized to 20 ng/μL. The DNA (200 ng) was digested with 8 U of high-fidelity PstI at 37°C for 2 h and heated to 65°C for 20 min to inactivate the enzyme. Six DNA libraries were constructed for GBS as described previously (Nam et al., 2021), pooled, and amplified by multiplex polymerase chain reaction (PCR). The products were purified using a QIAquick PCR Purification Kit (Qiagen, Hilden, Germany) and the distribution of fragment sizes was evaluated with a BioAnalyzer 2,100 instrument (Agilent Technologies, Santa Clara, CA, United States). The GBS libraries were sequenced with the Illumina NextSeq500 platform (San Diego, CA, United States) using 150-bp single reads in DNA Link, the authorized service provider. Library preparation has been described previously (Nam et al., 2021). A summary is provided in Supplementary Table S1.
Variant calling, SNP selection, and estimate association
Sequences were subjected to quality and adapter trimming with Trimmomatic 0.32 using the following parameter settings: leading, 5; trailing, 5; sliding window, 4:15; and min, 30 (Bolger et al., 2014). The processed reads were mapped to the abalone reference genome (Nam et al., 2017) using Bowtie2 v.2.2.8 (Langmead and Salzberg, 2012) and variant calling was performed with the Haplotype caller in the Genome Analysis Toolkit (GATK) (McKenna et al., 2010). SNPs were selected with GATK parameters, i.e., normalized quality score ≥2 and mapping quality ≥40. Additionally, the missing genotypes were input with Beagle method v.4.1 (Browning et al., 2021). SNPs were annotated with SnpEff v.4.2 (Cingolani et al., 2012). Finally, high-quality SNPs were selected with the following steps: 1) bi-allelic sites, 2) genotyping rate of the samples at each variable site ≥90%, 3) minor allele frequency (MAF) > 5%, and 4) Hardy-Weinberg equilibrium (HWE) < 0.001 using PLINK1.9 (Purcell et al., 2007). The selected SNPs were subjected to population stratification with the STRUCTURE algorithm (Pritchard et al., 2000) with a K range of 1–7 with 10,000 iterations. The association between genotype and phenotype was estimated as follows. The samples were classified as heat-resistant (Case) or -sensitive (Control). Dead abalones were considered sensitive and the rest were considered heat-resistant. Features such as fixation index (Fst) and genomic nucleotide diversity (π) were calculated with VCFtools v. 0.1.3 (Danecek et al., 2011). Reduction of diversity (ROD) was determined by the following equation [ROD = 1 − (π of case/π of control)] within a 10 kb window size. The region at which the ROD >0.8 was defined as the selective sweep. The SNP ran from −5 kb to +5 kb in the gene region until the end of the gene. SNPs were considered to be significantly associated with the trait when p < 0.01; analyses were done using PLINK with the—assoc function.
Targeted fluidigm chip design
A Fluidigm chip was designed with 96 SNP markers from the GBS dataset-assisted genome wide association study (GWAS). The target SNP genotyping chip was constructed with a Fluidigm 96.96 dynamic array integrated fluidic circuit (IFC) using an Adventa sample ID genotyping panel. For chip design, the primers for each SNP were selected as 100 bp of the flanking regions. Primers such as allele-specific primers (ASPs), locus-specific primers (LSPs), and specific target amplification primers (STAs) were designed using the Fluidigm SNP type assay protocol. The PCR cocktail was prepared with ASPs, LSPs, and STAs according to the manufacturer’s protocol along with the high-quality DNA prepared from each sample. The samples were loaded in a 96-well plate (12 columns x 8 rows) and subjected to the SNPtype 96X96 thermal cycling protocol with a FC1 PCR cycler to detect fluorescence by Biomark HD and processed with the SNP Trace™ Panel Analysis tool in the SNP genotyping analysis software. Genomic DNA (gDNA) quality was assessed by the concentration (ng/μLS) of each sample, which was measured with a Biotek Epoch spectrometer at 260/280 nm. All experimental protocols were performed by the TNT research service provider in Anyang, South Korea. A detailed summary is provided in Supplementary Table S2.
Machine learning approach to predict phenotype
To assess the classification potential of selected SNPs, 7 ML models were used: AdaBoost (AB), Bagged Tree (BT), Generalized Boosted Regression (GBR), Boosted Logistic Regression (BLR), partial least squares (PLS), Random Forest (RF), and support vector machine with Linear Kernel (SVM-LK). In this study, we generated three datasets (datasets 1–3). We developed the basic models using a pooled dataset with individual validation datasets to understand the association of SNPs with three different populations. ML was performed using the ‘train’ in function in the caret package in R software with tenfold cross validation (ver. 3.3; R Development Core Team, Vienna, Austria) (Zhao, 2014; Emir et al., 2016). Model assessments were performed using parameters such as accuracy, kappa, sensitivity, specificity, pro pred value, negative pred value, precision, recall, F1, prevalence, and balance accuracy as described previously (Malik et al., 2022).
Preliminary analysis report
In total, 185.9 GB of sequence was generated by GBS from 263 samples, and 81.59% of the reads were mapped to the abalone genome (Supplementary Figure S2). The mapped reads covered approximately 3% of the genome. From the mapped reads, 16,119 high-quality SNPs were obtained from 232,231 called SNPs, as illustrated in Supplementary Figure S3. Among those, 96 SNPs were selected with the metrics described in the Materials and Methods. In summary, 18 markers were selected using the ROD, 32 SNP markers were selected using the ODD score, and 46 markers were selected from the GWAS. The majority of the selected SNPs were present in intergenic regions and upstream of coding regions (Supplementary Figure S4). These SNPs were encoded with a Fluidigm chip for genotyping, and genotypes were generated from two other populations (Figure 1B) from different regions of Korea to include abalone diversity (Figure 1C). Detailed genotype summaries are provided in Supplementary Table S3. Finally, the 96 selected SNPs were subjected to the ML models to determine the predictive potential from. In this study, we used 7 ML methods (AB, BT, GBR, BLR, PLS, RF, and SVM-LK) with three combinations of genotyped datasets (Supplementary Table S4; Supplementary Figure S5, S6). We observed that an increase in the size of the dataset from different populations increased the ML prediction balanced accuracy (Table 1). Furthermore, the RF performed well in a pooled dataset (i.e., Type 3) with 0.714 balanced accuracy. Further, while optimizing the machine with 96 SNPs as features, we identified a subset of the SNPs (i.e., 38 SNPs) that contributed to the higher accuracy (Table 1; Supplementary Table S5). The features were selected with the same seven machines and the final features were selected from the seven machines with mean probabilities ≥0.1. This preliminary dataset could be a valuable asset to gain insight into heat resistance trait selection during abalone breeding. Detailed annotations of the SNPs are provided in Supplementary Table S5.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
The animal study was reviewed and approved by The studies involving animals were reviewed and approved by National Institute of Fisheries Science (Research Planning Division), 051-720-2871.
Author contributions
SS, SC, and YS: bioinformatic analysis; EN, YS, SS, and BN: manuscript preparation; EN, YK, CP, JL, and BN: sampling and sequencing; and YK and BN: funding and modelling.
Funding
This work was supported by the Collaborative Genome Program of the Korea Institute of Marine Science and Technology Promotion funded by the Ministry of Oceans and Fisheries (No. 20180430) and the National Institute of Fisheries Science, Ministry of Oceans and Fisheries, South Korea (R2022028).
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2023.1151427/full#supplementary-material
References
Adachi, K., and Okumura, S-I. (2012). Determination of genome size of Haliotis discus hannai and H. diversicolor aquatilis (Haliotidae) and phylogenetic examination of this family. Fish. Sci. 78, 849–852. doi:10.1007/s12562-012-0518-0
Bolger, A. M., Lohse, M., and Usadel, B. (2014). Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 30, 2114–2120. doi:10.1093/bioinformatics/btu170
Browning, B. L., Tian, X., Zhou, Y., and Browning, S. R. (2021). Fast two-stage phasing of large-scale sequence data. Am. J. Hum. Genet. 108, 1880–1890. doi:10.1016/j.ajhg.2021.08.005
Cingolani, P., Platts, A., Wang, Le L., Coon, M., Nguyen, T., Wang, L., et al. (2012). A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly. (Austin) 6, 80–92. doi:10.4161/fly.19695
Crosson, L. M., Lottsfeldt, N. S., Weavil-Abueg, M. E., and Friedman, C. S. (2020). Abalone withering syndrome disease dynamics: Infectious dose and temporal stability in seawater. J. Aquat. Anim. Health 32, 83–92. doi:10.1002/aah.10102
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, 2156–2158. doi:10.1093/bioinformatics/btr330
Frederick, A. R., Lee, A. M., Wehrle, B. A., Catabay, C. C., Rankins, D. R., Clements, K. D., et al. (2022). Abalone under moderate heat stress have elevated metabolic rates and changes to digestive enzyme activities. Comp. Biochem. Physiol. A Mol. Integr. Physiol. 270, 111230. doi:10.1016/j.cbpa.2022.111230
Im, J., and Kim, H. S. (2020). Genetic features of Haliotis discus hannai by infection of vibrio and virus. Genes Genomics 42, 117–125. doi:10.1007/s13258-019-00892-w
Kang, H. Y., Lee, Y. J., Song, W. Y., Kim, T. I., Lee, W. C., Kim, T. Y., et al. (2019). Physiological responses of the abalone Haliotis discus hannai to daily and seasonal temperature variations. Sci. Rep. 9, 8019. doi:10.1038/s41598-019-44526-3
Kim, C. H., Kim, E. J., Seo, C., Park, C. J., and Nam, Y. K. (2021). Transcriptome expression profiles between diploid and triploid Pacific abalone (Haliotis discus hannai) juveniles in response to acute heat-stress and hypoxia treatments. Mar. Genomics 57, 100820. doi:10.1016/j.margen.2020.100820
Kyeong, D., Kim, J., Shin, Y., Subramaniyam, S., Kang, B-C., Shin, E-H., et al. (2020). Expression of heat shock proteins in thermally challenged pacific abalone Haliotis discus hannai. Genes 11, 22. doi:10.3390/genes11010022
Langmead, B., and Salzberg, S. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359. doi:10.1038/nmeth.1923
Malik, A., Subramaniyam, S., Kim, C-B., and Manavalan, B. (2022). SortPred: The first machine learning based predictor to identify bacterial sortases and their classes using sequence-derived information. Comput. Struct. Biotechnol. J. 20, 165–174. doi:10.1016/j.csbj.2021.12.014
Mckenna, A., Hanna, M., Banks, E., Sivachenko, A., Cibulskis, K., Kernytsky, A., et al. (2010). The genome analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 20, 1297–1303. doi:10.1101/gr.107524.110
Nam, B-H., Jung, M., Subramaniyam, S., Yoo, S-I., Markkandan, K., Moon, J-Y., et al. (2016). Transcriptome analysis revealed changes of multiple genes involved in Haliotis discus hannai innate immunity during Vibrio parahemolyticus infection. PLOS ONE 11, e0153474. doi:10.1371/journal.pone.0153474
Nam, B-H., Kim, H., Seol, D., Kim, H., Noh, E. S., Kim, E. M., et al. (2021). Genotyping-by-Sequencing of the regional Pacific abalone (Haliotis discus) genomes reveals population structures and patterns of gene flow. PLOS ONE 16, e0247815. doi:10.1371/journal.pone.0247815
Nam, B. H., Kwak, W., Kim, Y. O., Kim, D. G., Kong, H. J., Kim, W. J., et al. (2017). Genome sequence of Pacific abalone (Haliotis discus hannai): The first draft genome in family haliotidae. Gigascience 6, 1–8. doi:10.1093/gigascience/gix014
Peter, A. C. (2016). Recent trends in worldwide abalone production. J. Shellfish Res. 35, 581–583. doi:10.2983/035.035.0302
Pritchard, J. K., Stephens, M., and Donnelly, P. (2000). Inference of population structure using multilocus genotype data. Genetics 155, 945–959. doi:10.1093/genetics/155.2.945
Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A., 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
Swezey, D. S., Boles, S. E., Aquilino, K. M., Stott, H. K., Bush, D., Whitehead, A., et al. (2020). Evolved differences in energy metabolism and growth dictate the impacts of ocean acidification on abalone aquaculture. Proc. Natl. Acad. Sci. U. S. A. 117, 26513–26519. doi:10.1073/pnas.2006910117
Tripp-Valdez, M. A., Harms, L., Pörtner, H. O., Sicard, M. T., and Lucassen, M. (2019). De novo transcriptome assembly and gene expression profile of thermally challenged green abalone (Haliotis fulgens: Gastropoda) under acute hypoxia and hypercapnia. Mar. Genomics 45, 48–56. doi:10.1016/j.margen.2019.01.007
Keywords: abalone, Haliotis discus hannai, heat resistant, genotype, aquaculture
Citation: Noh ES, Subramaniyam S, Cho S, Kim Y-O, Park C-J, Lee J-H, Nam B-H and Shin Y (2023) Genotyping of Haliotis discus hannai and machine learning models to predict the heat resistant phenotype based on genotype. Front. Genet. 14:1151427. doi: 10.3389/fgene.2023.1151427
Received: 26 January 2023; Accepted: 20 March 2023;
Published: 31 March 2023.
Edited by:
Jordan Poley, Scientist at the Center for Aquaculture Technologies in Canada, CanadaCopyright © 2023 Noh, Subramaniyam, Cho, Kim, Park, Lee, Nam and Shin. 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: Younhee Shin, yhshin@insilicogen.com;