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

Advances in biotechnology and breeding innovations in China's marine aquaculture

Biotechnology is the key driving force behind the sustainable development of aquaculture, as biological innovation would significantly improve the capabilities of aquatic breeding and achieve independent and controllable seeding sources to ensure food safety. In this article, we have analyzed the current status and existing problems of marine aquaculture in China. Based on these data, we have summarized the recent (especially the last 10 years) biotechnological innovation and breeding progress of marine aquaculture in China, including whole genome sequencing, sex-related marker screening, genomic selection, and genome editing, as well as progress of improved marine fish varieties in China. Finally, the perspectives in this field have been discussed, and three future countermeasures have been proposed.

Unveiling Gene Expression Dynamics during Early Embryogenesis in Cynoglossus semilaevis: A Transcriptomic Perspective

Commencing with sperm-egg fusion, the early stages of metazoan development include the cleavage and formation of blastula and gastrula. These early embryonic events play a crucial role in ontogeny and are accompanied by a dramatic remodeling of the gene network, particularly encompassing the maternal-to-zygotic transition. Nonetheless, the gene expression dynamics governing early embryogenesis remain unclear in most metazoan lineages. We conducted transcriptomic profiling on two types of gametes (oocytes and sperms) and early embryos (ranging from the four-cell to the gastrula stage) of an economically valuable flatfish-the Chinese tongue sole Cynoglossus semilaevis (Pleuronectiformes: Cynoglossidae). Comparative transcriptome analysis revealed that large-scale zygotic genome activation (ZGA) occurs in the blastula stage, aligning with previous findings in zebrafish. Through the comparison of the most abundant transcripts identified in each sample and the functional analysis of co-expression modules, we unveiled distinct functional enrichments across different gametes/developmental stages: actin- and immune-related functions in sperms; mitosis, transcription inhibition, and mitochondrial function in oocytes and in pre-ZGA embryos (four- to 1000-cell stage); and organ development in post-ZGA embryos (blastula and gastrula). These results provide insights into the intricate transcriptional regulation of early embryonic development in Cynoglossidae fish and expand our knowledge of developmental constraints in vertebrates.

De novo assembly and characterization of a highly degenerated ZW sex chromosome in the fish Megaleporinus macrocephalus

BACKGROUND

Megaleporinus macrocephalus (piauçu) is a Neotropical fish within Characoidei that presents a well-established heteromorphic ZZ/ZW sex determination system and thus constitutes a good model for studying W and Z chromosomes in fishes. We used PacBio reads and Hi-C to assemble a chromosome-level reference genome for M. macrocephalus. We generated family segregation information to construct a genetic map, pool sequencing of males and females to characterize its sex system, and RNA sequencing to highlight candidate genes of M. macrocephalus sex determination.

RESULTS

The reference genome of M. macrocephalus is 1,282,030,339 bp in length and has a contig and scaffold N50 of 5.0 Mb and 45.03 Mb, respectively. In the sex chromosome, based on patterns of recombination suppression, coverage, FST, and sex-specific SNPs, we distinguished a putative W-specific region that is highly differentiated, a region where Z and W still share some similarities and is undergoing degeneration, and the PAR. The sex chromosome gene repertoire includes genes from the TGF-β family (amhr2, bmp7) and the Wnt/β-catenin pathway (wnt4, wnt7a), some of which are differentially expressed.

CONCLUSIONS

The chromosome-level genome of piauçu exhibits high quality, establishing a valuable resource for advancing research within the group. Our discoveries offer insights into the evolutionary dynamics of Z and W sex chromosomes in fish, emphasizing ongoing degenerative processes and indicating complex interactions between Z and W sequences in specific genomic regions. Notably, amhr2 and bmp7 are potential candidate genes for sex determination in M. macrocephalus.

Microsatellite Genome-Wide Database Development for the Commercial Blackhead Seabream (Acanthopagrus schlegelii)

Simple sequence repeats (SSRs), the markers with the highest polymorphism and co-dominance degrees, offer a crucial genetic research resource. Limited SSR markers in blackhead seabream have been reported. The availability of the blackhead seabream genome assembly provided the opportunity to carry out genome-wide identification for all microsatellite markers, and bioinformatic analyses open the way for developing a microsatellite genome-wide database in blackhead seabream. In this study, a total of 412,381 SSRs were identified in the 688.08 Mb genome by Krait software. Whole-genome sequences (10×) of 42 samples were aligned against the reference genome and genotyped using the HipSTR tools by comparing and counting repeat number variation across the SSR loci. A total of 156,086 SSRs with a 2–4 bp repeat were genotyped by HipSTR tools, which accounted for 55.78% of the 2–4 bp SSRs in the reference genome. High accuracy of genotyping was observed by comparing HipSTR tools and PCR amplification. A set of 109,131 loci with a number of alleles ≥ 3 and with a number of genotyped individuals ≥ 6 were reserved to constitute the polymorphic SSR database. Fifty-one polymorphic SSR loci were identified through PCR amplification. This strategy to develop polymorphic SSR markers not only obtained a large set of polymorphic SSRs but also eliminated the need for laborious experimental screening. SSR markers developed in this study may facilitate blackhead seabream research, which lays a certain foundation for further gene tagging and genetic linkage analysis, such as marker-assisted selection, genetic mapping, as well as comparative genomic analysis.

Germplasm Resources and Strategy for Genetic Breeding of Lycium Species: A Review

Lycium species (goji), belonging to Solanaceae, are widely spread in the arid to semiarid environments of Eurasia, Africa, North and South America, among which most species have affinal drug and diet functions, resulting in their potential to be a superior healthy food. However, compared with other crop species, scientific research on breeding Lycium species lags behind. This review systematically introduces the present germplasm resources, cytological examination and molecular-assisted breeding progress in Lycium species. Introduction of the distribution of Lycium species around the world could facilitate germplasm collection for breeding. Karyotypes of different species could provide a feasibility analysis of fertility between species. The introduction of mapping technology has discussed strategies for quantitative trait locus (QTL) mapping in Lycium species according to different kinds of traits. Moreover, to extend the number of traits and standardize the protocols of trait detection, we also provide 1,145 potential traits (275 agronomic and 870 metabolic) in different organs based on different reference studies on Lycium, tomato and other Solanaceae species. Finally, perspectives on goji breeding research are discussed and concluded. This review will provide breeders with new insights into breeding Lycium 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.

Chromosome anchoring in Senegalese sole (Solea senegalensis) reveals sex-associated markers and genome rearrangements in flatfish

