Skip to main content

A high-resolution linkage map for comparative genome analysis and QTL fine mapping in Asian seabass, Lates calcarifer

Abstract

Background

High density linkage maps are essential for comparative analysis of synteny, fine mapping of quantitative trait loci (QTL), searching for candidate genes and facilitating genome sequence assembly. However, in most foodfish species, marker density is still low. We previously reported a first generation linkage map with 240 DNA markers and its application to preliminarily map QTL for growth traits in Asian seabass (Lates calcarifer). Here, we report a high-resolution linkage map with 790 microsatellites and SNPs, comparative analysis of synteny, fine-mapping of QTL and the identification of potential candidate genes for growth traits.

Results

A second generation linkage map of Asian seabass was developed with 790 microsatellite and SNP markers. The map spanned a genetic length of 2411.5 cM, with an average intermarker distance of 3.4 cM or 1.1 Mb. This high density map allowed for comparison of the map with Tetraodon nigroviridis genome, which revealed 16 synteny regions between the two species. Moreover, by employing this map we refined QTL to regions of 1.4 and 0.2 cM (or 400 and 50 kb) in linkage groups 2 and 3 in a population containing 380 progeny; potential candidate genes for growth traits in QTL regions were further identified using comparative genome analysis, whose effects on growth traits were investigated. Interestingly, a QTL cluster at Lca371 underlying growth traits of Asian seabass showed similarity to the cathepsin D gene of human, which is related to cancer and Alzheimer's disease.

Conclusions

We constructed a high resolution linkage map, carried out comparative mapping, refined the positions of QTL, identified candidate genes for growth traits and analyzed their effects on growth. Our study developed a framework that will be indispensable for further identification of genes and analysis of molecular variation within the refined QTL to enhance understanding of the molecular basis of growth and speed up genetic improvement of growth performance, and it also provides critical resource for future genome sequence assembly and comparative genomics studies on the evolution of fish genomes.

Background

Most economically important traits are quantitative in nature and are determined by many genes and gene complex where are described as quantitative trait loci (QTL) [1]. Traditional methods of genetic improvement of quantitative traits have relied mainly on phenotype and pedigree information [1], which are easily influenced by environmental factors. Genetic markers have made it possible to detect QTL that are significantly associated with traits [2], and made selection more effective. Genetic response can be improved by including the QTL in marker-assisted selection, which is a method of selection that makes use of phenotypic, genotypic and pedigree data [3].

Linkage maps are essential for mapping QTL [4]. In the past, genotyping of many markers was expensive, therefore, specific experimental designs were developed to reduce the impact of having fewer markers on statistical power [1]. More recently, however, high throughput methods have been developed to genotype markers such as microsatellites [5] and single nucleotide polymorphisms (SNP) [6], which have significantly reduced the cost. Linkage maps have been constructed for a number of foodfish species, such as salmon [7], rainbow trout [8, 9], catfish [10], tilapia [11], grass carp [12], common carp [13], Asian seabass [14], European seabass [15] and Japanese flounder [16], gilthead seabream Sparus aurata[17–19] using RAPD, AFLP and microsatellites. Only in a few species, linkage maps solely based on codominant DNA markers (microsatellites and SNPs) were constructed. Most linkage maps in food fish species are not dense in comparison to these linkage maps in model fish species (e.g. Zebrafish [20]), chicken [21], live stock species (e.g. cattle [22, 23], and pig [24]) and agronomic plant species (e.g. barley, soybean, grapevine [25–27]). QTL mapping in foodfish species is still in its infancy [6]. Only in a few species, such as Asian seabass [6], salmon [28], tilapia [29], Japanese flounder [30], rainbow trout [31] and European seabass [32], QTL for growth, meat quality, stress and disease resistance have been mapped in large genomic regions due to lack of a high resolution linkage map.

Linkage map with sequence-based markers is also a platform for comparative genome studies [33–36]. Recent comparative genome analyses based on genetic maps have already provided new insights into genome organization, evolution, and function across different organisms [12, 33, 35, 37]. For example, comparison of the Caenorhabditis briggsae genetic map and the Caenorhabditis elegans genome reveals extensive conservation of chromosome organization and synteny despite a very long divergence time (80 to 110 million years), suggesting that natural selection operates at the level of chromosomal organization [37]. In another study, a genetic linkage map of the blind Mexican cavefish Astyanax mexicanus has been successfully applied to predict candidate quantitative trait genes relating to rib number and eye size by anchoring cavefish QTLs to the zebrafish genome [36]. BLAST searches of sequences of mapped markers of grass carp against the whole genome sequence of zebrafish revealed substantial macrosynteny relationship and extensive colinearity of markers between grass carp and zebrafish [12]. Identification of conserved synteny blocks across fish genomes would help to unravel ancestral genome architecture of fish and transfer genome information from model fish species to non-model foodfish species.

Asian seabass, Lates calcarifer, also called Barramundi, is one of the important foodfish species. This species has been cultured for more than 20 years in brackish-water ponds and in recent years in floating cages. The global annual production of Asian seabass was currently 400,000 metric tons according to FAO statistics [38]. In the past few years, we started a breeding program for Asian seabass [39] and developed a number of genomic tools such as microsatellites [40–42], SNPs in genes [43], microRNA [44], a linkage map with 240 microsatellites [14], BAC and cDNA libraries [43, 45, 46] and a BAC-based physical map [44] to facilitate the selective breeding program. The linkage map has been used to map QTL for growth traits, and significant QTL for growth traits were mapped on linkage groups 2 and 3 [39, 47]. However, due to the lack of markers in the QTL regions, it is impossible to map the QTL in smaller chromosomal regions. For fine mapping QTL, comparative analysis of synteny and searching for candidate genes in QTL region, a high-density linkage map is essential.

In this report we present a second generation linkage map of Asian seabass. The current updated version of the Asian seabass consensus linkage map is a considerable improvement compared with the previous version [14]. This linkage map allowed for carrying out comparative mapping of synteny between Asian seabass and Tetraodon nigroviridis, and enabled fine mapping of QTL for growth traits. In addition, we identified potential candidate genes in QTL for growth traits, and defined the phenotypic consequences of alternative candidate gene alleles.

Results

Identification and genotyping of DNA markers

A total of 4300 clones collected from libraries enriched for CA-, GA-, CAA- GACA- and GATA-microsatellites were sequenced in both directions. Two thousand and eight hundred clones contained microsatellites, yielding 1520 unique sequences. Among the 1520 sequences, 1280 had enough flanking regions for primer design. The first set of 920 primer pairs was used to amplify three parents from two reference families [14] for linkage mapping. Six hundreds and sixty primers were selected to genotype the two reference families including 96 individuals (3 parents and 93 offspring) due to the fact that these primers could amplify easily scorable PCR products. Among the 660 microsatellite markers, 280 were tetranucleotide microsatellites, which could be more easily scored than the di-nucleotide microsatellites (see Additional file 1). Ten SNPs in nine genes were genotyped by direct sequencing each individual in the two reference families.

Linkage mapping

Of 851 informative markers including 240 markers mapped in the first generation linkage map [14], 822 markers were assigned to linkage groups by two-point linkage analysis with LOD scores >3.0 using CRIMAP (Green et al. 1990). From this, a total of 790 (97%) markers were mapped by multipoint linkage analysis. Among these markers, 53 initiated with LcaB were microsatellites isolated from BAC clones, 53 LcaE from ESTs and 10 SNPs from genes (Additional file 1). Details about primer sequences, GenBank accession number, annealing temperature for PCR, PCR product size, and locations of the 790 markers are summarized in Additional file 1.

