Transcriptome characterisation of Pinus tabuliformis and evolution of genes in the Pinusphylogeny
© Niu et al.; licensee BioMed Central Ltd. 2013
Received: 8 November 2012
Accepted: 15 April 2013
Published: 18 April 2013
The Chinese pine (Pinus tabuliformis) is an indigenous conifer species in northern China but is relatively underdeveloped as a genomic resource; thus, limiting gene discovery and breeding. Large-scale transcriptome data were obtained using a next-generation sequencing platform to compensate for the lack of P. tabuliformis genomic information.
The increasing amount of transcriptome data on Pinus provides an excellent resource for multi-gene phylogenetic analysis and studies on how conserved genes and functions are maintained in the face of species divergence. The first P. tabuliformis transcriptome from a normalised cDNA library of multiple tissues and individuals was sequenced in a full 454 GS-FLX run, producing 911,302 sequencing reads. The high quality overlapping expressed sequence tags (ESTs) were assembled into 46,584 putative transcripts, and more than 700 SSRs and 92,000 SNPs/InDels were characterised. Comparative analysis of the transcriptome of six conifer species yielded 191 orthologues, from which we inferred a phylogenetic tree, evolutionary patterns and calculated rates of gene diversion. We also identified 938 fast evolving sequences that may be useful for identifying genes that perhaps evolved in response to positive selection and might be responsible for speciation in the Pinus lineage.
A large collection of high-quality ESTs was obtained, de novo assembled and characterised, which represents a dramatic expansion of the current transcript catalogues of P. tabuliformis and which will gradually be applied in breeding programs of P. tabuliformis. Furthermore, these data will facilitate future studies of the comparative genomics of P. tabuliformis and other related species.
KeywordsPinus tabuliformis Carr 454 pyrosequencing SNPs SSRs Pinus phylogeny Comparative transcriptomics
Conifers are widely distributed globally as the largest and most diverse group of gymnosperms  that evolved independently from angiosperms >300 million years ago . Modern conifers are divided into eight families including 68 genera and 630 species, which form an integral part of the economy in many parts of the world . Chinese pine (Pinus tabuliformis Carr.) is a widespread indigenous conifer species and an economically and ecologically important hard pine in northern China [4, 5]. Because of its irreplaceable economic development and environmental protection status, a genetic improvement program for P. tabuliformis was initiated in the 1970s, and considerable progress has been made in many basic physiological aspects . The study of natural genetic variation in P. tabuliformis has traditionally been investigated using a common garden approach, whereas the pace of development of genomic resources has been slow, as only 288 P. tabuliformis entries are included in the NCBI database. Information regarding the genetic control of many important traits and fine-scale genetic variations is extremely limited, and more is needed given the renewed emphasis to accelerate the pace of P. tabuliformis breeding and shorten the breeding cycle.
Despite the economic and ecological importance of the genus Pinus, the progress of entire genome sequencing and associated marker development has been limited [6, 7]. Huge genomes with highly heterozygous and large amounts of repetitive DNA elements are the major obstacles towards sequencing the genomes of all Pinus spp. [8, 9]. The genome sizes of conifers are larger than those of most other plant species. The genome in all extant members of the genus Pinus is 18,000–40,000 Mbp . In contrast, several representative genera of angiosperm trees have genome sizes of 540–2,000 Mb . Therefore, researchers have focused on the transcribed part of the genome using dedicated technologies [6, 7]. Transcriptome analysis and construction of large-scale expressed sequence tag (EST) collections in pines are a promising means of providing genomic resources [2, 9, 11], as this technique produces expressed sequence portions of chromosomes at a fraction of the cost of sequencing the complete genome . It also facilitates the analysis of the transcribed part of the genome, which is not easy to predict from the entire genome . Next-generation sequencing is a viable and favourable alternative to Sanger sequencing and provides researchers with a relatively rapid and affordable option for developing genomic resources in non-model organisms [14–16]. The Roche 454 massively parallel pyrosequencing platform, GS FLX Titanium, can generate one million reads with an average read length of 400 bases at 99.5% accuracy per run [17, 18].
In addition to the discovery of new genes and investigations of gene expression, thousands of simple sequence repeats (SSRs), single nucleotide polymorphisms (SNPs) and insertions and deletions (Indels) have been detected in transcriptome data [6, 19]. It is possible to use these genome-wide and abundant markers to develop very dense genetic maps that can be applied to conduct marker-assisted selection breeding programs .
Moreover, the increasing availability of transcriptome data represents an excellent resource for comparative genomic analysis. Although there has been much work on the chloroplast DNA sequences (cpDNA) and mitochondria DNA sequences (mtDNA), based on phylogenetic analysis of Pinus[21–23], less emphasis has been placed on multi-gene phylogenetic analysis and on determination of how conserved genes and functions are maintained despite species divergence.
In the current study, we used the Roche 454 GS-FLX Titanium pyrosequencing platform to obtain a comprehensive transcriptome of P. tabuliformis from normalised cDNA libraries of adult trees (xylem, phloem, vascular cambium, needles, cones and strobili). As a result, thousands of molecular markers were characterised. Evolutionary studies based on these data and other shared transcriptome data of five pine species and one spruce species were conducted. These data provide compelling new insights into the transcriptome of P. tabuliformis and evolution of genes in the Pinus phylogeny.
Transcriptome sequencing and de novoassembly
Sequencing, assembly and data analysis
Raw results (after trimming)
Total number of reads
Total number of isotigs
Total read length (bp)
295 125 234
Total isotigs length (bp)
21 076 176
Minimum read length (bp)
Isotig N50 (bp)
Median read length (bp)
Maximum isotig length (bp)
Maximum read length (bp)
Mean read length (bp)
Number of singletons
GC content (%)
Total number of unigenes
The unigene coverage distribution revealed that most unigenes had a read-depth coverage <20-fold (Figure 1c, d). The steep decline in read-depth coverage suggests that cDNA normalisation was effective, which is typical for a normalised library . Isotig lengths were related to the number of sequences assembled into each isotig. The average unigene length exhibited a gradual increase with increasing read depth (Figure 1c, d).
Functional annotation of the transcriptome
Identification of SSRs, SNPs and Indels
Orthologue identification and functional characterisation between six conifer species
Large-scale transcriptome characterisations have been carried out for Pinus taeda, Pinus contorta, Pinus sylvestris and Pinus pinaster. The shared transcriptomes of Pinus in the PlantGDB and NCBI databases are valuable sources of information for multi-gene comparative and phylogenetic analyses .
Phylogenetic and speciation analysis
Evolutionary pattern of Pinusspp. genes
Transcriptome sequencing in Pinus spp.
Platform or approach
Reads or ESTs
GS XL R70
needles and conelets
Sanger and GS-FLX
3 968 794
Despite the fact that a large number of pine ESTs have been obtained from cDNA libraries based on traditional sequencing technology, the methods used were inefficient. Four P. taeda cDNA libraries were sequenced and yielded a total of 142,533 ESTs (Table 2); however, only one normalised cDNA library yielded 822,891 ESTs in P. tabuliformis (Table 1). Although traditional sequencing yields longer EST sequences, it has little advantage compared to new assembly technology based on the large-scale ESTs. Additionally, most previous studies used a cDNA library of one tissue (Table 2), whereas we used a normalised cDNA library comprising multiple tissues and individuals. These large-scale ESTs will provide more comprehensive pine transcriptome information and facilitate the assembly of Pinus spp. ESTs in the future.
Next-generation sequencing technology yields a large number of sequences at considerably lower costs compared to traditional sequencing methods, and, therefore, provides a valuable starting point to expedite analysis of less-studied species [18, 36, 37]. Normalised cDNA libraries were used to sample large numbers of transcripts to maximise sequence diversity. Next-generation sequencing of normalised libraries is more efficient than that of non-normalised libraries, particularly for rare transcripts . The capacity to deliver large numbers of gene-based markers from transcriptome sequencing projects is a major advantage of next-generation sequencing technology [18, 20, 36]. Because of cost and throughput, conventional markers such as restriction fragment length polymorphism and random amplified polymorphic DNA are being replaced with SSRs and SNPs . The genome-wide and abundant EST-based SSRs and SNPs/Indels markers obtained by next-generation sequencing represent an effective approach to marker discovery in many plant species, as these markers facilitate generation of dense genetic maps and have the advantage of higher cross-species transferability [6, 39–41]. However, relevant studies in Pinus spp. are limited. In this study, 724 distinct EST-SSR loci and more than 92,000 SNPs/InDels were identified. It is possible to use these markers in a broad range of applications, including genetic mapping, genotype identification, marker-assisted selection breeding, and molecular tagging of genes. Among the EST-derived SSRs, tri-nucleotide repeat units were predominant. Considering the importance of maintaining reading frames to generating a polypeptide within a partially or fully active, it is no surprise that this observation is common for tri-nucleotide expansions (or their multiples) within translated regions [6, 42, 43]. As usual, comparisons of P. tabuliformis transcriptome SNPs show an excess of transitional over transversional substitutions. Similarly, A and T were the most frequent insertion and deletion types of Indels. Part of this bias is due to the relatively high rate of mutation of methylated cytosines to thymines [44, 45].
Comparative phylogenetic analysis at the genome level dramatically improves the precision and sensitivity of evolutionary inference . However, comparative genomics in plants has been limited by the considerable phylogenetic distances between sequenced organisms . Transcriptome sequencing using massively parallel sequencing technologies provides an attractive approach to obtaining large-scale sequence data for non-model organisms necessary for comparative genomic analysis [24, 48]. Phylogenetic utility of transcriptome sequence data yields well-resolved and highly supported tree topologies for many groups of animals [49–51]; however, few such studies have been conducted with plant taxa . Phylogenetic analysis of the genus Pinus has been limited mostly to plastid genome (cpDNA and mtDNA) sequences [21–23]. The results of this study are consistent with previous data on plastid genome phylogeny [21, 22], but transcriptome analyses, producing more robust results, are presented for the first time. Given that this study was not limited to particular genes or motifs, the results presented here are more representative of Pinus evolution than previous studies.
Understanding the factors that affect the evolutionary patterns and rates of genes is central in many research fields . For the past 30 years, it was thought that the rate of gene evolution was determined by protein function . Studies on yeast and bacteria indicate that the expression level of a protein affects the evolution rate more than its functional category, at least in unicellular species [54, 55]. In this study, we have shown that sequence polymorphisms of the 191 putatively orthologous sets of ESTs of six Pinus species are widespread using GO terms. This suggests that selection of protein function does not contribute to the variation in the rates of gene evolution. However, most of the important factors are correlated with each other. More systematic analyses of genomic data are required to further demonstrate the effect of a range of factors on the evolutionary patterns and rates of genes.
This study is the first comprehensive sequencing effort and analysis of gene function in the transcriptome of P. tabuliformis and represents the most extensive expressed sequence resource available for P. tabuliformis to date. GO and KEGG analyses were carried out, and all unigenes were classified into functional categories so as to understand their functions and regulation pathways. An enormous number of SSR and SNP/Indel loci were detected. These data can be used to develop oligonucleotide microarrays or serve as a reference transcriptome for future RNA-seq experiments in large-scale gene expression assays. These data will accelerate our understanding of genetic variation in populations and the genetic control of important traits in P. tabuliformis. Additionally, the generation of such large-scale sequence data is a potentially invaluable scientific resource for mapping, marker-assisted breeding and conservation-genetic-oriented studies in P. tabuliformis and comparative evolutionary analysis of Pinus plants.
Sample collection, cDNA library creation and 454 sequencing
Samples used for sequencing
Samples in May
Number of individuals
Samples in July
Number of individuals
Cambium (stress side)
Cambium (stress side)
Cambium (tension side)
Cambium (tension side)
Needles (juvenile + mature)
Total RNA isolation from samples of all selected plant tissues, and cDNA library construction and normalisation were performed as described previously . The pooled library was sequenced in a full 454 plate run on the GS-FLX Titanium platform (Roche, Indianapolis, IN, USA).
Assembly and annotation
All generated ESTs were pre-screened to remove adaptor-ligated regions and contaminants by Seqclean and to trim low-quality regions by LUCY2 . Because no reference P. tabuliformis genome exists, cleaned and qualified reads were assembled de novo in Newbler 2.5.3, which performs best for restoring full-length transcripts [13, 58].
The assembled isotig and singleton sequences were combined and clustered with CD-HIT (version 4.0) [59, 60]. The sequences with similarity >95% were divided into one class, and the longest sequence of each class was treated as a unigene during later processing. Descriptive annotations and GO classifications were performed as described previously .
We simultaneously instituted a search for putative unigenes against the NCBI protein database using a BLASTx and annotated each sequence with GO terms using Blast2GO.
Identification of SSRs, SNPs and InDels
Assembled isotigs with coverage of at least four reads were screened for SSRs, SNPs and InDels using Misa and ssahaSNP software, respectively . Similar criteria for screening high-quality SNPs have been used in previous studies [20, 62]. Only perfect repeats of two to six nucleotide repeats were identified. The minimum repeat-unit size for di-nucleotides was set at six and at five for tri- to hexa-nucleotide repeats.
Identification of orthologues between six conifer species
The shared transcriptome data of five conifer species in the PlantGDB and NCBI databases were downloaded. The numbers of unigenes for each species were as follows: Picea glauca (48,619), P. contorta (13,570), P. pinaster (15,648), P. sylvestris (73,609) and P. taeda (77,540). Along with 46,584 unigenes of P. tabuliformis, clustering was carried out among the transcribed sequences using UCLUST software . Aligned sequences (at least 100 bp) showing 90% identity were defined as pairs of putative orthologues among six species. The best-hit sequence of each cluster was then used in subsequent analyses. Orthologues of P. taeda and P. tabuliformis were searched using the same approach. Sequences of P. tabuliformis were annotated with GO terms using Blast2GO.
Estimation at the level of synonymous substitution and non-synonymous substitution between orthologues
Because unigenes are derived from EST sequences, have no annotated open reading frames and may contain frame shift sequencing errors, each member of a pair of sequences was searched using BLASTX against all plant protein sequences available in GenBank. The approach used was as described previously . PAML software was used to estimate the non-synonymous substitutions per non-synonymous site (Ka) and the synonymous substitutions per synonymous site (Ks) .
Because the genus Pinus has a rich history of phylogenetic analysis and the relationships among the species in the genus are well understood [21–23], the precise topology is not critical for the purposes of this study. We chose to focus our analyses on the evolutionary pattern and rate of genes. The synonymous substitution and non-synonymous substitution between the orthologues of six conifer species were analysed as described previously. Phylograms were derived using pairwise substitution rates of orthologous transcripts as a distance metric with the neighbour-joining method . Picea glauca was used as an out-group to root trees.
The raw 454 EST data obtained in this study were deposited in the NCBI Sequence Read Archive (SRA) under the accession number SRA 056887.
This work was supported by grants from The Fundamental Research Funds for the Central Universities and the Special Fund for Scientific Forestry Research in the Public Interest (No.201104022).
- Ahuja MR, Neale DB: Evolution of genome size in conifers. Silvae Genet. 2005, 54 (3): 126-137.Google Scholar
- Fernandez-Pozo N, Canales J, Guerrero-Fernandez D, Villalobos DP, Diaz-Moreno SM, Bautista R, Flores-Monterroso A, Guevara MA, Perdiguero P, Collada C: EuroPineDB: a high-coverage web database for maritime pine transcriptome. BMC Genomics. 2011, 12: 366-10.1186/1471-2164-12-366.PubMed CentralView ArticlePubMedGoogle Scholar
- Cairney J, Zheng L, Cowels A, Hsiao J, Zismann V, Liu J, Ouyang S, Thibaud-Nissen F, Hamilton J, Childs K: Expressed sequence tags from loblolly pine embryos reveal similarities with angiosperm embryogenesis. Plant Mol Biol. 2006, 62 (4–5): 485-501.View ArticlePubMedGoogle Scholar
- Li W, Wang X, Li Y: Stability in and correlation between factors influencing genetic quality of seed lots in seed orchard of Pinus tabuliformis Carr. over a 12-year span. PLoS One. 2011, 6 (8): e23544-10.1371/journal.pone.0023544.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen K, Abbott RJ, Milne RI, Tian XM, Liu J: Phylogeography of Pinus tabulaeformis Carr. (Pinaceae), a dominant species of coniferous forest in northern China. Mol Ecol. 2008, 17 (19): 4276-4288. 10.1111/j.1365-294X.2008.03911.x.View ArticlePubMedGoogle Scholar
- Lesser MR, Parchman TL, Buerkle CA: Cross-species transferability of SSR loci developed from transciptome sequencing in lodgepole pine. Mol Ecol Resour. 2012, 12 (3): 448-455. 10.1111/j.1755-0998.2011.03102.x.View ArticlePubMedGoogle Scholar
- Ralph SG, Yueh H, Friedmann M, Aeschliman D, Zeznik JA, Nelson CC, Butterfield YS, Kirkpatrick R, Liu J, Jones SJ: Conifer defence against insects: microarray gene expression profiling of Sitka spruce (Picea sitchensis) induced by mechanical wounding or feeding by spruce budworms (Choristoneura occidentalis) or white pine weevils (Pissodes strobi) reveals large-scale changes of the host transcriptome. Plant Cell Environ. 2006, 29 (8): 1545-1570. 10.1111/j.1365-3040.2006.01532.x.View ArticlePubMedGoogle Scholar
- Valledor L, Jorrin JV, Rodriguez JL, Lenz C, Meijon M, Rodriguez R, Canal MJ: Combined proteomic and transcriptomic analysis identifies differentially expressed pathways associated to Pinus radiata needle maturation. J Proteome Res. 2010, 9 (8): 3954-3979. 10.1021/pr1001669.View ArticlePubMedGoogle Scholar
- Parchman TL, Geist KS, Grahnen JA, Benkman CW, Buerkle CA: Transcriptome sequencing in an ecologically important tree species: assembly, annotation, and marker discovery. BMC Genomics. 2010, 11: 180-10.1186/1471-2164-11-180.PubMed CentralView ArticlePubMedGoogle Scholar
- Morse AM, Peterson DG, Islam-Faridi MN, Smith KE, Magbanua Z, Garcia SA, Kubisiak TL, Amerson HV, Carlson JE, Nelson CD: Evolution of genome size and complexity in Pinus. PLoS One. 2009, 4 (2): e4332-10.1371/journal.pone.0004332.PubMed CentralView ArticlePubMedGoogle Scholar
- Neale DB: Genomics to tree breeding and forest health. Curr Opin Genet Dev. 2007, 17 (6): 539-544. 10.1016/j.gde.2007.10.002.View ArticlePubMedGoogle Scholar
- Emrich SJ, Barbazuk WB, Li L, Schnable PS: Gene discovery and annotation using LCM-454 transcriptome sequencing. Genome Res. 2007, 17 (1): 69-73.PubMed CentralView ArticlePubMedGoogle Scholar
- Mundry M, Bornberg-Bauer E, Sammeth M, Feulner PG: Evaluating characteristics of de novo assembly software on 454 transcriptome data: a simulation approach. PLoS One. 2012, 7 (2): e31410-10.1371/journal.pone.0031410.PubMed CentralView ArticlePubMedGoogle Scholar
- Hiremath PJ, Farmer A, Cannon SB, Woodward J, Kudapa H, Tuteja R, Kumar A, Bhanuprakash A, Mulaosmanovic B, Gujaria N: Large-scale transcriptome analysis in chickpea (Cicer arietinum L.), an orphan legume crop of the semi-arid tropics of Asia and Africa. Plant Biotechnol J. 2011, 9 (8): 922-931. 10.1111/j.1467-7652.2011.00625.x.PubMed CentralView ArticlePubMedGoogle 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-10.1186/1471-2164-12-202.PubMed CentralView ArticlePubMedGoogle Scholar
- Russell JR, Bayer M, Booth C, Cardle L, Hackett CA, Hedley PE, Jorgensen L, Morris JA, Brennan RM: Identification, utilisation and mapping of novel transcriptome-based markers from blackcurrant (Ribes nigrum). BMC Plant Biol. 2011, 11: 147-10.1186/1471-2229-11-147.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang Y, Zeng X, Iyer NJ, Bryant DW, Mockler TC, Mahalingam R: Exploring the switchgrass transcriptome using second-generation sequencing technology. PLoS One. 2012, 7 (3): e34225-10.1371/journal.pone.0034225.PubMed CentralView ArticlePubMedGoogle Scholar
- Kaur S, Pembleton LW, Cogan NO, Savin KW, Leonforte T, Paull J, Materne M, Forster JW: Transcriptome sequencing of field pea and faba bean for discovery and validation of SSR genetic markers. BMC Genomics. 2012, 13: 104-10.1186/1471-2164-13-104.PubMed CentralView ArticlePubMedGoogle Scholar
- Renaut S, Nolte AW, Bernatchez L: Mining transcriptome sequences towards identifying adaptive single nucleotide polymorphisms in Lake Whitefish species pairs (Coregonus spp. Salmonidae). Mol Ecol. 2010, 19 (Suppl 1): 115-131.View ArticlePubMedGoogle Scholar
- Barbazuk WB, Emrich SJ, Chen HD, Li L, Schnable PS: SNP discovery via 454 transcriptome sequencing. Plant J. 2007, 51 (5): 910-918. 10.1111/j.1365-313X.2007.03193.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Eckert AJ, Hall BD: Phylogeny, historical biogeography, and patterns of diversification for Pinus (Pinaceae): phylogenetic tests of fossil-based hypotheses. Mol Phylogenet Evol. 2006, 40 (1): 166-182. 10.1016/j.ympev.2006.03.009.View ArticlePubMedGoogle Scholar
- Gernandt DS, Lopez GG, Garcia SO, Liston A: Phylogeny and classification of Pinus. Taxon. 2005, 54 (1): 29-42. 10.2307/25065300.View ArticleGoogle Scholar
- Wang B, Mao JF, Gao J, Zhao W, Wang XR: Colonization of the Tibetan Plateau by the homoploid hybrid pine Pinus densata. Mol Ecol. 2011, 20 (18): 3796-3811. 10.1111/j.1365-294X.2011.05157.x.View ArticlePubMedGoogle Scholar
- Der JP, Barker MS, Wickett NJ, DePamphilis CW, Wolf PG: De novo characterization of the gametophyte transcriptome in bracken fern, Pteridium aquilinum. BMC Genomics. 2011, 12: 99-10.1186/1471-2164-12-99.PubMed CentralView ArticlePubMedGoogle Scholar
- DiGuistini S, Ralph SG, Lim YW, Holt R, Jones S, Bohlmann J, Breuil C: Generation and annotation of lodgepole pine and oleoresin-induced expressed sequences from the blue-stain fungus Ophiostoma clavigerum, a Mountain Pine Beetle-associated pathogen. FEMS Microbiol Lett. 2007, 267 (2): 151-158. 10.1111/j.1574-6968.2006.00565.x.View ArticlePubMedGoogle Scholar
- Heller G, Adomas A, Li G, Osborne J, van Zyl L, Sederoff R, Finlay RD, Stenlid J, Asiegbu FO: Transcriptional analysis of Pinus sylvestris roots challenged with the ectomycorrhizal fungus Laccaria bicolor. BMC Plant Biol. 2008, 8: 19-10.1186/1471-2229-8-19.PubMed CentralView ArticlePubMedGoogle Scholar
- Logacheva MD, Kasianov AS, Vinogradov DV, Samigullin TH, Gelfand MS, Makeev VJ, Penin AA: De novo sequencing and characterization of floral transcriptome in two species of buckwheat (Fagopyrum). BMC Genomics. 2011, 12: 30-10.1186/1471-2164-12-30.PubMed CentralView ArticlePubMedGoogle Scholar
- Buschiazzo E, Ritland C, Bohlmann J, Ritland K: Slow but not low: genomic comparisons reveals slower evolutionary rate and higher dN/dS in conifers compared to angiosperms. BMC Evol Biol. 2012, 12 (1): 8-10.1186/1471-2148-12-8.PubMed CentralView ArticlePubMedGoogle Scholar
- Li X, Wu HX, Dillon SK, Southerton SG: Generation and analysis of expressed sequence tags from six developing xylem libraries in Pinus radiata D. Don. BMC Genomics. 2009, 10: 41-10.1186/1471-2164-10-41.PubMed CentralView ArticlePubMedGoogle Scholar
- Li X, Wu HX, Southerton SG: Transcriptome profiling of Pinus radiata juvenile wood with contrasting stiffness identifies putative candidate genes involved in microfibril orientation and cell wall mechanics. BMC Genomics. 2011, 12: 480-10.1186/1471-2164-12-480.PubMed CentralView ArticlePubMedGoogle Scholar
- Li X, Wu HX, Southerton SG: Seasonal reorganization of the xylem transcriptome at different tree ages reveals novel insights into wood formation in Pinus radiata. New Phytol. 2010, 187 (3): 764-776. 10.1111/j.1469-8137.2010.03333.x.View ArticlePubMedGoogle Scholar
- Allona I, Quinn M, Shoop E, Swope K, St CS, Carlis J, Riedl J, Retzel E, Campbell MM, Sederoff R: Analysis of xylem formation in pine by cDNA sequencing. Proc Natl Acad Sci U S A. 1998, 95 (16): 9693-9698. 10.1073/pnas.95.16.9693.PubMed CentralView ArticlePubMedGoogle Scholar
- Kirst M, Johnson AF, Baucom C, Ulrich E, Hubbard K, Staggs R, Paule C, Retzel E, Whetten R, Sederoff R: Apparent homology of expressed genes from wood-forming tissues of loblolly pine (Pinus taeda L.) with Arabidopsis thaliana. Proc Natl Acad Sci U S A. 2003, 100 (12): 7383-7388. 10.1073/pnas.1132171100.PubMed CentralView ArticlePubMedGoogle Scholar
- Pavy N, Laroche J, Bousquet J, Mackay J: Large-scale statistical analysis of secondary xylem ESTs in pine. Plant Mol Biol. 2005, 57 (2): 203-224. 10.1007/s11103-004-6969-7.View ArticlePubMedGoogle Scholar
- Lorenz WW, Sun F, Liang C, Kolychev D, Wang H, Zhao X, Cordonnier-Pratt MM, Pratt LH, Dean JF: Water stress-responsive genes in loblolly pine (Pinus taeda) roots identified by analyses of expressed sequence tag libraries. Tree Physiol. 2006, 26 (1): 1-16. 10.1093/treephys/26.1.1.View ArticlePubMedGoogle Scholar
- Jhanwar S, Priya P, Garg R, Parida SK, Tyagi AK, Jain M: Transcriptome sequencing of wild chickpea as a rich resource for marker development. Plant Biotechnol J. 2012, 10 (6): 690-702. 10.1111/j.1467-7652.2012.00712.x.View ArticlePubMedGoogle Scholar
- Edwards CE, Parchman TL, Weekley CW: Assembly, gene annotation and marker development using 454 floral transcriptome sequences in Ziziphus celata (Rhamnaceae), a highly endangered, Florida endemic plant. DNA Res. 2012, 19 (1): 1-9. 10.1093/dnares/dsr037.PubMed CentralView ArticlePubMedGoogle Scholar
- Wall PK, Leebens-Mack J, Chanderbali AS, Barakat A, Wolcott E, Liang H, Landherr L, Tomsho LP, Hu Y, Carlson JE: Comparison of next generation sequencing technologies for transcriptome characterization. BMC Genomics. 2009, 10: 347-10.1186/1471-2164-10-347.PubMed CentralView ArticlePubMedGoogle Scholar
- Ellis JR, Burke JM: EST-SSRs as a resource for population genetic analyses. Heredity (Edinb). 2007, 99 (2): 125-132. 10.1038/sj.hdy.6801001.View ArticleGoogle Scholar
- Barbara T, Palma-Silva C, Paggi GM, Bered F, Fay MF, Lexer C: Cross-species transfer of nuclear microsatellite markers: potential and limitations. Mol Ecol. 2007, 16 (18): 3759-3767. 10.1111/j.1365-294X.2007.03439.x.View ArticlePubMedGoogle Scholar
- Bouck A, Vision T: The molecular ecologist's guide to expressed sequence tags. Mol Ecol. 2007, 16 (5): 907-924.View ArticlePubMedGoogle Scholar
- Rowland LJ, Alkharouf N, Darwish O, Ogden EL, Polashock JJ, Bassil NV, Main D: Generation and analysis of blueberry transcriptome sequences from leaves, developing fruit, and flower buds from cold acclimation through deacclimation. BMC Plant Biol. 2012, 12 (1): 46-10.1186/1471-2229-12-46.PubMed CentralView ArticlePubMedGoogle Scholar
- Blanca J, Canizares J, Roig C, Ziarsolo P, Nuez F, Pico B: Transcriptome characterization and high throughput SSRs and SNPs discovery in Cucurbita pepo (Cucurbitaceae). BMC Genomics. 2011, 12: 104-10.1186/1471-2164-12-104.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhao Z, Boerwinkle E: Neighboring-nucleotide effects on single nucleotide polymorphisms: a study of 2.6 million polymorphisms across the human genome. Genome Res. 2002, 12 (11): 1679-1686. 10.1101/gr.287302.PubMed CentralView ArticlePubMedGoogle Scholar
- Keller I, Bensasson D, Nichols RA: Transition-transversion bias is not universal: a counter example from grasshopper pseudogenes. PLoS Genet. 2007, 3 (2): e22-10.1371/journal.pgen.0030022.PubMed CentralView ArticlePubMedGoogle Scholar
- Clark AG, Eisen MB, Smith DR, Bergman CM, Oliver B, Markow TA, Kaufman TC, Kellis M, Gelbart W, Iyer VN: Evolution of genes and genomes on the Drosophila phylogeny. Nature. 2007, 450 (7167): 203-218. 10.1038/nature06341.View ArticlePubMedGoogle Scholar
- Kunstner A, Wolf JB, Backstrom N, Whitney O, Balakrishnan CN, Day L, Edwards SV, Janes DE, Schlinger BA, Wilson RK: Comparative genomics based on massive parallel transcriptome sequencing reveals patterns of substitution and selection across 10 bird species. Mol Ecol. 2010, 19 (Suppl 1): 266-276.PubMed CentralView ArticlePubMedGoogle Scholar
- Zakas C, Schult N, McHugh D, Jones KL, Wares JP: Transcriptome analysis and SNP development can resolve population differentiation of Streblospio benedicti, a developmentally dimorphic marine annelid. PLoS One. 2012, 7 (2): e31613-10.1371/journal.pone.0031613.PubMed CentralView ArticlePubMedGoogle Scholar
- Meusemann K, von Reumont BM, Simon S, Roeding F, Strauss S, Kuck P, Ebersberger I, Walzl M, Pass G, Breuers S: A phylogenomic approach to resolve the arthropod tree of life. Mol Biol Evol. 2010, 27 (11): 2451-2464. 10.1093/molbev/msq130.View ArticlePubMedGoogle Scholar
- Roeding F, Borner J, Kube M, Klages S, Reinhardt R, Burmester T: A 454 sequencing approach for large scale phylogenomic analysis of the common emperor scorpion (Pandinus imperator). Mol Phylogenet Evol. 2009, 53 (3): 826-834. 10.1016/j.ympev.2009.08.014.View ArticlePubMedGoogle Scholar
- Roeding F, Hagner-Holler S, Ruhberg H, Ebersberger I, von Haeseler A, Kube M, Reinhardt R, Burmester T: EST sequencing of Onychophora and phylogenomic analysis of Metazoa. Mol Phylogenet Evol. 2007, 45 (3): 942-951. 10.1016/j.ympev.2007.09.002.View ArticlePubMedGoogle Scholar
- Pal C, Papp B, Lercher MJ: An integrated view of protein evolution. Nat Rev Genet. 2006, 7 (5): 337-348. 10.1038/nrg1838.View ArticlePubMedGoogle Scholar
- McInerney JO: The causes of protein evolutionary rate variation. Trends Ecol Evol. 2006, 21 (5): 230-232. 10.1016/j.tree.2006.03.008.View ArticlePubMedGoogle Scholar
- Rocha EP: The quest for the universals of protein evolution. Trends Genet. 2006, 22 (8): 412-416. 10.1016/j.tig.2006.06.004.View ArticlePubMedGoogle Scholar
- Wall DP, Hirsh AE, Fraser HB, Kumm J, Giaever G, Eisen MB, Feldman MW: Functional genomic analysis of the rates of protein evolution. Proc Natl Acad Sci U S A. 2005, 102 (15): 5483-5488. 10.1073/pnas.0501761102.PubMed CentralView ArticlePubMedGoogle Scholar
- Garzon-Martinez GA, Zhu I, Landsman D, Barrero LS, Marino-Ramirez L: The Physalis peruviana leaf transcriptome: assembly, annotation and gene model prediction. BMC Genomics. 2012, 13 (1): 151-10.1186/1471-2164-13-151.PubMed CentralView ArticlePubMedGoogle Scholar
- Li S, Chou HH: LUCY2: an interactive DNA sequence quality trimming and vector removal tool. Bioinformatics. 2004, 20 (16): 2865-2866. 10.1093/bioinformatics/bth302.View ArticlePubMedGoogle Scholar
- Kumar S, Blaxter ML: Comparing de novo assemblers for 454 transcriptome data. BMC Genomics. 2010, 11: 571-10.1186/1471-2164-11-571.PubMed CentralView ArticlePubMedGoogle Scholar
- Huang Y, Niu B, Gao Y, Fu L, Li W: CD-HIT Suite: a web server for clustering and comparing biological sequences. Bioinformatics. 2010, 26 (5): 680-682. 10.1093/bioinformatics/btq003.PubMed CentralView ArticlePubMedGoogle Scholar
- Li W, Godzik A: Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006, 22 (13): 1658-1659. 10.1093/bioinformatics/btl158.View ArticlePubMedGoogle Scholar
- van Oeveren J, Janssen A: Mining SNPs from DNA sequence data; computational approaches to SNP discovery and analysis. Methods Mol Biol. 2009, 578: 73-91. 10.1007/978-1-60327-411-1_4.View ArticlePubMedGoogle Scholar
- Milano I, Babbucci M, Panitz F, Ogden R, Nielsen RO, Taylor MI, Helyar SJ, Carvalho GR, Espineira M, Atanassova M: Novel tools for conservation genomics: comparing two high-throughput approaches for SNP discovery in the transcriptome of the European hake. PLoS One. 2011, 6 (11): e28008-10.1371/journal.pone.0028008.PubMed CentralView ArticlePubMedGoogle Scholar
- Edgar RC: Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010, 26 (19): 2460-2461. 10.1093/bioinformatics/btq461.View ArticlePubMedGoogle Scholar
- Blanc G, Wolfe KH: Widespread paleopolyploidy in model plant species inferred from age distributions of duplicate genes. Plant Cell. 2004, 16 (7): 1667-1678. 10.1105/tpc.021345.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang Z: PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24 (8): 1586-1591. 10.1093/molbev/msm088.View ArticlePubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24 (8): 1596-1599. 10.1093/molbev/msm092.View ArticlePubMedGoogle Scholar
- Wan LC, Zhang H, Lu S, Zhang L, Qiu Z, Zhao Y, Zeng QY, Lin J: Transcriptome-wide identification and characterization of miRNAs from Pinus densata. BMC Genomics. 2012, 13: 132-10.1186/1471-2164-13-132.PubMed CentralView ArticlePubMedGoogle 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 cited.