- Research article
- Open Access
Construction of high-density genetic linkage maps for orange-spotted grouper Epinephelus coioidesusing multiplexed shotgun genotyping
BMC Genetics volume 14, Article number: 113 (2013)
Orange-spotted grouper, Epinephelus coioides, is one of the most valuable fish species in China. Commercial production of orange-spotted grouper could be increased by developing higher growth rates and improving commercially important traits. Information on genetic markers associated with quantitative trait loci (QTL) can be used in breeding programs to identify and select individuals carrying desired traits. A high-density genetic linkage map is the basis for QTL study, and multiplexed shotgun genotyping (MSG) facilitates the development of single nucleotide polymorphisms (SNPs) and genotyping. In this study, the first high-density genetic linkage maps for groupers were generated on the basis of the MSG method.
The sex-averaged map contained a total of 4,608 SNPs, which spanned 1581.7 cM, with a mean distance between SNPs of 0.34 cM. The 4,608 SNPs were located in 2,849 unique locations on the linkage map, with an average inter-location space at 0.56 cM. There were 2,516 SNPs on the female map, and the number of unique locus was 1,902. However, the male map contained more numbers of SNP (2,939) and unique locations (2,005). The total length of the female and male maps was 1,370.9 and 1,335.5 cM, respectively.
The high-resolution genetic linkage maps will be very useful for QTL analyses and marker-assisted selection (MAS) for economically important traits in molecular breeding of the orange-spotted grouper.
Grouper is the common name for a group of fishes belonging to a number of genera in the subfamily Epinephelinae of the family Serranidae, and they mainly inhabit in the Indian Ocean, the Red Sea, the Mediterranean Sea, the Northern Pacific Ocean, the Western Pacific Ocean and the Southeast Asia . The orange-spotted grouper, Epinephelus coioides, is a valuable and popular seafood fish, and is one of the major mariculture species in China. According to FAO fishery statistics, the global aquaculture production of the orange-spotted grouper was 152 t and the value had reached to 1,246,000 US dollars in 2011 . Currently in China, the orange-spotted grouper has become a major food in live fish markets and is an important cultured fish for commercial sale in Guangdong, Hainan and Fujian provinces of China. The vast market demands for orange-spotted grouper have driven efforts to breed families and populations with higher growth rates and lower food coefficients. The application of marker-assisted selection (MAS) to the orange-spotted grouper will become a promising strategy for improving growth traits .
A high-density genetic linkage map is the basis for quantitative trait loci (QTL) mapping, MAS, and functional gene positional cloning, and will be useful for functional genomics and genetic breeding studies. Genetic maps can also provide insights into genome organization and evolution through comparative studies. At present, high-density genetic linkage maps are already available for some important economic fishes, including Atlantic salmon , catfish , rainbow trout , tilapia , and Japanese flounder . Recently, a genetic linkage map of kelp grouper (Epinephelus bruneus) has been established based on microsatellite markers, which is the first linkage map in the subfamily Epinephelinae . However, construction of a high-density genetic linkage map for groupers has not been reported.
The density of genetic linkage map is mainly determined by the chosen genetic markers. Single nucleotide polymorphisms (SNPs) describes polymorphisms caused by point mutations at a given nucleotide position within a locus. SNPs are abundant and widely distributed throughout the genomes. As a result, SNPs have emerged as the most attractive genetic markers for construction of high-density genetic linkage map . With the advent of next-generation sequencing (NGS), there are several such approaches that are capable of discovering, sequencing and genotyping thousands of SNPs across almost any genome of interest in a single step, even in species in which little or even no genetic information is available . A particularly efficient protocol, termed as “multiplexed shotgun genotyping” (MSG), base on the NGS, has been used for high-throughput discovering of SNPs. In fact, MSG is similar in spirit to restriction site-associated DNA (RAD) sequencing  and whole-genome resequencing (WGD) , but combines the advantages of both techniques. MSG involves a highly simplified protocol for library preparation that requires only ~2 d of lab work to process 96 individuals or more. The technique is inexpensive because it requires only standard molecular laboratory equipment, uses only one set of bar-coded adapters, and does not require shearing and repairing of genomic DNAs. Because the approach does not depend on manual shearing, small amounts of DNA isolated from single individuals can be processed .
Here, we generated an orange-spotted grouper F1 full-sib family, which was genotyped by multiplexed shotgun sequencing, and at the first time reported the construction of high-density genetic linkage maps for groupers.
Family material and DNA isolation
An orange-spotted grouper F1 full-sib family was generated on the experimental base of Hainan Green Aquatic Breeding Co. Ltd., China, in May 2011. Fin clips of the parents were collected and kept in absolute ethanol, whereas the whole body of 142 offspring at the age of 35 days post hatch were soaked in absolute ethanol, and kept in a -20°C freezer. Genomic DNA was isolated using the standard phenol-chloroform protocol . DNA quality was evaluated by Qubit Fluorometer (Invitrogen, USA) and electrophoresis on a 0.6% agarose gel. All experiments were performed in accordance with the guidelines of the Animal Ethics Committee and were approved by the Institutional Review Board on Bioethics and Biosafety of BGI.
MSG library construction and sequencing
Genomic DNA of 144 individuals (142 offspring and 2 parents) were used to construct the library. We prepared the sequencing library using multiplexed shotgun sequencing method proposed by Andolfatto et al. , with some modifications. Genomic DNA of each sample (1 μg) was digested with 1 μl of FastDigest TaqI restriction enzyme (Thermo Scientific, USA) in FastDigest buffer for 10 minutes at 65°C in a total volume of 30 μl. Barcode adapters were designed and modified according to the standard Illumina adapters designed for paired-end read libraries. Unique barcode adapters (10 μmol) were added to each sample well. The ligation reaction was incubated for 1 hour at 22°C with 2 μl of T4 DNA ligase (Enzymatics, USA), 4 μl of 10× ligase buffer, 30 μl of digested products and 2.5 μl H2O. The T4 ligase was heat deactivated at 65°C for 20 minutes. Twenty-four ligation products of different samples were pooled in a single tube, and then 2 μl of chloroform was added to inactivate the restriction enzyme. DNA fragments were purified by excising a DNA fraction of 400–600 bp on a 2% agarose gel and then extracted by a QIA quick Gel Extraction Kit (Qiagen, Germany). All the products were amplified with 10 cycles of PCR in 50-μl reactions with 1 μl of common primer (10 μM), 1 μl of index primer and 25 μl of Phusion Master Mix (Finnzymes, Thermo Scientific, USA). The PCR profile was as follows: initial 30 s at 98°C, then 10 cycles, each with 30 s DNA denaturation at 98°C, 30 s at the annealing temperature 65°C and 30 s extension at 72°C, and a final extension of 5 min at 72°C. The PCR amplicons of 450–750 bases were size extracted using gel electrophoresis (0.5 × TAE, 2% agarose), subsequently purified using a QIAquick PCR Purification Kit (Qiagen), and quantified on Agilent 2100 Bioanalyzer (Agilent, USA). Sequencing was performed on the Illumina Hiseq 2000 platform (Illumina, USA) in 90-bp pair-end reads following standard protocols. Four lanes in the Illumina Hiseq 2000 were used.
Raw data filter and barcode reads split
Sequence reads from the Illumina runs were filtered as follows: reads with ambiguous base calls (N) more than ten percent, excessive low-quality positions (>40 positions with quality <2) were removed. The remaining trimmed, high-quality reads formed the basis for all subsequent analysis. Sequence reads from the same library were then sorted into individuals by barcode. Meanwhile reads with a barcode that did not match one of the expected barcodes (i.e. a sequencing error in the barcode) were discarded.
Retained reads were sorted into loci and genotyped using Stacks software (Version 0.9996) to analyze these maps cross data . The Stacks software package is freely downloadable at http://creskolab.uoregon.edu/stacks/. We used Stacks to identify loci in a set of individuals, de novo, and then genotype each locus. Stacks incorporates the likelihood-based SNP calling algorithm (maximum likelihood statistical model), which can evaluate each nucleotide position in every tag of all individuals, and then to identify sequence polymorphisms and distinguish them from sequencing errors [15, 16]. Some tag genotypes contained a single SNP, but others represented alleles that differed by multiple SNPs that were scored from these haplotypes.
Markers segregated in four different patterns. Type nn × np (1:1) was homozygous in the male and heterozygous in the female; lm × ll (segregating 1:1) was heterozygous in the male parent and homozygous in the female parent; hk × hk (1:2:1) was heterozygous in both parents with two shared alleles; and ef × eg (1:1:1:1) was heterozygous in both parents with two sex-specific alleles and one shared allele.
Selection and acquisition of high-quality SNP markers
As mentioned above, not every tag genotype contained a single SNP. In fact, some contained multiple SNPs. We selected SNPs for our further genetic map construction according to haplotypes identified by the Stacks genotypes program output using our custom scripts. The SNP markers were firstly filtered by excluded loci that genotype call rate less than 95%, and then a Chi-square goodness-of-fit test was used to assess the Mendelian segregation patterns. Those SNP markers showing significant segregation distortion (χ2 test, P < 0.001) were discarded.
Genetic linkage map construction
Linkage groups (LGs) were established by JoinMap v4.0 . The linkage map was built using the regression mapping algorithm, a recombination frequency smaller than 0.4, and an independence LOD threshold of 8.0. Recombination frequencies were converted to centiMorgans (cM) using Kosambi’s method for map-distance calculation, and the LGs maps were drawn and aligned using MapChart v2.2 .
Consensus sequences of the mapped MSG tag (84 bases in length) were aligned with the genomic sequences of five other fish species. The zebrafish Danio rerio (Zv9), three-spined stickleback Gasterosteus aculeatus (ver. 1.0), medaka Oryzias latipes (ver. 1.0), tilapia Oreochromis niloticus (ver. 1.0) and fugu Takifugu rubripes (ver. 5.0) genome sequences were downloaded, and blastn (BLAST + ver. 2.2.21)  searches with an e-value cutoff of 10-10 were conducted. In cases where the search of a query sequence hit two or more loci, a hit with the smallest e-value was considered significant. Significant hits on the chromosomes were used, including unoriented scaffolds assigned to chromosomes in the fugu genome. The Oxford grids  were constructed to study synteny and to compare positions of the homologous loci using Grid Map ver. 3.0a (http://cbr.jic.ac.uk/dicks/software/Grid_Map/).
Sequencing and genotyping
The Hiseq 2000 sequencing yielded 594,142,945 90-base reads. We generated more than 93 Gb of clean sequence data after excluded low quality raw reads. The sequences have been submitted to DDBJ under the accession number PRJDB1493. The average count of MSG tags per individual was 7,694,515. The mean coverage depth of MSG tags is 17.4 ×. MSG tags were aligned and clustered into 423,943 stacks, and 58,708 candidate MSG loci were inferred. For the analysis of the F1 mapping population, 25,892 SNPs were informative and were scored for sufficient numbers of progeny. Among them, 18,256 SNPs were retained after discarding those with a deviation from a Mendelian segregation pattern, and they were passed forward into the linkage map construction. The sequence dataset for this study was shown in the Additional file 1.
Genetic linkage map
Linkage analysis identified 24 linkage groups (LG1-LG24), which is consistent with the haploid chromosome number of the orange-spotted grouper . The sex-averaged map (Figure 1) contained a total of 4,608 SNPs, which spanned 1581.7 cM, with a mean distance between SNPs of 0.34 cM. The number of SNPs mapped on a linkage group of the sex-averaged map ranged from 115 (LG14) to 250 (LG5); the lengths of the linkage groups ranged from 43.8 (LG15) to 80.8 cM (LG4) (see details in Table 1). However, distribution of markers was uneven, likely reflecting differences in recombination frequency along the length of the grouper chromosomes. Clustered marker regions were observed in every linkage group of the sex-averaged map, especially in positions close to the centromeres and, less frequently, at the telomeres. To characterize marker clustering, we tabulated the number of markers sharing an identical genetic map location with other markers, as well as the number of unique location. The 4,608 SNPs were located in 2849 unique locations on the linkage map, with an average inter-location space of 0.56 cM (Table 1).
Sex specific maps (Additional file 2: Figure S1 and Additional file 3: Figure S2) were also constructed. There were 2,516 SNPs on the female map, and the number of unique locus was 1,902. The length of the female map was 1370.9 cM, with the average inter-SNP distance and inter-location space at 0.54 and 0.72 cM respectively. The male map contained 2,939 SNPs, and there were 2,005 unique locations. The length of the male map was 1335.5 cM, and the mean length between SNPs and locations was 0.45 and 0.67 cM respectively. The female map comprised linkage groups ranging in length from 45.3 to 72.4 cM while the male map contained linkage groups with a length ranging from 46.2 to 65.3 cM (Table 1).
The differences in recombination rates between the sexes are presented in Table 1. The overall recombination rate between female and male is 1.03, and the rates among different linkage groups range from 0.71 to 1.33.
Syntenies between different fish species
BLAST searches of the 5,194 mapped MSG tag consensus sequences from the orange-spotted grouper against the genome sequences of tilapia, stickleback, fugu, medaka and zebrafish indicated variation in the syntenic relationship between orange-spotted grouper and the respective species. Homology was most frequently inferred to the tilapia genome, with 225 tags being mapped to it. In contrast, the other four species yielded limited numbers of similarity hits. Only 177, 83, 72 and 15 tags mapped to the stickleback, fugu, medaka and zebrafish genome sequences, respectively. Owing to the tilapia genome sequence have not been integrated with linkage map, the syntenic relationship of linkage group between orange-spotted grouper and tilapia could not be evaluated. Detailed syntenic pairs between orange-spotted grouper linkage groups and four model fish species chromosomes were show on Additional file 4: Figure S3.
MSG and linkage mapping
In this study, 25,892 SNPs were informative and were scored for sufficient numbers of progeny, and then a Chi-square goodness-of-fit test (χ2 test, P < 0.001) was used to assess these SNPs whether deviate from the Mendelian segregation patterns. 18,256 SNPs were retained after discarding those with a deviation from a Mendelian segregation pattern, suggesting non-monogenetic inheritance of those discarded SNP markers. Technical artifacts may be responsible for the distorted markers, but biological processes known as transmission ratio distortion (TRD) also cause a deviation from Mendelian segregation . Both regression and ML mapping algorithm in JoinMap4.0 can be used to arrange markers. More SNPs can be mapped base on the ML mapping algorithm (data not shown). However, when we evaluated the results by aligning the SNPs that arranged on the maps to the genome of E. coioides (genome data not published), we found the order of SNPs based on the regression mapping algorithm was more reliable. So we choose the regression mapping algorithm. Finally, a high-density linkage map which contained a total of 4,608 SNPs and spanned 1581.7 cM, with a mean distance between SNPs of 0.34 cm was obtained. Such highly dense linkage maps contain rich information on the genomic structure of an organism and could be useful for studies involving comparative genomics and QTL mapping.
To date, microsatellite marker has been a popular option for linkage analyses in organisms without genomic information, especially for aquatic animals [23, 24]. Microsatellite markers are sequence-based, but they are costly and time-consuming if hundreds or thousands of markers are involved . In contrast to microsatellite markers, SNPs generated by MSG are abundant and highly suitable for cost effective high-throughput genotyping. These SNPs are sequence-based, allowing the practice of comparative genomics [22, 25], which aids in exploring candidate genes for traits of interest  and even assembling denovo genomic sequences . Moreover, allelic information on a large number of SNPs is readily available without prior curation and labourious experiments. The present study further demonstrates the utility of MSG in the genomic study of a non-model organism, yielding a wealth of genomic information without prior knowledge of the genome of an examined species.
The utility of genetic maps is correlated with the distribution of markers across linkage groups. Clustered markers in areas of minimal recombination, while allowing general linkage group assignment, often cannot be used for definitive fine mapping and positional cloning . In fish, clustering of DNA markers on genetic linkage maps has been observed in medaka , rainbow trout , channel catfish [5, 30], and Atlantic salmon . Although potential explanations for high levels of marker clustering are not completely understood [5, 32], large areas of repetitive DNA in teleost genome result from whole genome duplication could be the important reason for this phenomenon.
Differences in recombination rates between the male and female
It is common to find a difference in the recombination rate between the two sexes in any fish species with female map distances usually larger than those in male map. For instance, the female:male recombination rate was 1.37:1 in Atlantic salmon , 1.68:1 in rainbow trout , 1.43:1 in Japanese flounder  and 1.6:1 in catfish [5, 33]. In contrast to these previous reports, the linkage maps presented in this study (Table 1) show a smaller overall recombination rate between female and male (1.03:1). The reason for this remarkable difference in reported recombination rates is likely to be our maps improved marker coverage in telomeric regions. Since male recombination is often elevated in telomeres [34, 35], the more comprehensive coverage of these regions in our study has resulted in a more even recombination rate between sexes. Thus, it is quite likely that the male and female map lengths will converge further when marker density is increased.
Comparative genome analysis
The sequences of 225, 177, 83, 72 and 15 mapped markers had significant hits in the whole genome sequences of tilapia, stickleback, fugu, medaka and zebrafish, respectively, suggesting that tilapia is more closely related to grouper than the other four fishes.
Linkage map with sequence-based markers is a platform for comparative genome studies. It is well known that determination of gene functions is difficult in non-model species; functional genome analysis will have to rely heavily on the establishment of homologies from model species. Mapping more gene sequences on the linkage map of the orange-spotted grouper should enhance comparative mapping, thereby transferring genome information from model species to the orange-spotted grouper.
We constructed highly dense genetic linkage maps of the orange-spotted grouper using MSG. The number of linkage groups is consistent with the haploid chromosome number of the orange-spotted grouper.
This study presents the high density genetic linkage maps for the groupers. It was produced from genotypes of 142 F1 full-sib progeny and included 4608 SNP markers. The sex-averaged map revealed 24 linkage groups that covered 1581.7 cM, with an average intermarker distance of 0.34 cM. These high density maps allowed for comparison among them and those from model fish species, and will be useful for research on grouper genetics and molecular breeding.
Nelson JS: Fishes of the world. 1994, New York, NY, USA: John Wiley and Sons
FAO: Cultured Aquatic Species Information Programme, Epinephelus coioides. Available online: http://www.fao.org (accessed on 11 June 2013)
Wei Y, Huang H, Meng ZN, Zhang Y, Luo J, Chen GH, Lin HR: Single nucleotide polymorphisms in the leptin-a gene and associations with growth traits in the orange-spotted grouper (epinephelus coioides). Int J Mol Sci. 2013, 14: 8625-8637. 10.3390/ijms14048625.
Lien S, Gidskehaug L, Moen T, Hayes BJ, Berg PR, Davidson WS, Omholt SW, Kent MP: A dense SNP-based linkage map for Atlantic salmon (Salmo salar) reveals extended chromosome homeologies and striking differences in sex-specific recombination patterns. BMC Genomics. 2011, 12: 615-10.1186/1471-2164-12-615.
Ninwichian P, Peatman E, Liu H, Kucuktas H, Somridhivej B, Liu SK, Li P, Jiang YL, Sha ZX, Kaltenboeck L, Abernathy JW, Wang WQ, Chen F, Lee Y, Wong LL, Wang SL, Lu JG, Liu ZJ: Second-generation genetic linkage map of catfish and its integration with the BAC-based physical map. G3: Genes, Genomes, Genetics. 2012, 2: 1233-1241.
Rexroad CE, Palti Y, Gahr SA, Vallejo RL: A second generation genetic map for rainbow trout (Oncorhynchus mykiss). BMC Genet. 2008, 9: 74-10.1186/1471-2156-9-74.
Lee BY, Lee WJ, Streelman JT, Carleton KL, Howe AE, Hulata G, Slettan A, Stern JE, Terai Y, Kocher TD: A second-generation genetic linkage map of tilapia (Oreochromis spp.). Genetics. 2005, 170: 237-244. 10.1534/genetics.104.035022.
Castaño-Sánchez C, Fuji K, Ozaki A, Hasegawa O, Sakamoto T, Morishima K, Nakayama I, Fujiwara A, Masaoka T, Okamoto H, Hayashida K, Taganmi M, Kawai J, Hayashizaki Y, Okamoto N: A second generation genetic linkage map of Japanese flounder (Paralichthys olivaceus). BMC Genomics. 2010, 11: 554-10.1186/1471-2164-11-554.
Liu Q, Sakamoto T, Kubota S, Okamoto N, Yamashita H, Takagi M, Shigenobu Y, Sugaya T, Nakamura Y, Sano M, Wuthisuthimethavee S, Ozaki A: A genetic linkage map of kelp grouper (Epinephelus bruneus) based on microsatellite markers. Aquaculture. 2013, 414–415: 63-81.
Davey JW, Hohenlohe PA, Etter PD, Boone JQ, Catchen JM, Blaxter ML: Genome-wide genetic marker discovery and genotyping using next-generation sequencing. Nat Rev Genet. 2011, 12: 499-510. 10.1038/nrg3012.
Baird N, Etter P, Atwood T, Currey M, Shiver A, Lewis Z, Selker E, Cresko W, Johnson E: Rapid SNP discovery and genetic mapping using sequenced RAD markers. Plos One. 2008, 3 (10): e3376-10.1371/journal.pone.0003376.
Huang X, Feng Q, Qian Q, Zhao Q, Wang L, Wang A, Guan J, Fan D, Weng Q, Huang T, Dong GJ, Sang T, Han B: High-throughput genotyping by whole-genome resequencing. Genome Res. 2009, 19: 1068-1076. 10.1101/gr.089516.108.
Andolfatto P, Davison D, Erezyilmaz D, Hu TT, Mast J, Sunayama-Morita T, Stern DL: Multiplexed shotgun genotyping for rapid and efficient genetic mapping. Genome Res. 2011, 21: 610-617. 10.1101/gr.115402.110.
Sambrook J, Russell D: Molecular Cloning: A Laboratory Manual. 2011, Cold Spring Harbor, NY: Cold Spring Harbor Press, 3
Catchen JM, Amores A, Hohenlohe P, Cresko W, Postlethwait JH: Stacks: building and genotyping loci de novo from short-read sequences. G3: Genes, Genomes, Genetics. 2011, 1: 171-182.
Hohenlohe PA, Bassham S, Etter PD, Stiffler N, Johnson EA, Cresko WA: Population genomics of parallel adaptation in threespine stickleback using sequenced RAD tags. PloS Genetics. 2010, 6 (2): e1000862-10.1371/journal.pgen.1000862.
Van Ooijen JW: JoinMap® 4, Software for the calculation of genetic linkage maps in experimental populations. 2006, Wageningen: Kyazma BV
Voorrips RE: MapChart: software for the graphical presentation of linkage maps and QTLs. J Hered. 2002, 93 (1): 77-78. 10.1093/jhered/93.1.77.
Edwards JH: The Oxford gird. Ann Hum Genet. 1991, 55: 17-31. 10.1111/j.1469-1809.1991.tb00394.x.
Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, Madden T: BLAST+: architecture and applications. BMC Bioinforma. 2009, 10: 421-10.1186/1471-2105-10-421.
Ding SX, Wang SF, Wang DX, Wang J: Analysis of the karyotype of Epinephalus coioides. J Xiamen Univ (Nat Sci). 2004, 43 (3): 426-428.
Kakioka R, Kokita T, Kumada H, Watanabe K, Okuda N: A RAD-based linkage map and comparative genomics in the gudgeons (genus Gnathopogon, Cyprinidae). BMC Genomics. 2013, 14: 32-10.1186/1471-2164-14-32.
Wang CM, Bai ZY, He XP, Lin G, Xia JH, Sun F, Lo LC, Feng F, Zhu ZY, Yue GH: A high-resolution linkage map for comparative genome analysis and QTL fine mapping in Asian seabass, Lates calcarifer. BMC Genomics. 2011, 12: 174-10.1186/1471-2164-12-174.
Liu F, Sun F, Li J, Xia JH, Lin G, Tu RJ, Yue GH: A microsatellite-based linkage map of salt tolerant tilapia (Oreochromis mossambicus x Oreochromis spp.) and mapping of sex-determining loci. BMC Genomics. 2013, 14: 58-10.1186/1471-2164-14-58.
Amores A, Catchen J, Ferrara A, Fontenot Q, Postlethwait JH: Genome evolution and meiotic maps by massively parallel DNA sequencing: Spotted gar, an out group for the teleost genome duplication. Genetics. 2011, 188: 799-808. 10.1534/genetics.111.127324.
Yang HA, Tao Y, Zheng ZQ, Li CD, Sweetingham MW, Howieson JG: Application of next-generation sequencing for rapid marker development in molecular plant breeding: a case study on anthracnose disease resistance in Lupinus angustifolius L. BMC Genomics. 2012, 13: 318-10.1186/1471-2164-13-318.
Dasmahapatra KK, Walters JR, Briscoe AD, Davey JW, Whibley A, Nadeau NJ, Zimin AV, Hughes DS, Ferguson LC, The Heliconius Genome Consortium, et al: Butterfly genome reveals promiscuous exchange of mimicry adaptations among species. Nature. 2012, 487 (7405): 94-98.
Naruse K, Fukamachi S, Mitani H, Kondo M, Matsuoka T, Kondo S, Hanamura N, Morita Y, Hasegawa K, Nishigaki R, Shimada A, Wada H, Kusakabe T, Suzuki N, Kinoshata M, Kanamori A, Terado T, Kimura H, Nonoka M, Shima A: A detailed linkage map of medaka, Oryzias latipes: Comparative genomics and genome evolution. Genetics. 2000, 154: 1773-1784.
Robison BD, Wheeler PA, Sundin K, Sikka P, Thorgaard GH: Composite interval mapping reveals a major locus influencing embryonic development rate in rainbow trout (Oncorhynchus mykiss). J Hered. 2001, 92: 16-22. 10.1093/jhered/92.1.16.
Liu ZJ, Karsi A, Li P, Cao DF, Dunham R: An AFLP-based genetic linkage map of channel catfish (Ictalurus punctatus) constructed by using an interspecific hybrid resource family. Genetics. 2003, 165: 687-694.
Lorenz S, Brenna-Hansen S, Moen T, Roseth A, Davidson WS, Omholt SW, Lien S: BAC-based upgrading and physical integration of a genetic SNP map in Atlantic salmon. Anim Genet. 2010, 41: 48-54. 10.1111/j.1365-2052.2009.01963.x.
Watanabe T, Yoshida M, Nakajima M, Taniguchi N: Linkage mapping of AFLP and microsatellite DNA markers with the body color and sex-determining loci in the guppy (Poecilia reticulata). Zoolog Sci. 2005, 22: 883-889. 10.2108/zsj.22.883.
Kucuktas H, Wang SL, Li P, He CB C, Xu P, Sha ZX, Liu H, Jiang YL, Baoprasertkul P, Somridhivej B, Wang YP, Abernathy J, Guo XM, Liu L, Muir W, Liu ZJ: Construction of genetic linkage maps and comparative genome analysis of catfish using gene-associated markers. Genetics. 2009, 181: 1649-1660. 10.1534/genetics.108.098855.
Danzmann RG, Davidson EA, Ferguson MM, Gharbi K, Koop BF, Hoyheim B, Lien S, Lubieniecki KP, Moghadam HK, Park J, Phillips RB, Davidson WS: Distribution of ancestral proto-Actinopterygian chromosome arms within the genomes of 4R-derivative salmonid fishes (Rainbow trout and Atlantic salmon). BMC Genomics. 2008, 9: 557-10.1186/1471-2164-9-557.
Moen T, Hayes B, Baranski M, Berg PR, Kjoglum S, Koop BF, Davidson WS, Omholt SW, Lien S: A linkage map of the Atlantic salmon (Salmo salar) based on EST-derived SNP markers. BMC Genomics. 2008, 9: 223-10.1186/1471-2164-9-223.
This work was supported by China 863 project (No. 012AA10A407), Shenzhen Scientific R & D Grant (No. XB201108250095A), the National Natural Science Foundation (No.31370047) and Guangdong Marine Fisheries Promotion Program (No.A201100D01, A201100D03).
The authors declare that they have no competing interests.
QS, XY, YZ and HL initiated and supervised the study. LJ, GC, HZ and YZ prepared the F1 full-sib progeny of orange-spotted grouper. JC, JL, JB, QM, QX, QC and YC performed the experiments. LS, XY and SL conducted data analysis. XY, LS and QS prepared the manuscript. All authors read and approved the final manuscript.
Xinxin You, Liping Shu, Shuisheng Li contributed equally to this work.