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

Construction of Three High-Density Genetic Linkage Maps and Dynamic QTL Mapping of Growth Traits in Yellow River Carp (Cyprinus carpio haematopterus)

To provide the theoretical basis for researching growth, development, and molecular marker-assisted breeding of the economically important Yellow River carp (Cyprinus carpio haematopterus) using dynamic quantitative trait locus (QTL) mapping, we constructed three genetic linkage maps from 207 progeny using a new modified genotyping-by-sequencing method. The three maps contained 16,886, 16,548, and 7482 single nucleotide polymorphism markers, respectively, with an average interval of 0.36 cM, 0.45 cM, and 1.00 cM. We identified 148 QTLs related to four growth traits that were located on 25 chromosomes from three growth stages of Yellow River carp. A total of 32, 36, 43, and 37 QTLs were associated with body length, height, width, and weight, respectively. Among them, 47 QTLs were detected for only one growth trait in one stage, but all of the other QTLs were co-localized. Of the 14 main QTLs, 13 were located on chromosome 12, which suggests the presence of growth-related genes on this chromosome. We then detected 17 candidate genes within 50 K upstream and downstream of the 14 main QTLs. This is the first report of the dynamic QTL mapping of growth traits of Yellow River carp, and the results can be used in future studies of growth, development, and molecular-assisted breeding of this species.

A High-Density Genetic Linkage Map and Fine Mapping of QTL For Feed Conversion Efficiency in Common Carp (Cyprinus carpio)

Feed conversion efficiency (FCE) is an economically crucial trait in fish, however, little progress has been made in genetics and genomics for this trait because phenotypes of the trait are difficult to measure. In this study, we constructed a high-density and high-resolution genetic linkage map with 28,416 SNP markers for common carp (Cyprinus carpio) based on high throughput genotyping with the carp 250K single nucleotide polymorphism (SNP) array in a full-sib F1 family of mirror carp (Cyprinus carpio) consisting of 141 progenies. The linkage map contained 11,983 distinct loci and spanned 3,590.09 cM with an average locus interval of 0.33 cM. A total of 17 QTL for the FCE trait were detected on four LGs (LG9, LG20, LG28, and LG32), explaining 8.9-15.9% of the phenotypic variations. One major cluster containing eight QTL (qFCE1-28, qFCE2-28, qFCE3-28, qFCE4-28, qFCE5-28, qFCE6-28, qFCE7-28, and qFCE8-28) was detected on LG28. Two clusters consisting of four QTL (qFCE1-32, qFCE2-32, qFCE3-32, and qFCE4-32) and three QTL (qFCE1-20, qFCE2-20, and qFCE3-20) were detected on LG32 and LG20, respectively. Nine candidate genes (ACACA, SCAF4, SLC2A5, TNMD, PCDH1, FOXO, AGO1, FFAR3, and ARID1A) underlying the feed efficiency trait were also identified, the biological functions of which may be involved in lipid metabolism, carbohydrate metabolism, energy deposition, fat accumulation, digestion, growth regulation, and cell proliferation and differentiation according to GO (Gene Ontology). As an important tool, high-density and high-resolution genetic linkage maps play a crucial role in the QTL fine mapping of economically important traits. Our novel findings provided new insights that elucidate the genetic basis and molecular mechanism of feed efficiency and the subsequent marker-assisted selection breeding in common carp.

The stage-specific long non-coding RNAs and mRNAs identification and analysis during early development of common carp, Cyprinus carpio

Cyprinus carpio is considered an alternative vertebrate fish model to zebrafish. However, systemic times-series research on the lncRNAs and mRNAs during early development of C. carpio has not been reported yet. This study provides the first long non-coding RNA (lncRNA)-mRNA expression profiles during six main early development stages (2 h post-fertilization hpf, 6 hpf, 12 hpf, 20 hpf, 64 hpf and 1 day post-hatching). A total of 51,979 lncRNAs were identified. We screened the top 10 abundance lncRNAs and mRNAs and stage-specific lncRNAs and mRNAs (specificity measure SPM > 0.9). We identified significant differentially expressed lncRNAs and mRNAs (|log2 (fold change)| ≥ 1 and false discovery rate FDR of <0.05). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis identified numerous signaling pathways. Additionally, the lncRNA-mRNA co-regulated network analysis of two lncRNAs (lncrps25 and malat1) and two mRNAs (mitf and troponin T) were investigated. Our results provide new insight into the role of lncRNAs and mRNAs, and would advance the understanding of lncRNA-mediated mechanisms in early development of fish.

First High-Density Linkage Map and QTL Fine Mapping for Growth-Related Traits of Spotted Sea bass (Lateolabrax maculatus)

Possessing powerful adaptive capacity and a pleasant taste, spotted sea bass (Lateolabrax maculatus) has a broad natural distribution and is one of the most popular mariculture fish in China. However, the genetic improvement program for this fish is still in its infancy. Growth is the most economically important trait and is controlled by quantitative trait loci (QTL); thus, the identification of QTLs and genetic markers for growth-related traits is an essential step for the establishment of marker-assisted selection (MAS) breeding programs. In this study, we report the first high-density linkage map of spotted sea bass constructed by sequencing 333 F1 generation individuals in a full-sib family using 2b-RAD technology. A total of 6883 SNP markers were anchored onto 24 linkage groups, spanning 2189.96 cM with an average marker interval of 0.33 cM. Twenty-four growth-related QTLs, including 13 QTLs for body weight and 11 QTLs for body length, were successfully detected, with phenotypic variance explained (PVE) ranging from 5.1 to 8.6%. Thirty potential candidate growth-related genes surrounding the associated SNPs were involved in cell adhesion, cell proliferation, cytoskeleton reorganization, calcium channels, and neuromodulation. Notably, the fgfr4 gene was detected in the most significant QTL; this gene plays a pivotal role in myogenesis and bone growth. The results of this study may facilitate marker-assisted selection for breeding populations and establish the foundation for further genomic and genetic studies investigating spotted sea bass.

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.

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.

