Detection of Growth-Related Quantitative Trait Loci and High-Resolution Genetic Linkage Maps Using Simple Sequence Repeat Markers in the Kelp Grouper (Epinephelus bruneus)

To initiate breeding programs for kelp grouper (Epinephelus bruneus), the establishment of genetic linkage maps becomes essential accompanied by the search for quantitative trait loci that may be utilized in selection programs. We constructed a high-resolution genetic linkage map using 1055 simple sequence repeat (SSR) markers in an F1 family. Genome-wide and chromosome-wide significances of growth-related quantitative trait loci (QTLs) (body weight (BW) and total length (TL)) were detected using non-parametric mapping, Kruskal-Wallis (K-W) analysis, simple interval mapping (IM) and a permutation test (PT). Two stages and two families of fish were used to confirm the QTL regions. Ultimately, 714 SSR markers were matched that evenly covered the 24 linkage groups. In total, 509 and 512 markers were localized to the female and male maps, respectively. The genome lengths were approximately 1475.95 and 1370.39 cM and covered 84.68 and 83.21 % of the genome, with an average interval of 4.1 and 4.0 cM, in females and males, respectively. One major QTL affecting BW and TL was found on linkage group EBR 17F that identified for 1 % of the genome-wide significance and accounted for 14.6–18.9 and 14.7–18.5 % of the phenotypic variance, and several putative QTL with 5 % chromosome-wide significance were detected on eight linkage groups. Furthermore, the confirmed results of the regions harboring the major and putative QTLs showed consistent significant experiment-wide values of 1 and 5 % as well as a chromosome-wide value of 5 %. We identified growth-related QTLs that could be applied to find candidate genes for growth traits in further studies, and potentially useful in MAS breeding. Electronic supplementary material The online version of this article (doi:10.1007/s10126-015-9673-5) contains supplementary material, which is available to authorized users.


Introduction
The kelp grouper (Epinephelus bruneus) is a commercially important marine fish in East Asia. This species belongs to the subfamily Ephinephelinae, family Serranidae, and order Electronic supplementary material The online version of this article (doi:10.1007/s10126-015-9673-5) contains supplementary material, which is available to authorized users. Perciformes. Groupers, or Serranidae in general, are protogynous, which means they first start life as a female fish and then later switch into being males once they pass a certain size threshold or due to social cues (Lee et al. 2002;Tsuchihashi et al. 2003;Yeh et al. 2003) and matures at more than 6 years of age . The kelp grouper is a carnivorous fish that feeds on small fish and crustaceans. Generally, juvenile kelp groupers are found in shallow water estuaries (Heemstra and Randall 1995) and coastal areas, while the adult fish inhabit waters ranging 20-200 m of depth ) around coral reefs, rocky reefs, and mud bottom areas. The maximum size of a kelp grouper is reported to be about 128 cm in length and 33 kg in body weight (Tupper and Sheriff 2008). At least 16 species of grouper, including the kelp grouper, have been used successfully in aquaculture in many countries in East Asia (Tupper and Sheriff 2008).
In Japan, the kelp grouper has a high value because of its high market demand and low quantity of catch in natural waters (Mitcheson et al. 2003). Recently, this species was listed as a vulnerable species by the International Union for Conservation of Nature and Natural Resources (the IUCN Red List of Threatened Species) because of the rapid decrease in the natural population (Thierry et al. 2008). The kelp grouper is a target species for aquaculture in Japan (Fui et al. 2014); however, during artificial larval rearing, high mortality is frequent in the early life stages (Sawada et al. 1999). In addition, the kelp grouper grows slowly in farms and a prolonged farming period is required to reach a marketable size. To date, domestication of broodstock and a selective breeding program on a commercial scale for the kelp grouper in Japan have not yet been fully developed. Marker-assisted selection (MAS) based on quantitative trait loci (QTLs) is an effective method to improve quantitative traits (Max and Anatoly 2007) such as slow growth and high mortality in the larval stage of groupers.
SSR markers are highly polymorphic and show high inheritance and codominance of inheritance, making them suitable to identify homozygotes and heterozygotes. They are usually evenly distributed throughout the genome, and their results are simple to interpret,  *P≧0.05 normal distribution of phenotypic a This is the lower bound of the rue significance highly reproducible, and easily automated (Liu 2007). Thus, they are useful to construct a genetic linkage map. Nevertheless, the detection of a reasonable proportion of QTLs segregating in a population requires a large number of markers to increase the accuracy of QTL detection. Economic traits in aquaculture fish, especially growth-related quantitative traits, are the main goals  for improvement in a genetic breeding program. Growth-related traits have been measured and reported in several economically important marine fishes (Yue 2013). To study growth-related traits in fish and other species by molecular tools is complex, because growth-related traits are influenced not only by genetics, but also by the environment (Abraham et al. 2007;Molano et al. 2011).
In 2013, the first-generation genetic linkage map for t h e k e l p g r o u p e r w a s c o n s t r u c t e d u s i n g 2 2 2 microsatellite markers, covering 23 and 25 linkage groups in the male and female maps, with marker intervals of 5.0 and 6.7 cM, respectively ). In the present study, a high-resolution genetic linkage map and a genome scan for QTLs affecting growth-related traits (BW and TL) in F 1 progeny of kelp groupers were conducted. These results could be used to investigate candidate genes that will accelerate genetic improvement using MAS breeding programs in the kelp grouper.

