Microsatellite markers from an expressed sequence tag library of half-smooth tongue sole (Cynoglossus semilaevis) and their application in other related fish species

First haploid genetic map based on microsatellite markers in Senegalese sole (Solea senegalensis, Kaup 1858)

The Senegalese sole (Solea senegalensis, Kaup 1858) is a flatfish species of great value for aquaculture. In this study, we develop the first linkage map in this species based on microsatellite markers characterized from genomic DNA libraries and EST databases of Senegalese sole and from other flatfish species. Three reference gynogenetic families were obtained by chromosome-manipulation techniques: two haploid gynogenetics, used to assign and order microsatellites to linkage groups and another diploid gynogenetic family, used for estimating marker-centromere distances. The consensus map consists of 129 microsatellites distributed in 27 linkage groups (LG), with an average density of 4.7 markers per LG and comprising 1,004 centimorgans (cM). Additionally, 15 markers remained unlinked. Through half-tetrad analysis, we were able to estimate the centromere distance for 81 markers belonging to 24 LG, representing an average of 3 markers per LG. Comparative mapping was performed between flatfish species LG and model fish species chromosomes (stickleback, Tetraodon, medaka, fugu and zebrafish). The usefulness of microsatellite markers and the genetic map as tools for comparative mapping and evolution studies is discussed.

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.

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

High-density genetic linkage maps of half-smooth tongue sole were developed with 1007 microsatellite markers, two SCAR markers and an F1 family containing 94. The female map was composed of 828 markers in 21 linkage groups, covering a total of 1447.3 cM, with an average interval 1.83 cM between markers. The male map consisted of 794 markers in 21 linkage groups, spanning 1497.5 cM, with an average interval of 1.96 cM. The female and male maps had 812 and 785 unique positions, respectively. The genome length of half-smooth tongue sole was estimated to be 1527.7 cM for the females and 1582.1 cM for the males. Based on estimations of the map lengths, the female and male maps covered 94.74 and 94.65% of the genome, respectively. The consensus map was composed of 1007 microsatellite markers and two SCAR markers in 21 linkage groups, covering a total of 1624 cM with an average interval of 1.67 cM. Furthermore, 159 sex-linked SSR markers were identified. Five sex-linked microsatellite markers were confirmed in their association with sex in a large number of individuals selected from different families. These sex-linked markers were mapped on the female map LG1f with zero recombination. Two QTLs that were identified for body weight, designated as We-1 and We-2, accounted for 26.39% and 10.60% of the phenotypic variation. Two QTLs for body width, designated Wi-1 and Wi-2, were mapped in LG4f and accounted for 14.33% and 12.83% of the phenotypic variation, respectively. Seven sex-related loci were mapped in LG1f, LG14f and LG1m by CIM, accounting for 12.5–25.2% of the trait variation. The results should prove to be very useful for improving growth traits using molecular MAS.

Exploitation of a turbot (Scophthalmus maximus L.) immune-related expressed sequence tag (EST) database for microsatellite screening and validation

In this study, we identified and characterized 160 microsatellite loci from an expressed sequence tag (EST) database generated from immune-related organs of turbot (Scophthalmus maximus). A final set of 83 new polymorphic microsatellites were validated after the analysis of 40 individuals of Atlantic origin including both wild and farmed individuals. The allele number and the expected heterozygosity ranged from 2 to 18 and from 0.021 to 0.951, respectively. Evidences of null alleles at moderate-high frequencies were detected at six loci using population data. None of the analysed loci showed deviations from Mendelian segregation after the analysis of five full-sib families including approximately 92 individuals/family. The markers are used to consolidate the turbot genetic map, and because they are mostly EST-derived, they will be very useful for comparative genomic studies within flatfishes and with model fish species. Using an in silico approach, we detected significant homologies of microsatellite sequences with the EST databases of the flatfish species with highest genomic resources (Senegalese sole, Atlantic halibut, bastard halibut) in 31% of these turbot markers. The conservation of these microsatellites within Pleuronectiformes will pave the way for anchoring genetic maps of different species and identifying genomic regions related to productive traits.

Development and characterization of 60 novel EST-SSR markers in half-smooth tongue sole Cynoglossus semilaevis

