+Advanced Search
Volume 2 Issue 2
May  2020
Article Contents


High-resolution linkage and quantitative trait locus mapping using an interspecifc cross between Argopecten irradians irradians (♀) and A. purpuratus (♂)

  • Edited by Xin Yu.
  • Junxia Mao and Qifan Zeng contributed equally to this work.
  • Electronic supplementary material The online version of this article (https://doi.org/10.1007/s42995-020-00029-z) contains supplementary material, which is available to authorized users.
  • The bay scallop and Peruvian scallop are economically important species. Interspecific hybrids of these two scallops outperformed both of their parent species in multiple growth traits but exhibited decreased fertility, which provides good models for the study of heterosis and species divergence. Genetic mapping serves as a chromosomal-level framework to investigate the molecular mechanisms of hybridization and introgression. In this study, high-resolution linkage maps were constructed for the bay and Peruvian scallops with an interspecific hybrid family. The linkage map of the bay scallop covered over 98.9% of the whole genome with 2994 mapped markers and the average marker interval of 0.32 cM. For the Peruvian scallop, 1585 markers were mapped with the average maker interval of 0.51 cM, covering 97.7% of the genome. Both the two linkage maps have 16 linkage groups, corresponding to the haploid chromosome number of the two species. Approximately, 54.5% of markers exhibited significant deviation from the expected Mendelian ratio of segregation, lending insights into the intrinsic incompatibility between the two species. QTLs related to growth and shell coloration were detected, which could explain 13.1% and 74.9% of the phenotypic variance, respectively. This represents important information for further evaluation. These findings are an important addition to the genomic resources for scallop genetic studies, and are especially useful for investigations on genomic incompatibility for hybridization, genome evolution of closely related species, and genetic enhancement programs in aquaculture.
  • 加载中
  • Adamkewicz L, Castagna M (1988) Genetics of shell color and pattern in the bay scallop Argopecten irradians. J Hered 79:14-17 doi: 10.1093/oxfordjournals.jhered.a110436
    Barfield S, Aglyamova GV, Matz MV (2016) Evolutionary origins of germline segregation in Metazoa: evidence for a germ stem cell lineage in the coral Orbicella faveolata (Cnidaria, Anthozoa). Proc R Soc Biol Sci 283:20152128 doi: 10.1098/rspb.2015.2128
    Bath E, Bowden S, Peters C, Reddy A, Tobias J, Easton-Calabria E, Seddon N, Goodwin S, Wigby S (2017) Sperm and sex peptide stimulate aggression in female Drosophila. Nat Ecol Evol 1:0154 doi: 10.1038/s41559-017-0154
    Blanc G, Charcosset A, Mangin B, Gallais A, Moreau L (2006) Connected populations for detecting quantitative trait loci and testing for epistasis: an application in maize. Theor Appl Genet 113:206-224 doi: 10.1007/s00122-006-0287-1
    Brand A (2006) Scallop ecology: distributions and behaviour. In: Shumway SE, Parsons GJ (eds) Scallops: biology, ecology, and aquaculture. Elsevier, Amsterdam, pp 651-713
    Brennan A, Hiscock S, Abbott R (2014) Interspecific crossing and genetic mapping reveal intrinsic genomic incompatibility between two Senecio species that form a hybrid zone on Mount Etna, Sicily. Heredity 113:195-204 doi: 10.1038/hdy.2014.14
    Broman K, Sen S (2009) A guide to QTL mapping with R/qtl. Springer, New York
    Broman K, Hao W, Sen S, Churchill G (2003) R/qtl: QTL mapping in experimental crosses. Bioinformatics 19:889-890 doi: 10.1093/bioinformatics/btg112
    Buonamici S, Chakraborty S, Senyuk V, Nucifora G (2003) The role of EVI1 in normal and leukemic cells. Blood Cells Mol Dis 31:206-212 doi: 10.1016/S1079-9796(03)00159-1
    Chakravarti A, Lasher L, Reefer J (1991) A maximum likelihood method for estimating genome length using genetic linkage data. Genetics 128:175-182 
    Chanley PE (1961) Inheritance of shell markings and growth in the hard clam Venus mercenaria. Proc Natl Shellfish Assoc 50:163-169
    Cui Z, Hui M, Liu Y, Song C, Li X, Li Y, Liu L, Shi G, Wang S, Li F, Zhang X, Liu C, Xiang J, Chu KH (2015) High density linkage mapping aided by transcriptomics documents ZW sex determination system in the Chinese mitten crab Eriocheir sinensis. Heredity 115:206-215 doi: 10.1038/hdy.2015.26
    Ding J, Zhao L, Chang Y, Zhao W, Du Z, Hao Z (2015) Transcriptome sequencing and characterization of Japanese scallop Patinopecten yessoensis from different shell color lines. PLoS ONE 10:e0116406 doi: 10.1371/journal.pone.0116406
    Dixon GB, Davies SW, Aglyamova GA, Meyer E, Bay LK, Matz MV (2015) Genomic determinants of coral heat tolerance across latitudes. Science 348:1460-1462 doi: 10.1126/science.1261224
    Dobzhansky T (1947) Genetics of natural populations. XIV. A response of certain gene arrangements in the third chromosome of Drosophila pseudoobscura to natural selection. Genetics 32:142-160 
    Douglas R, Par K (2003) Common features of segregation distortion in plants and animals. Genetica 117:27-35 doi: 10.1023/A:1022308414864
    Elshire R, Glaubitz J, Sun Q, Poland JA, Kawamoto K, Buckler ES, Mitchell SE (2011) A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS ONE 6:e19379 doi: 10.1371/journal.pone.0019379
    FAO STAT (2016) Food and Agriculture Organization of the United Nations. Statistical Yearbook, Rome
    Fishman L, Kelly AJ, Morgan E, Willis J (2001) A genetic map in the Mimulus guttatus species complex reveals transmission ratio distortion due to heterospecific interactions. Genetics 159:1701-1706 
    Fu X, Dou J, Mao J, Su H, Jiao W, Zhang L, Hu X, Huang X, Wang S, Bao Z (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.0079960
    Garcia A, Garcia S, Melchinger A, Zeng Z (2008) QTL mapping and the genetic basis of heterosis in maize and rice. Genetics 180(3): 1707-1724 doi: 10.1534/genetics.107.082867
    Gardner C, Lonnquist J (1959) Linkage and the degree of dominance of genes controlling quantitative characters in maize. Agron J 51:524-528 doi: 10.2134/agronj1959.00021962005100090005x
    Gary FN (1980) Genetics of shell color in Mytilus edulisland the association of growth rate with shell color. J Exp Mar Biol Ecol 47:89-94 doi: 10.1016/0022-0981(80)90140-9
    Guo Y, Yuan H, Fang D, Song L, Liu Y, Liu Y, Wu L, Yu J, Li Z, Xu X, Zhang H (2014) An improved 2b-RAD approach (I2b-RAD) offering genotyping tested by a rice (Oryza sativa L.) F2 population. BMC Genom 15:956 doi: 10.1186/1471-2164-15-956
    Haley C, Knott S (1992) A simple regression method for mapping quantitative trait loci in line crosses using flanking markers. Heredity 69:315-324 doi: 10.1038/hdy.1992.131
    Hu L, Huang X, Mao J, Wang C, Bao Z (2013) Genomic characterization of interspecific hybrids between the scallops Argopecten purpuratus and A. irradians irradians. PLoS ONE 8:e62432 doi: 10.1371/journal.pone.0062432
    Hu L, Huang X, Sun Y, Mao J, Wang S, Wang C, Bao Z (2015) Molecular genetic analysis of heterosis in interspecific hybrids of Argopecten purpuratus × A. irradians irradians. Genet Mol Res 14:10692-10704 doi: 10.4238/2015.September.9.9
    Jiao W, Fu X, Dou J, Li H, Su H, Mao J, Yu Q, Zhang L, Hu X, Huang X, Wang Y, Wang S, Bao Z (2013) 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
    Kaeppler SM (2012) Heterosis: many genes, many mechanisms-end the search for an undiscovered unifying theory. ISRN Bot 2012:682824 
    Kim B, Jang S, Chu S, Bordiya Y, Akter M, Lee J, Chin J, Koh H (2014) Analysis of segregation distortion and its relationship to hybrid barriers in rice. Rice 7:3 doi: 10.1186/s12284-014-0003-8
    Kosambi D (1944) The estimation of map distances from recombination values. Ann Eugen 12:172-175 
    Krieger U, Lippman Z, Zamir D (2010) The flowering gene SINGLE FLOWER TRUSS drives heterosis for yield in tomato. Nat Genet 42:459-463 doi: 10.1038/ng.550
    Kurasawa Y, Hu H, Zhou Q, Li Z (2018) A trypanosome-specific protein cooperates with the CIF1 protein to promote cytokinesis in Trypanosoma brucei. J Biol Chem 293:10275 doi: 10.1074/jbc.RA118.003113
    Kusterer B, Muminovic J, Utz HF, Piepho H, Barth S, Heckenberger M, Meyer R, Altman T, Melchinger AE (2007) Analysis of a triple testcross design with recombinant inbred lines reveals a significant role of epistasis in heterosis for biomass-related traits in Arabidopsis. Genetics 175:2009-2017 doi: 10.1534/genetics.106.069005
    Layer G, Reichelt J, Jahn D, Heinz DW (2010) Structure and function of enzymes in heme biosynthesis. Protein Sci 19:1137-1161 doi: 10.1002/pro.405
    Li L, Guo X (2004) AFLP-based genetic linkage maps of the Pacific oyster Crassostrea gigas Thunberg. Mar Biotechnol 6:26-36 doi: 10.1007/s10126-003-0001-0
    Li X, Schimenti J (2007) Mouse pachytene checkpoint 2 (Trip13) is required for completing meiotic recombination but not synapsis. PLoS Genet 3:e130 doi: 10.1371/journal.pgen.0030130
    Li L, Xiang J, Liu X, Zhang Y, Dong B, Zhang X (2005) Construction of AFLP-based genetic linkage map for Zhikong scallop, Chlamys farreri Jones et Preston and mapping of sex-linked markers. Aquaculture 245:63-73 doi: 10.1016/j.aquaculture.2004.12.015
    Li H, Liu X, Zhang G (2012) A consensus microsatellite-based linkage map for the hermaphroditic bay scallop (Argopecten irradians) and its application in size-related QTL analysis. PLoS ONE 7:e46926 doi: 10.1371/journal.pone.0046926
    Li Y, Sun X, Hu X, Xun X, Zhang J, Guo X, Jiao W, Zhang L, Liu W, Wang J, Li J, Sun Y, Miao Y, Zhang X, Cheng T, Xu G, Fu X, Wang Y, Yu X, Huang X et al (2017) Scallop genome reveals molecular adaptations to semi-sessile life and neurotoxins. Nat Commun 8:1721 doi: 10.1038/s41467-017-01927-0
    Li C, Liu X, Liu B, Ma B, Liu F, Liu G, Shi Q, Wang C (2018) Draft genome of the Peruvian scallop Argopecten purpuratus. GigaScience 7:giy031
    Lindberg DR, Pearse JS (1990) Experimental manipulation of shell color and morphology of the limpets Lottia asmi (Middendorff) and Lottia digitalis (Rathke) (Mollusca: Patellogastropoda). J Exp Mar Biol Ecol 140:173-185 doi: 10.1016/0022-0981(90)90125-V
    Liu X, Wu F, Zhao H, Zhang G, Guo X (2009) A novel shell color variant of the Pacific abalone Haliotis discus hannai Ino subject to genetic control and dietary influence. J Shellfish Res 28:419-424 doi: 10.2983/035.028.0226
    Liu X, Guo L, You J, Liu X, He Y, Yuan J, Liu G, Feng Z (2010) Progress of segregation distortion in genetic mapping of plants. Res J Agron 4:78-83 doi: 10.3923/rjagr.2010.78.83
    Lowry DB, Hernandez K, Taylor SH, Meyer E, Logan TL, Barry KW, Chapman JA, Rokhsar DS, Schmutz J, Juenger TE (2015) The genetics of divergence and reproductive isolation between ecotypes of Panicum hallii. New Phytol 205:402-414 doi: 10.1111/nph.13027
    Lu Z, Yang A, Wang Q, Liu Z, Zhou L (2006) Assortative fertilization in Chlamys farreri and Patinopecten yessoensis and its implication in scallop hybridization. J Shellfish Res 25:509-514 doi: 10.2983/0730-8000(2006)25[509:AFICFA]2.0.CO;2
    Lyttle T (1991) Segregation distorters. Annu Rev Genet 25:511-557 doi: 10.1146/annurev.ge.25.120191.002455
    Maniatis T, Fritsch E, Sambrook J (1989) Molecular cloning: a laboratory manual. Cold Spring Harbor, New York
    Manichaikul A, Moon J, Sen S, Yandell B (2009) A model selection approach for the identification of quantitative trait loci in experimental crosses, allowing epistasis. Genetics 181:1077-1086 doi: 10.1534/genetics.108.094565
    Marín A, Fujimoto T, Arai K (2013) Genetic structure of the Peruvian scallop Argopecten purpuratus inferred from mitochondrial and nuclear DNA variation. Mar Genom 31:1-8 
    Miura O, Nishi S, Chiba S (2007) Temperature-related diversity of shell colour in the intertidal gastropod Batillaria. J Mollus Stud 73:235-240 doi: 10.1093/mollus/eym019
    Moll R, Lindsey M, Robinson H (1964) Estimates of genetic variances and level of dominance in maize. Genetics 49:411-423 
    Narum S, Buerkle C, Davey J, Miller MR, Hohenlohe PA (2013) Genotyping-by-sequencing in ecological and conservation genomics. Mol Ecol 22:2841-2847 doi: 10.1111/mec.12350
    Newkirk GF (1980) Genetics of shell color in Mytilus edulis L. and the association of growth rate with shell color. J Exp Mar Biol Ecol 47:89-94 doi: 10.1016/0022-0981(80)90140-9
    Ogut F, Bian Y, Bradbury PJ, Holland JB (2015) Joint-multiple family linkage analysis predicts within-family variation better than single-family analysis of the maize nested association mapping population. Heredity 114:552-563 doi: 10.1038/hdy.2014.123
    Ortiz-Barrientos D, Reiland J, Hey J, Noor MAF (2002) Recombination and the divergence of hybridizing species. Genetica 116:167-178 doi: 10.1023/A:1021296829109
    Ostberg C, Hauser L, Pritchard V, Garza J, Naish K (2013) Chromosome rearrangements, recombination suppression, and limited segregation distortion in hybrids between Yellowstone cutthroat trout (Oncorhynchus clarkiibouvieri) and rainbow trout (O. mykiss). BMC Genom 14:570 doi: 10.1186/1471-2164-14-570
    Panithanarak T, Hauffe H, Dallas J, Glover A, Ward R, Searle J (2004) Linkage-dependent gene flow in a house mouse chromosomal hybrid zone. Evolution 58:184-192 doi: 10.1111/j.0014-3820.2004.tb01585.x
    Pauletto M, Carraro L, Babbucci M, Lucchini R, Bargelloni L, Cardazzo B (2016) Extending RAD tag analysis to microbial ecology: a comparison between multi locus sequence typing (MLST) and 2b-RAD to investigate Listeria monocytogenes genetic structure. Mol Ecol Resour 16:823-835 doi: 10.1111/1755-0998.12495
    Petersen J, Baerwald M, Ibarra A, May B (2012) A first-generation linkage map of the Pacific lion-paw scallop (Nodipecten subnodosus): initial evidence of QTL for size traits and markers linked to orange shell color. Aquaculture 350:200-209
    Phadnis N, Orr H (2009) A single gene causes both male sterility and segregation distortion in Drosophila hybrids. Science 323:376-379 doi: 10.1126/science.1163934
    Poland J, Rife T (2012) Genotyping-by-sequencing for plant breeding and genetics. Plant Genome 5:92-102 doi: 10.3835/plantgenome2012.05.0005
    Qin Y, Liu X, Zhang H, Zhang G, Guo X (2007a) Identification and mapping of amplified fragment length polymorphism markers linked to shell color in bay scallop, Argopecten irradians irradians (Lamarck, 1819). Mar Biotechnol 9:66-73 doi: 10.1007/s10126-006-6076-7
    Qin Y, Liu X, Zhang H, Zhang G, Guo X (2007b) Genetic mapping of size-related quantitative trait loci (QTL) in the bay scallop (Argopecten irradians) using AFLP and microsatellite markers. Aquaculture 272:281-290 doi: 10.1016/j.aquaculture.2007.07.214
    Sastry AN (1963) Reproduction of the bay scallop, Aequipecten irradians Lamarck. Influence of temperature on maturation and spawning. Biol Bull 125:146-153
    Schwartz D, Laughner W (1969) A molecular basis for heterosis. Science 166:626-627 doi: 10.1126/science.166.3905.626
    Seetharam AS, Stuart GW (2013) Whole genome phylogeny for 21 Drosophila species using predicted 2b-RAD fragments. PeerJ 1:e226 doi: 10.7717/peerj.226
    Semel Y, Nissenbaum J, Menda N, Zinder M, Krieger U, Issman N, Pleban T, Lippman Z, Gur A, Zamir D (2006) Overdominant quantitative trait loci for yield and fitness in tomato. Proc Natl Acad Sci USA 103:12981-12986 doi: 10.1073/pnas.0604635103
    Shang L, Ma L, Wang Y, Su Y, Wang X, Li Y, Abduweli A, Cai S, Liu F, Wang K, Hua J (2016) Main effect QTL with dominance determines heterosis for dynamic plant height in upland cotton. G3(Bethesda) 3:3373-3379 
    Stam P (1993) Construction of integrated genetic linkage maps by means of a new computer package: Join Map. Plant J 3:739-744 doi: 10.1111/j.1365-313X.1993.00739.x
    Streiff R, Courtois B, Meusnier S, Bourguet D (2014) Genetic mapping of two components of reproductive isolation between two sibling species of moths Ostrinia nubilalis and O. scapulalis. Heredity 112:370-381 doi: 10.1038/hdy.2013.113
    Sturtevant AH (1917) Genetic factors affecting the strength of linkage in Drosophila. Proc Natl Acad Sci 3:555-558 doi: 10.1073/pnas.3.9.555
    Tian M, Li Y, Jing J, Mu C, Du H, Dou J, Mao J, Li X, Jiao W, Wang Y, Hu X, Wang S, Wang R, Bao Z (2015) Construction of a high-density genetic map and quantitative trait locus mapping in the sea cucumber Apostichopus japonicus. Sci Rep 5:14852 doi: 10.1038/srep14852
    Venkateswarlu K, Hanada T, Chishti A (2005) Centaurin-α1 interacts directly with kinesin motor protein KIF13B. J Cell Sci 118:2471-2484 doi: 10.1242/jcs.02369
    Voorrips R (2001) MapChart: software for the graphical presentation of linkage maps and QTLs. J Hered 93:77-78 
    Wada KT, Komaru A (1990) Inheritance of white coloration of the prismatic layer of shells in the Japanese pearl oyster Pincada fucata martensii and its importance in the pearl culture industry. Nippon Suisan Gakkaishi 56:1787-1790 doi: 10.2331/suisan.56.1787
    Wang Y, Guo X (2004) Chromosomal rearrangement in Pectinidae revealed by rRNA loci and implications for bivalve evolution. Biol Bull 207:247-256 doi: 10.2307/1543213
    Wang S, Guo X (2004) AFLP linkage map of an intraspecific cross in Chlamys farreri. J Shellfish Res 23:491-500
    Wang L, Song L, Chang Y, Xu W, Ni D, Guo X (2005) A preliminary genetic map of Zhikong scallop (Chlamys farreri Jones et Preston 1904). Aquacult Res 36:643-653 doi: 10.1111/j.1365-2109.2005.01268.x
    Wang L, Song L, Zhang H, Gao Q, Guo X (2007) Genetic linkage map of bay scallop, Argopecten irradians irradians (Lamarck 1819). Aquacult Res 38:409-419 doi: 10.1111/j.1365-2109.2007.01683.x
    Wang S, Zhang L, Hu J, Bao Z, Liu Z (2010) Molecular and cellular evidence for biased mitotic gene conversion in hybrid scallop. BMC Evol Biol 10:6 doi: 10.1186/1471-2148-10-6
    Wang C, Liu B, Li J, Liu S, Li J, Hu L, Fan X, Du H, Fang H (2011) Introduction of the Peruvian scallop and its hybridization with the bay scallop in China. Aquaculture 310:380-387 doi: 10.1016/j.aquaculture.2010.11.014
    Wang S, Meyer E, McKay JK, Matz MV (2012) 2b-RAD: a simple and flexible method for genome-wide genotyping. Nat Methods 9:808-810 doi: 10.1038/nmeth.2023
    Wang S, Liu P, Lv J, Li Y, Cheng T, Zhang L, Xia Y, Sun H, Hu X, Bao Z (2016) Serial sequencing of isolength RAD tags for cost-efficient genome-wide profiling of genetic and epigenetic variations. Nat Protoc 11:2189-2200 doi: 10.1038/nprot.2016.133
    Wang S, Zhang J, Jiao W, Li J, Xun X, Sun Y, Guo X, Huan P, Dong B, Zhang L, Hu X, Sun X, Wang J, Zhao C, Wang Y, Wang D, Huang X, Wang R, Lv J, Li Y et al (2017) Scallop genome provides insights into evolution of bilaterian karyotype and development. Nat Ecol Evol 1:0120 doi: 10.1038/s41559-017-0120
    Wijsman J, Troost K, Fang J, Roncarati A (2019) Global production of marine bivalves. Trends and challenges. In: Smaal A, Ferreira J, Grant J, Petersen J, Strand Ø (eds) Goods and services of marine bivalves. Springer, Cham, pp 7-26
    Williams ST (2017) Molluscan shell colour. Biol Rev 92:1039-1058 doi: 10.1111/brv.12268
    Williams ST, Lockyer AE, Dyal P, Nakano T, Churchill CK, Speiser DI (2017) Colorful seashells: identification of haem pathway genes associated with the synthesis of porphyrin shell color in marine snails. Ecol Evol 7:10379-10397 doi: 10.1002/ece3.3552
    Winkler F, Estévez B, Jollán L, Garrido J (2001) Inheritance of the general shell color in the Scallop Argopecten purpuratus (Bivalvia: Pectinidae). J Hered 92:521-525 doi: 10.1093/jhered/92.6.521
    Wolf D, Hallauer R (1997) Triple testcross analysis to detect epistasis in maize. Crop Sci 37:763-770 doi: 10.2135/cropsci1997.0011183X003700030012x
    Wolff M, Garrido J (1991) Comparative study of growth and survival of two color morphs of the Chilean scallop Argopecten purpuratus Lamarck (1819) in suspended culture. J Shellfish Res 10:47-53
    Wolff M, Mendo J (2000) Management of the Peruvian bay scallop (Argopecten purpuratus) metapopulation with regard to environmental change. Aquat Conserv 10:117-126 doi: 10.1002/(SICI)1099-0755(200003/04)10:2<117::AID-AQC399>3.0.CO;2-T
    Xiao J, Li J, Yuan L, Tanksley S (1995) Dominance is the major genetic basis of heterosis in rice as revealed by QTL analysis using molecular markers. Genetics 140:745-754 
    Yapici N, Kim Y, Ribeiro C, Dickson B (2008) A receptor that mediates the post-mating switch in Drosophila reproductive behaviour. Nature 451:33-37 doi: 10.1038/nature06483
    Yu S, Li J, Xu C, Tan Y, Gao Y, Li X, Zhang Q, Maroof M (1997) Importance of epistasis as the genetic basis of heterosis in an elite rice hybrid. Proc Natl Acad Sci 94:9226-9231 doi: 10.1073/pnas.94.17.9226
    Zhan A, Hu J, Hu X, Hui M, Wang M, Peng W, Huang X, Wang S, Lu W, Sun C, Bao Z (2009) Construction of microsatellite-based linkage maps and identification of size-related quantitative trait loci for Zhikong scallop (Chlamys farreri). Anim Genet 40:821-831 doi: 10.1111/j.1365-2052.2009.01920.x
    Zhang S, Li L, Wu F, Zhang G (2014) Yield trait improvement of bay scallops following complete diallel crosses between different scallop stocks. Chin J Oceanol Limnol 32:1-7 doi: 10.1007/s00343-014-3079-5
    Zhao L, Li Y, Li Y, Yu J, Liao H, Wang S, Lv J, Liang J, Huang X, Bao Z (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
    Zheng H, Zhang G, Liu X (2005) Comparison of growth and survival of larvae among different shell color stocks of bay scallop Argopecten irradians irradians (Lamarck 1819). Chin J Oceanol Limnol 23:183-188 doi: 10.1007/BF02894236
  • 加载中
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Figures(4) / Tables(4)

Article Metrics

Article views(599) PDF downloads(5) Cited by()

Proportional views

High-resolution linkage and quantitative trait locus mapping using an interspecifc cross between Argopecten irradians irradians (♀) and A. purpuratus (♂)

    Corresponding author: Chunde Wang, chundewang2007@163.com
    Corresponding author: Shi Wang, swang@ouc.edu.cn
  • 1. MOE Key Laboratory of Marine Genetics and Breeding, College of Marine Life Sciences, Ocean University of China, Qingdao 266003, China
  • 2. Key Laboratory of Mariculture and Stock Enhancement in North China's Sea, Ministry of Agriculture and Rural Afairs, Dalian Ocean University, Dalian 116023, China
  • 3. Laboratory for Marine Biology and Biotechnology, Pilot National Laboratory for Marine Science and Technology, Qingdao 266237, China
  • 4. Laboratory for Marine Fisheries Science and Food Production Processes, Pilot National Laboratory for Marine Science and Technology, Qingdao 266237, China
  • 5. Marine Science and Engineering College, Qingdao Agricultural University, Qingdao 266109, China

Abstract: The bay scallop and Peruvian scallop are economically important species. Interspecific hybrids of these two scallops outperformed both of their parent species in multiple growth traits but exhibited decreased fertility, which provides good models for the study of heterosis and species divergence. Genetic mapping serves as a chromosomal-level framework to investigate the molecular mechanisms of hybridization and introgression. In this study, high-resolution linkage maps were constructed for the bay and Peruvian scallops with an interspecific hybrid family. The linkage map of the bay scallop covered over 98.9% of the whole genome with 2994 mapped markers and the average marker interval of 0.32 cM. For the Peruvian scallop, 1585 markers were mapped with the average maker interval of 0.51 cM, covering 97.7% of the genome. Both the two linkage maps have 16 linkage groups, corresponding to the haploid chromosome number of the two species. Approximately, 54.5% of markers exhibited significant deviation from the expected Mendelian ratio of segregation, lending insights into the intrinsic incompatibility between the two species. QTLs related to growth and shell coloration were detected, which could explain 13.1% and 74.9% of the phenotypic variance, respectively. This represents important information for further evaluation. These findings are an important addition to the genomic resources for scallop genetic studies, and are especially useful for investigations on genomic incompatibility for hybridization, genome evolution of closely related species, and genetic enhancement programs in aquaculture.


  • Scallops represent a conspicuous group of bivalve molluscs, comprising over 300 extant species worldwide (Brand 2006). They are valuable models for studying early animal evolution and bivalve adaptation as they possess remarkable preservation of ancestral bilaterian linkage groups, and novel genomic features for adaptations to benthic life as semi-sessile filter feeders (Li et al. 2017; Wang et al. 2017). Also, scallops are popular aquaculture species and provide significant economic benefits. The production of maricultural bivalves is more than 13 million tonnes per year from 2010 to 2015 (Wijsman et al. 2019). The bay scallop, Argopecten irradians irradians, and the Peruvian scallop, A. purpuratus, are economically important species naturally distributed along the Atlantic coasts from Nova Scotia to the Gulf of Mexico and Pacific coasts from Peru to Chile, respectively (Wolff and Mendo 2000). Geographic isolation and climatic variability promoted their divergence and allopatric speciation (Brand 2006; Marín et al. 2013). Despite that, Argopecten scallops have distinct size, life span, temperature and hypoxia tolerance, these two species share similar life history and karyotypes (Hu et al. 2013, 2015; Wang and Guo 2004). Interspecific hybridization between the two species exhibited extraordinary increase in growth rates in both shell height and whole-body weight in F1 offspring (Wang et al. 2011). The exploitation of this heterosis could greatly benefit the scallop industry (Hu et al. 2013).

    Genetic linkage mapping is an efficient method to investigate the molecular mechanisms of hybridization and introgression between species. On the basis of genetic mapping, suppression of gene flow could be revealed by limited recombination between rearranged chromosomes (Ostberg et al. 2013; Panithanarak et al. 2004). Selections against particular hybrid genotypes result in serious deviations of observed genotypic frequencies from the expected Mendelian inheritance (Fishman et al. 2001; Kim et al. 2014). Therefore, segregation distortion markers have been extensively studied with interspecific populations to reveal the intrinsic incompatibilities (Brennan et al. 2014; Phadnis and Orr 2009; Streiff et al 2014).

    Moreover, genetic mapping serves as a chromosomal-level framework to detect quantitative trait loci (QTL) including that for heterosis and trace inheritance of various performance traits (Garcia et al. 2008; Shang et al. 2016). The dominant, over-dominant, epistatic and other effects for heterosis have been tested in plants and multiple livestock and aquaculture species (Gardner and Lonnquist 1959; Kaeppler 2012; Krieger et al. 2010; Kusterer et al. 2007; Moll et al. 1964; Schwartz and Laughner 1969; Semel et al. 2006; Wolf and Hallauer 1997; Xiao et al.1995; Yu et al. 1997). Limited by the number of developed markers and high-throughput genotyping platforms, genetic studies in scallops lack progress in uncovering the molecular mechanisms underlying various commercial traits. Most of the previously constructed genetic maps in scallops were with relatively low resolutions, limiting the sensitivity and accuracy in determining quantitative trait genes (Li et al. 2005, 2012; Wang et al. 2004, 2005, 2007; Zhan et al. 2009). Recent development of genotyping by sequencing (GBS) technologies has allowed high-throughput and efficient genotyping (Elshire et al. 2011; Narum et al. 2013; Poland and Rife 2012). Of these techniques, an optimized representative restriction site-associated DNA (RAD) method, 2b-RAD, has gained popularity owing to its simplicity in library preparation and ability to provide even and tunable genome coverage with a low requirement of sample quality (Wang et al. 2012, 2016). This method has been successfully applied for genetic and breeding studies, such as high-density linkage map construction and QTL mapping, in a wide range of species (Barfield et al. 2016; Cui et al. 2015; Dixon et al. 2015; Guo et al. 2014; Jiao et al. 2013; Lowry et al. 2015; Pauletto et al. 2016; Seetharam and Stuart 2013; Tian et al. 2015).

    In the present study, we conducted genetic linkage analysis using the 2b-RAD genotyping method based on a F1 interspecific family with the objectives of (1) construction of high-density genetic maps for the two scallops, (2) identification of segregation distortion loci to provide insights into genomic incompatibilities of the two species, and (3) genome scan of QTLs for growth-related traits to dissect the genetic architecture of interspecific heterosis.


    2b-RAD sequencing and genotyping of the mapping family

  • The sequencing data of the mapping family including two parents and 150 progenies (Fig. 1) are summarized in Supplementary Table S1. A total of 656 million reads were generated for the whole family, and 86.7% of the total reads passed quality trimming. For the parents, standard libraries were constructed and sequenced to a sufficient depth. A total of 33.1 and 32.3 million high-quality reads were generated for the parental bay scallop and Peruvian scallop, corresponding to the sequencing depth of 128.2 × and 122.8 ×, separately. For the hybrid progenies, RR libraries were used to reduce the total sequencing cost. Approximately, 3.4 million high-quality reads for each progeny were generated with an average sequencing depth of 22.3 ×, which was sufficient for accurate genotyping (Fu et al. 2013).

    Figure 1.  The interspecific cross between Argopecten irradians irradians (♀) and A. purpuratus (♂), and two shell color types of the hybrid progenies. Type Ⅰ: The left valves were orange and the right were purplish orange. Type Ⅱ: The left valves were purplish white with different degree coverage of dark blotches and the right were purplish white

    Clustering of the high-quality reads of the two parents generated 199, 962 representative reference sites, including 86, 544 parent-shared (~ 43.3%) and 113, 418 parent-specific (~ 56.4%) sites. After filtering sites with insufficient depth or from repetitive genomic regions, 73, 027 parent-shared and 111, 696 parent-specific sites were retained. For the parent-specific sites, 71, 692 (~ 64.2%) sites were from the bay scallop, and the remaining 40, 004 (~ 35.8%) sites were from the Peruvian scallop. Sequencing reads were then remapped to the constructed reference for genotyping. In total, 7653 polymorphic markers were identified in at least 80% progenies, which consisted of 1696 codominant markers and the rest 5957 dominant markers. The genotyping results for each segregation type are summarized in Table 1. The χ2 test was carried out to examine segregation distortion of the codominant markers and identified 924 markers that significantly deviated from the expected Mendelian ratio (P < 0.05), including 508 bay scallop-specific SNPs (Aa × aa type), 363 Peruvian scallop-specific SNPs (aa × Aa type), and 53 shared SNPs (Aa × Aa type) (Supplementary Figs. 1, 2). Notably, 93% of the shared codominant markers were identified with non-Mendelian inheritance.

    Marker type Segregation type Marker number Distorted marker number Distorted ratio
    Codominant Aa × aa 1018 508 49.9%
    aa × Aa 621 363 58.5%
    Aa × Aa 57 53 93.0%
    Sum 1696 924 54.5%
    Dominant A- × – 4047
    – × A- 1910
    Sum 5957
    Sum 7653

    Table 1.  Genotyping results and marker segregation analysis of the hybrid family

    Figure 2.  The high-resolution genetic linkage map for (a) A. irradians irradians containing 2994 markers in 16 linkage groups, and (b) A. purpuratus containing 1585 markers in 16 linkage groups. The white spaces are the gaps without markers mapped and black lines represent the positions of markers in the linkage maps

  • High-resolution linkage mapping for the bay and Peruvian scallops

  • In total, 6729 makers in accordance with Mendelian segregation (P ≥ 0.05) were used for the linkage mapping. These markers were partitioned into bay scallop heterozygous (Aa × aa and Aa × Aa type; 4561 markers) and Peruvian scallop heterozygous (aa × Aa and Aa × Aa type; 2172 markers) for constructing species-specific linkage maps. Using a LOD threshold of 10.0, markers were assigned into 16 linkage groups (LGs), which was consistent with the number of chromosomes of Argopecten haploid genome (Hu et al. 2013).

    The linkage map of the bay scallop contained 2994 markers and spanned 945.95 cM (Table 2, Fig. 2a). The average marker interval was 0.32 cM, ranging from 0.18 cM in LG1 to 0.63 cM in LG16. The expected map length was estimated to be 957.55 cM (Ge1) and 959.04 cM (Ge2) by two different methods. Genome coverage calculated using the average of the two estimates (958.30 cM) was as high as 98.9%. For the Peruvian scallop map, 1585 markers were mapped and spanned over 800.46 cM (Table 3, Fig. 2b). The average marker interval was 0.51 cM, ranging from 0.29 cM in LG1 to 1.35 cM in LG14. The constructed genome was estimated to be 819.01 cM in length by averaging Ge1 (816.78 cM) and Ge2 (821.24 cM) and covering about 97.7% of the whole genome. A few marker intervals were notably large in the genetic maps, which may represent recombination hotspots. Due to the large genetic differences between the bay and Peruvian scallops, insufficient shared markers were eligible to integrate the two maps. Therefore, consensus maps of the two species were not constructed.

    Linkage group Marker numbers Length (cM) Average marker interval (cM) ± standard deviation
    1 293 51.46 0.18 ± 0.73
    2 255 49.85 0.20 ± 0.91
    3 255 65.86 0.26 ± 1.97
    4 219 55.05 0.25 ± 1.34
    5 212 72.85 0.35 ± 2.13
    6 211 65.10 0.31 ± 1.49
    7 209 75.05 0.36 ± 1.42
    8 190 61.69 0.33 ± 1.54
    9 184 49.72 0.27 ± 1.71
    10 184 56.86 0.31 ± 1.79
    11 181 66.04 0.37 ± 1.77
    12 138 63.10 0.46 ± 1.43
    13 134 55.27 0.42 ± 1.53
    14 120 50.04 0.42 ± 1.50
    15 109 47.30 0.44 ± 1.56
    16 100 62.72 0.63 ± 2.92
    All 2994 947.95 0.32 ± 1.60

    Table 2.  Summary of the linkage map for A. irradians irradians

    Linkage group Marker numbers Length (cM) Average marker interval (cM)
    1 176 51.14 0.29 ± 0.80
    2 156 60.16 0.39 ± 1.19
    3 152 58.72 0.39 ± 1.31
    4 136 51.42 0.38 ± 1.24
    5 136 60.70 0.45 ± 1.70
    6 125 53.96 0.44 ± 1.11
    7 113 46.74 0.42 ± 1.33
    8 109 57.23 0.53 ± 1.69
    9 91 56.02 0.62 ± 1.89
    10 84 52.16 0.63 ± 1.68
    11 84 55.98 0.67 ± 2.30
    12 74 25.64 0.35 ± 0.92
    13 63 65.12 1.05 ± 2.35
    14 31 40.45 1.35 ± 2.48
    15 30 38.80 1.34 ± 3.04
    16 25 26.24 1.09 ± 2.03
    All 1585 800.46 0.51 ± 1.57

    Table 3.  Summary of the linkage map for A. purpuratus

  • QTL mapping for growth-related traits and shell color

  • The Kolmogorov–Smirnov test indicated that all the five growth-related traits followed normal distribution. The Pearson's correlation coefficients among the five growth traits ranged from 0.69 (shell width with shell weight) to 0.93 (shell height with shell length), indicating that they were highly correlated and possibly controlled by the same set of QTLs (Table 4, Supplementary Fig. S3). Therefore, we conducted a principle component analysis on the five traits. The first principle component (PC1), which explained over 80% of the total variation, was used as a composite trait for QTL mapping (Fig. 3). The full stepwise multiple QTL model identified a sole growth-related QTL (qgt) in LG2 (1.3–1.9 cM) of the Peruvian scallop, generated a LOD score of 4.7 and explained 13.1% of the variation in PC1 of the five traits (Fig. 4a). Evidence for epistatic effects has been found also among markers on LG1, LG2, and LG8, despite that the interactions do not meet the model selection significance threshold (Fig. 4b). Candidate genes, including SPR (sex peptide receptor), Trip13 (pachytene checkpoint protein 2 homolog), and KIF13B (kinesin-like protein KIF13B), were found to be located within the 95% credible interval of the QTL detected on LG2. These genes could be of causal importance that is worthy of further investigation. Notably, no LOD peaks above the threshold were detected in the bay scallop, which may be due to the non-polymorphism of candidate QTLs or associated markers in the bay scallop (Fig. 4c). Genomic scan of LOD profiles for five individual traits were also performed and revealed similar results (Supplementary Figs. S4, S5).

    SH 1.00 0.93 0.70 0.86 0.76
    SL 1.00 0.77 0.89 0.79
    SW 1.00 0.79 0.69
    BW 1.00 0.83
    AMW 1.00

    Table 4.  Pearson's correlation coefficients of the five growth traits

    Figure 3.  Principal component analysis of five growth traits in the hybrid scallops. a The PC1 explained about 85% of the total variance. b SW shell width, BW body weight, AMW adductor muscle weight, SL shell length, SH shell height

    Figure 4.  QTL mapping of growth traits and shell color in the hybrid scallops. The PC1 for five growth traits was used in QTL mapping with the linkage map of Peruvian scallop (a) and bay scallop (c). b LOD scores from a two-dimensional scan for growth traits. LODi measures the improvement in the fit of the full model over that of the additive model. It is displayed in the upper left triangle. LODf measures the improvement in the fit of the full two-locus model over the null model, indicating the evidence for at least one QTL, with allowance for interaction. It is displayed in the lower right triangle. QTL mapping for shell color with the linkage map of bay scallop (d) and Peruvian scallop (e). The solid and dashed lines indicate the LOD significance thresholds

    The two shell color types were checked to be consistent with the segregated ratio of 1:1 (P ≤ 0.05). Despite that, a single significant QTL (qsc) on LG2 in the bay scallop was identified with a LOD value of 45.1, no significant epistatic interactions were detected by the stepwise QTL selection model (Fig. 4d). The Bayesian 95% credible interval of qsc was only 0.83 cM, ranging from 13.2 to 14.1 cM. It explained ~ 74.9% of the variation in shell color. By mapping the markers to the bay scallop genome assembly, potential causal genes were identified, including PECR (peroxisomal trans-2-enoyl-CoA reductase), MECOM (MDS1 and EVI1 complex locus), and RTTN (Rotatin), while for the Peruvian scallop, no significant QTL signal was identified (Fig. 4e).

  • Scallops are important fishery and aquaculture species, with a production of over 2 million metric tonnes in 2016, accounting for 13% of the global mollusc production (FAO STAT 2016). Much of the research effort has been on evaluation and selection of crossbreeding and interspecific hybridization (Hu et al. 2015; Lu et al. 2006; Wang et al. 2010, 2011; Zhang et al. 2014). Understanding the genetic bases underlying genomic incompatibility and heterosis has been an important task to facilitate scallop breeding and genetic improvement. In this study, we performed genetic mapping with an interspecific hybrid resource family. Taking advantage of 2b-RAD method, high-resolution genetic linkage maps were constructed for the bay and Peruvian scallops by cost-effectively genotyping a full-sibling family of 152 scallops. The mapping density of the bay scallop was 0.32 cM with a 97.7% coverage, surpassing the resolution of all the previously constructed genetic maps (Li et al. 2012; Qin et al. 2007a, b; Wang et al. 2007). The linkage map for the Peruvian scallops constructed in this study is the first one for this species, with a resolution of 0.51 cM and a mapping coverage of 98.9%. Large differences in the polymorphic sites have been observed between the two species, with ~ 60.5% markers were species-specific polymorphic. A few marker intervals (i.e. 25.9 cM in LG16) were notably large in the current linkage map, which may represent hotspots of recombination. Therefore, more markers in these regions could be used for further study. Only ~ 3.4% codominant markers were shared between the two species, resulting in insufficient markers for construction of consensus maps.

    Markers deviating significantly from Mendelian inheritance were excluded from linkage mapping to avoid miscalculation. Despite that the level of segregation distortion varies greatly among species, ~ 54.5% markers showing segregation were identified in the hybrid progenies, which is much higher than the general segregation distortion ratio in marine bivalves, such as the Pacific oyster (Li and Guo 2004), the Zhikong scallop (Zhan et al. 2009), and the bay scallop (Li et al. 2012). In this study, over 93.0% of the SNPs shared by the two parents (i.e., Aa × Aa) deviated from the Mendelian segregation. Consistent with previous studies, the distorted markers tend to locate adjacently to form segregation distortion regions (SDRs), implying that inter-chromosomal rearrangements may occur within these regions (Douglas and Par 2003; Liu et al. 2010; Lyttle 1991). Chromosomal rearrangements may influence inheritance patterns by suppressing recombination in heterokaryotypes (Sturtevant 1917). Studies have reported that rearrangements, such as inversions, represent sets of coadapted alleles in a population that live in a constant environment. When two distinct populations are connected by migration, recombination loci involved in local adaptation regions are typically selected against (Dobzhansky 1947; Ortiz-Barrientos et al. 2002).

    Elucidation of the genetic basis underlying performance traits of scallop has been a central task for genetic improvement. QTLs for growth- and reproduction-related traits have been detected in the bay scallop (Li et al. 2012), the Zhikong scallop (Jiao et al. 2013), and the Pacific lion-paw scallop (Petersen et al. 2012). As revealed by Wang et al. (2011), the hybrids exhibited higher growth rates in both shell height and whole-body weight than the parental scallops, with a mid-parental heterosis of 125.9%–138.9% for whole-body weight and 145.4%–156.2% for adductor muscle weight. In this study, growth-related traits were highly correlated (Table 4), suggesting that they were possibly regulated by a same set of QTLs. Therefore, a principle component analysis was conducted to generate a composite trait for QTL mapping. By applying the newly constructed high-resolution linkage maps, a major QTL for the composite growth traits was mapped to LG2 of the Peruvian scallop. Despite that the gaps in genomic assembly may limit identification of the potential causal genes within the QTL, growth-related genes, including SPR, Trip13, and KIF13B were recognized within the loose window of this QTL (Supplementary Table S2). These genes are involved in the control of reproductive behaviors (Bath et al. 2017; Yapici et al. 2008), gamete generation (Li and Schimenti 2007), and cytokinesis (Kurasawa et al. 2018; Venkateswarlu et al. 2005). The list of candidate genes and SNPs associated with traits supplies potential targets in the future selective breeding programs to improve scallop aquaculture production, but more validation work is still needed in a broader population.

    In molluscs, shell coloration plays important roles in crypsis, mimicry and aposematism against visual predators including crabs, fish, and birds (reviewed in Williams 2017). Variations in shell colors and pigmentation patterns also associate with differences in growth, survival, and biological adaptations (Gary 1980; Wolff and Garrido 1991; Zheng et al. 2005). The mechanisms of shell coloration are complicated; they are affected both by environmental and genetic factors. Environmental effects including lights, climate, salinity and diets on shell color have been reported in previous studies (Lindberg and Pearse 1990; Liu et al. 2009; Miura et al. 2007). These modifications could be stable and long lasting once formed (Liu et al. 2009). Also, studies have revealed that shell colors are subject to strict genetic control. The general colorations and banding patterns of periostracum layer, raylike marks, and prismatic layer are controlled by single locus variations in the blue mussel (Newkirk 1980), hard clam (Chanley 1961), bay scallop (Adamkewicz and Castagna 1988), pearl oyster (Wada and Komaru 1990), and Yesso scallops (Zhao et al. 2017). Other coloration characteristics, such as pigmentation patterns, appear to be controlled by the epistatic effects of multiple nonallelic genes (Adamkewicz and Castagna 1988; Winkler et al. 2001). In this study, the hybrid scallop shells exhibited distinct background colors with a segregation ration of 1:1. A significant QTL signal on LG2 was detected in the bay scallop linkage map, explaining over 74.9% of the phenotypic variation. The markers located at the confidence intervals of this QTL constitute a group of valuable markers for further evaluation of their utility in marker-assisted selection. Candidate causal genes, and in particular, the MECOM gene, were identified within the QTL and focused our interest (Supplementary Table S3). MECOM is a transcriptional factor that is involved in regulation of myeloid differentiation of bone marrow progenitor cells to erythroid cells and granulocytes (Buonamici et al. 2003). Heme biosynthesis is induced during erythroid differentiation and is the main pathway for the formation of uroporphyrin pigments (Layer et al. 2010). Studies have revealed that uroporphyrin pigments are responsible for coloration of shells and soft tissues in several molluscs (reviewed in Williams 2017). In scallops and marine snails, several genes related to porphyrin metabolism were differentially expressed in mantles with different colors (Ding et al. 2015; Williams et al. 2017), whereas for MECOM, the difference is not significant in adult mantle tissues suggesting that it may express and function during early shell formation. In this study, the QTLs related to shell coloration and growth were revealed by mapping with one bi-parental family, which represents only a small cross section of the breeding germplasm. The expression of QTL depends on its presence and could be influenced by the genetic background (Blanc et al. 2006; Ogut et al. 2015). Therefore, further work would be needed to evaluate the QTL positions and effects with independent populations.

    In conclusion, we generated the first high-resolution linkage maps for A. purpuratus and A. irradians irradians. The marker density and mapping coverage were comparable to the highest level of linkage maps in molluscs allowing for comparative genome analysis using both genetic and genomic approaches. QTLs and candidate genes related to growth and shell coloration were detected, indicating important information for further evaluation. These findings are an important addition to the genomic resources for scallop genetic studies and are especially useful for investigations on genomic incompatibility for hybridization, genome evolution of closely related species, and genetic enhancement programs in aquaculture.

Materials and methods

    Mapping family

  • The mapping family was constructed using bay scallops from a farmed strain in Jiaonan, Qingdao, China, and Peruvian scallops from the progenies of wild stocks in Peru. As both species are hermaphrodite with external fertilization, inter-specific hybridization of A. irradians irradians (♀) and A. purpuratus (♂) was performed following the procedures previously described by Wang et al. (2011) to avoid self-fertilization. Briefly, sexually mature scallops with gonadal condition at stage Ⅳ (Sastry 1963) were collected and subjected to an air exposure of 30 min followed by a temperature shock from 21 to 23 ℃. Spawning scallops were separated in 5 L capacity polyethylene buckets containing 3 L of 30 μm filtered seawater to prevent contamination. The collected sperm were filtered through a 25 μm mesh screen to mitigate mixed eggs. The eggs were checked under a microscope 20 min post-spawning to verify that they were not self-fertilized. The unfertilized eggs of one bay scallop were mixed with the sperm of one Peruvian scallop for insemination. The mapping family was reared following a previous study (Wang et al. 2011). Briefly, larvae were raised at 21–23 ℃ with aeration, and fed with microalgae, i.e. Isochrysis galbana, at concentrations of 20, 000–100, 000 cells/ml depending on their stages. The density of larvae in the tank was maintained at 10 larvae/ml by adjusting water volume. Collectors were placed into the tanks as substrates for settlement when the eyespots became apparent. After 10 days of metamorphosis, the settled spat were transferred to the field and subsequently sorted into grow-out cages in the open sea of Jiaonan.

  • 2b-RAD sequencing of the mapping family

  • After six-months post fertilization, a total of 150 hybrid progenies were randomly taken for genotyping and subsequent analysis. The adductor muscles from the progenies and their parents were collected and stored in 70% ethanol for DNA extraction. Genomic DNA samples were prepared following standard phenol–chloroform DNA extraction procedure (Maniatis et al. 1989). 2b-RAD sequencing libraries were constructed using the protocol developed by Wang et al. (2012). The adaptors and PCR primers were modified to fit the Illumina HiSeq sequencing platforms (Supplementary Table S4). For the parents, standard sequencing libraries were constructed using adaptors with 5′-NNN-3′ overhangs to target all the BsaXI fragments, whereas for the hybrid progenies, reduced representation (RR) libraries were constructed by adaptors with 5′-NNA-3′ and 5′-NNT-3′ overhangs to capture a subset of all the fragments. All the libraries were labeled with unique barcodes and pooled together for single end sequencing (1 × 36 bp) utilizing a HiSeq-2000 system. All the sequencing data are deposited at the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) platform with accession no. PRJNA488174.

  • Data processing and de novo genotyping

  • RADtyping (version 1.0) was used to process the sequencing data and perform de novo dominant and codominant genotyping (Fu et al. 2013). Briefly, raw reads were first trimmed to remove adaptor sequences, low-quality reads (> 5 positions with quality < 10), uninformative reads (with no restriction sites, or containing ambiguous base calls), or reads with excessive homopolymer repeats (> 10 bp). Filtered reads of both parents were assembled into the exactly matching "allele" clusters, and then merged into "locus" clusters allowing two mismatches, which finally compromised the reference sites. Clusters with insufficient coverage (< 8 reads) or excessive depth (> 3000 reads) were excluded to avoid the interference of sequencing errors and repetitive elements. These clusters were further classified as parent-shared or parent-specific sites (sites present in one parent and absent in another due to the mutation of a recognition site) for subsequent dominant and codominant genotyping, respectively. Finally, filtered reads of the 150 progenies were mapped to the constructed reference for genotyping.

  • Genetic linkage map construction for the bay and Peruvian scallops

  • Markers were excluded from linkage analysis if they were not polymorphic in any of the parents or had low quality of genotype calls in more than 20% progenies. Mendelian segregation was examined through χ2 goodness-of-fit test, markers with significant segregation distortion were also filtered out (P < 0.05). Segregation-distorted regions were identified by aligning these distorted markers to the genome assemblies of the bay scallop (BioProject No. PRJNA428789) and the Peruvian scallop (Li et al. 2018) using Blastn (BLAST 2.2.28 +). Filtered markers were then used to construct the linkage map using JoinMap (version 4.0) (Stam 1993). Markers were categorized into three groups according to their phase types (lm × ll and nn × np as species-specific markers, and hk × hk as species-shared markers). Only species-specific markers were used for map construction. Markers were assigned to different linkage groups with logarithm of odds (LOD) cut-off of 10. The recombination frequencies were converted into map distances (centimorgans) through the Kosambi mapping function (Kosambi 1944). The expected genetic map lengths (Ge) were estimated according to Chakravarti et al. (1991) with the following formulas: Ge1 = Goa + 2S (Goa: the observed genetic map length; S: average marker interval); \(G_{{{\text{e}}1}} = \sum\nolimits_{i}^{n} {L_{i} \times (m_{i} + 1)/(m_{i} - 1)}\) (Li: the observed length of LGi; mi: the marker number of LGi; n: the number of LGs). Map coverage was measured as the ratio of observed genetic map length and the average of expected lengths. Finally, the genetic linkage map was graphically presented through MapChart (version 2.2) (Voorrips 2001).

  • QTL mapping for growth-related traits and shell color

  • The distributions of the five growth-related traits, including body weight (BW), adductor muscle weight (AMW), shell length (SL), shell height (SH), and shell width (SW), were checked whether they followed the normal distributions by the Kolmogorov–Smirnov test using MATLAB. The phenotypic correlations among the five growth traits were estimated by the Pearson's correlation coefficients. The progenies were characterized with two distinct types of shell coloration patterns (Fig. 1). The shell color was analyzed as a binary trait, and its segregation ratio was tested by χ2 test.

    QTL mapping was conducted for both growth-related traits and shell color with R/qtl (Broman et al. 2003) using the following procedures. The best multiple QTL model was determined using Haley–Knott regression (Haley and Knott 1992). Subsequently, automated stepwise model selection procedures were used to scan additive and epistatic QTL at each step (Manichaikul et al. 2009). Significant QTLs were identified with the LOD thresholds (experiment-wise α = 0.05) determined by 1000 permutations. Bayesian 95% credible intervals were calculated for each QTL as a measure of location uncertainties (Broman and Sen 2009). Also, single QTL and pairwise marker interactions were checked to confirm the selected model. The markers were mapped to the genome assemblies of the bay scallop (BioProject No. PRJNA428789) and the Peruvian scallop (Li et al. 2018) using Blastn (BLAST 2.2.28 +), genes located within the QTLs were identified as candidates. The full linkage maps of the two species, the sequences of markers, the phenotypic data, and the codes for QTL mapping have been deposited in https://www.figshare.com with the Digital Object Identifier (DOI) https://doi.org/10.6084/m9.figshare.8246867.

  • We acknowledge the Grant support from National Natural Science Foundation of China (U1706203, 31172404 and 31572618), Taishan Scholar Project Fund of Shandong Province of China, and Youth Talent Program Supported by Laboratory for Marine Fisheries Science and Food Production Processes, Pilot National Laboratory for Marine Science and Technology (Qingdao) (2018-MFS-T07).

Author contributions
  • SW and CW: conceived and designed the study; JM, ZY, HP, and LY: prepared the materials and constructed the sequencing libraries; JM and QZ: conducted genetic mapping and other data analysis; ZB: advised and coordinated the study. All the authors contributed to manuscript writing, reviewing, and approved the fnal version for submission.

Compliance and ethical standards

    Conflict of interest

  • The authors declare that they have no confict of interest.

  • Animal and human rights statement

  • We declare that all applicable international, national, and/or institutional guidelines for sampling, care, and experimental use of organisms for the study have been followed and all necessary approvals have been obtained.

Reference (99)



DownLoad:  Full-Size Img  PowerPoint