Reference Family and DNA Extraction
Paternal half-sib F 1 progeny from two families (families A and B) produced from two females and a single male were used. The fish were taken from recently derived wild broodstock of the kelp grouper maintained at the Ehime Fisheries Research Center, Japan. Fish were measured at two timepoints. Stage I at 5 months post-hatching (average total length of 150 mm), and stage II at 11 months post-hatching. At stage I, individual fish were tracked using an embedded passive integrated transponder (PIT) tag. This facilitated comparisons of fast growth phases that occur in the fish at this point in their development. All fish were measured for body weight (BW) and total length (TL). In total, 360 and 163 progeny in stages I and II of family A; and 112 and 45 progeny in stages I and II of family B were measured for BW and TL. Fin clip samples were collected and kept in absolute ethanol (99.9 % ethanol solution). DNA extraction was carried out from these samples using the Agincourt DNAdvance Genomic DNA Isolation Kit (Beckman Coulter, USA),  Map distances are shown in centimorgans (cM) LG linkage group, G e1 genome estimate size 1, G e2 genome estimate size 2 following the manufacturer's recommended protocol. The quality and quantity of the extracted DNA was quantified using a spectrophotometer (Untrospec 2100 pro, GE USA) and the DNA was diluted to 10 ng/μL for PCR. The high-resolution genetic linkage maps to find candidate growth-related QTL regions were constructed by using the parents and 90 F 1 progeny in stage II of family A. After that, all progeny in both stages of families A and B were used to confirm the candidate QTL regions.