Sixty novel simple sequence repeat (SSR) markers were developed from expressed sequence tags (EST) of half-smooth tongue sole Cynoglossus semilaevis exploited in the laboratory. The number of alleles, observed and expected heterozygosity per locus ranged from two to 16, from 0·0833 to 1·0000 and from 0·0816 to 0·913, respectively. Of these SSRs, 20 had significant homology to known genes by BLASTx (basic local alignment search tool x) search. For cross-species amplification, there are 53 positive amplifications in Japanese flounder Paralichthys olivaceus with 12 polymorphic loci and 51 positive amplifications in Senegalese sole Solea senegalensis with 11 polymorphic loci. These new EST-SSR markers will be useful for genetic studies and genome mapping of C. semilaevis and its closely related fishes.

Genomic resources for flatfish research and their applications

Flatfishes are a group of teleosts of high commercial and environmental interest, whose biology is still poorly understood. The recent rapid development of different 'omic' technologies is, however, enhancing the knowledge of the complex genetic control underlying different physiological processes of flatfishes. This review describes the different functional genomic approaches and resources currently available for flatfish research and summarizes different areas where microarray-based gene expression analysis has been applied. The increase in genome sequencing data has also allowed the construction of genetic linkage maps in different flatfish species; these maps are invaluable for investigating genome organization and identifying genetic traits of commercial interest. Despite the significant progress in this field, the genomic resources currently available for flatfish are still scarce. Further intensive research should be carried out to develop larger genomic sequence databases, high-density microarrays and, more detailed, complete linkage maps, using second-generation sequencing platforms. These tools will be crucial for further expanding the knowledge of flatfish physiology, and it is predicted that they will have important implications for wild fish population management, improved fish welfare and increased productivity in aquaculture.

Construction of a genetic linkage map and mapping of a female-specific DNA marker in half-smooth tongue sole (Cynoglossus semilaevis)

The half-smooth tongue sole (Cynoglossus semilaevis, hereafter, "tongue sole") is a marine flatfish with great commercial importance for fisheries and aquaculture in China. It has also been a promising model for the study of sex determination mechanisms in fish. Here, we report the construction of a genetic linkage map for the tongue sole, based on 137 markers including 103 AFLP markers, 33 microsatellite markers, and one female-specific DNA marker. Twenty-six linkage groups (LGs) were found. The total map length was 934.6 cM (Kosambi), with an average spacing of 8.4 cM, covering 64.4% of the estimated genome size. Furthermore, a female-specific SCAR marker, CseF-382, was mapped on LG5. This study represents the first genetic linkage map in the tongue sole. This map has great potential in the identification of quantitative traits loci and sex-related genes and marker-assisted selection in the tongue sole. Meanwhile, the new set of polymorphic microsatellite markers developed in this study is not only useful for genetic mapping but also of critical importance for studies on genetic diversity and broodstock management in tongue sole.

Fosmid library construction and initial analysis of end sequences in female half-smooth tongue sole (Cynoglossus semilaevis)

Half-smooth tongue sole (Cynoglossus semilaevis: Pleuronectiformes) is a commercially important cultured marine flatfish in China and forms an important fishery resource, but the research of its genome is underdeveloped. In this study, we constructed a female C. semilaevis fosmid library and analyzed the fosmid end sequences to provide a preliminary assessment of the genome. The library consists of 49,920 clones with an average insert size of about 39 kb, amounting to 3.23 genome equivalents. Fosmid stability assays indicate that female C. semilaevis DNA was stable during propagation in the fosmid system. Library screening with eight microsatellite markers yielded between two and five positive clones, and none of those tested was absent from the library. End-sequencing of both 5' and 3' ends of 1,152 individual clones generated 2,247 sequences after trimming, with an average sequence length of 855 bp. BLASTN searches of the nr and EST databases of GenBank and BLASTX searches of the nr database resulted in 259 (11.53%) and 287 (12.77%) significant hits (E < e (-5)), respectively. Repetitive sequences analysis resulted in 5.23% of base pairs masked using both the Fugu and Danio databases, repetitive elements were composed of retroelements, DNA transposons, satellites, simple repeats, and low-complexity sequences. The fosmid library, in conjunction with the fosmid end sequences, will serve as a useful resource for large-scale genome sequencing, physical mapping, and positional cloning, and provide a better understanding of female C. semilaevis genome.