In most regions, the order of the markers on the new map was consistent with the previous map [14], although some regions were rearranged through incorporating the new markers and correcting old marker data. The current sex-averaged map spanned 2411.5 cM of the Asian seabass genome (Figure 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12). In the map, the intermarker distance was 3.0 cM. These 790 markers were located in 501 unique locations on the linkage map of the Asian seabass, with an average inter-location space of 4.8 cM (Table 1).

Figure 1
figure 1

The second generation linkage map of Asian seabass (LGs 1-2). Female (F) and male (M) maps are shown on the left and right, respectively, and the consensus map is shown in the center. The same loci are connected with solid lines.

Figure 2
figure 2

The second generation linkage map of Asian seabass (LGs 3-4).

Figure 3
figure 3

The second generation linkage map of Asian seabass (LGs 5-6).

Figure 4
figure 4

The second generation linkage map of Asian seabass (LGs 7-8).

Figure 5
figure 5

The second generation linkage map of Asian seabass (LGs 9-10).

Figure 6
figure 6

The second generation linkage map of Asian seabass (LGs 11-12).

Figure 7
figure 7

The second generation linkage map of Asian seabass (LGs 13-14).

Figure 8
figure 8

The second generation linkage map of Asian seabass (LGs 15-16).

Figure 9
figure 9

The second generation linkage map of Asian seabass (LGs 17-18).

Figure 10
figure 10

The second generation linkage map of Asian seabass (LGs 19-20).

Figure 11
figure 11

The second generation linkage map of Asian seabass (LGs 21-21).

Figure 12
figure 12

The second generation linkage map of Asian seabass (LGs 23-24).

Table 1 Summary of the second generation linkage group of Asian seabass

Sex-specific maps were also constructed. The length of the male map was 2674.6 cM, with an average intermarker distance of 3.4 cM whereas the female map was 2294.8 cM long, with an average intermarker spacing of 2.9 cM. The male map comprised linkage groups ranging in length from 6.4 to 244.4 cM while the female map contained linkage groups with a length ranging from 22.3 to 225.6 cM.

The distances between LcaB201 and LcaTe0215 at the end region of LG1, are 0 and 47.6 cM on female and male map respectively, while the distances between Lca270 and LcaTe0382 in the near middle region of LG1, possible proximal to centromere region, are 39.2 and 0 cM on female and male map respectively. Similarly, in the end regions between LcaE142 and LcaB128 on LG2, Lca1508 and LcaTe0045 on LG3, LcaTe0194 and LcaE169 on LG5 and Lca114 and Lca454 on LG9, females had much lower recombination rates in telomeric regions than males. While in the near middle regions, possibly within regions proximal to the centromere between Lca359 and Lca137 on LG3, LcaB045 and Lca996 on LG9, Lca365 and LcaTe0191 on LG10, Lca1012 and Lca565 on LG17, Lca512 and Lca564 on LG23, LcaTe0485 and LcaTe0214 on LG24, recombination rates were much higher in females. It showed that females have much lower recombination rates in telomeric regions than males, while recombination rates were much higher in females within regions proximal to the centromere.

Comparative genome analysis

We used BLAT to identify homologs of the sequence-based markers of Lates calcarifer on the map. Sixty seven markers had a homolog and could be assigned to the chromosomes of Tetraodon nigroviridis. We compared the Lates calcarifer genetic map with assembled genome sequences of Tetraodon nigroviridis, identifying conserved synteny blocks in 16 of the 24 Lates calcarifer linkage groups, each of which contained 3 to 8 markers. The largest synteny block conserved between Lates calcarifer and Tetraodon nigroviridis was found in LG5 with 8 markers spanning 48 cM in the Lates calcarifer linkage group and their best matches spanning 4 Mb in chromosome 18 of the Tetraodon nigroviridis genome (Figure 13, 14, Additional file 2). Blasting sequences of DNA markers of Asian seabass against whole genome sequences of zebrafish and medaka detected only very few conserved synteny (data not shown).

Figure 13
figure 13

Lates calcarifer (Lca) linkage map and Tetraodon nigroviridis (Tni) synteny. Homologous Lca (empty) and Tni (dashed) chromosomes are shown with lines connecting homologous markers.

Figure 14
figure 14

Lates calcarifer (Lca) linkage map and Tetraodon nigroviridis (Tni) synteny (continued). Homologous Lca (empty) and Tni (dashed) chromosomes are shown with lines connecting homologous markers.

Fine mapping QTL

Additional 33 novel microsatellites located near QTL for growth traits on LGs 2 and 3 were selected and mapped with QTL panel containing 380 offspring. The linkage map was significantly improved for further QTL mapping as compared to the previous map by enhancing the density of markers for refinement of QTL positions (Figure 15). QTL analysis was carried out with genotype data of the markers and phenotypic data of all the progeny in the QTL mapping panel.

Figure 15
figure 15

QTL for growth traits identified on Asian seabass. See Table 2 for details about the effects of QTL. The position of the QTL is indicated on the right of the linkage groups. The QTL bars indicate experiment-wise LOD support confidence interval in which the inner line indicates position of maximum LOD score. The highlighted region on LGs shows QTL interval between two flanking markers. The peaks of qBW2-a and BW2-b were flanked by Lca182 and Lca287, and closely linked to Lca287 with a distance of 2 cM. The peaks of the other QTL (qBW2-b, c, d, e and qBW3) were detected near the positions of markers Lca562, Lca371, Te0359, Lca250 and Lca137.

QTL affecting body weight, total length and standard length were identified on an experiment-wise scale. The experiment-wise LOD significance thresholds were 3.6, 5.5 and 5.4 for body weight, total length and standard length, respectively, while the linkage-group-wise LOD significance thresholds varied from 2.3 to 5.5 (Table 2). Eleven QTL controlling body weight, total length and standard length were detected on LGs 2 and 3. Multiple QTL Model (MQM) mapping with initial QTL did not change the results.

Table 2 Location of QTL and magnitude of QTL effects on growth traits in Asian seabass

The three QTL qBW2-a, qTL2-a and qSL2-a near the marker Lca287 showed high percentage of phenotypic variance explained (PVE) of 30.2, 53.1 and 53.8%, respectively (Table 2). These QTL together with qBW2-b are consistent to the previous QTL mapping results and their positions were finely refined as discussed later on.

Due to the more markers integrated in the QTL region, five QTL, i.e. qBW2-c, qBW2-d, qBW2-e, qTL2-c and qTL2-d, were newly detected. Although LOD values of qSL2-a and qSL2-c were not higher than the significant threshold, we still listed qSL2-a and qSL2-c to Table 2 as they showed obvious peaks and the LOD scores were very close to the threshold.

On LG2, the peaks of qBW2-a and qBW2-b were flanked by Lca182 and Lca287 with distance 6.0 cM, and closely linked to Lca287 with a distance of 2 cM. Similarly, the peaks of the other QTL (qBW2-b, c, d and e) were located in a small region between flanking markers, with intervals of 6.1, 1.4, 0.6, 0.2 and 0.1 cM respectively. On LG3, qBW3 was located in an interval of 0.1 cM between Lca137 and Lca159. The LOD peaks were detected near the positions of markers Lca562, Lca371, LcaTe0359, Lca250 and Lca137 respectively (Table 2, Figure 15).

