BACKGROUND
the solanaceae, or nightshade family, is a dicot plant family that includes many economically important genera that are used in agriculture, horticulture, and other industries.
family members include the tuber bearing potato ; a large number of fruit-bearing vegetables, such as peppers , tomatoes , and eggplant ; leafy tobacco ; and ornamental flowers from the petunia and solanum genera.
tomato is generally considered to be a model crop plant species, for which many high-quality genetic and genomic resources are available, such as high-density molecular maps  <cit> , many well-characterized near-isogenic lines , and rich collections of ests and full-length cdnas  <cit> .
potato is the most important crop within the solanaceae, ranking fourth as a world food crop following wheat, maize and rice.
similar resources are available for potato, including an ultra-high density linkage map  <cit> , a collection of phenotype data  <cit> , and a large transcript database  <cit> .
like most other nightshades, tomato and potato both have a basic chromosome number of twelve, and there is genome-wide colinearity between their genomes  <cit> .
much effort is currently being invested to sequence the nuclear and organellar genomes of these organisms.
the international tomato genome sequencing project  <cit>  is sequencing the tomato  genome in the context of the family-wide solanaceae project .
rather than sequencing the complete genome, which is approximately  <dig> mb  <cit> , only the gene-rich euchromatic regions  are being sequenced using a bac-by-bac walking approach  <cit> .
the potato genome sequencing consortium   <cit>  aims to sequence the complete potato  genome of approximately  <dig> mb  <cit>  using a similar marker-anchored bac-by-bac sequencing strategy.
both sequencing projects rely heavily on bac libraries, of which three exist for tomato  and two exist for potato .
the tomato libraries are available through the sol genomics network   <cit>  and the potato libraries will soon by available at through the pgsc  <cit> .
all of these libraries have been end-sequenced to support bac-by-bac sequencing and extension, and to provide a base of genome-wide survey sequences to support studies such as the one presented here.
this paper describes the detailed sequence analysis of  <dig>  tomato bac end sequences , representing  <dig>  mb  of the tomato genome, and  <dig>  potato bess, corresponding to  <dig>  mb  of the potato genome .
this comparative genomics study aims to gain insight into the similarity between the tomato and potato genomes, both on the structural level through repeat and gene content analyses and on the functional level through gene function analyses.
furthermore, we investigate micro-syntenic relationships between these two solanaceous genomes, and several other sequenced plant genomes.
the sequence content of bess from a particular library is biased by which restriction enzyme was used to make the library.
to avoid comparing sequence sets with different biases, tomato-potato comparisons are made only between bess from libraries made with the same enzyme.
the sequences are subdivided into libraries, which are labeled with a three-letter code, with the corresponding restriction enzyme listed between brackets.
RESULTS
repeat density and categorization
based on similarity searches of the repeat database, between  <dig> % and  <dig> % of the nucleotides in the tomato bess were identified as belonging to a repeat .
the most common repeat families in the tomato libraries were the gypsy  and copia  classes of retrotransposons.
another prominent class of repeats comprised the ribosomal rna genes .
the tomato eco  library had the lowest repeat density at  <dig> %, which can be attributed to a lower amount of gypsy retrotransposons .
the highest repeat content was found in the tomato mbo  library , more than a third of which  consisted of ribosomal rna genes.
note that, since the repeat detection was based on sequence similarity, different segments in a bes could be assigned to more than one repeat family.
as a result, the sum of the repeat content per repeat type can be slightly larger than the total repeat content.
numbers represent percentages of nucleotides that show similarity to a repeat of the indicated category.
an 'x' represents the absence of a repeat family; ' <dig> ' indicates that the repeat is present, but at a frequency lower than  <dig>  % of the nucleotides in the bess.
species names have been abbreviated as follows: tom.
: tomato; pot.
: potato.
in contrast to the tomato bess, only between  <dig> % and  <dig> % of the nucleotides in the potato bess showed similarity to known magnoliaphytae repeats .
as in tomato, the majority of the repeats were found in the gypsy  and copia  retrotransposon families, whereas the fraction of ribosomal rna genes was small .
potato appeared to contain approximately two times as many line and sine elements as tomato , although the absolute percentages were low.
furthermore, a higher percentage of class ii dna transposons was observed in potato , the majority of which could not be classified.
in agreement with the differences observed between the tomato hba  and eco libraries, the potato ppt  library had an overall lower repeat content than the pot  library, and more specifically, a lower amount of gypsy retrotransposons .
the ppt library was also enriched in ribosomal rna genes in comparison to the pot library , just as was found comparing the eco library to the hba library in tomato.
since similarity-based repeat detection can be limited by the size and diversity of the repeat database, a self-comparison of the bess was performed in order to estimate the redundancy within the bess.
even with the stringent requirement that at least 50% of a given query sequence match another bes with at least 90% identity,  <dig> % of the nucleotides in the tomato bess had a match to one or more other tomato bess, and  <dig> % matched five or more other bess.
the redundancy in the potato bess was lower than in tomato;  <dig> % of the nucleotides in the potato bess had a hit to at least one other potato bess, and  <dig> % had a hit to five or more bess.
this difference could not be attributed solely to the larger number of tomato bess, compared to the number of potato bess; a self-comparison of the tomato hba library, which is of approximately the same size as the potato pot and ppt libraries combined, showed that  <dig> % of the nucleotides in this library matched at least one other hba bes, and  <dig> % matched five or more other hba bess.
the percentage of nucleotides in both species that matched five or more other bess was only slightly higher than the findings from the repeatmasker analysis , suggesting that the repeat database used in this study was sufficient to detect the majority of highly abundant repeats in these species.
these findings also confirm the observation from the similarity-based repeat detection that the tomato bess are more repetitive than the potato bess.
simple sequence repeats
a total of  <dig>  ssrs with a motif length between one and five nt, and a total length of at least  <dig> nt were detected in the tomato bess, representing one ssr per  <dig>  kb of genomic sequence.
the term 'motif length' is used here to describe the length of the motif that is repeated in the ssr; for example, an atatat repeat has a motif length of two .
the most abundant motif length was five nucleotides , followed by motif lengths of two , four , three , and lastly one .
in potato,  <dig>  ssrs were found, out of which  <dig>   belonged to class i .
thus, the potato bess had one ssr per  <dig>  kb of genomic sequence, which is higher than that in tomato .
as in tomato, the most abundant motif length in the potato ssrs was five nucleotides .
however, the next most abundant length was three , followed by motif lengths of two , four  and one .
considering only the class i ssrs, the most abundant ssr motifs in tomato and potato were at/ta  and a/t .
in tomato, the di-nucleotide motifs ac/gt  and ag/ct  were the most abundant after these two, whereas in potato the mononucleotide c/g  and tri-nucleotide aat/att  and aag/ctt  occurred at the second, third and fourth highest frequency, respectively.
this suggests that the differences in primary motif frequencies between tomato and potato also hold when considering only class i ssrs.
gene content
in the tomato bess, the percentage of nucleotides that matched by at least one database sequence ranged from  <dig> % for the eco library, to  <dig> % for the mbo library.
figure  <dig> presents a breakdown of these blast hits into three main categories , based on the keyword filtering described in materials and methods.
each category was then subdivided into 'masked' and 'unmasked' subcategories, with 'masked' indicating an overlap with repetitive sequences identified by repeatmasker, and 'unmasked' indicating a lack of such overlap.
in this way, the blast and repeatmasker results were combined in order to generate the best possible estimation of the percentage of putative protein-coding nucleotides in the bess.
the 'coding' category represents the percentage of nucleotides that matched one or more database sequences, and were not identified as repetitive by the keyword filtering.
after removing the overlap with repeats identified by repeatmasker, the percentage of coding nucleotides in the three libraries ranged from  <dig> % for the mbo library to  <dig> % for the hba library .
the mbo library had the highest percentage of the three libraries in the 'coding masked' category, which is likely the result of the high number of ribosomal repeat sequences in this library that have escaped the keyword filtering.
the 'repeats' category contains the blast matches to transposon and other repeat related sequences.
in all three libraries, there was a considerable fraction of nucleotides that the keyword filtering assigned to the 'repeats' category but that did not overlap with the repeats identified by repeatmasker .
this fraction ranged from  <dig> % in the eco library to  <dig> % in the hba library and may represent a combination of repeats that were missed by repeatmasker and true protein-coding genes that were miss-classified by the keyword filtering.
the final category in figure  <dig>  'other', represents all non-transposon-related repetitive sequences that were identified by the keyword filtering .
in the potato pot and ppt libraries,  <dig>  and  <dig> % of the nucleotides matched the protein database, respectively.
while these numbers were slightly lower than those for the tomato hba and eco libraries , the percentage of nucleotides assigned to the 'coding' category  was larger than those of the corresponding tomato libraries , suggesting that potato may have a larger gene repertoire than tomato.
furthermore, the number of transposon regions and other repeat-related regions that was found in this comparison to the protein database was more than  <dig> -fold higher for tomato than for potato.
this is consistent with the difference in transposon content that was found in the repeat analysis.
for the potato bess,  <dig> %  and  <dig> %  of the nucleotides had match in the potato est database, which is in fairly good agreement with the tomato hba and eco comparisons versus the tomato database .
fewer matches in the potato bess were 'masked' than in tomato, confirming the observation from the blastx comparison to the protein database that the potato bess have more protein coding nucleotides and lower repeat content.
functional annotation
a total of  <dig>  go terms, out of which  <dig> unique terms, were assigned to the tomato hba bess based matches in the pfam database .
although there were more than half as many eco bess as hba bess, only  <dig>  go terms  were assigned to them.
in potato,  <dig>  terms  were assigned to the pot library, whereas only  <dig>  terms  were assigned to the ppt library.
comparing the go annotations of tomato to those of potato  resulted in  <dig> significantly overrepresented terms between the hindiii digested libraries  and nine significantly overrepresented terms between the ecori digested libraries .
in both species, many of the terms that were overrepresented in the hindiii libraries compared to their ecori counterparts were related to retrotransposon activity, such as dna binding , dna integration , rna-directed dna polymerase activity , and chromatin-related terms .
furthermore, many of these transposon-related terms were significantly overrepresented in tomato, compared to potato .
this is consistent with the findings from the repeatmasker and blast analyses discussed above.
surprisingly, some terms that were overrepresented in both the ecori digested libraries could be linked to transcription factor genes.
in tomato, zinc ion binding , dna-dependent regulation of transcription , and transcription factor activity  were overrepresented in the eco library.
the potato ppt library was enriched for zinc ion binding , nucleic acid binding , and transcription factor activity .
analysis of the protein families identified by panther revealed similar trends for the number of matches, both within and between the tomato and potato libraries .
in tomato,  <dig>  distinct families were found in the hba bess for a total of  <dig>  hits, and  <dig>  hits representing  <dig> families were found in the eco bess.
analysis of the potato pot library revealed  <dig> distinct panther families for a total of  <dig>  hits; however, only  <dig>  hits to  <dig> families were found in the ppt bess.
two and three panther families were found to be overrepresented in the tomato hba and eco libraries, compared to eleven and five overrepresented families in the potato pot and ppt libraries, respectively.
consistent with the greater abundance of gypsy retrotransposons in the hindiii libraries of both tomato and potato, the gag/pol/env polyprotein  panther family was found to be overrepresented in both hindiii libraries, compared to the corresponding ecori libraries.
furthermore, the gag-pol-related retrotransposon  panther family was relatively more abundant in the ecori libraries, which also agrees with the difference in the gypsy:copia ratio between the hindiii and ecori libraries .
both of these retrotransposon-related terms were found to be significantly  overrepresented in tomato when compared to potato.
in the tomato eco library, transcription-factor related terms such as zinc finger cchc domain contain protein , zinc finger protein  and mads box protein  were significantly overrepresented , confirming the results from the go analysis.
no transcription-factor related panther families were significantly overrepresented in the potato ppt library.
between tomato and potato, the majority of the overrepresented terms in potato corresponded to important biological and biochemical processes.
for example, zinc finger cchc domain containing proteins  and general transcription factor 2-related zinc finger proteins  occurred with a significantly  higher frequency in potato pot than in tomato hba; the latter was also overrepresented in the potato ppt library.
this was also reflected in the go annotation through terms such as nucleic acid binding  and zinc ion binding .
the overrepresentation of these terms relative to tomato suggests an expansion of transcription factors or other genes for dna binding proteins in the potato genome.
another example is the cytochrome p <dig> superfamily , which was also found in the go analysis through terms such as iron ion binding  and mono-oxygenase activity .
cytochrome p <dig> proteins play important roles in the biosynthesis of secondary metabolites, and the overrepresentation of these proteins in potato could indicate an expanded network of pathways that synthesize secondary metabolites in potato.
a final example involves the large family of plant-type serine-threonine protein kinases , which are known to play important roles in disease resistance in various plant species .
in the panther database, this family consists of  <dig> different subfamilies,  <dig> of which were found in the tomato and potato bess.
out of these  <dig> subfamilies,  <dig> were found only in tomato, and five were unique to potato.
most of the subfamilies that were found in both species were overrepresented in potato, such as lrr receptor-like kinases  and lrr transmembrane kinases .
several subfamilies occurred at a higher frequency in tomato, including serine/threonine specific receptor-like protein kinases  and pto-like kinases .
thus, while the complement of serine-threonine protein kinases in potato exceeds that of tomato, several of the subfamilies have expanded specifically in tomato.
this may reflect an adaptation for resistance to different pathogens, or a difference in the dominant mechanism of pathogen resistance between these species.
comparative genome mapping
out of the  <dig>  pairs of tomato bess that were compared to the a. thaliana genome,  <dig>  pairs had one or more matches.
these matches were divided into five categories, as is shown in the last five columns of table  <dig>  the 'single end' category represents the bac end pairs from which only one of the two sequences had a match to the a. thaliana genome, and contained the majority of the matches .
paired end matches, in which the bess from the same bac each had a match to a different chromosome, were assigned to the 'non-linear' category.
the 'gapped' category contained  <dig>  bac end pairs that matched to the same a. thaliana chromosome with a distance between the paired matches that was either smaller than  <dig> kb or larger than  <dig> kb.
the final two categories represented the bacs from which both end sequences were matched to the genome within a distance of  <dig> to  <dig> kb of each other, either in the correct orientation with respect to each other , or rearranged with respect to each other .
out of the  <dig>  tomato bes pairs that hit to the same a. thaliana chromosome, three pairs fell into the 'colinear' category, and one pair fell into the 'rearranged' category, suggesting the presence of four putative micro-syntenic regions between tomato and a. thaliana.
potato had  <dig>  pairs of bess, out of which  <dig> pairs were mapped to the a. thaliana genome, with both bess of the pair matching the same chromosome.
two potato bacs displayed putative microsynteny based on the end sequence matches, one of which was colinear, whereas the other represented a possible rearrangement.
in comparison to tomato, potato had very few bacs that fell into the 'gapped' category, although the smaller ppt library had more than five times as many sequences in this category as the pot library.
interestingly, the large majority of the tomato bacs that fell into this category was from the eco and mbo libraries .
the ecori and mboi digested libraries were found to contain a high fraction of ribosomal rna genes in the repeatmasker analysis, and indeed more than 80% of the sequences from these libraries that fell into the 'gapped' category contained ribosomal rna genes.
repeating the same analysis against the p. trichocarpa genome, only  <dig> of the tomato bes pairs matched with both ends to the same chromosome .
it should be noted here that p. trichocarpa has both a larger number of chromosomes than a. thaliana  and approximately twenty-two thousand additional contig sequences that have not yet been integrated into the chromosome pseudomolecules.
based on these numbers alone, one would expect a smaller number of paired bess to map to the same chromosome or contig sequence.
even so, p. trichocarpa displayed more regions of micro-synteny with tomato than a. thaliana:  <dig> pairs of bess mapped within a distance between  <dig> and  <dig> kb of the other bes in the pair.
more than two-thirds of these matches  showed colinearity between tomato and p. trichocarpa, whereas the remaining  <dig> hits represented rearrangements in their respective regions of micro-synteny.
consistent with the difference between the tomato – a. thaliana and tomato – p. trichocarpa mappings, a smaller number of potato bes pairs  could be mapped with both ends to the same chromosome in p. trichocarpa, than in a. thaliana.
of these, there were  <dig> regions of potential microsynteny, out of which  <dig> were colinear.
compared to tomato, the 'non-linear' and to a lesser extent the 'gapped' categories were underrepresented in potato.
again these differences seem to originate from the fact that many of the bess in the eco and mbo libraries contain ribosomal rna genes.
the majority of these sequences fell into the 'non-linear' category in the p. trichocarpa comparison, rather than the 'gapped' category as was the case with a. thaliana, due to the ribosomal rna genes being contained in some of the unassembled contig sequences rather than in the chromosomal pseudomolecules.
discussion
sequence properties
based on the differences between the libraries in both tomato and potato, it seems unlikely that any of these partial digestion-based libraries represents an unbiased cross section of the genome.
for example, in tomato the mbo library has a higher gc percentage than the hba and eco libraries.
this difference is likely caused by the length and gc content of the restriction sites that were targeted in the digestion of the genome: both the hindiii and ecori sites  have a length of six nucleotides and a gc content of  <dig> %, whereas the mboi site  has a length of four nucleotides and a gc content of 50%.
the consequences of this are clearly visible in the results of the gene and repeat content analyses presented in this paper: results differ markedly among libraries made with different enzymes.
however, we think it reasonable to assume that tomato and potato libraries derived from digestion with the same restriction enzyme would have similar sequence bias.
using this assumption, we strive to minimize any effect of sequence bias on our results by maintaining logical separation of bess from different libraries, and only directly comparing data for bess from libraries constructed with the same restriction enzymes.
the tomato bess  are shorter than the potato bess on average.
the difference in average sequence length between the tomato hindiii and ecori libraries and their potato counterparts is approximately  <dig> nt for both libraries and is most likely the result of a difference in sequencing quality and equipment.
however, we think it reasonable to assume that a difference in sequence length on this scale would not influence the results of the similarity-based analyses that have been performed in this study.
repeat density and categorization
both the tomato and potato libraries vary in total repeat content and in ratios between repeat types.
for example, ribosomal dna sequences are overrepresented in the tomato mbo and eco, and the potato ppt libraries, relative to the tomato hba and potato pot library, respectively.
this phenomenon was also observed in a study of zea mays bess  <cit> , where it was attributed to the presence of many mboi sites in the z. mays ribosomal dna cluster, compared to one ecori site, and no hindiii sites.
by similar reasoning, the under-representation of gypsy retrotransposons in the eco and ppt libraries might result from a lower frequency of ecori sites in this element compared to hindiii and mboi sites.
the discrepancy between the repeats identified by repeatmasker  and blastx  indicates the need for tomato- and potato-specific repeat databases.
a repeat database had previously been generated from the tomato bess , however comparing the tomato bess to this database using repeatmasker resulted in approximately 60% of the tomato bess being annotated as repetitive .
the majority of these repeats could however not be assigned to a known repeat family.
thus, while the findings in this paper may present an underestimation of the actual repeat content of the tomato and potato bess, the findings from the repeatmasker and blastx analyses both clearly suggest a higher repeat content in the tomato bess than in the potato bess.
a correlation between genome size and retrotransposon content has previously been identified in the brassicaceae  <cit> .
there, it was found that the retrotransposon content increases with genome size, from approximately  <dig> to 10% in a. thaliana , to 14% in brassica rapa , to 20% in b. olacerea .
comparing this to cereal crops such as oryza sativa  suggests that while the actual retrotransposon content in cereals is higher than in brassicaceae, the correlation with genome size may be universally present in plants.
the data presented in this research indicate that genome expansion in the solanaceae is also associated with retrotransposon amplification; potato  has an estimated retrotransposon content between  <dig>   and  <dig> % , whereas that of tomato  is notably higher .
the ratio between gypsy and copia retrotransposon sequences in the tomato bess is between 1: <dig> and 2: <dig>  whereas this ratio in the potato bess is between 2: <dig> and 3: <dig>  while this ratio clearly differs within each species between libraries generated with a different restriction enzyme, the difference in ratios between tomato and potato is observed in both the hindiii and the ecori digested libraries .
in a. thaliana  <cit> , b. rapa  <cit> , carica papaya  <cit>  and z. mays  <cit> , this ratio is approximately 1: <dig>  the tomato and potato genomes appear more similar to the o. sativa genome in this respect, where the gypsy to copia ratio was found to be around 2: <dig>  <cit> .
the difference in the gypsy:copia ratio between tomato and potato suggests that the retrotransposon amplification associated with the genome expansion in tomato is predominantly the result of additional copia copies.
simple sequence repeats
the most abundant ssrs in all size categories for both tomato and potato were at-rich.
this is consistent with findings in other plant species, such as a. thaliana  <cit> , b. rapa  <cit> , c. papaya  <cit> , glycine max  <cit> , and musa acuminata  <cit> .
in both potato and tomato, penta-nucleotide repeats are the most common form of ssrs, and aaaat is the predominant repeat motif.
this is in sharp contrast to previously studied plant species, in which di- and penta-nucleotide repeats generally occur least frequently  <cit> .
in many plant species, such as a. thaliana, b. rapa  <cit> , and o. sativa  <cit> , tri-nucleotide repeats are the most abundant microsatellites.
however, bes analysis of c. papaya  <cit> , g. max  <cit>  and m. acuminata  <cit>  suggests that di-nucleotide repeats are more common in these plant species.
thus, both tomato and potato display a unique distribution of microsatellite frequencies compared to other studied plant species.
the tomato bess have a higher fraction of di- and tetra-nucleotide repeats compared to the potato bess.
this may be because one or more of the tomato bac end libraries are enriched for bacs that are derived from centromeric regions in the tomato genome, as these regions have previously been found to be enriched for long, class i di- and tetra-nucleotide repeats  <cit> .
however, the relative enrichment for di- and tetra-nucleotide repeats in tomato compared to potato is observed in all three tomato libraries; this would only be compatible with the hypothesis of enrichment for centromeric regions if these regions contain more hindiii, ecori and mboi sites than average for the tomato genome.
gene content
after repeat masking and keyword filtering, the percentage of nucleotides in the potato pot and ppt bess that have a match in the non-redundant protein database is  <dig> - to  <dig> -fold that of the tomato hba and eco bess, respectively.
both the percentage of nucleotides and the number of bess having a hit to the protein database after repeat masking and keyword filtering are higher in potato  than in tomato , supporting the hypothesis that potato has more putative protein-coding regions than tomato.
in the blastn comparison of the bess to the ests, a similar discrepancy between potato and tomato was observed, with potato having a  <dig> - to  <dig> -fold higher est coverage than tomato.
furthermore, cross-comparisons of the tomato bess to the potato ests and vice versa confirmed that the difference in est coverage of the bess was not caused by a difference in number of unique transcripts between the tomato and potato est collections .
the difference between the blast comparisons to the protein and transcript databases may be attributed to the presence of full-length cdna sequences in the tomato transcript data, whereas these are not present in the potato data, resulting in an overrepresentation in the tomato bess for the interior regions of coding sequences.
even if one assumes that this more conservative lower bound is correct, the results still suggest that potato has a larger gene repertoire than tomato since the tomato genome is only approximately  <dig>  times larger than the potato genome.
in both tomato and potato, a smaller percentage of nucleotides show similarity to the est database than to the protein database, while the percentage of non-repetitive coding sequence in the est database comparison  is higher than that in the protein database comparison .
surprisingly, the majority of the matches to the protein and transcript databases do not overlap.
for example, in the tomato hba library,  <dig> % and  <dig> % of the nucleotides have a match in the est and protein databases, respectively, while only  <dig> % have a match in both.
similarly, for the potato pot library, only  <dig> % of the nucleotides have a match in both the transcript and protein sequences, whereas the individual percentages of nucleotides that have a match in these databases are  <dig> % and  <dig> %, respectively.
on one hand, the matches to the est databases that do not overlap with matches to the protein database may represent unique, taxon- or species-specific protein-coding genes that are not represented in the non-redundant protein database, or transcribed but untranslated regions in these genomes.
on the other hand, matches to the protein database that do not overlap with matches in the est database may indicate either the presence of genes that were not sufficiently expressed in the tissues under the conditions that were sampled during est library construction, or mis-annotated or otherwise incorrect sequences in the protein database.
the est data likely provides the most reliable sampling of the true protein coding regions in these genomes, since it is based on experimental data that contain species-specific sequences not available in the protein database.
due to the selection for poly-a tails that is normally used in the construction of est libraries, the number of non-protein coding transcripts will be relatively small.
taking the nucleotides from the hba and eco libraries that match ests and do not overlap with repeats as a measure of coding sequences, the tomato genome  is estimated to contain between  <dig>  and  <dig>  mb of coding regions.
similarly, assuming a genome size of  <dig> mb, the total coding region length for potato would be between  <dig>  and  <dig>  mb.
these numbers set lower bounds on the estimated coding content of these genomes, as the est data is unlikely to represent the full complement of full-length protein-coding sequences in these genomes.
previous estimates put the gene content of tomato at  <dig>  genes, based on an analysis of  <dig>  unigenes and  <dig> bac sequences  <cit> .
if these  <dig>  genes are indeed represented by  <dig>  mb of coding sequence , then the average transcript length of tomato would be approximately  <dig>  kb.
this is longer than the average transcript length in a. thaliana, which is  <dig>  kb according to the tair <dig> genome annotation  <cit> .
assuming the same average transcript length, potato  would contain approximately  <dig>  genes, or  <dig>  more genes than tomato.
since the data presented here are based on similarity-searches on short genomic sequences only, this difference does not necessarily represent a difference in functional genes, but may also reflect a larger proportion of pseudogenes or otherwise non-functional alleles in potato.
functional annotation
the results from the go and panther analysis generally show a similar trend.
the tomato bess have more go terms and panther families associated to them than the potato bess do.
however, the potato bess have a larger number of unique terms associated to them.
this agrees with the results of the blastx comparison to the non-redundant protein database, where it was found that the tomato bess have a higher overall coverage of blast hits, but a lower percentage of putative protein coding regions .
in both the go term and panther family analyses, the majority of the terms occur at a relatively low frequency.
for example, in the tomato hba versus potato pot comparison, only  <dig> out of the  <dig> distinct go terms that were assigned to the bess occurred ten or more times in at least one of the species.
this group of  <dig> go terms contained all  <dig> of the terms that were significantly  overrepresented in one of the species in this comparison.
moreover,  <dig> out of these  <dig> terms were found at least  <dig> times in at least one species, and this subgroup contained  <dig> out of the  <dig> significantly overrepresented terms.
similarly, in the panther family analysis,  <dig> out of the  <dig>  distinct families that were found in the bess occurred at least ten times in at least one species, out of which  <dig> families were found at least  <dig> times.
the  <dig> families that were found at least ten times contained every one of the  <dig> families that were significantly overrepresented in one of the species; ten of these were counted more than  <dig> times in at least one species.
while only the tomato hba versus potato pot comparison is shown here, the other comparisons show a similar pattern, indicating that many of the highly abundant go terms and panther families are significantly overrepresented in either tomato or potato.
the majority of these overrepresented terms and families are most abundant in potato, and represent biologically important functions and processes.
in tomato, a smaller number of terms and families is overrepresented; these are primarily connected to structural genomic features such as retrotransposons.
the overrepresentation of transposon-related go terms and panther families in tomato was consistent with the results of the repeat analysis, confirming the observations that tomato is richer in retrotransposons than potato.
however, in the panther analysis, reverse transcriptases  were significantly overrepresented in potato.
at first glance, this does not correspond well with the overrepresentation of rna-directed dna polymerase activity  and rna-dependent dna replication  in tomato.
however, in both tomato and potato, the large majority of the reverse transcriptases originated from non-ltr retroelements , which in fact is consistent with the higher frequency of non-ltr retrotransposons in potato found in the repeatmasker analysis .
the cytochrome p <dig> mono-oxygenases represent a large gene superfamily in plants that are commonly associated with the biosynthesis of secondary metabolites.
in a. thaliana, at least  <dig> p <dig> genes have been found, representing approximately one percent of the gene complement of this species.
in o. sativa, this family is even larger, with  <dig> p <dig> genes identified so far  <cit> .
not all the p450s in these genomes represent true protein-coding sequences; in a. thaliana, 90% of the genes are truly protein coding, compared to 72% in o. sativa.
in total,  <dig> distinct families of p <dig> genes were identified in a. thaliana and o. sativa, several of which were found to be overrepresented in either of these species.
moreover, some families were present in one, but completely absent in the other species  <cit> .
in the hindiii and ecori libraries,  <dig> and  <dig> bess that could be associated with the cytochrome p <dig> panther family  were found in tomato and potato, respectively.
since these bac end sequences represent approximately 14% and 10% of their respective genomes, these data suggest an enormous expansion of p <dig> genes in the solanaceae.
this could be the result of an expansion of specific p <dig> families, but also of the evolution of species- or family-specific p450s.
for example, the allene oxide synthase has currently only been found in solanaceous species, including tomato and petunia inflata  <cit> .
the overrepresentation of p450s in potato compared to tomato may be another result of species-specific p <dig> families, but may also indicate expansion of families that are shared between these species.
comparative genome mapping
in this study, paired bac ends have been exploited to detect regions of microsynteny between the solanaceous species tomato and potato, and the model plant organisms a. thaliana and p. trichocarpa.
using similar approaches, microsynteny has been observed between a. thaliana and b. rapa  <cit> ; c. papaya and p. trichocarpa  <cit> ; and m. acuminata and o. sativa  <cit> .
a higher number of tomato and potato bacs display microsynteny to p. trichocarpa, than to a. thaliana.
the reduced level of microsynteny between tomato/potato and a. thaliana is not likely a difference in evolutionary distances between these species.
both a. thaliana and p. trichocarpa are part of the rosids clade, whereas tomato and potato belong to the asterids clade.
it may be the result of a recent duplication of the a. thaliana genome, followed by the loss of roughly 70% of the duplicated genes  <cit> .
assuming that this loss occurred randomly, the large majority of possible microsyntenic regions that existed before the duplication will have disappeared due to the major local expansions and contractions that would be associated with such a duplication and subsequent loss.
this hypothesis is strengthened by the observation that only approximately 1% of  <dig>  a. thaliana bes pairs could previously be mapped within  <dig> kb to the p. trichocarpa genome, indicating that the organization of these genomes is indeed vastly different  <cit> .
regions of microsynteny have previously been detected between tomato/potato and a. thaliana.
a  <dig> kb region of tomato chromosome  <dig> containing five genes was shown to be syntenic with a  <dig> kb region on a. thaliana chromosome  <dig>  although the order and orientation of the genes suggested two inversion events  <cit> .
in another study, a  <dig> kb bac sequence matched to four different segments on a. thaliana chromosomes  <dig>   <dig>   <dig>  and 5; however, each of the four regions in a. thaliana were shorter than their tomato counterpart  <cit> .
recently, five microsyntenic blocks were detected between a region on potato chromosome  <dig> harbouring a qtl for resistance to late blight and root cyst nematodes, and a. thaliana chromosomes  <dig>   <dig> and  <dig>  <cit> .
these syntenic blocks spanned between three and seven orfs, and were interrupted by non-syntenic blocks.
in each of these examples, the microsynteny between tomato/potato and a. thaliana involves shorter regions on the a. thaliana genome than the average tomato and potato bac sequence length.
furthermore, regions of synteny are often detected between coding sequences, whereas the fraction of coding sequences in the tomato and potato bess is relatively low , providing a good explanation for the reduced amount of microsynteny between these species observed here.
synteny between potato and a. thaliana has also been identified on a genome-wide level using a comparative mapping approach.
this revealed  <dig> putative syntenic blocks between potato and a. thaliana that cover 41% of the potato genetic map, and 50% of the a. thaliana physical map  <cit> .
these syntenic blocks were unevenly distributed over the potato genetic map, and redundant in respect to the number of areas on the a. thaliana genome that displayed synteny to most areas on the potato map.
the regions of microsynteny between tomato/potato and a. thaliana that were found with the bes-based approach described in this study could not be used to confirm or renounce any putative higher-order syntenic regions, due to the relatively short distances between the bac ends.
six paired tomato bac end matches cluster in the  <dig>  –  <dig>  mb interval of p. trichocarpa chromosome  <dig>  furthermore, seven pairs of potato bess map to the partially overlapping interval between  <dig>  –  <dig>  mb, indicating the presence of either a number of distinct microsyntenic regions, or possibly a single region of macrosynteny, between the tomato/potato and p. trichocarpa genomes.
these findings provide an interesting starting point for a detailed comparison between these species in this region, once more tomato and potato genomic sequences become available.
CONCLUSIONS
the large scale analysis of tomato and potato bess presented in this paper revealed many interesting structural and functional differences between the genomes of these closely related species.
we have shown that the tomato genome is not only more repetitive than the potato genome, but that these genomes also differ in their repeat composition, most importantly in the distribution of gypsy and copia retrotransposons.
in contrast to other studied plant genomes, we have shown that the tomato and potato genomes contain a large number of ssrs with a motif length of five, which may be a unique feature of solanaceous genomes.
comparative analysis of the putative protein coding regions in these bess revealed an enrichment of these regions in the potato genome.
moreover, several protein families were found to be overrepresented in potato compared to tomato, such as cytochrome p <dig> mono-oxygenases and serine-threonine protein kinases.
the p <dig> superfamily appears to have expanded dramatically in both species compared to a. thaliana, suggesting an expanded network of secondary metabolic pathways in the solanaceae.
both tomato and potato appear to have low microsynteny with a. thaliana, which is likely a result of this species' relatively recent genomic rearrangement.
a higher degree of synteny was observed with p. trichocarpa.
difference in evolutionary distances is not likely to be the reason for this increased microsynteny, since both a. thaliana and p. trichocarpa are part of the rosids clade, whereas tomato and potato belong to the asterids clade.
taken together, these findings present a first glimpse into the evolution of solanaceous genomes, both within the family and relative to other plant species.
when the complete genomic sequences of these species become available, whole-genome comparisons and protein- or repeat-family specific studies may shed more light on the intriguing observations made in this paper.

 the solanaceae, or nightshade family, is a dicot plant family that includes many economically important genera that are used in agriculture, horticulture, and other industries. family members include the tuber bearing potato; a large number of fruit-bearing vegetables, such as peppers, tomatoes, and eggplant; leafy tobacco; and ornamental flowers from the petunia and solanum genera. tomato is generally considered to be a model crop plant species, for which many high-quality genetic and genomic resources are available, such as high-density molecular maps, many well-characterized near-isogenic lines, and rich collections of ests and full-length cdnas. potato is the most important crop within the solanaceae, ranking fourth as a world food crop following wheat, maize and rice. similar resources are available for potato, including an ultra-high density linkage map, a collection of phenotype data, and a large transcript database. like most other nightshades, tomato and potato both have a basic chromosome number of twelve, and there is genome-wide colinearity between their genomes. much effort is currently being invested to sequence the nuclear and organellar genomes of these organisms. the international tomato genome sequencing project is sequencing the tomato genome in the context of the family-wide solanaceae project. rather than sequencing the complete genome, which is approximately mb, only the gene-rich euchromatic regions are being sequenced using a bac-by-bac walking approach. the potato genome sequencing consortium aims to sequence the complete potato genome of approximately mb using a similar marker-anchored bac-by-bac sequencing strategy. both sequencing projects rely heavily on bac libraries, of which three exist for tomato and two exist for potato. the tomato libraries are available through the sol genomics network and the potato libraries will soon by available at through the pgsc. all of these libraries have been end-sequenced to support bac-by-bac sequencing and extension, and to provide a base of genome-wide survey sequences to support studies such as the one presented here. this paper describes the detailed sequence analysis of tomato bac end sequences, representing mb of the tomato genome, and potato bess, corresponding to mb of the potato genome. this comparative genomics study aims to gain insight into the similarity between the tomato and potato genomes, both on the structural level through repeat and gene content analyses and on the functional level through gene function analyses. furthermore, we investigate micro-syntenic relationships between these two solanaceous genomes, and several other sequenced plant genomes. the sequence content of bess from a particular library is biased by which restriction enzyme was used to make the library. to avoid comparing sequence sets with different biases, tomato-potato comparisons are made only between bess from libraries made with the same enzyme. the sequences are subdivided into libraries, which are labeled with a three-letter code, with the corresponding restriction enzyme listed between brackets.  repeat density and categorization based on similarity searches of the repeat database, between % and % of the nucleotides in the tomato bess were identified as belonging to a repeat. the most common repeat families in the tomato libraries were the gypsy and copia classes of retrotransposons. another prominent class of repeats comprised the ribosomal rna genes. the tomato eco library had the lowest repeat density at %, which can be attributed to a lower amount of gypsy retrotransposons. the highest repeat content was found in the tomato mbo library, more than a third of which consisted of ribosomal rna genes. note that, since the repeat detection was based on sequence similarity, different segments in a bes could be assigned to more than one repeat family. as a result, the sum of the repeat content per repeat type can be slightly larger than the total repeat content. numbers represent percentages of nucleotides that show similarity to a repeat of the indicated category. an'x' represents the absence of a repeat family;'' indicates that the repeat is present, but at a frequency lower than % of the nucleotides in the bess. species names have been abbreviated as follows: tom.: tomato; pot.: potato. in contrast to the tomato bess, only between % and % of the nucleotides in the potato bess showed similarity to known magnoliaphytae repeats. as in tomato, the majority of the repeats were found in the gypsy and copia retrotransposon families, whereas the fraction of ribosomal rna genes was small. potato appeared to contain approximately two times as many line and sine elements as tomato, although the absolute percentages were low. furthermore, a higher percentage of class ii dna transposons was observed in potato, the majority of which could not be classified. in agreement with the differences observed between the tomato hba and eco libraries, the potato ppt library had an overall lower repeat content than the pot library, and more specifically, a lower amount of gypsy retrotransposons. the ppt library was also enriched in ribosomal rna genes in comparison to the pot library, just as was found comparing the eco library to the hba library in tomato. since similarity-based repeat detection can be limited by the size and diversity of the repeat database, a self-comparison of the bess was performed in order to estimate the redundancy within the bess. even with the stringent requirement that at least 50% of a given query sequence match another bes with at least 90% identity, % of the nucleotides in the tomato bess had a match to one or more other tomato bess, and % matched five or more other bess. the redundancy in the potato bess was lower than in tomato; % of the nucleotides in the potato bess had a hit to at least one other potato bess, and % had a hit to five or more bess. this difference could not be attributed solely to the larger number of tomato bess, compared to the number of potato bess; a self-comparison of the tomato hba library, which is of approximately the same size as the potato pot and ppt libraries combined, showed that % of the nucleotides in this library matched at least one other hba bes, and % matched five or more other hba bess. the percentage of nucleotides in both species that matched five or more other bess was only slightly higher than the findings from the repeatmasker analysis, suggesting that the repeat database used in this study was sufficient to detect the majority of highly abundant repeats in these species. these findings also confirm the observation from the similarity-based repeat detection that the tomato bess are more repetitive than the potato bess. simple sequence repeats a total of ssrs with a motif length between one and five nt, and a total length of at least nt were detected in the tomato bess, representing one ssr per kb of genomic sequence. the term'motif length' is used here to describe the length of the motif that is repeated in the ssr; for example, an atatat repeat has a motif length of two. the most abundant motif length was five nucleotides, followed by motif lengths of two, four, three, and lastly one. in potato, ssrs were found, out of which belonged to class i. thus, the potato bess had one ssr per kb of genomic sequence, which is higher than that in tomato. as in tomato, the most abundant motif length in the potato ssrs was five nucleotides. however, the next most abundant length was three, followed by motif lengths of two, four and one. considering only the class i ssrs, the most abundant ssr motifs in tomato and potato were at/ta and a/t. in tomato, the di-nucleotide motifs ac/gt and ag/ct were the most abundant after these two, whereas in potato the mononucleotide c/g and tri-nucleotide aat/att and aag/ctt occurred at the second, third and fourth highest frequency, respectively. this suggests that the differences in primary motif frequencies between tomato and potato also hold when considering only class i ssrs. gene content in the tomato bess, the percentage of nucleotides that matched by at least one database sequence ranged from % for the eco library, to % for the mbo library. figure presents a breakdown of these blast hits into three main categories, based on the keyword filtering described in materials and methods. each category was then subdivided into'masked' and'unmasked' subcategories, with'masked' indicating an overlap with repetitive sequences identified by repeatmasker, and'unmasked' indicating a lack of such overlap. in this way, the blast and repeatmasker results were combined in order to generate the best possible estimation of the percentage of putative protein-coding nucleotides in the bess. the'coding' category represents the percentage of nucleotides that matched one or more database sequences, and were not identified as repetitive by the keyword filtering. after removing the overlap with repeats identified by repeatmasker, the percentage of coding nucleotides in the three libraries ranged from % for the mbo library to % for the hba library. the mbo library had the highest percentage of the three libraries in the'coding masked' category, which is likely the result of the high number of ribosomal repeat sequences in this library that have escaped the keyword filtering. the'repeats' category contains the blast matches to transposon and other repeat related sequences. in all three libraries, there was a considerable fraction of nucleotides that the keyword filtering assigned to the'repeats' category but that did not overlap with the repeats identified by repeatmasker. this fraction ranged from % in the eco library to % in the hba library and may represent a combination of repeats that were missed by repeatmasker and true protein-coding genes that were miss-classified by the keyword filtering. the final category in figure'other', represents all non-transposon-related repetitive sequences that were identified by the keyword filtering. in the potato pot and ppt libraries, and % of the nucleotides matched the protein database, respectively. while these numbers were slightly lower than those for the tomato hba and eco libraries, the percentage of nucleotides assigned to the'coding' category was larger than those of the corresponding tomato libraries, suggesting that potato may have a larger gene repertoire than tomato. furthermore, the number of transposon regions and other repeat-related regions that was found in this comparison to the protein database was more than -fold higher for tomato than for potato. this is consistent with the difference in transposon content that was found in the repeat analysis. for the potato bess, % and % of the nucleotides had match in the potato est database, which is in fairly good agreement with the tomato hba and eco comparisons versus the tomato database. fewer matches in the potato bess were'masked' than in tomato, confirming the observation from the blastx comparison to the protein database that the potato bess have more protein coding nucleotides and lower repeat content. functional annotation a total of go terms, out of which unique terms, were assigned to the tomato hba bess based matches in the pfam database. although there were more than half as many eco bess as hba bess, only go terms were assigned to them. in potato, terms were assigned to the pot library, whereas only terms were assigned to the ppt library. comparing the go annotations of tomato to those of potato resulted in significantly overrepresented terms between the hindiii digested libraries and nine significantly overrepresented terms between the ecori digested libraries. in both species, many of the terms that were overrepresented in the hindiii libraries compared to their ecori counterparts were related to retrotransposon activity, such as dna binding, dna integration, rna-directed dna polymerase activity, and chromatin-related terms. furthermore, many of these transposon-related terms were significantly overrepresented in tomato, compared to potato. this is consistent with the findings from the repeatmasker and blast analyses discussed above. surprisingly, some terms that were overrepresented in both the ecori digested libraries could be linked to transcription factor genes. in tomato, zinc ion binding, dna-dependent regulation of transcription, and transcription factor activity were overrepresented in the eco library. the potato ppt library was enriched for zinc ion binding, nucleic acid binding, and transcription factor activity. analysis of the protein families identified by panther revealed similar trends for the number of matches, both within and between the tomato and potato libraries. in tomato, distinct families were found in the hba bess for a total of hits, and hits representing families were found in the eco bess. analysis of the potato pot library revealed distinct panther families for a total of hits; however, only hits to families were found in the ppt bess. two and three panther families were found to be overrepresented in the tomato hba and eco libraries, compared to eleven and five overrepresented families in the potato pot and ppt libraries, respectively. consistent with the greater abundance of gypsy retrotransposons in the hindiii libraries of both tomato and potato, the gag/pol/env polyprotein panther family was found to be overrepresented in both hindiii libraries, compared to the corresponding ecori libraries. furthermore, the gag-pol-related retrotransposon panther family was relatively more abundant in the ecori libraries, which also agrees with the difference in the gypsy:copia ratio between the hindiii and ecori libraries. both of these retrotransposon-related terms were found to be significantly overrepresented in tomato when compared to potato. in the tomato eco library, transcription-factor related terms such as zinc finger cchc domain contain protein, zinc finger protein and mads box protein were significantly overrepresented, confirming the results from the go analysis. no transcription-factor related panther families were significantly overrepresented in the potato ppt library. between tomato and potato, the majority of the overrepresented terms in potato corresponded to important biological and biochemical processes. for example, zinc finger cchc domain containing proteins and general transcription factor 2-related zinc finger proteins occurred with a significantly higher frequency in potato pot than in tomato hba; the latter was also overrepresented in the potato ppt library. this was also reflected in the go annotation through terms such as nucleic acid binding and zinc ion binding. the overrepresentation of these terms relative to tomato suggests an expansion of transcription factors or other genes for dna binding proteins in the potato genome. another example is the cytochrome p superfamily, which was also found in the go analysis through terms such as iron ion binding and mono-oxygenase activity. cytochrome p proteins play important roles in the biosynthesis of secondary metabolites, and the overrepresentation of these proteins in potato could indicate an expanded network of pathways that synthesize secondary metabolites in potato. a final example involves the large family of plant-type serine-threonine protein kinases, which are known to play important roles in disease resistance in various plant species. in the panther database, this family consists of different subfamilies, of which were found in the tomato and potato bess. out of these subfamilies, were found only in tomato, and five were unique to potato. most of the subfamilies that were found in both species were overrepresented in potato, such as lrr receptor-like kinases and lrr transmembrane kinases. several subfamilies occurred at a higher frequency in tomato, including serine/threonine specific receptor-like protein kinases and pto-like kinases. thus, while the complement of serine-threonine protein kinases in potato exceeds that of tomato, several of the subfamilies have expanded specifically in tomato. this may reflect an adaptation for resistance to different pathogens, or a difference in the dominant mechanism of pathogen resistance between these species. comparative genome mapping out of the pairs of tomato bess that were compared to the a. thaliana genome, pairs had one or more matches. these matches were divided into five categories, as is shown in the last five columns of table the'single end' category represents the bac end pairs from which only one of the two sequences had a match to the a. thaliana genome, and contained the majority of the matches. paired end matches, in which the bess from the same bac each had a match to a different chromosome, were assigned to the'non-linear' category. the'gapped' category contained bac end pairs that matched to the same a. thaliana chromosome with a distance between the paired matches that was either smaller than kb or larger than kb. the final two categories represented the bacs from which both end sequences were matched to the genome within a distance of to kb of each other, either in the correct orientation with respect to each other, or rearranged with respect to each other. out of the tomato bes pairs that hit to the same a. thaliana chromosome, three pairs fell into the'colinear' category, and one pair fell into the'rearranged' category, suggesting the presence of four putative micro-syntenic regions between tomato and a. thaliana. potato had pairs of bess, out of which pairs were mapped to the a. thaliana genome, with both bess of the pair matching the same chromosome. two potato bacs displayed putative microsynteny based on the end sequence matches, one of which was colinear, whereas the other represented a possible rearrangement. in comparison to tomato, potato had very few bacs that fell into the'gapped' category, although the smaller ppt library had more than five times as many sequences in this category as the pot library. interestingly, the large majority of the tomato bacs that fell into this category was from the eco and mbo libraries. the ecori and mboi digested libraries were found to contain a high fraction of ribosomal rna genes in the repeatmasker analysis, and indeed more than 80% of the sequences from these libraries that fell into the'gapped' category contained ribosomal rna genes. repeating the same analysis against the p. trichocarpa genome, only of the tomato bes pairs matched with both ends to the same chromosome. it should be noted here that p. trichocarpa has both a larger number of chromosomes than a. thaliana and approximately twenty-two thousand additional contig sequences that have not yet been integrated into the chromosome pseudomolecules. based on these numbers alone, one would expect a smaller number of paired bess to map to the same chromosome or contig sequence. even so, p. trichocarpa displayed more regions of micro-synteny with tomato than a. thaliana: pairs of bess mapped within a distance between and kb of the other bes in the pair. more than two-thirds of these matches showed colinearity between tomato and p. trichocarpa, whereas the remaining hits represented rearrangements in their respective regions of micro-synteny. consistent with the difference between the tomato – a. thaliana and tomato – p. trichocarpa mappings, a smaller number of potato bes pairs could be mapped with both ends to the same chromosome in p. trichocarpa, than in a. thaliana. of these, there were regions of potential microsynteny, out of which were colinear. compared to tomato, the'non-linear' and to a lesser extent the'gapped' categories were underrepresented in potato. again these differences seem to originate from the fact that many of the bess in the eco and mbo libraries contain ribosomal rna genes. the majority of these sequences fell into the'non-linear' category in the p. trichocarpa comparison, rather than the'gapped' category as was the case with a. thaliana, due to the ribosomal rna genes being contained in some of the unassembled contig sequences rather than in the chromosomal pseudomolecules. discussion sequence properties based on the differences between the libraries in both tomato and potato, it seems unlikely that any of these partial digestion-based libraries represents an unbiased cross section of the genome. for example, in tomato the mbo library has a higher gc percentage than the hba and eco libraries. this difference is likely caused by the length and gc content of the restriction sites that were targeted in the digestion of the genome: both the hindiii and ecori sites have a length of six nucleotides and a gc content of %, whereas the mboi site has a length of four nucleotides and a gc content of 50%. the consequences of this are clearly visible in the results of the gene and repeat content analyses presented in this paper: results differ markedly among libraries made with different enzymes. however, we think it reasonable to assume that tomato and potato libraries derived from digestion with the same restriction enzyme would have similar sequence bias. using this assumption, we strive to minimize any effect of sequence bias on our results by maintaining logical separation of bess from different libraries, and only directly comparing data for bess from libraries constructed with the same restriction enzymes. the tomato bess are shorter than the potato bess on average. the difference in average sequence length between the tomato hindiii and ecori libraries and their potato counterparts is approximately nt for both libraries and is most likely the result of a difference in sequencing quality and equipment. however, we think it reasonable to assume that a difference in sequence length on this scale would not influence the results of the similarity-based analyses that have been performed in this study. repeat density and categorization both the tomato and potato libraries vary in total repeat content and in ratios between repeat types. for example, ribosomal dna sequences are overrepresented in the tomato mbo and eco, and the potato ppt libraries, relative to the tomato hba and potato pot library, respectively. this phenomenon was also observed in a study of zea mays bess, where it was attributed to the presence of many mboi sites in the z. mays ribosomal dna cluster, compared to one ecori site, and no hindiii sites. by similar reasoning, the under-representation of gypsy retrotransposons in the eco and ppt libraries might result from a lower frequency of ecori sites in this element compared to hindiii and mboi sites. the discrepancy between the repeats identified by repeatmasker and blastx indicates the need for tomato- and potato-specific repeat databases. a repeat database had previously been generated from the tomato bess, however comparing the tomato bess to this database using repeatmasker resulted in approximately 60% of the tomato bess being annotated as repetitive. the majority of these repeats could however not be assigned to a known repeat family. thus, while the findings in this paper may present an underestimation of the actual repeat content of the tomato and potato bess, the findings from the repeatmasker and blastx analyses both clearly suggest a higher repeat content in the tomato bess than in the potato bess. a correlation between genome size and retrotransposon content has previously been identified in the brassicaceae. there, it was found that the retrotransposon content increases with genome size, from approximately to 10% in a. thaliana, to 14% in brassica rapa, to 20% in b. olacerea. comparing this to cereal crops such as oryza sativa suggests that while the actual retrotransposon content in cereals is higher than in brassicaceae, the correlation with genome size may be universally present in plants. the data presented in this research indicate that genome expansion in the solanaceae is also associated with retrotransposon amplification; potato has an estimated retrotransposon content between and %, whereas that of tomato is notably higher. the ratio between gypsy and copia retrotransposon sequences in the tomato bess is between 1: and 2: whereas this ratio in the potato bess is between 2: and 3: while this ratio clearly differs within each species between libraries generated with a different restriction enzyme, the difference in ratios between tomato and potato is observed in both the hindiii and the ecori digested libraries. in a. thaliana, b. rapa, carica papaya and z. mays, this ratio is approximately 1: the tomato and potato genomes appear more similar to the o. sativa genome in this respect, where the gypsy to copia ratio was found to be around 2:. the difference in the gypsy:copia ratio between tomato and potato suggests that the retrotransposon amplification associated with the genome expansion in tomato is predominantly the result of additional copia copies. simple sequence repeats the most abundant ssrs in all size categories for both tomato and potato were at-rich. this is consistent with findings in other plant species, such as a. thaliana, b. rapa, c. papaya, glycine max, and musa acuminata. in both potato and tomato, penta-nucleotide repeats are the most common form of ssrs, and aaaat is the predominant repeat motif. this is in sharp contrast to previously studied plant species, in which di- and penta-nucleotide repeats generally occur least frequently. in many plant species, such as a. thaliana, b. rapa, and o. sativa, tri-nucleotide repeats are the most abundant microsatellites. however, bes analysis of c. papaya, g. max and m. acuminata suggests that di-nucleotide repeats are more common in these plant species. thus, both tomato and potato display a unique distribution of microsatellite frequencies compared to other studied plant species. the tomato bess have a higher fraction of di- and tetra-nucleotide repeats compared to the potato bess. this may be because one or more of the tomato bac end libraries are enriched for bacs that are derived from centromeric regions in the tomato genome, as these regions have previously been found to be enriched for long, class i di- and tetra-nucleotide repeats. however, the relative enrichment for di- and tetra-nucleotide repeats in tomato compared to potato is observed in all three tomato libraries; this would only be compatible with the hypothesis of enrichment for centromeric regions if these regions contain more hindiii, ecori and mboi sites than average for the tomato genome. gene content after repeat masking and keyword filtering, the percentage of nucleotides in the potato pot and ppt bess that have a match in the non-redundant protein database is - to -fold that of the tomato hba and eco bess, respectively. both the percentage of nucleotides and the number of bess having a hit to the protein database after repeat masking and keyword filtering are higher in potato than in tomato, supporting the hypothesis that potato has more putative protein-coding regions than tomato. in the blastn comparison of the bess to the ests, a similar discrepancy between potato and tomato was observed, with potato having a - to -fold higher est coverage than tomato. furthermore, cross-comparisons of the tomato bess to the potato ests and vice versa confirmed that the difference in est coverage of the bess was not caused by a difference in number of unique transcripts between the tomato and potato est collections. the difference between the blast comparisons to the protein and transcript databases may be attributed to the presence of full-length cdna sequences in the tomato transcript data, whereas these are not present in the potato data, resulting in an overrepresentation in the tomato bess for the interior regions of coding sequences. even if one assumes that this more conservative lower bound is correct, the results still suggest that potato has a larger gene repertoire than tomato since the tomato genome is only approximately times larger than the potato genome. in both tomato and potato, a smaller percentage of nucleotides show similarity to the est database than to the protein database, while the percentage of non-repetitive coding sequence in the est database comparison is higher than that in the protein database comparison. surprisingly, the majority of the matches to the protein and transcript databases do not overlap. for example, in the tomato hba library, % and % of the nucleotides have a match in the est and protein databases, respectively, while only % have a match in both. similarly, for the potato pot library, only % of the nucleotides have a match in both the transcript and protein sequences, whereas the individual percentages of nucleotides that have a match in these databases are % and %, respectively. on one hand, the matches to the est databases that do not overlap with matches to the protein database may represent unique, taxon- or species-specific protein-coding genes that are not represented in the non-redundant protein database, or transcribed but untranslated regions in these genomes. on the other hand, matches to the protein database that do not overlap with matches in the est database may indicate either the presence of genes that were not sufficiently expressed in the tissues under the conditions that were sampled during est library construction, or mis-annotated or otherwise incorrect sequences in the protein database. the est data likely provides the most reliable sampling of the true protein coding regions in these genomes, since it is based on experimental data that contain species-specific sequences not available in the protein database. due to the selection for poly-a tails that is normally used in the construction of est libraries, the number of non-protein coding transcripts will be relatively small. taking the nucleotides from the hba and eco libraries that match ests and do not overlap with repeats as a measure of coding sequences, the tomato genome is estimated to contain between and mb of coding regions. similarly, assuming a genome size of mb, the total coding region length for potato would be between and mb. these numbers set lower bounds on the estimated coding content of these genomes, as the est data is unlikely to represent the full complement of full-length protein-coding sequences in these genomes. previous estimates put the gene content of tomato at genes, based on an analysis of unigenes and bac sequences. if these genes are indeed represented by mb of coding sequence, then the average transcript length of tomato would be approximately kb. this is longer than the average transcript length in a. thaliana, which is kb according to the tair genome annotation. assuming the same average transcript length, potato would contain approximately genes, or more genes than tomato. since the data presented here are based on similarity-searches on short genomic sequences only, this difference does not necessarily represent a difference in functional genes, but may also reflect a larger proportion of pseudogenes or otherwise non-functional alleles in potato. functional annotation the results from the go and panther analysis generally show a similar trend. the tomato bess have more go terms and panther families associated to them than the potato bess do. however, the potato bess have a larger number of unique terms associated to them. this agrees with the results of the blastx comparison to the non-redundant protein database, where it was found that the tomato bess have a higher overall coverage of blast hits, but a lower percentage of putative protein coding regions. in both the go term and panther family analyses, the majority of the terms occur at a relatively low frequency. for example, in the tomato hba versus potato pot comparison, only out of the distinct go terms that were assigned to the bess occurred ten or more times in at least one of the species. this group of go terms contained all of the terms that were significantly overrepresented in one of the species in this comparison. moreover, out of these terms were found at least times in at least one species, and this subgroup contained out of the significantly overrepresented terms. similarly, in the panther family analysis, out of the distinct families that were found in the bess occurred at least ten times in at least one species, out of which families were found at least times. the families that were found at least ten times contained every one of the families that were significantly overrepresented in one of the species; ten of these were counted more than times in at least one species. while only the tomato hba versus potato pot comparison is shown here, the other comparisons show a similar pattern, indicating that many of the highly abundant go terms and panther families are significantly overrepresented in either tomato or potato. the majority of these overrepresented terms and families are most abundant in potato, and represent biologically important functions and processes. in tomato, a smaller number of terms and families is overrepresented; these are primarily connected to structural genomic features such as retrotransposons. the overrepresentation of transposon-related go terms and panther families in tomato was consistent with the results of the repeat analysis, confirming the observations that tomato is richer in retrotransposons than potato. however, in the panther analysis, reverse transcriptases were significantly overrepresented in potato. at first glance, this does not correspond well with the overrepresentation of rna-directed dna polymerase activity and rna-dependent dna replication in tomato. however, in both tomato and potato, the large majority of the reverse transcriptases originated from non-ltr retroelements, which in fact is consistent with the higher frequency of non-ltr retrotransposons in potato found in the repeatmasker analysis. the cytochrome p mono-oxygenases represent a large gene superfamily in plants that are commonly associated with the biosynthesis of secondary metabolites. in a. thaliana, at least p genes have been found, representing approximately one percent of the gene complement of this species. in o. sativa, this family is even larger, with p genes identified so far. not all the p450s in these genomes represent true protein-coding sequences; in a. thaliana, 90% of the genes are truly protein coding, compared to 72% in o. sativa. in total, distinct families of p genes were identified in a. thaliana and o. sativa, several of which were found to be overrepresented in either of these species. moreover, some families were present in one, but completely absent in the other species. in the hindiii and ecori libraries, and bess that could be associated with the cytochrome p panther family were found in tomato and potato, respectively. since these bac end sequences represent approximately 14% and 10% of their respective genomes, these data suggest an enormous expansion of p genes in the solanaceae. this could be the result of an expansion of specific p families, but also of the evolution of species- or family-specific p450s. for example, the allene oxide synthase has currently only been found in solanaceous species, including tomato and petunia inflata. the overrepresentation of p450s in potato compared to tomato may be another result of species-specific p families, but may also indicate expansion of families that are shared between these species. comparative genome mapping in this study, paired bac ends have been exploited to detect regions of microsynteny between the solanaceous species tomato and potato, and the model plant organisms a. thaliana and p. trichocarpa. using similar approaches, microsynteny has been observed between a. thaliana and b. rapa; c. papaya and p. trichocarpa; and m. acuminata and o. sativa. a higher number of tomato and potato bacs display microsynteny to p. trichocarpa, than to a. thaliana. the reduced level of microsynteny between tomato/potato and a. thaliana is not likely a difference in evolutionary distances between these species. both a. thaliana and p. trichocarpa are part of the rosids clade, whereas tomato and potato belong to the asterids clade. it may be the result of a recent duplication of the a. thaliana genome, followed by the loss of roughly 70% of the duplicated genes. assuming that this loss occurred randomly, the large majority of possible microsyntenic regions that existed before the duplication will have disappeared due to the major local expansions and contractions that would be associated with such a duplication and subsequent loss. this hypothesis is strengthened by the observation that only approximately 1% of a. thaliana bes pairs could previously be mapped within kb to the p. trichocarpa genome, indicating that the organization of these genomes is indeed vastly different. regions of microsynteny have previously been detected between tomato/potato and a. thaliana. a kb region of tomato chromosome containing five genes was shown to be syntenic with a kb region on a. thaliana chromosome although the order and orientation of the genes suggested two inversion events. in another study, a kb bac sequence matched to four different segments on a. thaliana chromosomes and 5; however, each of the four regions in a. thaliana were shorter than their tomato counterpart. recently, five microsyntenic blocks were detected between a region on potato chromosome harbouring a qtl for resistance to late blight and root cyst nematodes, and a. thaliana chromosomes and. these syntenic blocks spanned between three and seven orfs, and were interrupted by non-syntenic blocks. in each of these examples, the microsynteny between tomato/potato and a. thaliana involves shorter regions on the a. thaliana genome than the average tomato and potato bac sequence length. furthermore, regions of synteny are often detected between coding sequences, whereas the fraction of coding sequences in the tomato and potato bess is relatively low, providing a good explanation for the reduced amount of microsynteny between these species observed here. synteny between potato and a. thaliana has also been identified on a genome-wide level using a comparative mapping approach. this revealed putative syntenic blocks between potato and a. thaliana that cover 41% of the potato genetic map, and 50% of the a. thaliana physical map. these syntenic blocks were unevenly distributed over the potato genetic map, and redundant in respect to the number of areas on the a. thaliana genome that displayed synteny to most areas on the potato map. the regions of microsynteny between tomato/potato and a. thaliana that were found with the bes-based approach described in this study could not be used to confirm or renounce any putative higher-order syntenic regions, due to the relatively short distances between the bac ends. six paired tomato bac end matches cluster in the – mb interval of p. trichocarpa chromosome furthermore, seven pairs of potato bess map to the partially overlapping interval between – mb, indicating the presence of either a number of distinct microsyntenic regions, or possibly a single region of macrosynteny, between the tomato/potato and p. trichocarpa genomes. these findings provide an interesting starting point for a detailed comparison between these species in this region, once more tomato and potato genomic sequences become available.  the large scale analysis of tomato and potato bess presented in this paper revealed many interesting structural and functional differences between the genomes of these closely related species. we have shown that the tomato genome is not only more repetitive than the potato genome, but that these genomes also differ in their repeat composition, most importantly in the distribution of gypsy and copia retrotransposons. in contrast to other studied plant genomes, we have shown that the tomato and potato genomes contain a large number of ssrs with a motif length of five, which may be a unique feature of solanaceous genomes. comparative analysis of the putative protein coding regions in these bess revealed an enrichment of these regions in the potato genome. moreover, several protein families were found to be overrepresented in potato compared to tomato, such as cytochrome p mono-oxygenases and serine-threonine protein kinases. the p superfamily appears to have expanded dramatically in both species compared to a. thaliana, suggesting an expanded network of secondary metabolic pathways in the solanaceae. both tomato and potato appear to have low microsynteny with a. thaliana, which is likely a result of this species' relatively recent genomic rearrangement. a higher degree of synteny was observed with p. trichocarpa. difference in evolutionary distances is not likely to be the reason for this increased microsynteny, since both a. thaliana and p. trichocarpa are part of the rosids clade, whereas tomato and potato belong to the asterids clade. taken together, these findings present a first glimpse into the evolution of solanaceous genomes, both within the family and relative to other plant species. when the complete genomic sequences of these species become available, whole-genome comparisons and protein- or repeat-family specific studies may shed more light on the intriguing observations made in this paper. 