Construction of a High-Density Linkage Map and QTL Fine Mapping for Growth- and Sex-Related Traits in Channel Catfish (Ictalurus punctatus)

A high-density genetic linkage map is of particular importance in the fine mapping for important economic traits and whole genome assembly in aquaculture species. The channel catfish (Ictalurus punctatus), a species native to North America, is one of the most important commercial freshwater fish in the world. Outside of the United States, China has become the major producer and consumer of channel catfish after experiencing rapid development in the past three decades. In this study, based on restriction site associated DNA sequencing (RAD-seq), a high-density genetic linkage map of channel catfish was constructed by using single nucleotide polymorphisms (SNPs) in a F₁ family composed of 156 offspring and their two parental individuals. A total of 4,768 SNPs were assigned to 29 linkage groups (LGs), and the length of the linkage map reached 2,480.25 centiMorgans (cM) with an average distance of 0.55 cM between loci. Based on this genetic linkage map, 223 genomic scaffolds were anchored to the 29 LGs of channel catfish, and a total length of 704.66 Mb was assembled. Quantitative trait locus (QTL) mapping and genome-wide association analysis identified 10 QTLs of sex-related and six QTLs of growth-related traits at LG17 and LG28, respectively. Candidate genes associated with sex dimorphism, including spata2, spata5, sf3, zbtb38, and fox, were identified within QTL intervals on the LG17. A sex-linked marker with simple sequence repeats (SSR) in zbtb38 gene of the LG17 was validated for practical verification of sex in the channel catfish. Thus, the LG17 was considered as a sex-related LG. Potential growth-related genes were also identified, including important regulators such as megf9, npffr1, and gas1. In a word, we constructed the high-density genetic linkage map and developed the sex-linked marker in channel catfish, which are important genetic resources for future marker-assisted selection (MAS) of this economically important teleost.

Genome-Scale Association Study of Abnormal Scale Pattern in Yellow River Carp Identified Previously Known Causative Gene in European Mirror Carp

Common carp (Cyprinus carpio) is one of the most widely studied fish species due to its great economic value and strong environmental adaptability. Scattered scale, a typical phenotype of the mirror carp that is derived from Europe, has never been observed in the Yellow River carp previously. We recently identified approximately one fourth of the F1 progenies displaying scattered scale in a full-sib Yellow River carp family in our breeding program, despite both parents that showed wild type with normal scale patterns. This family provides us unique materials to investigate the genetic basis underlying the abnormal scale mutant in Yellow River carp population. Genome-wide association study (GWAS) and association mapping were performed based on genome-wide single nucleotide polymorphisms (SNP) genotyped with common carp 250 K SNP genotyping array in 82 samples of the Yellow River carp family. We identified a 1.4 Mb genome region that was significantly associated with abnormal scattered scale patterns. We further identified a deletion mutation in fibroblast growth factor receptor 1 a1 (fgfr1a1) gene within this genome region. Amplification and sequencing analysis of this gene revealed a 311-bp deletion in intron 10 and exon 11, which proved that fgfr1a1 could be the causal gene responsible for abnormal scattered scale in the Yellow River carp family. Since similar fragment mutation with 306-bp and 310-bp deletions had been previously reported as causal mutation of scattered scale patterns in the mirror carp, we speculate that either the deletion mutation was introduced from Europe-derived mirror carp or the deletion independently occurred in the mutation hotspot in fgfr1a1 gene. The results provided insights into the genetic basis of scale pattern mutant in Yellow River carp population, which would help us to eliminate the recessive allele of the abnormal scale patterns in Yellow River carp population by molecular marker-assisted breeding.

Genome-Wide Association Analysis With a 50K Transcribed Gene SNP-Chip Identifies QTL Affecting Muscle Yield in Rainbow Trout

Detection of coding/functional SNPs that change the biological function of a gene may lead to identification of putative causative alleles within QTL regions and discovery of genetic markers with large effects on phenotypes. This study has two-fold objectives, first to develop, and validate a 50K transcribed gene SNP-chip using RNA-Seq data. To achieve this objective, two bioinformatics pipelines, GATK and SAMtools, were used to identify ~21K transcribed SNPs with allelic imbalances associated with important aquaculture production traits including body weight, muscle yield, muscle fat content, shear force, and whiteness in addition to resistance/susceptibility to bacterial cold-water disease (BCWD). SNPs ere identified from pooled RNA-Seq data collected from ~620 fish, representing 98 families from growth- and 54 families from BCWD-selected lines with divergent phenotypes. In addition, ~29K transcribed SNPs without allelic-imbalances were strategically added to build a 50K Affymetrix SNP-chip. SNPs selected included two SNPs per gene from 14K genes and ~5K non-synonymous SNPs. The SNP-chip was used to genotype 1728 fish. The average SNP calling-rate for samples passing quality control (QC; 1,641 fish) was ≥ 98.5%. The second objective of this study was to test the feasibility of using the new SNP-chip in GWA (Genome-wide association) analysis to identify QTL explaining muscle yield variance. GWA study on 878 fish (representing 197 families from 2 consecutive generations) with muscle yield phenotypes and genotyped for 35K polymorphic markers (passing QC) identified several QTL regions explaining together up to 28.40% of the additive genetic variance for muscle yield in this rainbow trout population. The most significant QTLs were on chromosomes 14 and 16 with 12.71 and 10.49% of the genetic variance, respectively. Many of the annotated genes in the QTL regions were previously reported as important regulators of muscle development and cell signaling. No major QTLs were identified in a previous GWA study using a 57K genomic SNP chip on the same fish population. These results indicate improved detection power of the transcribed gene SNP-chip in the target trait and population, allowing identification of large-effect QTLs for important traits in rainbow trout.

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.

