Next Article in Journal
Nucleotide Excision Repair and Vitamin D—Relevance for Skin Cancer Therapy
Next Article in Special Issue
Prediction of Toxin Genes from Chinese Yellow Catfish Based on Transcriptomic and Proteomic Sequencing
Previous Article in Journal
Present Advances and Future Perspectives of Molecular Targeted Therapy for Osteosarcoma
Previous Article in Special Issue
Characterization of Virulence Properties of Aeromonas veronii Isolated from Diseased Gibel Carp (Carassius gibelio)
Article

Genome-Wide Mapping of Growth-Related Quantitative Trait Loci in Orange-Spotted Grouper (Epinephelus coioides) Using Double Digest Restriction-Site Associated DNA Sequencing (ddRADseq)

by 1,2,†, 2,3,†, 1,2, 2, 3, 3, 4, 3,5,*, 3,4,5,* and 1,2,6,*
1
BGI Education Center, University of Chinese Academy of Sciences, Shenzhen 518083, China
2
Shenzhen Key Lab of Marine Genomics, Guangdong Provincial Key Lab of Molecular Breeding in Marine Economic Animals, State Key Laboratory of Agricultural Genomics, BGI, Shenzhen 518083, China
3
State Key Laboratory of Biocontrol, Institute of Aquatic Economic Animals and Guangdong Provincial Key Laboratory for Aquatic Economic Animals, School of Life Sciences, Sun Yat-Sen University, Guangzhou 510275, China
4
Marine Fisheries Development Center of Guangdong Province, Huizhou 510610, China
5
South China Sea Bio-Resource Exploitation and Utilization Collaborative Innovation Center, Guangzhou 510275, China
6
College of Life Sciences, Shenzhen University, Shenzhen 518060, China
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Academic Editors: Jun Li and Li Lin
Int. J. Mol. Sci. 2016, 17(4), 501; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms17040501
Received: 1 March 2016 / Revised: 28 March 2016 / Accepted: 29 March 2016 / Published: 6 April 2016
(This article belongs to the Special Issue Fish Molecular Biology)

Abstract

Mapping of quantitative trait loci (QTL) is essential for the discovery of genetic structures that related to complex quantitative traits. In this study, we identified 264,072 raw SNPs (single-nucleotide polymorphisms) by double digest restriction site associated DNA sequencing (ddRADseq), and utilized 3029 of these SNPs to construct a genetic linkage map in orange-spotted grouper (Epinephelus coioides) using a regression mapping algorithm. The genetic map contained 24 linkage groups (LGs) spanning a total genetic distance of 1231.98 cM. Twenty-seven significant growth-related QTLs were identified. Furthermore, we identified 17 genes (fez2, alg3, ece2, arvcf, sla27a4, sgk223, camk2, prrc2b, mchr1, sardh, pappa, syk, tert, wdrcp91, ftz-f1, mate1 and notch1) including three (tert, ftz-f1 and notch1) that have been reported to be involved in fish growth. To summarize, we mapped growth-related QTLs in the orange-spotted grouper. These QTLs will be useful in marker-assisted selection (MAS) efforts to improve growth-related traits in this economically important fish.
Keywords: quantitative trait loci (QTL); genetic linkage map; growth-related genes; ddRADseq; orange-spotted grouper (Epinephelus coioides) quantitative trait loci (QTL); genetic linkage map; growth-related genes; ddRADseq; orange-spotted grouper (Epinephelus coioides)

1. Introduction

Orange-spotted grouper, Epinephelus coioides (Epinephelinae, Serranidae), a protogynous hermaphrodite species, is mainly distributed in the Indo-West-Pacific region [1]. It is an economically important aquaculture species in many Asian countries, especially in China, due to its desirable taste and as a source of nutrition [2]. Orange-spotted grouper has become an important edible fish species in live fish markets and is an important cultured fish for sale in markets in southeast China. Market demand for the orange-spotted grouper has prompted the development of fish families and populations characterized by lower food consumption and higher growth rates. The application of marker-assisted selection (MAS) to orange-spotted grouper will be useful to improve important traits, such as feed conversion rate, meat quality and disease resistance. These traits have a major effect on productivity and profitability.
High-density genetic linkage maps and quantitative trait loci (QTL) mapping provide a framework for the MAS program. At present, genetic linkage maps have been constructed in over 28 fish species and economic traits have been mapped in at least 13 fish species [3]. These include the Nile tilapia (orechromis niloticus) [4], Atlantic salmon (Salmo salar) [5], rainbow trout (Oncorhynchus mykiss) [6], and channel catfish (Ictalurus punctatus) [7]. Microsatellite-based linkage maps have been reported for the white grouper (Epinephelus aeneus) [8] and kelp grouper (Epinephelus bruneus) [9]. The density of the genetic linkage maps is typically determined by the marker type. Single-nucleotide polymorphisms (SNPs), the most stable and abundant form of genetic marker, are the ideal marker type for construction of high-density genetic linkage maps [10,11,12]. The advent of next-generation sequencing (NGS) has made it possible to discover thousands of SNPs dispersed throughout the genome in a single procedure, even if little genetic information is available on the species [13]. A SNP-based high-density genetic linkage map has been constructed for orange-spotted grouper by our lab [14], but QTL mapping was not reported.
This is the first report of a SNP-based high-resolution genetic linkage map of the orange-spotted grouper developed using double digest restriction site associated DNA sequencing (ddRADseq). We also report the identification of 27 significant growth-related QTLs and 17 corresponding genes in this economically important fish.

2. Results

2.1. Sequencing and Genotyping

In total, 8.2 × 108 raw 90 bp reads were generated and 7.3 × 108 clean reads were retained after removal of low-quality raw reads. These high-quality reads were partitioned into ddRAD tags. The average number of ddRAD tags in each individual was 10,428,571. For the analysis of the F1 mapping population, a total of 264,072 candidate ddRAD loci were inferred from 68 individuals. If both parents were homozygous, the SNPs were eliminated. The SNPs that were significant segregation distortion based on Mendelism by Chi-square analysis were also eliminated. Ultimately, 26,661 SNPs were obtained (see detailed information in Table S1).

