Construction of a genetic linkage map for silver carp (Hypophthalmichthys molitrix)

Construction of a high-density genetic linkage map and mapping of quantitative trait loci for growth-related traits in silver carp (Hypophthalmichthys molitrix)

High-density genetic map and quantitative trait loci (QTL) mapping are powerful tools for identifying genomic regions that may be responsible for such polygenic trait as growth. A high-density genetic linkage map was constructed by sequencing 198 individuals in a F₁ family of silver carp (Hypophthalmichthys molitrix) in this study. This genetic map spans a length of 2,721.07 cM with 3,134 SNPs distributed on 24 linkage groups (LGs). Comparative genomic mapping presented a high level of syntenic relationship between silver carp and zebrafish. We detected one major and nineteen suggestive QTL for 4 growth-related traits (body length, body height, head length and body weight) at 6, 12 and 18 months post hatch (mph), explaining 10.2~19.5% of phenotypic variation. All six QTL for growth traits of 12 mph generally overlapped with QTL for 6 mph, while the majority of QTL for 18 mph were identified on two additional LGs, which may reveal a different genetic modulation during early and late muscle growth stages. Four potential candidate genes were identified from the QTL regions by homology searching of marker sequences against zebrafish genome. Hepcidin, a potential candidate gene identified from a QTL interval on LG16, was significantly associated with growth traits in the analyses of both phenotype-SNP association and mRNA expression between small-size and large-size groups of silver carp. These results provide a basis for elucidating the genetic mechanisms for growth and body formation in silver carp, a world aquaculture fish.

Comparative analysis of spleen transcriptome detects differences in evolutionary adaptation of immune defense functions in bighead carp and silver carp

The evolutionary divergence of the immune defense functions in bighead carp (A. nobilis) and silver carp (H. molitrix) is still not understood at the molecular level. Here, we obtained 48,821,754 and 55,054,480 clean reads from spleen tissue libraries prepared for bighead carp and silver carp using Illumina paired-end sequencing technology, respectively, and identified 4976 orthologous genes from the transcriptome data sets by comparative analysis. Adaptive evolutionary analysis showed that 212 orthologous genes and 255 Gene Ontology (GO) terms were subjected to positive selection(Ka/Ks values > 1) only in bighead carp, and 195 orthologous genes and 309 GO terms only in silver carp. Among immune defense functions with significant evolutionary divergence, the positively selected biological processes in bighead carp mainly included B cell-mediated immunity, chemokine-mediated signaling pathway, and immunoglobulin mediated immune response, whereas those in silver carp mainly included the antigen processing and presentation, defense response to fungus, and detection of bacteria. Moreover, we found 2974 genes expressed only in spleen of bighead carp and 3494 genes expressed only in spleen of silver carp, where these genes were mostly enriched in the same biological processes or pathways. These results provide a better understanding of the differences in resistance to some diseases by bighead carp and silver carp, and also facilitate the identification of candidate genes related to disease resistance.

Comparative transcriptomic analyses of two bighead carp (Hypophthalmichthys nobilis) groups with different growth rates

Growth is one of the most important and desired economic traits for aquaculture species, and identification of loci controlling growth is a difficult task without genomic sequences. In this study, the liver transcriptomes of two groups (F and S) of bighead carp (Hypophthalmichthys nobilis) within a full-sib family with significant differences in growth rate were sequenced. Following de novo assembly of the combined reads from the two groups, a total of 410 differentially expressed genes were identified. Functional annotation and analysis of these genes indicated that some of these were involved in regulation of glucose levels and lipid metabolism, particularly fatty acid oxidation and transport. In addition to the differences on expression levels between the two groups, we also identified many non-synonymous coding single-nucleotide polymorphisms (SNPs) that were specific to each group, including SNPs from 4 genes involved in the lipid metabolism process (GO: 0006629). These differences in gene expression and DNA sequences may in part comprise the genetic background for the regulation of early growth rate in bighead carp.

Construction of a high density SNP linkage map of kelp (Saccharina japonica) by sequencing Taq I site associated DNA and mapping of a sex determining locus

BACKGROUND

Kelp (Saccharina japonica) has been intensively cultured in China for almost a century. Its genetic improvement is comparable with that of rice. However, the development of its molecular tools is extremely limited, thus its genes, genetics and genomics. Kelp performs an alternative life cycle during which sporophyte generation alternates with gametophyte generation. The gametophytes of kelp can be cloned and crossed. Due to these characteristics, kelp may serve as a reference for the biological and genetic studies of Volvox, mosses and ferns.

