Skip to main content
  • Research article
  • Open access
  • Published:

Construction of a high-coverage bacterial artificial chromosome library and comprehensive genetic linkage map of yellowtail Seriola quinqueradiata



Japanese amberjack/yellowtail (Seriola quinqueradiata) is a commonly cultured marine fish in Japan. For cost effective fish production, a breeding program that increases commercially important traits is one of the major solutions. In selective breeding, information of genetic markers is useful and sufficient to identify individuals carrying advantageous traits but if the aim is to determine the genetic basis of the trait, large insert genomic DNA libraries are essential. In this study, toward prospective understanding of genetic basis of several economically important traits, we constructed a high-coverage bacterial artificial chromosome (BAC) library, obtained sequences from the BAC-end, and constructed comprehensive female and male linkage maps of yellowtail using Simple Sequence Repeat (SSR) markers developed from the BAC-end sequences and a yellowtail genomic library.


The total insert length of the BAC library we constructed here was estimated to be approximately 11 Gb and hence 16-times larger than the yellowtail genome. Sequencing of the BAC-ends showed a low fraction of repetitive sequences comparable to that in Tetraodon and fugu. A total of 837 SSR markers developed here were distributed among 24 linkage groups spanning 1,026.70 and 1,057.83 cM with an average interval of 4.96 and 4.32 cM in female and male map respectively without any segregation distortion. Oxford grids suggested conserved synteny between yellowtail and stickleback.


In addition to characteristics of yellowtail genome such as low repetitive sequences and conserved synteny with stickleback, our genomic and genetic resources constructed and revealed here will be powerful tools for the yellowtail breeding program and also for studies regarding the genetic basis of traits.


Species of yellowtail (family Carangidae) are widely distributed in the world’s ocean and are major target species for fisheries and aquaculture. The Japanese amberjack/yellowtail (Seriola quinqueradiata) is one of the most popular fish for consumption in Japan, where about 150,000 tons of farmed fish are produced each year. Although there is a huge market demand, seeds of this fish mostly rely on wild catch and hence artificial seed production is required for stable cultivation and breeding as well as reducing the negative effects of large-scale sampling of seed fish on natural stock. It is well known that using cultured brood fish for seed production reduces the environmental impact and allows the selection of commercially important traits and in such a case, marker-assisted selection (MAS) breeding based on studies regarding quantitative trait locus (QTL) is powerful and cost effective choice. Indeed, QTL studies have been performed in several fishes so far to improve production and life-history traits such as disease resistance and enhance growth rate [1]. To enable the QTL studies, linkage maps are required. In yellowtail, although a female linkage map has been constructed with 180 microsatellite markers [2, 3], the number of markers is not sufficient for fine QTL mapping and/or MAS in yellowtail breeding programs. Therefore, a higher-density linkage map is still required.

To isolate simple sequence repeats (SSRs) such as microsatellites and to further investigate the genetic basis of the traits, genomic information is essential. The sequences were isolated from genomic library and of the genomic library, one using bacterial artificial chromosome (BAC) system, called the BAC library, has been frequently used such as to generate whole-genome physical maps by DNA fingerprinting [4], to develop sequence-tagged connectors [5], and to sequence the genome itself [6] because of their insert size capacity, reproducibility and stability as the DNA sample [7]. By integrating BAC clones into linkage maps using BAC-derived sequences such as BAC-end sequences (BESs), BAC library also play important roles in genetic studies and subsequent positional cloning [8]. BAC libraries have been developed in several domestic animals, e.g. cattle [9], pig [10] and sheep [11], and in fishes, salmon [12], catfish [13], rainbow trout, carp, tilapia [14, 15], European sea bass [16, 17] and barramundi [18] but in yellowtails a BAC library has not yet been constructed.

In this study, to advance yellowtail genomic and genetic resources and for understanding of the genetic basis of several traits, we constructed a high-coverage BAC library, obtained BESs for preliminary survey of the genomic content and constructed comprehensive genetic linkage map of yellowtail.

Results and discussion

BAC library construction and BAC end sequencing