2.2. Genetic Linkage Map

After construction of the genetic linkage map using the 26,661 SNPs, we obtained 24 linkage groups (LGs), which is consistent with the haploid chromosome number of the orange-spotted grouper [15]. A total of 3029 SNPs were successfully mapped to the 24 LGs (Table S2 and Figure S1) after discarding contradictory SNPs. These LGs spanned a total genetic distance of 1231.98 cM with the length of each LG ranging from 25.43 cM (LG18) to 111.47 cM (LG8). The average genetic length of the LGs was 51.20 cM. The number of SNPs in the different LGs ranged from 94 (LG12) to 150 (LG22) with an average number of 126. Detailed information about loci and SNPs is summarized in Table S3.

2.3. Growth Trait–Associated Quantitative Trait Loci (QTL) and Related Genes

In total, 27 significant QTLs for growth traits were found to be distributed on LG1, LG5, LG7, LG21 and LG24 of the orange-spotted grouper (Table 1 and Figure 1). These included 18 QTLs for body weight (BW) and nine QTLs for body length (BL). The majority of these QTLs clustered together on their respective LGs. Especially, QTLs (qLG5-1, qLG5-2, qLG5-3, qLG5-4, qLG5-5, qLG5-6, qLG5-7, qLG5-8, qLG5-9, qLG5-10, qLG5-11, qLG5-12, qLG5-13, qLG5-14, qLG5-15, qLG5-16, qLG5-17, qLG5-18, qLG5-19, qLG5-20 and qLG5-21) were found to be clustered in a narrow region (15.1–33.3 cM) on the LG5.
Among these, the logarithm of odds (LOD) value of qLG5-1 for body weight, which accounts for 24.3% of the phenotypic variation, located at 15.1–15.6 cM is the highest (4.7), while qLG5-16 located at 15.1–15.4 cM for the trait of body length had the LOD value of 4.4 and it explains 23.3% of the phenotypic variation. The other QTLs on LG5 for BW were detected at positions 17.0–17.5, 17.8–18.0, 18.7–19.4, 19.9–20.2, 20.4–20.8, 27.4–27.9, 28.1–28.3, 29.4–29.7, 32.2–32.4 and 33.0–33.3 cM, with LOD values of 2.6–4.2, accounting for 14.6%–29.1% of the phenotypic variation. The other BL QTLs on LG5 were detected at positions 17.0–17.1, 19.2–19.3, 19.9–20.2, 20.4–20.8 cM, with LOD values of 2.6–3.8. These accounted for 14.5%–22.2% of the phenotypic variation. The QTLs for BW and BL were distributed across all 24 LGs. However, there was only one QTL for BW on LG1 and one QTL for BL on LG2.
Seventeen genes, including fasciculation and elongation protein ζ-2-like (fez2), Dol-P-Man: Man(5)GlcNAc(2)-PP-Dol α-1,3-mannosyltransferase (alg3), endothelin converting enzyme 2 (ece2), armadillo repeat gene deleted in velocardiofacial syndrome (arvcf), solute carrier family 27 (fatty acid transporter), member 4 (sla27a4), tyrosine-protein kinase SgK223 (sgk223), calcium/calmodulin-dependent protein kinase 2 (camk2), proline-rich coiled-coil 2B (prrc2b), melanin-concentrating hormone receptor 1 (mchr1), sarcosine dehydrogenase (sardh), pregnancy-associated plasma protein-A (papp-a), spleen tyrosine kinase (syk), telomerase reverse transcriptase (tert), WD repeat-containing protein 91-like (wdrcp91), ftz transcription factor 1 (ftz-f1), multidrug and toxin extrusion protein 1-like (mate1) and neurogenic locus notch homolog protein 1-like (notch1), were identified (Table 2). Among these, 12 genes (arvcf, sla27a4, sgk223, camk2, prrc2b, mchr1, sardh, papp-a, syk, tert, wdrcp9 and ftz-f1), corresponding to nine QTLs, were clustered on LG5. Two genes (alg3 and ece2), representing a single QTL, were distributed on LG2. Three genes including fez2, mate4 and notch1 mapped as a single QTL to LG1, LG7 and LG21, respectively.

3. Discussion