The integration of physical and high-density genetic maps is a very useful approach to achieve chromosome-level genome assemblies. Here, the genome of a male Senegalese sole (Solea senegalensis) was de novo assembled and the contigs were anchored to a high-quality genetic map for chromosome-level scaffolding. Hybrid assembled genome was 609.3 Mb long and contained 3403 contigs with a N50 of 513 kb. The linkage map was constructed using 16,287 informative SNPs derived from ddRAD sequencing in 327 sole individuals from five families. Markers were assigned to 21 linkage groups with an average number of 21.9 markers per megabase. The anchoring of the physical to the genetic map positioned 1563 contigs into 21 pseudo-chromosomes covering 548.6 Mb. Comparison of genetic and physical distances indicated that the average genome-wide recombination rate was 0.23 cM/Mb and the female-to-male ratio 1.49 (female map length: 2,698.4 cM, male: 2,036.6 cM). Genomic recombination landscapes were different between sexes with crossovers mainly concentrated toward the telomeres in males while they were more uniformly distributed in females. A GWAS analysis using seven families identified 30 significant sex-associated SNP markers located in linkage group 18. The follicle-stimulating hormone receptor appeared as the most promising locus associated with sex within a region with very low recombination rates. An incomplete penetrance of sex markers with males as the heterogametic sex was determined. An interspecific comparison with other Pleuronectiformes genomes identified a high sequence similarity between homologous chromosomes, and several chromosomal rearrangements including a lineage-specific Robertsonian fusion in S. senegalensis.

A Chromosome-Level Genome Assembly of the Mandarin Fish (Siniperca chuatsi)

The mandarin fish, Siniperca chuatsi, is an economically important perciform species with widespread aquaculture practices in China. Its special feeding habit, acceptance of only live prey fishes, contributes to its delicious meat. However, little is currently known about related genetic mechanisms. Here, we performed whole-genome sequencing and assembled a 758.78 Mb genome assembly of the mandarin fish, with the scaffold and contig N50 values reaching 2.64 Mb and 46.11 Kb, respectively. Approximately 92.8% of the scaffolds were ordered onto 24 chromosomes (Chrs) with the assistance of a previously established genetic linkage map. The chromosome-level genome contained 19,904 protein-coding genes, of which 19,059 (95.75%) genes were functionally annotated. The special feeding behavior of mandarin fish could be attributable to the interaction of a variety of sense organs (such as vision, smell, and endocrine organs). Through comparative genomics analysis, some interesting results were found. For example, olfactory receptor (OR) genes (especially the beta and delta types) underwent a significant expansion, and endocrinology/vision related npy, spexin, and opsin genes presented various functional mutations. These may contribute to the special feeding habit of the mandarin fish by strengthening the olfactory and visual systems. Meanwhile, previously identified sex-related genes and quantitative trait locis (QTLs) were localized on the Chr14 and Chr17, respectively. 155 toxin proteins were predicted from mandarin fish genome. In summary, the high-quality genome assembly of the mandarin fish provides novel insights into the feeding habit of live prey and offers a valuable genetic resource for the quality improvement of this freshwater fish.

A simple and rapid method for fish sex identification based on recombinase-aided amplification and its use in Cynoglossus semilaevis

Fish sex identification is a basic technique of great importance for both fish genetic studies and fisheries. Due to the sexual reversal phenomenon in many fish species, a simple and rapid molecular identification method for fish genetic sex is urgently needed to suit versatile detection scenarios, such as point-of-need applications. In this study, we took Cynoglossus semilaevis as an example, established a recombinase-aided amplification (RAA)-based method for sex identification, and combined the RAA-detection with two result visualization approaches with distinct features, capillary electrophoresis (CE) and lateral flow dipstick (LFD). Specific primers and probe were designed to specifically detect the sex chromosome W of C. semilaevis in order to distinguish the genetic sex between males, pseudo-males and females. To evaluate the performance of our methods, the genetic sex for twenty-eight males, sixty-eight pseudo-males and fifty-four females were examined with the RAA-based method and classical PCR-based genotyping method, demonstrating the consistent results of sex identification between both methods. The RAA-LFD method is operationally simple, rapid (~ 30 min) and holds great potential for point-of-need applications of fish sex identification, including fishery fields. The method presented here could be effective for identifying fish gender with the ZW karyotype.

A high-density genetic map construction and sex-related loci identification in Chinese Giant salamander

BACKGROUND

The Chinese giant salamander Andrias davidianus is an important amphibian species in China because of its increasing economic value, protection status and special evolutionary position from aquatic to terrestrial animal. Its large genome presents challenges to genetic research. Genetic linkage mapping is an important tool for genome assembly and determination of phenotype-related loci.

RESULTS

In this study, we constructed a high-density genetic linkage map using ddRAD sequencing technology to obtain SNP genotyping data of members from an full-sib family which sex had been determined. A total of 10,896 markers were grouped and oriented into 30 linkage groups, representing 30 chromosomes of A. davidianus. The genetic length of LGs ranged from 17.61 cM (LG30) to 280.81 cM (LG1), with a mean inter-locus distance ranging from 0.11(LG3) to 0.48 cM (LG26). The total genetic map length was 2643.10 cM with an average inter-locus distance of 0.24 cM. Three sex-related loci and four sex-related markers were found on LG6 and LG23, respectively.

CONCLUSION

We constructed the first High-density genetic linkage map and identified three sex-related loci in the Chinese giant salamander. Current results are expected to be a useful tool for future genomic studies aiming at the marker-assisted breeding of the species.

Transcriptome Profiling Reveals the Sexual Dimorphism of Gene Expression Patterns during Gonad Differentiation in the Half-Smooth Tongue Sole (Cynoglossus semilaevis)

The half-smooth tongue sole (Cynoglossus semilaevis), one of the most economically-important fish species in China, exhibits sexually dimorphic growth. An understanding of sex-related gene expression patterns in the tongue sole may inform sex regulation and breeding processes that increase fish production. However, the gene expression patterns during gonad development in the tongue sole remain unknown. In this study, transcriptome sequencing analyses were performed during gonad differentiation in the tongue sole, namely, at 62 days post-hatching (dph), 100 dph, 120 dph, and 150 dph. Differentially expressed genes associated with sex differentiation and gonad development were identified at each time point. Trend analysis showed that gene expression patterns varied over time. These expression patterns either explained common, non-sexually-dimorphic features or indicated significant sexual dimorphism. Transcript structure analyses identified both sex and time differences among samples. This study investigated the time-dependent expression patterns of several sex-related genes, including Dmrt1, Amh, Foxl2, aromatase encoding gene, Esr, and the Sox gene family, during gonad differentiation in the tongue sole. These results might clarify the significant sexual differences during early development in the tongue sole and might provide insight into the mechanisms controlling sex differentiation and development.