SSR Markers and Genotyping
A total of 2348 microsatellite-enriched segments from the kelp grouper were developed using next-generation sequencing (NGS) by the GS FLX system (Roche, Switzerland) (denoted as the EBR series) (Kubota et al. 2014) and 889 simple tandem repeats (STR) markers were obtained from the NCBI database of a cross section of species in the subfamily Epinephelinae (denoted as the STR series) (Chapman et al. 1999;Dong et al. 2008;Liu et al. 2008;Lo and Yue 2007;Mokhtar et al. 2011;Ramirez et al. 2006;Renshaw et al. 2010;Rivera et al. 2003;Zeng et al. 2008;Zhao et al. 2009a, b;Zhu et al. 2005). In total, 1867 SSR markers (1466 EBR markers and 401 STR markers) were designed using the TROLL program at http://wsmartins.net/websat/ (Martins et al. 2009) under the default settings and considering a product size of 100-250 bp. For the SSR markers, the forward primers were labeled with tetrachloro-6-carboxy-fluorescein (TET) fluorescent dye at the 5′-end. Polymerase chain reactions (PCR) were performed in 11 μl volumes containing 50 ng of genomic DNA, 1× Ex Taq buffer (Mg 2+ free), 2.0 mM MgCl 2 , 0. 2 mM dNTP, 1 % BSA, 0.025 U of Taq polymerase (Takara: Ex-Taq™ (Mg 2+ free buffer)), 0.5 pmol/μL of the reverse primer, and 0.05 pmol/μL of the forward primer. Cycle amplification was performed on an MJ PTC-100 (Bio-Rad, USA), with the program conditions of 95°C for 5 min for initial denaturation; followed by 36 cycles of 30 s at 95°C, 1 min at the annealing temperature 56°C, and 1 min at 72°C, and a final extension at 72°C for 10 min. The amplified products were mixed with an equal volume of loading buffer (98 % formaldehyde, 10 mM EDTA, and 0.05 % bromophenol blue), heated for 10 min at 95°C and then immediately cooled on ice. The samples were separated on a 6 % polyacrylamide gel containing 7 M urea and 0.5× Trizma base/Boric Acid/EDTA-2Na (TBE) buffer and 40 % Page-plus (Amresco, USA) with a 500-bp DNA ladder (GeneScan™-500 TAMRA™). Electrophoresis was performed using 0.5× TBE buffer at a constant voltage of 1800 V for 1.5 h. After electrophoresis, the gel was scanned and imaged using an FMBIO III Multi-View fluorescence image analyzer (Hitachi-soft, Japan).

Linkage Analysis
Linkage analysis was performed using LINKMFEX version 2.3 (Danzmann 2006). This application can separate alleles that originated from males or females. To avoid errors during genotyping, the accuracy of genotypes in their progeny was checked from parental male and female alleles. Genotype data were converted to a backcross format even though the grandparent genotype was unknown. Pairwise analysis was performed, and markers were sorted into linkage groups at a logarithm of odds (LOD) threshold of 4.0. Linkage phases were determined retrospectively by examining the assortment of alleles among linked markers. The goodness of fit of the chi-square analysis ( χ 2 ) was used to test for Mendelian segregation distortion of the locus. Therefore, the distance of the marker was estimated on each linkage group, assuming the Kosambi mapping function. Double recombination was checked using the application in Map Manager QTX (Manly et al. 2001). Graphical representation of the linkage groups was performed using MAPCHART version 2.1 (Voorrips 2002). In addition, a consensus linkage map was constructed using JoinMap version 4 (Ooijen 2006) and the module of the combined group The recombination rate female/male (1.12:1) cM centimorgan, G oa observed genome length, G e1 genome estimate size 1, G e2 genome estimate size 2, Ge genome length estimate, C f genome coverage of (%)

Estimation of Genome Size and Coverage
A sex-specific map of genome length was estimated by two different calculation methods. First, genome estimation size 1 (G e1 ) was calculated by adding 2s, where s is the average framework marker spacing that was calculated by dividing the summed length of all the genetic linkage groups by the number of intervals (number of markers minus the number of genetic linkage groups) to the length of each genetic linkage group, accounting for chromosome ends beyond the terminal markers coverage. Second, genome estimation size 2 (G e2 ) was calculated by multiplying the length of each genetic linkage group by a factor (m + 1)/(m − 1). Where m is the number of framework markers for each genetic linkage group (Chakravarti et al. 1990). The estimated genome length (G e ) for each sex was used as an average of the two estimates (Fishman et al. 2001;Sanchez et al. 2010). The genome coverage for each sex was calculated as the  observed genome length (G oa ) divided by the estimated genome length (G e) (Song et al. 2013) while the observed genome length (G oa ) was taken as the combination of total length in all linkage group.