In this study, we applied the technology of ddRADseq, an extension of RADseq [16]. This technique has the advantage of providing improved efficiency and robustness by utilizing two restriction enzymes (a frequently-cutting enzyme and a rare-cutting enzyme) [16,17,18,19]. Under-sampling in read count was a constant question resulting from biased read representation in pooled sequencing experiments among individual samples. ddRADseq increases its sturdiness compared to RADseq [20,21,22]. ddRAD-based genetic maps of Nile tilapia (Oreochromis niloticus) [23], Japanese eel (Anguilla japonica) [17], Midas cichlids (Amphilophus spp.) [24,25], and Eurasian perch (Perca fluviatilis L.) [19] have been successfully constructed. Using ddRADseq, 264,072 candidate RAD loci were inferred, and 3029 high-quality SNPs were retained after a series of filtering. The number of mapped genetic markers in the present study was less than that in our previous report [14]. However, ddRAD-seq was better than MSG (multiplexed shotgun genotyping) in the detection of SNPs because MSG with a stringent methodology may overlook some significant loci [26]. There were fewer markers mainly due to the use of a small F1 full-sib population.
QTL mapping is considered an efficient strategy for analyzing complex quantitative traits in a variety of fish species. For example, 17 QTLs were detected for the traits of body size (body length and weight) of Chinook salmon (Oncorhynchus tshawytscha) [27] and 11 QTLs were found to be related to the body shape of lake cichlid fishes [25]. Additionally, growth-related QTLs had been identified in Atlantic salmon (Salmo salar) [5] and Arctic charr (Saivelinus alpinus) [28,29]. Recently, a genetic linkage map of kelp grouper had been constructed using simple sequence repeat (SSR) markers, and growth-related (body weight and total length) QTL analysis was performed [30]. The high-resolution genetic map consisted of 714 SSR markers. One major growth-related QTL and several putative QTLs were detected.
The previous data were from the family without phenotype information and were used to construct a SNP-based high-density genetic linkage map. However, our current mapping data from ddRADseq were applied to construct a genetic map and perform QTL analysis based on the available phenotypic information. Twenty-seven QTLs associated with growth (body weight or body length) were identified and found to be distributed on six LGs (LG1, LG2, LG5, LG7, LG21 and LG24). Interestingly, 21 out of the 27 QTLs were concentrated on LG5 within a narrow region (15.1–33.3 cM). These had the highest LOD value (4.7), accounting for 14.5%–29.1% of the phenotypic variation. The small physical and genetic distance among QTLs within the cluster suggested that the cluster would be highly effective for future marker-assisted selection.
Target genes associated with the growth traits of orange-spotted grouper have been previously studied [31]. We compared our QTLs with the scaffold assembly and annotation of the orange-spotted grouper reference genome [32], and discovered a total of 17 genes within the QTL regions. Three of these genes (notch1, ftz-f1 and tert) have been reported to be involved in fish growth [33,34,35]. The notch1 signaling gene plays a role in the notochord development of zebrafish [32]. The development of the notochord is strongly linked with body length. Growth was correlated with gonadal development and sex change in groupers [36]. In tilapia, ftz-f1 is involved in the development of adrenal-gonadal and sex determination and its transcripts were only expressed in the gonads and kidneys [34]. Phylogenetic analysis of the tilapia ftz-f1 indicated that it was highly conserved among other teleosts. In addition, the orange-spotted grouper is a protogynous hermaphroditic and the expression of ftz-f1 in gonads is influenced by the development of the gonad [37]. The gene of tert is associated with the physiological aging of teleosts and plays a major part in the cell process of proliferation, differentiation and tumorigenesis [35]. Expression of tert was also significantly correlated with muscle telomerase activity (TA) in the skeletal muscle of many fish species [38]. TA has a major effect on the growth rate of orange-spotted grouper. Fourteen additional genes (Table S4) were also mapped to growth-related QTLs, but their exact roles with regard to fish growth remain to be determined.

4. Materials and Methods

4.1. Sample Preparation

Parent fishes (orange-spotted grouper) were captured from the South China Sea in Hainan. One female and one male, with desirable properties such as the coefficient of mature (stage IV) gonad, genetically diverse and highly heterozygous, high vigor, etc., were selected. An F1 full-sib family was generated by crossing the female and male fish at the Daya Bay Seawater Fish Farm in Huizhou, Guangdong Province, China. The F1 progeny were raised under a natural photoperiod, at the water temperature 28 °C. They were fed according to management practices of the fishery. Fin clips of the parents and 68 offspring (four months of age) were collected and stored in absolute ethanol at −20 °C before use. Phenotypic information on the F1 progeny is summarized in Table S5. Phenotypic correlations were tested with the Pearson correlation coefficient, and the results indicated a significant correlation (p < 0.01) between body weight and body length (Table S6). Genomic DNA was extracted using a standard phenol-chloroform protocol [39]. The quality of all DNA samples was evaluated by Qubit Fluorometer (Invitron, Waltham, MA, USA). Electrophoresis was conducted on 0.6% agarose gels. All the experiments were carried out in accordance with the guidelines of the Animal Ethics Committee and were approved by the Institutional Review Board on Bioethics and Biosafety of BGI (No. FT14015, 12 March 2014).

4.2. ddRAD Library Construction and Sequencing

The ddRADseq library was prepared using a previously reported protocol [16] with some modifications. Briefly, the workflow included double-digestion, ligation reaction, pooling, purification, and amplification. The double-digestion and ligation reactions were prepared as 30- and 40-µL reactions (see more details in Tables S7 and S8), respectively. One µL of the ligation products from each sample was pooled 24 samples in a new centrifuge tube, respectively. The DNA fragments were purified (300–700 bp) were isolated on a 3% agarose gel, and then extracted with the QIA quick Gel Extraction Kit (Qiagen, Hilden, Germany). All the samples were amplified by 12 cycles of PCR in the 50-µL reaction volumes containing 2-µL purified of the PCR products, 20-µL Master Mix, 2-µL primers (1-µL each) and 26-µL ddH2O. The amplified products were purified using the QIA quick PCR Purification Kit (Qiagen) and checked the quality (concentration and purity) of DNA products via an Agilent 2100 Bioanalyzer (Agilent, Santa Clara, CA, USA). Finally, library sequencing was performed using a Hiseq 2000 platform (Illumina, San Diego, CA, USA) with 90-bp pair-end reads.

4.3. SNP (Single-Nucleotide Polymorphism) Calling and Genotyping

We discarded the Illumina short reads without sample-specific barcodes and restriction enzyme motifs. Subsequently, we used Soapsnp [40] to sort retained reads into loci and to genotype them. Lastly, we aligned the clean reads of the two parents to the previously assembled orange-spotted grouper genome [32] to eliminate monomorphic DNA sequences and then added the reads containing SNPs. SNPs were selected as described previously [41]. Regions with these putative SNPs were defined as reference SNP regions, and then the clean reads of progeny samples were aligned to the reference SNP regions. As noted earlier, the parental genotypes determined the genotypes of the individual offspring [42].

4.4. High-Density Genetic Map Construction