Fine-mapping of genes determining extrafusal fiber properties in murine soleus muscle

Muscle fiber cross-sectional area (CSA) and proportion of different fiber types are important determinants of muscle function and overall metabolism. Genetic variation plays a substantial role in phenotypic variation of these traits; however, the underlying genes remain poorly understood. This study aimed to map quantitative trait loci (QTL) affecting differences in soleus muscle fiber traits between the LG/J and SM/J mouse strains. Fiber number, CSA, and proportion of oxidative type I fibers were assessed in the soleus of 334 genotyped female and male mice of the F₃₄ generation of advanced intercross lines (AIL) derived from the LG/J and SM/J strains. To increase the QTL detection power, these data were combined with 94 soleus samples from the F₂ intercross of the same strains. Transcriptome of the soleus muscle of LG/J and SM/J females was analyzed by microarray. Genome-wide association analysis mapped four QTL (genome-wide P < 0.05) affecting the properties of muscle fibers to chromosome 2, 3, 4, and 11. A 1.5-LOD QTL support interval ranged between 2.36 and 4.67 Mb. On the basis of the genomic sequence information and functional and transcriptome data, we identified candidate genes for each of these QTL. The combination of analyses in F₂ and F₃₄ AIL populations with transcriptome and genomic sequence data in the parental strains is an effective strategy for refining QTL and nomination of the candidate genes.

Genome-Wide Association Study for Muscle Fat Content and Abdominal Fat Traits in Common Carp (Cyprinus carpio)

Muscle fat content is an important phenotypic trait in fish, as it affects the nutritional, technical and sensory qualities of flesh. To identify loci and candidate genes associated with muscle fat content and abdominal fat traits, we performed a genome-wide association study (GWAS) using the common carp 250 K SNP assay in a common carp F₂ resource population. A total of 18 loci surpassing the genome-wide suggestive significance level were detected for 4 traits: fat content in dorsal muscle (MFdo), fat content in abdominal muscle (MFab), abdominal fat weight (AbFW), and AbFW as a percentage of eviscerated weight (AbFP). Among them, one SNP (carp089419) affecting both AbFW and AbFP reached the genome-wide significance level. Ten of those loci were harbored in or near known genes. Furthermore, relative expressions of 5 genes related to MFdo were compared using dorsal muscle samples with high and low phenotypic values. The results showed that 4 genes were differentially expressed between the high and low phenotypic groups. These genes are, therefore, prospective candidate genes for muscle fat content: ankyrin repeat domain 10a (ankrd10a), tetratricopeptide repeat, ankyrin repeat and coiled-coil containing 2 (tanc2), and four jointed box 1 (fjx1) and choline kinase alpha (chka). These results offer valuable insights into the complex genetic basis of fat metabolism and deposition.

Quantitative trait loci for morphometric traits in multiple families of common carp (Cyprinus carpio)

Quantitative trait locus (QTL) mapping is frequently used to understand the genetic architecture of quantitative traits. Herein, we performed a genome scan for QTL affecting the morphometric characters in eight full-sib families containing 522 individuals using different statistical methods (Sib-pair and half-sib model). A total of 194 QTLs were detected in 25 different regions on 10 linkage groups (LGs). Among them, 37 QTLs on five LGs (eight, 13, 24, 40 and 45) were significant (5% genome-wide level), while the remaining 40 (1% chromosome-wide level) and 117 (5% chromosome-wide level) indicated suggestive effect on those traits. Heritabilities for most morphometric traits were moderate to high, ranging from 0.21 to 0.66, with generally strong phenotypic and genetic correlations between the traits. A large number of QTLs for morphometric traits were co-located, consistent with their high correlations, and may reflect pleiotropic effect on the same genes. Biological pathways were mapped for possible candidate genes on QTL regions. One significantly enriched pathway was identified on LG45, which had a P-value of 0.04 and corresponded to the "regulation of actin cytoskeleton pathway". The results are expected to be useful in marker-assisted selection (MAS) and provide valuable information for the study of gene pathway for morphometric and growth traits of the common carp.

An ultra-high density linkage map and QTL mapping for sex and growth-related traits of common carp (Cyprinus carpio)

High density genetic linkage maps are essential for QTL fine mapping, comparative genomics and high quality genome sequence assembly. In this study, we constructed a high-density and high-resolution genetic linkage map with 28,194 SNP markers on 14,146 distinct loci for common carp based on high-throughput genotyping with the carp 250 K single nucleotide polymorphism (SNP) array in a mapping family. The genetic length of the consensus map was 10,595.94 cM with an average locus interval of 0.75 cM and an average marker interval of 0.38 cM. Comparative genomic analysis revealed high level of conserved syntenies between common carp and the closely related model species zebrafish and medaka. The genome scaffolds were anchored to the high-density linkage map, spanning 1,357 Mb of common carp reference genome. QTL mapping and association analysis identified 22 QTLs for growth-related traits and 7 QTLs for sex dimorphism. Candidate genes underlying growth-related traits were identified, including important regulators such as KISS2, IGF1, SMTLB, NPFFR1 and CPE. Candidate genes associated with sex dimorphism were also identified including 3KSR and DMRT2b. The high-density and high-resolution genetic linkage map provides an important tool for QTL fine mapping and positional cloning of economically important traits, and improving common carp genome assembly.

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.

