- Research article
- Open Access
Construction of a high-coverage bacterial artificial chromosome library and comprehensive genetic linkage map of yellowtail Seriola quinqueradiata
BMC Research Notes volume 7, Article number: 200 (2014)
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 . 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 , to develop sequence-tagged connectors , and to sequence the genome itself  because of their insert size capacity, reproducibility and stability as the DNA sample . 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 . BAC libraries have been developed in several domestic animals, e.g. cattle , pig  and sheep , and in fishes, salmon , catfish , rainbow trout, carp, tilapia [14, 15], European sea bass [16, 17] and barramundi  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 . 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 .
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%) (http://esper.lab.nig.ac.jp/genome-composition-database/).
Preliminary survey of genomic content using BESs
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%) , common carp (17.3%) , channel catfish (11.9%)  and Nile tilapia (14.0%)  and comparable to that in Tetraodon (6.2%) and in fugu (4.3%) .
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).
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 . 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].
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 . 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 .
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.
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 . 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.
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 . 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 .
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 http://www.repeatmasker.org/ 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: Xipmac18.104.22.168, stickleback: BROADS1.70, fugu: FUGU4.70 and zebrafish: Zv9.70) were obtained from Ensembl (ftp://ftp.ensembl.org/). 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 (http://hgdownload.soe.ucsc.edu/) 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 (http://wsmartins.net/websat/) 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].
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.
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 (http://www.uoguelph.ca/~rdanzman/software.htm). 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 . The resultant genetic map was visualized using MapChart version 2.2 .
The genome length L was estimated using two different methods following Fishman et al. . 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 (http://cbr.jic.ac.uk/dicks/software/Grid_Map/).
Availability of supporting data
All the supporting data are included as additional files.
Kocher TD, Kole C: Genome Mapping and Genomics in Fishes and Aquatic Animals. 2008, Berlin: Springer
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.
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.
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.
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.
International Human Genome Sequencing Consortium: Initial sequencing and analysis of the human genome. Nature. 2001, 409: 860-921. 10.1038/35057062.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Nelson JS: Fishes of the World. 2006, New Jersey: Wiley, 4
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.
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.
Ojima Y, Yamamoto K: Cellular DNA contents of fishes determined by flow cytometry. La Kromosomo II. 1990, 57: 1871-1888.
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.
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.
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.
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.
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.
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.
Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998, 8: 186-194.
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.
Manly KF, Olson JM: Overview of QTL mapping software and introduction to Map manager QT. Mamm Genome. 1999, 10: 327-334. 10.1007/s003359900997.
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.
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.
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.
The authors declare they have no competing interests.
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)
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) doi:10.1186/1756-0500-7-200
- BAC library
- Linkage map