The SNPs that at least one parent was heterozygous and had high quality of genotype calls were retained for genetic linkage analysis. Those SNPs with significant segregation distortion (p < 0.01, x2 test) were discarded. We used JoinMap 4.1 (Kyazma B.V., Wageningen, Gelderland, Netherlands) to conduct the genetic map construction with the genotypic data from the F1 mapping population. The logarithm of odds (LOD) threshold was set as from two to 15. It was adopted as indicator to cluster analysis using regression mapping algorithm. Map distances were calculated using the Kosambi mapping function in centiMorgans (cM). All the SNPs were arranged and grouped into 24 linkage groups (LGs). Graphics of the linkage groups were generated using custom perl script.

4.5. QTL Mapping

We used WinQTLCart2.5 software (Wang S., Basten C.J. and Zeng Z.-B., Raleigh N.C., USA) to perform the QTL analyses with the method of composite interval mapping (CIM) [43]. The CIM method was executed using Model 6 with parameters of five control markers, a 10-cM window size and backward regression. We used 1000 permutations and the p-value (0.05) of whole genome-wide significance to determine the threshold level of likelihood ratio (LR) (11.5) or LOD (2.5). The LR peak location and surrounding district of every QTL was determined. Phenotypic variation resulting from growth-related QTLs was also calculated by the software.

4.6. Identification of the Genes on Growth Traits Associated QTLs

The genes on growth-associated QTLs were identified via following steps. First we identified the orange-spotted grouper genome scaffolds [32] from these QTL regions. Subsequently, we mapped the corresponding gene identities (IDs) based on their position on the scaffolds. Finally, we retrieved the homologous genes from the orange-spotted grouper gene annotation file [32].

5. Conclusions

A high-density genetic linkage map with a total genetic distance of 1231.98 cM was constructed for the orange-spotted grouper using the ddRADseq method. Twenty-seven growth-related QTLs were identified by QTL mapping and corresponding 17 genes were discovered. These results will provide a valuable resource in future marker-assisted molecular breeding studies of this economically important fish species.

Supplementary Materials

Supplementary materials can be found at https://0-www-mdpi-com.brum.beds.ac.uk/1422-0067/17/4/501/s1.

Acknowledgments

This work was financially supported by China 863 project (No. 2012AA10A407), Nature Science Foundation of China (31370047), Special Project on the Integration of Industry, Education and Research of Guangdong Province (No. 2013B090800017), Shenzhen Special Program for Future Industrial Development (No. JSGG20141020113728803), Shenzhen Special Program for Strategic Emerging Industry Development (HY20130205008), Shenzhen Dapeng Special Program for Industrial Development (No. QT20140101), Guangdong Provincial Science and Technology Program (2014B090905008), Special Fund for Fisheries-Scientific Research of Guangdong Province (A201400A01, A201501A03, A201501A09), and State Key Laboratory of Agricultural Genomics, BGI-Shenzhen, Shenzhen, China (No. 2011DQ782025).

Author Contributions

