A genome scan for quantitative trait loci affecting growth-related traits in an F1 family of Asian seabass (Lates calcarifer)

Chromosome-scale assembly and quantitative trait locus mapping for major economic traits of the Culter alburnus genome using Illumina and PacBio sequencing with Hi-C mapping information

Topmouth culter (Culter alburnus) is an economically important freshwater fish with high nutritional value. However, its potential genetic advantages have not been fully exploited. Therefore, we aimed to determine the genome sequence of C. alburnus and examine quantitative trait loci (QTLs) related to major economic traits. The results showed that 24 pseudochromosomes were anchored by 914.74 Mb of the C. alburnus genome sequence. De novo sequencing identified 31,279 protein-coding genes with an average length of 8507 bp and average coding sequ ence of 1115 bp. In addition, a high-density genetic linkage map consisting of 24 linkage groups was constructed based on 353,532 high-quality single nucleotide polymorphisms and 4,710 bin markers. A total of 28 QTLs corresponding to 11 genes, 26 QTLs corresponding to 11 genes, and 12 QTLs corresponding to 5 genes were identified for sex, intermuscular spine number and body weight traits, respectively. In this study, we assembled an accurate and nearly complete genome of C. alburnus by combining Illumina, PacBio, and high-throughput Chromosome conformation capture (Hi-C) technologies. In addition, we identified QTLs that explained variances in intermuscular spine number, body weight, and sex differences in C. alburnus. These genetic markers or candidate genes associated with growth traits provide a basis for marker-assisted selection in C. alburnus.

A High-Density Genetic Map Enables Genome Synteny and QTL Mapping of Vegetative Growth and Leaf Traits in Gardenia

The gardenia is a traditional medicinal horticultural plant in China, but its molecular genetic research has been largely hysteretic. Here, we constructed an F₁ population with 200 true hybrid individuals. Using the genotyping-by-sequencing method, a high-density sex-average genetic map was generated that contained 4,249 SNPs with a total length of 1956.28 cM and an average genetic distance of 0.46 cM. We developed 17 SNP-based Kompetitive Allele-Specific PCR markers and found that 15 SNPs were successfully genotyped, of which 13 single-nucleotide polymorphism genotypings of 96 F₁ individuals showed genotypes consistent with GBS-mined genotypes. A genomic collinearity analysis between gardenia and the Rubiaceae species Coffea arabica, Coffea canephora and Ophiorrhiza pumila showed the relativity strong conservation of LG11 with NC_039,919.1, HG974438.1 and Bliw01000011.1, respectively. Lastly, a quantitative trait loci analysis at three phenotyping time points (2019, 2020, and 2021) yielded 18 QTLs for growth-related traits and 31 QTLs for leaf-related traits, of which qBSBN7-1, qCD8 and qLNP2-1 could be repeatably detected. Five QTL regions (qCD8 and qSBD8, qBSBN7 and qSI7, qCD4-1 and qLLLS4, qLNP10 and qSLWS10-2, qSBD10 and qLLLS10) with potential pleiotropic effects were also observed. This study provides novel insight into molecular genetic research and could be helpful for further gene cloning and marker-assisted selection for early growth and development traits in the gardenia.

Genome-Wide Association Study of Body Shape-Related Traits in Large Yellow Croaker (Larimichthys crocea)

Large yellow croaker (Larimichthys crocea) is one of the most important cultured marine fish on the southeast coast of China. Its body shape is important for the aquaculture industry since it affects the behavior such as swimming, ingesting, and evading, as well as customer preference. Due to the greater consumer demand of small head, slender body large yellow croaker, selecting and breeding of slender individuals with the assistance of genetic markers will benefit the industry quickly. In this study, several traits were employed to represent body shape, including body depth/body length (BD/BL), body thickness/body length (BT/BL), caudal peduncle depth/caudal peduncle length (CPDLR), tail length/body length (TL/BL), and body area/head area (BA/HA). Genome-wide association study was conducted with a panmictic population of 280 individuals to identify SNP and genes potentially associated with body shape. A set of 20 SNPs on 12 chromosomes were identified to be significantly associated with body shape-related traits. Besides, 5 SNPs were identified to be suggestive associated with CPDLR and BT/BL. Surrounding these SNPs, we found some body shape-related candidate genes, including fabp1, acrv1, bcor, mstn, bambi, and neo1, which involved in lipid metabolism, TGF-β signaling, and BMP pathway and other important regulatory pathways. These results will be useful for the understanding of the genetic basis of body shape formation and helpful for body shape controlling of large yellow croaker by using marker-assisted selection.

Transcriptome analysis of grass carp (Ctenopharyngodon idella) between fast- and slow-growing fish

Grass carp is one of the most important freshwater aquaculture species in China. However, the mechanisms underlying the growth of muscle tissue in the fish are unclear. High-throughput RNA-Seq was used to analyze the transcriptome of grass carp muscle tissue between fast- and slow-growing fish family groups. Twenty-four individuals each from 4 fast-growing families and 4 slow-growing families were used to reduce background noise. 71 up-regulated and 35 down-regulated genes were identified in the differentially expressed genes (DEGs). GO and KEGG enrichment analyses revealed the DEGs were involved in the GH/IGF axis, calcium metabolism, protein and glycogen synthesis, oxygen transport, cytoskeletal and myofibrillar components. IGFBP1 was up-regulated in big fish while GHR2 was down-regulated. Glutamic pyruvate transaminase 2, an indicator of liver tissue damage, was down-regulated in big grass carp, which indicates that the fish was better adapted to an artificially formulated diet. GAPDH, the rate-limiting enzyme in glycolytic flux was highly expressed in fast-growing grass carp, reflecting enhanced carbohydrate metabolism. Higher expression of ALAS2 and myoglobin 1 in big grass carp, related to oxygen transport might promote aerobic exercise along with food intake and muscle growth. Genes for cytoskeletal and myofibrillar components such as tropomyosin, meromyosin, and troponin I were also up-regulated in big grass carp. These results provide valuable information about the key genes for use as biomarkers of growth in selective breeding programs for grass carp and contribute to our understanding of the molecular mechanisms and regulative pathways regulating growth in fish.

High-Density Linkage Map and Mapping for Sex and Growth-Related Traits of Largemouth Bass (Micropterus salmoides)

The largemouth bass is an important species, and its culture has risen sharply with the surge in fish aquaculture in China. Due to the lack of selective breeding technology for the largemouth bass, the growth rate and disease resistance are low, its sexual maturation is slow, and other serious problems are contributing to a sharp decline in the safety and quality of largemouth bass products in recent decades. Therefore, comprehensive breeding programs to improve the economic performance and promote the modern industrial development of largemouth bass must be considered a priority. Here, a total of 152 adult largemouth bass, including two parents and 150 progenies, were selected to produce the genetic mapping family. Then, a high-density linkage map was constructed based on restriction site–associated DNA sequencing using 6,917 single-nucleotide polymorphisms (SNPs) located in 24 linkage groups (LGs). The total genetic length of the linkage map was 1,261.96 cM, and the length of each LG varied from 24.72 cM for LG02 to 117.53 cM for LG16, with an average length of 52.58 cM and an average SNP number of 286. Thirteen significant quantitative trait loci (QTLs) for sex determination were located on LG04, LG05, LG08, LG12, LG15, LG21, and LG23. An informative QTL cluster that included six QTLs was detected on LG12. However, one notable QTL, which accounted for 71.48% of the total phenotypic variation, was located in the region of 1.85 cM on LG05. In addition, 32 identified QTLs were related to growth, including body weight, body length, body height, and head length. The QTLs for these growth-related traits are located in 13 LG regions and have little effect on phenotypic variation. This high-density genetic linkage map will enable the fine-mapping of economic traits and support the future genome assembly of the largemouth bass. Additionally, our study will be useful for future selective culture of largemouth bass and could potentially be used in molecular-assisted breeding of largemouth bass for aquaculture.

Mapping and Validating QTL for Fatty Acid Compositions and Growth Traits in Asian Seabass

Asian seabass is an important food fish species. While improving growth, increasing the nutritional value is important, omega-3 fatty acids are indispensable to human health. Identifying and validating DNA markers associated with traits is the first step towards marker-assisted selection (MAS). We quantified 13 different fatty acids and three growth traits in 213 F₂ Asian seabass from a family at the age 270 days post hatch, and screened QTL for these traits. The content of total fatty acids in 100 g flesh was 2.57 ± 0.80 g, while the proportions of docosahexaenoic acid (DHA) and eicosapentaenoic acid (EPA) were 16.96 ± 2.20% and 5.42 ± 0.90%, respectively. A linkage map with 2424 SNPs was constructed and used for QTL mapping. For fatty acid compositions, 14 significant QTL were identified on three linkage groups (LG5, LG11 and LG14), with phenotypic variance explained (PVE) from 12.8 to 24.6%. Thirty-nine suggestive QTL were detected on 16 LGs. Two significant QTL for EPA were identified on LG5 and LG14, with PVE of 15.2% and 15.1%, respectively. No significant QTL was identified for DHA. For growth traits, six significant and 13 suggestive QTL were identified on two and seven LGs, respectively. Only a few significant QTL for fatty acids overlapped with previously mapped QTL for these traits, suggesting that most QTL detected in a family are family-specific and could only be used in MAS in the family per se. To facilitate population-wide molecular breeding, more powerful methods (e.g. GWAS) should be used to identify SNPs for genomic selection.

QTL Mapping for Red Blotches in Malaysia Red Tilapia (Oreochromis spp.)

Body color is an interesting economic trait in fish. Red tilapia with red blotches may decrease its commercial values. Conventional selection of pure red color lines is a time-consuming and labor-intensive process. To accelerate selection of pure lines through marker-assisted selection, in this study, double-digest restriction site-associated DNA sequencing (ddRAD-seq) technology was applied to genotype a full-sib mapping family of Malaysia red tilapia (Oreochromis spp.) (N = 192). Genome-wide significant quantitative trait locus (QTL)-controlling red blotches were mapped onto two chromosomes (chrLG5 and chrLG15) explaining 9.7% and 8.2% of phenotypic variances by a genome-wide association study (GWAS) and linkage-based QTL mapping. Six SNPs from the chromosome chrLG5 (four), chrLG15 (one), and unplaced supercontig GL831288-1 (one) were significantly associated to the red blotch trait in GWAS analysis. We developed nine microsatellite markers and validated significant correlations between genotypes and blotch data (p < 0.05). Our study laid a foundation for exploring a genetic mechanism of body colors and carrying out genetic improvement for color quality in tilapia.

A Genome-Wide Association Study Reveals That Genes with Functions for Bone Development Are Associated with Body Conformation in Catfish

Body conformation is of great scientific and commercial interest for aquaculture fish species because it affects biological adaptation of the organism to environments, and is of economic importance to the aquaculture industry considering its direct effect on fillet yield. Catfish is the primary aquaculture species in the USA. Two major species used in the aquaculture industry, channel catfish and blue catfish, differ in body shape and therefore the backcross progenies serve as a good model for quantitative trait locus (QTL) analysis. Here, a genome-wide association study (GWAS) with hybrid catfish was conducted to identify the QTL for body conformation, including deheaded body length (DBL), body length (BL), body depth (BD), and body breadth (BB), which were all standardized by cubic root of body weight. Overall, the results indicate that the traits are polygenic. For DBL, linkage group (LG) 2 and LG 24 contain significant QTL, and LG 13 and LG 26 contain suggestively associated QTL (-log₁₀(P value) > 4.5). Compared with DBL, additional SNPs were identified to be associated with body length on LG 2, LG 7, and LG 18. Although no significant QTL for body depth was found, three suggestively associated QTLs were identified on LG 5, LG 13, and LG 14. No SNP for body breadth reached the threshold for suggestive association. Genes close to the associated SNPs were determined, many of which are known to be involved in bone development. This work therefore provides the basis for future identification of causal genes for the control of body conformation.

QTL Mapping for Resistance to Iridovirus in Asian Seabass Using Genotyping-by-Sequencing

Identifying quantitative trait loci (QTL) for viral disease resistance is of particular importance in selective breeding programs of fish species. Genetic markers linked to QTL can be useful in marker-assisted selection (MAS) for elites resistant to specific pathogens. Here, we conducted a genome scan for QTL associated with Singapore grouper iridovirus (SGIV) resistance in an Asian seabass (Lates calcarifer) family, using a high-density linkage map generated with genotyping-by-sequencing. One genome-wide significant and three suggestive QTL were detected at LG21, LG6, LG13, and LG15, respectively. The phenotypic variation explained (PVE) by the four QTL ranged from 7.5 to 15.6%. The position of the most significant QTL at LG21 was located between 31.88 and 36.81 cM. The SNP marker (SNP130416) nearest to the peak of this QTL was significantly associated with SGIV resistance in an unrelated multifamily population. One candidate gene, MECOM, close to the peak of this QTL region, was predicted. Evidence of alternative splicing was observed for MECOM and one specific category of splicing variants was differentially expressed at 5 days post-SGIV infection. The QTL detected in this study are valuable resources and can be used in the selective breeding programs of Asian seabass with regard to resistance to SGIV.