QTL Analysis
First, the normality of the phenotypes (BW and TL) was tested using the Kolmogorov-Smirnov test (N>50) and Shapiro-Wilk test (N<50), implemented in SPSS 16.0 package. The data were converted to Z scores before analysis using MapQTL software. QTL analysis was carried out using MapQTL 5 software (Ooijen 2004). Ninety F 1 progeny from stage II of family A were used to find candidate QTLs. A non-parametric Kruskal-Wallis analysis was used to determine the significance level of all marker loci associated with the growth-related traits (BW and TL). Meanwhile, simple interval mapping was used to detect significant associations with growth-related traits and marker loci in the data sets under the significant threshold of genome-wide (P value < 0.01 and P value < 0.05) and chromosome-wide (P value<0.05) analyses. A minimum LOD threshold of 4.0 was used for determining a significant QTL and the percentage of phenotypic variance of each QTL. Permutation tests were performed (1000 replicates) to determine the LOD threshold by type one error. The significant thresholds derived from the permutation tests was estimated by dividing the nominal P value by the total number of chromosomes (Churchill and Doerge 1994;Ozaki et al. 2013). A graphical representation of the significant QTLs was constructed using MAPCHART version 2.1 and MapQTL 5. The results of the growth-related QTL regions of stage II family A were confirmed to be reproducible in the other stage and family.

Correlation of Phenotypes and Growth-Related Traits in Families A and B
The correlation of phenotypes was tested using Pearson's correlation coefficient. The results showed a high correlation between BW and TL in both stages of the two families (Table 1). The normal distribution of the phenotype was tested by a Kolmogorov-Smirnov test or Shapiro-Wilk test depending on the number of samples (Table 2, Additional file 1). The high correlation between BW and TL and normal distribution of phenotypes in stage II of family A led us to select family A to construct the high genetic linkage map and to screen candidate QTL regions.

High-Resolution Genetic Linkage Map and Genome Coverage
A total of 1867 SSR markers were designed. Of them, approximately 1050 SSR markers were polymorphic (56.2 %), and composed 905 EBR and 145 STR SSR markers. Ultimately, 714 SSR markers were used to construct a linkage map with reference species. The list of SSR markers used for mapping is given in additional file 2. Twenty-four genetic linkage groups (LG1-LG24) were identified. The female linkage map contained 509 markers distributed in 24 linkage groups (EBR 1F-EBR 24F) (Fig. 1). The total genome size of the female map was estimated as 1249.8 cM. The number of markers per linkage group varied from 5 to 29, with an average of 21; the longest linkage group of the female map extended to 65.4 cM (EBR 7F). Meanwhile, 512 markers were distributed in 24 linkage groups of the male map (EBR 1M-EBR 24M) (Fig. 1). The total genome was estimated at 1140.3 cM. The longest linkage group of the male extended to 58.0 cM (EBR 1), while the average number of markers per linkage group was 21, and varying from 9 to 31. The framework interval in each group was estimated based on the distance between clusters or markers, because some markers located on the same cluster. The female and male linkage maps comprised 305 and 285 framework, respectively, and the average interval between markers was 4.1 and 4.0 cM, respectively (Tables 3 and 4).
Recombination rate between the sex-specific genetic linkages were estimated by co-segregation markers. At least two SSR markers shared loci in the female and male maps and could be used to calculate the recombination rate among adjacently paired markers. The total length of genetic distance obtained from 24 genetic linkage groups (LGs) were 1249.8 and 1140.3 cM in female and male maps, respectively. The relative recombination ratio between females and males in these pairs was 1.12:1, which indicated that female LGs had a higher recombination rate than male except for LGs 1,8,11,12,14,and 19 (Table 5).
Genome length (G e ) was estimated as approximately 1475.95 and 1370.39 cM in the female and male maps, respectively. The female map was 1.07 times longer than the male map. Only nine LGs (1,6,11,12,14,19,21,23, and 24) on the male map were longer than the female map. The genome coverages of the female and male maps were estimated at 84.68 and 83.21 %, respectively (Table 4).