Targeted disruption of sp7 and myostatin with CRISPR-Cas9 results in severe bone defects and more muscular cells in common carp

The common carp (Cyprinus carpio) as one of the most important aquaculture fishes produces over 3 million metric tones annually, approximately 10% the annual production of the all farmed freshwater fish worldwide. However, the tetraploidy genome and long generation-time of the common carp have made its breeding and genetic studies extremely difficult. Here, TALEN and CRISPR-Cas9, two versatile genome-editing tools, are employed to target common carp bone-related genes sp7, runx2, bmp2a, spp1, opg, and muscle suppressor gene mstn. TALEN were shown to induce mutations in the target coding sites of sp7, runx2, spp1 and mstn. With CRISPR-Cas9, the two common carp sp7 genes, sp7a and sp7b, were mutated individually, all resulting in severe bone defects; while mstnba mutated fish have grown significantly more muscle cells. We also employed CRISPR-Cas9 to generate double mutant fish of sp7a;mstnba with high efficiencies in a single step. These results demonstrate that both TALEN and CRISPR-Cas9 are highly efficient tools for modifying the common carp genome, and open avenues for facilitating common carp genetic studies and breeding.

Construction of the BAC Library of Small Abalone (Haliotis diversicolor) for Gene Screening and Genome Characterization

The small abalone (Haliotis diversicolor) is one of the most important aquaculture species in East Asia. To facilitate gene cloning and characterization, genome analysis, and genetic breeding of it, we constructed a large-insert bacterial artificial chromosome (BAC) library, which is an important genetic tool for advanced genetics and genomics research. The small abalone BAC library includes 92,610 clones with an average insert size of 120 Kb, equivalent to approximately 7.6× of the small abalone genome. We set up three-dimensional pools and super pools of 18,432 BAC clones for target gene screening using PCR method. To assess the approach, we screened 12 target genes in these 18,432 BAC clones and identified 16 positive BAC clones. Eight positive BAC clones were then sequenced and assembled with the next generation sequencing platform. The assembled contigs representing these 8 BAC clones spanned 928 Kb of the small abalone genome, providing the first batch of genome sequences for genome evaluation and characterization. The average GC content of small abalone genome was estimated as 40.33%. A total of 21 protein-coding genes, including 7 target genes, were annotated into the 8 BACs, which proved the feasibility of PCR screening approach with three-dimensional pools in small abalone BAC library. One hundred fifty microsatellite loci were also identified from the sequences for marker development in the future. The BAC library and clone pools provided valuable resources and tools for genetic breeding and conservation of H. diversicolor.

De novo assembly and characterization of the spleen transcriptome of common carp (Cyprinus carpio) using Illumina paired-end sequencing

Common carp (Cyprinus carpio) is one of the most important aquacultured species of the family Cyprinidae, and breeding this species for disease resistance is becoming more and more important. However, at the genome or transcriptome levels, study of the immunogenetics of disease resistance in the common carp is lacking. In this study, 60,316,906 and 75,200,328 paired-end clean reads were obtained from two cDNA libraries of the common carp spleen by Illumina paired-end sequencing technology. Totally, 130,293 unique transcript fragments (unigenes) were assembled, with an average length of 1400.57 bp. Approximately 105,612 (81.06%) unigenes could be annotated according to their homology with matches in the Nr, Nt, Swiss-Prot, COG, GO, or KEGG databases, and they were found to represent 46,747 non-redundant genes. Comparative analysis showed that 59.82% of the unigenes have significant similarity to zebrafish Refseq proteins. Gene expression comparison revealed that 10,432 and 6889 annotated unigenes were, respectively, up- and down-regulated with at least twofold changes between two developmental stages of the common carp spleen. Gene ontology and KEGG analysis were performed to classify all unigenes into functional categories for understanding gene functions and regulation pathways. In addition, 46,847 simple sequence repeats (SSRs) were detected from 35,618 unigenes, and a large number of single nucleotide polymorphism (SNP) and insertion/deletion (INDEL) sites were identified in the spleen transcriptome of common carp. This study has characterized the spleen transcriptome of the common carp for the first time, providing a valuable resource for a better understanding of the common carp immune system and defense mechanisms. This knowledge will also facilitate future functional studies on common carp immunogenetics that may eventually be applied in breeding programs.

Comparative mapping for bighead carp (Aristichthys nobilis) against model and non-model fishes provides insights into the genomic evolution of cyprinids

Comparative mapping provides an efficient method to connect genomes of non-model and model fishes. In this study, we used flanking sequences of the 659 microsatellites on a genetic map of bighead carp (Aristichthys nobilis) to comprehensively study syntenic relationships between bighead carp and nine model and non-model fishes. Of the five model and two food fishes with whole genome data, Cyprinus carpio showed the highest rate of positive BLAST hits (95.3 %) with bighead carp map, followed by Danio rerio (70.9 %), Oreochromis niloticus (21.7 %), Tetraodon nigroviridis (6.4 %), Gasterosteus aculeatus (5.2 %), Oryzias latipes (4.7 %) and Fugu rubripes (3.5 %). Chromosomal syntenic analyses showed that inversion was the basic chromosomal rearrangement during genomic evolution of cyprinids, and the extent of inversions and translocations was found to be positively correlated with evolutionary relationships among fishes studied. Among the five investigated cyprinids, linkage groups (LGs) of bighead carp, Hypophthalmichthys molitrix and Ctenopharyngodon idella exhibited a one-to-one relationship. Besides, LG 9 of bighead carp and homologous LGs of silver carp and grass carp all corresponded to the chromosomes 10 and 22 of zebrafish, suggesting that chromosomal fission may have occurred in the ancestor of zebrafish. On the other hand, LGs of bighead carp and common carp showed an approximate one-to-two relationship with extensive translocations, confirming the occurrence of a 4th whole genome duplication in common carp. This study provides insights into the understanding of genome evolution among cyprinids and would aid in transferring positional and functional information of genes from model fish like zebrafish to non-model fish like bighead carp.