Construction of a high-density linkage map and mapping of sex determination and growth-related loci in the mandarin fish (Siniperca chuatsi)

Background

The mandarin fish (Siniperca chuatsi) is an important and widely cultured fish in China. However, the lack of selective breeding of mandarin fish in previous decades has resulted in a decline in the growth rate of pond-cultured fish, a shortened period of sexual maturity, and reduced disease resistance; these issues seriously affect the quality and safety of the fish products. Therefore, it is necessary to establish a selective breeding program for the mandarin fish to improve the economical traits of the fish and to sustain the development of the mandarin fish industry.

Results

We constructed a high-density linkage map for it based on double digest restriction site associated DNA sequencing (ddRAD-Sequencing). This map contained 3283 dimorphic single nucleotide polymorphism markers and 24 linkage groups (LGs). The total map-length was 1972.01 cM, with an average interlocus distance of 0.61 cM. One significant quantitative trait locus (QTL) for sex determination trait was detected on LG23, which was supported by five markers, clustered between 60.27 and 68.71 cM. The highest logarithm of odds value (17.73) was located at 60.27 cM, near the marker r1_73194, accounting for 53.3% of the phenotypic variance. Genotypes of all the male fish on r1_33008 were homozygous, whereas those of all females were heterozygous. Thus, LG23 was considered a sex-related linkage group. Eleven significant QTLs, for three growth traits, at two growth stages and the increased values were distributed on four LGs; their contributions to the phenotypic variation were quite low (12.4–17.2%), suggesting that multiple genes affected the growth traits.

Conclusion

This high-resolution genetic map provides a valuable resource for fine-mapping of important traits and for identification of sex-related markers that should facilitate breeding of all-female mandarin fish for aquaculture and mechanistic studies on sex determination.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3830-3) contains supplementary material, which is available to authorized users.

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.

Aquaculture genomics, genetics and breeding in the United States: current status, challenges, and priorities for future research

Advancing the production efficiency and profitability of aquaculture is dependent upon the ability to utilize a diverse array of genetic resources. The ultimate goals of aquaculture genomics, genetics and breeding research are to enhance aquaculture production efficiency, sustainability, product quality, and profitability in support of the commercial sector and for the benefit of consumers. In order to achieve these goals, it is important to understand the genomic structure and organization of aquaculture species, and their genomic and phenomic variations, as well as the genetic basis of traits and their interrelationships. In addition, it is also important to understand the mechanisms of regulation and evolutionary conservation at the levels of genome, transcriptome, proteome, epigenome, and systems biology. With genomic information and information between the genomes and phenomes, technologies for marker/causal mutation-assisted selection, genome selection, and genome editing can be developed for applications in aquaculture. A set of genomic tools and resources must be made available including reference genome sequences and their annotations (including coding and non-coding regulatory elements), genome-wide polymorphic markers, efficient genotyping platforms, high-density and high-resolution linkage maps, and transcriptome resources including non-coding transcripts. Genomic and genetic control of important performance and production traits, such as disease resistance, feed conversion efficiency, growth rate, processing yield, behaviour, reproductive characteristics, and tolerance to environmental stressors like low dissolved oxygen, high or low water temperature and salinity, must be understood. QTL need to be identified, validated across strains, lines and populations, and their mechanisms of control understood. Causal gene(s) need to be identified. Genetic and epigenetic regulation of important aquaculture traits need to be determined, and technologies for marker-assisted selection, causal gene/mutation-assisted selection, genome selection, and genome editing using CRISPR and other technologies must be developed, demonstrated with applicability, and application to aquaculture industries.Major progress has been made in aquaculture genomics for dozens of fish and shellfish species including the development of genetic linkage maps, physical maps, microarrays, single nucleotide polymorphism (SNP) arrays, transcriptome databases and various stages of genome reference sequences. This paper provides a general review of the current status, challenges and future research needs of aquaculture genomics, genetics, and breeding, with a focus on major aquaculture species in the United States: catfish, rainbow trout, Atlantic salmon, tilapia, striped bass, oysters, and shrimp. While the overall research priorities and the practical goals are similar across various aquaculture species, the current status in each species should dictate the next priority areas within the species. This paper is an output of the USDA Workshop for Aquaculture Genomics, Genetics, and Breeding held in late March 2016 in Auburn, Alabama, with participants from all parts of the United States.

Transcriptome assembly and identification of genes and SNPs associated with growth traits in largemouth bass (Micropterus salmoides)

Growth is one of the most crucial economic traits of all aquaculture species, but the molecular mechanisms involved in growth of largemouth bass (Micropterus salmoides) are poorly understood. The objective of this study was to screen growth-related genes of M. salmoides by RNA sequencing and identify growth-related single-nucleotide polymorphism (SNP) markers through a growth association study. The muscle transcriptomes of fast- and slow-growing largemouth bass were obtained using the RNA-Seq technique. A total of 54,058,178 and 54,742,444 qualified Illumina read pairs were obtained for the fast-growing and slow-growing groups, respectively, giving rise to 4,865,236,020 and 4,926,819,960 total clean bases, respectively. Gene expression profiling showed that 3,530 unigenes were differentially expressed between the fast-growing and slow-growing phenotypes (false discovery rate ≤0.001, the absolute value of log₂ (fold change) ≥1), including 1,441 up-regulated and 2,889 down-regulated unigenes in the fast-growing largemouth bass. Analysis of these genes revealed that several signalling pathways, including the growth hormone-insulin-like growth factor 1 axis and signalling pathway, the glycolysis pathway, and the myostatin/transforming growth factor beta signalling pathway, as well as heat shock protein, cytoskeleton, and myofibril component genes might be associated with muscle growth. From these genes, 10 genes with putative SNPs were selected, and 17 SNPs were genotyped successfully. Marker-trait analysis in 340 individuals of Youlu No. 1 largemouth bass revealed three SNPs associated with growth in key genes (phosphoenolpyruvate carboxykinase 1, FOXO3b, and heat shock protein beta-1). This research provides information about key genes and SNPs related to growth, providing new clues to understanding the molecular basis of largemouth bass growth.

Characterization of a novel disease resistance gene rtp3 and its association with VNN disease resistance in Asian seabass

