Mapping QTL for an adaptive trait: the length of caudal fin in Lates calcarifer

Genome-Wide Association Study Identifies Loci Associated with Resistance to Viral Nervous Necrosis Disease in Asian Seabass

Viral nervous necrosis disease (VNN), caused by nervous necrosis virus (NNV), is one major threat to mariculture. Identifying loci and understanding the mechanisms associated with resistance to VNN are important in selective breeding programs. We performed a genome-wide association study (GWAS) using genotyping-by-sequencing (GBS) to study the genomic architecture of resistance to NNV infection in Asian seabass. We genotyped 986 individuals from 43 families produced by 15 founders with 44498 bi-allelic genetic variants using GBS. The GWAS identified three genome-wide significant loci on chromosomes 16, 19, and 20, respectively, and six suggestive loci on chromosomes 1, 8, 14, 15, 21, and 24, respectively, associated with resistance to NNV infection measured as binary and quantitative traits. Using the 500 most significant markers in combination with a training population of 800 samples could reach a genomic prediction accuracy of 0.7. Candidate genes significantly associated with resistance to NNV, including lysine-specific demethylase 2A, beta-defensin 1, and cystatin-B, which play important roles in immune responses against virus infection, were identified. Almost all the candidate genes were differentially expressed in different tissues against NNV infection. The significant genetic variants can be used in genomic selection and help understand the mechanism of resistance to VNN. Future studies should use populations of large effective size and whole genome resequencing to identify more useful genetic variants.

RAD-QTL Mapping Reveals Both Genome-Level Parallelism and Different Genetic Architecture Underlying the Evolution of Body Shape in Lake Whitefish (Coregonus clupeaformis) Species Pairs

Parallel changes in body shape may evolve in response to similar environmental conditions, but whether such parallel phenotypic changes share a common genetic basis is still debated. The goal of this study was to assess whether parallel phenotypic changes could be explained by genetic parallelism, multiple genetic routes, or both. We first provide evidence for parallelism in fish shape by using geometric morphometrics among 300 fish representing five species pairs of Lake Whitefish. Using a genetic map comprising 3438 restriction site−associated DNA sequencing single-nucleotide polymorphisms, we then identified quantitative trait loci underlying body shape traits in a backcross family reared in the laboratory. A total of 138 body shape quantitative trait loci were identified in this cross, thus revealing a highly polygenic architecture of body shape in Lake Whitefish. Third, we tested for evidence of genetic parallelism among independent wild populations using both a single-locus method (outlier analysis) and a polygenic approach (analysis of covariation among markers). The single-locus approach provided limited evidence for genetic parallelism. However, the polygenic analysis revealed genetic parallelism for three of the five lakes, which differed from the two other lakes. These results provide evidence for both genetic parallelism and multiple genetic routes underlying parallel phenotypic evolution in fish shape among populations occupying similar ecological niches.

Mapping QTLs for swimming ability related traits in Cyprinus carpio L

Body height (BH), head length (HL), snout length (SL), and tail length (TL) are important traits related with swimming ability of fish. Therefore, improving these traits will increase the production which is the basic goal of aquaculture breeding. To understand the genetic basis of swimming ability related traits in Cyprinus carpio L., a high-density linkage map spanning 3,301 cM in 50 linkage groups was utilized for quantitative trait locus (QTL) mapping. Mapping family comprised 190 offspring and 627 molecular markers were genotyped with average distance of 5.6 cM. A total of 15 QTLs including four (qBH13, qBH30, qBH33, qBH48) for BH, four (qHL10, qHL18, qHL29, qHL48) for HL, three (qSL24, qSL27, qSL45) for SL, and four (qTL15, qTL17, qTL18, qTL44) for TL were detected on 13 linkage groups LG10, LG13, LG15, LG17, LG18, LG24, LG27, LG29, LG30, LG33, LG44, LG45, and LG48. Each LG consisted on single QTL except LG18 and LG48. LG18 was found with two QTLs associated with HL and TL. While LG48 was comprised, the QTLs related with BH and HL. The phenotype variance was recorded from 12.6 to 40.6 %. Five QTLs, qHL48, qSL45, qTL15, qTL18, and qTL44, explained phenotype variance of >20 % with a significant levels of 0.047, 0.049, 0.037, 0.025, and 0.023, respectively. The neighbored loci of these QTLs were considered as main region of chromosomes controlling the traits. These identified genetic regions will be the main source of discovering gene(s) associated with swimming ability related traits in C. carpio L.

Molecular characterization and mapping of Fgf21 gene in a foodfish species asian seabass