A novel 31-bp indel in the paired box 7 (PAX7) gene is associated with chicken performance traits

1. A novel 31-bp indel polymorphism in intron 3 of the chicken paired box 7 (PAX7) gene was identified and genotyped in an F2 resource population of Gushi chicken crossed with Anka broiler to analyse its associations with chicken growth, carcass and meat quality traits. 2. Results showed that the 31-bp indel was significantly associated with body weight at 4, 6, 8, 10, 12 weeks of age and body size indices including shank length, shank girth, body slanting length at 8 and 12 weeks of age. Significant associations were found for carcass weight, semi-evisceration weight, evisceration weight, breast muscle fibre diameter, leg muscle fibre diameter, breast muscle fibre density, while no significant association with leg muscle fibre density was observed. 3. It was concluded that the 31-bp indel in intron 3 of the PAX7 gene was associated with chicken growth, carcass and meat quality traits where the 31-bp deletion had a negative effect on chicken growth and carcass traits and positive effect on meat quality traits.

Novel single nucleotide polymorphisms of the insulin-like growth factor-I gene and their associations with growth traits in common carp (Cyprinus carpio L.)

Insulin-like growth factor-I (IGF-I) plays an important role in the growth and development of vertebrates. To study polymorphisms of IGF-I, we screened a total of 4555 bp of genomic sequences in four exons and partial introns for the discovery of single nucleotide polymorphism (SNP) in common carp (Cyprinus carpio). Three SNPs (g.3759T>G, g.7627T>A and g.7722T>C) in intron 2 and a nonsynonymous SNP (g.7892C>T) in exon 3 were identified in a pilot population including random parents and their progenies. 289 progenies were further genotyped for studying possible associations between genotypes or combined genotypes and growth traits. The results showed that the locus g.7627T>A was significantly associated with body weight and body length, and fish with genotype AA had a mean body weight 5.9% higher than those with genotype TT. No significant associations were observed between genotypes of other loci and growth traits. However, when both g.7627T>A and g.7722T>C were considered, the combined genotype TT/TT was extremely associated with the lowest values of body length and body weight and the highest K value in comparison with other diplotypes (p < 0.01). These results suggest that genotype AA at g.7627T>A and its combined genotypes with alleles from another locus have positive effects on growth traits, which would be a candidate molecular marker for further studies in marker-assisted selection in common carp.

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.

Screening of SSR markers associated with scale cover pattern and mapped to a genetic linkage map of common carp (Cyprinus carpio L.)

Fish scale is an attractive model in bone physiology research and is also a crucial character for breeding new varieties. Thus, it is important to identify loci in the genome associated with scale formation. In this study, 290 microsatellite markers in common carp (Cyprinus carpio L.) were selected and tested for their segregation in a full-sib mapping panel containing 96 individuals (population 1). Association analysis identified seven simple sequence repeats (SSRs) (HLJ2509, HLJ3227, HLJ3675, HLJ3766, HLJ3863, FGFR1, FGFR7) that showed significant correlation with scale cover pattern in population 1. When the seven SSRs were investigated in two other populations, seven and five SSRs were significantly correlated with scale cover pattern in population 2 (116 individuals) and population 3 (57 individuals), respectively. The exceptions were FGFR1 and HLJ3227. A genetic linkage map was constructed using the 290 SSRs and 241 SSRs were mapped into 47 linkage groups (LGs), with 2-15 markers per LG. The map spanned 2,241.7 cM, with LG sizes that varied from 1.1 to 124.9 cM. All seven markers associated with scale cover mapped into LG3. We considered that a gene cluster that affected the scale cover pattern possibly existed in LG3. By aligning the seven markers with the zebrafish (Danio rerio) genome, we identified six candidate genes (atoh1a, ptch1, bmp1a, fgfr1a, fgf17, wnt5a) that may be associated with scale formation. We propose that the seven markers could be used with marker-assisted selection to breed a new variety of common carp, and the six candidate genes could help in understanding the scale cover mechanism.

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.

A linkage map of transcribed single nucleotide polymorphisms in rohu (Labeo rohita) and QTL associated with resistance to Aeromonas hydrophila

BACKGROUND

Production of carp dominates world aquaculture. More than 1.1 million tonnes of rohu carp, Labeo rohita (Hamilton), were produced in 2010. Aeromonas hydrophila is a bacterial pathogen causing aeromoniasis in rohu, and is a major problem for carp production worldwide. There is a need to better understand the genetic mechanisms affecting resistance to this disease, and to develop tools that can be used with selective breeding to improve resistance. Here we use a 6 K SNP array to genotype 21 full-sibling families of L. rohita that were experimentally challenged intra-peritoneally with a virulent strain of A. hydrophila to scan the genome for quantitative trait loci associated with disease resistance.

RESULTS