Asian seabass, an important food fish in Southeast Asia, has suffered from nervous necrosis virus (NNV) infection, resulting in massive mortality of Asian seabass larvae and enormous economic losses. Identification and characterization of disease resistance genes is important. Previous transcriptome analysis of Asians seabass epithelial cells after NNV infection revealed a highly inducible gene, receptor-transporting protein 3 (rtp3), indicating it could play an important role in Asian seabass - NNV interaction. To characterize this gene, we determined its expression pattern and subcellular localization. The rtp3 was highly induced in most examined tissues and organs of Asian seabass after NNV infection, and protein Rtp3 was localized in cytoplasm. Further association study in multiple families revealed that a microsatellite marker, (GT)ntt(GT)n, in the 3' UTR of rtp3 was significantly associated with VNN disease resistance in Asian seabass. Our results imply that rtp3 may be a novel disease resistance gene in Asian seabass. This data could improve our understanding of molecular interaction between Asian seabass and NNV, and has the potential to be applied in marker-assisted selection for disease resistance breeding in Asian seabass.

Construction of high-resolution recombination maps in Asian seabass

Background

A high-density genetic map is essential for de novo genome assembly, fine mapping QTL for important complex traits, comparative genomic studies and understanding the mechanisms of genome evolution. Although a number of genomic resources are available in Asian seabass (Lates calcarifer), a high-density linkage map is still lacking. To facilitate QTL mapping for marker-assisted selection and genome assembly, and to understand the genome-wide recombination rates, we constructed high density linkage maps using three families and genotyping by sequencing.

Results

A high-density consensus linkage map consisting of 8, 274 markers was constructed based on sex-averaged genetic maps. The genetic maps were then aligned and integrated with the current genome assembly of Asian seabass. More than 90% of the genome contig sequences were anchored onto the consensus genetic map. Evidence of assembly errors in the current genome assembly was identified. A fragment of up to 2.5 Mb belonging to LG14 was assembled into Chr15. The length of family-specific sex-averaged maps ranged from 1348.96 to 1624.65 cM. Female maps were slightly longer than male maps using common markers. Female-to-male ratios were highly variable both across chromosomes within each family and throughout three families for each chromosome. However, the distribution patterns of recombination along chromosomes were similar between sexes across the whole genome. The overall recombination rates were significantly correlated with genome-wide GC content and the correlations were revealed to be stronger in females than in males.

Conclusions

These high-density genetic maps provide not only essential tools for facilitating de novo genome assembly and comparative genomic studies in teleosts, but also critical resources for fine mapping QTL and genome-wide association mapping for economically important traits in Asian seabass.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3462-z) contains supplementary material, which is available to authorized users.

Fine mapping QTL for resistance to VNN disease using a high-density linkage map in Asian seabass

Asian seabass has suffered from viral nervous necrosis (VNN) disease. Our previous study has mapped quantitative trait loci (QTL) for resistance to VNN disease. To fine map these QTL and identify causative genes, we identified 6425 single nucleotide polymorphisms (SNPs) from 85 dead and 94 surviving individuals. Combined with 155 microsatellites, we constructed a genetic map consisting of 24 linkage groups (LGs) containing 3000 markers, with an average interval of 1.27 cM. We mapped one significant and three suggestive QTL with phenotypic variation explained (PVE) of 8.3 to 11.0%, two significant and two suggestive QTL with PVE of 7.8 to 10.9%, for resistance in three LGs and survival time in four LGs, respectively. Further analysis one QTL with the largest effect identified protocadherin alpha-C 2-like (Pcdhac2) as the possible candidate gene. Association study in 43 families with 1127 individuals revealed a 6 bp insertion-deletion was significantly associated with disease resistance. qRT-PCR showed the expression of Pcdhac2 was significantly induced in the brain, muscle and skin after nervous necrosis virus (NNV) infection. Our results could facilitate marker-assisted selection (MAS) for resistance to NNV in Asian seabass and set up the basis for functional analysis of the potential causative gene for resistance.

Morpho-histological characterisation of the alimentary canal of an important food fish, Asian seabass (Lates calcarifer)

Asian seabass (Lates calcarifer) is a food fish of increasing aquaculture importance. In order to improve our understanding on the digestive system and feeding of this species, morphological and histological features of the gut were studied. Morphologically, the Asian seabass gut is defined by a short and muscular esophagus, well-developed stomach and comparatively short intestine. Mucous secreting goblet cells reactive to PAS (Periodic Acid Schiff) and AB (Alcian Blue) stain were present throughout the esophagus. The stomach was sac-like and could be distinguished into the cardiac, fundic and pyloric regions. Gastric glands and mucus cells were predominately present in the cardiac and fundic regions. Five finger-like pyloric caeca were present between the stomach and intestine. The intestine was a short, tubular structure with no morphological differences between the various regions. Histologically, the intestinal regions were similar, the main difference being in the number of goblet cells that increased from anterior to posterior intestine, with 114 ± 9, 153 ± 7 and 317 ± 21 goblet cells in the anterior, mid and posterior regions, respectively. The intestinal epithelium stained positively for PAS, but the staining was stronger for acidic glycoproteins. The rectum was similar to intestine, except for increased goblet cell numbers (anterior rectum: 529 ± 26; posterior rectum: 745 ± 29). Gut morpho-histology did not respond to salinity changes, however, there was a significant reduction of mucosal height, goblet cell numbers and muscularis thickness upon food deprivation.

Transcriptome analysis of genes responding to NNV infection in Asian seabass epithelial cells

Asian seabass is an important food fish in Southeast Asia. Viral nervous necrosis (VNN) disease, triggered by nervous necrosis virus (NNV) infection, has caused mass mortality of Asian seabass larvae, resulting in enormous economic losses in the Asian seabass industry. In order to better understand the complex molecular interaction between Asian seabass and NNV, we investigated the transcriptome profiles of Asian seabass epithelial cells, which play an essential role in immune regulation, after NNV infection. Using the next generation sequencing (NGS) technology, we sequenced mRNA from eight samples (6, 12, 24, 48 h post-inoculation) of mock and NNV-infected Asian seabass epithelial cell line, respectively. Clean reads were de novo assembled into a transcriptome consisting of 89026 transcripts with a N50 of 2617 bp. Furthermore, 251 differentially expressed genes (DEGs) in response to NNV infection were identified. Top DEGs include protein asteroid homolog 1-like (ASTE1), receptor-transporting protein 3 (RTP3), heat shock proteins 30 (HSP30) and 70 (HSP70), Viperin, interferon regulatory factor 3 (IRF3) and other genes related to innate immunity. Our data suggest that abundant and diverse genes corresponding to NNV infection. The results of this study could also offer vital information not only for identification of novel genes involved in Asian seabass-NNV interaction, but also for our understanding of the molecular mechanism of Asian seabass' response to viral infection. In addition, 24807 simple sequence repeats (SSRs) were detected in the assembled transcriptome, providing valuable resources for studying genetic variations and accelerating quantitative trait loci (QTL) mapping for disease resistance in Asian seabass in the future.

QTL variations for growth-related traits in eight distinct families of common carp (Cyprinus carpio)

