Genome-wide characterization of the biggest grass, bamboo, based on 10,608 putative full-length cDNA sequences
- Zhenhua Peng†1, 3,
- Tingting Lu†2,
- Lubin Li†1,
- Xiaohui Liu†2,
- Zhimin Gao3,
- Tao Hu1,
- Xuewen Yang3,
- Qi Feng2,
- Jianping Guan2,
- Qijun Weng2,
- Danlin Fan2,
- Chuanrang Zhu2,
- Ying Lu2,
- Bin Han2, 4Email author and
- Zehui Jiang1, 3Email author
© Peng et al; licensee BioMed Central Ltd. 2010
Received: 30 September 2009
Accepted: 18 June 2010
Published: 18 June 2010
With the availability of rice and sorghum genome sequences and ongoing efforts to sequence genomes of other cereal and energy crops, the grass family (Poaceae) has become a model system for comparative genomics and for better understanding gene and genome evolution that underlies phenotypic and ecological divergence of plants. While the genomic resources have accumulated rapidly for almost all major lineages of grasses, bamboo remains the only large subfamily of Poaceae with little genomic information available in databases, which seriously hampers our ability to take a full advantage of the wealth of grass genomic data for effective comparative studies.
Here we report the cloning and sequencing of 10,608 putative full length cDNAs (FL-cDNAs) primarily from Moso bamboo, Phyllostachys heterocycla cv. pubescens, a large woody bamboo with the highest ecological and economic values of all bamboos. This represents the third largest FL-cDNA collection to date of all plant species, and provides the first insight into the gene and genome structures of bamboos. We developed a Moso bamboo genomic resource database that so far contained the sequences of 10,608 putative FL-cDNAs and nearly 38,000 expressed sequence tags (ESTs) generated in this study.
Analysis of FL-cDNA sequences show that bamboo diverged from its close relatives such as rice, wheat, and barley through an adaptive radiation. A comparative analysis of the lignin biosynthesis pathway between bamboo and rice suggested that genes encoding caffeoyl-CoA O-methyltransferase may serve as targets for genetic manipulation of lignin content to reduce pollutants generated from bamboo pulping.
We rely on grasses more than any other groups of plants for food and potential renewable energy. With the number of genome sequences growing much more rapidly than any other plant family, the grass family (Poaceae) becomes an ideal system for comparative studies of gene and genome structure and function [1–6]. To develop the system, it is critical to accumulate genomic resources for all major lineages of the grass family. The majority of large subfamilies of grasses have already had a great deal of genomic or expressional data available primarily because they possess crop species. The only exception is the subfamily of bamboos, Bambusoideae, which contains more than 1,000 species but has little data available in DNA or protein sequence databases [4, 7]. By January of 2009, the number of ESTs deposited in the GenBank ranged from 436,535 to 2,018,337 for rice, wheat, maize, barley, sorghum, sugarcane, and switchgrass, but only 3,087 for bamboo. This creates a major missing link in the grass family for comparative genomics.
More serious than the missing link for comparative analyses is the lack of genomic resources of bamboo. This hampers biological investigations of this group of morphologically and physiologically unique and ecologically and economically important grasses. Unlike the majority of ~10,000 grass species that are herbaceous and occupy open habitats such as grassland, bamboo represents the only major lineage of grasses that lives exclusively in forests and grows large woody culms up to 30 cm in diameter and 12 m in height . In addition to remarkable sizes and woodiness, bamboo has rather striking life history characterized by a prolonged vegetative phase lasting up to more than 100 years before flowering.
With these unique features, bamboos are important components of tropical and subtropical forest ecosystems, especially in Asia, where they have had a long history of being utilized as garden ornamentals and forest products for making construction material, paper pulp, and furniture. With the realization that bamboo produces high-quality fibers and can be harvested repeatedly without severe destruction of the ecosystems, it becomes an increasingly valuable forest product that could replace a substantial portion of tree-based timber and paper pulp plantation. This highlights another important economic value of grasses in addition to food and renewable energy.
In this study, we cloned and sequenced more than ten thousand putatively unique FL-cDNA derived primarily from vegetative tissues of Moso bamboo, Phyllostachys heterocycla cv. Pubescens, a large woody bamboo with the highest ecological, economic, and cultural values of all bamboos in Asia and accounting for ~70% of total area of bamboo growth and 5 billion US dollars of annual forest production in China  (Additional file 1). The sequences reported in this study represent the third largest collection of FL-cDNA sequences of all plant species, only smaller in number than those of Arabidopsis and rice. It provides the first large sequence dataset for studying the structure and function of a substantial portion of bamboo genes, and fills the gap in the grass family for comparative genomics. A multiple-gene phylogeny inferred from these data shed light on the evolutionary relationships within the grass family. Comparative analyses of the bamboo sequence data with those of rice, barley, and wheat yielded new insights into gene evolution associated with rapid and marked phenotypic and ecological divergence between bamboo and closely related grasses. Comparison of bamboo cDNAs and rice genes involved in lignin biosynthesis suggested that genes encoding caffeoyl-CoA O-methyltransferase may serve as effective targets for genetic manipulation of lignin content to reduce pollutants generated from bamboo pulping.
Structure of bamboo cDNA sequences
From screening and sequencing five cDNA libraries, a total of 37,797 5'-end sequences comprised of at least 100 consecutive nucleotides with a Phred score above 20 were obtained and assembled into 10,669 singletons and 3,373 contigs. The average and maximal sizes of contigs were 8 and 288 cDNA clones, respectively. All singletons and one representative of each contig were fully sequenced and further filtered for redundancy, which yielded 10,608 putatively unique FL-cDNA sequences. The average length of these FL-cDNAs was 1,092 bp, and 8,695 (82%) cDNAs had open reading frames (ORFs) longer than 100 amino acids, similar to 87% in soybean  and 91% in poplar .
Mean percentage of GC content of bamboo, rice, and Arabidopsis FL-cDNAs.
Alternative splicing, natural antisense transcripts, and transposable elements
With the absence of bamboo genome sequences, we identified putative alternative splicing (AS) patterns based on six major AS types and the GT-AG rule of exon/intron boundaries (Additional file 4A). A total of 304 cDNA groups were predicted to have AS, of which 138 were putative intron retention (IntronR) or Alternative position (AltP) types, including 6 with GC-AG exon/intron boundaries, and 92, 29, 15 and 14 were putative ExonO, ExonS, AltA and AltD types, respectively (Additional file 4B). Of these, 16 cDNA groups have more than one AS patterns.
Natural antisense transcripts (NATs) overlap partially the sequences of other endogenous sense transcripts in the opposite direction. Overlapping genes involving exon regions are known as cis-encoded natural sense-antisense gene pairs [12, 13]. Of the 10,608 FL-cDNAs, 25 were identified to be sense-antisense transcripts (Additional file 5). The average overlapping length was 467 bp, with a range from 73 to 960 bp. Thirteen had predicted ORFs longer than100 amino acid and 10 had predicted ORFs between 30 and 100 amino acid, while one failed to have predicted ORFs.
BLASTX searches of the bamboo cDNAs against plant transposable element (TE) database revealed that 53 of 10,608 cDNA or 0.5% showed significant match to TEs, of which 32, 10, and 10 were homologous to polyproteins, transposases of class-I DNA elements (Mutators), and class-II RNA elements (non-long terminal repeat, non-LTR), respectively. One matched to a TE in the reverse orientation, which might be a result of genomic contamination or read-through from neighboring retrotransposons .
Functional annotation and genes involved in lignin biosynthesis
A total of 9,496 (89.5% of ) bamboo FL-cDNAs matched proteins of known or unknown functions based on BLASTX against NCBI non-redundant database (nrDB). A total of 5,404 (50.9%) cDNAs showed similarity with PFAM protein families in the InterPro database, of which those belonging to the top 20 most abundant families are shown in Additional file 6. A total of 387 (3.6%) cDNAs were predicted to be putative transcription factors according to nrDB, of which those in the top 10 most abundant classes are shown (Additional file 7). The InterPro Gene Ontology (GO) assignment identified 4,594 cDNA with GO terms associated with molecular functions in 11 categories (Additional file 8).
Number of bamboo FL-cDNAs and number of genes found in the rice genome that encode nine key enzymes in the lignin biosynthesis pathway.
4-coumarate-CoA ligase (4CL)
Caffeoyl caffeoyl-CoA O-methyltransferase (CCoAOMT)
Cinnamoyl-CoA reductase (CCR)
Caffeic acid O-methyltransferase (COMT)
Cinnamoyl alcohol dehydrogenase (CAD)
5-hydroxyconiferyl aldehyde O-methyltransferase (AldOMT)
3-deoxy-D-arabino-heptulosonate 7-phosphate synthase (DAHPS)
Comparative and phylogenetic analyses with other grasses
Sequence similarity search was conducted between bamboo and eight other grasses with the largest amount data available in TIGR Plant Transcript Assemblies. They represent three of the four large subfamilies with more than 100 species: rice (Oryza sativa) from the subfamily Ehrhartoideae with >100 species; wheat (Triticum aestivum), barley (Hordeum vulgare), and Brachypodium distachyon from Pooideae with >3,000 species; and maize (Zea mays), sorghum (Sorghum bicolor), sugarcane (Saccharum officinarum), and switchgrass (Panicum virgatum) from Panicoideae with >3,000 species. As bamboo represents the subfamily Bambusoideae with >1,000 species, Chloridoideae with >1,000 species is the only large subfamily not represented in this study due to the lack of sufficient amount of genomic data for informative comparative analysis (maximal ~20,000 ESTs for a single species).
Search for nucleotide sequence similarity with a relatively high stringency (E-value < 1e-10 in BLASTn) showed that 49.2% of bamboo cDNAs had similarity hits to rice transcripts, the highest among the grasses (Additional file 10). This compares to 81.2% and 58.1% of cDNAs that had similarity hits to rice and Arabidopsis, respectively, when the search was conducted for amino acid sequences under a less stringent condition (E-value < 1e-6 in tBLASTx) . The higher proportion of bamboo cDNAs had hits to Arabidopsis than Branchypodium and switchgrass was due to the relatively small number of transcripts available for these two grasses. Similarity search was also conducted between the bamboo FL-cDNAs and the genome sequences of rice, sorghum, Arabidopsis, and poplar with both nucleotide and amino acid sequences under different stringencies. The proportions of similarity hits were somewhat lower than those hit the transcript database (Additional file 11).
We then compared the functional classification between two groups of bamboo cDNAs, one including shared homologs with other two subfamilies and the other being unique to bamboo (Figure 3B). It is noteworthy that shared cDNAs are relatively abundant in catalytic molecules, structural molecules, and transporters. In contrast, cDNAs unique to bamboo are relatively abundant in groups with binding activities and regulatory functions.
Databasing FL-cDNAs and ESTs of Moso bamboo and EST redundancy analysis
To facilitate the access and utilization of the bamboo cDNA sequences, we developed a Moso bamboo cDNA database (MBCD), including 10,608 putative FL-cDNAs and nearly 38,000 ESTs. The database is hosted on a REDHAT LINUX server accessible at the URL http://www.ncgr.ac.cn/MBCD/ (Additional file 15). The source code is available at http://www.ncgr.ac.cn/mbcd/code/ or http://www.ncgr.ac.cn/mbcd/res/down.php The database has the following features: 1) FL-cDNA clone list search. Given cDNA clone ID or NCBI accession number, MBCD will display the cDNA sequences with predicted ORF information, library (tissue) sources, clone ID, Genbank accession numbers, and corresponding ESTs; 2) Sequence search. Blast is provided for sequence alignment search against the cDNA database. Conserved domain search is also available when putative function information (Pfam) is entered; 3) FTP download links. The cDNA sequences, annotations, and corresponding ESTs can be downloaded as compressed files of gzip or bzip; and 4) Other useful web links. Other database websites and research institutes related to bamboo research are conveniently accessible through this website.
Taking advantage of the database, we took a snapshot on the genome redundancy of Moso bamboo by identifying and analyzing EST clusters. BLASTN searches of all ESTs of Moso bamboo were conducted against the FL-cDNAs. We used stringent criteria to identify EST clusters that contained at least 6 homologous ESTs with each sequence polymorphism supported by at least two ESTs. A total of 2,728 clusters were identified, which included 21,497 ESTs.
Bamboo genes and genome
In comparison with a large number of plant ESTs deposited in public database (NCBI ESTs of plants >100 millions), only a small fraction of high-quality sequence-finished FL-cDNAs have been generated. Prior to this study, rice and Arabidopsis were the only plant species with more than 10,000 fully sequenced FL-cDNAs available in the public databases [17–19]. Thus, here we report the third largest FL-cDNA collection of a plant species. Given the estimate that the fully-sequenced genomes of rice and sorghum had approximately 35,000 - 38,000 genes [20, 21], the 10,608 putative FL-cDNAs are likely to represent more than a quarter of bamboo genes. One should note this cDNA collection was from various vegetative tissues while the floral tissues were not available for this study because Moso bamboo flowers every 67 years on average. Nevertheless, these data provide a valuable resource for the first glance at bamboo gene diversity and for comparative genomic studies among grasses.
GC content is an important characteristic of genomic composition, reflecting various genome features, such as gene density and structure, codon usage, repeat element distribution, and methylation patterns [15, 22, 23]. Here we found that bamboo had the highest overall GC content as well as the highest GC content at the third codon position in comparison to rice and Arabidopsis (Table 1; [24, 25]). Further, the third codon position of bamboo and rice both displays the highest GC content, followed by the first and then the second positions, whereas Arabidopsis has the highest GC content at the first codon position. With regard to the distribution of GC content in the 5' UTR, ORF, and 3' UTR, we found similarities between the two grasses that differed from Arabidopsis (Additional file 2). The distributions of GC content in the 5' UTR and ORF are rather flat or slightly bimodal for bamboo and rice, but clearly unimodal for Arabidopsis.
The codon usages of bamboo and rice also share considerable similarities in comparison to that of Arabidopsis (Additional file 2). Given the much higher GC content of bamboo and rice cDNAs, it is not surprising that they use G and C more frequently whereas A and T appear at a higher frequency for the same codons in Arabidopsis. There has been evidence that high GC content and corresponding codon usage bias of rice might have influenced the level of gene expression and the mechanisms of gene evolution . The even higher GC content found in Bamboo could have played a role in the genetic and genomic divergence of the large, woody grass.
The type and distribution of SSRs of bamboo are also more similar to those of rice than Arabidopsis. The tri-nucleotide SSRs are more abundant than mono- and di-nucleotide SSRs in bamboo and rice, whereas mono-nucleotide SSRs are the most abundant in Arabidopsis (Figure 1). Of tri-nucleotide SSRs of bamboo and rice that occur predominantly in ORFs, the CCG/CGG type is much more abundant than others. This is perhaps correlated with the higher GC content of the grasses, and may have allowed frequently insertion/deletion of certain amino acids without causing frame shift mutations . Interestingly, bamboo has higher frequency of all three major SSRs in 5' UTRs from the start codon to 150 bp upstream than rice and Arabidopsis, implying that SSRs may have played an important role in regulating gene expression in bamboo.
Alternative splicing may be involved in posttranscriptional regulations that increase protein diversity and affect mRNA stability [28, 29]. According to the previous classification of AS , we found two types of AS, IntronR and single exon overlapping (sExonO), were more abundant than others. We calculated the length between two continuous exons of each cDNA pair in the IntronR type, which gave an estimate of mean intron size of 218 bp, comparing to 461 bp and 158 bp of the mean intron sizes of rice and Arabidopsis, respectively [19, 31].
After the BLASTX searching against NCBI nrDB, 672 bamboo FL-cDNAs had no hit in the database, some of which may be new genes. Functional annotation of transcription factors indicated that ERF was the most abundant among bamboo cDNAs, which was followed by MYB, Zinc finger, WRKY, homeobox, bZIP, bHLH, and NAC. This differs from what was found in a rice cDNA collection where Zinc finger was much more abundant than other transcription factors, whereas ERF was ranked as the third abundant, fewer than Myb . The difference could be due to differential amplification of various families of transcription factors during the evolutionary divergence between bamboo and rice. But it is possible that at least a part of the difference is due to an artifact of cDNA isolation from different tissues at different developmental stages between the two studies.
Because we filtered gene redundancy by selecting one cDNA clone from a contig for full-length sequencing, the sequences of FL cDNAs could not be used to evaluate the genome redundancy. Because the Moso bamboo was considered to be a tetraploid , we went back to analyze the ESTs to take a snapshot at the genome redundancy. From a total of 100 sampled EST clusters, the divergence times estimated between two homologous EST types within clusters had an interesting distribution (Figure 5). The abrupt increase in redundancy around 15 million years ago (mya) or more recent is probably a result of genome duplication through polyploid formation. With the EST data, we were unable to determine whether this was an autopolyploidy or allopolyploidy event. In any case, it seems like that the two diploid genomes have diverged for less than 15 million years according to this estimation. The actual divergence time, however, could be longer given that the Moso Bamboo has a longer generation time than most of grasses from which the synonymous substitution rate was calculated.
Genes involved in lignin biosynthesis
The lignin content of bamboo is comparable to that of woody plants and higher than most herbaceous plants, which has contributed to high culm rigidity that allows bamboo to grow into a large forest . The lignin content of mature culms of Moso bamboo, estimated at about 23% of dry biomass , is nearly twice as high as that of rice straw (12.5%, ). The high lignin content, however, has a negative impact on paper production from bamboos because removal of lignin during pulping results in hazardous waste. A better under understanding of lignin biosynthesis is important for manipulation lignin content or composition through genetic engineering.
On average, bamboo cDNAs for lignin biosynthesis enzymes isolated in this study are one fifth of the total number of rice genes encoding lignin biosynthesis enzymes. This is not surprising because the cDNAs may represent only one third to one fourth of bamboo genes. From these ratios as well as the phylogenetic inference, it seems likely that there has not been a large-scale duplication of lignin biosynthesis genes in bamboo after it diverged from rice (Additional file 9). Whereas gene duplication might have played a role in altering certain physiological properties, such as drought tolerance in sorghum , it was probably not a major reason for the increased lignin content of bamboo in comparison to rice.
Nevertheless, it is noteworthy that the number of cDNAs found for caffeoyl-CoA O-methyltransferase (CCoAOMT) was nearly the same as that of rice genes. This much higher than average ratio could be a result of duplication of genes encoding this particular enzyme in bamboo. It is also possible that the majority of genes for CCoAOMT had high expression levels so that a larger proportion of the cDNA population was captured for sequencing. Either case could support an important role of CCoAOMT in lignin biosynthesis of Moso bamboo. Interestingly, CCoAOMT was identified as key enzyme that determined lignin content in woody trees such as poplar . Furthermore, the repression of CCoAOMT by antisense led to significant reduction of lignin content of transgenic popular but did not affect normal growth of the transgenic plants . These findings together suggested that CCoAOMT may serve as an effective target for genetic manipulation to reduce pollutants generated from bamboo pulping.
Phylogeny and evolutionary divergence of bamboo from other grasses
The cloning and sequencing of more than ten thousand bamboo putative FL-cDNAs generated the largest set of unlinked nuclear loci to date for phylogenetic analysis of major grass lineages. Even with 43 sequences, the relationship of bamboo was not unambiguously resolved. Different methods of phylogenetic analyses disagreed whether bamboo was more closely related to rice or barley and wheat. While the disagreement was relatively well supported in the respective phylogenies, the branch leading to either grouping was short (Figure 2), suggesting that the three subfamilies, Bambusoideae, Pooideae, and Ehrhartoideae, diverged rapidly from each other, possibly through a process known as adaptive radiation.
In the similar cases where adaptive radiation was involved, lineage sorting posed a major problem to correct reconstruction of the species phylogeny from individual gene sequences, and as a consequence a large number of independent molecular markers were required to resolve the phylogenetic relationships . In the previous phylogenetic studies of the grass family based on only a few chloroplast and nuclear genes, there seemed to have been a consensus that Bambusoideae and Ehrhartoideae were sister groups, i.e., a closer relationship between bamboo and rice [8, 38], until a recent study suggesting that Bamboideae and Pooideae (wheat and barley) were more closely related . Our analysis of the 43 putative orthologous genes demonstrates that the relationship among these three subfamilies remains unresolved and this is mostly likely due to rapid diversification of these grasses. The finding raises an intriguing question of what evolutionary forces drove the rapid diversification of bamboo from other grasses and its adaptation to a very different habitat.
Although a direct answer to this question requires long-term research using integrative approaches, some insights may be gained by analyzing the patterns and rates of gene evolution between bamboo and its close relatives. We found that the number of FL-cDNAs with and without homologs in rice, wheat, and barley is nearly equal (Figure 3A). Interestingly, the shared genes encode mostly enzymes, transporters, and structure proteins, whereas genes unique to bamboo encode more frequently binding and regulatory factors. This suggests that regulatory factors evolved more rapidly than enzymes and structural proteins following the divergence of bamboo from other grasses. It has been widely observed that rapid phenotypic evolution of plants at population and species levels was controlled primarily by regulatory genes [40–42]. The characteristics of gene evolution between bamboo and related subfamilies may represent a signature of this evolutionary mechanism at a higher taxonomic level. We thus hypothesize that rapid evolution of regulatory genes could have played an important role during the rapid adaptive diversification of bamboo from other grasses.
For nearly 6,000 bamboo cDNAs with homologs in rice, there are on ly a few pairs that had a K A /K S ratio larger than 1, a threshold measuring positive selection. The lack of gene with high K A /K S values is likely a result of sufficiently long divergence time between bamboo and rice that has obscured the signal of positive selection driving adaptation. When the K A /K S ratio was compared in a finer scale, we found that genes encoding enzymes, transporters, and structural proteins were more abundant in the categories with the K A /K S ratio smaller than 0.2. In contrast, a larger number of genes encoding binding and regulatory factors had the K A /K S ratio between 0.2 and 0.3 (Additional file 14B). This is consistent with the finding when the overall level of divergence was considered, together indicating that regulatory factors not only evolved more rapidly but were more likely to have evolved in response to natural selection driving rapid adaptive divergence of these grasses. It was probably the rapid divergence of regulatory genes that allowed drastic morphological and physiological evolution of bamboos and rapid adaptation to the new forest habitats.
This study generated the first large collection of bamboo FL-cDNAs and ESTs and the first genomic resource database for bamboos. This is the initial effort to sequence and analyze the genome of Moso bamboo, which has gained increasingly important ecological and economical values. The analyses of yielded putative FL-cDNAs provided the first glance at the structural and functional features of bamboo genomes. Phylogenetic analyses of the bamboo sequences with those of rice, barley, and wheat yielded suggested that bamboo has diverged from these morphologically and physiologically different relatives possibly through an adaptive radiation. Comparative analyses of bamboo cDNAs and rice genes involved in lignin biosynthesis indicated that genes encoding caffeoyl-CoA O-methyltransferase may serve as effective targets for genetic manipulation of lignin content to reduce pollutants resulting from bamboo pulping. The sequences of the FL-cDNAs and ESTs generated in this study will close a critical gap existing in grass comparative genomics and consequently allow the more efficient development of the grass system for evolutionary and functional studies of plant genes and genomes.
As similar study in other plants, we here provided an enriched putative full-length cDNA data collection and corresponding amino acid sequence information as a start of genomic research. More precise transcripts would be confirmed in coming gene cloning and functional research in bamboo.
Plant materials and cDNA library construction
Five cDNA libraries were constructed for Phyllostachys heterocycla cv. pubescens from shoots just breaking out from ground, shoots reaching a height of ~40 cm, young leaves, and shoots and roots removed from germinating seeds (Additional file 1). For the first three libraries, shoots and leaves were collected in April from plants naturally grown in the Tianmu-Mountain National Nature Reserve in Zhejiang province of China and immediately preserved in liquid nitrogen and then stored at -80°C in the lab until RNA extraction. Shoots were dug out, and non-lignified tissues at the top of shoots with their roots were sampled. From germinating seeds, 1-2 cm shoots and roots were cut and pooled for library construction. Another cDNA library was constructed from leaves of Bashania fangiana, collected in May from the Wolong National Nature Reserve in Sichuan province of China. A Cap-Tagging method was used to construct putative full-length enriched cDNA libraries .
Sequencing and assembling FL-cDNA clones
Approximately 50,000 cDNA clones were randomly selected from the five libraries for 5' end single-pass sequencing using BigDye Terminator Cycle sequencing V2.0 Ready Reaction (Applied Biosystems). PHRED and PHRAP [43, 44] were used to generate and assemble raw data. Vector sequences were filtered automatically and low-quality bases (quality value <20) were removed. All 5' tagged ESTs were clustered into contigs or singletons using TIGR Gene Indices clustering tools [45, 46]. If sequences share >95% similarity over 80 consecutive bases, they were clustered into one contig. All singletons and one randomly chosen clone from each contig were fully sequenced in both directions, with internal sequencing primers designed for clones longer than 1 kb. To minimize sequence errors, all assembled sequences were manually checked. Programs CAP3 and BLAST were then used to filter redundant sequences judged by having 99% or higher identity .
ORFs were predicted using the "getorf" program of EMBOSS package , with the longest ORF extracted for each FL-cDNA. ORFs longer than 100 amino acids were analyzed further for codon usage and SSR using the CUSP program implemented in EMBOSS and a perl script, MISA http://pgrc.ipk-gatersleben.de/misa/, respectively. For comparison, cDNA data of Arabidopsis and rice were downloaded from public databases PlantGDB http://www.plantgdb.org/AtGDB/ and NCBI ftp://ftp.ncbi.nih.gov, respectively.
To identify putative alternative splicing variants, program CAP3  was first used to filter singletons with default parameters. The remaining sequences assembled to contigs were subjected to a BLASTN analysis against themselves. Those with sequence identity higher than 99% for at least 100 bp were extracted. The GT-AG rule of exon/intron boundary was then applied using program BLAST2 followed by manual check. To identify putative sense-antisense transcript pairs, all sequences were searched against themselves using BLASTN. Those with plus and minus matching pairs that were at least 50 nucleotide long and had less than 2 bp mismatches were recognized.
A database of plant transposable element peptide sequences was retrieved from the GenBank non-redundant protein database using keywords, including transpose, retro, non-LTR, en/spm, ac/ds, gypsy, copia, polyprotein, mutator, and mudr [24, 49]. The bamboo FL-cDNA sequences were then searched against this database using BLASTx with E-value lower than 1e-20.
All bamboo putative full-length cDNA sequences were searched against NCBI nrDB using BLASTX (E-value < 1e-10). The InterPro database was used for identifying putative protein domains . Functional classification followed PFAM  and GO terms attached to the InterPro domain names.
Comparative and phylogenetic analyses with other grasses
Similarity searches were performed with BLASTn and tBLASTx (version 2.2.14)  against the public database TIGR Plant Transcript Assemblies http://plantta.tigr.org/; ) (including Oryza sativa, Triticum aestivum, Zea mays, Zea mays B73, Hordeum vulgare, Sorghum bicolor, Saccharum officinarum, Brachypodium distachyon, and Panicum virgatum, and Arabidopsis thaliana), with an E-value cutoff of 10-10 (BLASTn) and 10-6 (tBLASTx), respectively. The bamboo cDNAs were also compared with currently available genome sequences of Arabidopsis http://www.plantgdb.org/AtGDB/, rice (IRGSP version 4.0: http://rgp.dna.affrc.go.jp/IRGSP/, Sorghum bicolor http://www.gramene.org/Sorghum_bicolor/index.html, and Populus trichocarpa http://www.ornl.gov/sci/ipgc/ with BLASTn and tBLASTn, with the E-value cutoff of 10-1 (BLASTn) and 10-7 (tBLASTn), respectively. Finally, the bamboo cDNAs were searched against NCBI ntDB ftp://ftp.ncbi.nih.gov using BLASTN (E-value < 1e-10). For nucleotide similarity search, homologs were recognized as with >75% sequence identity for >50% bamboo FL-cDNA sequences. For tBLASTn and tBLASTx search of amino acid sequences, homologs were recognized as with >60% sequence identity and for >50 consecutive amino acids.
To identify putative orthologs for phylogenetic analyses, we applied a stringent criterion requiring that there was only one homolog identified in any other species when a bamboo cDNA was used as the query in the BLAST search and vice versa. Identification of orthologs in other species was carried out in three steps. 1) Three hits with the best E-value in BLAST were selected. Their DNA sequences were re retrieved from public database by their accession number. 2) The selected sequences were put into ORF prediction again because we found some of them do not carry entire ORF. 3) The predicted amino acid sequences were aligned to the corresponding sequence from bamboo. Only ONE sequence with the highest identity in alignment was considered as a putative ortholog.
For phylogenetic analysis, DNA sequences were aligned using program ClustalW with the default settings . Phylogenetic trees were inferred by maximum likelihood (ML), Neighbour joining (NJ), and Bayesian inference (BI). ML and NJ were implemented with PAUP 4.0b10  and the branch-and-bound algorithm was used for tree searching. BI was employed for phylogenetic reconstruction using the MrBayes software of version 3.1 . For sequence divergence, synonymous (K S ) and nonsynonymous (K A ) distances were estimated using program PAL2NAL .
Analysis of EST clusters for gene/genome duplication
ESTs of Moso bamboo were aligned with the putative FL-cDNAs using the ungapped BLASTN search at the highest E-value of 0.0. For those FL-cDNAs with at least six EST counter parts, the corresponding ESTs were aligned with CLUSTAL 2.0.11 http://software.informer.com/getfree-clustalx2-cite/; ). For the resulting EST clusters, we retained those with any given sequence polymorphism supported by at least two members of the ESTs to avoid artifacts of sequence errors. Only those EST clusters that covered the entire predicted coding regions of the FL-cDNA were sampled for estimating sequence divergence. Synonymous substitutions (K S ) between homologs were calculated [59, 60]. Divergence time were estimated based on T = K S /2r using the divergence rate r = 6.5 × 10-9 mutations per synonymous site per year .
Databasing ESTs and FL-cDNAs of Moso bamboo
To facilitate the access and utilization of putative FL-cDNA and EST sequences generated for Moso bamboo, we developed a database to store and search these sequences and their annotations. The web interface of this database was constructed in PHP script and the data were stored in a relational database management system, MYSQL. Search function and display were built with a combination of SQL commands and PHP script.
Sequence data from this article can be found in the GenBank under accession numbers FP091249-FP101855.
We thank Xinhui Zou and Xuehui Huang for assistance with phylogenetic and statistical analyses. We thank Tao Sang for critical reading the manuscript and correction. This work was supported by the Forestry Project of the Ministry of Science and Technology of the People's Republic of China (Grant No. 200704001) and by the Chinese Academy of Sciences.
- Gale MD, Devos KM: Comparative genetics in the grasses. Proc Natl Acad Sci USA. 1998, 95: 1971-1974. 10.1073/pnas.95.5.1971.PubMedPubMed CentralView ArticleGoogle Scholar
- Gaut BS: Evolutionary dynamics of grass genomes. New Phytologist. 2001, 154: 15-28. 10.1046/j.1469-8137.2002.00352.x.View ArticleGoogle Scholar
- Bennetzen JL: Patterns in grass genome evolution. Curr Opin in Plant Biol. 2007, 10: 176-181. 10.1016/j.pbi.2007.01.010.View ArticleGoogle Scholar
- Buell CR: Poaceae genomes: Going from unattainable to becoming a model clade for comparative plant genomics. Plant Physiol. 2009, 149: 111-116. 10.1104/pp.108.128926.PubMedPubMed CentralView ArticleGoogle Scholar
- Messing J: Synergy of two reference genomes for the grass family. Plant Physiol. 2009, 149: 117-124. 10.1104/pp.108.128520.PubMedPubMed CentralView ArticleGoogle Scholar
- Paterson AH, Bowers JE, Feltus FA, Tang H, Lin L, Wang X: Comparative genomics of grasses promises a bountiful harvest. Plant Physiol. 2009, 49: 125-131. 10.1104/pp.108.129262.View ArticleGoogle Scholar
- Das M, Bhattacharya S, Singh P, Filgueiras T, Pal A: Bamboo taxonomy and diversity in the era of molecular markers. Adv Bot Res. 2008, 47: 225-267. 10.1016/S0065-2296(08)00005-0.View ArticleGoogle Scholar
- Barker NP, Clark LG, Davis JI, Duvall MR, Guala GF, Hsiao C, Kellogg EA, Linder PH: Phylogeny and subfamilial classification of the grasses (Poaceae). Ann Missouri Bot Garden. 2001, 88: 373-457. 10.2307/3298585.View ArticleGoogle Scholar
- Jiang ZH: Bamboo and rattan in the world. LiaoNing Science and Technology Published House. 2002Google Scholar
- Umezawa T, Sakurai T, Totoki Y, Toyoda A, Seki M, Ishiwata A, Akiyama K, Kurotani A, Yoshida T, Mochida K, Kasuga M, Todaka D, Maruyama K, Nakashima K, Enju A, Mizukado S, Ahmed S, Yoshiwara K, Harada K, Tsubokura Y, Hayashi M, Sato S, Anai T, Ishimoto M, Funatsuki H, Teraishi M, Osaki M, Shinano T, Akashi R, Sakaki Y, Yamaguchi-Shinozaki K, Shinozaki K: Sequencing and analysis of approximately 40,000 soybean cDNA clones from a full-length-enriched cDNA library. DNA Res. 2008, 15: 333-346. 10.1093/dnares/dsn024.PubMedPubMed CentralView ArticleGoogle Scholar
- Ralph SG, Chun HJ, Cooper D, Kirkpatrick R, Kolosova N, Gunter L, Tuskan GA, Douglas CJ, Holt RA, Jones SJ, Marra MA, Bohlmann J: Analysis of 4,664 high-quality sequence-finished poplar full-length cDNA clones and their utility for the discovery of genes responding to insect feeding. BMC Genomics. 2008, 9: 57-10.1186/1471-2164-9-57.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang Y, Liu XS, Liu QR, Wei L: Genome-wide in silico identification and analysis of cis natural antisense transcripts (cis-NATs) in ten species. Nucleic Acids Res. 2006, 34: 3465-3475. 10.1093/nar/gkl473.PubMedPubMed CentralView ArticleGoogle Scholar
- Chen J, Sun M, Kent WJ, Huang X, Xie H, Wang W, Zhou G, Shi RZ, Rowley JD: Over 20% of human transcripts might form sense-antisense pairs. Nucleic Acids Res. 2004, 32: 4812-4820. 10.1093/nar/gkh818.PubMedPubMed CentralView ArticleGoogle Scholar
- Elrouby N, Bureau TE: A novel hybrid open reading frame formed by multiple cellular gene transductions by a plant long terminal repeat retroelement. J Biol Chem. 2001, 276: 41963-41968. 10.1074/jbc.M105850200.PubMedView ArticleGoogle Scholar
- Yu J, Hu S, Wang J, Wong GK, Li S, Liu B, Deng Y, Dai L, Zhou Y, Zhang X, Cao M, Liu J, Sun J, Tang J, Chen Y, Huang X, Lin W, Ye C, Tong W, Cong L, Geng J, Han Y, Li L, Li W, Hu G, Huang X, Li W, Li J, Liu Z, Li L, Liu J, Qi Q, Liu J, Li L, Li T, Wang X, Lu H, Wu T, Zhu M, Ni P, Han H, Dong W, Ren X, Feng X, Cui P, Li X, Wang H, Xu X, Zhai W, Xu Z, Zhang J, He S, Zhang J, Xu J, Zhang K, Zheng X, Dong J, Zeng W, Tao L, Ye J, Tan J, Ren X, Chen X, He J, Liu D, Tian W, Tian C, Xia H, Bao Q, Li G, Gao H, Cao T, Wang J, Zhao W, Li P, Chen W, Wang X, Zhang Y, Hu J, Wang J, Liu S, Yang J, Zhang G, Xiong Y, Li Z, Mao L, Zhou C, Zhu Z, Chen R, Hao B, Zheng W, Chen S, Guo W, Li G, Liu S, Tao M, Wang J, Zhu L, Yuan L, Yang H: A draft sequence of the rice genome (Oryza sativa L. ssp. indica). Science. 2002, 296: 79-92. 10.1126/science.1068037.PubMedView ArticleGoogle Scholar
- Gaut BS, Morton BR, McCaig BC, Clegg MT: Substitution rate comparisons between grasses and palms: synonymous rate differences at the nuclear gene Adh parallel rate differences at the plastid gene rbcL. Proc Natl Acad Sci USA. 1996, 19: 10274-10279. 10.1073/pnas.93.19.10274.View ArticleGoogle Scholar
- Seki M, Narusaka M, Kamiya A, Ishida J, Satou M, Sakurai T, Nakajima M, Enju A, Akiyama K, Oono Y, Muramatsu M, Hayashizaki Y, Kawai J, Carninci P, Itoh M, Ishii Y, Arakawa T, Shibata K, Shinagawa A, Shinozaki K: Functional annotation of a full-length Arabidopsis cDNA collection. Science. 2002, 296: 141-145. 10.1126/science.1071006.PubMedView ArticleGoogle Scholar
- The Rice Full-Length cDNA Consortium: Collection, mapping, and annotation of over 28,000 cDNA clones from japonica rice. Science. 2003, 301: 376-379. 10.1126/science.1081288.View ArticleGoogle Scholar
- Liu X, Lu T, Yu S, Li Y, Huang Y, Huang T, Zhang L, Zhu J, Zhao Q, Fan D, Mu J, Shangguan Y, Feng Q, Guan J, Ying K, Zhang Y, Lin Z, Sun Z, Qian Q, Lu Y, Han B: A collection of 10,096 indica rice full-length cDNAs reveals highly expressed sequence divergence between Oryza sativa indica and japonica subspecies. Plant Mol Biol. 2007, 65: 403-415. 10.1007/s11103-007-9174-7.PubMedView ArticleGoogle Scholar
- International Rice Genome Sequencing Project: The map-based sequence of the rice genome. Nature. 2005, 436: 793-800. 10.1038/nature03895.View ArticleGoogle Scholar
- Paterson AH, Bowers JE, Bruggmann R, Dubchak I, Grimwood J, Gundlach H, Haberer G, Hellsten U, Mitros T, Poliakov A, Schmutz J, Spannagl M, Tang H, Wang X, Wicker T, Bharti AK, Chapman J, Feltus FA, Gowik U, Grigoriev IV, Lyons E, Maher CA, Martis M, Narechania A, Otillar RP, Penning BW, Salamov AA, Wang Y, Zhang L, Carpita NC, Freeling M, Gingle AR, Hash CT, Keller B, Klein P, Kresovich S, McCann MC, Ming R, Peterson DG, Mehboob-ur-Rahman Ware D, Westhoff P, Mayer KF, Messing J, Rokhsar DS: The Sorghum bicolor genome and the diversification of grasses. Nature. 2009, 457: 551-556. 10.1038/nature07723.PubMedView ArticleGoogle Scholar
- Galtier N, Piganeau G, Mouchiroud D, Duret L: GC- content evolution in mammalian genomes: the biased gene conversion hypothesis. Genetics. 2001, 159: 907-911.PubMedPubMed CentralGoogle Scholar
- Ingvarsson PK: Molecular evolution of synonymous codon usage in Populus. BMC Evol Biol. 2008, 8: 307-10.1186/1471-2148-8-307.PubMedPubMed CentralView ArticleGoogle Scholar
- Kuhl JC, Cheung F, Yuan Q, Martin W, Zewdie Y, McCallum J, Catanach A, Rutherford P, Sink KC, Jenderek M, Prince JP, Town CD, Havey MJ: A unique set of 11,008 onion expressed sequence tags reveals expressed sequence and genomic differences between the monocot orders Asparagales and Poales. Plant Cell. 2004, 16: 114-125. 10.1105/tpc.017202.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang HC: Hickey DARapid divergence of codon usage patterns within the rice genome. BMC Evol Biol. 2007, 7: S6-10.1186/1471-2148-7-S1-S6.PubMedPubMed CentralView ArticleGoogle Scholar
- Guo X, Bao J, Fan L: Evidence of selectively driven codon usage in rice. Implications for GC content evolution of Gramineae genes. FEBS Letters. 2007, 581: 1015-1021. 10.1016/j.febslet.2007.01.088.PubMedView ArticleGoogle Scholar
- Lawson MJ, Zhang L: Distinct patterns of SSR distribution in the Arabidopsis thaliana and rice genomes. Genome Biol. 2006, 7: R14-10.1186/gb-2006-7-2-r14.PubMedPubMed CentralView ArticleGoogle Scholar
- Lareau LF, Green RE, Bhatnagar RS, Brenner SE: The evolving roles of alternative splicing. Curr Opin Struct Biol. 2004, 14: 273-282. 10.1016/j.sbi.2004.05.002.PubMedView ArticleGoogle Scholar
- Stamm S, Ben-Ari S, Rafalska I, Tang Y, Zhang Z, Toiber D, Thanaraj TA, Soreq H: Function of alternative splicing. Gene. 2005, 344: 1-20. 10.1016/j.gene.2004.10.022.PubMedView ArticleGoogle Scholar
- Wang BB, Brendel V: Genomewide comparative analysis of alternative splicing in plants. Proc Natl Acad Sci USA. 2006, 103: 7175-7180. 10.1073/pnas.0602039103.PubMedPubMed CentralView ArticleGoogle Scholar
- Hong X, Scofield DG, Lynch M: Intron Size, Abundance, and Distribution within Untranslated Regions of Genes. Mol Biol Evol. 2006, 23: 2392-2404. 10.1093/molbev/msl111.PubMedView ArticleGoogle Scholar
- Gui Y, Wang S, Quan L, Zhou C, Long S, Zheng H, Jin L, Zhang X, Ma N, Fan L: Genome size and sequence composition of moso bamboo: a comparative study. Sci China C Life Sci. 2007, 50: 700-705. 10.1007/s11427-007-0081-6.PubMedView ArticleGoogle Scholar
- Scurlock JMO, Dayton DC, Hames B: Bamboo: an overlooked biomass resource?. Biomass Bioenergy. 2000, 19: 229-244. 10.1016/S0961-9534(00)00038-6.View ArticleGoogle Scholar
- Li XB, Shupe TF, Peter GF, Hse CY, Eberhardt TL: Chemical changes with maturation of the bamboo species Phyllostrachys pubescens. J Trop Forest Sci. 2007, 19: 6-12.Google Scholar
- Pauly M, Keegstra K: Cell-wall carbohydrates and their modification as a resource for biofuels. Plant J. 2008, 54: 559-568. 10.1111/j.1365-313X.2008.03463.x.PubMedView ArticleGoogle Scholar
- Zhong R, Morrison WH, Himmelsbach DS, Poole FL, Ye ZH: Essential role of caffeoyl coenzyme A O-methyltransferase in lignin biosynthesis in woody poplar plants. Plant Physiol. 2000, 124: 563-578. 10.1104/pp.124.2.563.PubMedPubMed CentralView ArticleGoogle Scholar
- Zou XH, Zhang FM, Zhang JG, Zang LL, Tang L, Wang J, Sang T, Ge S: Analysis of 142 genes resolves the rapid diversification of the rice genus. Genome Biol. 2008, 9: R49-10.1186/gb-2008-9-3-r49.PubMedPubMed CentralView ArticleGoogle Scholar
- Kellogg EA: Evolutionary history of the grasses. Plant Physiol. 2001, 125: 1198-1205. 10.1104/pp.125.3.1198.PubMedPubMed CentralView ArticleGoogle Scholar
- Bouchenak-Khelladi Y, Salamin N, Savolainen V, Forest F, Van Der Bank M, Chase MW, Hodkinson TR: Large multi-gene phylogenetic trees of the grasses (Poaceae): Progress towards complete tribal and generic level sampling. Mol Phylogent Evl. 2008, 47: 488-505. 10.1016/j.ympev.2008.01.035.View ArticleGoogle Scholar
- Doebley JF, Lukens L: Transcriptional regulators and the evolution of plant form. Plant Cell. 1998, 10: 1075-1082. 10.1105/tpc.10.7.1075.PubMedPubMed CentralView ArticleGoogle Scholar
- Doebley JF, Gaut BS, Smith BD: The molecular genetics of crop domestication. Cell. 2006, 127: 1309-1321. 10.1016/j.cell.2006.12.006.PubMedView ArticleGoogle Scholar
- Kellogg EA: Progress and challenges in studies of the evolution of development. J Exp Bot. 2006, 57: 3505-3516. 10.1093/jxb/erl132.PubMedView ArticleGoogle Scholar
- Ewing B, Green P: Base-calling of automated sequencer traces using PHRED. II. Error probabilities. Genome Res. 1998, 8: 186-194.PubMedView ArticleGoogle Scholar
- Ewing B, Hillier L, Wendl MC, Green P: Base-calling of automated sequencer traces using PHRED. I. Accuracy assessment. Genome Res. 1998, 8: 175-185.PubMedView ArticleGoogle Scholar
- Chou HH, Holmes MH: DNA sequence quality trimming and vector removal. Bioinformatics. 2001, 17: 1093-1104. 10.1093/bioinformatics/17.12.1093.PubMedView ArticleGoogle Scholar
- Pertea G, Huang X, Liang F, Antonescu V, Sultana R, Karamycheva S, Lee Y, White J, Cheung F, Parvizi B, Tsai J, Quackenbush J: TIGR gene indices clustering tools (TGICL): a software system for fast clustering of large EST datasets. Bioinformatics. 2003, 19: 651-652. 10.1093/bioinformatics/btg034.PubMedView ArticleGoogle Scholar
- Huang X, Madan A: CAP3: A DNA sequence assembly program. Genome Res. 1999, 9: 868-877. 10.1101/gr.9.9.868.PubMedPubMed CentralView ArticleGoogle Scholar
- Rice P, Longden I, Bleasby A: EMBOSS: the European molecular biology open software suite. Trends Genet. 2000, 16: 276-277. 10.1016/S0168-9525(00)02024-2.PubMedView ArticleGoogle Scholar
- Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Wheeler DL: GenBank. Nucleic Acids Res. 2003, 31: 23-27. 10.1093/nar/gkg057.PubMedPubMed CentralView ArticleGoogle Scholar
- Apweiler R, Attwood TK, Bairoch A, Bateman A, Birney E, Biswas M, Bucher P, Cerutti L, Corpet F, Croning MD, Durbin R, Falquet L, Fleischmann W, Gouzy J, Hermjakob H, Hulo N, Jonassen I, Kahn D, Kanapin A, Karavidopoulou Y, Lopez R, Marx B, Mulder NJ, Oinn TM, Pagni M, Servant F, Sigrist CJ, Zdobnov EM: The InterPro database, an integrated documentation resource for protein families, domains and functional sites. Nucleic Acids Res. 2001, 29: 37-41. 10.1093/nar/29.1.37.PubMedPubMed CentralView ArticleGoogle Scholar
- Bateman A, Coin L, Durbin R, Finn RD, Hollich V, Griffiths-Jones S, Khanna A, Marshall M, Moxon S, Sonnhammer EL, Studholme DJ, Yeats C, Eddy SR: The Pfam protein families database. Nucleic Acids Res. 2004, 32: 138-141. 10.1093/nar/gkh121.View ArticleGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang JH, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search Programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMedPubMed CentralView ArticleGoogle Scholar
- Childs KL, Hamilton JP, Zhu W, Ly E, Cheung F, Wu H, Rabinowicz PD, Town CD, Buell CR, Chan AP: The TIGR Plant Transcript Assemblies database. Nucleic Acids Res. 2007, 35: D846-D851. 10.1093/nar/gkl785.PubMedPubMed CentralView ArticleGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: Improving the sensitivity of progressive multiple sequence alignment through sequence weighting position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.PubMedPubMed CentralView ArticleGoogle Scholar
- Swofford DL: PAUP*: Phylogenetic Analysis Using Parsimony (*and Other Methods). Sinauer, Sunderland, MA; 2002.Google Scholar
- Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogeny. Bioinformatics. 2001, 17: 754-755. 10.1093/bioinformatics/17.8.754.PubMedView ArticleGoogle Scholar
- Suyama M, Torrents D, Bork P: PAL2NAL: robust conversion of protein sequence alignments into the corresponding codon alignments. Nucleic Acids Res. 2006, 34: W609-W612. 10.1093/nar/gkl315.PubMedPubMed CentralView ArticleGoogle Scholar
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The ClustalX windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25: 4876-4882. 10.1093/nar/25.24.4876.PubMedPubMed CentralView ArticleGoogle Scholar
- Yang Z, Nielsen R: Estimating synonymous and nonsynonymous substitution rates under realistic evolutionary models. Mol Biol Evol. 2000, 17: 32-43.PubMedView ArticleGoogle Scholar
- Zhang Z, Li J, Zhao XQ, Wang J, Wong GK, Yu J: KaKs_Calculator: calculating Ka and Ks through model selection and model averaging. Genomics Proteomics Bioinformatics. 2006, 4: 259-263. 10.1016/S1672-0229(07)60007-2.PubMedView 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 cited.