Screening Candidate QTL Regions
Screening for candidate QTL of BW using the Kruskal-Wallis analysis of stage II family A (90 progeny) identified 5, 23, and 6 of the 34 total markers were significant (P<0.01) on three linkage groups corresponding to chromosomes EBR 13F, EBR 17F, and EBR 18M (data not shown). The results of simple interval mapping and the permutation test showed a significant major QTL (qBW17f) at the 1 and 5 % genomewide level on the linkage group EBR 17F (Fig. 2). The LOD score of qBW17f (LOD=4.09) was higher than the genomewide LOD significance threshold of 3.7. This candidate major QTL region was detected close to the SSR markers Ebr00153FRA, Ebr00702FRA, Ebr00314FRA, and EguSTR119DB, and accounted for 14.6-18.9 % of the phenotypic variance with 1.00-1.13 of the additive effect (Table 6). While seven putative QTLs (qBW5f, qBW10m, qBW13f, qBW15m, qBW18m, qBW19f, and qBW21f) were significant at a 5 % chromosome-wide level on linkage groups EBR 5F, EBR 13F, EBR 19F, and EBR 21F of the female map, and linkage groups EBR 10M, EBR 15M, and EBR 18M of the male map (Fig. 3) and could explain 7.5-12 % of the phenotypic variance with 0.70-0.92 of the additive effect (Table 6). Due to a high significance level (P<0.01) in the nonparametric K-W analysis and the LOD score of a candidate major and two putative QTLs exceeded the high significance thresholds (P < 0.01 and P < 0.05) of genome-and chromosome-wide after simple interval mapping and permutation analysis. Then, 35 representative microsatellite markers (Additional file 3) related with these QTLs from three linkage groups (EBR 13F , EBR 17F, and EBR 18M) were used to confirm the candidate QTL regions in all progeny in the two stages of families A and B.

Confirmation of the Candidate QTL Regions
Thirty-five marker loci from three candidate QTL regions of three linkage groups affecting BW in stage II family A were used to confirm the QTL region in the other stage of the same family and in the other family by collecting genotype data in both stages of the two families. In the case of the stage II family analysis, the number of progeny analyzed for the trait analysis increased from 90 to 163 progeny. For family A, the K-W test results showed that eight markers from linkage groups EBR 13F and EBR 17F of the female map showed consistently significant results (P<0.001) in stage II. Of them, three markers (Ebr00254FRA, Ebr00314FRA, and EguSTR119DB) showed the highest consistently significant results (P<0.0005), while only two markers (ElaSTR366DB and Ebr00443FRA) showed consistently significant results (P<0.005) in the male map (Table 7). Simple interval mapping on a chromosome-wide basis was then performed in each stage. The results showed only three QTLs (qBW13f, qBW17f, and qBW18m) in stage II were still significant. However, the results of interval mapping in stage II showed decreasing LOD scores (4.09 to 3.17) from the genome-wide analysis, with an LOD experimental-wide significance threshold of 2.0 ( Fig. 4a) with the LOD maximum locus (qBW17f) could explain phenotypic variance ranging 5.9-8.6 % with 0.49-0.59 of the additive effect of the BW traits. In contrast, for two candidate QTLs (qBW13f, qBW18m) on linkage groups EBR 13F and EBR 18M, their LOD scores increased from 2.5 to 3.38, and from 2.47 to 2.9, respectively, under the experiment-wide analysis. LOD significant threshold of 2.0 and 2.0 (Fig. 4b, c) with the region of LOD maximum locus (qBW13f and qBW18m) could explain phenotypic variance ranging 4.2-9.1 and 5-7.9 % with 0.42-0.62 and 0.44-0.56 of the additive effect of BW traits (Table 8). Nevertheless, we could not find any consistently significant results in stage I of family A.
In family B, the results showed only one marker (Ebr00702FRA) on linkage group EBR 17M in stage I, which presented consistent highly significant results (P< 0.001), was a putative QTL (qBW17m-1). It had a LOD score of 2.65, which was higher than the chromosomewide LOD significance threshold of 2.0, with a range of 10.3 % of the phenotypic variance with 0.64 of the additive effect (Table 8, Fig. 5). Meanwhile, other significant regions in all linkages were rejected as QTLs in stage II of family B.