BACKGROUND

Comparing QTL analyses of multiple pair-mating families can provide a better understanding of important allelic variations and distributions. However, most QTL mapping studies in common carp have been based on analyses of individual families. In order to improve our understanding of heredity and variation of QTLs in different families and identify important QTLs, we performed QTL analysis of growth-related traits in multiple segregating families.

RESULTS

We completed a genome scan for QTLs that affect body weight (BW), total length (TL), and body thickness (BT) of 522 individuals from eight full-sib families using 250 microsatellites evenly distributed across 50 chromosomes. Sib-pair and half-sib model mapping identified 165 QTLs on 30 linkage groups. Among them, 10 (genome-wide P <0.01 or P < 0.05) and 28 (chromosome-wide P < 0.01) QTLs exhibited significant evidence of linkage, while the remaining 127 exhibited a suggestive effect on the above three traits at a chromosome-wide (P < 0.05) level. Multiple QTLs obtained from different families affect BW, TL, and BT and locate at close or identical positions. It suggests that same genetic factors may control variability in these traits. Furthermore, the results of the comparative QTL analysis of multiple families showed that one QTL was common in four of the eight families, nine QTLs were detected in three of the eight families, and 26 QTLs were found common to two of the eight families. These common QTLs are valuable candidates in marker-assisted selection.

CONCLUSION

A large number of QTLs were detected in the common carp genome and associated with growth-related traits. Some of the QTLs of different growth-related traits were identified at similar chromosomal regions, suggesting a role for pleiotropy and/or tight linkage and demonstrating a common genetic basis of growth trait variations. The results have set up an example for comparing QTLs in common carp and provided insights into variations in the identified QTLs affecting body growth. Discovery of these common QTLs between families and growth-related traits represents an important step towards understanding of quantitative genetic variation in common carp.

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.

Identification of Quantitative Trait Loci Associated with the Skeletal Deformity LSK complex in Gilthead Seabream (Sparus aurata L.)

Morphological abnormalities, especially skeletal deformities, are some of the most important problems affecting gilthead seabream (Sparus aurata L.) aquaculture industry. In this study, a QTL analysis for LSK complex deformity in gilthead seabream is reported. LSK complex is a severe deformity consisting of a consecutive repetition of three vertebral deformities: lordosis, scoliosis, and kyphosis. Seventy-eight offspring from six breeders from a mass-spawning were analyzed: five full-sibling families, three maternal, and two paternal half-sibling families. They had shown a significant association with the LSK complex prevalence in a previous segregation analysis. Fish were genotyped using a set of multiplex PCRs (ReMsa1-13), which includes 106 microsatellite markers. Two methods were used to perform the QTL analysis: a linear regression with the GridQTL software and a linear mixed model with the Qxpak software. A total of 18 QTL were identified. Four of them (QTLSK3, 6, 12, and 14), located in LG5, 8, 17, and 20, respectively, were the most solid ones. These QTL were significant at genome level and showed an extremely large effect (>35%) with both methods. Markers close to the identified QTL showed a strong association with phenotype. Two of these molecular markers (DId-03-T and Bt-14-F) were considered as potential linked-to-this-deformity markers. The detection of these QTL supposes a critical step in the implementation of marker-assisted selection in this species, which could decrease the incidence of this deformity and other related deformities. The identification of these QTL also represents a major step towards the study of the etiology of skeletal deformities in this species.

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 (QTLs) 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 QTLs (body weight: BW and total length: TL) were detected using non-parametric mapping, Kruskal-Wallis 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 17 F 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 marker assisted selection (MAS) breeding.

Construction of a high-density linkage map and fine mapping of QTL for growth in Asian seabass

A high-density genetic map is essential for comparative genomic studies and fine mapping of QTL, and can also facilitate genome sequence assembly. Here, a high density genetic map of Asian seabass was constructed with 3321 SNPs generated by sequencing 144 individuals in a F₂ family. The length of the map was 1577.67 cM with an average marker interval of 0.52 cM. A high level of genomic synteny among Asian seabass, European seabass, Nile tilapia and stickleback was detected. Using this map, one genome-wide significant and five suggestive QTL for growth traits were detected in six linkage groups (i.e. LG4, LG5, LG11, LG13, LG14 and LG15). These QTL explained 10.5–16.0% of phenotypic variance. A candidate gene, ACOX1 within the significant QTL on LG5 was identified. The gene was differentially expressed between fast- and slow-growing Asian seabass. The high-density SNP-based map provides an important tool for fine mapping QTL in molecular breeding and comparative genome analysis.

Genetic and genomic analyses for economically important traits and their applications in molecular breeding of cultured fish

The traits of cultured fish must continually be genetically improved to supply high-quality animal protein for human consumption. Economically important fish traits are controlled by multiple gene quantitative trait loci (QTL), most of which have minor effects, but a few genes may have major effects useful for molecular breeding. In this review, we chose relevant studies on some of the most intensively cultured fish and concisely summarize progress on identifying and verifying QTLs for such traits as growth, disease and stress resistance and sex in recent decades. The potential applications of these major-effect genes and their associated markers in marker-assisted selection and molecular breeding, as well as future research directions are also discussed. These genetic and genomic analyses will be valuable for elucidating the mechanisms modulating economically important traits and to establish more effective molecular breeding techniques in fish.

Genetic mapping and QTL analysis of growth traits in the large yellow croaker Larimichthys crocea

Large yellow croaker (Larimichthys crocea) is an important maricultured species in China. A genetic linkage map of the large yellow croaker was constructed using type II microsatellites and expressed sequence tag (EST)-derived microsatellites in two half-sib families (two females and one male). A total of 289 microsatellite markers (contained 93 EST-SSRs) were integrated into 24 linkage groups, which agreed with the haploid chromosome number. The map spanned a length of 1,430.8 cm with an average interval of 5.4 cm, covering 83.9 % of the estimated genome size (1,704.8 cm). A total of seven quantitative trait locis (QTLs) were detected for growth traits on five linkage groups, including two 1 % and five 5 % chromosome-wide significant QTLs, and explained from 2.33 to 5.31 % of the trait variation. The identified QTLs can be applied in marker-assisted selection programs to improve the growth traits.

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.

A genome scan for quantitative trait loci associated with Vibrio anguillarum infection resistance in Japanese flounder (Paralichthys olivaceus) by bulked segregant analysis