RESULTS

We constructed a high density single nucleotide polymorphism (SNP) linkage map for kelp by restriction site associated DNA (RAD) sequencing. In total, 4,994 SNP-containing physical (tag-defined) RAD loci were mapped on 31 linkage groups. The map expanded a total genetic distance of 1,782.75 cM, covering 98.66% of the expected (1,806.94 cM). The length of RAD tags (85 bp) was extended to 400-500 bp with Miseq method, offering us an easiness of developing SNP chips and shifting SNP genotyping to a high throughput track. The number of linkage groups was in accordance with the documented with cytological methods. In addition, we identified a set of microsatellites (99 in total) from the extended RAD tags. A gametophyte sex determining locus was mapped on linkage group 2 in a window about 9.0 cM in width, which was 2.66 cM up to marker_40567 and 6.42 cM down to marker_23595.

CONCLUSIONS

A high density SNP linkage map was constructed for kelp, an intensively cultured brown alga in China. The RAD tags were also extended so that a SNP chip could be developed. In addition, a set of microsatellites were identified among mapped loci, and a gametophyte sex determining locus was mapped. This map will facilitate the genetic studies of kelp including for example the evaluation of germplasm and the decipherment of the genetic bases of economic traits.

A second-generation genetic linkage map for bighead carp (Aristichthys nobilis) based on microsatellite markers

Bighead carp (Aristichthys nobilis) is an important aquaculture fish worldwide. Genetic linkage maps for the species were previously reported, but map resolution remained to be improved. In this study, a second-generation genetic linkage map was constructed for bighead carp through a pseudo-testcross strategy using interspecific hybrids between bighead carp and silver carp. Of the 754 microsatellites genotyped in two interspecific mapping families (with 77 progenies for each family), 659 markers were assigned to 24 linkage groups, which were equal to the chromosome numbers of the haploid genome. The consensus map spanned 1917.3 cM covering 92.8% of the estimated bighead carp genome with an average marker interval of 2.9 cM. The length of linkage groups ranged from 52.2 to 133.5 cM with an average of 79.9 cM. The number of markers per linkage group varied from 11 to 55 with an average of 27.5 per linkage group. Normality tests on interval distances of the map showed a non-normal marker distribution; however, significant correlation was found between the length of linkage group and the number of markers below the 0.01 significance level (two-tailed). The length of the female map was 1.12 times that of the male map, and the average recombination ratio of female to male was 1.10:1. Visual inspection showed that distorted markers gathered in some linkage groups and in certain regions of the male and female maps. This well-defined genetic linkage map will provide a basic framework for further genome mapping of quantitative traits, comparative mapping and marker-assisted breeding in bighead carp.

Incorporation of genetic technologies associated with applied reproductive technologies to enhance world food production

Animal breeding and reproductive physiology have been closely related throughout the history of animal production science. Artificial insemination provides the best method of increasing the influence of sires with superior genetics to improve production traits. Multiple ovulation embryo transfer (MOET) provides some ability to increase the genetic influence of the maternal line as well. The addition of genetic technologies to this paradigm allows for improved methods of selecting sires and dams carrying the best genes for production and yield of edible products and resistance to diseases and parasites. However, decreasing the number of influential parents within a population also increases the risk of propagating a recessive gene that could negatively impact the species (Reprod Domest Anim 44:792-796, 2009; BMC Genomics 11:337, 2010). Furthermore, antagonistic genotypic relationships between production traits and fertility (Anim Prod Sci 49:399-412, 2009; Anim Genet 43:442-446, 2012) suggest that care must be taken to ensure that increasing the frequency of genes with a positive influence on production does not negatively impact the fertility of the replacement females entering the herd.

Massively parallel sequencing of Chikso (Korean brindle cattle) to discover genome-wide SNPs and InDels