Cloning, expression and functional analysis of the desert hedgehog (dhh) gene in Chinese tongue sole (Cynoglossus semilaevis)

Desert hedgehog (dhh) is a gene that is crucial for spermatogenesis and Leydig cell differentiation, but little is known regarding its influence on gonadal differentiation and development in fish. To understand its function, we cloned and characterized the dhh gene from Cynoglossus semilaevis (csdhh). The full length csdhh cDNA was 2473 bp, including a 1386 bp open reading frame (ORF), a 475 bp 5'-UTR, and a 612 bp 3'-UTR, encoding a predicted protein of 461 amino acid residues. Phylogenetic analysis showed that the putative protein belongs to the hedgehog (HH) family, and contains typical HH-N and HH-C domains. Amino acid sequence analysis revealed that CsDhh shares many features with Dhh analogues in other teleost species. Real-time quantitative PCR showed that csdhh was detected in eight different tissues in male and female tongue sole. During early embryonic development, the relative expression of the csdhh was significantly higher in the neural stage than in other embryonic developmental stages (P < 0.05). csdhh was detected at 20 days after hatching (dah) and at the critical period of male gonadal differentiation (80-95 dah), the relative expression of the csdhh was significantly higher in the male gonads than the female gonads. In 5, 8, and 12 month old gonads, the relative expression of the csdhh was significantly higher in male and pseudo-male than in female fish. The in situ hybridization (ISH) results showed that the hybridization signal was strongly expressed in primary and secondary spermatocytes, spermatids, and sertoli cells of the 1-year-old fish testis, with only weak signal expression in the corresponding ovarian tissue. These results suggest that csdhh is highly conserved in evolution and plays an important role in spermatogenesis in males and pseudo-males.

A first genetic linage map construction and QTL mapping for growth traits in Larimichthys polyactis

Larimichthys polyactis is a commercially important marine fish species in Eastern Asia, yet very few genetic resources exist. In particular, genetic linkage maps are critical tools for genetic breeding. In this study, we generated a high resolution linkage map from a family of 110 individuals and their parents by resequencing the individuals. 3,802 effective SNPs were mapped to 24 linkage groups (LGs). The map spanned 2,567.39 cm, with an average marker interval of 0.66 cm. We used the map to conduct QTL analysis for growth traits, and found 31 markers were significantly associated with growth-related traits. Specifically, three SNPs were identified for total length, nineteen SNPs for body length, and nine SNPs for body weight. The identified SNPs could explain 15.2-22.6% of the phenotypic variation. SNPs associated with growth traits were distributed on LG6 and LG11, and candidate genes included, kif26b, bat1, gna1, gbgt1, and amfr, which may regulate growth. The linkage map and mapped QTLs would be useful for improving the quality of L. polyactis via marker-assisted selection.

QTL Fine Mapping for Sex Determination Region in Bighead Carp (Hypophthalmichthys nobilis) and Comparison with Silver Carp (Hypophthalmichthys molitrix)

Bighead carp (Hypophthalmichthys nobilis) and silver carp (Hypophthalmichthys molitrix) are genetically close aquaculture fish in the Cyprinidae, which have been confirmed to hold XX/XY sex determination. However, genomic locations of potential sex-related loci in these two fishes are still unknown. In this study, a high-resolution genetic linkage map was constructed by using 2976 SNP and 924 microsatellite markers in a F1 full-sib family of bighead carp, the length of which spanned 2022.34 cM with an average inter-marker distance of 0.52 cM. Comparative genomics revealed a high level of genomic synteny between bighead carp and zebrafish as well as grass carp. QTL fine mapping for sex trait was performed based on this linkage map of bighead carp and an unpublished linkage map of silver carp. A map distance of 3.863 cM (69.787-73.650 cM) on LG19 of bighead carp and 4.705 cM (79.096-83.801 cM) on LG21 of silver carp was significantly associated with sex phenotypes, and these two LGs are homologous between two fish species. Fourteen markers harboring in these regions were in strong linkage disequilibrium with the sex phenotype variance explained (PVE) varying from 89 to 100%. Two common markers were mapped on the QTL regions of bighead carp and silver carp, suggesting that these two carp species may have similar genetic bases for sex determination. Eleven potentially sex-related genes were identified within or near the sex QTL markers in two species. This study provided insights into elucidating mechanisms and evolution of sex determination in cyprinid fishes.

Construction of High-Density Genetic Map and Mapping of Sex-Related Loci in the Yellow Catfish (Pelteobagrus fulvidraco)

The yellow catfish (Pelteobagrus fulvidraco) is a very important aquaculture species distributed in freshwater area of China. All-male yellow catfish is very popular in aquaculture because of their significant sex dimorphism phenomena. The males grow much faster than females in full-sibling family. However, the sex dimorphism mechanism is still unclear in yellow catfish. In order to better understand the genetic basis of yellow catfish sexual dimorphism, it is vital to map the sex-related traits and localize the candidate genes across yellow catfish whole genome. Here, we constructed a high-density linkage map of yellow catfish using genotyping-by-sequencing (GBS) strategy. A total of 5705 single-nucleotide polymorphism (SNP) markers were mapped to 26 different linkage groups (LGs) using 184 F1 offspring. The total genetic map length was 3071.59 cM, with an average interlocus distance of 0.54 cM. Eleven significant sex-related QTLs in yellow catfish were identified. Six sex-related genes were identified from the region of reference genome near these QTLs including amh, gnrhr, vasa, lnnr1, foxl2, and bmp15. The high-density genetic linkage map provides valuable resources for yellow catfish molecular assistant breeding and elucidating sex differentiation process. Moreover, the comparative genomic study was analyzed among yellow catfish, channel catfish, and zebrafish. It revealed highly conserved chromosomal distribution between yellow catfish and channel catfish.

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.

Genome survey and high-resolution backcross genetic linkage map construction of the ridgetail white prawn Exopalaemon carinicauda applications to QTL mapping of growth traits

BACKGROUND

High-resolution genetic linkage map is critical for QTL mapping, genome sequence assembly and marker-assisted selection in aquaculture species. The ridgetail white prawn Exopalaemon carinicauda is one of the most economic shrimp species naturally distributed in the coasts of eastern China and western Korea. However, quite limited genomics and genetics information have been exploited for genetic improvement of economic traits in this species.

RESULTS

In the present study, we conducted genome survey and constructed high-resolution genetic linkage maps of the ridgetail white prawn with reciprocal-cross mapping family genotyped using next-generation sequencing approaches. The estimated genome size was 9.33 Gb with a heterozygosity of 0.26% and a repeat sequence ratio of 76.62%. 65,772 protein-coding genes were identified by genome annotation. A total of 10,384 SNPs were used to high-throughput genotyping and assigned to 45 linkage groups (LGs) from reciprocal backcross families of E. carinicauda, and the average marker distances were 0.73 cM and 0.55 cM, respectively. Based on the high-resolution linkage map, twenty-three QTLs related to five growth traits were detected. All QTLs could explain 8.8-15.7% of the total growth-traits variation.