A recent genetic linkage map was employed to detect quantitative trait loci (QTLs) associated with Vibrio anguillarum resistance in Japanese flounder. An F1 family established and challenged with V. anguillarum in 2009 was used for QTL mapping. Of the 221 simple sequence repeat (SSR) markers used to detect polymorphisms in the parents of F1, 170 were confirmed to be polymorphic. The average distance between the markers was 10.6 cM. Equal amounts of genomic DNA from 15 fry that died early and from 15 survivors were pooled separately to constitute susceptible bulk and resistance bulk DNA. Bulked segregant analysis and QTL mapping were combined to detect candidate SSR markers and regions associated with the disease. A genome scan identified four polymorphic SSR markers, two of which were significantly different between susceptible and resistance bulk (P=0.008). These two markers were located in linkage group (LG) 7; therefore, all the SSR markers in LG7 were genotyped in all the challenged fry by single marker analysis. Using two different models, 11-17 SSR markers were detected with different levels of significance. To confirm the associations of these markers with the disease, composite interval mapping was employed to genotype all the challenged individuals. One and three QTLs, which explained more than 60 % of the phenotypic variance, were detected by the two models. Two of the QTLs were located at 48.6 cM. The common QTL may therefore be a major candidate region for disease resistance against V. anguillarum infection.

QTL mapping for economically important traits of common carp (Cyprinus carpio L.)

Quantitative trait loci (QTL) were analyzed for three economically important traits, i.e., body weight (BW), body length (BL), and body thickness (BT), in an F1 family of common carp holding the 190 progeny. A genetic linkage map spanning 3,301 cM in 50 linkage groups with 627 markers and an average distance of 5.6 cM was utilized for QTL mapping. Sixteen QTLs associated with all three growth-related traits were scattered across ten linkage groups, LG6, LG10, LG17, LG19, LG25, LG27, LG28, LG29, LG30, and LG39. Six QTLs for BW and five each for BL and BT explained phenotypic variance in the range 17.0-32.1%. All the nearest markers of QTLs were found to be significantly (p ≤ 0.05) related with the trait. Among these QTLs, a total of four, two (qBW30 and qBW39) related with BW, one (qBL39) associated with BL, and one (qBT29) related to BT, were found to be the major QTLs with a phenotypic variance of >20%. qBW30 and qBW39 with the nearest markers HLJ1691 and HLJ1843, respectively, show significant values of 0.0038 and 0.0031, correspondingly. QTLs qBL39 and qBT29 were found to have significant values of 0.0047 and 0.0015, respectively. Three QTLs (qBW27, qBW30, qBW39) of BW, two for BL (qBL19, qBL39), and two for BT (qBT6, qBT25) found in this study were similar to populations with different genetic backgrounds. In this study, the genomic region controlling economically important traits were located. These genomic regions will be the major sources for the discovery of important genes and pathways associated with growth-related traits in common carp.

Mapping quantitative trait loci (QTL) for body weight, length and condition factor traits in backcross (BC1) family of Common carp (Cyprinus carpio L.)

Body weight and length are economical important traits in aquaculture species influenced by quantitative trait loci (QTL) and environmental factors. In this study, a backcross (BC1) common carp family, with 86 progeny, was utilized to construct genetic map for preliminary QTL mapping. The genetic map was constructed with 366 markers, including 191 SNP from gene, coverage 50 linkage groups with an average marker distance of 18.5 cM. A total of fourteen QTLs associated with body weight (BW), body length (BL) and condition factor (K) were detected on ten linkage groups (LGs). Among these QTLs detected, three (qBW8, qBL8 and qK8) were associated with BW, BL and K respectively, were mapped on LG8. qBW8 and qK8 were identified on similar interval neared locus HLJ2394 explained 14.9 and 20.9 % of the phenotype variance, while qBL8 was identified on separate nearby locus HLJ571 with 30.8 % of phenotype variance. Two QTLs, qBW13 and qK13, related with BW and K separately, were found on LG13 at different locus with phenotype variance of 25.3 and 20.9 %. Other two QTLs, qBW19 and qBL19, associated to BW and BL were mapped on same region near SNP0626 on LG19, and explained 10.3 and 15.6 % of phenotype variance. While other seven QTLs related with BW and BL were located on different LGs. Confidential interval was ranged from 1.1 to 10 cM in the present study. These markers, with lower QTL interval, have great influence on the body weight and length. Therefore, these QTLs will be helpful to find out the genes related with specific trait.

Development of molecular resources for an intertidal clam, Sinonovacula constricta, using 454 transcriptome sequencing

Background

The razor clam Sinonovacula constricta is a benthic intertidal bivalve species with important commercial value. Despite its economic importance, knowledge of its transcriptome is scarce. Next generation sequencing technologies offer rapid and efficient tools for generating large numbers of sequences, which can be used to characterize the transcriptome, to develop effective molecular markers and to identify genes associated with growth, a key breeding trait.

Results

Total RNA was isolated from the mantle, gill, liver, siphon, gonad and muscular foot tissues. High-throughput deep sequencing of S. constricta using 454 pyrosequencing technology yielded 859,313 high-quality reads with an average read length of 489 bp. Clustering and assembly of these reads produced 16,323 contigs and 131,346 singletons with average lengths of 1,376 bp and 458 bp, respectively. Based on transcriptome sequencing, 14,615 sequences had significant matches with known genes encoding 147,669 predicted proteins. Subsequently, previously unknown growth-related genes were identified. A total of 13,563 microsatellites (SSRs) and 13,634 high-confidence single nucleotide polymorphism loci (SNPs) were discovered, of which almost half were validated.

Conclusion

De novo sequencing of the razor clam S. constricta transcriptome on the 454 GS FLX platform generated a large number of ESTs. Candidate growth factors and a large number of SSRs and SNPs were identified. These results will impact genetic studies of S. constricta.

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.

Primary genome scan for complex body shape-related traits in the common carp Cyprinus carpio

To identify quantitative trait loci (QTL) that affect body shape in common carp Cyprinus carpio, a linkage map, 2159·23 cM long, was constructed with a total of 307 markers covering 51 linkage groups (LG). The map included 167 new single nucleotide polymorphism (SNP) markers derived from expressed sequence tags (EST) together with 140 microsatellite markers reported earlier. A primary genome scan was conducted for QTL for standard length (L(S)), head length (L(H)), body height (H(B)), body width (W(B)) and tail length (L(TAIL)) in an F1 line containing 92 offspring. A total of 15 suggestive QTL on six LGs were found to associate with L(S), L(H), H(B), W(B) and L(TAIL) which explained 10·7-17·4% of the variance. Five significant QTL were detected for body-shape related traits and located for LGs (lg1, 12 and 20). These QTL included: one associated with L(S) (21·1% variance explained), three for H(B) (almost 20% variance explained) and one for W(B) (20·7% variance explained).

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.

Molecular parentage analysis is essential in breeding Asian seabass