Potential candidate genes in refined QTL

As positions of some microsatellites collocating with QTL peaks, sequences of these DNA markers were blasted against the whole genome sequences of zebrafish and medaka. We detected some potential candidate genes for growth. Among these markers, a QTL cluster including qBW2-c, qTL2-c and qSL2-c were located in between markers Lca371 and Lca480 with distance of 1.4 cM, and QTL peaks collocated with the position of Lca371, whose sequence hit to genome sequences at 16.325 Mb on the chromosome 18 (Dr-18) in zebrafish, and at 2.311 Mb chromosome 3 (Ol-3) in medaka (Figure 16). When looking closely, we found that the sequence of the marker Lca371 showed very high similarity of the cathepsin D gene of medaka and zebrafish. We also found that the Lca371 showed similarity to the cathepsin D gene (Genbank accession nos: TC168775 and TC156330) of rainbow trout. Comparing the sequences of Lca371 to all 24,000 EST sequences generated by our group (unpublished data); we found that Lca371 was located in the cathepsin D gene of Asian seabass.

Figure 16
figure 16

Mapping of QTL for growth traits. a, QTL for body weight on LG2; b, QTL for body weight on LG3; c, QTL for total length on LG2; d, QTL for total length on LG3; e, QTL for standard length on LG2; f, QTL for standard length on LG3. The lines were drawn by plotting the LOD scores at each marker as well as at 1.0 cM intervals along the linkage group. The arrow shows the position of genes and their markers where overlapped with QTL peaks.

The peak of qBW2-e was located between the markers Lca568/LcaTe0265 and Lca250 with distance of 0.2 cM (Figure 16). The sequence of the marker Lca568 hit to genome sequences at 29.898 Mb on the chromosome 3 (Ol-3) in medaka, and some of the sequences of the marker Lca568 showed very high similarity of the KCTD15 gene of medaka. The sequences of LcaTe0265 hit to genome sequences at 35.626 Mb on the chromosome 19 (Dr-19) in zebrafish, and some of the sequences of the marker LcaTe0265 showed very high similarity of the csmd2 gene of zebrafish.

For Lca562, LcaTe0359 and Lca137, there were no hits to genes in medaka and zebrafish genomes.

Effects of gene allele substitution

A two-way ANOVA was performed on the 380 progeny using four allelic combinations (m1f1, m1f2, m2f1 and m2f2) from markers in the three candidate genes of cathepsin D, KCTD15 and csmd2 in order to investigate associations between traits and genotypes of these genes. The phenotype values of each allelic combination of these three genes are listed in Table 2 and Figure 17. Significant differences of phenotype means among different allelic combinations were identified, revealing the effects of alternate gene alleles inherited from the parents. Progeny with m2f2 genotype at the marker Lca371 located in cathepsin D, showed the highest phenotype values. Similarly, progeny with m1f2 genotype at the marker Lca568 located in KCTD15 and LcaTe0265 in csmd2, showed the highest phenotype values. These results suggested the effect of these genes on growth related traits.

Figure 17
figure 17

The effects of genotypes of cathepsin D (left), KCTD15 and csmd2 (right). The letters on the top of bars indicate the level of difference in phenotypic value between genotypes. The same letter indicates that the difference is statistically insignificant, whereas different letters represent statistically significant (Bonferroni T tests at α < 0.01) difference.

Discussion

Map density and recombination rate