CONCLUSIONS

The genome size of E. carinicauda was estimated more accurately by genome survey analysis, which revealed basic genomic architecture. The first high-resolution backcross genetic linkage map and QTLs related to growth traits will provide important information for QTL fine mapping, genome assembly and genetic improvement of E. carinicauda and other palaemon shrimps.

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.

A high-density genetic linkage map and QTL mapping for growth and sex of yellow drum (Nibea albiflora)

A high-density genetic linkage map is essential for the studies of comparative genomics and gene mapping, and can facilitate assembly of reference genome. Herein, we constructed a high-density genetic linkage map with 8,094 SNPs selected from 113 sequenced fish of a F1 family. Ultimately, the consensus map spanned 3818.24 cM and covered nearly the whole genome (99.4%) with a resolution of 0.47 cM. 1,457 scaffolds spanning 435.15 Mb were anchored onto 24 linkage groups, accounting for 80.7% of the draft genome assembly of the yellow drum. Comparative genomic analyses with medaka and zebrafish genomes showed superb chromosome-scale synteny between yellow drum and medaka. QTL mapping and association analysis congruously revealed 22 QTLs for growth-related traits and 13 QTLs for sex dimorphism. Some important candidate genes such as PLA2G4A, BRINP3 and P2RY1 were identified from these growth-related QTL regions. A gene family including DMRT1, DMRT2 and DMRT3 was identified from these sex-related QTL regions on the linkage group LG9. We demonstrate that this linkage map can facilitate the ongoing marker-assisted selection and genomic and genetic studies for yellow drum.

Genomic Selection Using BayesCπ and GBLUP for Resistance Against Edwardsiella tarda in Japanese Flounder (Paralichthys olivaceus)

The Japanese flounder is one of the most widely farmed economic flatfish species throughout eastern Asia including China, Korea, and Japan. Edwardsiella tarda is a major species of pathogenic bacteria that causes ascites disease and, consequently, a huge economy loss for Japanese flounder farming. After generation selection, traditional breeding methods can hardly improve the E. tarda resistance effectively. Genomic selection is an effective way to predict the breeding potential of parents and has rarely been used in aquatic breeding. In this study, we chose 931 individuals from 90 families, challenged by E. tarda from 2013 to 2015 as a reference population and 71 parents of these families as selection candidates. 1,934,475 markers were detected via genome sequencing and applied in this study. Two different methods, BayesCπ and GBLUP, were used for genomic prediction. In the reference population, two methods led to the same accuracy (0.946) and Pearson's correlation results between phenotype and genomic estimated breeding value (GEBV) of BayesCπ and GBLUP were 0.912 and 0.761, respectively. In selection candidates, GEBVs from two methods were highly similar (0.980). A comparison of GEBV with the survival rate of families that were structured by selection candidates showed correlations of 0.662 and 0.665, respectively. This study established a genomic selection method for the Japanese flounder and for the first time applied this to E. tarda resistance breeding.

Gonadal transcriptome sequencing of the critically endangered Acipenser dabryanus to discover candidate sex-related genes

Background

Acipenser dabryanus, an endemic Chinese species, has been listed as a first-class protected animal in China. Sturgeons are among the oldest and most primitive group of existing fish in the world and occupy a special place in the evolutionary history of fish. Thus, a study of the reproduction and sex differentiation of sturgeon will be of great value for fish as well as the whole vertebrate group.

Methods

In this study, we conducted comparative analysis of the testes and ovaries transcriptomes of A. dabryanus to screen for sex-differentiation and sexual development-related genes.

Results

The transcriptome sequencing of six cDNA libraries generated 265 million clean reads, encompassing 79 Gb of sequences. The N50 and mean length of the identified 91,375 unigenes were 1,718 and 989 bp, respectively. A total of 6,306, 9,961, 13,170, 15,484, and 23,588 unigenes were annotated in the clusters of orthologous groups, gene ontology categories, Kyoto Encyclopedia of Genes and Genomes Pathway, euKaryotic orthologous groups, and NCBI non-redundant protein databases, respectively. A total of 5,396 differentially expressed genes were found between the two sexes, with 1,938 predicted to be up-regulated in ovaries and 3,458 in testes. A total of 73 candidate genes known to be involved in sex differentiation and sexual development were searched in the transcriptome of A. dabryanus of which 52 showed significant similarity. We highlighted six genes that are differentially expressed between the two sexes and may play important roles in sex differentiation and gonad maintenance. In addition, 24,271 simple sequence repeats (SSRs) and 550,519 single-nucleotide polymorphisms (SNPs) were detected.

Discussion

This work represents the first transcriptome study comparing the ovary and testis in A. dabryanus. The putative differentially expressed genes between the gonads provide an important source of information for further study of the sex-differentiation related genes and the sex-differentiation mechanism in sturgeons. The SSRs or SNPs identified in this study will be helpful in the discovery of sex-related markers in A. dabryanus.

A high-resolution genetic linkage map and QTL fine mapping for growth-related traits and sex in the Yangtze River common carp (Cyprinus carpio haematopterus)

Background

A high-density genetic linkage map is essential for QTL fine mapping, comparative genome analysis, identification of candidate genes and marker-assisted selection for economic traits in aquaculture species. The Yangtze River common carp (Cyprinus carpio haematopterus) is one of the most important aquacultured strains in China. However, quite limited genetics and genomics resources have been developed for genetic improvement of economic traits in such strain.

Results

A high-resolution genetic linkage map was constructed by using 7820 2b-RAD (2b-restriction site-associated DNA) and 295 microsatellite markers in a F2 family of the Yangtze River common carp (C. c. haematopterus). The length of the map was 4586.56 cM with an average marker interval of 0.57 cM. Comparative genome mapping revealed that a high proportion (70%) of markers with disagreed chromosome location was observed between C. c. haematopterus and another common carp strain (subspecies) C. c. carpio. A clear 2:1 relationship was observed between C. c. haematopterus linkage groups (LGs) and zebrafish (Danio rerio) chromosomes. Based on the genetic map, 21 QTLs for growth-related traits were detected on 12 LGs, and contributed values of phenotypic variance explained (PVE) ranging from 16.3 to 38.6%, with LOD scores ranging from 4.02 to 11.13. A genome-wide significant QTL (LOD = 10.83) and three chromosome-wide significant QTLs (mean LOD = 4.84) for sex were mapped on LG50 and LG24, respectively. A 1.4 cM confidence interval of QTL for all growth-related traits showed conserved synteny with a 2.06 M segment on chromosome 14 of D. rerio. Five potential candidate genes were identified by blast search in this genomic region, including a well-studied multi-functional growth related gene, Apelin.