In aquaculture species, maintaining pedigree information and genetic variation in each generation is essential, but very difficult. In this study, we used nine microsatellites to genotype 2,520 offspring from four independent full-factorial crosses (10 males ×10 females) of Asian seabass to reconstruct pedigree and monitor the change of genetic variations. In all four crosses, over 96.8% of the offspring could be assigned to their parents, indicating the high power of the nine microsatellites for parentage assignment. This study revealed several interesting results: (1). In all four crosses, the contribution of parents to offspring was significantly uneven, and some dominant breeding fishes (i.e. brooders) were found; (2). In two mass crosses where the brooders were carefully checked for reproductive status, a majority (≥90%) of brooders contributed to offspring, whereas in another two crosses, where the brooders were randomly picked without checking reproductive status, only a few brooders (40.0–45.0%) produced offspring; (3). Females had more problems in successful spawning compared to males; and (4). In the two crosses where a few brooders produced offspring, there was a substantial loss in allelic (24.1–34.3%) and gene (20.5–25.7%) diversities in offspring, while in the other two crosses, the majority of allelic (96.8–97.0%) and gene diversities (94.8–97.1%) were maintained. These observations suggest that a routine molecular parentage analysis is required to maintain both allelic and gene diversity in breeding Asian seabass.

RNA-Seq identifies SNP markers for growth traits in rainbow trout

Fast growth is an important and highly desired trait, which affects the profitability of food animal production, with feed costs accounting for the largest proportion of production costs. Traditional phenotype-based selection is typically used to select for growth traits; however, genetic improvement is slow over generations. Single nucleotide polymorphisms (SNPs) explain 90% of the genetic differences between individuals; therefore, they are most suitable for genetic evaluation and strategies that employ molecular genetics for selective breeding. SNPs found within or near a coding sequence are of particular interest because they are more likely to alter the biological function of a protein. We aimed to use SNPs to identify markers and genes associated with genetic variation in growth. RNA-Seq whole-transcriptome analysis of pooled cDNA samples from a population of rainbow trout selected for improved growth versus unselected genetic cohorts (10 fish from 1 full-sib family each) identified SNP markers associated with growth-rate. The allelic imbalances (the ratio between the allele frequencies of the fast growing sample and that of the slow growing sample) were considered at scores >5.0 as an amplification and <0.2 as loss of heterozygosity. A subset of SNPs (n = 54) were validated and evaluated for association with growth traits in 778 individuals of a three-generation parent/offspring panel representing 40 families. Twenty-two SNP markers and one mitochondrial haplotype were significantly associated with growth traits. Polymorphism of 48 of the markers was confirmed in other commercially important aquaculture stocks. Many markers were clustered into genes of metabolic energy production pathways and are suitable candidates for genetic selection. The study demonstrates that RNA-Seq at low sequence coverage of divergent populations is a fast and effective means of identifying SNPs, with allelic imbalances between phenotypes. This technique is suitable for marker development in non-model species lacking complete and well-annotated genome reference sequences.

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.

Detection of growth-related QTL in turbot (Scophthalmus maximus)

Background

The turbot (Scophthalmus maximus) is a highly appreciated European aquaculture species. Growth related traits constitute the main goal of the ongoing genetic breeding programs of this species. The recent construction of a consensus linkage map in this species has allowed the selection of a panel of 100 homogeneously distributed markers covering the 26 linkage groups (LG) suitable for QTL search. In this study we addressed the detection of QTL with effect on body weight, length and Fulton's condition factor.

Results

Eight families from two genetic breeding programs comprising 814 individuals were used to search for growth related QTL using the panel of microsatellites available for QTL screening. Two different approaches, maximum likelihood and regression interval mapping, were used in order to search for QTL. Up to eleven significant QTL were detected with both methods in at least one family: four for weight on LGs 5, 14, 15 and 16; five for length on LGs 5, 6, 12, 14 and 15; and two for Fulton's condition factor on LGs 3 and 16. In these LGs an association analysis was performed to ascertain the microsatellite marker with the highest apparent effect on the trait, in order to test the possibility of using them for marker assisted selection.

Conclusions

The use of regression interval mapping and maximum likelihood methods for QTL detection provided consistent results in many cases, although the high variation observed for traits mean among families made it difficult to evaluate QTL effects. Finer mapping of detected QTL, looking for tightly linked markers to the causative mutation, and comparative genomics are suggested to deepen in the analysis of QTL in turbot so they can be applied in marker assisted selection programs.

Genetic linkage mapping and analysis of muscle fiber-related QTLs in common carp (Cyprinus carpio L.)

A genetic linkage map of common carp (Cyprinus carpio L.) was constructed using Type I and Type II microsatellite markers and a pseudo-testcross mapping strategy. The microsatellite markers were isolated from microsatellite-enriched genomic libraries and tested for their segregation in a full-sib mapping panel containing 92 individuals. A total of 161 microsatellite loci were mapped into 54 linkage groups. The total lengths of the female, male and consensus maps were 2,000, 946, and 1,852 cM, with an average marker spacing of approximately 13, 7, and 11 cM, respectively. Muscle fiber-related traits, including muscle fiber cross-section area and muscle fiber density, were mapped to the genetic map. Three QTLs for muscle fiber cross-section area and two QTLs for muscle fiber density were identified when considering both significant and suggestive QTL effects. The QTLs with largest effects for muscle fiber cross-section area and muscle fiber density were 21.9% and 18.9%, and they were located in LG3, respectively.

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

Background

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

Results

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

Conclusions

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

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

The caudal fin represents a fundamental design feature of fishes and plays an important role in locomotor dynamics in fishes. The shape of caudal is an important parameter in traditional systematics. However, little is known about genes involved in the development of different forms of caudal fins. This study was conducted to identify and map quantitative trait loci (QTL) affecting the length of caudal fin and the ratio between tail length and standard body length in Asian seabass (Lates calcarifer). One F1 family containing 380 offspring was generated by crossing two unrelated individuals. One hundred and seventeen microsatellites almost evenly distributed along the whole genome were genotyped. Length of caudal fin at 90 days post-hatch was measured. QTL analysis detected six significant (genome-wide significant) and two suggestive (linkage-group-wide significant) QTL on seven linkage groups. The six significant QTL explained 5.5-16.6% of the phenotypic variance, suggesting these traits were controlled by multiple genes. Comparative genomics analysis identified several potential candidate genes for the length of caudal fin. The QTL for the length of caudal fin detected for the first time in marine fish may provide a starting point for the future identification of genes involved in the development of different forms of caudal fins in fishes.

Quantitative trait loci involved in sex determination and body growth in the gilthead sea bream (Sparus aurata L.) through targeted genome scan