In all, 3193 SNPs were found to be informative and were used to create a linkage map and to scan for QTL affecting resistance to A. hydrophila. The linkage map consisted of 25 linkage groups, corresponding to the number of haploid chromosomes in L. rohita. Male and female linkage maps were similar in terms of order, coverage (1384 and 1393 cM, respectively) and average interval distances (1.32 and 1.35 cM, respectively). Forty-one percent of the SNPs were annotated with gene identity using BLAST (cut off E-score of 0.001). Twenty-one SNPs mapping to ten linkage groups showed significant associations with the traits hours of survival and dead or alive (P <0.05 after Bonferroni correction). Of the SNPs showing significant or suggestive associations with the traits, several were homologous to genes of known immune function or were in close linkage to such genes. Genes of interest included heat shock proteins (70, 60, 105 and "small heat shock proteins"), mucin (5b precursor and 2), lectin (receptor and CD22), tributyltin-binding protein, major histocompatibility loci (I and II), complement protein component c7-1, perforin 1, ubiquitin (ligase, factor e4b isoform 2 and conjugation enzyme e2 c), proteasome subunit, T-cell antigen receptor and lymphocyte specific protein tyrosine kinase.

CONCLUSIONS

A panel of markers has been identified that will be validated for use with both genomic and marker-assisted selection to improve resistance of L. rohita to A. hydrophila.

Development and evaluation of the first high-throughput SNP array for common carp (Cyprinus carpio)

Background

A large number of single nucleotide polymorphisms (SNPs) have been identified in common carp (Cyprinus carpio) but, as yet, no high-throughput genotyping platform is available for this species. C. carpio is an important aquaculture species that accounts for nearly 14% of freshwater aquaculture production worldwide. We have developed an array for C. carpio with 250,000 SNPs and evaluated its performance using samples from various strains of C. carpio.

Results

The SNPs used on the array were selected from two resources: the transcribed sequences from RNA-seq data of four strains of C. carpio, and the genome re-sequencing data of five strains of C. carpio. The 250,000 SNPs on the resulting array are distributed evenly across the reference C.carpio genome with an average spacing of 6.6 kb. To evaluate the SNP array, 1,072 C. carpio samples were collected and tested. Of the 250,000 SNPs on the array, 185,150 (74.06%) were found to be polymorphic sites. Genotyping accuracy was checked using genotyping data from a group of full-siblings and their parents, and over 99.8% of the qualified SNPs were found to be reliable. Analysis of the linkage disequilibrium on all samples and on three domestic C.carpio strains revealed that the latter had the longer haplotype blocks. We also evaluated our SNP array on 80 samples from eight species related to C. carpio, with from 53,526 to 71,984 polymorphic SNPs. An identity by state analysis divided all the samples into three clusters; most of the C. carpio strains formed the largest cluster.

Conclusions

The Carp SNP array described here is the first high-throughput genotyping platform for C. carpio. Our evaluation of this array indicates that it will be valuable for farmed carp and for genetic and population biology studies in C. carpio and related species.

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

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

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.

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.

A consensus linkage map provides insights on genome character and evolution in common carp (Cyprinus carpio L.)

Common carp (Cyprinus carpio L.) is cultured worldwide and is a major contributor to the world's aquaculture production. The common carp has a complex tetraploidized genome, which may historically experience additional whole genome duplication than most other Cyprinids. Fine maps for female and male carp were constructed using a mapping panel containing one F1 family with 190 progeny. A total of 1,025 polymorphic markers were used to construct genetic maps. For the female map, 559 microsatellite markers in 50 linkage groups cover 3,468 cM of the genome. For the male map, 383 markers in 49 linkage groups cover 1,811 cM of the genome. The consensus map was constructed by integrating the new map with two published linkage maps, containing 732 markers and spanning 3,278 cM in 50 linkage groups. The number of consensus linkage groups corresponds to the number of common carp chromosomes. A significant difference on sex recombinant rate was observed that the ratio of female and male recombination rates was 4.2:1. Comparative analysis was performed between linkage map of common carp and genome of zebrafish (Danio rerio), which revealed clear 2:1 relationship of common carp linkage groups and zebrafish chromosomes. The results provided evidence that common carp did experienced a specific whole genome duplication event comparing with most other Cyprinids. The consensus linkage map provides an important tool for genetic and genome study of common carp and facilitates genetic selection and breeding for common carp industry.

A dense genetic linkage map for common carp and its integration with a BAC-based physical map

Background

Common carp (Cyprinus carpio) is one of the most important aquaculture species with an annual global production of 3.4 million metric tons. It is also an important ornamental species as well as an important model species for aquaculture research. To improve the economically important traits of this fish, a number of genomic resources and genetic tools have been developed, including several genetic maps and a bacterial artificial chromosome (BAC)-based physical map. However, integrated genetic and physical maps are not available to study quantitative trait loci (QTL) and assist with fine mapping, positional cloning and whole genome sequencing and assembly. The objective of this study was to integrate the currently available BAC-based physical and genetic maps.

Results

The genetic map was updated with 592 novel markers, including 312 BAC-anchored microsatellites and 130 SNP markers, and contained 1,209 genetic markers on 50 linkage groups, spanning 3,565.9 cM in the common carp genome. An integrated genetic and physical map of the common carp genome was then constructed, which was composed of 463 physical map contigs and 88 single BACs. Combined lengths of the contigs and single BACs covered a physical length of 498.75 Mb, or around 30% of the common carp genome. Comparative analysis between common carp and zebrafish genomes was performed based on the integrated map, providing more insights into the common carp specific whole genome duplication and segmental rearrangements in the genome.

Conclusion

We integrated a BAC-based physical map to a genetic linkage map of common carp by anchoring BAC-associated genetic markers. The density of the genetic linkage map was significantly increased. The integrated map provides a tool for both genetic and genomic studies of common carp, which will help us to understand the genomic architecture of common carp and facilitate fine mapping and positional cloning of economically important traits for genetic improvement and modification.

A consensus linkage map of common carp (Cyprinus carpio L.) to compare the distribution and variation of QTLs associated with growth traits