Conclusions

We mapped a set of suggestive and significant QTLs for growth-related traits and sex based on a high-density genetic linkage map using SNP and microsatellite markers for Yangtze River common carp. Several candidate growth genes were also identified from the QTL regions by comparative mapping. This genetic map would provide a basis for genome assembly and comparative genomics studies, and those QTL-derived candidate genes and genetic markers are useful genomic resources for marker-assisted selection (MAS) of growth-related traits in the Yangtze River common carp.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4613-1) contains supplementary material, which is available to authorized users.

Microsatellites Associated with Growth Performance and Analysis of Resistance to Aeromonas hydrophila in Tambaqui Colossoma macropomum

Tambaqui, Colossoma macropomum, is the main native fish species produced in Brazil, and is an important species for genetic improvement in aquaculture. In addition, breeding studies on this species can be optimized with the use of molecular markers associated with productive phenotypes. The objective of the present study was to test the performance of growth traits and resistance to the bacteria, Aeromonas hydrophila, in association with microsatellite markers in C. macropomum. In this study, three full-sib families were subjected to bacterial challenge and morphometric growth assessments. Tambaqui families subjected to the bacterial challenge differed significantly in death time and mortality rate. There was, however, no association between resistance to bacteria and microsatellite markers. In relation to growth traits, we observed a marker/phenotype association in two microsatellites. The marker in the 6b isoform x5 gene (TNCRC6b) was associated with length, whereas an anonymous marker was associated with height. The present study highlighted the evaluation of molecular markers associated with growth traits, and can serve as the basis for future marker-assisted selection (MAS) of tambaqui.

Characterization of a low-density lipoprotein receptor, Lrp13, in Chinese tongue sole (Cynoglossus semilaevis) and medaka (Oryzias latipes)

As an important economic marine species cultured in China, Chinese tongue sole (Cynoglossus semilaevis) has interested us due to its sexual dimorphism and ZW/ZZ sex determination system. In a previous study, dmrt1 was identified as a dosage-dependent male-determining gene. In the present study, a female-specific expressed gene, cse0440, initially annotated as lrp1b-like, was identified from chromosome W of C. semilaevis. In view of the differences between cse0440 and lrp1b in terms of expression pattern, a phylogenetic analysis containing 85 LRP proteins was constructed and provided an evidence to re-annotate cse0440 as cseLRP13. In addition, two orthologues of cseLRP13 were separately identified from W and Z chromosomes: cseLRP13-W and cseLRP13-Z. The subsequent multiple sequence alignment and syntenic arrangements of LRP13 in C. semilaevis, Japanese medaka (Oryzias latipes), large yellow croaker (Larimichthys crocea), striped bass (Morone saxatilis), white perch (Morone americana) and Fugu rubripes (Takifugu rubripes) further supported this re-annotation. RT-PCR and in situ hybridization revealed that cselrp13 was exclusively expressed in the oocytes and follicles of ovaries. These results suggested that lrp13 may play important roles in female reproduction. In future, with the advancement of micromanipulation in flatfish, the detailed function of two lrp13 orthologues in C. semilaevis will be elucidated.

Locus Mapping, Molecular Cloning, and Expression Analysis of rps6kb2, a Novel Metamorphosis-Related Gene in Chinese Tongue Sole (Cynoglossus semilaevis)

Flatfish metamorphosis denotes the extraordinary transformation of a symmetric pelagic larva into an asymmetric benthic juvenile. This unique process involves eye migration, a 90° rotation in posture, and asymmetrical pigmentation for adaptation to a benthic lifestyle. In the present study, we used genetics to map a metamorphosis-related locus (q-10M) in the male linkage group (LG10M), a small interval of 0.9 cM corresponding to a 1.8 M-bp physical area in chromosome 9 in the Chinese tongue sole (Cynoglossus semilaevis). Combined with single-marker analysis, ribosomal protein S6 kinase 2 (rps6kb2) a member of the family of AGC kinases was identified as a novel metamorphosis-related candidate gene. Its expression pattern during metamorphosis was determined by quantitative RT-PCR and whole-mount in situ hybridization analysis. rps6kb2 gene was significantly expressed in metamorphic climax stage larvae and distributed in all the tissues transforming during metamorphosis, including tail, jaw, eye and skin of larvae. The results suggest that rps6kb2 has a general role in tissue transformations during flatfish metamorphosis including tail changes, skull remodeling, eye migration, and asymmetrical pigmentation.

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

Background

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

Results

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

Conclusion

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

Electronic supplementary material

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

Genome editing reveals dmrt1 as an essential male sex-determining gene in Chinese tongue sole (Cynoglossus semilaevis)

Chinese tongue sole is a marine fish with ZW sex determination. Genome sequencing suggested that the Z-linked dmrt1 is a putative male determination gene, but direct genetic evidence is still lacking. Here we show that TALEN of dmrt1 efficiently induced mutations of this gene. The ZZ dmrt1 mutant fish developed ovary-like testis, and the spermatogenesis was disrupted. The female-related genes foxl2 and cyp19a1a were significantly increased in the gonad of the ZZ dmrt1 mutant. Conversely, the male-related genes Sox9a and Amh were significantly decreased. The dmrt1 deficient ZZ fish grew much faster than ZZ male control. Notably, we obtained an intersex ZW fish with a testis on one side and an ovary on the other side. This fish was chimeric for a dmrt1 mutation in the ovary, and wild-type dmrt1 in the testis. Our data provide the first functional evidence that dmrt1 is a male determining gene in tongue sole.

Single Locus Maintains Large Variation of Sex Reversal in Half-Smooth Tongue Sole (Cynoglossus semilaevis)

Sex determination is a fundamental biological process for individual sex development and population sex ratios. However, for some species, the primary sex might be altered during development, and individuals can develop into the opposite sex. Sex reversal may happen in insects, reptiles, amphibians, and fishes. In half-smooth tongue sole (Cynoglossus semilaevis), some genetically female fish irreversibly reverse to pseudomales, resulting in higher costs in aquaculture owing to a lower growth rate of male fish during a 2-yr growth period. Here, we identified a locus with large controlling effect on sex reversal in the half-smooth tongue sole through genome-wide association study with high-density single nucleotide polymorphisms (SNPs). This SNP is located at the third intron of the F-box and leucine rich repeat protein 17 (FBXL17) gene on the Z chromosome, and it has two alleles, A and T. Genetic females with Z^AW genotypes will never reverse into phenotypic males, but those with Z^TW genotypes can sometimes undergo sex reversal. This SNP explains 82.7% of the genetic variation, or 58.4% of the phenotypic variation. Based on our results, a reproductive management program could be developed to improve the phenotypic female ratio in aquaculture, and elucidate the mechanism of sex reversal in half-smooth tongue sole. We expect that these findings will have a substantial impact on the population management in many harvested species where sex reversal occurs.