Among vertebrates, teleost fish exhibit a considerably wide range of sex determination patterns that may be influenced by extrinsic parameters. However even for model fish species like the zebrafish Danio rerio the precise mechanisms involved in primary sex determination have not been studied extensively. The zebrafish, a gonochoristic species, is lacking discernible sex chromosomes and the sex of juvenile fish is difficult to determine. Sequential protandrous hermaphrodite species provide distinct determination of the gender and allow studying the sex determination process by looking at the mechanism of sex reversal. This is the first attempt to understand the genetic basis of phenotypic variation for sex determination and body weight in a sequential protandrous hermaphrodite species, the gilthead sea bream (Sparus aurata). This work demonstrates a fast and efficient strategy for Quantitative Trait Loci (QTL) detection in the gilthead sea bream, a non-model but target hermaphrodite fish species. Therefore a comparative mapping approach was performed to query syntenies against two other Perciformes, the European sea bass (Dicentrarchus labrax), a gonochoristic species and the Asian sea bass (Lates calcarifer) a protandrous hermaphrodite. In this manner two significant QTLs, one QTL affecting both body weight and sex and one QTL affecting sex, were detected on the same linkage group. The co-segregation of the two QTLs provides a genomic base to the observed genetic correlation between these two traits in sea bream as well as in other teleosts. The identification of QTLs linked to sex reversal and growth, will contribute significantly to a better understanding of the complex nature of sex determination in S. aurata where most individuals reverse to the female sex at the age of two years through development and maturation of the ovarian portion of the gonad and regression of the testicular area. [Genomic sequences reported in this manuscript have been submitted to GenBank under accession numbers HQ021443-HQ021749.].

A first generation BAC-based physical map of the Asian seabass (Lates calcarifer)

Background

The Asian seabass (Lates calcarifer) is an important marine foodfish species in Southeast Asia and Australia. Genetic improvement of this species has been achieved to some extent through selective breeding programs since 1990s. Several genomic tools such as DNA markers, a linkage map, cDNA and BAC libraries have been developed to assist selective breeding. A physical map is still lacking, although it is essential for positional cloning of genes located in quantitative trait loci (QTL) and assembly of whole genome sequences.

Methodology/Principal Findings

A genome-wide physical map of the Asian seabass was constructed by restriction fingerprinting of 38,208 BAC clones with SNaPshot HICF FPC technique. A total of 30,454 were assembled into 2,865 contigs. The physical length of the assembled contigs summed up to 665 Mb. Analyses of some contigs using different methods demonstrated the reliability of the assembly.

Conclusions/Significance

The present physical map is the first physical map for Asian seabass. This physical map will facilitate the fine mapping of QTL for economically important traits and the positional cloning of genes located in QTL. It will also be useful for the whole genome sequencing and assembly. Detailed information about BAC-contigs and BAC clones are available upon request.

Mapping of quantitative trait loci for flesh colour and growth traits in Atlantic salmon (Salmo salar)

BACKGROUND

Flesh colour and growth related traits in salmonids are both commercially important and of great interest from a physiological and evolutionary perspective. The aim of this study was to identify quantitative trait loci (QTL) affecting flesh colour and growth related traits in an F2 population derived from an isolated, landlocked wild population in Norway (Byglands Bleke) and a commercial production population.

METHODS

One hundred and twenty-eight informative microsatellite loci distributed across all 29 linkage groups in Atlantic salmon were genotyped in individuals from four F2 families that were selected from the ends of the flesh colour distribution. Genotyping of 23 additional loci and two additional families was performed on a number of linkage groups harbouring putative QTL. QTL analysis was performed using a line-cross model assuming fixation of alternate QTL alleles and a half-sib model with no assumptions about the number and frequency of QTL alleles in the founder populations.

RESULTS

A moderate to strong phenotypic correlation was found between colour, length and weight traits. In total, 13 genome-wide significant QTL were detected for all traits using the line-cross model, including three genome-wide significant QTL for flesh colour (Chr 6, Chr 26 and Chr 4). In addition, 32 suggestive QTL were detected (chromosome-wide P < 0.05). Using the half-sib model, six genome-wide significant QTL were detected for all traits, including two for flesh colour (Chr 26 and Chr 4) and 41 suggestive QTL were detected (chromosome-wide P < 0.05). Based on the half-sib analysis, these two genome-wide significant QTL for flesh colour explained 24% of the phenotypic variance for this trait.

CONCLUSIONS

A large number of significant and suggestive QTL for flesh colour and growth traits were found in an F2 population of Atlantic salmon. Chr 26 and Chr 4 presented the strongest evidence for significant QTL affecting flesh colour, while Chr 10, Chr 5, and Chr 4 presented the strongest evidence for significant QTL affecting growth traits (length and weight). These QTL could be strong candidates for use in marker-assisted selection and provide a starting point for further characterisation of the genetic components underlying flesh colour and growth.

Malaria parasite genome scan: insights into antimalarial resistance

Genome scan and genotype-phenotype association study offer excellent opportunities to unearth drug/vaccine targets in human pathogens including malaria parasites. A recently conducted such study in worldwide isolates in the most devastating malaria parasite Plasmodium falciparum has reported important genomic information on genetic basis of antimalarial resistance. Several unknown genes were also found to be under strong influence of natural selection. The findings provide important insights into the malaria parasite genome evolution in general and to use this information to develop more focused malaria control strategies, in particular.

Construction of a BAC library and mapping BAC clones to the linkage map of Barramundi, Lates calcarifer

BACKGROUND

Barramundi (Lates calcarifer) is an important farmed marine food fish species. Its first generation linkage map has been applied to map QTL for growth traits. To identify genes located in QTL responsible for specific traits, genomic large insert libraries are of crucial importance. We reported herein a bacterial artificial chromosome (BAC) library and the mapping of BAC clones to the linkage map.

RESULTS

This BAC library consisted of 49,152 clones with an average insert size of 98 kb, representing 6.9-fold haploid genome coverage. Screening the library with 24 microsatellites and 15 ESTs/genes demonstrated that the library had good genome coverage. In addition, 62 novel microsatellites each isolated from 62 BAC clones were mapped onto the first generation linkage map. A total of 86 BAC clones were anchored on the linkage map with at least one BAC clone on each linkage group.

CONCLUSION

We have constructed the first BAC library for L. calcarifer and mapped 86 BAC clones to the first generation linkage map. This BAC library and the improved linkage map with 302 DNA markers not only supply an indispensable tool to the integration of physical and linkage maps, the fine mapping of QTL and map based cloning genes located in QTL of commercial importance, but also contribute to comparative genomic studies and eventually whole genome sequencing.