this study presents a first genome-wide analysis of these two species, based on two large collections of bac end sequences representing approximately 19% of the tomato genome and 10% of the potato genome. on the other hand, simple sequence repeats are more abundant in potato than in tomato. when the complete genome sequences of these species become available, whole-genome comparisons and protein- or repeat-family specific studies may shed more light on the observations made here. the tomato genome has a higher repeat content than the potato genome, primarily due to a higher number of retrotransposon insertions in the tomato genome. the findings in this paper present a first glimpse into the evolution of solanaceous genomes, both within the family and relative to other plant species. moreover, the p superfamily appears to have expanded spectacularly in both species compared to arabidopsis thaliana, suggesting an expanded network of secondary metabolic pathways in the solanaceae. major gene families such as cytochrome p mono-oxygenases and serine-threonine protein kinases are significantly overrepresented in potato, compared to tomato. both tomato and potato appear to have a low level of microsynteny with a. thaliana.
 tomato and potato are two economically important crop species, the genomes of which are currently being sequenced. this study presents a first genome-wide analysis of these two species, based on two large collections of bac end sequences representing approximately 19% of the tomato genome and 10% of the potato genome.  the tomato genome has a higher repeat content than the potato genome, primarily due to a higher number of retrotransposon insertions in the tomato genome. on the other hand, simple sequence repeats are more abundant in potato than in tomato. the two genomes also differ in the frequency distribution of ssr motifs. based on est and protein alignments, potato appears to contain up to more putative coding regions than tomato. major gene families such as cytochrome p mono-oxygenases and serine-threonine protein kinases are significantly overrepresented in potato, compared to tomato. moreover, the p superfamily appears to have expanded spectacularly in both species compared to arabidopsis thaliana, suggesting an expanded network of secondary metabolic pathways in the solanaceae. both tomato and potato appear to have a low level of microsynteny with a. thaliana. a higher degree of synteny was observed with populus trichocarpa, specifically in the region between and mb on p. trichocarpa chromosome  the findings in this paper present a first glimpse into the evolution of solanaceous genomes, both within the family and relative to other plant species. when the complete genome sequences of these species become available, whole-genome comparisons and protein- or repeat-family specific studies may shed more light on the observations made here. 