Figures
Abstract
Background
Body color and coloration patterns are important phenotypic traits to maintain survival and reproduction activities. The Oujiang color varieties of common carp (Cyprinus carpio var. color), with a narrow distribution in Zhejiang Province of China and a history of aquaculture for over 1,200 years, consistently exhibit a variety of body color patterns. The molecular mechanism underlying diverse color patterns in these variants is unknown. To the practical end, it is essential to develop molecular markers that can distinguish different phenotypes and assist selective breeding.
Methodology/Principal Findings
In this exploratory study, we conducted Roche 454 transcriptome sequencing of two pooled skin tissue samples of Oujiang common carp, which correspond to distinct color patterns, red with big black spots (RB) and whole white (WW), and a total of 737,525 sequence reads were generated. The reads obtained in this study were co-assembled jointly with common carp Roche 454 sequencing reads downloaded from NCBI SRA database, resulting in 43,923 isotigs and 546,676 singletons. Over 31 thousand (31,445; 71.6%) isotigs were found with significant BLAST matches (E<1e-10) to the nr protein database, which corresponds to 12,597 annotated zebrafish genes. A total of 70,947 isotigs and singletons (transcripts) were annotated with Gene Ontology, and 60,221 transcripts were found with corresponding EC numbers. Out of 145 zebrafish pigmentation genes, orthologs for 117 were recovered in Oujiang color carp transcriptome, including 18 found only among singletons. Our transcriptome analysis revealed over 52,902 SNPs in Oujiang common carp, and identified 63 SNP markers that are putatively unique either for RB or WW.
Conclusions
The transcriptome of Oujiang color varieties of common carp obtained through this study, along with the pigmentation genes recovered and the color pattern-specific molecular markers developed, will facilitate future research on the molecular mechanism of color patterns and promote aquaculture of Oujiang color varieties of common carp through molecular marker assisted-selective breeding.
Citation: Wang C, Wachholtz M, Wang J, Liao X, Lu G (2014) Analysis of the Skin Transcriptome in Two Oujiang Color Varieties of Common Carp. PLoS ONE 9(3): e90074. https://doi.org/10.1371/journal.pone.0090074
Editor: László Orbán, Temasek Life Sciences Laboratory, Singapore
Received: June 6, 2013; Accepted: January 29, 2014; Published: March 6, 2014
Copyright: © 2014 Wang et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: Funding was provided by grants from the Special Fund for Agro-Scientific Research in the Public Interest (No. 200903045), National Natural Science Foundation of China (No. 30972250), and by the Shanghai University Knowledge Service Platform (ZF1206). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Body color and coloration patterns are important phenotypic traits associated with the survival and reproduction activities in many organisms, including their camouflage, mimicry, social communication and selective mating [1], [2], [3]. Coloration is determined mainly by pigments synthesized by chromatophores, i.e., pigment cells. So far, six types of pigment cells, including melanocytes (black, dark brown), xanthophores (yellow), erythrophores (red/orange), iridophores (reflecting), leucophores (white), and cyanophores (blue), have been reported in vertebrates [4]. Mammals and birds were found to have only melanocytes in their skin, whereas reptiles also possess xanthophores, erythrophores and iridophores. Interestingly, only teleost fishes were found to have all six types of pigment cells [5], [6], [7].
The molecular mechanism of melanin biosynthesis has been xextensively studied due to its biomedical significance; the melanogenesis pathways have been found conserved in vertebrates [8], [9], [10]. Melanin synthesis takes place within the melanosomes of melanocytes and involves the tyrosine metabolism pathway. The tyrosinases (TYR), tyrosinase-related protein 1 (TYRP1) and dopachrome tautomerase (DCT) are important enzymes in melanin synthesis [8], [11]. The most important regulator of melanogenesis is the MC1 receptor (MC1R) with its melanocortic peptides [10]. A single mutation of this gene was found to be associated with the beach mouse color pattern [12]. MC1R activates the cyclic AMP (cAMP) response-element binding protein (CREB), an important cellular transcription factor. Increased expression of microphthalmia-associated transcription factor (MITF) and its activation by phosphorylation (P) stimulate the transcription of TYR, TYRP1 and DCT [7], [8], [13], [14]. The biosynthesis of other types of pigments may involve different pathways and regulatory networks, but are in general less studied [6].
The common carp (Cyprinus carpio L.), a cyprinid belonging to Cyprinidae of ray-finned fishes (teleosts), is considered the most widely distributed and important aquaculture species in the world [15]. During the long history of domestication, numerous strains or variants of common carp have been established through artificial selection from the wild ancestral populations, of which many have diverse colors and coloration patterns [16], [17]. The molecular mechanisms underlying these diverse color patterns in common carp remain less well-understood.
Oujiang color varieties of common carp (Cyprinus carpio var. color), a unique variant of common carp, has been cultured in paddies and backyard ponds for about 1,200 years in the Oujiang river basin of Zhejiang province in China [18]. It has been found that several body colors and coloration patterns consistently exist, including whole red (WR), red with scattered black spots (RB), whole white (WW), and white with scattered black spots (WB), which provide an excellent model to explore the molecular mechanisms of pigment formation and development [16].
In a mating experiment involving whole red (WR) and whole white (WW) Oujiang color common carp varieties, we observed a ratio of three red to one white, providing evidence of Mendelian inheritance in the red and white body colors [19]. Although the black color showed dominant inheritance over red and white, it does not seem to follow the Mendelian law, suggesting that the regulation of this pattern is more complicated than the others (Wang et al., unpublished). In a study of red and black color inheritance with Koi, an ornamental variety developed from the common carp, three phenotypes were found associated with the red color and two phenotypes with the black one [20].
In our previous study, significant genetic differences were revealed among the different varieties of Oujiang color common carp through the use of microsatellite markers [21]. In addition, our field growth experiment demonstrated that RB grew faster than other body color types of Oujiang color common carp, while WW grew at the slowest rate [22]. Developing more genetic markers will not only facilitate the studies of population genetics of Oujiang color varieties of common carp, but it will also improve its aquaculture production through the identification of coloration patterns at early developmental stages.
In recent years, the transcriptome analysis via next-generation sequencing technology has brought new insight into the knowledge of whole transcriptomes in many organisms [23], [24], [25], [26], [27], [28]. This powerful approach allows the detection of significantly expressed genes and important pathways involved in different developmental stages and/or physiological processes [25], [27], [29], [30], [31], [32]. In addition, the large amount of transcriptomic sequences provides a rich resource for the discovery of different types of molecular markers [25], [26], [33], [34], [35], [36], [37].
Here, we conducted transcriptome analysis of Oujiang color common carp skin and obtained the first transcriptomic repertoire of this unique common carp strain. We also investigated the pigmentation genes expressed in the skin of Oujiang color varieties of common carp and found a set of candidate SNPs markers that can be used for selective molecular breeding in the future.
Materials and Methods
Ethics statement
This study was approved by the Institutional Animal Care and Use Committee (IACUS) of Shanghai Ocean University, Shanghai, China. The collection of fish samples was permitted by the Provincial Farm of Oujiang Color Common Carp in Zhejiang Province, China, and all sampling procedures were complied with the guidelines of IACUS on the care and use of animals for scientific purposes.
Sample collection
Oujiang color common carp samples (tetraploid) were collected from six siblings, including three whole white (WW) and three red with big black spots (RB), from the Provincial Farm of Oujiang Color Common Carp in Zhejiang Province, China. Fresh skin tissues (∼3 g) were collected and preserved immediately in liquid nitrogen for RNA extraction.
RNA extraction and mRNA purification
Each frozen sample was ground in mortars with liquid nitrogen, and total RNA was isolated using TRIzol reagent (Invitrogen) following the standard protocol. The concentration of total RNA dissolved in 200 µL RNase-free water was determined by NanoDrop (Thermo Scientific, USA), whereas the RNA integrity value (RIN) was estimated using RNA 6000 Pico LabChip of Agilent 2100 Bioanalyzer (Agilent, USA). Total RNA of each pooled sample of three individuals from the same color pattern was incubated with 10 U DNase I (Ambion) at 37°C for 1 hr, and then nuclease-free water was added to bring the sample volume to 250 µL. Messenger RNA was further purified with MicroPoly(A) Purist Kit (Ambion) as per the recommended protocol by the manufacturer. The mRNA was dissolved in 100 µL THE RNA Storage Solution, and the final concentration was determined using NanoDrop.
cDNA synthesis
Double-stranded cDNA was synthesized from mRNA according to Ng's full-length cDNA synthesis protocol with minor modifications [38]. A GsuI-oligo dT primer was used to prime first-strand cDNA synthesis from 10 µg of mRNA, using 1000 units of Superscript II reverse transcriptase (Invitrogen). After incubation at 42°C for 1 hr, the 5′-CAP structure of mRNA was oxidized by NaIO4 (Sigma) and ligated to biotin hydrazide, which was used to select complete mRNA/cDNA by binding Dynal M280 beads (Invitrogen). After the second strand cDNA synthesis, the polyA and 5′ adaptor were removed through GsuI digestion.
cDNA sequencing
cDNA size fractionation was performed using a cDNA size fractionation column (Agencourt). Each cDNA fraction larger than 800 bp was sonicated to the range of 300–800 bp, and then pooled together with the other cDNA samples within the same range of lengths. The prepared cDNA of each pooled sample of three individuals from the same color pattern was transformed into single-stranded template DNA (sstDNA) libraries by using the GS DNA Library Preparation Kit (Roche Applied Science). The sstDNA libraries were clonally amplified in a bead-immobilized form by using the GS emPCR Kit (Roche Applied Science) and sequenced on the Roche 454 Genome Sequencer FLX instrument.
Sequence assembly and annotation
Sequence assembly, downstream analysis and annotation were conducted following the pipeline illustrated in Figure 1. The raw sequencing reads generated in this study have been deposited to the NCBI Short Read Archive (SRA) database (SRA 291286). The reads from Oujiang color common carp (RB and WW) along with the gynogenic common carp (CC) originating from Roche 454 sequencing reads downloaded from NCBI (SRA050545) were assembled using Newbler 2.6 [39].
The pictures of RB, WW and CC are shown at the top panel. CC raw sequencing reads were downloaded from NCBI SRA database (SRA 050545).
The assembled transcripts (isotigs and singletons) were BLASTed against NCBI nr protein database as well as zebrafish protein sequences, downloaded from ENSEMBL database, with a cutoff E-value smaller than 1e-10 [40]. Gene Ontology analysis was performed with the BLAST2GO program [41]. The BLAST program was also used to identify fish pigmentation genes [4], [6] in Oujiang color common carp transcriptome.
Mapping and SNP Identification
All reads were mapped to the de novo common carp transcriptome assembled in this study to identify candidate SNP markers using Newbler GS mapper [39]. To find a consensus set of candidate SNPs markers, we also used BWA 0.6.1 [42] for mapping and FaSD [43] for SNPs detection. SNPs detected by the above two approaches and with at least three reads mapped were considered as putative markers specific to RB or WW.
Results and Discussion
Transcriptome of common carp and Oujiang color varieties
The statistics of common carp sequencing and assembly are summarized in Table 1. Altogether, 737,525 raw reads were generated via Roche 454 sequencing technology with an average length of 252 bp and maximal length of 680 bp. The reads generated in this study and the over two million Roche 454 sequencing reads from common carp available at NCBI SRA database (SRA050545) were assembled into 43,923 isotigs, with 546,676 reads unassembled (singletons). A total of 344,319 (63.0%) singletons were from the two color varieties of common carp and 75,260 singletons were found to have significant Blast hits. The N50 value, the average length and the largest length of isotigs were 1,237, 1,000 and 11,663 bp, respectively. The BLAST search resulted in 31,445 isotigs (71.6%) with significant matches to the NCBI nr database and 30,304 isotigs (68.9%) to ENSEMBL zebrafish proteins (E-value<1e-10). In addition, 12,597 zebrafish genes (about half of the zebrafish transcriptome) were found with significant BLAST matches to the assembled common carp transcriptome, indicating that approximately half of the common carp transcriptome can be annotated. For singletons, 128,507 out of 546,676 (23.5%) were found to have significant BLAST hits.
Further analysis of Gene Ontology (GO) showed that 25,564 isotigs and 45,383 singletons were successfully mapped to Gene Ontology categories - molecular functions, biological processes, and cellular components (Figure 2). For molecular functions, the majority of the transcripts (including both isotigs and singletons) were associated with binding, protein binding, and catalytic activity. For cellular component, the final products of these transcripts (proteins) were located mainly within cell, which is probably not surprising because the skin tissue was sampled for RNA extraction. For biological processes, most of these transcripts were related to cellular process and metabolic process, which was in agreement with most of other studies [24], [44]. A total of 134 KEGG pathways were found in the common carp transcriptome with the top 20 pathways shown in Table 2. The pathways with most genes represented were purine metabolism, glycolysis and gluconeogenesis, and oxidative phosphorylation. Tyrosine metabolism, a pigmentation related pathway, was found with 37 isotigs and 43 singletons that correspond to 23 proteins.
The transcriptome was assembled from Oujiang color varieties of common carp (RB - red with big black spots and WW – whole white) and gynogenic common carp (CC).
In our study, approximately 30% of the isotigs and 80% of the singletons were found without significant BLAST matches to known proteins in public databases. In fact, most transcriptome studies showed a portion of isotig and singleton sequences that cannot be annotated [24], [25], [26], [34], [35], [45]. One possible reason for this is that these non-model organisms possess numerous putative novel genes or transcripts, of which the sequences are not available in public databases. Further analysis and annotation of these putative novel transcripts are necessary, but appear to be challenging. With the availability of transcriptomic sequences in many organisms, the comparative analysis of closely related species may assist in the identification of novel orthologous genes. Besides the widely used BLAST algorithm, the profile- or model-based approaches such as PRINTS [46], PROSITE [47], and Pfam [48] can be used to predict the functions of isotig/singleton sequences.
Pigmentation genes in Oujiang color varieties of common carp
After a BLAST search with 145 zebrafish-derived pigmentation genes against the common carp transcriptome, 593 isotigs were found with hits to 99 pigmentation genes and 501 singletons with matches to 95 pigmentation genes (E-value<1e-10). Of the 501 singletons, 291 were from WW and 210 from RB. A total of 77 pigmentation genes were commonly found in both isotigs and singletons, with 22 occurring solely in isotigs and 18 found only in singletons. When comparing the mapped reads of these pigmentation genes in RB and WW, 187 isotigs were found with at least a two-fold change, including 106 higher in WW and 81 in RB. In addition, 134 RB isotigs were not found in WW whereas 100 WW isotigs were not found in RB (Table S1).
Isolation of SNP markers specific to the RB or WW color varieties
Reads from Oujiang color varieties of common carp and from NCBI SRA database were mapped to the assembled de novo common carp transcriptome. Over 11,000 isotigs were mapped with reads from RB, WW, and CC (Figure 3). There were about 1,403 and 635 isotigs uniquely mapped with RB and WW reads, respectively, suggesting that more specific genes were involved in the formation of RB phenotype compared with WW. Many more unique isotigs were mapped with reads only from CC, which may be due to the much larger number of original reads from CC (Figures 1 and 3).
RB - red with big black spots, WW - whole white.
A total of 52,902 single nucleotide polymorphism (SNP) markers were found in 9,379 isotigs, including 3,721 single nucleotide insertions or deletions (indels). The frequency distribution of different substitutions or indels is shown in Figure 4. The highest SNP pattern was G/A, accounting for 16% of the total. Compared to other substitution types, the substitution between C and G was found to be the least frequent.
* indicates single nucleotide polymorphisms other than single nucleotide substitution, including deletion.
To discover reliable markers that can be used to distinguish RB and WW, we set up the following three criteria: 1) detected with two algorithms; 2) homozygous in either RB or WW; and 3) more than two mapped reads. When these criteria were applied, 63 putative SNPs markers were found, including 21 SNPs with RB showing the reference nucleotide, while WW showing the variant one and 42 SNPs with the opposite (Table 3). In addition, the analysis of pigmentation-related transcripts revealed only 10 SNPs that may be associated with different coloration phenotypes. One of our future research directions is to validate the predicted markers and evaluate their accuracy in identifying specific Oujiang color varieties of common carp with distinct color patterns [36], [37], [49].
Conclusions
The de novo skin transcriptome sequencing of Oujiang color common carp varieties were conducted, resulting in 737,525 reads assembled to 43,923 isotigs. Over 10,000 genes were annotated using BLAST, Gene Ontology, and KEGG pathways. The orthologs of approximately 80% (117 out of 145) zebrafish pigmentation genes were recovered. In addition, 52,902 SNPs were detected, with 63 candidate markers predicted. Future research directions include expanding this study to analyze more samples of other color patterns and validating the results on potentially differentially expressed genes between the RB and WW phenotypes. We are also planning to validate color pattern-specific molecular markers, which will presumably lead to a better understanding of the molecular mechanisms underlying diverse color patterns and promote aquaculture of Oujiang color varieties of common carp through the application of molecular markers in selective breeding.
Supporting Information
Table S1.
Pigmentation genes found in Oujiang color varieties of common carp (RB - red with big black spots and WW - whole white).
https://doi.org/10.1371/journal.pone.0090074.s001
(XLSX)
Table S2.
The complete list of KEGG pathways identified in the common carp skin transcriptome involving two Oujiang color common carp variants.
https://doi.org/10.1371/journal.pone.0090074.s002
(DOCX)
Table S3.
The complete list of putative SNP markers found in Oujiang color varieties of common carp (RB - red with big black spots, WW - whole white).
https://doi.org/10.1371/journal.pone.0090074.s003
(DOCX)
Acknowledgments
We acknowledge Mr. Songping Xiang and Jian Wang for the assistance in fish culture and sample collection. In addition, we want to thank Mr. Xinxin Yang with DNA content measurement, library construction and data analysis. Finally, we are grateful to Mary Christman and Julian Egger for proofreading the manuscript.
Author Contributions
Conceived and designed the experiments: CW. Performed the experiments: CW. Analyzed the data: MW JW XL. Contributed reagents/materials/analysis tools: CW GL. Wrote the paper: CW GL JW.
References
- 1. Protas ME, Patel NH (2008) Evolution of coloration patterns. Annu Rev Cell Dev Biol 24: 425–446.
- 2. Rodgers GM, Kelley JL, Morrell LJ (2010) Colour change and assortment in the western rainbowfish. Animal Behaviour 79: 1025–1030.
- 3. Steiner CC, Weber JN, Hoekstra HE (2007) Adaptive variation in beach mice produced by two interacting pigmentation genes. PLoS biology 5: 1–e219.
- 4. Braasch I, Schartl M, Volff JN (2007) Evolution of pigment synthesis pathways by gene and genome duplication in fish. BMC Evol Biol 7: 74.
- 5. Bagnara JT, Fernandez PJ, Fujii R (2007) On the blue coloration of vertebrates. Pigment Cell Res 20: 14–26.
- 6. Braasch I, Brunet F, Volff JN, Schartl M (2009) Pigmentation pathway evolution after whole-genome duplication in fish. Genome Biol Evol 1: 479–493.
- 7. Mellgren EM, Johnson SL (2002) The evolution of morphological complexity in zebrafish stripes. Trends Genet 18: 128–134.
- 8. Li S, Wang C, Yu W, Zhao S, Gong Y (2012) Identification of genes related to white and black plumage formation by RNA-Seq from white and black feather bulbs in ducks. PLoS One 7: e36592.
- 9. Logan DW, Burn SF, Jackson IJ (2006) Regulation of pigmentation in zebrafish melanophores. Pigment Cell Res 19: 206–213.
- 10. Richardson J, Lundegaard PR, Reynolds NL, Dorin JR, Porteous DJ, et al. (2008) mc1r Pathway regulation of zebrafish melanosome dispersion. Zebrafish 5: 289–295.
- 11. Braasch I, Liedtke D, Volff JN, Schartl M (2009) Pigmentary function and evolution of tyrp1 gene duplicates in fish. Pigment Cell Melanoma Res 22: 839–850.
- 12. Hoekstra HE, Hirschmann RJ, Bundey RA, Insel PA, Crossland JP (2006) A single amino acid mutation contributes to adaptive beach mouse color pattern. Science 313: 101–104.
- 13. Sugie A, Terai Y, Ota R, Okada N (2004) The evolution of genes for pigmentation in African cichlid fishes. Gene 343: 337–346.
- 14. Altschmied J, Delfgaauw J, Wilde B, Duschl J, Bouneau L, et al. (2002) Subfunctionalization of duplicate mitf genes associated with differential degeneration of alternative exons in fish. Genetics 161: 259–267.
- 15. Balon EK (1995) Origin and domestication of the wild carp,Cyprinus carpio: from Roman gourmets to the swimming flowers. Aquaculture 129: 3–48.
- 16. Wang C, Li S, Nagy ZT, Lehoczky I, Huang L, et al. (2010) Molecular genetic structure and relationship of Chinese and Hungarian common carp (Cyprinus carpio L.) strains based on mitochondrial sequence. Aquaculture Research 41: 1339–1347.
- 17. David L, Rothbard S, Rubinstein I, Katzman H, Hulata G, et al. (2004) Aspects of red and black color inheritance in the Japanese ornamental (Koi) carp (Cyprinus carpio L.). Aquaculture 233: 129–147.
- 18. Li S (2003) Aquaculture research and its relation to development in China. Agricultural development and the opportunities for aquatic resources research in China 65: 17.
- 19. Wang CH, Xiang SP, Lv YP, Wang J, Pan ZY, et al. (2008) Preliminary study on inheritance of red and white color phenotypes in Oujiang color carp, C. carpio var. color. Journal of Shanghai Fisheries University 4: 004.
- 20. David L, Blum S, Feldman MW, Lavi U, Hillel J (2003) Recent Duplication of the Common Carp (Cyprinus carpio L.) Genome as Revealed by Analyses of Microsatellite Loci. Molecular Biology and Evolution 20: 1425–1434.
- 21. Ma Y, Wang C, Wang J, Yang X, Bi X, et al. (2011) Genetic differentiation of wild and hatchery Oujiang color common carp: potential application to the identification of escapees. Fisheries Science 77: 591–597.
- 22. Zhu L, Ma Y, Xiang S, Bi X, Wang J, et al. (2013) Observation and analysis of growth dynamics in Oujiang color common carp with different pigmentation types. Journal of Shanghai Fisheries University In Press.
- 23. Wang Z, Gerstein M, Snyder M (2009) RNA-Seq: a revolutionary tool for transcriptomics. Nature Reviews Genetics 10: 57–63.
- 24. Ji P, Liu G, Xu J, Wang X, Li J, et al. (2012) Characterization of common carp transcriptome: sequencing, de novo assembly, annotation and comparative genomics. PLoS One 7: e35152.
- 25. Liao X, Cheng L, Xu P, Lu G, Wachholtz M, et al. (2013) Transcriptome Analysis of Crucian Carp (Carassius auratus), an Important Aquaculture and Hypoxia-Tolerant Species. PLoS One 8: e62308.
- 26. Salem M, Rexroad CE 3rd, Wang J, Thorgaard GH, Yao J (2010) Characterization of the rainbow trout transcriptome using Sanger and 454-pyrosequencing approaches. BMC Genomics 11: 564.
- 27. Samanta MP, Tongprasit W, Istrail S, Cameron RA, Tu Q, et al. (2006) The transcriptome of the sea urchin embryo. Science 314: 960–962.
- 28. Shen GM, Dou W, Niu JZ, Jiang HB, Yang WJ, et al. (2011) Transcriptome analysis of the oriental fruit fly (Bactrocera dorsalis). PLoS One 6: e29127.
- 29. Huan P, Wang H, Liu B (2012) Transcriptomic analysis of the clam Meretrix meretrix on different larval stages. Marine Biotechnology 14: 69–78.
- 30. Mu Y, Ding F, Cui P, Ao J, Hu S, et al. (2010) Transcriptome and expression profiling analysis revealed changes of multiple signaling pathways involved in immunity in the large yellow croaker during Aeromonas hydrophila infection. BMC Genomics 11: 506.
- 31. Visser WE, Swagemakers SM, Őzgűr Z, Schot R, Verheijen FW, et al. (2010) Transcriptional profiling of fibroblasts from patients with mutations in MCT8 and comparative analysis with the human brain transcriptome. Human molecular genetics 19: 4189–4200.
- 32. Li C, Zhang Y, Wang R, Lu J, Nandi S, et al. (2012) RNA-seq analysis of mucosal immune responses reveals signatures of intestinal barrier disruption and pathogen entry following Edwardsiella ictaluri infection in channel catfish, Ictalurus punctatus. Fish & Shellfish Immunology 32: 816–827.
- 33. Chouvarine P, Cooksey AM, McCarthy FM, Ray DA, Baldwin BS, et al. (2012) Transcriptome-based differentiation of closely-related Miscanthus lines. PLoS One 7: e29850.
- 34. Ma K, Qiu G, Feng J, Li J (2012) Transcriptome analysis of the oriental river prawn, Macrobrachium nipponense using 454 pyrosequencing for discovery of genes and markers. PLoS One 7: e39727.
- 35. Pereiro P, Balseiro P, Romero A, Dios S, Forn-Cuni G, et al. (2012) High-throughput sequence analysis of turbot (Scophthalmus maximus) transcriptome using 454-pyrosequencing for the discovery of antiviral immune genes. PLoS One 7: e35369.
- 36. Renaut S, Nolte AW, Bernatchez L (2010) Mining transcriptome sequences towards identifying adaptive single nucleotide polymorphisms in lake whitefish species pairs (Coregonus spp. Salmonidae). Mol Ecol 19 Suppl 1: 115–131.
- 37. Ji P, Zhang Y, Li C, Zhao Z, Wang J, et al. (2012) High throughput mining and characterization of microsatellites from common carp genome. Int J Mol Sci 13: 9798–9807.
- 38. Ng P, Wei CL, Sung WK, Chiu KP, Lipovich L, et al. (2005) Gene identification signature (GIS) analysis for transcriptome characterization and genome annotation. Nat Methods 2: 105–111.
- 39. Margulies M, Egholm M, Altman WE, Attiya S, Bader JS, et al. (2005) Genome sequencing in microfabricated high-density picolitre reactors. Nature 437: 376–380.
- 40. Flicek P, Amode MR, Barrell D, Beal K, Brent S, et al. (2012) Ensembl 2012. Nucleic acids research 40: D84–D90.
- 41. Conesa A, Gotz S, Garcia-Gomez JM, Terol J, Talon M, et al. (2005) Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21: 3674–3676.
- 42. Li H, Durbin R (2010) Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics 26: 589–595.
- 43. Xu F, Wang W, Wang P, Jun Li M, Chung Sham P, et al. (2012) A fast and accurate SNP detection algorithm for next-generation sequencing data. Nat Commun 3: 1258.
- 44. Bai X, Mamidala P, Rajarapu SP, Jones SC, Mittapalli O (2011) Transcriptomics of the bed bug (Cimex lectularius). PLoS One 6: e16336.
- 45. Calduch-Giner JA, Bermejo-Nogales A, Benedito-Palos L, Estensoro I, Ballester-Lozano G, et al. (2013) Deep sequencing for de novo construction of a marine fish (Sparus aurata) transcriptome database with a large coverage of protein-coding transcripts. BMC Genomics 14: 178.
- 46. Attwood TK, Blythe MJ, Flower DR, Gaulton A, Mabey J, et al. (2002) PRINTS and PRINTS-S shed light on protein ancestry. Nucleic acids research 30: 239–241.
- 47. Sigrist CJ, Cerutti L, de Castro E, Langendijk-Genevaux PS, Bulliard V, et al. (2010) PROSITE, a protein domain database for functional characterization and annotation. Nucleic acids research 38: D161–D166.
- 48. Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, et al. (2012) The Pfam protein families database. Nucleic acids research 40: D290–D301.
- 49. Salem M, Vallejo RL, Leeds TD, Palti Y, Liu S, et al. (2012) RNA-Seq identifies SNP markers for growth traits in rainbow trout. PLoS One 7: e36264.