The linkage map presented in this paper represents one of the most dense linkage maps based on microsatellites and SNPs for foodfish species. This second generation linkage map of Asian seabass contained 3.3 times as many sequence-based co-dominant loci as did the previous linkage map [14], and more sequence-based co-dominant markers than the linkage maps of major foodfish species, such as salmon [7], tilapia [11], common carp [13], grass carp [12], Japanese flounder [16], catfish [10] and European seabass [15]. However, in comparison to the linkage maps of model fish species (e.g. zebrafish [20]), livestock (e.g. cattle [22, 23], pig [24]), chicken [21]) and agronomical plant species (e.g. rice, barley, soybean, grapevine, [25–27, 48], the marker density of the Asian seabass linkage map is still lower. This is mainly because the whole genome sequence and large number of SNPs in Asian seabass are still not available. Fortunately, sequencing the whole genome of Asian seabass is in agenda. The third generation of linkage map could be based on large number of SNPs. The current Asian seabass map spanned 2411.5 cM, and it is estimated to span much more of the Asian seabass genome in comparison with the previous linkage groups [14]. It had a resolution of 3.4 cM, which is sufficient for fine mapping QTL for future marker-assisted selection.

The average recombination rate across all linkage groups is approximately 3.4 cM/Mb in Asian seabass, which is higher than that in zebrafish (1.35 cM/Mb) [20], catfish (1.65 cM/Mb) [10], tilapia (1.3 cM/Mb) [11], grass carp (1.2 cM/Mb) [12], human (1.20 cM/Mb), mouse (0.5 cM/Mb) [49] and lower than the plant Arabidopsis thaliana (5 cM/Mb; calculated based on data from The Arabidopsis Information Resource website). Based on the genome size of 700 Mb [14, 50] we estimated the average intermarker distance to be approximately 0.88 Mb in the new linkage map which consisting 790 markers. This suggests that QTL, if identified, can be narrowed down to rather small genomic regions. The enhanced map will be invaluable not only for QTL and gene mapping but also for comparative genome analysis and assembling the Asian seabass genome once the whole genome sequencing available.

Sex-specific patterns of recombination

Different recombination rates between the sexes have been reported in mammals [8] and other fish species [15, 17, 51] with female map distances usually greater than those in male maps [8, 52]. In Japanese flounder, the recombination rate was unusually higher in males (7.4 times) compared to females [16]. Different to the previous map [14], the overall sizes of the male and female maps are comparable (female: 2294.0 cM; male: 2674.6 cM). Females had much lower recombination rates in telomeric regions than males, while recombination rates were much higher in females within regions proximal to the centromere. Females have much lower recombination rates in telomeric regions than males on LGs 1, 2, 3, 5 and 9, while recombination rates are much higher in females within regions possibly proximal to the centromere on LGs 1, 3, 9, 10, 17, 23 and 24. This is similar to linkage map of rainbow trout [8].

QTL fine mapping

To improve the utility of the QTL in MAS, and to move toward the positional cloning of candidate genes, fine mapping of the QTL to a smaller region of the chromosome is necessary [53, 54]. In foodfish species, although QTL mapping has been conducted in a few species, such rainbow trout [31, 55], salmon [56], Japanese flounder [30] and tilapia [29], QTL were mapped in large genomic regions, usually bigger than 10 cM. In Asian seabass, the number of available markers on our previous map [14] limited the possibility for fine mapping QTL in Asian seabass; therefore the development of additional markers in the QTL regions is an important target in this study. As compared to the previous QTL mapping results, the positions of QTL were much refined with the markers from the enhanced map (Figure 15). The peaks of qBW2-a and BW2-b were flanked by Lca182 and Lca287 with distance reduced from 11.3 cM on previous map to 6.0 cM, and closely linked to Lca287 with a distance of 2.0 cM. Similarly, the peaks of the other QTL (qBW2-b, c, d and e) was refined from 8.3 cM on previous map to a small region between flanking markers, with intervals 6.1, 1.4, 0.6, 0.2 and 0.1 cM respectively. On LG3, qBW3 was from 9.6 cM to 0.1 cM. The refined QTL supply a basis for identifying potential candidate genes located in these refined QTL through comparative genome analysis. Fortunately, the complete genome sequences of several mode fish species (i.e. zebrafish, medaka, freshwater pufferfish, fugu and stickleback) [57–60] are available, and the sequencing of genomes of several important foodfish species (common carp, salmon and tilapia) are in progress.

Comparative genome analysis

Methods to map sequence based markers and to compare maps with relevant model systems are crucial to extend genomic-level analysis to non-model species [54]. According to the phylogenetic tree including over 30 fish species, constructed using 12 mitochondrial genes in our previous study, Asian seabass is more closely related to T. nigroviridis than to zebrafish and medaka [61]. However, in our previous comparative mapping, no synteny block could be identified between Aisan seabass and T. nigroviridis although flanking sequences of 55 microsatellites showed high similarity to known genomic DNA sequences of T. nigroviridis[14]. Herein, we developed our second generation linkage map with high density and demonstrated that, 16 syteny blocks of Asian seabass were syntenic to the 16 counterparts of T. nigroviridis chromosomes. The result showed certain colinearity for the 16 syntenic chromosome/linkage pairs between the two foodfish species, T. nigroviridis and Asian seabass. The conserved syntenies identified here between the Asian seabass and T. nigroviridis should facilitate studies on genome evolution and analysis of structural genome, but more importantly should facilitate functional inference of genes in Asian seabass. 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 ESTs or gene sequences on the linkage map of Asian seabass should enhance comparative mapping, thereby transferring genome information from model species to Asian seabass.

Identification of candidate genes in QTL

Maps with sequence-based markers are useful not only for comparative genomics but also resolving mapped genomic regions to a tractable number of candidate genes, especially if there is synteny with related model species. After refining the positions of QTL, the regions of candidate genes were delimited more precisely. The LOD peaks were detected near the positions of markers Lca562, Lca371, LcaTe0359, Lca250 and Lca137 respectively, therefore we examined regions of the sequenced zebrafish and medaka genome corresponding to the QTL regions on Asian seabass, using the markers overlapped on QTL LOD peaks. Such an approach allows us to prioritize our research effort on candidate genes that are found collocating with QTL LOD peaks. Candidate gene approaches have been successfully used for identification of QTL [62, 63]. Three putative growth-related genes, cathepsin D, KCTD15 and csmd2, were found to be positioned directly at QTL peaks between small region of 0.2-1.4 cM, making them the strongest candidate genes for growth traits. Further two-way ANOVA revealed that allelic substitution at these two genes showed significant effects on growth-related traits. Certain allelic combinations showed significantly higher values of the growth traits. The candidate genes are rare examples of QTL fine mapping in foodfish species. Interestingly, the QTL cluster at Lca371 underlying growth traits of Asian seabass showed similarity to the cathepsin D gene of human, which is related to cancer [64] and Alzheimer's disease [65] and has become a hotspot of human genetics study. The genotypes at these genes may be useful for growth improvement through marker-assisted selection, gene cloning and functional analysis. In humans, KCTD15 may be associated with obesity [66], while csmd2 may be an oligodendroglioma suppressor [67]. Therefore it is interesting to further study the mechanisms underlying the associations between polymorphisms in these two genes with growth traits in Asian seabass.

Conclusions

We constructed a second generation genetic linkage map and carried out comparative mapping of synteny between Asian seabass and T. nigroviridis. Moreover, we applied this map for refinement of QTL for growth traits using extensive progeny testing with defined recombination within the QTL region. We estimated their locations within short intervals, identified potential candidate genes, and further defined the phenotypic consequences of alternative candidate gene alleles. This second generation linkage map should facilitate the advancement of genetic studies for a wide variety of complex traits in foodfish species. In the future, more SNP markers should be identified by using next generation sequencing [68] to enable whole genome association studies [69] to facilitate genomic selection [70] and to understand the genetic basis of phenotypic variation of important traits [71, 72].

Methods

Reference families for linkage mapping and QTL analysis

The reference families used for the construction of the first generation linkage map [14] were used for the construction of the second generation linkage map. Briefly, a whole broodstock containing 94 brooders, including 48 males and 46 females collected from the wild in Southeast Asia four years ago, were genotyped with nine polymorphic microsatellites as described previously [14, 41]. One female and two male brooders were selected for constructing a mapping panel because of their high allelic diversity and genetic differences. By crossing the female and two male brooders, millions of eggs were produced. A total of 47 and 46 full-sib progeny were randomly collected from the two full-sib families, respectively. The reference family including one parental pair and 380 offspring used for preliminary QTL analysis for growth traits was used in the fine mapping of QTL for growth traits.

Fin clips of the parents were collected and kept in absolute ethanol, whereas the whole body of each offspring at the age of 90 days post hatch (dph) was cut into small pieces, soaked in absolute ethanol, and kept in a -80° freezer. DNA was isolated and arrayed into 96-well PCR plates as described preciously [14].

Identification and genotyping of DNA markers

Partial genomic DNA libraries enriched for CA-, GA-, CAA-, GACA- and GATA- repeats were constructed as described previously [73]. Repeat-enriched DNA fragments of 400-1200 bp in length were cloned into pGEM-T vector (Promega, San Luis Obispo, CA), and transformed into XL-10 blue supercompetent cells (Stratagene, La Jolla, CA). The libraries were arrayed into 96- well plates for bidirectional sequencing on an ABI3730xl DNA sequencer (ABI, Foster City, CA) using the BigDye V3.0 kit and M13 and M13 reverse primers. Redundant and overlapping sequences were grouped using Sequencher (Gene Codes, Ann Arbor, MI). Unique sequences were compared to known microsatellite sequences of Asian seabass prior to primer design to remove redundancy. Microsatellites in ESTs were isolated using a method described previously [74]. Sequences containing CA>7, GA>7, CAA>6, GACA>5 and GATA>5 were subjected to primer design using PrimerSelect (DNASTAR, Brighton, MA), targeting a product size between 100 and 400 bp. Primers were designed for each unique sequence using PrimerSelect (DNASTAR, Madison, WI). One primer of each pair was labeled with FAM or HEX fluorescent dyes at the 5'-end. PCR for genotyping was carried out as described previously [14]. Products were analyzed using the DNA sequencer ABI3730xl, and genotyping was carried out to determine fragment size against the size standard GS-ROX-500 (Applied Biosystems, Foster City, CA) with software GeneMapper V3.5 (Applied Biosystems).

SNPs in genes were detected by PCR amplification of DNA of three parents and sequencing of PCR products. Briefly, genes of Asian seabass were aligned with genomic sequence data from zebrafish on GenBank. Primer sites in conserved exon regions were identified and primer pairs allowing PCR amplification of an intron-spanning fragment were developed. Amplified intron-spanning fragments which were sequenced as described above to detect SNP in the three parents from the reference families for linkage mapping. Ten SNPs were detected in 9 genes (Additional file 1). The genotyping for the 10 polymorphic SNP markers were performed by direct sequencing PCR products using ABI 3730xl Genetic Analyzer (Applied Biosystems). SNP genotypes were scored using Sequencher 4.9 (Genecodes).

Linkage Analysis

CRIMAP 2.4 software [75] was used for linkage analysis. Markers were placed into linkage groups based on twopoint LOD scores (>3.0) with markers from the previous map. Ordering the markers within a linkage group began with the previous map. A new marker was inserted into the map by evaluating the LOD at every possible location. The marker was then inserted at the location with the best LOD, and the change in length of the linkage group was evaluated. If a switched pair improved the LOD, the pair was switched and the process repeated until no better LOD could be found. The process of switching the order of marker pairs was repeated to finalize marker order. The final maps represent the most likely marker order identified with the complete data set. Once the most likely order had been derived, sex-average and sex-specific linkage distances were estimated for each using the Kosambi function. MapChart 2.2 software was used for graphical visualization of the linkage groups [76]. The total length of the linkage map was calculated by summing up the length of all 24 linkage groups.

Comparative genome analysis

As our previous study showed that Asian seabass is more closely related to T. nigroviridis than to zebrafish and medaka [61], comparison of flanking sequences of each Asian seabass marker on the map with the assembled genomic DNA sequences of the T. nigroviridis was conducted using BLAT http://www.genoscope.cns.fr/blat-server/cgi-bin/tetraodon/webBlat. BLAT searching was performed with a score above 80 and an alignment length of more than 50 bp as recommended [77].

QTL fine mapping

In our previous studies [39, 47], we identified significant QTL for growth traits on LGs (linkage groups) 2 and 3. To map the QTL on LGs 2 and 3 with more precision, 33 additional microsatellites located near the QTL on these two linkage groups were genotyped. With genotype data of the markers on the QTL regions and phenotypic data of the 380 progeny, QTL analysis was carried out using the program MapQTL 5.0 [78]. Interval mapping and multiple QTL model (MQM) mapping were utilized to detect any significant association between growth-related traits and marker loci in the data sets. Cofactors for MQM analyses were automatically selected with a p-value of 0.02. The LOD score significance thresholds were calculated by permutation tests in MapQTL 5.0, with a experiment-wise significance level of α < 0.05, n = 1000 for significant linkages. Calculation of the percentage of phenotypic variance explained (PVE) by a QTL was performed in MapQTL 5 on the basis of the population variance found within the progeny of the cross.

Potential candidate genes and their effects on growth traits

For those microsatellites, positions of which were detected well-overlapping with QTL peaks, we did a blast search of the sequences of microsatellites to identify potential candidate genes in these loci. Blast searches were done against the whole genome sequences of zebrafish and medaka in the current ENSEMBL release version http://www.ensembl.org and against all known sequences in GenBank.

A two-way ANOVA was performed on the 380 progeny using four allelic combinations (m1f1, m1f2, m2f1 and m2f2) from markers in the three candidate genes in order to investigate associations between phenotypic traits and genotypes of these genes. Mean phenotypic values of each trait were calculated for those progeny with the alternate alleles of the microsatellite markers, inherited from the male parent (m1 or m2), alleles inherited from the female parent (f1 or f2). This was conducted by using the general linear model (GLM) procedure of SAS (SAS Institute) and the Bonferroni method of multiple comparisons with α < 0.01.

Abbreviations

QTL:

Quantitative trait loci

SNP:

Single nucleotide polymorphism

cM:

Centimorgan

RAPD:

Rapid amplified polymorphic DNA

AFLP:

Amplified fragment length polymorphism.

References

  1. Falconer DS, Mackay TFC: Introduction to Quantitative Genetics. 1996, Harlow, Essex, UK: Addison Wesley Longman, Fourth

    Google Scholar 

  2. Paterson AH, DeVerna JW, Lanini B, Tanksley SD: Fine mapping of quantitative trait loci using selected overlapping recombinant chromosomes, in an interspecies cross of tomato. Genetics. 1990, 124 (3): 735-742.

    CAS  PubMed  PubMed Central  Google Scholar 

  3. Varshney RK, Tuberosa R, (eds): Genomics-Assisted Crop Improvement- Vol 2: Genomics Applications in Crops. 2007, Netherlands: Springer

  4. Lander ES, Green P, Abrahamson J, Barlow A, Daly MJ, Lincoln SE, Newburg L: MAPMAKER: An interactive computer package for constructing primary genetic linkage maps of experimental and natural populations. Genomics. 1987, 1: 174-181. 10.1016/0888-7543(87)90010-3.

    CAS  PubMed  Google Scholar 

  5. Tautz D: Hypervariabflity of simple sequences as a general source for polymorphic DNA markers. Nucl Acids Res. 1989, 17 (16): 6463-10.1093/nar/17.16.6463.

    CAS  PubMed  Google Scholar 

  6. Liu ZJ, Cordes JF: DNA marker technologies and their applications in aquaculture genetics. Aquaculture. 2004, 238 (1-4): 1-37. 10.1016/j.aquaculture.2004.05.027.

    CAS  Google Scholar 

  7. Gilbey J, Verspoor E, McLay A, Houlihan D: A microsatellite linkage map for Atlantic salmon (Salmo salar). Anim Genet. 2004, 35 (2): 98-105. 10.1111/j.1365-2052.2004.01091.x.

    CAS  PubMed  Google Scholar 

  8. Sakamoto T, Danzmann RG, Gharbi K, Howard P, Ozaki A, Khoo SK, Woram RA, Okamoto N, Ferguson MM, Holm LE, et al: A microsatellite linkage map of rainbow trout (Oncorhynchus mykiss) characterized by large sex-specific differences in recombination rates. Genetics. 2000, 155 (3): 1331-1345.

    CAS  PubMed  PubMed Central  Google Scholar 

  9. Guyomard R, Mauger S, Tabet-Canale K, Martineau S, Genet C, Krieg F, Quillet E: A type I and type II microsatellite linkage map of rainbow trout (Oncorhynchus mykiss) with presumptive coverage of all chromosome arms. BMC Genomics. 2006, 7: 302-10.1186/1471-2164-7-302.

    PubMed  PubMed Central  Google Scholar 

  10. Kucuktas H, Wang S, Li P, He C, Xu P, Sha Z, Liu H, Jiang Y, Baoprasertkul P, Somridhivej B, et al: Construction of genetic linkage maps and comparative genome analysis of catfish using gene-associated markers. Genetics. 2009, 181 (4): 1649-1660. 10.1534/genetics.108.098855.

    CAS  PubMed  PubMed Central  Google Scholar 

  11. 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 (1): 237-244. 10.1534/genetics.104.035022.

    CAS  PubMed  PubMed Central  Google Scholar 

  12. Xia JH, Feng L, Zhu ZY, Fu JJ, Feng JB, Li JL, Yue GH: A consensus linkage map of the grass carp (Ctenopharyngodon idella) based on microsatellites and SNPs. BMC Genomics. 2010, 11:

    Google Scholar 

  13. Sun X, Liang L: A genetic linkage map of common carp (Cyprinus carpio L.) and mapping of a locus associated with cold tolerance. Aquaculture. 2004, 238 (1-4): 165-172. 10.1016/S0044-8486(03)00445-9.

    CAS  Google Scholar 

  14. Wang CM, Zhu ZY, Lo LC, Feng F, Lin G, Yang WT, Li J, Yue GH: A microsatellite linkage map of Barramundi, Lates calcarifer. Genetics. 2007, 175 (2): 907-915. 10.1534/genetics.106.059972.

    CAS  PubMed  PubMed Central  Google Scholar 

  15. Chistiakov DA, Hellemans B, Haley CS, Law AS, Tsigenopoulos CS, Kotoulas G, Bertotto D, Libertini A, Volckaert FAM: A microsatellite linkage map of the European sea bass Dicentrarchus labrax L. Genetics. 2005, 170 (4): 1821-1826. 10.1534/genetics.104.039719.

    CAS  PubMed  PubMed Central  Google Scholar 

  16. Coimbra MRM, Kobayashi K, Koretsugu S, Hasegawa O, Ohara E, Ozaki A, Sakamoto T, Naruse K, Okamoto N: A genetic linkage map of the Japanese flounder, Paralichthys olivaceus. Aquaculture. 2003, 220 (1-4): 203-218. 10.1016/S0044-8486(02)00353-8.

    CAS  Google Scholar 

  17. Franch R, Louro B, Tsalavouta M, Chatziplis D, Tsigenopoulos CS, Sarropoulou E, Antonello J, Magoulas A, Mylonas CC, Babbucci M: A genetic linkage map of the hermaphrodite teleost fish Sparus aurata L. Genetics. 2006, 174 (2): 851-10.1534/genetics.106.059014.

    CAS  PubMed  PubMed Central  Google Scholar 

  18. Senger F, Priat C, Hitte C, Sarropoulou E, Franch R, Geisler R, Bargelloni L, Power D, Galibert F: The first radiation hybrid map of a perch-like fish: the gilthead seabream (Sparus aurata L). Genomics. 2006, 87 (6): 793-800. 10.1016/j.ygeno.2005.11.019.

    CAS  PubMed  Google Scholar 

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

    CAS  PubMed  Google Scholar 

  20. Shimoda N, Knapik EW, Ziniti J, Sim C, Yamada E, Kaplan S, Jackson D, de Sauvage F, Jacob H, Fishman MC: Zebrafish genetic map with 2000 microsatellite markers. Genomics. 1999, 58 (3): 219-232. 10.1006/geno.1999.5824.

    CAS  PubMed  Google Scholar 

  21. Groenen MA, Wahlberg P, Foglio M, Cheng HH, Megens HJ, Crooijmans RP, Besnier F, Lathrop M, Muir WM, Wong GK, et al: A high-density SNP-based linkage map of the chicken genome reveals sequence features correlated with recombination rate. Genome Res. 2009, 19 (3): 510-519. 10.1101/gr.086538.108.

    CAS  PubMed  PubMed Central  Google Scholar 

  22. Ihara N, Takasuga A, Mizoshita K, Takeda H, Sugimoto M, Mizoguchi Y, Hirano T, Itoh T, Watanabe T, Reed K: A comprehensive genetic map of the cattle genome based on 3802 microsatellites. Genome Research. 2004, 14 (10a): 1987-10.1101/gr.2741704.

    CAS  PubMed  PubMed Central  Google Scholar 

  23. Arias JA, Keehan M, Fisher P, Coppieters W, Spelman R: A high density linkage map of the bovine genome. BMC Genet. 2009, 10: 18-10.1186/1471-2156-10-18.

    PubMed  PubMed Central  Google Scholar 

  24. Vingborg RK, Gregersen VR, Zhan B, Panitz F, Hoj A, Sorensen KK, Madsen LB, Larsen K, Hornshoj H, Wang X, et al: A robust linkage map of the porcine autosomes based on gene-associated SNPs. BMC Genomics. 2009, 10: 134-10.1186/1471-2164-10-134.

    PubMed  PubMed Central  Google Scholar 

  25. Hwang TY, Sayama T, Takahashi M, Takada Y, Nakamoto Y, Funatsuki H, Hisano H, Sasamoto S, Sato S, Tabata S, et al: High-density integrated linkage map based on SSR markers in soybean. DNA Res. 2009, 16 (4): 213-225. 10.1093/dnares/dsp010.

    CAS  PubMed  PubMed Central  Google Scholar 

  26. Sato K, Nankaku N, Takeda K: A high-density transcript linkage map of barley derived from a single population. Heredity. 2009, 103 (2): 110-117. 10.1038/hdy.2009.57.

    CAS  PubMed  Google Scholar 

  27. Troggio M, Malacarne G, Coppola G, Segala C, Cartwright DA, Pindo M, Stefanini M, Mank R, Moroldo M, Morgante M, et al: A dense single-nucleotide polymorphism-based genetic linkage map of grapevine (Vitis vinifera L.) anchoring Pinot Noir bacterial artificial chromosome contigs. Genetics. 2007, 176 (4): 2637-2650. 10.1534/genetics.106.067462.

    CAS  PubMed  PubMed Central  Google Scholar 

  28. Reid D, Szanto A, Glebe B, Danzmann R, Ferguson M: QTL for body weight and condition factor in Atlantic salmon (Salmo salar): comparative analysis with rainbow trout (Oncorhynchus mykiss) and Arctic charr (Salvelinus alpinus). Heredity. 2004, 94 (2): 166-172. 10.1038/sj.hdy.6800590.

    Google Scholar 

  29. Cnaani A, Hallerman EM, Ron M, Weller JI, Indelman M, Kashi Y, Gall GAE, Hulata G: Detection of a chromosomal region with two quantitative trait loci, affecting cold tolerance and fish size, in an F-2 tilapia hybrid. Aquaculture. 2003, 223 (1-4): 117-128. 10.1016/S0044-8486(03)00163-7.

    CAS  Google Scholar 

  30. Fuji K, Kobayashi K, Hasegawa O, Coimbra M, Sakamoto T, Okamoto N: Identification of a single major genetic locus controlling the resistance to lymphocystis disease in Japanese flounder (Paralichthys olivaceus). Aquaculture. 2006, 254 (1-4): 203-210. 10.1016/j.aquaculture.2005.11.024.

    CAS  Google Scholar 

  31. Ozaki A, Sakamoto T, Khoo S, Nakamura K, Coimbra M, Akutsu T, Okamoto N: Quantitative trait loci (QTLs) associated with resistance/susceptibility to infectious pancreatic necrosis virus (IPNV) in rainbow trout (Oncorhynchus mykiss). Mol Genet Genom. 2001, 265 (1): 23-31. 10.1007/s004380000392.

    CAS  Google Scholar 

  32. Massault C, Hellemans B, Louro B, Batargias C, Van Houdt JKJ, Canario A, Volckaert FAM, Bovenhuis H, Haley C, De Koning DJ: QTL for body weight, morphometric traits and stress response in European sea bass Dicentrarchus labrax. Anim Genet. 2010, 41: 337-345.

    CAS  PubMed  Google Scholar 

  33. Wang S, Zhang L, Meyer E, Matz MV: Construction of a high-resolution genetic linkage map and comparative genome analysis for the reef-building coral Acropora millepora. Genom Biol. 2009, 10 (11): R126-10.1186/gb-2009-10-11-r126.

    Google Scholar 

  34. Criscione CD, Valentim CL, Hirai H, LoVerde PT, Anderson TJ: Genomic linkage map of the human blood fluke Schistosoma mansoni. Genom Biol. 2009, 10 (6): R71-10.1186/gb-2009-10-6-r71.

    Google Scholar 

  35. Yamamoto K, Nohata J, Kadono-Okuda K, Narukawa J, Sasanuma M, Sasanuma S, Minami H, Shimomura M, Suetsugu Y, Banno Y, et al: A BAC-based integrated linkage map of the silkworm Bombyx mori. Genom Biol. 2008, 9 (1): R21-10.1186/gb-2008-9-1-r21.

    Google Scholar 

  36. Gross JB, Protas M, Conrad M, Scheid PE, Vidal O, Jeffery WR, Borowsky R, Tabin CJ: Synteny and candidate gene prediction using an anchored linkage map of Astyanax mexicanus. Proc Natl Acad Sci USA. 2008, 105 (51): 20106-20111. 10.1073/pnas.0806238105.

    CAS  PubMed  Google Scholar 

  37. Hillier LW, Miller RD, Baird SE, Chinwalla A, Fulton LA, Koboldt DC, Waterston RH: Comparison of C. elegans and C. briggsae genome sequences reveals extensive conservation of chromosome organization and synteny. PLoS Biol. 2007, 5 (7): e167-10.1371/journal.pbio.0050167.

    PubMed  PubMed Central  Google Scholar 

  38. FAO: Food and Agriculture Organization of the United Nations Statistics. 2008

    Google Scholar 

  39. Wang CM, Lo LC, Zhu ZY, Yue GH: A genome scan for quantitative trait loci affecting growth-related traits in an F1 family of Asian seabass (Lates calcarifer). BMC Genomics. 2006, 7: 274-10.1186/1471-2164-7-274.

    PubMed  PubMed Central  Google Scholar 

  40. Yue GH, Li Y, Chao TM, Chou R, Orban L: Novel microsatellites from Asian seabass (Lates calcarifer) and their application to broodstock analysis. Mar Biotechnol. 2002, 4 (5): 503-511.

    CAS  PubMed  Google Scholar 

  41. Zhu ZY, Lin G, LC Lo, YX Xu, Renee C, et al: Genetic analyses of Asian seabass stocks using novel polymorphic microsatellites. Aquaculture. 2006, 256: 167-173. 10.1016/j.aquaculture.2006.02.033.

    CAS  Google Scholar 

  42. Zhu ZY, Wang CM, Lo LC, Feng F, Lin G, Yue GH: Isolation, characterization, and linkage analyses of 74 novel microsatellites in Barramundi (Lates calcarifer). Genome. 2006, 49 (8): 969-976. 10.1139/G06-080.

    CAS  PubMed  Google Scholar 

  43. Xu YX, Zhu ZY, Lo LC, Wang CM, Lin G, Feng F, Yue GH: Characterization of two parvalbumin genes and their association with growth traits in Asian seabass (Lates calcarifer). Anim Genet. 2006, 37: 266-268. 10.1111/j.1365-2052.2006.01423.x.

    CAS  PubMed  Google Scholar 

  44. Xia JH, He XP, Bai ZY, Yue GH: Identification and characterization of 63 microRNAs in the Asian seabass Lates calcarifer. Plos ONE. 2011, 6: e17537-10.1371/journal.pone.0017537.

    CAS  PubMed  PubMed Central  Google Scholar 

  45. 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.

    PubMed  PubMed Central  Google Scholar 

  46. Xia JH, Yue GH: Identification and analysis of immune-related transcriptome in Asian seabass Lates calcarifer. BMC Genomics. 2010, 11: 356-10.1186/1471-2164-11-356.

    PubMed  PubMed Central  Google Scholar 

  47. Wang C, Lo L, Feng F, Zhu Z, Yue G: Identification and verification of QTL associated with growth traits in two genetic backgrounds of Barramundi (Lates calcarifer). Anim Genet. 2008, 39 (1): 34-39. 10.1111/j.1365-2052.2007.01672.x.

    PubMed  Google Scholar 

  48. Harushima Y, Yano M, Shomura A, Sato M, Shimano T, Kuboki Y, Yamamoto T, Lin SY, Antonio BA, Parco A, et al: A high-density rice genetic linkage map with 2275 markers using a single F2 population. Genetics. 1998, 148 (1): 479-494.

    CAS  PubMed  PubMed Central  Google Scholar 

  49. Jensen-Seaman MI, Furey TS, Payseur BA, Lu Y, Roskin KM, Chen CF, Thomas MA, Haussler D, Jacob HJ: Comparative recombination rates in the rat, mouse, and human genomes. Genome Res. 2004, 14 (4): 528-538. 10.1101/gr.1970304.

    CAS  PubMed  PubMed Central  Google Scholar 

  50. Hardie DC, Hebert PDN: Genome-size evolution in fishes. Can J Fish Aqua Sci. 2004, 61: 1636-1646. 10.1139/f04-106.

    Google Scholar 

  51. Woram RA, McGowan C, Stout JA, Gharbi K, Ferguson MM, Hoyheim B, Davidson EA, Davidson WS, Rexroad C, Danzmann RG: A genetic linkage map for Arctic char (Salvelinus alpinus): evidence for higher recombination rates and segregation distortion in hybrid versus pure strain mapping parents. Genome. 2004, 47 (2): 304-315. 10.1139/g03-127.

    CAS  PubMed  Google Scholar 

  52. Dib C, Faure S, Fizames C, Samson D, Drouot N, Vignal A, Millasseau P, Marc S, Hazan J, Seboun E, et al: A comprehensive genetic map of the human genome based on 5,264 microsatellites. Nature. 1996, 380 (6570): 152-154. 10.1038/380152a0.

    CAS  PubMed  Google Scholar 

  53. Mott R, Talbot C, Turri M, Collins A, Flint J: A method for fine mapping quantitative trait loci in outbred animal stocks. Proc Natl Acad Sci USA. 2000, 97 (23): 12649-10.1073/pnas.230304397.

    PubMed  Google Scholar 

  54. Liu ZJ, (ed): Aquaculture Genome Technologies. 2007, Wiley-Blackwell

  55. Haidle L, Janssen JE, Gharbi K, Moghadam HK, Ferguson MM, Danzmann RG: Determination of quantitative trait loci (QTL) for early maturation in rainbow trout (Oncorhynchus mykiss). Mar Biotechnol. 2008, 10 (5): 579-592. 10.1007/s10126-008-9098-5.

    CAS  PubMed  PubMed Central  Google Scholar 

  56. Moen T, Baranski M, Sonesson AK, Kjoglum S: Confirmation and fine-mapping of a major QTL for resistance to infectious pancreatic necrosis in Atlantic salmon (Salmo salar): population-level associations between markers and trait. BMC Genomics. 2009, 10: 368-10.1186/1471-2164-10-368.

    PubMed  PubMed Central  Google Scholar 

  57. Aparicio S, Chapman J, Stupka E, Putnam N, Chia J, Dehal P, Christoffels A, Rash S, Hoon S, Smit A: Whole-genome shotgun assembly and analysis of the genome of Fugu rubripes. Science. 2002, 297 (5585): 1301-10.1126/science.1072104.

    CAS  PubMed  Google Scholar 

  58. Kitano J, Ross J, Mori S, Kume M, Jones F, Chan Y, Absher D, Grimwood J, Schmutz J, Myers R: A role for a neo-sex chromosome in stickleback speciation. Nature. 2009, 461 (7267): 1079-1083. 10.1038/nature08441.

    CAS  PubMed  PubMed Central  Google Scholar 

  59. Jaillon O, Aury J, Brunet F, Petit J, Stange-Thomann N, Mauceli E, Bouneau L, Fischer C, Ozouf-Costaz C, Bernot A: Genome duplication in the teleost fish Tetraodon nigroviridis reveals the early vertebrate proto-karyotype. Nature. 2004, 431 (7011): 946-957. 10.1038/nature03025.

    PubMed  Google Scholar 

  60. Kasahara M, Naruse K, Sasaki S, Nakatani Y, Qu W, Ahsan B, Yamada T, Nagayasu Y, Doi K, Kasai Y: The medaka draft genome and insights into vertebrate genome evolution. Nature. 2007, 447 (7145): 714-719. 10.1038/nature05846.

    CAS  PubMed  Google Scholar 

  61. Lin G, Lo LC, Zhu ZY, Feng F, Chou R, Yue GH: The complete mitochondrial genome sequence and characterization of single-nucleotide polymorphisms in the control region of the Asian seabass (Lates calcarifer). Mar Biotechnol. 2006, 8 (1): 71-79. 10.1007/s10126-005-5051-z.

    CAS  PubMed  PubMed Central  Google Scholar 

  62. Rothschild M, Jacobson C, Vaske D, Tuggle C, Wang L, Short T, Eckardt G, Sasaki S, Vincent A, McLaren D, et al: The estrogen receptor locus is associated with a major gene influencing litter size in pigs. Proc Natl Acad Sci USA. 1996, 93 (1): 201-205. 10.1073/pnas.93.1.201.

    CAS  PubMed  Google Scholar 

  63. Doust AN, Devos KM, Gadberry MD, Gale MD, Kellogg EA: Genetic control of branching in foxtail millet. Proc Natl Acad Sci USA. 2004, 101 (24): 9045-9050. 10.1073/pnas.0402892101.

    CAS  PubMed  Google Scholar 

  64. Albayrak Ö, Tirniceriu A, Riemenschneider M, Kurz A, Scherag A, Egensperger R: The Cathepsin D (224C/T) polymorphism confers an increased risk to develop Alzheimer's disease in men. J Geront: Series A. 2010, 65 (3): 219-10.1093/gerona/glp209.

    CAS  Google Scholar 

  65. Kontorovich T, Levy A, Korostishevsky M, Nir U, Friedman E: Single nucleotide polymorphisms in miRNA binding sites and miRNA genes as breast/ovarian cancer risk modifiers in Jewish high-risk women. Inter J Cancer. 2009, 127 (3): 589-597. 10.1002/ijc.25065.

    Google Scholar 

  66. Strausberg RL, Feingold EA, Grouse LH, Derge JG, Klausner RD, Collins FS, Wagner L, Shenmen CM, Schuler GD, Altschul SF, et al: Generation and initial analysis of more than 15,000 full-length human and mouse cDNA sequences. Proc Natl Acad Sci USA. 2002, 99 (26): 16899-16903. 10.1073/pnas.242603899.

    PubMed  Google Scholar 

  67. Lau WL, Scholnick SB: Identification of two new members of the CSMD gene family. Genomics. 2003, 82 (3): 412-415. 10.1016/S0888-7543(03)00149-6.

    CAS  PubMed  Google Scholar 

  68. Schuster SC: Next-generation sequencing transforms today's biology. Nature. 2007, 200: 8-

    Google Scholar 

  69. Carlson CS, Eberle MA, Kruglyak L, Nickerson DA: Mapping complex disease loci in whole-genome association studies. Nature. 2004, 429 (6990): 446-452. 10.1038/nature02623.

    CAS  PubMed  Google Scholar 

  70. Goddard M, Hayes B: Genomic selection. J Anim Breed Genet. 2007, 124 (6): 323-330. 10.1111/j.1439-0388.2007.00702.x.

    CAS  PubMed  Google Scholar 

  71. Lazzaro BP, Sceurman BK, Clark AG: Genetic basis of natural variation in D. melanogaster antibacterial immunity. Science. 2004, 303 (5665): 1873-10.1126/science.1092447.

    CAS  PubMed  Google Scholar 

  72. Barton NH, Keightley PD: Understanding quantitative genetic variation. Nat Review Genet. 2002, 3 (1): 11-21. 10.1038/nrg700.

    CAS  Google Scholar 

  73. Yue GH, Chen F, Orban L: Rapid isolation and characterization of microsatellites from the genome of Asian arowana (Scleropages formosus, Osteoglossidae, Pisces). Mol Ecol. 2000, 9 (7): 1007-1009. 10.1046/j.1365-294x.2000.00939-10.x.

    CAS  PubMed  Google Scholar 

  74. Yue GH, Zhu ZY, Wang CM, Xia JH: A simple and efficient method for isolating polymorphic microsatellites from cDNA. BMC Genomics. 2009, 10: 125-10.1186/1471-2164-10-125.

    PubMed  PubMed Central  Google Scholar 

  75. Green P, Falls K, Crooks S: Documentation for CRI-MAP version 2.4. Washington University School of Medicine. 1990, St. Louis, Mo., USA

    Google Scholar 

  76. 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.

    CAS  PubMed  Google Scholar 

  77. Kent WJ: BLAT--the BLAST-like alignment tool. Genom Res. 2002, 12 (4): 656-

    CAS  Google Scholar 

  78. Van Ooijen JW: MapQTL 5, Software for the Mapping of Quantitative Trait Loci in Experimental Populations. 2004, Kyazma B.V., Wageningen, The Netherlands

    Google Scholar 

Download references

Acknowledgements

This study is funded by Ministry of National Development for the project "To ensure self-sufficiency of safe seafood for Singapore through the development of aquaculture genomic tools for marker-assisted selective breeding of tropical marine foodfish". We thank staffs of Marine Aquaculture Center, AVA for technical supports.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Gen Hua Yue.

Additional information

Authors' contributions

GHY initiated and overviewed the Asian seabass project and finalized the manuscript. WCM and YGH designed the study. WCM, BZY, HXP, XJH, SF, LLC, FF and ZZY conducted the experiments. WCM conducted the data analysis and wrote the manuscript. All authors have read and approved the final manuscript.

Electronic supplementary material

Additional file 1:790 DNA markers mapped to the second generation linkage map of Asian seabass. (XLS 214 KB)

12864_2010_10042_MOESM2_ESM.DOC

Additional file 2:Comparative mapping of markers of Lates calcarifer through BLAT search against Tetraodon nigroviridis genome. (DOC 132 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

Authors’ original file for figure 3

Authors’ original file for figure 4

Authors’ original file for figure 5

Authors’ original file for figure 6

Authors’ original file for figure 7

Authors’ original file for figure 8

Authors’ original file for figure 9

Authors’ original file for figure 10

Authors’ original file for figure 11

Authors’ original file for figure 12

Authors’ original file for figure 13

Authors’ original file for figure 14

Authors’ original file for figure 15

Authors’ original file for figure 16

Authors’ original file for figure 17

Rights and permissions

This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Cite this article

Wang, C.M., Bai, Z.Y., He, X.P. et al. A high-resolution linkage map for comparative genome analysis and QTL fine mapping in Asian seabass, Lates calcarifer. BMC Genomics 12, 174 (2011). https://doi.org/10.1186/1471-2164-12-174

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-12-174

Keywords