De novo sequencing and characterization of floral transcriptome in two species of buckwheat (Fagopyrum)
© Logacheva et al; licensee BioMed Central Ltd. 2011
Received: 18 October 2010
Accepted: 13 January 2011
Published: 13 January 2011
Transcriptome sequencing data has become an integral component of modern genetics, genomics and evolutionary biology. However, despite advances in the technologies of DNA sequencing, such data are lacking for many groups of living organisms, in particular, many plant taxa. We present here the results of transcriptome sequencing for two closely related plant species. These species, Fagopyrum esculentum and F. tataricum, belong to the order Caryophyllales - a large group of flowering plants with uncertain evolutionary relationships. F. esculentum (common buckwheat) is also an important food crop. Despite these practical and evolutionary considerations Fagopyrum species have not been the subject of large-scale sequencing projects.
Normalized cDNA corresponding to genes expressed in flowers and inflorescences of F. esculentum and F. tataricum was sequenced using the 454 pyrosequencing technology. This resulted in 267 (for F. esculentum) and 229 (F. tataricum) thousands of reads with average length of 341-349 nucleotides. De novo assembly of the reads produced about 25 thousands of contigs for each species, with 7.5-8.2× coverage. Comparative analysis of two transcriptomes demonstrated their overall similarity but also revealed genes that are presumably differentially expressed. Among them are retrotransposon genes and genes involved in sugar biosynthesis and metabolism. Thirteen single-copy genes were used for phylogenetic analysis; the resulting trees are largely consistent with those inferred from multigenic plastid datasets. The sister relationships of the Caryophyllales and asterids now gained high support from nuclear gene sequences.
454 transcriptome sequencing and de novo assembly was performed for two congeneric flowering plant species, F. esculentum and F. tataricum. As a result, a large set of cDNA sequences that represent orthologs of known plant genes as well as potential new genes was generated.
Transcriptome sequencing is a convenient way to rapidly obtain information on the expressed fraction of genome. With the advent of next-generation sequencing transcriptomic data for many species became available. However many important taxonomic groups of living organisms remain underrepresented. This is especially pressing problem for plants since few of them have been a priority for large-scale sequencing projects. The present study aimed at filling the gap in the genomic data for the genus Fagopyrum, a group of plants important for both practical and theoretical reasons.
Fagopyrum (buckwheat) belongs to the eudicot family Polygonaceae. The genus Fagopyrum comprises about 17 species; one of the species, F. esculentum (common buckwheat), is a crop and honey-producing plant widely cultivated in several countries, in particular, Canada, China, Japan, Russia and Ukraine. Recently, an extensive search for the wild ancestor of the cultivated buckwheat identified the ancestral group, called F. esculentum subspecies ancestrale[1, 2]. Two more distant relatives of the cultivated buckwheat are F. homotropicum and F. tataricum. Both of them have an important favorable trait absent in the common buckwheat, the ability to self-pollination (self-compatibility). Thus they have a potential for use in breeding for the development of self-compatible cultivars.
Despite its economic importance, molecular studies on buckwheat are few. They are mainly confined to the development of molecular markers based on anonymous DNA sequences (RAPD, SSR, AFLP) and the characterization of several proteins [4–6]. Also, the molecular systematics of the genus Fagopyrum was extensively studied using nuclear and plastid sequences [7–9], and a complete sequence of the chloroplast genome of the common buckwheat was reported recently . These data, however, have only limited applicability for the buckwheat breeding and improvement. The identification of genes expressed during flower development can enable the search for candidate genes responsible for agriculturally important traits.
The utility of large-scale buckwheat gene sequencing data is not limited to their potential practical applications. They will also contribute to the plant comparative and evolutionary genomics. Not only the buckwheat, but the entire family Polygonaceae was out of scope of molecular genetic studies. The Polygonaceae is a middle-sized family (ca. 1000 species) characterized by several morphological features that are not typical for other eudicots. Among them is an unusual floral structure with flowers having a uniseriate perianth not differentiated into sepals and petals. This feature, as well as the trimerous perianth found in several Polygonaceae genera and considered as ancestral state for the family, is common for monocots and basal angiosperms but rarely found in eudicots.
The genetic control of the floral development has been well studied in several model plants, primarily Arabidopsis thaliana and some gene interactions have been found to be conserved in a variety of species . However, it is widely accepted that for understanding of the plant morphological evolution new model systems representing different lineages of the angiosperms should be selected [12, 13]. Given the morphological peculiarities of Polygonaceae, Fagopyrum, and specifically, F. esculentum, a fast growing annual weed with high seed production, is a plausible candidate for such a system.
We expect that characterization of buckwheat transcriptomes will contribute not only to the genetics of this genus. The order Caryophyllales is characterized by a striking diversity of the floral structure and the perianth evolution that can hardly be explained based only on the knowledge about classic model species . Few transcriptome datasets of the Caryophylalles are available to date [17, 18]. Publicly available, annotated set of buckwheat genes expressed in the flower and inflorescence resulting from this study will facilitate comparative and evolutionary studies within the Caryophylalles, important for the plant evo-devo field.
454 sequencing and assembly
Pyrosequencing of normalized cDNA libraries resulted in 266782 reads for F. esculentum and 229031 for F. tataricum with the average length of 349 and 341 nucleotides respectively. The raw data were deposited in the NCBI Sequence Read Archive (SRA) under the accession number SRA023408.
Characteristics of raw data and contigs
Number of reads
Length, mean (min-max)
Number of contigs
Length, mean (min-max)
Number of reads per contig, mean (min-max)
Since several full-length coding sequences of buckwheat genes were available in GenBank, we used this information to assess the assembly quality by mapping the contigs onto these sequences. Of 34 available sequences of F. esculentum, 30 were found in our assembly. In the latter, the average fraction of coding nucleotides covered by at least one read was 73.4%; while for 11 genes the coding sequence was covered completely (Additional file 3). This demonstrates that the assembled contigs can be a useful source of complete or near-complete coding sequences for F. esculentum.
Summary of BLAST search and annotation
Number of contigs with BLAST match
Number of contigs with GO annotation
Number of contigs without BLAST matches
Number of contigs with predicted ORF/mean length of ORF, bp
25345/368 (30-1637) 25115 > 90
24984/390 (39-1983) 24743 > 90
Species provided most of all best BLAST hits
Total number of singletons
Number of singletons with BLAST match
Number of singletons with GO annotation
Number of singletons without BLAST matches
Species provided most of all best BLAST hits
For the singletons, the fraction of sequences that had BLAST matches was lower than for the contigs (Table 2), as expected given their smaller length . The taxonomic distribution of species that provided most top hits was the same as for the contigs.
An automatic ORF prediction revealed open reading frames with length more than 90 bp in about 98% contigs in both buckwheat species. This length, though being much less than the average length of protein-coding sequence in plants, is higher than expected to arise by chance (assessed by random shuffling of the contigs with subsequent ORF prediction) suggesting that most of the contigs correspond to protein-coding genes. This is consistent with the method used for cDNA synthesis, as it is based on the presence of polyA-tail on the 3'-end of most mRNAs, the library is highly enriched in mRNA sequences. Interestingly, most contigs that had no BLAST matches, also contained predicted ORFs with length exceeding the one expected for a random sequence of a similar GC-content (data not shown).
The contigs that had no BLASTX hits but contained an ORF with a length exceeding 90 nucleotides were searched against the Pfam database (e-value threshold 1). A similar number of translated contigs had significant match in the database for both species, 1795 for F. esculentum and 1775 for F. tataricum, with 1108 and 1156, respectively, having a match in pfamA families. Among the latter, about 50% (47% for F. esculentum and 50% for F. tataricum) had hits to families containing at least one protein from plants.
Putative orthologs of buckwheat genes involved in flower development (defined according to Alvarez-Buylla et al. 2008) identified by bidirectional BLAST search
WIGGUM/ENHANCED RESPONSE TO ABSCISIC ACID1
AUXIN RESPONSE FACTOR 5/MONOPTEROS
AUXIN RESPONSE FACTOR 6
AUXIN RESPONSE FACTOR 8
BARELY ANY MERISTEM 1
CORONATINE INSENSITIVE 1
DELAYED DEHISCENCE1/OPDA-REDUCTASE 3
FATTY ACID DESATURASE 7
GA INSENSITIVE DWARF1A
GA INSENSITIVE DWARF1B
GA INSENSITIVE DWARF1C
RECEPTOR-LIKE PROTEIN KINASE 2
NO TRANSMITTING TRACT
Since this type of search is not able to detect duplications, we performed a manual inspection of the results of BLAST search for several gene families. For the B-class MADS-box gene, APETALA3 (AP3), we were able to find two highly similar sequences. The similarity between them is lower than that expected for alleles, with multiple indels and non-synonymous substitutions. Thus we assume that they represent paralogs arisen from a recent duplication. In another species from the same family, Rumex acetosa, two genes with high similarity to AP3 were also found . Hence, this duplication likely occured in the common ancestor of the Polygonaceae.
Comparative analysis of the buckwheat transcriptomes
We addressed the question whether these transcriptomes completely consist of orthologous genes or there are genes expressed in one species but with orthologs that are not represented (absent or not expressed) in the other species. Such transcripts are plausible candidates for the determinants of morphological differences between the species.
Characteristics of potential differentially expressed genes (PDEG) fraction of buckwheat transcriptomes
total number of PDEG
without BLAST matches
with GO annotation
number of contigs covered by more than 10 reads
In addition to the direct comparison of the PDEG sets, we addressed the question whether they represent random subsets of the overall transcriptomic datasets. To test this, the enrichment analysis was performed, with the PDEG set taken as a test group and the overall set taken as a reference group (p = 0.05). No enrichment by any specific category was found for F. tataricum; while several categories were found to be overrepresented in the F. esculentum PDEG set (Table 4), in particular categories obviously related to retrotransposons (DNA integration, RNA-directed DNA polymerase activity) and to processes inherent to plants such as sugar metabolism, signaling and amino acid transport.
Expected number of differentially expressed genes in the PDEG sets
number of reads covering a contig
observed number of PDEG
expected number of differentially expressed genes
Shared single-copy genes in the 454 transcriptome assemblies and their phylogenetic analysis
An interesting result is the position of the Caryophyllales. With the addition of Fagopyrum and Silene, the sister relationships of the Caryophyllales gained high support (PP = 0.96 for nucleotide data and 1 for amino acid) in the Bayesian analysis. This is consistent with our previous data on the plastid genome phylogeny , but in the nuclear gene analysis this topology is present as highly supported for the first time. In the dataset of Duarte et al. , where the only caryophyllid species, Beta vulgaris was present, it was also grouped with asterids in the maximum likelihood tree but with very a low support (< 50%).
The supertree resulting from SDM is in most points congruent with the Bayesian trees and also supports sister relationships of the Caryophyllales and asterids (Additional file 5).
Utility of 454 sequencing for gene discovery in non-model species
Despite the economical importance of buckwheat and the great advances in DNA sequencing technologies no genomic data were available to date for any of the species of Fagopyrum, with the exception of the complete chloroplast genome sequence of F. esculentum ssp. ancestrale. Prior to this study, the total number of sequences in GenBank was 149 for F. esculentum and 121 (including 83 unannotated EST) for F. tataricum. Sequencing and assembly of genes expressed in flower and inflorescence presented in this study results in 25 thousand contigs for each species. Given that the average number of genes encoded in a plant nuclear genome is about 30 thousands (as estimated from seven completely sequenced genomes), we assume that our dataset represents a substantial fraction of the buckwheat genes. Among them are putative orthologs of genes that play key roles in the flower and inflorescence development, including those that are expressed at low level and are not widely represented in other transcriptome assemblies. These results are highly consistent with previous experimental and simulation-based studies [31–33] and provide further support for the use of combination of cDNA normalization and 454 sequencing for fast transcriptome characterization in non-model species. Application of the floral transcriptome sequence data for gene expression studies in F. esculentum mutants with altered flower development [34–36] will constitute the basis for the study of the genetic control of flower development in this species.
A fraction of contigs that did not produce any BLAST hits may also be an integral part of genomic data. The existence of such "non-blastable" sequences is reported in virtually all plant transcriptomes characterized to date, with their fraction varying from about 30 to 70% [23–25], depending on the species, the depth of sequencing and the parameters of BLAST search. There are two major groups of effects that may be responsible for this. The first group comprises the technical issues such as low quality of raw data, contamination by genomic DNA, inaccurate contig assembly and wrong choice of search parameters (e.g. too stringent e-value cut-off). The second group of effects, biological ones, is related to the characteristics of the data that are inherent to the species: the existence of non-coding cDNA fraction (primarily 5' and 3'UTR), lineage-specific genes (i.e. genes that are present in the genome being considered but absent from genomes represented in the databases) and fast-evolving genes, those that have orthologs in the other species but with a high level of sequence divergence precluding efficient recognition of orthologs.
In our study we tried to minimize the influence of technical effects. DNAse treatment of RNA was performed before cDNA synthesis to preclude genomic DNA contamination; also low quality data were excluded from the assembly process. As for the assembly accuracy, in the absence of a reference genome the possibility of assembly errors can not be completely ruled out. However given the relatively large length of 454-generated reads even chimeric combination of reads corresponding to different genes in one contig is unlikely to hamper efficient BLAST search. Indeed, assembly errors cause problems with similarity search only if short reads corresponding to small and weakly conserved fragments of different genes are combined into one contig.
Thus we assume that most of non-blastable contigs in our data are due to biological reasons. Some of them might indeed represent non-translated regions, but since the average length of predicted ORF in such contigs is significantly higher than in artificial sequences produced by random shuffling, we suggest that most of them correspond to protein-coding sequences, either lineage specific or highly variable. Currently there are no means to select one of these two options. Addition of more genomic and/or transcriptomic data for the Polygonaceae and other Caryophyllales is expected to improve the annotation and gene prediction in buckwheat.
The favourable effect of using the data from closely related species for annotation was reported for the analysis of chestnut (Castanea, Fagaceae) transcriptomes . Several genes that could not be annotated using the information from the Arabidopsis proteome (Castanea and Arabidopsis belong to different clades of rosids - eurosids I and eurosids II) were successfully annotated when compared with Populus (belonging to eurosids I) proteome. This emphasizes that despite the availability of large-scale sequence data for model plant species and the increased rate of data acquisition for non-model species denser taxon sampling will further improve plant gene annotation. Given that the Caryophylalles is a group of exceptional interest from the point of view of evolutionary developmental biology we expect that more transcriptomic data will be available for this group.
Comparison of two buckwheat transcriptomes: finding potential differentially expressed genes
Gene expression is often compared within one species at different developmental stages or in different conditions [22, 37], and several efficient tools have been developed to analyze the enrichment in specific sets of genes (reviewed in ). The interspecific comparisons are rare and confined to the cases when it is possible to use sequence information from a phylogenetically close model species [39, 40]. The possibility to use de novo sequencing data for the comparison of gene expression (even on the qualitative level, that is, to identify the presence or absence of expression) in different species has, to our knowledge, never been addressed in plant genomics.
Here we used an approach that is similar in its concept to the digital transcriptome subtraction . The latter is based on the similarity search of sequences from one species against the other and selection of sequences that do not have significant matches. As expected, for most sequences from both datasets the sequences with high similarity were found in another species, either in contigs, or in singletons. However, about four thousands contigs in each species were found to be "species-specific". It is difficult to directly associate these contigs with the observed differences between F. esculentum and F. tataricum, the more so because about two thirds of them do not have significant BLAST hits in the nr database, and thus lack any functional annotation.
Still, some patterns are likely caused by species-specific expression. Firstly, the F. esculentum PDEG set shows a strong prevalence of contigs highly similar to retrotransposons. This suggests that either the F. esculentum genome has more retrotransposons than the genome of F. tataricum, or more of them are expressed. The former explanation seems to be more plausible, considering the fact that the genome size of F. esculentum is three-fold larger than that of F. tataricum. While such differences in the genome size in plants are often caused by the polyploidy, it is not the case for the buckwheat since both species have the same chromosome number [42, 43].
Secondly, the F. esculentum PDEG set is enriched in genes related to the disaccharide metabolism. Mono- and disaccharides, primarily glucose, fructose and sucrose, are a major component of the F. esculentum nectar . F. tataricum, being a self-pollinated species, does not produce nectar. So the disaccharide metabolism genes indeed might be differentially expressed due to the difference in the nectar production. However, any conclusion about the differential expression should be treated with caution and verified by more precise methods such as quantitative RT-PCR. The set of "unique" contigs identified in our study is a mixture of genes whose orthologs are actually not expressed in the other species and those genes that are expressed in both species but are absent in the sequence data from one of the species due to some random fluctuations on different stages of RNA and cDNA processing. Indeed, our statistical analysis indicates that about 25-30% of PDEG are really differentially expressed. One of the possible issues responsible for this is the normalization because in normalized libraries more transcripts are represented and thus the probability of being sequenced is lower for each transcript than in case of non-normalized library.
The probability that a gene is expressed, given the observations of no reads from the other genome, is proportional to its coverage. It is very low for low-covered genes but among genes represented by a contig covered by more than 7 at least 95% of genes are expected to be differentially expressed (Table 5). Thus, further investigation of the genetic basis of specific differences between F. esculentum and F. tataricum may be conducted in two ways. The first one is the survey of the expression of PDEGs from the categories that are enriched in either species and that are likely to have species-specific expression for biological reasons (retrotransposons, sugar metabolism genes). The second one is the analysis of those genes that have a high number of reads in one transcriptome but are absent in thе other. A combination of these approaches may provide insights in the genetic and evolutionary mechanisms underlying morphological and physiological differentiation in the genus Fagopyrum. Further insights, including quantitative analysis of gene expression levels are also expected to result from the sequencing of non-normalized cDNA. This approach demonstrated its utility on several plant species [37, 45].
Phylogenetic utility of transcriptome sequence data
Besides its potential utility for functional genetic studies, the data from large-scale transcriptome sequencing is an important source of information for phylogenetic analysis. Phylogenetic studies based on transcriptome sequences yielded well-resolved and highly supported tree topologies for many groups of animals . In plants, however, the phylogenomic approach has been mostly limited to plastid genome sequences . There are two reasons for that. The first one is the relative scarcity of plant genome-scale data. Secondly, the starting point of any phylogenetic analysis, identification of orthologous genes, is a major problem for plants due to abundant genome-wide and small scale duplications.
Recently a set of single-copy genes was proposed for phylogenetic analysis on different taxonomic level. Though a number of articles reporting de novo plant transcriptome sequence data were published [25, 33] none of them reported the use of these data for phylogenetic analysis. Here we attempted to validate the use of these genes for inferring angiosperm phylogeny by including in the analysis the data from four 454 transcriptome assemblies, with focus on Caryophyllales and asterids. Orthologs of most of these genes were found in all these transcriptomes. The tree inferred from the Bayesian analysis of these genes is similar to the ML tree from Duarte et al.  but has a higher resolution and node support. We suggest the increased resolution and node support reflects favourable effect of improved taxon sampling.
A major problem observed in our phylogenetic trees is the incongruence between different methods. The Bayesian inference and supertree (SDM) approach yield congruent topologies but these topologies are not congruent with those inferred from the maximum parsimony analysis. Also, MP analysis results in a tree with very low bootstrap support values. This may indicate that the MP is not applicable for transcriptomic data, where high variation is combined with lot of missing data. It is known that several types of data may require additional adjustments of procedures used for phylogenetic analysis (primarily alignment). For example, for non-coding sequences of plastid genomes special rules of alignment were developed and this increased resolution and reliability of phylogenetic trees inferred from these sequences [47, 48]. The evolutionary peculiarities of plastid sequences are well studied [47, 49]. On the contrary, little is known about the nuclear genes used for the phylogenetic analysis in this study. It is probable that the angiosperm-wide survey of the patterns of evolution of these genes will identify the source of incongruence between different methods of the analysis and, as the result, improve their utility for plant phylogenetics.
However, even under limited taxon sampling these genes provide valuable information on phylogenetic relationships within the angiosperms. The novel result inferred from these data is the strong support of sister relationships of Caryophyllales (including buckwheat) and asterids. Earlier we have shown that these relationships are strongly supported by the chloroplast genome sequence data . It is well known however that the results of phylogenetic analysis of chloroplast gene datasets may be misleading and thus require verification by independent data [50, 51]; the present study reports such verification for the phylogenetic position of Caryophyllales.
The results presented here emphasize that even shallow-coverage transcriptome sequence data are an important source of information for phylogenetic applications. The addition of such data from various lineages of flowering plants is expected to improve greatly the resolution, support and reliability of phylogenetic trees and to provide novel insights into the evolution of angiosperms.
454 pyrosequencing of normalized cDNA libraries produced a large set of cDNA sequences for two congeneric plant species, Fagopyrum esculentum and F. tataricum. These sequences are an important resource for the evolutionary and developmental genetics in these species and in the order Caryophyllales. Analysis of single copy genes derived from transcriptome sequence data have great potential for inferring angiosperm phylogeny, especially with increased taxonomic sampling.
Total RNA was extracted from developing inflorescences of Fagopyrum esculentum ssp. ancestrale and F. tataricum using the RNEasy Plant Mini kit (Qiagen) with following modifications: 50 μl of Plant RNA Isolation Aid reagent (Ambion) and 10 μl of b-mercaptoethanol were added to the lysis buffer RLT prior to homogenization and the homogenization was performed without liquid nitrogen. All other steps were performed according to the manufacturer's instructions. To avoid genomic DNA contamination RNA was treated with RNase-free DNase I (Qiagen). The RNA integrity was assessed by the agarose gel electrophoresis (1% gel with addition of SYBRGreen dye).
cDNA synthesis, amplification and normalization and 454 sequencing
Total RNA was used for double-stranded cDNA synthesis using the SMART approach . To increase the representation of weakly expressed genes, SMART-prepared amplified cDNA was then normalized using the DSN normalization method . Normalization included cDNA denaturation/reassociation, treatment by a duplex-specific nuclease (DSN, ) and amplification of the normalized fraction. cDNA synthesis and normalization were performed by the Evrogen company, the detailed protocols are available in Additional file 6.
The efficiency of normalization was tested by real-time PCR. Fragments of GAPDH (a housekeeping gene expressed at a high level) and of the WUS ortholog (encodes a transcription factor expressed at a low level) were amplified from the normalized and non-normalized libraries. The 32-fold decrease of GAPDH expression was observed in the normalized library. Then 7 μg of normalized cDNA were fractionated and sequenced using 454 GS-FLX sequencer with 454 Titanium chemistry. Sequencing was performed by The Centre for Applied Genomics, The Hospital for Sick Children, Toronto, Canada. Each sample was run on a half of a picotiter plate.
Assembly, annotation and comparative analysis of transcriptomes
SeqClean  was used to remove polyA sequences from the raw data. We did not use SeqClean to remove adapter sequences since it works only with 100-92% identity segments and in our case adapters contained many sequencing errors caused by the homopolymer runs. To remove the adapters with less than 92% identity we developed two scripts, which are available upon request. MIRA assembler  version 3.0.0rc4 was used for the 454 data assembly. The assembler was run in its 'accurate' mode with the assembly type set as 'EST'.
The resulting contigs were subject to the BLASTX search against the non-redundant protein database (nr) with the e-value threshold of 10-6 and the HSP length cut-off of 50, as implemented in the BLAST2Go program . The ORF Predictor tool was used to identify ORFs (http://proteomics.ysu.edu/tools/OrfPredictor.html, ). Searches against the Arabidopsis proteome (e-value threshold of 10-6) and the searches between F. esculentum and F. tataricum (e-value of 10-15) were performed using BLASTX/TBLASTN and BLASTN, respectively, as implemented in BioEdit 7.09 . The Gene Ontology annotation  was done using BLAST2Go with the annotation cut-off of 10-5. We used simplified GOSlim (plant) annotations (developed by TAIR, ) to select genes expressed in the flower and in the inflorescence. For the identification of genes that are represented in the transcriptome of one species of buckwheat but are not represented in the other (here and further called potential differentially expressed genes, PDEG), a Perl script was developed that processes the BLAST output and extracts sequences that do not have significant matches. The enrichment analysis was performed using Fisher's exact test as implemented in Gossip/BLAST2Go . The statistical assessment of the number of false positives in the PDEG sets was performed as follows. First, the probability p(i,j,s) to observe i reads for the gene in one species (s) given j reads for its ortholog in the other species was calculated. For each possible j and s, the observed number of orthologous pairs with j reads for gene in species s fitted well the Poisson model after filtering out a small number of outliers. The model parameters λ(j,s) were estimated from these data using the maximum-likelihood method, taking into account that pairs with zero number of reads for any of the genes were not observed for obvious reasons. Finally, the number of false positives was calculated as Σ k(j,s) p(0,j,s) where k(j,s) is the total number of genes with j reads in species s.
For the phylogenetic analysis, a set of 13 single copy genes from 69 taxa, from the  study was used. The orthologs of these genes were found by the BLAST search of corresponding Arabidopsis genes in three 454-sequenced transcriptomic datasets of Caryophylalles, F. esculentum, F. tataricum, and Silene latifolia (Filatov D.A., personal communication) and two transcriptomes of asterids: Artemisia annua and Craterostigma plantagineum. Since the orthologs from the two buckwheat species are highly similar (> 95%) the sequences from only one species (F. tataricum) used. Sequences of F. tataricum, S. latifolia, A. annua and C. plantagineum were added to the alignment from  then the whole set was checked for frameshifts and the latter, if present, were corrected by inserting N to recover in frame translation. Sequences of transcripts were translated into amino acid sequences by BioEdit and further aligned using MUSCLE ver. 3.6 . The nucleotide sequence alignment was overlaid on the amino acid sequence alignment. Gap-rich positions and positions of ambiguous alignment were excluded from analyses (789 in total). The phylogenetic analysis using the maximum parsimony (MP) method was performed using PAUP* ver. 4.0b8 . The MP analysis involved a heuristic search using TBR branch swapping and 500 random addition replicates. The non-parametric bootstrap analysis  was performed with 1000 replicates with TBR branch swapping. The Bayesian inference of phylogeny was explored using the MrBayes program ver. 3.1.2 . The Bayesian analysis was applied both to the amino acid and nucleotide data with two runs with three chains in each. Each chain started with a random tree, 1,000,000 replicates for amino acid data and 5,000,000 replicates for nucleotide data were generated, trees were sampled every 100 generations. The number of discarded trees was determined using the cold-chain log-likelihood examination. The JTT model of amino acid substitutions was determined by the Bayesian information criterion in Modelgenerator ver. 0.43 , the GTR+I+Γ model of nucleotide substitution was determined by the Akaike information criterion in Modeltest ver. 3.7 , the data were treated as a single partition. As an alternative to the combined analysis a super distance matrix (SDM) method  was applied. The SDM-algorithm combines the evolutionary distances obtained from each gene into a single distance supermatrix, therefore for each gene a corresponding model of evolution was determined in Modeltest, then the maximum likelihood distances according to the selected model were computed with PAUP*. A supertree was constructed using Fitch program from PHYLIP package .
List of abbreviations
amplified fragment length polymorphism
polymerase chain reaction
potential differentially expressed genes
random amplified polymorphic DNA
reverse transcription polymerase chain reaction
super distance matrix
simple sequence repeat
The authors are grateful to Dmitry Shagin (Evrogen Lab, Moscow, Russia) for cDNA normalization and Sergio Pereira (The Centre for Applied Genomics, The Hospital for Sick Children, Toronto, Canada) for 454 sequencing, to Andrei Mironov (Moscow State University, Moscow, Russia) for useful discussion, to Dmitry Filatov (University of Oxford, Oxford, UK) for providing unpublished data on Silene latifolia genes used for phylogenetic analysis and to two anonymous reviewers for their helpful comments. The study is supported by the Russian Foundation for Basic Research (project №09-04-01363), the Ministry of Education and Science of Russian Federation (contracts P913 and P315 and project "Phylogenetic analysis of complex selection in molecular evolution"), the President's program for support of young scientists and the Russian Academy of Science Presidium Program "Molecular and Cell Biology".
- Ohnishi O: Discovery of the wild ancestor of common buckwheat. Fagopyrum. 1991, 11: 5-10.
- Ohnishi O: Search for the wild ancestor of buckwheat. I. Description of new Fagopyrum (Polygonaceae) species and their distribution in China and Himalayan hills. Fagopyrum. 1998, 15: 18-28.
- Ohnishi O, Matsuoka Y: Search for the wild ancestor of buckwheat. II. Taxonomy of Fagopyrum (Polygonaceae) species based on morphology, isozymes and cpDNA variability. Genes and Genetic Systems. 1996, 71: 383-390. 10.1266/ggs.71.383.View Article
- Yasui Y, Wang Y, Ohnishi O, Campbell CG: Amplified fragment length polymorphism linkage analysis of common buckwheat (Fagopyrum esculentum) and its wild self-pollinated relative Fagopyrum homotropicum. Genome. 2004, 47: 345-351. 10.1139/g03-126.PubMedView Article
- Sharma R, Jana S: Species relationships in Fagopyrum revealed by PCR-based DNA fingerprinting. Theoretical and Applied Genetics. 2002, 105: 306-312. 10.1007/s00122-002-0938-9.PubMedView Article
- Licen M, Kreft I: Buckwheat (Fagopyrum esculentum Moench) low molecular weight seed proteins are restricted to the embryo and are not detectable in the endosperm. Plant physiology and biochemistry. 2005, 43: 862-865. 10.1016/j.plaphy.2005.08.002.PubMedView Article
- Yasui Y, Ohnishi O: Phylogenetic relationships among Fagopyrum species revealed by the nucleotide sequences of the ITS region of the nuclear rRNA gene. Genes and Genetic Systems. 1998, 73: 201-210. 10.1266/ggs.73.201.PubMedView Article
- Ohsako T, Ohnishi O: Intra- and interspecific phylogeny of wild Fagopyrum (Polygonaceae) species based on nucleotide sequences of noncoding regions in chloroplast DNA. American Journal of Botany. 2000, 87: 573-582. 10.2307/2656601.PubMedView Article
- Nishimoto Y, Ohnishi O, Hasegawa M: Topological incongruence between nuclear and chloroplast DNA trees suggesting hybridization in the urophyllum group of the genus Fagopyrum (Polygonaceae). Genes and Genetic Systems. 2003, 78: 139-153. 10.1266/ggs.78.139.PubMedView Article
- Logacheva MD, Samigullin TH, Dhingra A, Penin AA: Comparative chloroplast genomics and phylogenetics of Fagopyrum esculentum ssp. ancestrale -a wild ancestor of cultivated buckwheat. BMC Plant Biology. 2008, 8: 59-10.1186/1471-2229-8-59.PubMed CentralPubMedView Article
- Ma H, dePamphilis CW: The ABCs of floral evolution. Cell. 2000, 101: 5-8. 10.1016/S0092-8674(00)80618-2.PubMedView Article
- Baum DA, Doebley J, Irish VF, Kramer EM: Response: Missing links: the genetic architecture of flower and floral diversification. Trends in Plant Science. 2002, 7: 31-34. 10.1016/S1360-1385(01)02181-1.View Article
- Nutt P, Ziermann J, Hintz M, Neuffer B, Theißen G: Capsella as a model system to study the evolutionary relevance of floral homeotic mutants. Plant Systematics and Evolution. 2006, 259: 217-235. 10.1007/s00606-006-0420-2.View Article
- Albert VA, Soltis DE, Carlson JE, Farmerie WG, Wall PK, Ilut DC, Solow TM, Mueller LA, Landherr LL, Hu Y, Buzgo M, Kim S, Yoo MJ, Frohlich MW, Perl-Treves R, Schlarbaum SE, Bliss BJ, Zhang X, Tanksley SD, Oppenheimer DG, Soltis PS, Ma H, DePamphilis CW, Leebens-Mack JH: Floral gene resources from basal angiosperms for comparative genomics research. BMC Plant Biology. 2005, 5: 5-10.1186/1471-2229-5-5.PubMed CentralPubMedView Article
- Samimy C, Bjorkman T, Siritunga D, Blanchard L: Overcoming the barrier to interspecific hybridization of Fagopyrum esculentum with Fagopyrum tataricum. Euphytica. 1996, 91: 323-330. 10.1007/BF00033094.View Article
- Brockington SF, Roolse A, Ramdial J, Moore MJ, Crawley S, Dhingra A, Hilu K, Soltis DE, Soltis PS: Phylogeny of the Caryophyllales sensu lato: revisiting hypotheses on pollination biology and perianth differentiation in the core Caryophyllales. International Journal of Plant Science. 2009, 170: 627-643. 10.1086/597785.View Article
- DFCI Ice_plant (Mesembryanthemum crystallinum) Gene Index. [http://compbio.dfci.harvard.edu/cgi-bin/tgi/gimain.pl?gudb=ice_plant]
- Silene vulgaris Transcriptome Project. [http://faculty.virginia.edu/silenetranscriptome/]
- Jaillon O, Aury JM, Noel B, Policriti A, Clepet C, Casagrande A, Choisne N, Aubourg S, Vitulo N, Jubin C, Vezzi A, Legeai F, Hugueney P, Dasilva C, Horner D, Mica E, Jublot D, Poulain J, Bruyère C, Billault A, Segurens B, Gouyvenoux M, Ugarte E, Cattonaro F, Anthouard V, Vico V, Del Fabbro C, Alaux M, Di Gaspero G, Dumas V, Felice N, Paillard S, Juman I, Moroldo M, Scalabrin S, Canaguier A, Le Clainche I, Malacrida G, Durand E, Pesole G, Laucou V, Chatelet P, Merdinoglu D, Delledonne M, Pezzotti M, Lecharny A, Scarpelli C, Artiguenave F, Pè ME, Valle G, Morgante M, Caboche M, Adam-Blondon AF, Weissenbach J, Quétier F, Wincker P: The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007, 449: 463-467. 10.1038/nature06148.PubMedView Article
- Tuskan GA, Difazio S, Jansson S, Bohlmann J, Grigoriev I, Hellsten U, Putnam N, Ralph S, Rombauts S, Salamov A, Schein J, Sterck L, Aerts A, Bhalerao RR, Bhalerao RP, Blaudez D, Boerjan W, Brun A, Brunner A, Busov V, Campbell M, Carlson J, Chalot M, Chapman J, Chen GL, Cooper D, Coutinho PM, Couturier J, Covert S, Cronk Q, Cunningham R, Davis J, Degroeve S, Dejardin A, Depamphilis C, Detter J, Dirks B, Dubchak I, Duplessis S, Ehlting J, Ellis B, Gendler K, Goodstein D, Gribskov M, Grimwood J, Groover A, Gunter L, Hamberger B, Heinze B, Helariutta Y, Henrissat B, Holligan D, Holt R, Huang W, Islam-Faridi N, Jones S, Jones-Rhoades M, Jorgensen R, Joshi C, Kangasjarvi J, Karlsson J, Kelleher C, Kirkpatrick R, Kirst M, Kohler A, Kalluri U, Larimer F, Leebens-Mack J, Leple JC, Locascio P, Lou Y, Lucas S, Martin F, Montanini B, Napoli C, Nelson DR, Nelson C, Nieminen K, Nilsson O, Pereda V, Peter G, Philippe R, Pilate G, Poliakov A, Razumovskaya J, Richardson P, Rinaldi C, Ritland K, Rouze P, Ryaboy D, Schmutz J, Schrader J, Segerman B, Shin H, Siddiqui A, Sterky F, Terry A, Tsai CJ, Uberbacher E, Unneberg P, Vahala J, Wall K, Wessler S, Yang G, Yin T, Douglas C, Marra M, Sandberg G, Peer Y, Rokhsar D: The genome of black cottonwood, Populus trichocarpa (Torr. & Gray). Science. 2006, 313: 1596-1604. 10.1126/science.1128691.PubMedView Article
- Ricinus communis genome project. [http://gsc.jcvi.org/projects/msc/ricinus_communis/]
- Rodriguez MC, Edsgard D, Sarfraz HS, Alquezar D, Rasmussen M, Gilbert T, Nielsen BH, Bartels D, Mundy J: Transcriptomes of the desiccation tolerant resurrection plant Craterostigma plantagineum. Plant Journal. 2010, 63: 212-228. 10.1111/j.1365-313X.2010.04243.x.PubMedView Article
- 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 CentralPubMedView Article
- Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nature Genetics. 2000, 25: 25-29. 10.1038/75556.PubMed CentralPubMedView Article
- Sun C, Li Y, Wu Q, Luo H, Sun Y, Song J, Lui EM, Chen S: De novo sequencing and analysis of the American ginseng root transcriptome using a GS FLX Titanium platform to discover putative genes involved in ginsenoside biosynthesis. BMC Genomics. 2010, 11: 262-10.1186/1471-2164-11-262.PubMed CentralPubMedView Article
- Ainsworth C, Thangavelu M, Crossley S, Buchanan-Wollaston V, Parker J: Male and female flowers from the dioecious plant Rumex acetosa show different patterns of MADS-box gene expression. Plant Cell. 1995, 7: 1583-1598. 10.1105/tpc.7.10.1583.PubMed CentralPubMedView Article
- Duarte JM, Wall PK, Edger PP, Landherr LL, Ma H, Pires JC, Leebens-Mack J, dePamphilis CW: Identification of shared single copy nuclear genes in Arabidopsis, Populus, Vitis and Oryza and their phylogenetic utility across various taxonomic levels. BMC Evolutionary Biology. 2010, 10: 61-10.1186/1471-2148-10-61.PubMed CentralPubMedView Article
- Moore MJ, Soltis PS, Bell CD, Burleigh JG, Soltis DE: Phylogenetic analysis of 83 plastid genes further resolves the early diversification of eudicots. Proceedings of the National Academy of Sciences USA. 2010, 107: 4623-4628. 10.1073/pnas.0907801107.View Article
- Graham IA, Besser K, Blumer S, Branigan CA, Czechowski T, Elias L, Guterman I, Harvey D, Isaac PG, Khan AM, Larson TR, Li Y, Pawson T, Penfield T, Rae AM, Rathbone DA, Reid S, Ross J, Smallwood MF, Segura V, Townsend T, Vyas D, Winzer T, Bowles D: The genetic map of Artemisia annua L. identifies loci affecting yield of the antimalarial drug artemisinin. Science. 2010, 327: 328-331. 10.1126/science.1182612.PubMedView Article
- Bailey CD, Carr TG, Harris SA, Hughes CE: Characterization of angiosperm nrDNA polymorphism, paralogy, and pseudogenes. Molecular Phylogenetics and Evolution. 2003, 29: 435-455. 10.1016/j.ympev.2003.08.021.PubMedView Article
- Wall PK, Leebens-Mack J, Chanderbali AS, Barakat A, Wolcott E, Liang H, Landherr L, Tomsho LP, Hu Y, Carlson JE, Ma H, Schuster SC, Soltis DE, Soltis PS, Altman N, dePamphilis CW: Comparison of next generation sequencing technologies for transcriptome characterization. BMC Genomics. 2009, 10: 347-10.1186/1471-2164-10-347.PubMed CentralPubMedView Article
- Barakat A, DiLoreto DS, Zhang Y, Smith C, Baier K, Powell WA, Wheeler N, Sederoff R, Carlson JE: Comparison of the transcriptomes of American chestnut (Castanea dentata) and Chinese chestnut (Castanea mollissima) in response to the chestnut blight infection. BMC Plant Biology. 2009, 9: 51-10.1186/1471-2229-9-51.PubMed CentralPubMedView Article
- Zeng S, Xiao G, Guo J, Fei Z, Xu Y, Roe BA, Wang Y: Development of a EST dataset and characterization of EST-SSRs in a traditional Chinese medicinal plant, Epimedium sagittatum (Sieb. Et Zucc.) Maxim. BMC Genomics. 2010, 11: 94-10.1186/1471-2164-11-94.PubMed CentralPubMedView Article
- Fesenko AN, Fesenko IN, Logacheva MD, Penin AA: The role of gene TEPAL-LIKE BRACT (TLB) in the separation between the bracts and perianth in Fagopyrum esculentum Moench. Genetika. 2005, 41: 1644-1649.PubMed
- Logacheva MD, Fesenko IN, Fesenko AN, Penin AA: Genetic and morphological analysis of floral homeotic mutants tepal-like bract and fagopyrum apetala of Fagopyrum esculentum. Botany. 2008, 86: 367-375. 10.1139/B08-010.View Article
- Penin AA, Fesenko AN, Fesenko IN, Logacheva MD: Some characteristics of genetic control of Fagopyrum esculentum flower development. Wulfenia. 2009, 16: 1-11.
- Alagna F, D'Agostino N, Torchia L, Servili M, Rao R, Pietrella M, Giuliano G, Chiusano ML, Baldoni L, Perrotta G: Comparative 454 pyrosequencing of transcripts from two olive genotypes during fruit development. BMC Genomics. 2009, 10: 399-10.1186/1471-2164-10-399.PubMed CentralPubMedView Article
- Huang W, Sherman BT, Lempicki RA: Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Research. 2009, 37: 1-13. 10.1093/nar/gkn923.PubMed CentralView Article
- Enard W, Khaitovich P, Klose J, Zöllner S, Heissig F, Giavalisco P, Nieselt-Struwe K, Muchmore E, Varki A, Ravid R, Doxiadis GM, Bontrop RE, Pääbo S: Intra- and interspecific variation in primate gene expression patterns. Science. 2002, 296: 340-343. 10.1126/science.1068996.PubMedView Article
- Sharbel TF, Voigt ML, Corral JM, Galla G, Kumlehn J, Klukas C, Schreiber F, Vogel H, Rotter B: Apomictic and sexual ovules of Boechera display heterochronic global gene expression patterns. Plant Cell. 2010, 22: 655-671. 10.1105/tpc.109.072223.PubMed CentralPubMedView Article
- Feng H, Taylor JL, Benos PV, Newton R, Waddell K, Lucas SB, Chang Y, Moore PS: Human transcriptome subtraction using short sequence tags to search for tumor viruses in conjunctival carcinoma. Journal of Virology. 2007, 81: 11332-11340. 10.1128/JVI.00875-07.PubMed CentralPubMedView Article
- Nagano M, Aii J, Campbell C, Kawasaki S, Adachi T: Genome size analysis of the genus Fagopyrum. Fagopyrum. 2000, 17: 35-39.
- Quisenberry KS: Chromosome numbers in buckwheat species. Botanical Gazette. 1927, 83: 85-88. 10.1086/333697.View Article
- Vattala HD, Wratten SD, Phillips CB, Wäckers FL: The influence of flower morphology and nectar quality on the longevity of a parasitoid biological control agent. Biological Control. 2006, 39: 179-185. 10.1016/j.biocontrol.2006.06.003.View Article
- Yang P, Li X, Shipp MJ, Shockey JM, Cahoon EB: Mining the bitter melon (Momordica charantia L.) seed transcriptome by 454 analysis of non-normalized and normalized cDNA populations for conjugated fatty acid metabolism-related genes. BMC Plant Biology. 2010, 10: 250-10.1186/1471-2229-10-250.PubMed CentralPubMedView Article
- 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). Molecular Phylogenetics and Evolution. 2009, 53: 826-834. 10.1016/j.ympev.2009.08.014.PubMedView Article
- Kelchner SA: The evolution of non-coding chloroplast DNA and its application in plant systematics. Annals of the Missouri Botanical Garden. 2000, 87: 482-498. 10.2307/2666142.View Article
- Löhne C, Borsch T: Phylogenetic utility and molecular evolution of the petD group II intron in basal angiosperms. Molecular Biology and Evolution. 2005, 22: 317-332.PubMedView Article
- Learn GH, Shore JS, Furnier GR, Zurawski G, Clegg MT: Constraints on the evolution of plastid introns: the group II intron in the gene encoding tRNA-Val(UAC). Molecular Biology and Evolution. 1992, 9: 856-871.PubMed
- Degtjareva GV, Samigullin TH, Sokoloff DD, Valiejo-Roman CM: Gene sampling versus taxon sampling: is Amborella (Amborellaceae) a sister group to all other extant angiosperms?. Botanicheskiy Zhurnal. 2004, 89: 896-907.
- Soltis DE, Albert VA, Savolainen V, Hilu K, Qiu YL, Chase MW, Farris JS, Stefanovic S, Rice DW, Palmer JD, Soltis PS: Genome-scale data, angiosperm relationships, and 'ending incongruence': a cautionary tale in phylogenetics. Trends in Plant Science. 2004, 9: 477-483. 10.1016/j.tplants.2004.08.008.PubMedView Article
- Zhu YY, Machleder EM, Chenchik A, Li R, Siebert PD: Reverse transcriptase template switching, a SMART approach for full-length cDNA library construction. Biotechniques. 2001, 30: 892-897.PubMed
- Zhulidov PA, Bogdanova EA, Shcheglov AS, Vagner LL, Khaspekov GL, Kozhemyako VB, Matz MV, Meleshkevitch E, Moroz LL, Lukyanov SA, Shagin DA: Simple cDNA normalization using kamchatka crab duplex-specific nuclease. Nucleic Acids Research. 2004, 32: e37-10.1093/nar/gnh031.PubMed CentralPubMedView Article
- Shagin DA, Rebrikov DV, Kozhemyako VB, Altshuler IM, Shcheglov AS, Zhulidov PA, Bogdanova EA, Staroverov DB, Rasskazov VA, Lukyanov S: A novel method for SNP detection using a new duplex-specific nuclease from crab hepatopancreas. Genome Research. 2002, 12: 1935-1942. 10.1101/gr.547002.PubMed CentralPubMedView Article
- DFCI Gene Indices Software Tools. [http://compbio.dfci.harvard.edu/tgi/software/]
- Chevreux B, Wetter T, Suhai S: Genome sequence assembly using trace signals and additional sequence information. Computer Science and Biology: Proceedings of the German Conference on Bioinformatics. 1999, 99: 45-56.
- Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M: Blast2GO: A universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21: 3674-3676. 10.1093/bioinformatics/bti610.PubMedView Article
- Min XJ, Butler G, Storms R, Tsang A: OrfPredictor: predicting protein-coding regions in EST-derived sequences. Nucleic Acids Research. 2005, W677-W680. 10.1093/nar/gki394.
- Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symposium Series. 1999, 41: 95-98.
- Rhee SY, Beavis W, Berardini TZ, Chen G, Dixon D, Doyle A, Garcia-Hernandez M, Huala E, Lander G, Montoya M, Miller N, Mueller LA, Mundodi S, Reiser L, Tacklind J, Weems DC, Wu Y, Xu I, Yoo D, Yoon J, Zhang P: The Arabidopsis Information Resource (TAIR): a model organism database providing a centralized, curated gateway to Arabidopsis biology, research materials and community. Nucleic Acids Research. 2003, 31: 224-228. 10.1093/nar/gkg076.PubMedView Article
- Blüthgen N, Brand K, Cajavec B, Swat M, Herzel H, Beule D: Biological profiling of gene groups utilizing Gene Ontology. Genome Informatics. 2005, 16 (1): 106-15.PubMed
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Research. 2004, 32: 1792-1797. 10.1093/nar/gkh340.PubMed CentralPubMedView Article
- Swofford DL: PAUP*: Phylogenetic analysis using parsimony (*and other methods), ver. 4.0. 2003, Sunderland MA: Sinauer Associates
- Felsenstein J: Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985, 39: 783-791. 10.2307/2408678.View Article
- Ronquist F, Huelsenbeck JP: MRBAYES 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003, 19: 1572-1574. 10.1093/bioinformatics/btg180.PubMedView Article
- Keane TM, Creevey CJ, Pentony MM, Naughton TJ, McInerney JO: Assessment of methods for amino acid matrix selection and their use on empirical data shows that ad hoc assumptions for choice of matrix are not justified. BMC Evolutionary Biology. 2006, 6: 29-10.1186/1471-2148-6-29.PubMed CentralPubMedView Article
- Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.PubMedView Article
- Criscuolo A, Berry V, Douzery EJ, Gascuel O: SDM: a fast distance-based approach for (super) tree building in phylogenomics. Systematic Biology. 2006, 55: 740-755. 10.1080/10635150600969872.PubMedView Article
- Felsenstein J: PHYLIP (Phylogeny Inference Package) version 3.6. 2004, Department of Genome Sciences, University of Washington, Seattle
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.