A second generation SNP and SSR integrated linkage map and QTL mapping for the Chinese mitten crab Eriocheir sinensis

The Chinese mitten crab Eriocheir sinensis is the most economically important cultivated crab species in China, and its genome has a high number of chromosomes (2n = 146). To obtain sufficient markers for construction of a dense genetic map for this species, we employed the recently developed specific-locus amplified fragment sequencing (SLAF-seq) method for large-scale SNPs screening and genotyping in a F1 full-sib family of 149 individuals. SLAF-seq generated 127,677 polymorphic SNP markers, of which 20,803 valid markers were assigned into five segregation types and were used together with previous SSR markers for linkage map construction. The final integrated genetic map included 17,680 SNP and 629 SSR markers on the 73 linkage groups (LG), and spanned 14,894.9 cM with an average marker interval of 0.81 cM. QTL mapping localized three significant growth-related QTL to a 1.2 cM region in LG53 as well as 146 sex-linked markers in LG48. Genome-wide QTL-association analysis further identified four growth-related QTL genes named LNX2, PAK2, FMRFamide and octopamine receptors. These genes are involved in a variety of different signaling pathways including cell proliferation and growth. The map and SNP markers described here will be a valuable resource for the E. sinensis genome project and selective breeding programs.

Integrating genomic resources of flatfish (Pleuronectiformes) to boost aquaculture production

Flatfish have a high market acceptance thus representing a profitable aquaculture production. The main farmed species is the turbot (Scophthalmus maximus) followed by Japanese flounder (Paralichthys olivaceous) and tongue sole (Cynoglossus semilaevis), but other species like Atlantic halibut (Hippoglossus hippoglossus), Senegalese sole (Solea senegalensis) and common sole (Solea solea) also register an important production and are very promising for farming. Important genomic resources are available for most of these species including whole genome sequencing projects, genetic maps and transcriptomes. In this work, we integrate all available genomic information of these species within a common framework, taking as reference the whole assembled genomes of turbot and tongue sole (>210× coverage). New insights related to the genetic basis of productive traits and new data useful to understand the evolutionary origin and diversification of this group were obtained. Despite a general 1:1 chromosome syntenic relationship between species, the comparison of turbot and tongue sole genomes showed huge intrachromosomic reorganizations. The integration of available mapping information supported specific chromosome fusions along flatfish evolution and facilitated the comparison between species of previously reported genetic associations for productive traits. When comparing transcriptomic resources of the six species, a common set of ~2500 othologues and ~150 common miRNAs were identified, and specific sets of putative missing genes were detected in flatfish transcriptomes, likely reflecting their evolutionary diversification.

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.

Bayesian Analysis for Genetic Architectures of Body Weights and Morphological Traits Using Distorted Markers in Japanese Flounder Paralichthys olivaceus

A doubled haploids (DH) population of 160 individuals is constructed by using mitotic gynogenetics for the study of genetic architectures of the body weight and morphological traits in Japanese flounder. Each DH individual is genotyped for 458 SSR markers, 222 of which segregate distortionally. By modifying conditional probabilities of quantitative trait locus (QTL) genotypes on the distorted flanking markers, Bayesian model selection is used to dissect genetic architectures for the traits. As a result, we identify 42 main-effect QTLs on chromosomes 5, 6, 7, 8, 9, 10, 15, 20, 21, 22, and 59 pairs of interacting QTLs. Among these detected QTLs, the largest interacting QTL is between chromosome 6 and chromosome 9 and accounts for 25.196 % of phenotypic variance for body weights and in a similar trend. Also, many QTLs show pleiotropic effects. The QTL on chromosome 9 simultaneously governs seven traits, BL, BH, FL, HL, PFL, HW, and CW. As compared to method using the uncorrected conditional probabilities of QTL genotypes, our method using corrected conditional probabilities can detect more interacting QTLs for the traits.

Recent Major Advances of Biotechnology and Sustainable Aquaculture in China

Background:

Global aquaculture production has increased continuously over the last five decades, and particularly in China. Its aquaculture has become the fastest growing and most efficient agri-sector, with production accounting for more than 70% of the world’s aquaculture output. In the new century, with serious challenges regarding population, resources and the environment, China has been working to develop high-quality, effective, healthy, and sustainable blue agriculture through the application of modern biotechnology. Sound knowledge related to the biology and ecology of aquatic organisms has laid a solid foundation and provided the innovation and technology for rapid development of the aquaculture industry. Marine biotechnology, which is enabling solutions for ocean productivity and sustainability, has been promoted since the last decades of the 20th Century in China.

Objective:

In this article, priority areas of research, mainly genetic breeding, omics studies, novel production systems, biosecurity, bioprocesses and biorefinery, as well as the major progress of marine biotechnology R&D in China are reviewed.

Conclusion:

Current innovative achievements in China are not enough and the level and frequency of academic advancements must be improved. International cooperation and assistance remain crucial for the success of marine biotechnology.

The First High-Density Genetic Map Construction in Tree Peony (Paeonia Sect. Moutan) using Genotyping by Specific-Locus Amplified Fragment Sequencing

Genetic linkage maps, permitting the elucidation of genome structure, are one of most powerful genomic tools to accelerate marker-assisted breeding. However, due to a lack of sufficient user-friendly molecular markers, no genetic linkage map has been developed for tree peonies (Paeonia Sect. Moutan), a group of important horticultural plants worldwide. Specific-locus amplified fragment sequencing (SLAF-seq) is a recent molecular marker development technology that enable the large-scale discovery and genotyping of sequence-based marker in genome-wide. In this study, we performed SLAF sequencing of an F1 population, derived from the cross P. ostti ‘FenDanBai’ × P. × suffruticosa ‘HongQiao’, to identify sufficient high-quality markers for the construction of high-density genetic linkage map in tree peonies. After SLAF sequencing, a total of 78 Gb sequencing data and 285,403,225 pair-end reads were generated. We detected 309,198 high-quality SLAFs from these data, of which 85,124 (27.5%) were polymorphic. Subsequently, 3518 of the polymorphic markers, which were successfully encoded in to Mendelian segregation types, and were in conformity with the criteria of high-quality markers, were defined as effective markers and used for genetic linkage mapping. Finally, we constructed an integrated genetic map, which comprised 1189 markers on the five linkage groups, and spanned 920.699 centiMorgans (cM) with an average inter-marker distance of 0.774 cM. There were 1115 ‘SNP-only’ markers, 18 ‘InDel-only’ markers, and 56 ‘SNP&InDel’ markers on the map. Among these markers, 450 (37.85%) showed significant segregation distortion (P < 0.05). In conclusion, this investigation reported the first large-scale marker development and high-density linkage map construction for tree peony. The results of this study will serve as a solid foundation not only for marker-assisted breeding, but also for genome sequence assembly for tree peony.

Functional marker detection and analysis on a comprehensive transcriptome of large yellow croaker by next generation sequencing

Large yellow croaker (Larimichthys crocea) is an important economic fish in China and Eastern Asia. Because of the exhaustive fishing and overdense aquaculture, the wild population and the mariculture of the species are facing serious challenges on germplasm degeneration and susceptibility to infectious disease agents. However, a comprehensive transcriptome from multi-tissues of the species has not been reported and functional molecular markers have not yet been detected and analyzed. In this work, we applied RNA-seq with the Illumina Hiseq2000 platform for a multi-tissue sample of large yellow croaker and assembled the transcriptome into 88,103 transcripts. Of them, 52,782 transcripts have been successfully annotated by nt/nr, InterPro, GO and KEGG database. Comparing with public fish proteins, we have found that 34,576 protein coding transcripts are shared in large yellow croaker with zebrafish, medaka, pufferfish, and stickleback. For functional markers, we have discovered 1,276 polymorphic SSRs and 261, 000 SNPs. The functional impact analysis of SNPs showed that the majority (~75%) of small variants cause synonymous mutations in proteins, followed by variations in 3' UTR region. The functional enrichment analysis illuminated that transcripts involved in DNA bindings, enzyme activities, and signal pathways prominently exhibit less single-nucleotide variants but genes for the constituent of the muscular tissue, the cytoskeleton, and the immunity system contain more frequent SNP mutations, which may reflect the structural and functional selections of the translated proteins. This is the first work for the high-throughput detection and analysis of functional polymorphic SSR and SNP markers in a comprehensive transcriptome of large yellow croaker. Our study provides valuable transcript sequence and functional marker resources for the quantitative trait locus (QTL) identification and molecular selection of the species in the research community.

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.

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

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

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

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

Construction of a high-density, high-resolution genetic map and its integration with BAC-based physical map in channel catfish

Construction of genetic linkage map is essential for genetic and genomic studies. Recent advances in sequencing and genotyping technologies made it possible to generate high-density and high-resolution genetic linkage maps, especially for the organisms lacking extensive genomic resources. In the present work, we constructed a high-density and high-resolution genetic map for channel catfish with three large resource families genotyped using the catfish 250K single-nucleotide polymorphism (SNP) array. A total of 54,342 SNPs were placed on the linkage map, which to our knowledge had the highest marker density among aquaculture species. The estimated genetic size was 3,505.4 cM with a resolution of 0.22 cM for sex-averaged genetic map. The sex-specific linkage maps spanned a total of 4,495.1 cM in females and 2,593.7 cM in males, presenting a ratio of 1.7 : 1 between female and male in recombination fraction. After integration with the previously established physical map, over 87% of physical map contigs were anchored to the linkage groups that covered a physical length of 867 Mb, accounting for ∼90% of the catfish genome. The integrated map provides a valuable tool for validating and improving the catfish whole-genome assembly and facilitates fine-scale QTL mapping and positional cloning of genes responsible for economically important traits.

Characterization and development of EST-SSR markers derived from transcriptome of yellow catfish

Yellow catfish (Pelteobagrus fulvidraco) is one of the most important freshwater fish due to its delicious flesh and high nutritional value. However, lack of sufficient simple sequence repeat (SSR) markers has hampered the progress of genetic selection breeding and molecular research for yellow catfish. To this end, we aimed to develop and characterize polymorphic expressed sequence tag (EST)–SSRs from the 454 pyrosequencing transcriptome of yellow catfish. Totally, 82,794 potential EST-SSR markers were identified and distributed in the coding and non-coding regions. Di-nucleotide (53,933) is the most abundant motif type, and AC/GT, AAT/ATT, AAAT/ATTT are respective the most frequent di-, tri-, tetra-nucleotide repeats. We designed primer pairs for all of the identified EST-SSRs and randomly selected 300 of these pairs for further validation. Finally, 263 primer pairs were successfully amplified and 57 primer pairs were found to be consistently polymorphic when four populations of 48 individuals were tested. The number of alleles for the 57 loci ranged from 2 to 17, with an average of 8.23. The observed heterozygosity (H_O), expected heterozygosity (H_E), polymorphism information content (PIC) and fixation index (FIS) values ranged from 0.04 to 1.00, 0.12 to 0.92, 0.12 to 0.91 and −0.83 to 0.93, respectively. These EST-SSR markers generated in this study could greatly facilitate future studies of genetic diversity and molecular breeding in yellow catfish.

Microsatellite-centromere mapping in common carp through half-tetrad analysis in diploid meiogynogenetic families

Gene-centromere (G-C) mapping provides insights into the understanding of the composition, structure, and evolution of vertebrate genomes. Common carp (Cyprinus carpio) is an important aquaculture fish and has been proposed to undertake tetraploidization. In this study, we selected 214 informative microsatellite markers across 50 linkage groups of a common carp genetic map to perform gene-centromere mapping using half-tetrad analysis. A total of 199 microsatellites were segregated under the Mendelian expectations in at least one of the three gynogenetic families and were used for G-C distance estimation. The G-C recombination frequency (y) ranged from 0 to 0.99 (0.43 on average), corresponding to a fixation index (F) of 0.57 after one generation of gynogenesis. Large y values for some loci together with significant correlation between G-C distances and genetic linkage map distances suggested the presence of high interference in common carp. Under the assumption of complete interference, 50 centromeres were localized onto corresponding linkage groups (LGs) of common carp, with G-C distances of centromere-linked markers per LG ranging from 0 to 10.3 cM (2.9 cM on average). Based on the information for centromere positions, we proposed a chromosome formula of 2n = 100 = 58 m/sm + 42 t/st with 158 chromosome arms for common carp, which was similar to a study observed by cytogenetic method. The examination of crossover distributions along 10 LGs revealed that the proportion of crossover chromatids was overall higher than that of non-crossover chromatids in gynogenetic progenies, indicating high recombination levels across most LGs. Comparative genomics analyses suggested that the chromosomes of common carp have undergone extensive rearrangement after genome duplication. This study would be valuable to elucidate the mechanism of genome evolution and integrate physical and genetic maps in common carp.

Construction of a microsatellites-based linkage map for the white grouper (Epinephelus aeneus)

The white grouper (Epinephelus aeneus) is a promising candidate for domestication and aquaculture due to its fast growth, excellent taste, and high market price. A linkage map is an essential framework for mapping quantitative trait loci for economic traits and the study of genome evolution. DNA of a single individual was deep-sequenced, and microsatellite markers were identified in 177 of the largest scaffolds of the sequence assembly. The success rate of developing polymorphic homologous markers was 94.9% compared with 63.1% of heterologous markers from other grouper species. Of the 12 adult mature fish present in the broodstock tank, two males and two females were identified as parents of the assigned offspring by parenthood analysis using 34 heterologous markers. A single full-sib family of 48 individuals was established for the construction of first-generation linkage maps based on genotyping data of 222 microsatellites. The markers were assigned to 24 linkage groups in accordance to the 24 chromosomal pairs. The female and male maps consisting of 203 and 202 markers spanned 1053 and 886 cM, with an average intermarker distance of 5.8 and 5.0 cM, respectively. Mapping of markers to linkage groups ends was enriched by using markers originating from scaffolds harboring telomeric repeat-containing RNA. Comparative mapping showed high synteny relationships among the white grouper, kelp grouper (E. bruneus), orange-spotted grouper (E. coioides), and Nile tilapia (Oreochromis niloticus). Thus, it would be useful to integrate the markers that were developed for different groupers, depending on sharing of sequence data, into a comprehensive consensus map.

A first generation BAC-based physical map of the half-smooth tongue sole (Cynoglossus semilaevis) genome

Background

Half-smooth tongue sole (Cynoglossus semilaevis Günther) has been exploited as a commercially important cultured marine flatfish, and female grows 2–3 times faster than male. Genetic studies, especially on the chromosomal sex-determining system of this species, have been carried out in the last decade. Although the genome of half-smooth tongue sole was relatively small (626.9 Mb), there are still some difficulties in the high-quality assembly of the next generation genome sequencing reads without the assistance of a physical map, especially for the W chromosome of this fish due to abundance of repetitive sequences. The objective of this study is to construct a bacterial artificial chromosome (BAC)-based physical map for half-smooth tongue sole with the method of high information content fingerprinting (HICF).

Results

A physical map of half-smooth tongue sole was constructed with 30, 294 valid fingerprints (7.5 × genome coverage) with a tolerance of 4 and an initial cutoff of 1e-60. A total of 29,709 clones were assembled into 1,485 contigs with an average length of 539 kb and a N50 length of 664 kb. There were 394 contigs longer than the N50 length, and these contigs will be a useful resource for future integration with linkage map and whole genome sequence assembly. The estimated physical length of the assembled contigs was 797 Mb, representing approximately 1.27 coverage of the half-smooth tongue sole genome. The largest contig contained 410 BAC clones with a physical length of 3.48 Mb. Almost all of the 676 BAC clones (99.9%) in the 21 randomly selected contigs were positively validated by PCR assays, thereby confirming the reliability of the assembly.

Conclusions

A first generation BAC-based physical map of half-smooth tongue sole was constructed with high reliability. The map will promote genetic improvement programs of this fish, especially integration of physical and genetic maps, fine-mappings of important gene and/or QTL, comparative and evolutionary genomics studies, as well as whole genome sequence assembly.

Identification of transcriptome-derived microsatellite markers and their association with the growth performance of the mud crab (Scylla paramamosain)

Microsatellite markers from a transcriptome sequence library were initially isolated, and their genetic variation was characterized in a wild population of the mud crab (Scylla paramamosain). We then tested the association between these microsatellite markers and the growth performance of S. paramamosain. A total of 129 polymorphic microsatellite markers were identified, with an observed heterozygosity ranging from 0.19 to 1.00 per locus, an expected heterozygosity ranging from 0.23 to 0.96 per locus, and a polymorphism information content (PIC) ranging from 0.21 to 0.95 per locus. Of these microsatellite markers, 30 showed polymorphism in 96 full-sib individuals of a first generation family. Statistical analysis indicated that three microsatellite markers were significantly associated with 12 growth traits of S. paramamosain. Of these three markers, locus Scpa36 was significantly associated with eight growth traits, namely, carapace length, abdomen width (AW), body height (BH), fixed finger length of the claw, fixed finger width of the claw, fixed finger height of the claw, meropodite length of pereiopod 2, and meropodite length of pereiopod 3 (MLP3) (P<0.05). Locus Scpa75 was significantly associated with five growth traits, namely, internal carapace width, AW, carapace width at spine 8, distance between lateral spine 2 (DLS2), and MLP3 (P<0.05). Locus Spm30 was significantly associated with BH, DLS2, and body weight (P<0.05). Further analysis suggested a set of genotypes (BC at Scpa36, BC and BD at Scpa75, and AC at Spm30) that have great potential in the selection of S. paramamosain for growth traits. These findings will facilitate the development of population conservation genetics and molecular marker-assisted selective breeding of S. paramamosain and other closely related species.

Centromere localization for Bighead Carp (Aristichthys nobilis) through half-tetrad analysis in diploid gynogenetic families

Gene-centromere (G-C) mapping provides insights into structural and behavioural properties of chromosomes. In this study, G-C mapping using microsatellite markers and meiogynogenetic (meiotic gynogenetic) families were performed in bighead carp (Aristichthys nobilis, 2N = 48), which belongs to Cyprinidae. A total of 218 microsatellites were selected across 24 linkage groups (LGs) of a recently well-defined genetic linkage map for bighead carp, with 151 being heterozygous in at least one of six dams in diploid meiogynogenetic families. After tests for Mendelian segregation in two diploid control families, 103 microsatellites were used for G-C distance calculation in 383 gynogens. The second division segregation frequency (y) was computed through half-tetrad analyses, and the values ranged from 0 to 0.97 (mean 0.40). High G-C recombination frequencies (over 0.667) were observed in 18 (17.5%) of the loci examined, which revealed a low level of chiasma interferences compared with other fishes studied previously. Distribution of G-C distances across LGs ranged from 0 cM to 48.5 cM (mean 20 cM) under the assumption of complete interference. All 24 centromeres were localized according to their closest-related microsatellites at 95% confident intervals. The average distance between centromeres and their closest-linked markers was 6.1 cM with 15 out of 24 LGs having a distance below 5 cM. Based on the centromere positions in this study, we proposed a formula of 24 m/sm+24 t/st chromosomes with 92 arms for bighead carp, which was mostly in accordance with a previously reported karyotype for bighead carp (24 m/sm+24 st). These results of centromere localization provide a basic framework and important resources for genetics and comparative genomics studies in bighead carp and its closely-related cyprinid species.