Fgf21 is a newly discovered fibroblast growth factor. It is typically induced by fasting and plays important roles in the regulation of glucose and lipid metabolisms and energy balance in mammals, whereas potential functions of this gene in teleosts are still unknown. We identified the Fgf21 gene and studied its functions in Asian seabass (Lates calcarifer). The cDNA of the Fgf21 encoded a protein with 206 amino acids. Analysis of DNA and amino acid sequences of Fgf21 genes revealed that the sequences and structure of the Fgf21 genes were highly conserved in vertebrates. Real-time PCR revealed that Fgf21 was exclusively expressed in the intestine and kidney, which was different from the expression profiles of mammals. Fgf21 was down-regulated under fasting, whereas it was significantly increased during the LPS challenge. Exogenous recombinant FGF21 significantly suppressed the appetite of Asian seabass. Our data suggest that Fgf21 plays a role in energy regulation and acute phase response in Asian seabass, and may have different functions in fish and mammals. In addition, we identified one SNP in Fgf21. By using this SNP, the gene was mapped on the linkage group 23, where a suggestive QTL for growth was mapped previously. Association mapping identified significant associations between Fgf21 genotypes at the SNP and growth traits. These results not only provide important information of the functions of Fgf21, but also suggest that the SNP in this gene can be used as a marker in selecting fast-growing individuals of Asian seabass.

Mapping quantitative trait loci for omega-3 fatty acids in Asian seabass

Omega-3 fatty acids are essential fatty acids for human health. Therefore, increasing both percentage of omega-3 and a better fatty acid profile in fish fillets is one of the breeding goals in aquaculture. However, it is difficult to increase the omega-3 content in fish fillets, as the phenotypic selection of these traits is not easily feasible. To facilitate the genetic improvement of the Asian seabass for optimal fatty acid profiles, a genome-wide scan for quantitative trait loci (QTL) affecting fatty acid level in the flesh of the Asian seabass was performed on an F2 family containing 314 offspring. All family members were genotyped using 123 informative microsatellites and 22 SNPs. High percentages of n-3 polyunsaturated fatty acids (PUFA), especially C22:6 (DHA 16.48 ± 3.09 %) and C20:5 (EPA 7.19 ± 0.86 %) were detected in the flesh. One significant and 54 suggestive QTL for different fatty acids and a water content trait were detected on the whole genome. QTL for C18:0b was located on linkage groups (LG) 5. QTL for total n-3 PUFA content in flesh were mapped onto LG6 and LG23 with the phenotypic variance explained ranging from 3.8 to 6.3 %. Four QTL for C22:6 were detected on LG6, LG23, and LG24, explaining 3.9 to 4.9 % of the phenotypic variance, respectively. Mapping of QTL for contents of different fatty acids is the first step towards improving the omega-3 content in the fillets of fish by using marker-assisted selection and is important for understanding the biology of fatty acid deposition.

Identification of quantitative trait locus (QTL) linked to dorsal fin length from preliminary linkage map of molly fish, Poecilia sp

A preliminary linkage map was constructed by applying backcross and testcross strategy using microsatellite (SSR) markers developed for Xiphophorus and Poecilia reticulata in ornamental fish, molly Poecilia sp. The linkage map having 18 SSR loci consisted of four linkage groups that spanned a map size of 516.1cM. Association between genotypes and phenotypes was tested in a random fashion and QTL for dorsal fin length was found to be linked to locus Msb069 on linkage group 2. Coincidentally, locus Msb069 was also reported as putative homologue primer pairs containing SSRs repeat motif which encoded hSMP-1, a sex determining locus. Dorsal fin length particularly in males of Poecilia latipinna is an important feature during courtship display. Therefore, we speculate that both dorsal fin length and putative hSMP-1 gene formed a close proximity to male sexual characteristics.

Quantitative trait loci for morphometric body measurements of the hybrids of silver carp (Hypophthalmichthys molitrix) and bighead carp (H. nobilis)

Quantitative trait loci (QTL) for 11 morphometric body measurements of the hybrids of silver (Hypophthalmichthys molitrix) and bighead carp (H. nobilis) including body weight (BW), standard length (SL), body depth (BD), body thickness (BT), head length (HL), head depth (HD), length of ventral keel (LVK), length of pectoral fin (Lpec), length of pelvic fin (Lpel), length of caudal fin (Lcau) and space between pectoral and pelvic fins (SPP) were located on the sex average microsatellite linkage map constructed using the hybrids of a female bighead and a male silver carp, on which 15 microsatellites were newly mapped. One locus was found to be responsible for BW, LV K and SPP, respectively. As many as 6 loci were found to be responsible for HD. The variances of remaining traits were partitioned by different numbers of loci varying between 2 and 5. The variance explained each locus ranged from 9.1% to 23.8% of the total. The variance explained by all loci responsible for each measurement ranged from 17.7% to 75.1%. It was noted that multiple measurements were mapped on the same locus. For example, a region bounded by Hym435 and Hym145 was found to be responsible for all the measurements analyzed.

Whole genome scanning and association mapping identified a significant association between growth and a SNP in the IFABP-a gene of the Asian seabass