Association of Growth-Related Trait QTL Regions and TL
In this study, we also measured another phenotype, TL, which was highly correlated with BW of fish (Pearson correlation coefficient test P<0.01), particularly in stage II of both families. For stage II of family A, the results of the K-W analysis and simple interval mapping showed significant loci in eight linkage groups (EBR 5F, EBR 7F, EBR 8F, EBR 13F, EBR 10M, EBR 17F, EBR 18M, and EBR 22M). The LOD score of a major QTL (qTL17f) effected to TL in linkage group EBR 17F was 4.0. with genome-wide significance (P<0.01). This QTL region could explain 14.7-18.5 % of the phenotypic variance and 0.99-1.12 of the additive effect of TL trait. Meanwhile, another region with an LOD maximum locus in the other linkage group had a value that exceeded the chromosome-wide value and could explain 7.0-11.3 % of the phenotypic variance and 0.69-0.89 of the additive effect of the TL trait (Table 9). Moreover, we confirmed all the Fig. 3 Localization of putative QTLs for body weight traits in the female map of family A. EBR (linkage group) F marker distance on the female map; EBR (linkage group) M marker distance on the male map. a qBW5f: QTL for body weight on EBR5F. b qBW13f: QTL for body weight on EBR13F. c qBW19f: QTL for body weight on EBR19F. d qBW21f: QTL for body weight on EBR 21F. e qBW10m: QTL for body weight on EBR 10M. f qBW15m: QTL for body weight on EBR 15M. g qBW18m: QTL for body weight on EBR 18M. Map positions and LOD scores were based on simple interval mapping, QTL analysis was performed using the software MapQTL 5. LOD limit of detection (significance threshold), P c chromosome-wide significance threshold      (Table 10). LOD analysis showed a decreasing LOD score from 4.00 to 3.25 at an LOD experimental-wide significance threshold of 2.0, in the candidate major QTL (qBW17f) on the linkage group EBR 17F. By contrast, the confirmation of two candidate putative QTL regions (qTL13f and qTL18m) demonstrated LOD scores that increased from 2.34 to 3.24 and 2.32 to 2.46, respectively, on the experiment-wide scale. LOD significant threshold of 2.0 and 2.0. The region of the LOD maximum locus (qTL13f and qTL18m) could explain phenotypic variance ranging from 3.9-8 .to 4.8-6.7 % of the phenotypic variance and 0.40-0.61 and 0.43-0.52 of the additive effect of the TL trait. As with the results for BW, we could not find any consistently significant values for stage I of family A or for both stages of family B (Table 8).