Qiong Shi, Xinxin You, Yong Zhang and Hao-Ran Lin contributed study design; Ling Xiao, Zining Meng, Haifa Zhang and Yong Zhang contributed to the fish and sample preparation; Hui Yu, Jia Li, Hankui Liu and Xinxin You performed the bioinformatics analysis; Hui Yu, Xinxin You and Qiong Shi wrote the paper. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nelson, J.S. Fishes of the World; John Wiley and Sons, Inc.: New York, NY, USA, 1994. [Google Scholar]
  2. Wang, Y.D.; Huang, S.J.; Chou, H.N.; Liao, W.L.; Gong, H.Y.; Chen, J.Y. Transcriptome analysis of the effect of Vibrio alginolyticus infection on the innate immunity-related complement pathway in Epinephelus coioides. BMC Genom. 2014, 15. [Google Scholar] [CrossRef] [PubMed]
  3. Yue, G.H. Recent advances of genome mapping and marker-assisted selection in aquaculture. Fish Fish. 2014, 15, 376–396. [Google Scholar] [CrossRef]
  4. Eshel, O.; Shirak, A.; Weller, J.I.; Hulata, G.; Ron, M. Linkage and physical mapping of sex region on LG23 of nile tilapia (Oreochromis niloticus). G3 Bethesda 2012, 2, 35–42. [Google Scholar] [CrossRef] [PubMed]
  5. Tsai, H.Y.; Hamilton, A.; Guy, D.R.; Tinch, A.E.; Bishop, S.C.; Houston, R.D. The genetic architecture of growth and fillet traits in farmed Atlantic salmon (Salmo salar). BMC Genet. 2015, 16. [Google Scholar] [CrossRef] [PubMed]
  6. Liu, S.; Vallejo, R.L.; Gao, G.; Palti, Y.; Weber, G.M.; Hernandez, A.; Rexroad, C.E., 3rd. Identification of single-nucleotide polymorphism markers associated with cortisol response to crowding in rainbow trout. Mar. Biotechnol. 2015, 17, 328–337. [Google Scholar] [CrossRef] [PubMed]
  7. Geng, X.; Sha, J.; Liu, S.; Bao, L.; Zhang, J.; Wang, R.; Yao, J.; Li, C.; Feng, J.; Sun, F.; et al. A genome-wide association study in catfish reveals the presence of functional hubs of related genes within QTLs for columnaris disease resistance. BMC Genom. 2015, 16. [Google Scholar] [CrossRef] [PubMed]
  8. Liu, Q.; Sakamoto, T.; Kubota, S.; Okamoto, N.; Yamashita, H.; Takagi, M.; Shigenobu, Y.; Sugaya, T.; Nakamura, Y.; Sano, M.; et al. A genetic linkage map of kelp grouper (Epinephelus bruneus) based on microsatellite markers. Aquaculture 2013, 414–415, 63–81. [Google Scholar] [CrossRef]
  9. Dor, L.; Shirak, A.; Gorshkov, S.; Band, M.R.; Korol, A.; Ronin, Y.; Curzon, A.; Hulata, G.; Seroussi, E.; Ron, M. Construction of a microsatellites-based linkage map for the white grouper (Epinephelus aeneus). G3 Bethesda 2014, 4, 1455–1464. [Google Scholar] [CrossRef] [PubMed]
  10. Bourgeois, Y.X.; Lhuillier, E.; Cezard, T.; Bertrand, J.A.; Delahaie, B.; Cornuault, J.; Duval, T.; Bouchez, O.; Mila, B.; Thebaud, C. Mass production of SNP markers in a nonmodel passerine bird through RAD sequencing and contig mapping to the zebra finch genome. Mol. Ecol. Resour. 2013, 13, 899–907. [Google Scholar] [CrossRef] [PubMed]
  11. Stolting, K.N.; Nipper, R.; Lindtke, D.; Caseys, C.; Waeber, S.; Castiglione, S.; Lexer, C. Genomic scan for single nucleotide polymorphisms reveals patterns of divergence and gene flow between ecologically divergent species. Mol. Ecol. 2013, 22, 842–855. [Google Scholar] [CrossRef] [PubMed]
  12. Wang, J.; Luo, M.C.; Chen, Z.; You, F.M.; Wei, Y.; Zheng, Y.; Dvorak, J. Aegilops tauschii single nucleotide polymorphisms shed light on the origins of wheat D-genome genetic diversity and pinpoint the geographic origin of hexaploid wheat. New Phytol. 2013, 198, 925–937. [Google Scholar] [CrossRef] [PubMed]
  13. Davey, J.W.; Hohenlohe, P.A.; Etter, P.D.; Boone, J.Q.; Catchen, J.M.; Blaxter, M.L. Genome-wide genetic marker discovery and genotyping using next-generation sequencing. Nat. Rev. Genet. 2011, 12, 499–510. [Google Scholar] [CrossRef] [PubMed]
  14. You, X.; Shu, L.; Li, S.; Chen, J.; Luo, J.; Lu, J.; Mu, Q.; Bai, J.; Xia, Q.; Chen, Q.; et al. Construction of high-density genetic linkage maps for orange-spotted grouper Epinephelus coioides using multiplexed shotgun genotyping. BMC Genet. 2013, 14. [Google Scholar] [CrossRef] [PubMed]
  15. Ding, S.-X.; Wang, S.; Wang, D.; Wang, J. Analysis of the karyotype of Epinephalus coioides. J. Xiamen Univ. Nat. Sci. 2004, 43, 428. [Google Scholar]
  16. Peterson, B.K.; Weber, J.N.; Kay, E.H.; Fisher, H.S.; Hoekstra, H.E. Double digest RADseq: An inexpensive method for de novo SNP discovery and genotyping in model and non-model species. PLoS ONE 2012, 7, e37135. [Google Scholar] [CrossRef] [PubMed]
  17. Kai, W.; Nomura, K.; Fujiwara, A.; Nakamura, Y.; Yasuike, M.; Ojima, N.; Masaoka, T.; Ozaki, A.; Kazeto, Y.; Gen, K.; et al. A ddRAD-based genetic map and its integration with the genome assembly of Japanese eel (Anguilla japonica) provides insights into genome evolution after the teleost-specific genome duplication. BMC Genom. 2014, 15. [Google Scholar] [CrossRef] [PubMed]
  18. Zhou, X.; Xia, Y.; Ren, X.; Chen, Y.; Huang, L.; Huang, S.; Liao, B.; Lei, Y.; Yan, L.; Jiang, H. Construction of a SNP-based genetic linkage map in cultivated peanut based on large scale marker development using next-generation double-digest restriction-site-associated DNA sequencing (ddRADseq). BMC Genom. 2014, 15. [Google Scholar] [CrossRef] [PubMed]
  19. Pukk, L.; Ahmad, F.; Hasan, S.; Kisand, V.; Gross, R.; Vasemagi, A. Less is more: Extreme genome complexity reduction with ddRAD using Ion Torrent semiconductor technology. Mol. Ecol. Resour. 2015, 15, 1145–1152. [Google Scholar] [CrossRef] [PubMed]
  20. Alon, S.; Vigneault, F.; Eminaga, S.; Christodoulou, D.C.; Seidman, J.G.; Church, G.M.; Eisenberg, E. Barcoding bias in high-throughput multiplex sequencing of miRNA. Genome Res. 2011, 21, 1506–1511. [Google Scholar] [CrossRef] [PubMed][Green Version]
  21. Craig, D.W.; Pearson, J.V.; Szelinger, S.; Sekar, A.; Redman, M.; Corneveaux, J.J.; Pawlowski, T.L.; Laub, T.; Nunn, G.; Stephan, D.A.; et al. Identification of genetic variants using bar-coded multiplexed sequencing. Nat. Methods 2008, 5, 887–893. [Google Scholar] [CrossRef] [PubMed]
  22. Andolfatto, P.; Davison, D.; Erezyilmaz, D.; Hu, T.T.; Mast, J.; Sunayama-Morita, T.; Stern, D.L. Multiplexed shotgun genotyping for rapid and efficient genetic mapping. Genome Res. 2011, 21, 610–617. [Google Scholar] [CrossRef] [PubMed]
  23. Palaiokostas, C.; Bekaert, M.; Khan, M.G.; Taggart, J.B.; Gharbi, K.; McAndrew, B.J.; Penman, D.J. A novel sex-determining QTL in Nile tilapia (Oreochromis niloticus). BMC Genom. 2015, 16. [Google Scholar] [CrossRef] [PubMed]
  24. Recknagel, H.; Elmer, K.R.; Meyer, A. A hybrid genetic linkage map of two ecologically and morphologically divergent Midas cichlid fishes (Amphilophus spp.) obtained by massively parallel DNA sequencing (ddRADSeq). G3 Bethesda 2013, 3, 65–74. [Google Scholar] [CrossRef] [PubMed]
  25. Franchini, P.; Fruciano, C.; Spreitzer, M.L.; Jones, J.C.; Elmer, K.R.; Henning, F.; Meyer, A. Genomic architecture of ecologically divergent body shape in a pair of sympatric crater lake cichlid fishes. Mol. Ecol. 2014, 23, 1828–1845. [Google Scholar] [CrossRef] [PubMed]
  26. Sharma, R.; Goossens, B.; Kun-Rodrigues, C.; Teixeira, T.; Othman, N.; Boone, J.Q.; Jue, N.K.; Obergfell, C.; O’Neill, R.J.; Chikhi, L. Two different high throughput sequencing approaches identify thousands of de novo genomic markers for the genetically depleted Bornean elephant. PLoS ONE 2012, 7, e49533. [Google Scholar] [CrossRef] [PubMed]
  27. Everett, M.V.; Seeb, J.E. Detection and mapping of QTL for temperature tolerance and body size in Chinook salmon (Oncorhynchus tshawytscha) using genotyping by sequencing. Evol. Appl. 2014, 7, 480–492. [Google Scholar] [CrossRef] [PubMed]
  28. Moghadam, H.K.; Poissant, J.; Fotherby, H.; Haidle, L.; Ferguson, M.M.; Danzmann, R.G. 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. Genom. 2007, 277, 647–661. [Google Scholar] [CrossRef] [PubMed]
  29. Kuttner, E.; Moghadam, H.K.; Skulason, S.; Danzmann, R.G.; Ferguson, M.M. Genetic architecture of body weight, condition factor and age of sexual maturation in Icelandic Arctic charr (Salvelinus alpinus). Mol. Genet. Genom. 2011, 286, 67–79. [Google Scholar] [CrossRef] [PubMed]
  30. Kessuwan, K.; Kubota, S.; Liu, Q.; Sano, M.; Okamoto, N.; Sakamoto, T.; Yamashita, H.; Nakamura, Y.; Ozaki, A. Detection of growth-related quantitative trait loci and high-resolution genetic linkage maps using simple sequence repeat markers in the kelp grouper (Epinephelus bruneus). Mar. Biotechnol. 2016, 18, 57–84. [Google Scholar] [CrossRef] [PubMed]
  31. Guo, L.; Xia, J.; Yang, S.; Li, M.; You, X.; Meng, Z.; Lin, H. GHRH, PRP-PACAP and GHRHR target sequencing via an Ion Torrent personal genome machine reveals an association with growth in orange-spotted grouper (Epinephelus coioides). Int. J. Mol. Sci. 2015, 16, 26137–26150. [Google Scholar] [CrossRef] [PubMed]
  32. Zhang, Y. Sun Yat-Sen University. Personal data. 2016.
  33. Akerberg, A.A.; Stewart, S.; Stankunas, K. Spatial and temporal control of transgene expression in zebrafish. PLoS ONE 2014, 9, e92217. [Google Scholar] [CrossRef] [PubMed]
  34. Cao, J.; Chen, J.; Jiang, Z.; Luo, Y.; Gan, X. Molecular cloning and expression analysis of FTZ-F1 in the GIFT tilapia, Oreochromis niloticus. Zool. Res. 2012, 33, E40–E46. [Google Scholar] [PubMed]
  35. Lopez de Abechuco, E.; Bilbao, E.; Soto, M.; Diez, G. Molecular cloning and measurement of telomerase reverse transcriptase (TERT) transcription patterns in tissues of European hake (Merluccius merluccius) and Atlantic cod (Gadus morhua) during aging. Gene 2014, 541, 8–18. [Google Scholar] [CrossRef] [PubMed]
  36. Renan, X.; Trejo-Martinez, J.; Caballero-Arango, D.; Brule, T. Growth stanzas in an Epinephelidae–Lutjanidae complex: Considerations to length-weight relationships. Rev. Biol. Trop. 2015, 63, 175–187. [Google Scholar] [CrossRef] [PubMed]
  37. Zhang, W.; Zhang, Y.; Zhang, L.; Zhao, H.; Li, X.; Huang, H.; Lin, H. The mRNA expression of P450 aromatase, gonadotropin beta-subunits and FTZ-F1 in the orange-spotted grouper (Epinephelus Coioides) during 17α-methyltestosterone-induced precocious sex change. Mol. Reprod. Dev. 2007, 74, 665–673. [Google Scholar] [CrossRef] [PubMed]
  38. Peterson, D.R.; Mok, H.O.; Au, D.W. Modulation of telomerase activity in fish muscle by biological and environmental factors. Comp. Biochem. Physiol. Toxicol. Pharmacol. 2015, 178, 51–59. [Google Scholar] [CrossRef] [PubMed]
  39. Song, W.; Pang, R.; Niu, Y.; Gao, F.; Zhao, Y.; Zhang, J.; Sun, J.; Shao, C.; Liao, X.; Wang, L.; et al. Construction of high-density genetic linkage maps and mapping of growth-related quantitative trail loci in the Japanese flounder (Paralichthys olivaceus). PLoS ONE 2012, 7, e50404. [Google Scholar] [CrossRef] [PubMed]
  40. Li, R.; Li, Y.; Fang, X.; Yang, H.; Wang, J.; Kristiansen, K.; Wang, J. SNP detection for massively parallel whole-genome resequencing. Genome Res. 2009, 19, 1124–1132. [Google Scholar] [CrossRef] [PubMed][Green Version]
  41. Ao, J.; Li, J.; You, X.; Mu, Y.; Ding, Y.; Mao, K.; Bian, C.; Mu, P.; Shi, Q.; Chen, X. Construction of the high-density genetic linkage map and chromosome map of large yellow croaker (Larimichthys crocea). Int. J. Mol. Sci. 2015, 16, 26237–26248. [Google Scholar] [CrossRef] [PubMed]
  42. Shao, C.; Niu, Y.; Rastas, P.; Liu, Y.; Xie, Z.; Li, H.; Wang, L.; Jiang, Y.; Tai, S.; Tian, Y.; et al. Genome-wide SNP identification for the construction of a high-resolution genetic map of Japanese flounder (Paralichthys olivaceus): Applications to QTL mapping of Vibrio anguillarum disease resistance and comparative genomic analysis. DNA Res. 2015, 22, 161–170. [Google Scholar] [CrossRef] [PubMed]
  43. Silva Lda, C.; Wang, S.; Zeng, Z.B. Composite interval mapping and multiple interval mapping: Procedures and guidelines for using Windows QTL Cartographer. Methods Mol. Biol. 2012, 871, 75–119. [Google Scholar] [PubMed]
  44. Mei, Q.; Liu, J.; Liu, Y.; Li, C.; Wang, H.; Li, H.; Chen, X.; Lan, X. Expression of proline-rich coiled-coil 2B protein in developing rat brains. Neurosci. Lett. 2013. [Google Scholar] [CrossRef] [PubMed]
  45. Su, J.; Ekman, C.; Oskolkov, N.; Lahti, L.; Strom, K.; Brazma, A.; Groop, L.; Rung, J.; Hansson, O. A novel atlas of gene expression in human skeletal muscle reveals molecular changes associated with aging. Skelet. Muscle 2015, 5. [Google Scholar] [CrossRef] [PubMed]
  46. Dai, Z.; Aryal, U.K.; Shukla, A.; Qian, W.J.; Smith, R.D.; Magnuson, J.K.; Adney, W.S.; Beckham, G.T.; Brunecky, R.; Himmel, M.E.; et al. Impact of alg3 gene deletion on growth, development, pigment production, protein secretion, and functions of recombinant Trichoderma reesei cellobiohydrolases in Aspergillus niger. Fungal Genet. Biol. 2013, 61, 120–132. [Google Scholar] [CrossRef] [PubMed]
  47. Pan, H.; Mzhavia, N.; Devi, L.A. Endothelin converting enzyme-2: A processing enzyme involved in the generation of novel neuropeptides. Protein Pept. Lett. 2004, 11, 461–469. [Google Scholar] [CrossRef] [PubMed]
  48. Miners, J.S.; Palmer, J.C.; Love, S. Pathophysiology of hypoperfusion of the precuneus in early Alzheimer’s disease. Brain Pathol. 2015. [Google Scholar] [CrossRef] [PubMed]
  49. Yanagisawa, H.; Hammer, R.E.; Richardson, J.A.; Emoto, N.; Williams, S.C.; Takeda, S.; Clouthier, D.E.; Yanagisawa, M. Disruption of ECE-1 and ECE-2 reveals a role for endothelin-converting enzyme-2 in murine cardiac development. J. Clin. Investig. 2000, 105, 1373–1382. [Google Scholar] [CrossRef] [PubMed]
  50. Sim, K.; Chan, W.Y.; Woon, P.S.; Low, H.Q.; Lim, L.; Yang, G.L.; Lee, J.; Chong, S.A.; Sitoh, Y.Y.; Chan, Y.H.; et al. ARVCF genetic influences on neurocognitive and neuroanatomical intermediate phenotypes in Chinese patients with schizophrenia. J. Clin. Psychiatry 2012, 73, 320–326. [Google Scholar] [CrossRef] [PubMed]
  51. Xu, Z.Y.; Xiong, Y.Z.; Lei, M.G.; Li, F.E.; Zuo, B. Genetic polymorphisms and preliminary association analysis with production traits of the porcine SLC27A4 gene. Mol. Biol. Rep. 2009, 36, 1427–1432. [Google Scholar] [CrossRef] [PubMed]
  52. Tactacan, C.M.; Phua, Y.W.; Liu, L.; Zhang, L.; Humphrey, E.S.; Cowley, M.; Pinese, M.; Biankin, A.V.; Daly, R.J. The pseudokinase SgK223 promotes invasion of pancreatic ductal epithelial cells through JAK1/Stat3 signaling. Mol. Cancer 2015, 14. [Google Scholar] [CrossRef] [PubMed]
  53. Sugawara, K.; Fujikawa, M.; Yoshida, M. Screening of protein kinase inhibitors and knockdown experiments identified four kinases that affect mitochondrial ATP synthesis activity. FEBS Lett. 2013, 587, 3843–3847. [Google Scholar] [CrossRef] [PubMed]
  54. Chaki, S.; Shimazaki, T.; Nishiguchi, M.; Funakoshi, T.; Iijima, M.; Ito, A.; Kanuma, K.; Sekiguchi, Y. Antidepressant/anxiolytic potential and adverse effect liabilities of melanin-concentrating hormone receptor 1 antagonists in animal models. Pharmacol. Biochem. Behav. 2015, 135, 154–168. [Google Scholar] [CrossRef] [PubMed]
  55. Sherwood, A.; Holland, P.C.; Adamantidis, A.; Johnson, A.W. Deletion of Melanin Concentrating Hormone Receptor-1 disrupts overeating in the presence of food cues. Physiol. Behav. 2015, 152, 402–407. [Google Scholar] [CrossRef] [PubMed]
  56. Green, T.; Chen, X.; Ryan, S.; Asch, A.S.; Ruiz-Echevarria, M.J. TMEFF2 and SARDH cooperate to modulate one-carbon metabolism and invasion of prostate cancer cells. Prostate 2013, 73, 1561–1575. [Google Scholar] [CrossRef] [PubMed]
  57. Conover, C.A.; Bale, L.K.; Marler, R.J. Pregnancy-associated plasma protein-A deficiency improves survival of mice on a high fat diet. Exp. Gerontol. 2015, 70, 131–134. [Google Scholar] [CrossRef] [PubMed]
  58. Cheuk, Q.K.; Lo, T.K.; Wong, S.F.; Lee, C.P. Association between pregnancy-associated plasma protein-A levels in the first trimester and gestational diabetes mellitus in Chinese women. Hong Kong Med. J. 2016, 22, 30–38. [Google Scholar] [PubMed]
  59. Hancerliogullari, N.; Aktulay, A.; Engin-Ustun, Y.; Ozkan, M.S.; Oksuzoglu, A.; Danisman, N. Pregnancy-associated plasma protein a levels are decreased in obstetric cholestasis. Clin. Exp. Obstet. Gynecol. 2015, 42, 617–618. [Google Scholar] [PubMed]
  60. Wu, N.L.; Huang, D.Y.; Wang, L.F.; Kannagi, R.; Fan, Y.C.; Lin, W.W. Spleen tyrosine kinase mediates EGFR signaling to regulate keratinocyte terminal differentiation. J. Investig. Dermatol. 2016, 136, 192–201. [Google Scholar] [CrossRef] [PubMed]
  61. Lorenz, J.; Waldschmidt, J.; Wider, D.; Follo, M.; Ihorst, G.; Chatterjee, M.; May, A.M.; Duyster, J.; Rosenwald, A.; Wasch, R.; et al. From CLL to Multiple Myeloma—Spleen Tyrosine Kinase (SYK) influences multiple myeloma cell survival and migration. Br. J. Haematol. 2015. [Google Scholar] [CrossRef] [PubMed]
  62. Nakano, K.; Ando, H.; Kurokawa, S.; Hosohata, K.; Ushijima, K.; Takada, M.; Tateishi, M.; Yonezawa, A.; Masuda, S.; Matsubara, K.; et al. Association of decreased mRNA expression of multidrug and toxin extrusion protein 1 in peripheral blood cells with the development of flutamide-induced liver injury. Cancer Chemother. Pharmacol. 2015, 75, 1191–1197. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The genetic location of growth-related (body weight and body length) quantitative trait loci (QTLs) on the 24 linkage groups (LGs) of orange-spotted grouper. The green orizontal lines indicate the cutting threshold of logarithm of odds (LOD) at 2.5. The red and blue coloration represents the traits of body weight and body length, respectively.
