Sequencing and characterization of the transcriptome of half-smooth tongue sole (Cynoglossus semilaevis)
© Wang et al.; licensee BioMed Central Ltd. 2014
Received: 18 March 2013
Accepted: 4 June 2014
Published: 13 June 2014
Half-smooth tongue sole (Cynoglossus semilaevis) is a valuable fish for aquaculture in China. This fish exhibits sexual dimorphism, particularly different growth rates and body sizes between two genders. Thus, C. semilaevis is a good model that can be used to investigate mechanisms responsible for such dimorphism, this model can also be utilized to answer fundamental questions in evolution and applied fields of aquaculture. Hence, advances in second-generation sequencing technology, such as 454 pyrosequencing, could provide a robust tool to study the genome characteristics of non-model species.
In this study, C. semilaevis was subjected to de novo transcriptome sequencing and characterization. A total of 749,954 reads were generated using a single 454 sequencing run in a full PicoTiter plate. These reads were then assembled into 62,632 contigs with a 10-fold average sequencing coverage. A total of 26,589 sequences were successfully annotated based on sequence similarities; among these sequences, 3,451 transcripts exhibited gene ontology terms and 2,362 showed enzyme commissions associated with 186 pathways from Kyoto Encyclopedia of Gene and Genomes pathways. A search of repetitive elements was performed, and 1,898 transposable elements were identified. Approximately 7,800 simple-sequence repeats and 21,234 single-nucleotide polymorphisms were also detected.
Our data provided an integrated and comprehensive transcriptome resource for C. semilaevis. These data could be used for further research in population genetics, gene function, and tissue-specific gene expressions.
KeywordsCynoglossus semilaevis Half-smooth tongue sole Transcriptome SNP SSR
Half-smooth tongue sole (Cynoglossus semilaevis) is a large flatfish species naturally distributes in East Asia . This species has been considered as one of the most flavored species, which shows special advantages in a mild taste. Owing to limited wild resources, this flatfish has become an important farmed fish in China . Studies on this valuable aquatic fish have focused on the development of genetic markers, construction of genetic maps, and characterization of functional genes involved in growth, reproduction, stress, and immunity [3–6]. Recently, whole genome sequencing of C. semilaevis has been completed which made the study of species into a new phase .
Advances in high-throughput sequencing technologies have facilitated the studies on the genome and transcriptome of a non-model organism, such as C. semilaevis. Massively parallel 454 pyrosequencing, which shows its long-read characteristics (>400 bp), can be performed for de novo transcriptome analysis . Using 454 pyrosequencing, researchers sequenced and characterized the transcriptomes of many species [9–12].
The present study aimed to characterize the transcriptome of C. semilaevis. A multi-tissue and multi-individual library was constructed and sequenced using a 454 GS FLX titanium platform. Sequence assembly, gene annotation, transposable element (TE) analysis, and marker identification were performed.
Results and discussion
Sequence analysis and assembly
Summary of 454 transcriptome sequencing and assembly for C. semilaevis
Average length (bp)
Raw sequencing reads
The assembly of the trimmed and size-selected reads produced 62,632 contigs and 98,262 singletons. The length of the contigs ranged from 100 bp to 1,665 bp with an average length of 272 bp and an N50 of 303 bp (Figure 1B; Table 1). The average sequencing coverage, which was determined as the number of reads assembled into a particular contig, was 10.2. A positive relationship between the length of a specific contig and its coverage was observed in a randomly fragmented transcriptome (Figure 1C) . Using cd-hit , we performed cluster analysis with a standard such that the sequences with a similarity to this transcriptome of >95% were clustered into one class; the longest sequence of each class was used as a representative sequence. A total of 150,039 representative sequences, which may be unigenes, were obtained.
The complete genome sequence of C. semilaevis has been obtained, in which a 477 M genome with a scaffold N50 size of 867 kb was assembled . To assess our trancriptome assembly, we mapped the unigenes to C. semilaevis genome. Approximately 93.2% of unigenes exhibited significant hits on the genome. The mapping rates were 94.2% and 91.6% for isotigs and singletons, respectively. A high mapping rate indicated that our established assembly was assembly.
Compared with other fish transcriptomes using 454-pyrosequencing
Average length of raw reads
Numbers of raw reads
Total bases (M)
Average length of contigs
Number of contigs
Total bases (M)
1, 665, 609
1, 416, 404
Several complementary approaches were used to annotate the assembled sequences. The unigenes were initially compared with those in public protein databases by using BLASTX . With this procedure, the gene names of 26,589 (17.7%) sequences were successfully assigned (Additional file 1). A low annotation rate was mainly attributed to a short unigene because the significance of the BLAST comparison partially depends on the length of query sequence; thus, short reads obtained from sequencing rarely match known genes . In the present study, the annotation rate of short sequences (<300 bp, 14.5%) was lower than that of long sequences (>300 bp, 30.8%). Another reason accounted for low assignment percentage was the lack of information regarding flatfish. For example, only 349 sequences accounting for 1.31% of the total annotated sequences were annotated using the known information on flatfish. The three species with the most hits to Pleuronectiformes were Paralichthys olivaceus (129, 0.49%), C. semilaevis (97, 0.36%), and Solea senegalensis (57, 0.21%).
The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway approach used to perform high-order functional annotation was implemented using the web tool KASS . A total of 2,362 unigenes were mapped to 186 different pathways (Additional file 3). Among these pathways, the three highest maps were metabolic pathways (KO01100), biosynthesis of secondary metabolites (KO01110), and microbial metabolism in diverse environments (KO01120).
Functional genes involved in growth and immunity
As a valuable aquaculture fish, half-smooth tongue sole exhibits economical traits that are of great importance to aquaculturists. The sequences obtained by pyrosequencing and the annotation information provided significant data to determine economically important traits, including growth, reproduction, stress, and immunity.
The transcripts responsible for growth (GO: 0040007) were identified in our dataset. Seven sequences were identified as growth-related hormones. Among these seven sequences, three correspond to growth hormone-inducible transmembrane proteins, two correspond to growth hormone receptor 2, one corresponds to a predicted potassium channel subfamily K member, and the remaining sequence corresponds to a growth hormone-releasing hormone receptor.
MHC, a group of genes that encode for major histocompatibility antigens, functions in the immune response of vertebrates. Approximately 28 sequences have been annotated as MHC genes . This finding is consistent with that of the variability of MHC genes.
Cytochrome P450 (CYP) is a part of a large family of heme enzymes that catalyze diverse chemical reactions, including epoxidation, hydroxylation, and heteroatom oxidation . Steroidogenic enzyme is a member of the family P450 aromatase (P450arom); this enzyme is responsible for the conversion of androgen to estrogen and suppression of P450 arom gene expression; as a result, phenotypic sex reversal in fish occurs [23, 24]. The gender ratio of the cultivated half-smooth tongue sole was not 1:1 because the number of males was higher than that of the females. We identified 49 sequences annotated as P450; GO terms and KEGG pathways of these sequences were also detected.
Transposable elements identification
Previous studies reported that TEs are involved in the formation of the sex chromosomes of fish [28–30]. C. semilaevis contains a pair of well-differentiated sex chromosomes, and a large W chromosome is easily recognizable [31, 32]. However, the mechanism by which the W chromosome is formed in C. semilaevis remains unclear. In another study (unpublished), C. semilaevis contains a large number of TEs accumulated in the heterochromatic regions of the W chromosome. Our data sets are important to understand this phenomenon.
SSR and SNP discovery
A total of 7,869 SSRs located in 6,997 sequences were identified in our 454 pyrosequencing data set. The most common repeat motifs were dinucleotides, which accounted for 64.3% of the total SSRs; other repeat motifs were trinucleotides (31.1%), tetranucleotides (3.5%), pentanucleotides (0.7%), and hexanucleotides (0.4%). On the basis of the distribution of SSR motifs, we found that AC motifs were the most common (20.2%); CAG and AAAC repeat motifs were the most abundant motifs in tri- and tetra-nucleotides, respectively.
Molecular markers should be developed to facilitate marker-assisted selection to optimize commercial species and improve the available genetic resources of species. Transcriptome sequencing is a rapid and effective method to identify SSR and SNP markers. These markers identified by transcriptome sequencing are possibly associated with specific traits. In our trancriptome sequencing, 7,869 SSRs were identified, but the effectiveness should be validated. Previous studies may provide references for validation. For instance, Gao et al. identified 4,952 putative SSRs in the transcriptome of blunt snout bream (Megalobrama amblycephala) by 454 pyrosequencing, and 116 of 160 (72.5%) SSRs have been validated by PCR. In addition, 71 (44.4%) of these SSRs were polymorphic across a panel of 40 individuals. The C. semilaevis transcriptomes in this study and our study showed that dinucleotide repeats accounted for the highest proportions (64.3% and 60.5%, respectively). However, the proportion of dinucleotide repeats is higher (5,291/6,501, 81.4%) in Clupea harengus transcriptome ; this result indicated that dinucleotide repeats are the most common SSR type.
PCR and sequencing errors resulted in false positives to detect SNP. We implemented strict criteria to reduce the effects of sequencing errors. These criteria were listed as follows: (1) Q value of the bases was restricted at <23; (2) minor allele frequency was >15%; and (3) the minimum number of minor allele reads was set at 2. Hence, a total of 21,234 putative SNPs and 13,370 putative single-nucleotide indels were identified. These putative SNPs included 14,333 transitions and 6,901 transversions. The overall frequency of all of the SNP types and indels in the transcriptome was 1 per 491 bp.
The relationship of number putative SNPs and indels and number of minor allele reads
We performed the de novo transcriptome sequencing of half-smooth tongue sole (C. semilaevis) by using a 454 FLX titanium platform. Our results revealed a large number of candidate genes potentially involved in growth, reproduction, and stress/immunity response. Putative TEs were detected and analyzed. Moreover, numerous SNPs and SSRs were identified and prepared to perform marker development. Our data set provided a useful resource for future genetic or genomic studies on this species.
All of the experimental animal programs involved in this study were approved by the Ocean University of China’s Animal Care and Use Committee and performed in accordance with the experimental basic principles.
Biological materials and RNA extraction
Six adult C. semilaevis (three males and three females) individuals were obtained from the Yellow Sea Aquatic Product Co., Ltd. and temporarily housed in the laboratory without feeding for 2 d. These individuals were euthanized using 300 mg/L tricaine methanesulfonate (MS222). Afterward, different tissues, including brain, gill, heart, liver, kidney, spleen, intestines (the end part), muscle, testis, and ovary, were removed, immediately frozen in liquid nitrogen, and stored at -80°C. Total RNA was extracted from these tissues by using Trizol reagent (Invitrogen, USA). The quantity and quality of the total RNA were determined by spectrophotometry (Bio-Spec-mini, Shimadzu, Japan) and gel electrophoresis. Equal quantities of high-quality RNA of each tissue from different individual fish were pooled for cDNA synthesis.
cDNA library construction and pyrosequencing
cDNA library was constructed following the protocol described by Meyer et al. . In brief, the RNAs from the same tissue of the six C. semilaevis individuals were pooled before the first-strand cDNA was synthesized using a SMART PCR cDNA synthesis kit (Clontech, USA). The second-strand cDNA was amplified using an Advantage® 2 PCR kit (Clontech, USA). The PCR products were purified using a QIA quick® PCR purification kit (QIAGEN, Germany).
cDNA samples were normalized using a TRIMMER cDNA normalization kit (Evrogen, Russia) to balance the number of transcripts with high and low expressions. After normalization was performed, cDNA samples were sheared by sonication using an ultrasonic crusher (JY-92IIDN, Ningbo Xinzhi Biotech, China) to produce fragments with a length ranging from 300 bp to 800 bp, which is the appropriate fragment size range in 454 pyrosequencing.
Oligonucleotide adaptors contain a barcode sequence to discriminate samples from different tissues. Before sequencing was performed, all of the libraries were combined into a single pool and 5 μg of the mixed cDNA sample was used in a full PicoTiter plate of 454 GS FLX titanium (Roche, Switzerland) according to the standard manufacturer’s instructions (Shanghai Oebiotech Co., Ltd., China).
Sequence data analysis and assembly
The obtained raw reads were initially pre-processed by removing the adaptors and the primers using SeqClean (latest x86_64, Dana-Farber Cancer Institute) and Newbler (version 2.5.3, Roche). Low-quality reads were removed using Lucy (version 1.20p, -m50 –e 0.03 0.03 –w 30 0.03 10 0.1 –b 4 0.03) . High-quality reads with a length of >50 bp were maintained. The trimmed and size-selected reads were then assembled using Newbler with the following parameters: “use duplicate reads”; “extend low-depth overlaps”; “minimum read length = 45 bp”; “read limited to one contig”; and “single ACE file”. After the assembly was constructed, contigs and singletons were subjected to cluster analysis using cd-hit (version 4.0) with the following parameters: -r 1 and –c 0.95. The longest sequence of each class was used as a representative sequence. The alignment between unigenes and C. semilaevis genome was performed using BLASTN with e-value < 1e–10.
The unigenes were compared with the NCBI non-redundant protein database by using BLASTX. The alignment was processed using the following parameters: (1) e-value < 1e–5; (2) identity > 50%; and (3) number of matched amino acids >50. If the length of a unigene was <150 bp, the number of matched amino acids should be >70% of its own unigene. Gene function was assigned to each unigene based on the most common BLAST hit.
GO annotation was performed using Blast2GO, a software package that retrieves GO terms; this software package can be used to determine and compare gene functions [36–38]. KEGG pathways were assigned to the unigenes by using the online KEGG Automatic Annotation Server (KAAS; http://www.genome.jp/kegg/kaas/).
Transposable elements identification
Putative TEs were identified based on a homology search. Our data sets were compared with RepBase 17.09 using tBLASTx with a threshold of 1e–5 . The outputs were manually inspected, and significant matches to Simple Repeat, Pseudogene, and Integrated Virus were excluded.
SSR and SNP discovery
SSRs were identified using the MISA program. The minimum number of repetitions of dinucleotide was set at six. By comparison, trinucleotide, tetranucleotide, pentanucleotide, and hexanucleotide were set at five. The maximum difference between two SSRs was 100 bp.
Putative SNPs and single-nucleotide indels were detected using a ssahaSNP program . A putative SNP and indel site should satisfy the following conditions: (1) Q value of bases > 22; (2) minor allele frequency > 15%; and (3) minor allele was supported by at least two reads.
We thank the Yellow Sea Aquatic Product Co., Ltd. for providing the half-smooth tongue sole materials. We also thank Mr. Guoliang Wang for TEs analysis, and Mr. Linhe Huo for artwork preparation. Financial support for this work was provided by the National High-tech Research and Development Program (2012AA10A402 and 2012AA10A401).
Availability of supporting data
The raw reads of C. semilaevis transcriptome sequencing can be available at NCBI Short Read Archive with the accession number SRP020479 (http://www.ncbi.nlm.nih.gov/Traces/sra/sra.cgi?study=SRP020479).
- Li SZ, Wang HM: Chinese Fauna Osteichthyes Pleuronectiformes. 1995, Beijing: Academic PressGoogle Scholar
- Ma A, Liu X, Xu Y, Liang Y, Zhuang Z: Feeding rhythm and growth of the tongue sole, Cynoglossus semilaevis Günther, during its early life stages. Aquaculture Res. 2006, 37: 586-593. 10.1111/j.1365-2109.2006.01466.x.View ArticleGoogle Scholar
- Wang Y, Liu X, Ma L, Yu Y, Yu H, Mohammed S, Chu G, Mu L, Zhang Q: Identification and characterization of a hepcidin from half-smooth tongue sole Cynoglossus semilaevis. Fish Shellfish Immunol. 2012, 33: 213-219. 10.1016/j.fsi.2012.04.011.PubMedView ArticleGoogle Scholar
- Wang X, Zhang Q, Sun X, Chen Y, Zhai T, Zhuang W, Qi J, Wang Z: Fosmid library construction and initial analysis of end sequences in female half-smooth tongue sole (Cynoglossus semilaevis). Mar Biotechnol (NY). 2009, 11: 236-242. 10.1007/s10126-008-9137-2.View ArticleGoogle Scholar
- Song W, Pang R, Niu Y, Gao F, Zhao Y, Zhang J, Sun J, Shao C, Liao X, Wang L, Tian Y, Chen S: Construction of High-Density Genetic Linkage Maps and Mapping of Growth-Related Quantitative Trail Loci in the Japanese Flounder (Paralichthys olivaceus). PLoS One. 2012, 7 (11): e50404-10.1371/journal.pone.0050404.PubMed CentralPubMedView ArticleGoogle Scholar
- Yu Y, Zhong Q, Li C, Jiang L, Sun Y, Wang X, Zhang Q: Molecular cloning and characterization of interleukin-1beta in half-smooth tongue sole Cynoglossus semilaevis. Vet Immunol Immunopathol. 2012, 146: 270-276. 10.1016/j.vetimm.2012.02.011.PubMedView ArticleGoogle Scholar
- Chen S, Zhang G, Shao C, Huang Q, Liu G, Zhang P, Song W, An N, Chalopin D, Volff J-N: Whole-genome sequence of a flatfish provides insights into ZW sex chromosome evolution and adaptation to a benthic lifestyle. Nat Genet. 2014, 46 (3): 253-262. 10.1038/ng.2890.PubMedView ArticleGoogle Scholar
- Rothberg JM, Leamon JH: The development and impact of 454 sequencing. Nat Biotechnol. 2008, 26: 1117-1124. 10.1038/nbt1485.PubMedView ArticleGoogle Scholar
- Meyer E, Aglyamova GV, Wang S, Buchanan-Carter J, Abrego D, Colbourne JK, Willis BL, Matz MV: Sequencing and de novo analysis of a coral larval transcriptome using 454 GSFlx. BMC Genomics. 2009, 10: 219-236. 10.1186/1471-2164-10-219.PubMed CentralPubMedView ArticleGoogle Scholar
- Clark MS, Thorne MAS, Toullec J-Y, Meng Y, Guan LL, Peck LS, Moore S: Antarctic Krill 454 Pyrosequencing Reveals Chaperone and Stress Transcriptome. PLoS One. 2011, 6 (1): e15919-10.1371/journal.pone.0015919.PubMed CentralPubMedView ArticleGoogle Scholar
- Coppe A, Pujolar JM, Maes GE, Larsen PF, Hansen MM, Bernatchez L, Zane L, Bortoluzzi S: Sequencing, de novo annotation and analysis of the first Anguilla anguilla transcriptome: EeelBase opens new perspectives for the study of the critically endangered European eel. BMC Genomics. 2010, 11: 635-647. 10.1186/1471-2164-11-635.PubMed CentralPubMedView ArticleGoogle Scholar
- Hou R, Bao Z, Wang S, Su H, Li Y, Du H, Hu J, Wang S, Hu X: Transcriptome sequencing and de novo analysis for Yesso scallop (Patinopecten yessoensis) using 454 GS FLX. PLoS One. 2011, 6 (6): e21560-10.1371/journal.pone.0021560.PubMed CentralPubMedView ArticleGoogle Scholar
- Li W, Godzik A: Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006, 22: 1658-1659. 10.1093/bioinformatics/btl158.PubMedView ArticleGoogle Scholar
- Fraser BA, Weadick CJ, Janowitz I, Rodd FH, Hughes KA: Sequencing and characterization of the guppy (Poecilia reticulata) transcriptome. BMC Genomics. 2011, 12: 202-215. 10.1186/1471-2164-12-202.PubMed CentralPubMedView ArticleGoogle Scholar
- Salem M, Rexroad CE, Wang J, Thorgaard GH, Yao J: Characterization of the rainbow trout transcriptome using Sanger and 454-pyrosequencing approaches. BMC Genomics. 2010, 11: 564-573. 10.1186/1471-2164-11-564.PubMed CentralPubMedView ArticleGoogle Scholar
- Vera JC, Wheat CW, Fescemyer HW, Frilander MJ, Crawford DL, Hanski I, Marden JH: Rapid transcriptome characterization for a nonmodel organism using 454 pyrosequencing. Mol Ecol. 2008, 17: 1636-1647. 10.1111/j.1365-294X.2008.03666.x.PubMedView ArticleGoogle Scholar
- Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, Madden TL: BLAST+: architecture and applications. BMC Bioinformatics. 2009, 10: 421-429. 10.1186/1471-2105-10-421.PubMed CentralPubMedView ArticleGoogle Scholar
- Novaes E, Drost DR, Farmerie WG, Pappas GJ, Grattapaglia D, Sederoff RR, Kirst M: High-throughput gene and SNP discovery in Eucalyptus grandis, an uncharacterized genome. BMC Genomics. 2008, 9: 312-325. 10.1186/1471-2164-9-312.PubMed CentralPubMedView ArticleGoogle Scholar
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. Nat Genet. 2000, 25: 25-29. 10.1038/75556.PubMed CentralPubMedView ArticleGoogle Scholar
- Moriya Y, Itoh M, Okuda S, Yoshizawa AC, Kanehisa M: KAAS: an automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 2007, 35: W182-W185. 10.1093/nar/gkm321.PubMed CentralPubMedView ArticleGoogle Scholar
- Li C, Zhang Q, Yu Y, Li S, Zhong Q, Sun Y, Wang Z, Qi J, Zhai J, Wang X: Sequence polymorphism of two major histocompatibility (MH) class II B genes and their association with Vibrio anguillarum infection in half-smooth tongue sole (Cynoglossus semilaevis). Chin J Oceanol Limnol. 2011, 29: 1275-1286. 10.1007/s00343-011-0278-1.View ArticleGoogle Scholar
- Bistolas N, Wollenberger U, Jung C, Scheller FW: Cytochrome P450 biosensors–a review. Biosens Bioelectron. 2005, 20: 2408-2423. 10.1016/j.bios.2004.11.023.PubMedView ArticleGoogle Scholar
- Kitano T, Takamune K, Nagahama Y, Abe SI: Aromatase inhibitor and 17alpha-methyltestosterone cause sex-reversal from genetical females to phenotypic males and suppression of P450 aromatase gene expression in Japanese flounder (Paralichthys olivaceus). Mol Reprod Dev. 2000, 56: 1-5. 10.1002/(SICI)1098-2795(200005)56:1<1::AID-MRD1>3.0.CO;2-3.PubMedView ArticleGoogle Scholar
- Kitano T, Takamune K, Kobayashi T, Nagahama Y, Abe SI: Suppression of P450 aromatase gene expression in sex-reversed males produced by rearing genetically female larvae at a high water temperature during a period of sex differentiation in the Japanese flounder (Paralichthys olivaceus). J Mol Endocrinol. 1999, 23: 167-176. 10.1677/jme.0.0230167.PubMedView ArticleGoogle Scholar
- Finnegan DJ: Eukaryotic transposable elements and genome evolution. Trends Genet. 1989, 5: 103-107.PubMedView ArticleGoogle Scholar
- Wicker T, Sabot F, Hua-Van A, Bennetzen JL, Capy P, Chalhoub B, Flavell A, Leroy P, Morgante M, Panaud O, Paux E, SanMiguel P, Schulman AH: A unified classification system for eukaryotic transposable elements. Nat Rev Genet. 2007, 8: 973-982. 10.1038/nrg2165.PubMedView ArticleGoogle Scholar
- Kapitonov VV, Jurka J: A universal classification of eukaryotic transposable elements implemented in Repbase. Nat Rev Genet. 2008, 9: 411-412. 10.1038/nrg2165-c1.PubMedView ArticleGoogle Scholar
- de Bello CM, Sánchez A, Marchal JA, Kosyakova N, Liehr T, Trifonov V, Bertollo LAC: Whole chromosome painting reveals independent origin of sex chromosomes in closely related forms of a fish species. Genetica. 2011, 139: 1065-1072. 10.1007/s10709-011-9610-0.View ArticleGoogle Scholar
- Cioffi M, Molina W, Moreira-Filho O, Bertollo L: Chromosomal distribution of repetitive DNA sequences highlights the independent differentiation of multiple sex chromosomes in two closely related fish species. Cytogenet Genome Res. 2011, 134: 295-302. 10.1159/000329481.PubMedView ArticleGoogle Scholar
- Cabral-de-Mello D, Valente G, Nakajima R, Martins C: Genomic organization and comparative chromosome mapping of the U1 snRNA gene in cichlid fish, with an emphasis in Oreochromis niloticus. Chromosome Res. 2012, 20: 279-292. 10.1007/s10577-011-9271-y.PubMedView ArticleGoogle Scholar
- Zhou L, Yang A, Liu X, Du W, Zhuang Z: The karyotype of the tonguefish Cynoglossus semilaevis. J Fisheries Chin. 2005, 3: 417-419.Google Scholar
- Wang X, Zhang Q, Ren J, Jiang Z, Wang C, Zhuang W, Zhai T: The preparation of sex-[hromosome-specific painting probes and construction of sex chromosome DNA library in half-smooth tongue sole (Cynoglossus semilaevis). Aquaculture. 2009, 297: 78-84. 10.1016/j.aquaculture.2009.09.020.View ArticleGoogle Scholar
- Helyar SJ, Limborg MT, Bekkevold D, Babbucci M, Van Houdt J, Maes GE, Bargelloni L, Nielsen RO, Taylor MI, Ogden R: SNP discovery using next generation transcriptomic sequencing in Atlantic Herring (Clupea harengus). PLoS One. 2012, 7 (8): e42089-10.1371/journal.pone.0042089.PubMed CentralPubMedView ArticleGoogle Scholar
- Gao Z, Luo W, Liu H, Zeng C, Liu X, Yi S, Wang W: Transcriptome analysis and SSR/SNP markers information of the blunt snout bream (Megalobrama amblycephala). PLoS One. 2012, 7 (8): e42637-10.1371/journal.pone.0042637.PubMed CentralPubMedView ArticleGoogle Scholar
- Li S, Chou HH: LUCY2: an interactive DNA sequence quality trimming and vector removal tool. Bioinformatics. 2004, 20: 2865-2866. 10.1093/bioinformatics/bth302.PubMedView ArticleGoogle Scholar
- Conesa A, Götz S: Blast2GO: a comprehensive suite for functional analysis in plant genomics. Int J Plant Genomics. 2008, doi:10.1155/2008/619832Google Scholar
- Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21: 3674-3676. 10.1093/bioinformatics/bti610.PubMedView ArticleGoogle Scholar
- Götz S, García-Gómez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, Robles M, Talón M, Dopazo J, Conesa A: High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008, 36: 3420-3435. 10.1093/nar/gkn176.PubMed CentralPubMedView ArticleGoogle Scholar
- Jurka J, Kapitonov V, Pavlicek A, Klonowski P, Kohany O, Walichiewicz J: Repbase Update, a database of eukaryotic repetitive elements. Cytogenet Genome Res. 2005, 110: 462-467. 10.1159/000084979.PubMedView ArticleGoogle Scholar
- Ning Z, Caccamo M, Mullikin JC: ssahaSNP-a polymorphism detection tool on a whole genome scale. Computational Systems Bioinformatics Conference, 2005 Workshops and Poster Abstracts IEEE: 2005. 2005, Stanford: IEEE, 251-252.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited.