Search

Background Genetic diversity within a species reflects population evolution, ecology, and ability to adapt. Genome-wide population surveys of both natural and introduced populations provide insights into genetic diversity, the evolutionary processes and the genetic basis underlying local adaptation. Grass carp is the most important freshwater foodfish species for food and water weed control. However, there is as yet no overall picture on genetic variations and population structure of this species, which is important for its aquaculture. Results We used 43,310 SNPs to infer the population structure, evidence of local adaptation and sources of introduction. The overall genetic differentiation of this species was low. The native populations were differentiated into three genetic clusters, corresponding to the Yangtze, Pearl and Heilongjiang River Systems, respectively. The populations in Malaysia, India and Nepal were introduced from both the Yangtze and Pearl River Systems. Loci and genes involved in putative local selection for native locations were identified. Evidence of both positive and balancing selection was found in the introduced locations. Genes associated with loci under putative selection were involved in many biological functions. Outlier loci were grouped into clusters as genomic islands within some specific genomic regions, which likely agrees with the divergence hitchhiking scenario of divergence-with-gene-flow. Conclusions This study, for the first time, sheds novel insights on the population differentiation of the grass carp, genetics of its strong ability in adaption to diverse environments and sources of some introduced grass carp populations. Our data also suggests that the natural populations of the grass carp have been affected by the aquaculture besides neutral and adaptive forces. Electronic supplementary material The online version of this article (10.1186/s12864-019-5872-1) contains supplementary material, which is available to authorized users.

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 F2 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.

Background Tilapia is the common name for a group of cichlid fishes and is one of the most important aquacultured freshwater food fish. Mozambique tilapia and its hybrids, including red tilapia are main representatives of salt tolerant tilapias. A linkage map is an essential framework for mapping QTL for important traits, positional cloning of genes and understanding of genome evolution. Results We constructed a consensus linkage map of Mozambique tilapia and red tilapia using 95 individuals from two F1 families and 401 microsatellites including 282 EST-derived markers. In addition, we conducted comparative mapping and searched for sex-determining loci on the whole genome. These 401 microsatellites were assigned to 22 linkage groups. The map spanned 1067.6 cM with an average inter-marker distance of 3.3 cM. Comparative mapping between tilapia and stickleback, medaka, pufferfish and zebrafish revealed clear homologous relationships between chromosomes from different species. We found evidence for the fusion of two sets of two independent chromosomes forming two new chromosome pairs, leading to a reduction of 24 chromosome pairs in their ancestor to 22 pairs in tilapias. The XY sex determination locus in Mozambique tilapia was mapped on LG1, and verified in five families containing 549 individuals. The major XY sex determination locus in red tilapia was located on LG22, and verified in two families containing 275 individuals. Conclusions A first-generation linkage map of salt tolerant tilapia was constructed using 401 microsatellites. Two separate fusions of two sets of two independent chromosomes may lead to a reduction of 24 chromosome pairs in their ancestor to 22 pairs in tilapias. The XY sex-determining loci from Mozambique tilapia and red tilapia were mapped on LG1 and LG22, respectively. This map provides a useful resource for QTL mapping for important traits and comparative genome studies. The DNA markers linked to the sex-determining loci could be used in the selection of YY males for breeding all-male populations of salt tolerant tilapia, as well as in studies on mechanisms of sex determination in fish.

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.

Background Grass carp (Ctenopharyngodon idella) belongs to the family Cyprinidae which includes more than 2000 fish species. It is one of the most important freshwater food fish species in world aquaculture. A linkage map is an essential framework for mapping traits of interest and is often the first step towards understanding genome evolution. The aim of this study is to construct a first generation genetic map of grass carp using microsatellites and SNPs to generate a new resource for mapping QTL for economically important traits and to conduct a comparative mapping analysis to shed new insights into the evolution of fish genomes. Results We constructed a first generation linkage map of grass carp with a mapping panel containing two F1 families including 192 progenies. Sixteen SNPs in genes and 263 microsatellite markers were mapped to twenty-four linkage groups (LGs). The number of LGs was corresponding to the haploid chromosome number of grass carp. The sex-specific map was 1149.4 and 888.8 cM long in females and males respectively whereas the sex-averaged map spanned 1176.1 cM. The average resolution of the map was 4.2 cM/locus. BLAST searches of sequences of mapped markers of grass carp against the whole genome sequence of zebrafish revealed substantial macrosynteny relationship and extensive colinearity of markers between grass carp and zebrafish. Conclusions The linkage map of grass carp presented here is the first linkage map of a food fish species based on co-dominant markers in the family Cyprinidae. This map provides a valuable resource for mapping phenotypic variations and serves as a reference to approach comparative genomics and understand the evolution of fish genomes and could be complementary to grass carp genome sequencing project.

Background Body weight and length are economically important traits in foodfish species influenced by quantitative trait loci (QTL) and environmental factors. It is usually difficult to dissect the genetic and environmental effects. Asian seabass (Lates calcarifer) is an important marine foodfish species with a compact genome (~700 Mb). The recent construction of a first generation linkage map of Asian seabass with 240 microsatellites provides a good opportunity to determine the number and position of QTL, and the magnitude of QTL effects with a genome scan. Results We conducted a genome scan for QTL affecting body weight, standard length and condition factors in an F1 family containing 380 full-sib individuals from a breeding stock by using 97 microsatellites evenly covering 24 chromosomes. Interval mapping and multiple QTL model mapping detected five significant and 27 suggestive QTL on ten linkage groups (LGs). Among the five significant QTL detected, three (qBW2-a, qTL2-a and qSL2-a) controlling body weight, total and standard length respectively, were mapped on the same region near Lca287 on LG2, and explained 28.8, 58.9 and 59.7% of the phenotypic variance. The other two QTL affecting body weight, qBW2-b and qBW3, were located on LG2 and 3, and accounted for 6.4 and 8.8% of the phenotypic variance. Suggestive QTL associated with condition factors are located on six different LGs. Conclusion This study presents the first example of QTL detection for growth-related traits in an F1 family of a marine foodfish species. The results presented here will enable further fine-mapping of these QTL for marker-assisted selection of the Asian seabass, eventually identifying individual genes responsible for growth-related traits.

Background Mitochondrial DNA-derived sequences have become popular markers for evolutionary studies, as their comparison may yield significant insights into the evolution of both the organisms and their genomes. From the more than 24,000 teleost species, only 254 complete mtDNA sequences are available (GenBank status on 06 Sep 2006). In this paper, we report the complete mitochondrial genome sequence of Asian arowana, a basal bonytongue fish species, which belongs to the order of Osteoglossiformes. Results The complete mitochondrial genomic sequence (mtDNA) of Asian arowana (Scleropages formosus) was determined by using shotgun sequencing method. The length of Asian arowana mtDNA is ca. 16,650 bp (its variation is due to polymorphic repeats in the control region), containing 13 protein-coding genes, 22 tRNA and 2 rRNA genes. Twelve of the thirteen protein coding genes were found to be encoded by the heavy strand in the order typically observed for vertebrate mitochondrial genomes, whereas only nad6 was located on the light strand. An interesting feature of Asian arowana mitogenome is that two different repeat arrays were identified in the control region: a 37 bp tandem repeat at the 5' end and an AT-type dinucleotide microsatellite at the 3' end. Both repeats show polymorphism among the six individuals tested; moreover the former one is present in the mitochondrial genomes of several other teleost groups. The TACAT motif described earlier only from mammals and lungfish was found in the tandem repeat of several osteoglossid and eel species. Phylogenetic analysis of fish species representing Actinopterygii and Sarcopterygii taxa has shown that the Asian arowana is located near the baseline of the teleost tree, confirming its status among the ancestral teleost lineages. Conclusion The mitogenome of Asian arowana is very similar to the typical vertebrate mitochondrial genome in terms of gene arrangements, codon usage and base composition. However its control region contains two different types of repeat units at both ends, an interesting feature that to our knowledge has never been reported before for other vertebrate mitochondrial control regions. Phylogenetic analysis using the complete mtDNA sequence of Asian arowana confirmed that it belongs to an ancestral teleost lineage.

Background Fish diseases caused by pathogens are limiting their production and trade, affecting the economy generated by aquaculture. Innate immunity system is the first line of host defense in opposing pathogenic organisms or any other foreign material. For identification of immune-related genes in Asian seabass Lates calcarifer, an important marine foodfish species, we injected bacterial lipopolysaccharide (LPS), a commonly used elicitor of innate immune responses to eight individuals at the age of 35 days post-hatch and applied the suppression subtractive hybridization (SSH) technique to selectively amplify spleen cDNA of differentially expressed genes. Results Sequencing and bioinformatic analysis of 3351 ESTs from two SSH libraries yielded 1692 unique transcripts. Of which, 618 transcripts were unknown/novel genes and the remaining 1074 were similar to 743 known genes and 105 unannotated mRNA sequences available in public databases. A total of 161 transcripts were classified to the category "response to stimulus" and 115 to "immune system process". We identified 25 significantly up-regulated genes (including 2 unknown transcripts) and 4 down-regulated genes associated with immune-related processes upon challenge with LPS. Quantitative real-time PCR confirmed the differential expression of these genes after LPS challenge. Conclusions The present study identified 1692 unique transcripts upon LPS challenge for the first time in Asian seabass by using SSH, sequencing and bioinformatic analysis. Some of the identified transcripts are vertebrate homologues and others are hitherto unreported putative defence proteins. The obtained immune-related genes may allow for a better understanding of immunity in Asian seabass, carrying out detailed functional analysis of these genes and developing strategies for efficient immune protection against infections in Asian seabass.

Background Microsatellites in cDNA are useful as molecular markers because they represent transcribed genes and can be used as anchor markers for linkage and comparative mapping, as well as for studying genome evolution. Microsatellites in cDNA can be detected in existing ESTs by data mining. However, in most fish species, no ESTs are available or the number of ESTs is limited, although fishes represent half of the vertebrates on the earth. We developed a simple and efficient method for isolation of microsatellites from cDNA in fish. Results The method included normalization of 150 ng cDNA using 0.5 U duplex-specific nuclease (DSN) at 65°C for 30 min, enrichment of microsatellites using biotinylated oligonucleotides and magnetic field, and directional cloning of cDNA into a vector. We tested this method to enrich CA- and GA-microsatellites from cDNA of Asian seabass, and demonstrated that enrichment of microsatellites from normalized cDNA could increased the efficiency of microsatellite isolation over 30 times as compared to direct sequencing of clones from cDNA libraries. One hundred and thirty-nine (36.2%) out of 384 clones from normalized cDNA contained microsatellites. Unique microsatellite sequences accounted for 23.6% (91/384) of sequenced clones. Sixty microsatellites isolated from cDNA were characterized, and 41 were polymorphic. The average allele number of the 41 microsatellites was 4.85 ± 0.54, while the expected heterozygosity was 0.56 ± 0.03. All the isolated microsatellites inherited in a Mendelian pattern. Conclusion Normalization of cDNA substantially increased the efficiency of enrichment of microsatellites from cDNA. The described method for isolation of microsatellites from cDNA has the potential to be applied to a wide range of fish species. The microsatellites isolated from cDNA could be useful for linkage and comparative mapping, as well as for studying genome evolution.

Background Starvation not only affects the nutritional and health status of the animals, but also the microbial composition in the host’s intestine. Next-generation sequencing provides a unique opportunity to explore gut microbial communities and their interactions with hosts. However, studies on gut microbiomes have been conducted predominantly in humans and land animals. Not much is known on gut microbiomes of aquatic animals and their changes under changing environmental conditions. To address this shortcoming, we determined the microbial gene catalogue, and investigated changes in the microbial composition and host-microbe interactions in the intestine of Asian seabass in response to starvation. Results We found 33 phyla, 66 classes, 130 orders and 278 families in the intestinal microbiome. Proteobacteria (48.8%), Firmicutes (15.3%) and Bacteroidetes (8.2%) were the three most abundant bacteria taxa. Comparative analyses of the microbiome revealed shifts in bacteria communities, with dramatic enrichment of Bacteroidetes, but significant depletion of Betaproteobacteria in starved intestines. In addition, significant differences in clusters of orthologous groups (COG) functional categories and orthologous groups were observed. Genes related to antibiotic activity in the microbiome were significantly enriched in response to starvation, and host genes related to the immune response were generally up-regulated. Conclusions This study provides the first insights into the fish intestinal microbiome and its changes under starvation. Further detailed study on interactions between intestinal microbiomes and hosts under dynamic conditions will shed new light on how the hosts and microbes respond to the changing environment.