The yellowtail genomic DNA content, represented as C-value, was estimated to be 0.7 pg/cell (data not shown) using flow cytometric analysis and hence the genome size was calculated to be approximately 685 Mb. Of 100 randomly selected BAC clones, 71 (71%) contained inserts, indicating that approximately 78,520 (71% of 110,592 clones) clones had an insert. The size distribution of the 71 clones with inserts was from about 20 kb to 220 kb and average insert length was 140.7 kb (data not shown). Therefore, it is estimated that the total length of the yellowtail BAC library insert DNA was approximately 11 Gb and was 16-times larger than the yellowtail genome. It is known that a minimum of 5-10 × coverage across the entire genome is required for a BAC library to be useful for positional cloning, physical mapping, and genome sequencing [19]. Therefore, the yellowtail BAC library is sufficient for further genomic/genetic analysis except for studies regarding W-linked genes because of our ZZ male derived DNA source [20].

By sequencing both ends of randomly-selected 2,960 BAC clones, a total of 5,920 raw reads were obtained, and of those reads, 4,956 reads (2,471 in the SP6 side and 2,485 in the T7 side) were qualified for subsequent repeat identification and BLAST search (GA867436 - GA872391). Total length of the qualified BESs was 3,074,133 bp with an average size of 620 bp, representing approximately 0.45% of the yellowtail genome. The GC content was estimated to be 41.36%, which is almost the same as other fishes (Takifugu rubripes: 45.46%; Gasterosteus aculeatus: 44.60%; Oryzias latipes: 40.46%; Tetraodon nigroviridis: 46.43%) (

Preliminary survey of genomic content using BESs

Repeat content

The repeat elements were searched and screened from the qualified BESs. Total 211,184 bp (6.87%) of the qualified BESs are assigned to the repeat elements, of which 60,380 bp (1.96%), 40,242 bp (1.31%), 21,563 bp (0.70%) and 76,925 bp (2.50%) were classified as retroelements, DNA transposons, small RNA and simple repeats, respectively (Table 1). Assuming that the BAC library represents the genome of target specie, the abundance of the repetitive sequence in yellowtail genome is lower than the majority of teleost fishes studied so far such as rainbow trout (59.5%) [21], common carp (17.3%) [15], channel catfish (11.9%) [22] and Nile tilapia (14.0%) [23] and comparable to that in Tetraodon (6.2%) and in fugu (4.3%) [24].

Table 1 Repeat content of the yellowtail BESs

A total of 1,845 simple sequence repeats (SSRs) were identified from the BESs (Table 1). Of the SSRs, di-nucleotide repeats, particularly AC/GT repeats including CA/TG repeats, were the most abundant (Table 2).

Table 2 SSR distribution in the yellowtail BESs

Homology to other teleost genomes

To identify the homology between yellowtail and other fishes, the yellowtail qualified BESs were subjected to BLASTx and BLASTn searches against eight teleost proteomes and genomes respectively. The highest number of top hits, highest average bit score and % identity were observed in yellowtail-Nile tilapia in both BLAST results (Table 3). Total length of the queries in BLASTx hits between yellowtail and Nile tilapia was estimated to be 198,090 bp indicating that 6.4% of the qualifed BESs was protein coding sequence. The high sequence similarity between yellowtail and Nile tilapia can be explained by their phylogenic positions where they are both assigned in the order Perciformes [25]. In the BLASTn result, the second-highest number of top hit was observed in the yellowtail-stickleback comparison (Table 3). The high sequence similarity between stickleback and species in Perciformes such as striped bass and gilthead seabream has been reported and therefore our data is consistent with the previous observations [26, 27].

Table 3 Summary of BLAST searches of the yellowtail qualified BESs against eight fish genomes and proteomes

Genetic linkage map

Out of the 743 primer pairs designed from the qualified BESs, 373 primer pairs (27 mononucleotide repeats, 285 dinucleotide repeats, 31 trinucleotide repeats, 26 tetranucleotide repeats, 3 pentanucleotide repeats and 1 hexanucleotide repeat) produced amplicons. In addition to the 464 microsatellite markers derived from the genomic library Ohara et al. developed [2, 3], 837 markers in total were included in the yellowtail genetic linkage maps (Additional file 1). No segregation distortion was observed in any markers and hence lethal allele-linked markers were not included in our marker set.

Resultant yellowtail female and male genetic map consists of 715 and 702 markers including 232 and 271 framework markers, spanning 1,026.65 and 1,057.83 cM Kosambi with an average interval 4.96 and 4.32 cM on 24 linkage groups respectively (Table 4, Figure 1). The number of chromosomes in yellowtail has been reported to be 2n = 48 and hence the SSR markers we developed are distributed throughout the yellowtail genome [28]. The “gaps” observed in Squ21 and 24 in male and both map respectively might be caused by “recombination hot-spots” where recombination occurs frequently (Figure 1). The genome length was estimated to be 1,274.64 (L 1 ) and 1,284.34 (L 2 ) cM in the female and 1,282.35 (L 1 ) and 1,285.45 (L 2 ) cM in the male map by the two different methods respectively (see Materials and Methods). Using formula c = 1 – e-2dn/L and estimated genome length L, coverage of the female and male map is estimated to be 83.3 to 83.9% respectively (Table 4). Considering the average interval less than 10 cM and the genome coverage, we concluded that the yellowtail genetic map was sufficient for further QTL studies [29].

Table 4 Summary of the yellowtail genetic map
Figure 1
figure 1

Yellowtail female (left) and male (right) maps for linkage groups Squ1- Squ24. Total lengths of linkage groups are expressed in Kosambi cM. BES-derived SSR markers are coded “BAC” after a number, and microsatellite markers developed from genomic library are coded “TUF”.

Identification orthologous chromosomes with other fishes

In addition to the BAC or whole genome sequence, comparative genome analysis especially conserved synteny would be helpful for fine-scale QTL analyses and/or understanding the genetic basis of the traits [30, 31]. BLAST searches of the 818 mapped yellowtail loci against medaka, Tetraodon, stickleback, fugu and zebrafish proved that 25.7, 23.0, 42.2, 24.4 and 9.4% of the loci were mapped to each genome sequence. Oxford grids showed that eighteen linkage group pairs between yellowtail and stickleback retained a one-to-one relationship, and another three stickleback and six yellowtail linkage groups had a one-to-two relationship, implying that chromosomal fusions or breakages occurred after divergence from ancestor of both species (Figure 2). Nevertheless, the result suggests conserved synteny between yellowtail and stickleback and hence the stickleback genome data would be useful as a reference of yellowtail genome.

Figure 2
figure 2

Oxford grids between yellowtail and five model fish genomes. Numbers in boxes indicate the number of orthologous gene pairs. Boxes containing more than ten, seven and five orthologous gene pairs are highlighted in red, yellow and blue respectively.


We herein constructed a high-coverage BAC library and comprehensive genetic linkage map including BES-derived SSR markers of yellowtail (Seriola quinqueradiata). A survey of BESs showed a low frequency of repetitive sequences as much as that of Tetraodon and fugu. BLAST searches and Oxford grids against five fish genomes clearly showed conservation between yellowtail and stickleback genome. Generally, a high repetitive sequence frequency hampers chromosome walking and makes the positional cloning difficult [32]. A low frequency of repetitive sequences and relatively small genome size suggest that yellowtail would be an ideal species to study the genetic basis of economically important traits. In addition, conserved genome architecture with stickleback would be helpful for synteny-based identification of new genetic markers and genes in the target genomic segments. We have already started studies regarding several traits such as sex determination and disease resistance [20, 33]. We anticipate that the genomic and genetic resources we constructed will be powerful tools for further studies of these traits.


Ethics statements

Field permits are not required for this species in Japan. Since all fish treatments were performed in Goto Branch of Seikai National Fisheries Research Institute of Fisheries Research Agency, fish handling, husbandry and sampling methods were approved by Institutional Animal Care and Use Committee of National Research Institute of Aquaculture (IACUC-NRIA No. 03).

BAC construction and BAC-end sequencing

The BAC library was constructed according to Katagiri et al. with some modifications [14]. Briefly, at first, approximately 5 × 107 frozen sperm cells taken from one male yellowtail were embedded in agarose plugs, digested with proteinase K overnight at 37°C and stored in 0.5 M EDTA following proteinase K inhibitor treatment until use. The plugs were dialyzed in 0.5 × TE, partially digested with Mbo I and size fractionated by pulse-field electrophoresis. The fraction containing 150 to 250 kb genomic DNA was excised from the gel and was recovered as high molecular weight (HMW) genomic DNA. The HMW genomic DNA was then integrated into Bam HI site of pBACe3.6 vector and reactions were transfected to E. coli DH10B strain. Finally, a total of 110,592 recombinant BAC clones were picked and stored in 288 384-well microtiter plates. The length of the insert DNAs was estimated by analyzing 100 BAC inserts digested with Not I.

The BESs were obtained from eight 384-well plates containing 3,072 clones. The BAC DNAs extracted by conventional alkaline lysis method were sequenced from SP6 and T7 sides with BigDye Terminator v3.1 Cycle Sequencing Kit (Life Technologies) following the manufacturer’s instructions and reactions were electrophoresed with Applied Biosystems 3730 DNA Analyzer (Life Technologies). All raw reads were processed using PHRED software with default parameters except for the trimming error probability was set at 0.01 [34, 35], and vector and bacterial sequences were masked by CROSS_MATCH implemented in PHRAP software. The masked BESs of more than 100 bp in length, hereafter called “qualified BES”, were extracted using our in-house perl script. The GC content of the extracted BESs was estimated using the geecee program included in the EMBOSS package [36].

Sequence data analysis

Repetitive DNA elements in the qualified BESs, such as transposable elements and SSRs, were identified and masked using Crossmatch search engine (v1.090518), “teleostei” repeat database implemented in Repbase RepeatMasker Edition (20120418) and RepeatMasker program (see for details).

The repeat-masked qualified BESs were subjected to homology search. The eight fish proteome data sets (Atlantic cod: gadMor1.70, Tetraodon: TETRAODON8.70, medaka: MEDAKA1.70, Nile tilapia: Orenil1.0.70, platyfish: Xipmac4.4.2.70, stickleback: BROADS1.70, fugu: FUGU4.70 and zebrafish: Zv9.70) were obtained from Ensembl ( Only the longest protein for each gene was extracted and used for in-house database construction. For construction of genomic sequence database, the genomic sequences of stickleback (gasAcu1) and fugu (fr3) were downloaded from UCSC genome browser ( and others from Ensembl.

BLAST searches were performed with qualified BESs as query with cut-off e-value e−9. The top hit query-subject pairs were extracted using in-house perl script under the criteria in which if multiple query-subject pairs were observed and were overlapped each other, only the most significant pair was considered significant.

Development of SSR markers

The SSR motifs and primer pairs in the qualified BESs were searched by WebSat online application ( with default settings except for product size, which was set to 100–200 bp. For each BAC clone, the SP6-side BES was at first analyzed and if no SSR motif or primer-binding site was found in the sequence, T7-side BES of the same clone was alternatively used. In both sides, the SSR motifs containing over six repeats were considered as real SSRs. In the case where more than two SSR motifs were found in one read, the longest one was used as a representative. Finally, 743 SSRs (89 mono-, 550 di-, 68 tri-, 31 tetra-, 4 penta- and 1 hexa-nucleotide repeats) were selected for primer design.

In addition to the BES-derived SSRs, we also developed microsatellite markers from the genomic library constructed by Ohara et al. [2, 3]. The microsatellites containing CA/GT repeat motifs were isolated according to the protocol of Ohara et al. and primers were designed as described above [2, 3].

Mapping panel

The mapping panel consists of ninety progenies produced by artificial fertilization. Parent fish were caught off Goto Island, Nagasaki Prefecture, Japan and reared in a sea cage until they were matured with approximate body weight 7 kg. Human chorionic gonadotropin (ASKA Pharmaceutical) was intramuscularly administered to the parent fish at 600 IU/kg body weight and eggs and sperm were taken at 45 hours after administration. Fertilized eggs were kept in 500 L seawater at approximately 19°C with 0.5-1 L/min aeration until hatching. The juvenile fish were reared in 500 L seawater at 20-25°C until their body length reached 10 cm. The caudal fin was partially clipped from each progeny as the DNA source and kept in absolute ethanol until use. Genomic DNA of each fish was extracted using DNeasy Blood and Tissue Kit (Qiagen) according to manufacturer’s instructions.

Data acquisition

Genotyping was performed in an 11 μl reaction volume containing 0.5 pmol/μl of unlabelled primer, 0.05 pmol/μl of fluorescence-end-labeled primer with [5’-TET], 1 × buffer, 2.0 mM MgCl2, 0.2 mM dNTP, 1.1 μg of BSA, 0.025 U of EX Taq DNA polymerase (Takara) and 25 ng of template DNA. PCR was performed on a GeneAmp® PCR System 9700 (Applied Biosystems), and the program conditions were 95°C for 2 min for initial denaturation, followed by 30 cycles of 30 sec at 95°C, 1 min at the annealing temperature (52-55°C), 1 min at 72°C and 10 min at 72°C for final extension. Amplification products were mixed with an equal volume of loading buffer (98% formamide, 10 mM EDTA, 0.05 w/v% bromophenol blue), heated for 10 min at 95°C and then immediately cooled on ice. 2 μl of each sample was loaded onto a 6% PAGE-PLUS gel (Amresco) containing 8 M urea and 0.5 × TBE buffer. Electrophoresis was performed in 0.5 × TBE buffer, and after electrophoresis, the gel was scanned and imaged using a FLA-9000 image scanner (GE Healthcare).

Linkage map construction

Genotype data obtained above were subjected for linkage analysis for the male and the female meiosis independently. Marker genotypes were analyzed with LINKMFEX ver. 2.3 ( Linkage analysis was performed using genotype data converted to a backcross format. As the grandparent genotypes were unknown, pairwise analyses were performed, and markers were sorted in linkage groups at a minimum LOD score of 4.0. A goodness-of-fit for Mendelian segregation distortion was tested for all alleles using the chi-square test (p < 0.05, d.f. = 1). Finally, the marker order was determined and double recombination events were checked with MapManagerQTX version 2.0 [37]. The resultant genetic map was visualized using MapChart version 2.2 [38].

The genome length L was estimated using two different methods following Fishman et al. [39]. In the first method (L 1 ), average marker interval was estimated by dividing the summed length of all linkage groups by the number of intervals, and twice the average marker interval was added to each linkage group. In the second method (L 2 ), the length of the each linkage group was multiplied by the factor (m + 1)/(m - 1), where m is the number of framework markers on the linkage groups. Finally, genome coverage c of the linkage map was estimated by calculating c = 1 - e-2dn/L, where d is the average interval of markers, n is the number of markers, and L is the genome length estimated above.

Identification orthologous chromosomes with other fishes

The flanking sequences obtained from all SSR markers assigned to the yellowtail linkage groups were used for the BLASTn search against genomic sequences of medaka, Tetraodon, stickleback, fugu and zebrafish with a cut-off e-value of 0.01. The top hit query-subject was extracted using in-house perl script. In the case where multiple hits were obtained, we defined orthology as follows; let us consider only the first, second and third top hit, if query position of the first and second/third top hit is overlapped each other and quotient of e-value of the first hit divided by that of the second/third hit is greater than 10−3, the hit is considered to be an unclear orthologous pair and rejected. The substantial hits were processed for constructing Oxford grid using Grid Map ver. 3.0a (

Availability of supporting data

All the supporting data are included as additional files.


  1. Kocher TD, Kole C: Genome Mapping and Genomics in Fishes and Aquatic Animals. 2008, Berlin: Springer

    Book  Google Scholar 

  2. Ohara E, Nishimura T, Nagakura Y, Sakamoto T, Mushiake K, Okamoto N: Genetic linkage maps of two yellowtails (Seriola quinqueradiata and Seriola lalandi). Aquaculture. 2005, 244: 41-48. 10.1016/j.aquaculture.2004.10.022.

    Article  CAS  Google Scholar 

  3. Ohara E, Nishimura T, Sakamoto T, Nagakura Y, Mushiake K, Okamoto N: Isolation and characterization of microsatellite loci from yellowtail Seriola quinqueradiata and cross-species amplification within the genus Seriola. Mol Eco Notes. 2003, 3: 390-391. 10.1046/j.1471-8286.2003.00460.x.

    Article  CAS  Google Scholar 

  4. Chen M, Presting G, Barbazuk WB, Goicoechea JL, Blackmon B, Fang G, Kim H, Frisch D, Yu Y, Sun S, Higingbottom S, Phimphilai J, Phimphilai D, Thurmond S, Gaudette B, Li P, Liu J, Hatfield J, Main D, Farrar K, Henderson C, Barnett L, Costa R, Williams B, Walser S, Atkins M, Hall C, Budiman MA, Tomkins JP, Luo M: An integrated physical and genetic map of the rice genome. Plant Cell. 2002, 14: 537-545. 10.1105/tpc.010485.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Budiman MA, Mao L, Wood TC, Wing RA: A deep-coverage tomato BAC library and prospects toward development of an STC framework for genome sequencing. Genome Res. 2000, 10: 129-136.

    PubMed  CAS  PubMed Central  Google Scholar 

  6. International Human Genome Sequencing Consortium: Initial sequencing and analysis of the human genome. Nature. 2001, 409: 860-921. 10.1038/35057062.

    Article  Google Scholar 

  7. Shizuya H, Birren B, Kim UJ, Mancino V, Slepak T, Tachiiri Y, Simon M: Cloning and stable maintenance of 300-kilobase-pair fragments of human DNA in Escherichia coli using an F-factor-based vector. Proc Natl Acad Sci U S A. 1992, 89: 8794-8797. 10.1073/pnas.89.18.8794.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  8. Myosho T, Otake H, Masuyama H, Matsuda M, Kuroki Y, Fujiyama A, Naruse K, Hamaguchi S, Sakaizumi M: Tracing the emergence of a novel sex-determining gene in medaka, Oryzias luzonensis. Genetics. 2012, 191: 163-170. 10.1534/genetics.111.137497.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  9. Buitkamp J, Kollers S, Durstewitz G, Fries R, Welzel K, Schäfer K, Kellermann A, Lehrach H: Construction and characterization of a gridded cattle BAC library. Anim Genet. 2000, 31: 347-351. 10.1046/j.1365-2052.2000.00675.x.

    Article  PubMed  CAS  Google Scholar 

  10. Fahrenkrug SC, Rohrer GA, Freking BA, Smith TP, Osoegawa K, Shu CL, Catanese JJ, de Jong PJ: A porcine BAC library with tenfold genome coverage: a resource for physical and genetic map integration. Mamm Genome. 2001, 12: 472-474. 10.1007/s003350020015.

    Article  PubMed  CAS  Google Scholar 

  11. Vaiman D, Billault A, Tabel-Aoul K, Schibler L, Vilette D, Oustry-Vaiman A, Soravito C: Construction and characterization of a sheep BAC library of three genome equivalents. Mamm Genome. 1999, 10: 585-587. 10.1007/s003359901049.

    Article  PubMed  CAS  Google Scholar 

  12. Thorsen J, Zhu BL, Frengen E, Osoegawa K, de Jong PJ, Koop BF, Davidson WS, Hoyheim B: A highly redundant BAC library of Atlantic salmon (Salmo salar): an important tool for salmon projects. BMC Genomics. 2005, 6: 50-10.1186/1471-2164-6-50.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Liu H, Jiang Y, Wang S, Ninwichian P, Somridhivej B, Xu P, Abernathy J, Kucuktas H, Liu Z: Comparative analysis of catfish BAC end sequences with the zebrafish genome. BMC Genomics. 2009, 10: 592-10.1186/1471-2164-10-592.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Katagiri T, Asakawa S, Minagawa S, Shimizu N, Hirono I, Aoki T: Construction and characterization of BAC libraries for three fish species; rainbow trout, carp and tilapia. Anim Genet. 2001, 32: 200-204. 10.1046/j.1365-2052.2001.00764.x.

    Article  PubMed  CAS  Google Scholar 

  15. Xu P, Li J, Li Y, Cui R, Wang J, Wang J, Zhang Y, Zhao Z, Sun X: Genomic insight into the common carp (Cyprinus carpio) genome by sequencing analysis of BAC-end sequences. BMC Genomics. 2011, 12: 188-10.1186/1471-2164-12-188.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  16. Whitaker HA, McAndrew BJ, Taggart JB: Construction and characterization of a BAC library for the European sea bass Dicentrarchus labrax. Anim Genet. 2006, 37: 526-10.1111/j.1365-2052.2006.01514.x.

    Article  PubMed  CAS  Google Scholar 

  17. Kuhl H, Beck A, Wozniak G, Canario AV, Volckaert FA, Reinhardt R: The European sea bass Dicentrarchus labrax genome puzzle: comparative BAC-mapping and low coverage shotgun sequencing. BMC Genomics. 2010, 11: 68-10.1186/1471-2164-11-68.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Wang CM, Lo LC, Feng F, Gong P, Li J, Zhu ZY, Lin G, Yue GH: Construction of a BAC library and mapping BAC clones to the linkage map of Barramundi, Lates calcarifer. BMC Genomics. 2008, 9: 139-10.1186/1471-2164-9-139.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Ammiraju JS, Luo M, Goicoechea JL, Wang W, Kudrna D, Mueller C, Talag J, Kim H, Sisneros NB, Blackmon B, Fang E, Tomkins JB, Brar D, MacKill D, McCouch S, Kurata N, Lambert G, Galbraith DW, Arumuganathan K, Rao K, Walling JG, Gill N, Yu Y, SanMiguel P, Soderlund C, Jackson S, Wing RA: The Oryza bacterial artificial chromosome library resource: construction and analysis of 12 deep-coverage large-insert BAC libraries that represent the 10 genome types of the genus Oryza. Genome Res. 2006, 16: 140-147.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Fuji K, Yoshida K, Hattori K, Ozaki A, Araki K, Okauchi M, Kubota S, Okamoto N, Sakamoto T: Identification of the sex-linked locus in yellowtail, Seriola quinqueradiata. Aquaculture. 2010, 308: S51-S55.

    Article  CAS  Google Scholar 

  21. Genet C, Dehais P, Palti Y, Gao G, Gavory F, Wincker P, Quillet E, Boussaha M: Analysis of BAC-end sequences in rainbow trout: content characterization and assessment of synteny between trout and other fish genomes. BMC Genomics. 2011, 12: 314-10.1186/1471-2164-12-314.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  22. Xu P, Wang S, Liu L, Peatman E, Somridhivej B, Thimmapuram J, Gong G, Liu Z: Channel catfish BAC-end sequences for marker development and assessment of syntenic conservation with other fish species. Anim Genet. 2006, 37: 321-326. 10.1111/j.1365-2052.2006.01453.x.

    Article  PubMed  CAS  Google Scholar 

  23. Shirak A, Grabherr M, Di Palma F, Lindblad-Toh K, Hulata G, Ron M, Kocher TD, Seroussi E: Identification of repetitive elements in the genome of Oreochromis niloticus: tilapia repeat masker. Mar Biotechnol (NY). 2010, 12: 121-125. 10.1007/s10126-009-9236-8.

    Article  CAS  Google Scholar 

  24. Roest Crollius H, Jaillon O, Dasilva C, Ozouf-Costaz C, Fizames C, Fischer C, Bouneau L, Billault A, Quetier F, Saurin W, Bernot A, Weissenbach J: Characterization and repeat analysis of the compact genome of the freshwater pufferfish Tetraodon nigroviridis. Genome Res. 2000, 10: 939-949. 10.1101/gr.10.7.939.

    Article  PubMed  CAS  Google Scholar 

  25. Nelson JS: Fishes of the World. 2006, New Jersey: Wiley, 4

    Google Scholar 

  26. Sarropoulou E, Nousdili D, Magoulas A, Kotoulas G: Linking the genomes of nonmodel teleosts through comparative genomics. Mar Biotechnol (NY). 2008, 10: 227-233. 10.1007/s10126-007-9066-5.

    Article  CAS  Google Scholar 

  27. Liu S, Rexroad CE, Couch CR, Cordes JF, Reece KS, Sullivan CV: A microsatellite linkage map of striped bass (Morone saxatilis) reveals conserved synteny with the three-spined stickleback (Gasterosteus aculeatus). Mar Biotechnol (NY). 2012, 14: 237-244. 10.1007/s10126-011-9407-2.

    Article  CAS  Google Scholar 

  28. Ojima Y, Yamamoto K: Cellular DNA contents of fishes determined by flow cytometry. La Kromosomo II. 1990, 57: 1871-1888.

    Google Scholar 

  29. Massault C, Bovenhuis H, Haley C, de Koning DJ: QTL mapping designs for aquaculture. Aquaculture. 2008, 285: 23-29. 10.1016/j.aquaculture.2008.06.040.

    Article  CAS  Google Scholar 

  30. Streelman JT, Albertson RC, Kocher TD: Genome mapping of the orange blotch colour pattern in cichlid fishes. Mol Ecol. 2003, 12: 2465-2471. 10.1046/j.1365-294X.2003.01920.x.

    Article  PubMed  CAS  Google Scholar 

  31. Roberts RB, Ser JR, Kocher TD: Sexual conflict resolved by invasion of a novel sex determiner in Lake Malawi cichlid fishes. Science. 2009, 326: 998-1001. 10.1126/science.1174705.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  32. Locke J, Podemski L, Aippersbach N, Kemp H, Hodgetts R: A physical map of the polytenized region (101EF-102F) of chromosome 4 in Drosophila melanogaster. Genetics. 2000, 155: 1175-1183.

    PubMed  CAS  PubMed Central  Google Scholar 

  33. Ozaki A, Yoshida K, Fuji K, Kubota S, Kai W, Aoki J, Kawabata Y, Suzuki J, Akita K, Koyama T, Nakagawa M, Hotta T, Tsuzaki T, Okamoto N, Araki K, Sakamoto T: Quantitative trait loci (QTL) associated with resistance to a parasitic disease (Benedenia seriolae) in yellowtail (Seriola quinqueradiata) through genome wide analysis. PLOS ONE. 2013, 8: e64987-10.1371/journal.pone.0064987.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  34. Ewing B, Hillier L, Wendl MC, Green P: Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998, 8: 175-185. 10.1101/gr.8.3.175.

    Article  PubMed  CAS  Google Scholar 

  35. Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998, 8: 186-194.

    Article  PubMed  CAS  Google Scholar 

  36. Rice P, Longden I, Bleasby A: EMBOSS: the European molecular biology open software Suite. Trends Genet. 2000, 16: 276-277. 10.1016/S0168-9525(00)02024-2.

    Article  PubMed  CAS  Google Scholar 

  37. Manly KF, Olson JM: Overview of QTL mapping software and introduction to Map manager QT. Mamm Genome. 1999, 10: 327-334. 10.1007/s003359900997.

    Article  PubMed  CAS  Google Scholar 

  38. Voorips RE: Mapchart software for the graphical presentation of linkage maps and QTLs. J Hered. 2002, 93: 77-78. 10.1093/jhered/93.1.77.

    Article  Google Scholar 

  39. Fishman L, Kelly AJ, Morgan E, Willis JH: A genetic map in the Mimulus guttatus species complex reveals transmission ratio distortion due to heterospecific interactions. Genetics. 2001, 159: 1701-1716.

    PubMed  CAS  PubMed Central  Google Scholar 

Download references


We thank Misato Fujii for providing technical help in SSR genotyping, Kazuharu Nomura in NRIFS for his support with the analysis of measure an estimated genome size of yellowtail. This research was supported by the Program for Promotion of Basic and Applied Researches for Innovations in Bio-oriented Industry (BRAIN) and by grants-in-aid for Scientific Research “New Technology of Fish Breeding” from Fisheries Research Agency of Japan.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Takashi Sakamoto.

Additional information

Competing interests

The authors declare they have no competing interests.

Authors’ contributions

KF performed all molecular experiments, constructed the linkage map and drafted the manuscript. TK analyzed all data and drafted the manuscript. WK contributed to data analysis. SK helped SSR marker genotyping. KY performed fish breeding. AO took part in developing the linkage analysis pipeline. JA, YK, KA, and TT provided laboratory facilities and support this program in all steps. NO and TS conceived and overlooked the project and reviewed the manuscript. All authors read and approved the final manuscript.

Kanako Fuji, Takashi Koyama contributed equally to this work.

Electronic supplementary material


Additional file 1:SSR markers in the yellowtail map. Marker name, linkage group, polymorphic information (1: polymorphic in both female and male, 2: polymorphic in only female, 3: polymorphic in only male), genotype of dam and sire, kind of repeat motif number, repeat type, primer sequences, annealing temperature, PCR product size, and GenBank accession number are shown. (XLS 303 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( ) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Cite this article

Fuji, K., Koyama, T., Kai, W. et al. Construction of a high-coverage bacterial artificial chromosome library and comprehensive genetic linkage map of yellowtail Seriola quinqueradiata. BMC Res Notes 7, 200 (2014).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: