- Research article
- Open Access
Whole genome scanning and association mapping identified a significant association between growth and a SNP in the IFABP-a gene of the Asian seabass
BMC Genomics volume 14, Article number: 295 (2013)
Aquaculture is the quickest growing sector in agriculture. However, QTL for important traits have been only identified in a few aquaculture species. We conducted QTL mapping for growth traits in an Asian seabass F2 family with 359 individuals using 123 microsatellites and 22 SNPs, and performed association mapping in four populations with 881 individuals.
Twelve and nine significant QTL, as well as 14 and 10 suggestive QTL were detected for growth traits at six and nine months post hatch, respectively. These QTL explained 0.9-12.0% of the phenotypic variance. For body weight, two QTL intervals at two stages were overlapped while the others were mapped onto different positions. The IFABP-a gene located in a significant QTL interval for growth on LG5 was cloned and characterized. A SNP in exon 3 of the gene was significantly associated with growth traits in different populations.
The results of QTL mapping for growth traits suggest that growth at different stages was controlled by some common QTL and some different QTL. Positional candidate genes and association mapping suggest that the IFABP-a is a strong candidate gene for growth. Our data supply a basis for fine mapping QTL, marker-assisted selection and further detailed analysis of the functions of the IFABP-a gene in fish growth.
Growth traits are the major interest in aquaculture. Growth traits are complex traits that regulated by many loci that have quantitative effects on the phenotype. Due to the ease of measurement, the economic importance and intermediate heritabilities for growth traits, phenotypic selection has been commonly and effectively applied in fish breeding programs. However, the process of developing new varieties or lines is labor intensive, time consuming and costly by using traditional breeding methods. Marker-assisted selection (MAS) by using markers tightly linked to quantitative trait loci (QTL) as a substitute to assist phenotypic screening, can overcome the shortcomings in traditional breeding, thus increasing the accuracy and efficiency of selection . This technique is especially valuable for traits with low to moderate heritability, which are difficult to be improved by traditional selection.
Identification of QTL related to growth is the basis for application of genetic markers in fish breeding. A key step in QTL analysis of complex traits is the establishment of large collections of genetic markers . A few studies are available on linkage and QTL mapping for growth traits in fish, such as salmonids [3–6], sea bream [7–11], rainbow trout , tilapia [13, 14], Asian and European seabass [15–20] and turbot . However, the identification of QTL underlying complex traits have proved to be very challenging due to gene-gene and gene-environment interactions, low proportion of explained phenotypic variation for each locus [22, 23], the large number of QTL, and many additional sources of variation . Linkage-based QTL mapping and linkage disequilibrium (LD)-based genetic association mapping are two major approaches to locate genes that control phenotypes of interest. Combined linkage and association analysis showed significant improvement in detection of QTL and have shown the ability to unravel the secrets of complex traits in animal breeding . Few such reports exist for fish to date.
Many traits with importance in agriculture change with time or other independent variables . In animal, many factors such as skeletal structure, organ, muscle, and fat mass contribute to overall body size  and family-based, age-related and strain-specific growth differences are likely to co-exist [28–30]. Genetic factors affecting growth are expected to vary among individuals, families and different stages of development. The power of QTL analysis is therefore seriously limited, since only segregating QTL in one or both parents that were used in constructing mapping families can be detected . Although some family/population-specific QTL have been reported in the Asian seabass [15–17] and other fishes [28–30], the molecular and genetic basis of growth traits remains fragmentary, highlighting the need for additional studies of these important trait in fish.
The Asian seabass Lates calcarifer is an important farmed marine foodfish species, widely distributed in the Indo-West Pacific region (FishBase: http://www.fishbase.org). Currently a second generation linkage map for the Asian seabass was developed with 790 markers . Based on the genetic map information, some QTL for growth traits were mapped [15–17]. Significant QTL for growth traits measured at the age of 3 months were mapped onto linkage groups 2 and 3 [16, 17] using F1 families. Fine mapping of QTL on linkage groups 2 and 3 identified a QTL cluster underlying growth traits in the Asian seabass . These studies provided preliminary data for application of MAS selection in the Asian seabass. However, due to the low analytic power of F1 families for QTL mapping, some QTL for growth might be missed.
The potentials of MAS to improve traits are increasingly evident in aquaculture. Identification of QTL-related economic traits is a major aim of fish genome research. In Southeast Asia, the harvesting time for the aquacultured Asian seabass is starting from about nine months of age. In order to map QTL for growth traits at different farming stages of the species, we conducted a whole genome scan for growth traits at the age of six and nine months post hatch by using an F2 family including 359 individuals, and genotyping with 123 microsatellites and 22 SNPs. We identified significant and suggestive QTL that affect growth-related traits (body weight, standard length and total length measured at the age of 6 and 9 months) in the Asian seabass. In addition, we detected SNPs in 19 candidate genes related to growth, nutritional condition or environmental adaptation and mapped them in the linkage map of Asian seabass. The gene IFABP-a located on LG5 was significantly associated with growth. Our study provides a basis for identification of genes controlling growth traits and conducting of MAS in the Asian seabass.
Results and discussion
Phenotypic data in the F2mapping family
In the Asian seabass, QTL mapping for growth traits that measured at the age of 3 months have been reported [15–17]. Significant QTL for growth traits were mapped onto linkage groups 2 and 3 [16, 17]. In these previous studies, only F1 families were used in QTL mapping. Due to the low power of F1 families for QTL mapping, some QTL might not be detected. To identify QTL that have not been identified in previous studies, we carried out a QTL analysis for growth traits that measured at the age of 6 and 9 months in an F2 full sib mapping family (F2S2; N = 359) of the Asian seabass. Body weight, standard length and total length data (named as BW6M, BW9M, SL6M, SL9M, TL6M and TL9M, respectively) were measured at two different stages of development. The means (with standard deviation, sd) for the 6 growth traits in the F2 population are provided in Table 1. Significant phenotypic differences were observed among individuals of the family. The observed phenotypic values ranged from 44 g to 442 g (~10 fold differences) with an average of 227.6 (± 62.0 sd) g for the trait BW6M, and the values ranged from 117 g to 672 g (~6 fold differences) for the trait BW9M. The phenotypic variance within a population is the result of genetic sources and/or environmental sources. Under the same conditions, a larger phenotypic variance in a population suggests a higher degree of genetic variation. The population F2S2 showed substantial levels of phenotypic variation, therefore, had enough statistical power for QTL detection. Strong positive phenotypic correlations between pairwise traits were observed. The coefficients ranged from 0.73 between BW9M and SL6M to 0.96 between TL9M and SL9M, underlying a common genetic basis for these traits. Similar results were reported in tilapia  and turbot .
Linkage and QTL mapping
One hundred and twenty-three microsatellite markers showed informativeness in the F2 mapping family. These markers distributes evenly on the 24 linkage groups that defined by Wang et al. . In addition, 22 SNP markers that developed from genes related to growth, nutritional condition, or environmental adaptation were also applied in map construction. Primer information for the SNP markers was listed in Additional file 1: Table S1. Using linkage analysis and LOD ≥ 3 as threshold for mapping data, all loci were finally assigned to 24 linkage groups (LGs), designated as LG1 – LG24, corresponding to the linkage group names of the high-resolution linkage maps [15, 16]. The largest linkage group LG17 obtained in this study containing 8 markers and had a length of 93.5 cM, and the shortest linkage group LG20 consisted of 5 markers with a length of 2.3 cM. The total length for the map was 812 cM with an average linkage group length of 33.8 (± 22.1 sd) cM. The average spacing between adjacent markers was 7.8 (± 5.5 sd) cM, therefore meeting the minimum requirements for QTL analysis . Most markers mapped in the same linear order as in the high-resolution linkage maps [15, 16]. The differences between two maps could be caused from the differences in mapping population sizes and marker numbers, genotyping errors and the nature of the population, e.g., differences in family structures and female/male specific recombination rates. Summary of linkage map information and the averaged linkage map are presented in Additional file 2: Table S2 and Additional file 3: Figure S1, respectively.
We analyzed six types of growth traits. Genome-wide (GW) and linkage group (LG)-wide LOD threshold (p < 0.01 and 0.05) were calculated from 10000 permutations of the quantitative trait data. Forty-five suggestive QTL including 21 significant QTL (GW threshold: p < 0.05) were detected across all traits, which might be representative of the genetic architecture of growth traits in the Asian seabass. Two QTL intervals were detected on each of the 4 linkage groups (LG5, 9, 18 and 24, respectively). No significant association was found in 11 linkage groups (LG1, 4, 8, 10, 12, 14, 16, 17, 19, 20 and 22), probably due to the low informativeness of involved markers and the nature of the genetic architectures of these traits. For each trait, 5 (BW9M) to 9 (TL6M and SL6M) suggestive QTL were detected, suggesting that these traits are highly polygenic. One QTL confidence interval on linkage group LG15 for SL6M and SL9M was large, covering nearly the entire chromosomes. Of the 21 significant QTL (GW threshold: p < 0.05), 15 were detected at the genome-wide significant level (GW threshold: p < 0.01). QTL that surpassed the suggestive LOD threshold are summarized in Table 1 and shown in Figure 1.
The proportions of phenotypic variance explained (PVE) by all the suggestive QTL were generally small. The highest PVE was 12.0% for the suggestive QTL tl6m_i, which were located on LG24 (8.07 cM). The PVE for the eight QTL of the trait BW6M ranged from 2.5% (bw6m_e) to 7.2% (bw6m_b); and for the nine QTL of trait SL6M, the PVE ranged from 3% (sl6m_e) to 10.3% (sl6m_i). The small PVEs of the QTL that identified per trait suggest that the expression of these traits is under the control of a number of genes. Similar results were obtained in turbot by Sánchez-Molano , in which, the PVE of the tightest linked markers to the significant QTL varied from 8.0% to 13.6% for weight and from 8.2% to 12.9% for body length. Some potential candidate genes were detected near or on some of the QTL intervals, e.g. IG2F2BP3 on LG3, IFABP-a on LG5, 2-TNNL2 on LG6, GCDH on LG11, SP on LG13 and APOCI on LG15.
Significant differences in QTL positions for different growth traits were also found, such as sl6m_a on LG3, sl9m_b on LG5 and sl9m_c on LG6, suggesting that different sets of genes are involved in the control of these traits. In a previous study of the Asian seabass [15, 16], five significant QTL controlling growths were detected on LG2 and LG3, and explained 6.4 - 59.7% of the observed phenotypic variance, respectively. We also identified two QTL regions on these linkage groups. But, low PVEs and only suggestive significance level (p < 0.05) were found for these QTL. These differences could result from the growth variations at different development stages among different families, since we measured traits of fishes at the age of 6 and 9 months, but those used by Wang et al. [15–17] were measured at the 3 months. Specific QTL for body weight, standard length and total length that detected at different sampling time period suggest that there were dynamic expression for the traits, as in salmonids , in which some QTL were differentially expressed at different developing stages.
Co-localization of the detected QTL among traits is evident, indicating that pleiotropic effects may be playing a role in growth of the Asian seabass. This result was expected since high genetic correlations among growth traits were found in the seabass and other fish species. All suggestive QTL were located in 17 different genomic regions of 13 linkage groups. QTL regions having similar LOD profiles for multiple traits were detected on 10 linkage groups (LG2, 5, 7, 9, 11, 13, 15, 18, 21 and 24). Interestingly, one QTL region on linkage group LG5 (0~1 cM) was significantly associated to all traits (GW threshold: p < 0.01). Given that the high and sharp LOD peak, one common locus might be responsible for the different traits. The PVE of traits for this QTL ranged from 5.4% (tl9m_b) to 10.5% (sl9m_b). Near or on this QTL, eight markers including two SNP markers were located within 1.2 cM regions (Additional file 3: Figure S1 and Figure 1). Interestingly we found IFABP-a-SNP1245 marker showed a significant association to the traits. This result suggests that IFABP-a gene is of great importance in growth of fish and hence, warrant further studies.
Association of heterozygosity with growth traits in the full sib mapping family
The degree of heterozygosity was used to test for direct effects of the markers on the growth traits. To compare the distribution of genotypes between pairs of extreme size classes for each trait (N = 30 for each class), the degree of heterozygosity for used loci in two extreme size groups of each trait were calculated and tested by using Fisher's exact test. Pairwise tests revealed significant heterogeneity of genotypes among two extreme size groups in multiple markers (Additional file 4: Table S3). The degrees of heterozygosity for 4 loci (SL9M) to 12 loci (SL6M, TL6M) were significantly different between two extreme size groups of each trait. Three SNP markers (IFABP-a-SNP1245, ECHS1, APOE-SNP1) showed significant heterogeneity of genotypes among size-group pairs. The heterozygosities for four of the six loci for BW6M and five of the eight loci for BW9M were highly present in the large size group (80-97% and 37-93%, respectively). The remaining loci show high heterozygosities in the small size group. IFABP-a-SNP1245 show high heterozygosity (CG genotype) in the large size group (77-90%) and show low heterozygosity in the small size group (20-30%) for all six traits (p < 0.001) (Figure 2). This result confirmed the putative association of IFABP-a-SNP1245 with growth traits that obtained in QTL analysis. The heterozygosities of marker genotypes that significantly associated with positive effects in traits suggest possible heterosis and potential utility for MAS in the Asian seabass. For example, selection with heterozygote would favor growth in the population. The significant associations between the degree of heterozygosity of loci and traits are presented in Additional file 4: Table S3.
Analysis of associations between candidate markers and growth traits in 4 populations
The large impact of markers/genes detected in the QTL model does not necessarily mean that the polymorphisms are causative mutations. LD-based association mapping represents a much greater precision than QTL mapping . To verify the association of the markers to body weight that were identified by QTL mapping and display their applications in other families in breeding programs, six markers on or near QTL regions were selected and applied in association analysis of 4 populations produced by mass crosses. These markers were distributed on 5 linkage groups. Two SNP markers (IFABP-a-SNP1245 and IFABP-a-SNP1550) were located on the same gene of linkage group LG5, the other 4 markers located on separate linkage groups. Except the microsatellite marker Lca949, all the remaining five markers were located on QTL intervals at genome wide significance level. For each population, 142 to 248 individuals were used. Substantial phenotypic variation among populations ranging from 40.2 g (F2S5 and F2S6) to 65.2 g (F2S1) was observed at the age of 3 months (BW3M). This variation can result from uncontrolled environmental factors; each population was therefore treated separately during association analysis. Summary description of the association study in the four populations using selected markers is presented in Table 2.
The LD association was tested between the homozygous/heterozygous genotypes and haplotypes of the selected markers and phenotype data (BW3M). For each selected marker, an association analysis was performed. An important concordance was observed between QTL mapping and association studies. Except the marker Lca949, the remaining five markers show significant association with body weight trait in one to two populations. Five markers show significant association in population F2S5. The remaining three populations were associated to one or two markers under different models. This may result from the limited population sizes, age differences, and some genotypic combinations of the two loci are unfavorable under some specific conditions. In the four populations, the observed phenotypic data for the heterozygosity genotype (CG) for marker IFABP-a-SNP1245 were generally larger than those for CC genotypes (Figure 3), indicating the SNP has a significant additive effect on BW3M. IFABP-a-SNP1245 haplotypes and genotypes were significantly associated to BW3M in two populations (F2S5 and F2S6). However, to our surprise, when compared the trait data for three genotypes (CC, CG and GG) in F2S6 and F2S7, the GG genotype showed negative effect in F2S6, but showed positive effect in F2S7. These results indicated there were unfavorable factors that associated to G haplotype in F2S6 populations, which also compromised the CG heterozygosity effect in the population. Our study showed that the markers located on the QTL intervals at the genome wide significance level were quite reliable for application in MAS. Since family/population-specific QTL exist, further studies using more markers near QTL and a larger number of families with high performances are required to identify marker-trait associations.
Characterization of IFABP genes in the Asian seabass
Both QTL and association mapping revealed the IFABP-a SNP1245 marker was significantly associated to growth traits, which explained 5.4-10.5% of phenotypic variance. Therefore, IFABP-a may be an important candidate gene for growth. To further validate its functions in growth, we studied the IFABP gene by cloning and characterization. In fish, IFABP gene was characterized only in few species, e.g., Paralichthys olivaceus and zebrafish [35–39]. The IFABP gene shows binding specificity for long-chain fatty acids and is proposed to be involved in uptake of dietary fatty acids and their intracellular transport , similar to the roles in mammals. We found that IFABP gene has two isoforms, IFABP-a and IFABP-b, in the Asian seabass. These two isoforms locate on LG5 and LG14, respectively. We cloned the complete ORFs and partial UTRs for both genes. The two isoforms showed similar gene structures, including 4 exons and 3 introns with a deduced polypeptide sequence of 132 amino acids in each gene, similar to that of P. olivaceus. However, the two genes are far apart with a genetic distance of 0.278 between the amino acid sequences of two isoforms. The phylogenetic tree analysis suggests that the two seabass IFABP protein genes belong to two distinct groups (Additional file 5: Figure S2). The patterns of diversification suggest a long evolutionary history for the two isoforms. The SNP mutation IFABP-a-SNP1245 showing significant association to growth was located at exon 3 of the gene (Figure 4). In the QTL mapping and association populations, two alleles (C and G) were detected at SNP1245. This C to G transversion caused the amino acid Threonine (Thr/T) changing to Arginine (Arg/R) in the polypeptide sequence. This mutation might have influence on growth in fish. Another SNP IFABP-a-SNP1550 was mutated at the 3rd intron between exon 3 and exon 4 of the gene without significant association to growth (Figure 4).
Expressions of the two IFABP homologs in the fish intestines in response to fast challenge
qPCR revealed that both isoforms were abundantly expressed in the intestines but weakly or barely detectable in other tissues (P < 0.001; Figure 5-a). IFABP-a was more highly expressed than IFABP-b in the control intestines (~7 fold). Temporal expression of the two genes in the intestines showed substantial down-regulation at 3, 6, 12 days post fast (5.9-8.8 fold for IFABP-a, and 3.2-10.8 fold for IFABP-b) (Figure 5-b), underlining the functional importance in nutritional digestion and absorption-related pathways of the fish. The IFABP-a gene shows significant down-regulation in samples at 6 days post fast (P < 0.001; Figure 5-b), and the IFABP-b shows significant down-regulation in samples at 12 days post fast (P < 0.001; Figure 5-b), suggesting different functions for these isoforms. In mammals, significant associations among fatness, the abundance of IFABP and growth [40, 41] and common variants showing an important effect on fatness and growth [42–44] were reported. These studies suggested that the growth and fatness shared partially genetic architecture. In our study, QTL mapping and association analysis suggest that IFABP-a gene might regulate the growth in fish. Although direct functional evidence relating the variant to growth was lacking, analyses of IFABP in the Asian seabass by qPCR gave additional support to our QTL findings. Since IFABP gene plays key roles in metabolism though the digestion and absorption of fatty acids in animals, it is reasonable to consider that IFABP-a may be a strong candidate for the important growth genes in fish.
A genome-wide QTL study was performed for growth traits at age of 6 and 9 months post hatch in the Asian seabass. Twenty-one significant QTL and 24 suggestive QTL were detected across all traits, which might be representative of the genetic architecture of growth traits in the Asian seabass. QTL mapping and association analysis revealed that IFABP-a-SNP1245 was significantly associated to growth. Characterization of IFABP gene gives additional support to our QTL findings and suggests that IFABP-a is strong candidate gene for growth in fish. Our data supply a basis for fine mapping QTL, marker-assisted selection and further detailed analysis of the functions of the IFABP gene in fish growth.
Constructing mapping populations and phenotyping
To construct an F2 population for mapping QTL affecting growth traits, an F1 population was first constructed by mass crossing of 50 wild Asian seabass that captured from Southeast Asia (Thailand, Indonesia, Malaysia and Singapore). Two F1 individuals selected as F2 parents based on their performance difference in growth traits were then intercrossed to generate the F2 mapping family (F2S2), which contained 359 full sib individuals. The F2 individuals were raised communally in a 20-ton indoor tank and maintained on standard feeding regimes until 270 days post hatch. All individuals were tagged at the age of 3 months with microchips (Opulent, Singapore). Phenotypic data including body weight, standard length and total length were recorded on each fish of the family at the age of 6 and 9 months (Table 1), using methods described in Fishbase (http://www.fishbase.org/Identification/Morphometrics/centimeters/Index.php). Fin samples were collected from each fish and kept at −80°C for DNA isolation.
For LD-based association mapping, four F2 populations (F2S1, F2S5, F2S6, F2S7) were used. Each population was established from mating of unrelated 25 F1 sires with 25 F1 dams, whose parents originated from natural populations in Southeast Asian. The genetic relatedness and genetic distances among individuals were evaluated by using software PAPA  and MEGA 4 . Only unrelated individuals were used in mass crosses. Parentage assignment for each population was conducted to reconstruct pedigrees. Each population containing around 2000 fishes was transferred to a 7-ton indoor tank. The fishes were fed with commercial dry feed. After three months of growing in the tanks, two extreme trait groups (large size group: 60–200 fishes, and small size group: 60–100 fishes) from each tank were sampled and phenotypic data from each fish were measured. All experiments in this work have been reviewed and approved by the IACUC Committee on Bioethics of the Temasek Life Sciences Laboratory, Singapore.
Microsatellite and SNP genotyping
Based on the Asian seabass linkage maps [15, 16], 123 microsatellite markers shown polymorphism in parents of the mapping family (F2S2) were selected for genotyping analysis. These microsatellite markers distributed evenly on 24 linkage groups and were genotyped on a 3730xl DNA analyzer (Applied Biosystems) after PCR amplification using fluorescently labeled primers. The two parents and their offspring (N = 359) were genotyped using microsatellites as described previously .
Based on our Asian seabass transcriptome database (unpublished data), 19 ESTs showing high similarity to genes related to growth, nutritional condition or environmental adaptation were aligned with genomic sequence data from Stickleback (Gasterosteus aculeatus; http://www.ensembl.org/Gasterosteus_aculeatus/Info/Index). Primer sites in conserved exon regions were identified. One or two primer pairs for each gene (see Additional file 1: Table S1) were developed based on the EST sequences using the program PrimerSelect (DNASTAR, Wilmington, DE). Genotypes of SNPs in candidate genes were detected by sequencing of PCR products. SNP genotype was analyzed by using software Sequencher 4.9 (Gene Codes, Ann arbor, USA). To design primers for real-time PCR, the exon/intron boundaries and divergence regions from two IFABP homologs were first identified by using the software Sequencher 4.9. One primer pair that spans exon/intron boundary on the mRNA for each homolog was designed using the program PrimerSelect (DNASTAR, Wilmington, DE). This primer design allows differentiating the two homlogs and doesn’t amplify genomic DNA.
Fasting experiment and sampling
Twenty-eight fishes (at ~ 1 year old) that maintained at the fish facility of Temasek life sciences laboratory were used for functional analysis. Prior to fast experiment the fishes were housed in one tank containing 2000 L of freshwater. On day 0 (starvation day 0), seven tissues (kidney, heart, muscle, intestine, spleen, eye and liver) from four fishes were taken and kept in Trizol reagent. The remaining twenty-four fishes were then transferred to two tanks holding 1000 L of freshwater. The twelve fishes in control tank were fed with commercial dry feed twice per day. In test tank, a fast regime was imposed on twelve fishes for up to 12 additional days under the same conditions as the control except without feeding. Four fishes from each tank were sacrificed on days 3, 6 and 12 of the fast period. Intestine samples were kept in Trizol reagent (Invitrogen, CA, USA) for RNA extraction.
Cloning and gene expression analysis of the IFABP gene by quantitative real-time PCR (qPCR)
The full-length cDNA libraries described as in Xia et al. [47, 48] was used for cloning and characterization of the IFABP gene. For analysis of gene expression patterns 10 times dilution of the cDNA were assayed as template by qPCR using EF1A and RPL13A as control (see Additional file 1: Table S1), as suggested by Tang et al.  that the EF1A and Rpl13A genes are more suitable as a reference gene panel for zebrafish tissue analysis. PCR reactions were performed with the iQ SYBR Green Supermix (Bio-Rad, Hercules, CA, USA) as described by the manufacturer in an iQ™5 Real Time PCR Detection Systems (Bio-Rad, Hercules, CA, USA). PCRs were performed in four biological replicates and three technical replicates. For the analysis of the changes of gene expression, the values of qPCR reactions were normalized to EF1A and RPL13A gene expression, calculated by the ΔΔCT method  as implemented by the software Bio-Rad iQ5 (2.0) (Bio-Rad, Hercules, CA, USA). P-value significant difference was calculated using the student’s T-Test module installed in the Microsoft Office Excel 2008 program with parameters setting as two-tailed distribution and two-sample equal variance.
Linkage and QTL analysis
A sex averaged genetic linkage map was generated for the full sib family containing 359 individuals using JoinMap linkage analysis software packages . A LOD score of 3.0 was used as a threshold value for defining linkage groups. Marker distances were evaluated according to the Kosambi formula (in cM). QTL analyses were performed with MapQTL version 5.0  by using the multiple QTL model mapping analysis with selected cofactor markers. The LOD profiles were generated at 1 cM intervals along each linkage group to identify the most likely QTL position. The linkage group wide and genome wide significance thresholds for LOD scores were determined by implementing a bootstrapping method with 10000 permutation iterations at p = 0.01 and 0.05. QTL that exceeded the linkage group-wide LOD threshold at p < 0.05 was reported as suggestive QTL, while exceeding a genome-wide LOD threshold of p < 0.05 was considered evidence for a significant QTL effect, and exceeding a genome-wide LOD threshold of p < 0.01 was considered evidence for a very significant QTL effect. QTL affecting different traits were considered to be located in the same region if the chromosomal positions of these QTL were overlapped. Linkage map and QTL regions affecting growth were visualized using MapChart software (ver. 2.1)  and shown in Additional file 3: Figure S1.
LD-based association analysis and heterozygosity effect
In association analysis, each of the 4 populations from mass crosses was analyzed separately, as there are differences in the maintenance of each population in separated tanks, which could distort the estimation of association effects if they were analyzed together. The genetic structure of each population was subdivided by using a model-based approach in the STRUCTURE software . Kinship matrices were calculated using the SPAGeDi software package . Association analysis between the haplotypes and genotypes of the markers/genes and growth traits was then carried out using general linear model (GLM) function and mixed linear model (MLM) function that implemented in the software TASSEL 2.0.1 .
In order to investigate associations between phenotypic traits and molecular genotypes, two extreme trait groups (N = 30/group) for each trait were selected from 359 full sib fishes from the family F2S2. Excluding the loci with heterozygous genotypes (ab×cd) that segregated in two parents (with 4 alleles), the observed heterozygosities and homozygosities of remaining loci that used in the construction of the linkage groups were calculated for each group. The one-tailed probability for frequency distribution differences of heterozygosities between two extreme groups of a trait was calculated using on line Fisher's exact test calculator (Ver.3; http://www.danielsoper.com/statcalc3/calc.aspx?id=29).
Sequence analysis and phylogenetic analysis
The Asian seabass IFABP gene sequences (GenBank accession no.: JX678843-JX678844) were used to perform BLAST analyses for EST annotation. Twenty one IFABP protein sequences, sixteen of which from fish databases, were retrieved from the NCBI database (http://mirrors.vbi.vt.edu/mirrors/ftp.ncbi.nih.gov/blast/db/). The protein sequences were aligned using the Clustal X 1.83 program . Phylogenetic analyses was carried out using the MEGA 4 package  with neighbour-joining bootstrap tests (1000 times). Genetic distances among protein sequences were evaluated with Poisson correction by MEGA 4. Evolutionary relationships of IFABP protein genes were shown in Additional file 5: Figure S2.
Descriptive statistics and correlation analysis
Basic descriptive statistics, including number of observations (N), minimum values, maximum values, means and standard deviations (sd) and Pearson's correlation coefficient were calculated using the ‘data analysis’ module of the software Excel.
Quantitative trait loci
Single nucleotide polymorphism
Fatty-acid-binding protein gene
Intestine fatty-acid-binding protein alpha gene.
Tanksley SD, Young ND, Paterson AH, Bonierbale MW: RFLP mapping in plant breeding: new tools for an old science. Nat Biotech. 1989, 7: 257-264. 10.1038/nbt0389-257.
Doerge RW: Mapping and analysis of quantitative trait loci in experimental populations. Nat Rev Genet. 2002, 3: 43-
Reid DP, Szanto A, Glebe B, Danzmann RG, Ferguson MM: 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. 2005, 94 (2): 166-172. 10.1038/sj.hdy.6800590.
Moghadam H, Poissant J, Fotherby H, Haidle L, Ferguson M, Danzmann R: Quantitative trait loci for body weight, condition factor and age at sexual maturation in Arctic charr (Salvelinus alpinus): comparative analysis with rainbow trout (Oncorhynchus mykiss) and Atlantic salmon (Salmo salar). Mol Genet Genomics. 2007, 277 (6): 647-661. 10.1007/s00438-007-0215-3.
Houston RD, Bishop SC, Hamilton A, Guy DR, Tinch AE, Taggart JB, Derayat A, McAndrew BJ, Haley CS: Detection of QTL affecting harvest traits in a commercial Atlantic salmon population. Anim Genet. 2009, 40 (5): 753-755. 10.1111/j.1365-2052.2009.01883.x.
Baranski M, Moen T, Våge D: Mapping of quantitative trait loci for flesh colour and growth traits in Atlantic salmon (Salmo salar). Genet Sel Evol. 2010, 42: 17-10.1186/1297-9686-42-17.
Loukovitis D, Sarropoulou E, Batargias C, Apostolidis AP, Kotoulas G, Tsigenopoulos CS, Chatziplis D: Quantitative trait loci for body growth and sex determination in the hermaphrodite teleost fish Sparus aurata L. Anim Genet. 2012, 43 (6): 753-759. 10.1111/j.1365-2052.2012.02346.x.
Loukovitis D, Sarropoulou E, Tsigenopoulos CS, Batargias C, Magoulas A, Apostolidis AP, Chatziplis D, Kotoulas G: Quantitative trait loci involved in sex determination and body growth in the gilthead sea bream (Sparus aurata L.) through targeted genome scan. PLoS One. 2011, 6 (1): e16599-10.1371/journal.pone.0016599.
Antonello J, Massault C, Franch R, Haley C, Pellizzari C, Bovo G, Patarnello T, de Koning DJ, Bargelloni L: Estimates of heritability and genetic correlation for body length and resistance to fish pasteurellosis in the gilthead sea bream (Sparus aurata L.). Aquaculture. 2009, 298 (1–2): 29-35.
Boulton K, Massault C, Houston RD, de Koning DJ, Haley CS, Bovenhuis H, Batargias C, Canario AVM, Kotoulas G, Tsigenopoulos CS: QTL affecting morphometric traits and stress response in the gilthead seabream (Sparus aurata). Aquaculture. 2011, 319 (1–2): 58-66.
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 (4): 337-345.
Wringe BF, Devlin RH, Ferguson MM, Moghadam HK, Sakhrani D, Danzmann RG: Growth-related quantitative trait loci in domestic and wild rainbow trout (Oncorhynchus mykiss). BMC Genet. 2010, 11: 63-
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.
Cnaani A, Zilberman N, Tinman S, Hulata G, Ron M: Genome-scan analysis for quantitative trait loci in an F-2 tilapia hybrid. Mol Genet Genomics. 2004, 272 (2): 162-172.
Wang CM, Bai ZY, He XP, Lin G, Xia JH, Sun F, Lo LC, Feng F, Zhu ZY, Yue GH: A high-resolution linkage map for comparative genome analysis and QTL fine mapping in Asian seabass, Lates calcarifer. BMC Genomics. 2011, 12: 174-10.1186/1471-2164-12-174.
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.
Wang CM, Lo LC, Feng F, Zhu ZY, Yue GH: 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.
Chatziplis D, Batargias C, Tsigenopoulos CS, Magoulas A, Kollias S, Kotoulas G, Volckaert FAM, Haley CS: Mapping quantitative trait loci in European sea bass (Dicentrarchus labrax): The BASSMAP pilot study. Aquaculture. 2007, 272 (Suppl 1): 172-182.
Massault C, Hellemans J, Louro B, Batargias C, Houdt IS, Canario A, Volckaert FAM, Bovenhuis H, Haley CS, Koning DJ: QTL for body weight, morphometric traits and stress response in European sea bass Dicentrarchus labrax. Anim Genet. 2010, 41: 4-
Wang CM, Lo LC, Zhu ZY, Pang HY, Liu HM, Tan J, Lim HS, Chou R, Orban L, Yue GH: Mapping QTL for an adaptive trait: the length of caudal fin in Lates calcarifer. Mar Biotechnol. 2011, 13 (1): 74-82. 10.1007/s10126-010-9271-5.
Sánchez-Molano E, Cerna A, Toro MA, Bouza C, Hermida M, Pardo BG, Cabaleiro S, Fernández J, Martínez P: Detection of growth-related QTL in turbot (Scophthalmus maximus). BMC Genomics. 2011, 12: 473-10.1186/1471-2164-12-473.
Glazier AM, Nadeau JH, Aitman TJ: Finding genes that underlie complex traits. Science. 2002, 298: 2345-2349. 10.1126/science.1076641.
Darvasi A, Pisante-Shalom A: Complexities in the genetic dissection of quantitative trait loci. Trends Genet. 2002, 18: 489-491. 10.1016/S0168-9525(02)02767-1.
MacKay J, Boyle B, Kayal WE, Namroud M, Doerksen T, Cooke J, Isabel N, Beaulieu J, Rigault P, Bicho P: Gene mapping in white spruce (P. glauca): QTL and association studies integrating population and expression data. BMC Proc. 2011, 5 (Suppl 7): I6-10.1186/1753-6561-5-S7-I6.
Li GH, Cheung C, Xiao S, Lau K, Gao Y, Bow C, Huang Q, Sham P, Kung AW: Identification of QTL genes for BMD variation using both linkage and gene-based association approaches. Hum Genet. 2011, 130 (4): 539-546. 10.1007/s00439-011-0972-2.
Atchley WR: Ontogeny, timing of development, and genetic variancecovariance structure. Am Nat. 1984, 123: 519-540. 10.1086/284220.
Leamy LJ, Pomp D, Eisen EJ, Cheverud JM: Pleiotropy of quantitative trait loci for organ weights and limb bone lengths in mice. Physiol Genomics. 2002, 10: 21-29.
Larsen PF, Nielsen EE, Koed A, Thomsen DS, Olsvik PA, Loeschcke V: Interpopulation differences in expression of candidate genes for salinity tolerance in winter migrating anadromous brown trout (Salmo trutta L.). BMC Genet. 2008, 9: 12-
Nilsen TO, Ebbesson LOE, Madsen SS, McCormick SD, Andersson E, Bjornsson BT, Prunet P, Stefansson SO: Differential expression of gill Na+, K+−ATPase alpha- and beta-subunits, Na+, K+,2Cl(−) cotransporter and CFTR anion channel in juvenile anadromous and landlocked Atlantic salmon Salmo salar. J Exp Biol. 2007, 210 (16): 2885-2896. 10.1242/jeb.002873.
Norman JD, Danzmann RG, Glebe B, Ferguson MM: The genetic basis of salinity tolerance traits in Arctic charr (Salvelinus alpinus). BMC Genet. 2011, 12: 81-
Dekkers JCM, Hospital F: The use of molecular genetics in the improvement of agricultural populations. Nat Rev Genet. 2002, 3 (1): 22-32. 10.1038/nrg701.
Martyniuk CJ, Perry GML, Mogahadam HK, Ferguson MM, Danzmann RG: The genetic architecture of correlations among growth-related traits and male age at maturation in rainbow trout. J Fish Biol. 2003, 63 (3): 746-764. 10.1046/j.1095-8649.2003.00188.x.
Thornsberry JM, Goodman MM, Doebley J, Kresovich S, Nielsen D, Buckler ES: Dwarf8 polymorphisms associate with variation in flowering time. Nat Genet. 2001, 28: 286-289. 10.1038/90135.
Chen X, Shi Z: Sequence analysis of the full-length cDNA and protein structure homology modeling of FABP2 from Paralichthys Olivaceus. Bioinform Biol Insights. 2009, 3: 29-35.
Sharma MK, Denovan-Wright EM, Degrave A, Thisse C, Thisse B, Wright JM: Sequence, linkage mapping and early developmental expression of the intestinal-type fatty acid-binding protein gene (fabp2) from zebrafish (Danio rerio). Comp Biochem Physiol B Biochem Mol Biol. 2004, 138 (4): 391-398. 10.1016/j.cbpc.2004.05.009.
Her GM, Yeh YH, Wu JL: Functional conserved elements mediate intestinal-type fatty acid binding protein (I-FABP) expression in the gut epithelia of zebrafish larvae. Dev Dynam. 2004, 230 (4): 734-742. 10.1002/dvdy.20081.
Her GM, Chiang CC, Wu JL: Zebrafish intestinal fatty acid binding protein (I-FABP) gene promoter drives gut-specific expression in stable transgenic fish. Genesis. 2004, 38 (1): 26-31. 10.1002/gene.10248.
André M, Ando S, Ballagny C, Durliat M, Poupard G, Briançon C, Babin PJ: Intestinal fatty acid binding protein gene expression reveals the cephalocaudal patterning during zebrafish gut morphogenesis. Intl J Dev Biol. 2000, 44 (2): 249-252.
Pierce M, Wang Y, Denovan-Wright EM, Wright JM: Nucleotide sequence of a cDNA clone coding for an intestinal-type fatty acid binding protein and its tissue-specific expression in zebrafish (Danio rerio). Biochim Biophys Acta. 2000, 1490 (1–2): 175-183.
Woudstra TD, Drozdowski LA, Wild GE, Clandinin MT, Agellon LB, Thomson AB: The age-related decline in intestinal lipid uptake is associated with a reduced abundance of fatty acid-binding protein. Lipids. 2004, 39 (7): 603-610. 10.1007/s11745-004-1272-9.
Fawcett GL, Roseman CC, Jarvis JP, Wang B, Wolf JB, Cheverud JM: Genetic architecture of adiposity and organ weight using combined generation QTL analysis. Obesity. 2008, 16: 1861-1868. 10.1038/oby.2008.300.
Andersson L, Haley CS, Ellegren H, Knott SA, Johansson M, Andersson K, Andersson-Eklund L, Edfors-Lilja I, Fredholm M, Hansson I: Genetic mapping of quantitative trait loci for growth and fatness in pigs. Science. 1994, 263: 1771-1774. 10.1126/science.8134840.
Marklund L, Nystrom PE, Stern S, Andersson-Eklund L, Andersson L: Confirmed quantitative trait loci for fatness and growth on pig chromosome 4. Heredity. 1999, 82: 134-141. 10.1038/sj.hdy.6884630.
Loos RJF, Lindgren CM, Li S, Wheeler E, Zhao JH, Prokopenko I, Inouye M, Freathy RM, Attwood AP, Beckmann JS: Common variants near MC4R are associated with fat mass, weight and risk of obesity. Nat Genet. 2008, 40 (6): 768-775. 10.1038/ng.140.
Duchese P, Godbout M, Bernatchez L: PAPA (package for the analysis of parental allocation): a computer program for simulated and real parental allocation. Mol Ecol Notes. 2002, 2: 191-193. 10.1046/j.1471-8286.2002.00164.x.
Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.
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.
Xia JH, He XP, Bai ZY, Lin G, Yue GH: Analysis of the Asian seabass transcriptome based on expressed sequence tags. DNA Res. 2011, 18 (6): 513-522. 10.1093/dnares/dsr036.
Tang RY, Dodd A, Lai D, McNabb WC, Love DR: Validation of zebrafish (Danio rerio) reference genes for quantitative real-time RT-PCR normalization. Acta Bioch Bioph Sin. 2007, 39 (5): 384-390. 10.1111/j.1745-7270.2007.00283.x.
Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Δ Δ C(T)) method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.
Van Ooijen JW: Multipoint maximum likelihood mapping in a full-sib family of an outbreeding species. Genet Res. 2011, 93: 343-349. 10.1017/S0016672311000279.
Van Ooijen JW, Kyazma BV: MapQTL® 6. 2009, Wageningen, The Netherlands: Software for the mapping of quantitative trait loci in experimental populations of diploid species
Voorrips RE: MapChart: software for the graphical presentation of linkage maps and QTLs. J Hered. 2002, 93: 77-78. 10.1093/jhered/93.1.77.
Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005, 14: 2611-2620. 10.1111/j.1365-294X.2005.02553.x.
Hardy OJ, Vekemans X: SPAGEDi: a versatile computer program to analyse spatial genetic structure at the individual or population levels. Mol Ecol Notes. 2002, 2: 618-620. 10.1046/j.1471-8286.2002.00305.x.
Bradbury PJ, Zhang Z, Kroon DE, Casstevens TM, Ramdoss Y, Buckler ES: TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics. 2007, 23 (19): 2633-2635. 10.1093/bioinformatics/btm308.
Sievers F, Wilm A, Dineen DG, Gibson TJ, Karplus K, Li W, Lopez R, McWilliam H, Remmert M, Söding J: Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol. 2011, 7: 539-
We thank MAC, AVA for long-time co-operation on Asian seabass breeding, and our lab members for assisting in tissue sample collection. We are very grateful to our colleague Ms May Lee for English editing. This study is funded by the National Research Foundation Singapore under AQUAMARINOME 2016: Increasing Singapore’s Food Fish Production Through Aquaculture Genomics R & D (Award No: NRF2010NRF-CRP002-001). We thank MAC, AVA for long-time co-operation on breeding Asian seabass.
The authors declare that they have no competing interests.
GHY initiated and overviewed the Asian seabass project and finalized the manuscript. JHX and GHY designed the study. JHX, GL, XPH, PL, FS, FL and RJT conducted the experiments. JHX conducted the data analysis and wrote the manuscript. All authors have read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Table S1: Primer information for the SNP markers used in linkage mapping and gene expression analysis by quantitative RT-PCR. (XLS 52 KB)
Additional file 2: Table S2: Summary of the Asian seabass linkage map based on microsatellite and SNP markers. (XLS 23 KB)
Additional file 3: Figure S1: The averaged linkage map and the positions of QTL affecting growth traits in the F2S2 family of the Asian seabass. (DOC 47 KB)
Additional file 4: Table S3: Significant heterogeneity of genotypes between two extreme size groups in the full sib mapping family. (XLS 32 KB)
Additional file 5: Figure S2: Evolutionary relationships of IFABP protein genes. The evolutionary history is inferred using the neighbor-joining method. Bootstrap confidence values are listed to the left of the nodes. The GenBank accession no. and species name for each IFABP protein sequence is given to the right. The two IFABP isoforms that cloned from the Asian seabass are underlined. A scale near the bottom relates the length of a branch to the number of the amino acid substitutions that have taken place along the branch (evolutionary distances). (TIFF 1 MB)
About this article
Cite this article
Xia, J.H., Lin, G., He, X. et al. Whole genome scanning and association mapping identified a significant association between growth and a SNP in the IFABP-a gene of the Asian seabass. BMC Genomics 14, 295 (2013). https://doi.org/10.1186/1471-2164-14-295
- Single nucleotide polymorphism
- Growth trait
- Candidate gene
- Quantitative trait locus