The ability to detect and identify quantitative trait loci (QTLs) in a single population is often limited. Analyzing multiple populations in QTL analysis improves the power of detecting QTLs and provides a better understanding of their functional allelic variation and distribution. In this study, a consensus map of the common carp was constructed, based on four populations, to compare the distribution and variation of QTLs. The consensus map spans 2371.6 cM across the 42 linkage groups and comprises 257 microsatellites and 421 SNPs, with a mean marker interval of 3.7 cM/marker. Sixty-seven QTLs affecting four growth traits from the four populations were mapped to the consensus map. Only one QTL was common to three populations, and nine QTLs were detected in two populations. However, no QTL was common to all four populations. The results of the QTL comparison suggest that the QTLs are responsible for the phenotypic variability observed for these traits in a broad array of common carp germplasms. The study also reveals the different genetic performances between major and minor genes in different populations.

Genetic positioning of centromeres through half-tetrad analysis in gynogenetic diploid families of the Zhikong scallop (Chlamys farreri)

Centromere mapping is a powerful tool for improving linkage maps, investigating crossover events, and understanding chiasma interference during meiosis. Ninety microsatellite markers selected across all linkage groups (LGs) from a previous Chlamys farreri genetic map were studied in three artificially induced meiogynogenetic families for centromere mapping by half-tetrad analysis. Inheritance analyses showed that all 90 microsatellite loci conformed to Mendelian inheritance in the control crosses, while 4.4 % of the microsatellite loci showed segregation departures from an expected 1:1 ratio of two homozygote classes in meiogynogenetic progeny. The second division segregation frequency (y) of the microsatellites ranged from 0.033 to 0.778 with a mean of 0.332, confirming the occurrence of partial chiasma interference in this species. Heterogeneity of y is observed in one of 42 cases in which markers were typed in more than one family, suggesting variation in gene-centromere recombination among families. Centromere location was mostly in accordance with the C. farreri karyotype, but differences in marker order between linkage and centromere maps occurred. Overall, this study makes the genetic linkage map a more complete and informative tool for genomic studies and it will also facilitate future research of the structure and function of the scallop centromeres.

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

Two novel SNPs of the 3-hydroxy-3-methylglutaryl coenzyme A reductase gene associated with growth and meat quality traits in the chicken

The enzyme 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGCR) is rate-limiting for metabolism of cholesterol; it plays an important role in endogenous cholesterol biosynthesis. We used DNA sequencing technology and created restriction site PCR-RFLP to detect HMGCR SNPs in an F(2) resource population of Gushi chicken and Anka broilers. We found a G/T mutation (Gln/His) in exon 17 and a T/C mutation (Pro/Pro) in exon 18. Based on association analysis of these HMGCR polymorphisms in 864 Gushi/Anka F(2) hybrids, these two mutations have significant effects on growth, carcass, meat quality, and lipid concentration.

High throughput mining and characterization of microsatellites from common carp genome