Since the completion of the bovine sequencing projects, a substantial number of genetic variations such as single nucleotide polymorphisms have become available across the cattle genome. Recently, cataloguing such genetic variations has been accelerated using massively parallel sequencing technology. However, most of the recent studies have been concentrated on European Bos taurus cattle breeds, resulting in a severe lack of knowledge for valuable native cattle genetic resources worldwide. Here, we present the first whole-genome sequencing results for an endangered Korean native cattle breed, Chikso, using the Illumina HiSeq 2,000 sequencing platform. The genome of a Chikso bull was sequenced to approximately 25.3-fold coverage with 98.8% of the bovine reference genome sequence (UMD 3.1) covered. In total, 5,874,026 single nucleotide polymorphisms and 551,363 insertion/deletions were identified across all 29 autosomes and the X-chromosome, of which 45% and 75% were previously unknown, respectively. Most of the variations (92.7% of single nucleotide polymorphisms and 92.9% of insertion/deletions) were located in intergenic and intron regions. A total of 16,273 single nucleotide polymorphisms causing missense mutations were detected in 7,111 genes throughout the genome, which could potentially contribute to variation in economically important traits in Chikso. This study provides a valuable resource for further investigations of the genetic mechanisms underlying traits of interest in cattle, and for the development of improved genomics-based breeding tools.

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

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

Transcriptome analysis of silver carp (Hypophthalmichthys molitrix) by paired-end RNA sequencing

The silver carp (Hypophthalmichthys molitrix) is among the most intensively pond-cultured fish species and is used in the wild to counteract water bloom in China. However, little genomic information is available for this species, especially regarding its ability to grow rapidly in water, even water contaminated with high concentrations of poisonous microcystin. In this study, we performed de novo transcriptome assembly and analysis of the 17.10 million short-read sequences produced by the Illumina paired-end sequencing technology. Using an improved multiple k-mer contig assembly method coupled with further scaffolding, 85,759 sequences were obtained. There were 23,044 sequences annotated with 3423 gene ontology terms for 104 196 term occurrences and the three corresponding organizing principles. A total of 38,200 assembled sequences were involved in 218 predicted Kyoto Encyclopedia of Genes and Genomes metabolic pathways. We also recovered 41 of 44 genes involved in the biosynthesis of glutathione. Of these, five genes were identified as experienced positive selection between silver carp and zebrafish, as determined by the likelihood ratio test. This report is the first annotated review of the silver carp transcriptome. These data will be of interest to researchers investigating the evolution and biological processes of the silver carp. This work also provides an archive for future studies of recent speciation and evolution of Cyprinidae fishes and can be used in comparative studies of other fishes.

Association of an ACSL1 gene variant with polyunsaturated fatty acids in bovine skeletal muscle

Background

The intramuscular fat deposition and the fatty acid profiles of beef affect meat quality. High proportions of unsaturated fatty acids are related to beef flavor and are beneficial for the nutritional value of meat. Moreover, a variety of clinical and epidemiologic studies showed that particularly long-chain omega-3 fatty acids from animal sources have a positive impact on human health and disease.

Results

To screen for genetic factors affecting fatty acid profiles in beef, we initially performed a microsatellite-based genome scan in a F₂ Charolais × German Holstein resource population and identified a quantitative trait locus (QTL) for fatty acid composition in a region on bovine chromosome 27 where previously QTL affecting marbling score had been detected in beef cattle populations. The long-chain acyl-CoA synthetase 1 (ACSL1) gene was identified as the most plausible functional and positional candidate gene in the QTL interval due to its direct impact on fatty acid metabolism and its position in the QTL interval. ACSL1 is necessary for synthesis of long-chain acyl-CoA esters, fatty acid degradation and phospholipid remodeling. We validated the genomic annotation of the bovine ACSL1 gene by in silico comparative sequence analysis and experimental verification. Re-sequencing of the complete coding, exon-flanking intronic sequences, 3' untranslated region (3'UTR) and partial promoter region of the ACSL1 gene revealed three synonymous mutations in exons 6, 7, and 20, six noncoding intronic gene variants, six polymorphisms in the promoter region, and four variants in the 3' UTR region. The association analysis identified the gene variant in intron 5 of the ACSL1 gene (c.481-233A>G) to be significantly associated with the relative content of distinct fractions and ratios of fatty acids (e.g., n-3 fatty acids, polyunsaturated, n-3 long-chain polyunsaturated fatty acids, trans vaccenic acid) in skeletal muscle. A tentative association of the ACSL1 gene variant with intramuscular fat content indicated that an indirect effect on fatty acid composition via modulation of total fat content of skeletal muscle cannot be excluded.

Conclusions

The initial QTL analysis suggested the ACSL1 gene as a positional and functional candidate gene for fatty acid composition in bovine skeletal muscle. The findings of subsequent association analyses indicate that ACSL1 or a separate gene in close proximity might play a functional role in mediating the lipid composition of beef.