Figure 1. The genetic location of growth-related (body weight and body length) quantitative trait loci (QTLs) on the 24 linkage groups (LGs) of orange-spotted grouper. The green orizontal lines indicate the cutting threshold of logarithm of odds (LOD) at 2.5. The red and blue coloration represents the traits of body weight and body length, respectively.
Ijms 17 00501 g001
Table 1. The information of growth-related quantitative trait loci (QTLs) for body weight (BW) and body length (BL) in orange-spotted grouper.
Table 1. The information of growth-related quantitative trait loci (QTLs) for body weight (BW) and body length (BL) in orange-spotted grouper.
QTLTraitGenetic Position(cM)Logarithm of Odds (LOD)R2 (%)
qLG1BW30.8–31.03.319.6
qLG2BL29.8–30.23.321
qLG5_1BW15.1–15.64.724.3
qLG5_2BW17.0–17.53.117.1
qLG5_3BW17.0–17.52.816.5
qLG5_4BW17.8–18.02.915.8
qLG5_5BW18.7–19.43.219.1
qLG5_6BW18.7–19.44.224.6
qLG5_7BW19.9–20.23.520.8
qLG5_8BW20.4–20.83.420
qLG5_9BW27.4–27.92.621.5
qLG5_10BW27.4–27.92.824.5
qLG5_11BW28.1–28.32.828.3
qLG5_12BW29.4–29.72.714.6
qLG5_13BW29.4–29.7329.1
qLG5_14BW32.2–32.43.217.6
qLG5_15BW33.0–33.32.915.6
qLG5_16BL15.1–15.44.423.3
qLG5_17BL17.0–17.12.614.5
qLG5_18BL19.2–19.33.820.5
qLG5_19BL19.9–20.23.420.3
qLG5_20BL19.9–20.23.822.2
qLG5_21BL20.4–20.83.520.8
qLG7_1BW14.4–14.52.816.8
qLG7_2BL14.4–14.52.615.1
qLG21BL16.0–16.52.915.5
qLG24BW48.8–49.02.614.4
Table 2. Genes in the growth-related QTL of orange-spotted grouper.
Table 2. Genes in the growth-related QTL of orange-spotted grouper.
QTLGene IDGene DescriptionGene SymbolPhysical Position
qLG1Eco_gene_10012941fasciculation and elongation protein ζ-2-likefez2scaffold416
qLG2Eco_gene_10013595Dol-P-Man: Man(5)GlcNAc(2)-PP-Dolα-1,3-mannosyltransferasealg3scaffold438
Eco_gene_10013597endothelin converting enzyme 2ece2scaffold438
qLG5_1Eco_gene_10014238armadillo repeat gene deleted in velocardiofacial syndromearvcfscaffold468
qLG5_3Eco_gene_10004531solute carrier family 27 (fatty acid transporter), member 4slc27a4scaffold1517
Eco_gene_10005552tyrosine-protein kinase SgK223sgk223scaffold1706
qLG5_4Eco_gene_10002069calcium/calmodulin-dependent protein kinase 2camk2scaffold1180
qLG5_5Eco_gene_10005555proline-rich coiled-coil 2Bprrc2bscaffold1706
Eco_gene_10000893melanin-concentrating hormone receptor 1mchr1scaffold1067
qLG5_19Eco_gene_10005043sarcosine dehydrogenasesardhscaffold1600
qLG5_8Eco_gene_10018102pregnancy-associated plasma protein Apappascaffold675
qLG5_9Eco_gene_10014809spleen tyrosine kinasesykscaffold491
qLG5_11Eco_gene_10002109telomerase reverse transcriptasetertscaffold1187
qLG5_12Eco_gene_10005999WD repeat-containing protein 91-likewdrcp91scaffold1804
Eco_gene_10000470ftz transcription factor 1ftz-f1scaffold1023
qLG7Eco_gene_10021611multidrug and toxin extrusion protein 1-likemate1scaffold919
qLG21Eco_gene_10003094neurogenic locus notch homolog protein 1-likenotch1scaffold130
Back to TopTop