BACKGROUND

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 F(2) family with 359 individuals using 123 microsatellites and 22 SNPs, and performed association mapping in four populations with 881 individuals.

RESULTS

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.

CONCLUSIONS

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.

Construction of a high-density microsatellite genetic linkage map and mapping of sexual and growth-related traits in half-smooth tongue sole (Cynoglossus semilaevis)

High-density genetic linkage maps of half-smooth tongue sole were developed with 1007 microsatellite markers, two SCAR markers and an F1 family containing 94. The female map was composed of 828 markers in 21 linkage groups, covering a total of 1447.3 cM, with an average interval 1.83 cM between markers. The male map consisted of 794 markers in 21 linkage groups, spanning 1497.5 cM, with an average interval of 1.96 cM. The female and male maps had 812 and 785 unique positions, respectively. The genome length of half-smooth tongue sole was estimated to be 1527.7 cM for the females and 1582.1 cM for the males. Based on estimations of the map lengths, the female and male maps covered 94.74 and 94.65% of the genome, respectively. The consensus map was composed of 1007 microsatellite markers and two SCAR markers in 21 linkage groups, covering a total of 1624 cM with an average interval of 1.67 cM. Furthermore, 159 sex-linked SSR markers were identified. Five sex-linked microsatellite markers were confirmed in their association with sex in a large number of individuals selected from different families. These sex-linked markers were mapped on the female map LG1f with zero recombination. Two QTLs that were identified for body weight, designated as We-1 and We-2, accounted for 26.39% and 10.60% of the phenotypic variation. Two QTLs for body width, designated Wi-1 and Wi-2, were mapped in LG4f and accounted for 14.33% and 12.83% of the phenotypic variation, respectively. Seven sex-related loci were mapped in LG1f, LG14f and LG1m by CIM, accounting for 12.5–25.2% of the trait variation. The results should prove to be very useful for improving growth traits using molecular MAS.

Tracing Asian seabass individuals to single fish farms using microsatellites

Traceability through physical labels is well established, but it is not highly reliable as physical labels can be easily changed or lost. Application of DNA markers to the traceability of food plays an increasingly important role for consumer protection and confidence building. In this study, we tested the efficiency of 16 polymorphic microsatellites and their combinations for tracing 368 fish to four populations where they originated. Using the maximum likelihood and Bayesian methods, three most efficient microsatellites were required to assign over 95% of fish to the correct populations. Selection of markers based on the assignment score estimated with the software WHICHLOCI was most effective in choosing markers for individual assignment, followed by the selection based on the allele number of individual markers. By combining rapid DNA extraction, and high-throughput genotyping of selected microsatellites, it is possible to conduct routine genetic traceability with high accuracy in Asian seabass.

Construction of an integrated map of Haliotis diversicolor using microsatellite markers

Small abalone, Haliotis diversicolor, is naturally distributed along the coastal waters of East Asia from Japan to the Philippines. It is an economically important maricultured species in southern China and Taiwan. Genetic linkage maps for small abalone were constructed using a total of 308 simple sequence repeat markers including 297 novel markers. Segregation data on 96 progeny were genotyped using a pseudo-testcross strategy. Sixteen linkage groups were identified in both female and male maps, consistent with the haploid chromosome number. The female linkage map covered 758.3 cM, with an average interval of 5.2 cM. The male linkage map spanned a total genetic distance of 676.2 cM, with an average interval of 4.5 cM. An integrated map was constructed by incorporating the homologous parental linkage groups, resulting in 16 linkage groups with a total of 762.1 cM. Genome coverage of the integrated linkage map was approximately 80.7%. The genetic linkage maps of small abalone may facilitate marker-assisted selection and quantitative trait loci mapping.

Analysis of the Asian seabass transcriptome based on expressed sequence tags

Analysis of transcriptomes is of great importance in genomic studies. Asian seabass is an important fish species. A number of genomic tools in it were developed, while large expressed sequence tag (EST) data are lacking. We sequenced ESTs from nine normalized cDNA libraries and obtained 11 431 high-quality ESTs. We retrieved 8524 ESTs from dbEST database and analyzed all 19 975 ESTs using bioinformatics tools. After clustering, we obtained 8837 unique sequences (2838 contigs and 5999 singletons). The average contig length was 574 bp. Annotation of these unique sequences revealed that 48.9% of them showed significant homology to RNA sequences in GenBank. Functional classification of the unique ESTs identified a broad range of genes involved in different functions. We identified 6114 putative single-nucleotide polymorphisms and 634 microsatellites in ESTs. We discovered different temporal and spatial expression patterns of some immune-related genes in the Asian seabass after challenging with a pathogen Vibrio harveyi. The unique EST sequences are being used in developing a cDNA microarray to examine global gene expression and will also facilitate future whole-genome sequence assembly and annotation of Asian seabass and comparative genomics.