In order to supply sufficient microsatellite loci for high-density linkage mapping, whole genome shotgun (WGS) sequences of the common carp (Cyprinus carpio) were assembled and surveyed for microsatellite identification. A total of 79,014 microsatellites were collected which were harbored in 68,827 distinct contig sequences. These microsatellites were characterized in the common carp genome. Information of all microsatellites, including previously published BAC-based microsatellites, was then stored in a MySQL database, and a web-based database interface (http://genomics.cafs.ac.cn/ssrdb) was built for public access and download. A total of 3,110 microsatellites, including 1,845 from WGS and 1,265 from BAC end sequences (BES), were tested and genotyped on a mapping family with 192 individuals. A total of 963 microsatellites markers were validated with polymorphism in the mapping family. They will soon be used for high-density linkage mapping with a vast number of polymorphic SNP markers.

High genetic diversity and low differentiation in mud crab (Scylla paramamosain) along the southeastern coast of China revealed by microsatellite markers

The mud crab (Scylla paramamosain) is a carnivorous portunid crab, mainly distributed along the southeastern coast of China. Mitochondrial DNA analysis in a previous study indicated a high level of genetic diversity and a low level of genetic differentiation. In this study, population genetic diversity and differentiation of S. paramamosain were investigated using nine microsatellite markers. In total, 397 wild specimens from 11 locations on the southeastern coast of China were sampled and genotyped. A high level of genetic diversity was observed, with the number of alleles, and the observed and expected heterozygosity per location in the range 7.8-9.6, 0.62-0.77 and 0.66-0.76, respectively. AMOVA analysis indicated a low level of genetic differentiation among the 11 locations, despite the fact that a statistically significant fixation index (F(ST)) value was found (F(ST)=0.0183, P<0.05). Out of 55 pairwise location comparisons, 39 showed significant F(ST) values (P<0.05), but all of them were lower than 0.05, except for one between Sanmen and Shantou locations. No significant deficiency of heterozygotes (inbreeding coefficient F(IS)=0.0007, P>0.05) was detected for all locations except Sanmen and Zhanjiang. Cluster analysis using UPGMA showed that all locations fell into one group except Sanmen. Significant association was found between genetic differentiation in terms of F(ST)/(1-F(ST)) and the natural logarithm of geographical distance (r(2)=0.1139, P=0.02), indicating that the genetic variation pattern closely resembled an isolation by distance model. This study supports the proposal of high genetic diversity and low genetic differentiation in S. paramamosain along the southeastern coast of China.

Characterization of common carp transcriptome: sequencing, de novo assembly, annotation and comparative genomics

BACKGROUND

Common carp (Cyprinus carpio) is one of the most important aquaculture species of Cyprinidae with an annual global production of 3.4 million tons, accounting for nearly 14% of the freshwater aquaculture production in the world. Due to the economical and ecological importance of common carp, genomic data are eagerly needed for genetic improvement purpose. However, there is still no sufficient transcriptome data available. The objective of the project is to sequence transcriptome deeply and provide well-assembled transcriptome sequences to common carp research community.

RESULT

Transcriptome sequencing of common carp was performed using Roche 454 platform. A total of 1,418,591 clean ESTs were collected and assembled into 36,811 cDNA contigs, with average length of 888 bp and N50 length of 1,002 bp. Annotation was performed and a total of 19,165 unique proteins were identified from assembled contigs. Gene ontology and KEGG analysis were performed and classified all contigs into functional categories for understanding gene functions and regulation pathways. Open Reading Frames (ORFs) were detected from 29,869 (81.1%) contigs with an average ORF length of 763 bp. From these contigs, 9,625 full-length cDNAs were identified with sequence length from 201 bp to 9,956 bp. Comparative analysis revealed that 27,693(75.2%) contigs have significant similarity to zebrafish Refseq proteins, and 24,371(66.2%), 24,501(66.5%) and 25,025(70.0%) to teraodon, medaka and three-spined stickleback refseq proteins. A total of 2,064 microsatellites were initially identified from 1,730 contigs, and 1,639 unique sequences had sufficient flanking sequences on both sides for primer design.

CONCLUSION

The transcriptome of common carp had been deep sequenced, de novo assembled and characterized, providing the valuable resource for better understanding of common carp genome. The transcriptome data will facilitate future functional studies on common carp genome, and gradually apply in breeding programs of common carp, as well as closely related other Cyprinids.

Generation of the first BAC-based physical map of the common carp genome

Background

Common carp (Cyprinus carpio), a member of Cyprinidae, is the third most important aquaculture species in the world with an annual global production of 3.4 million metric tons, accounting for nearly 14% of the all freshwater aquaculture production in the world. Apparently genomic resources are needed for this species in order to study its performance and production traits. In spite of much progress, no physical maps have been available for common carp. The objective of this project was to generate a BAC-based physical map using fluorescent restriction fingerprinting.

Result

The first generation of common carp physical map was constructed using four- color High Information Content Fingerprinting (HICF). A total of 72,158 BAC clones were analyzed that generated 67,493 valid fingerprints (5.5 × genome coverage). These BAC clones were assembled into 3,696 contigs with the average length of 476 kb and a N50 length of 688 kb, representing approximately 1.76 Gb of the common carp genome. The largest contig contained 171 BAC clones with the physical length of 3.12 Mb. There are 761 contigs longer than the N50, and these contigs should be the most useful resource for future integrations with linkage map and whole genome sequence assembly. The common carp physical map is available at http://genomics.cafs.ac.cn/fpc/WebAGCoL/Carp/WebFPC/.

Conclusion

The reported common carp physical map is the first physical map of the common carp genome. It should be a valuable genome resource facilitating whole genome sequence assembly and characterization of position-based genes important for aquaculture traits.

A genetic linkage map and comparative genome analysis of common carp (Cyprinus carpio L.) using microsatellites and SNPs

A genetic linkage map is a powerful research tool for mapping traits of interest and is essential to understanding genome evolution. The aim of this study is to provide an expanded genetic linkage map of common carp to effectively carry out quantitative trait loci analysis and conduct comparative mapping analysis between lineages. Here, we constructed a genetic linkage map of common carp (Cyprinus carpio L.) using microsatellite and single-nucleotide polymorphism (SNP) markers in a 159 sibling family. A total of 246 microsatellites and 306 SNP polymorphic markers were genotyped in this family. Linkage analysis using JoinMap 4.0 organized 427 markers (186 microsatellites and 241 SNPs) to 50 linkage groups, ranging in size from 1.4 to 130.1 cM. Each group contained 2-30 markers. The linkage map covered a genetic distance of 2,039.2 cM and the average interval for markers within the linkage groups was approximately 6.4 cM. In addition, comparative genome analysis within five model teleost fish revealed a high percentage (74.7%) of conserved loci corresponding to zebrafish chromosomes. In most cases, each zebrafish chromosome comprised two common carp linkage groups. The comparative analysis also revealed independent chromosome rearrangements in common carp and zebrafish. The linkage map will be of great assistance in mapping genes of interest and serve as a reference to approach comparative mapping and enable further insights into the comprehensive investigations of genome evolution of common carp.

Molecular characterization of CART, AgRP, and MC4R genes and their expression with fasting and re-feeding in common carp (Cyprinus carpio)

Cocaine and amphetamine-regulated transcript (CART), agouti-related proteins (AgRP) and Melanocortin 4 Receptor (MC4R) involves in the control of appetite. The genes were cloned and characterized, and their regulation was studied in common carp. The CARTI and CARTII genes encode 117- and 120-amino acids, respectively. The AgRP-1 and AgRP-2 genes encode 128- and 136-amino acids, respectively. CARTI was principally expressed in the brain, eye and ovary, while CARTII was highly expressed in the brain. AgRP-1 was strongly expressed in the brain, intestine, testis and eye, while AgRP-2 was highly expressed only in the gill and eye. The MC4R gene, encoding 326-amino acids, was mainly expressed in the brain testis, pituitary and eye. Phylogenetic analysis had been conducted which implied that both CARTI/CARTII and AgRP-1/AgRP-2 might derived from gene duplication events during genome evolution of common carp. CART, AgRP and MC4R gene expression in brain were decreased after fasting treatment and increased sharply after refeeding comparing with normal fed controls, which suggested that CART, AgRP and MC4R are involved in appetite regulation in common carp.