Discussion
The high-resolution genetic linkage maps of the kelp grouper produced in this study greatly enhanced the previous genetic linkage map for the kelp grouper which was developed by using 222 microsatellite markers. The previous female and male map consisted of 25 and 23 linkage groups with 67.2 and 67.8 % of genome coverage and 1.5:1 of average recombination ratio ). In the new genetic linkage map, 714 SSR markers were mapped in the 24 linkage groups, which is consistent with the diploid chromosome number of the kelp grouper (2N=48) (Lan 2009). About 509 and 512 markers were identified and evenly covered the 24 linkage groups of the female and male maps, respectively. Only 10 of 714 markers remained as single markers. All of the microsatellite markers used in the previous genetic linkage map were also included and were consistently assigned in the same order and linkage groups in the present study, except for six markers. Of these, three markers (EguStr125DB, MiniSTR267DB, and Ebr00025FRA) and three other markers (MiniSTR266DB, Ebr00270FRA, and Ebr00253FRA) in linkage group EBR 24 and EBR 25 of the female map were moved to linkage groups EBR 23 and EBR 5 in the new female map, respectively. In addition, the genome coverage and average ratio of recombination between female and male maps were about 84.68, 83.21, and 1.12:1, respectively. This result revealed a large number of markers in the F 1 progeny that filled several gaps of the new linkage map, which led to a reduction in the average mapping interval and an increase of the genome coverage. Considering the average interval and the genome coverage, we conclude that the high-resolution genetic linkage map of the kelp grouper of this study offers a sufficient marker density to permit a preliminary genome-wide scan for QTLs for growth-related traits (Massault et al. 2008). In addition, markers from other grouper species could speed up the construction and completion of a genetic linkage map of the kelp grouper in the near future. The recombination rate of a gene located on a chromosome (autosomal) is different between females and males because of the number of crossing-over events that occur during meiosis I. Differences in recombination rates between sexes have been identified in many species; for example, humans (Dib et al. 1996), dogs (Wong et al. 2010), crocodiles (Miles et al. 2009), and fish. In fish, recombination rates have generally been reported to be higher in females compared to males ranging from 3.25:1 in rainbow trout (Sakamoto et al. 2000), 7.4:1 in the Japanese flounder (Coimbra et al. 2003), 1.37:1 in Atlantic salmon (Lien et al. 2011), 2.2:1 in the silver carp (Guo et al. 2013), 2:1 in the Atlantic halibut (Reid et al. 2007), 1.5:1 in the kelp grouper ), 1.03:1 in the orange-spotted grouper (You et al. 2013), and 1.19:1 in the white grouper (Dor et al. 2014). In this study, the recombination rate ratio between females and males was 1.12:1, which was lower than previous reports. This may reflect the increased number of markers linked to the male map rather than the female map, which would affect not only the density of the markers but also the recombination rate in all linkage groups. In the present study, we found that markers in the female and male maps were irregularly distributed and showed high clustering of markers in all linkage groups. These markers tended to be compressed in the telomeric and centromeric regions of the female and male maps. A higher rate of recombination in the female and male maps probably occurred near the EBR (linkage group) M marker distance on the male map. qBW17m-1: QTL number 1 for body weight on EBR 17M; Map positions and LOD score based on simple interval mapping. QTL analysis was performed using the software MapQTL 5. LOD limit of detection (significance threshold), P e experiment-wide significance threshold, P c ; chromosome-wide significance threshold centromeric and the telomeric regions (You et al. 2013). This could be explained by the higher frequency of recombination in females near the centromeric regions during oogenesis. Similarly, more frequent recombination in males was also found near the telomeres during meiosis (Strachan and Read 2011;You et al. 2013). For indicating the centromeric or telomeric region in female and male maps, these two regions were observed by the map distance between markers. In the case of high recombination, the maps will present high distance between markers or clusters. The distances between markers in the centrometric region were assessed to be larger than other sites (telemetric). Similar to the male map, the markers or clusters in telemetric regions were estimated to have a larger distance than the centrometric region. The difference in sex recombination is an important factor in the implementation of marker-assisted selection using QTL-associated mapping.
The growth-related quantitative trait QTLs in this study were identified using F 1 progeny of the kelp grouper. This was different from other studies that performed QTL mapping using F 2 generation from F 1 crosses in a genetically different line or F 2 back-cross (Hayashi and Awata 2004), such as the Pacific white leg shrimp (Andriantahina et al. 2013). Kelp groupers are protogynous hermaphrodites and it would take a long time to produce an F 2 generation. This type of reproductive system takes a longer time for the sex reversal from male to female when they exceed a certain age or body size. In the kelp grouper, it takes more than 6 years of culture for the fish to reach maturity (before the first maturation and spawning). This is too long to create an F 2 generation. This explains our choice of producing F 1 progeny for the QTL    study. In the past decade, the analysis of QTLs using F 1 prog-