In silico mining and characterization of simple sequence repeats from gilthead sea bream (Sparus aurata) expressed sequence tags (EST-SSRs); PCR amplification, polymorphism evaluation and multiplexing and cross-species assays

We screened for simple sequence repeats (SSRs) found in ESTs derived from an EST-database development project ('Marine Genomics Europe' Network of Excellence). Different motifs of di-, tri-, tetra-, penta- and hexanucleotide SSRs were evaluated for variation in length and position in the expressed sequences, relative abundance and distribution in gilthead sea bream (Sparus aurata). We found 899 ESTs that harbor 997 SSRs (4.94%). On average, one SSR was found per 2.95 kb of EST sequence and the dinucleotide SSRs are the most abundant accounting for 47.6% of the total number. EST-SSRs were used as template for primer design. 664 primer pairs could be successfully identified and a subset of 206 pairs of primers was synthesized, PCR-tested and visualized on ethidium bromide stained agarose gels. The main objective was to further assess the potential of EST-SSRs as informative markers and investigate their cross-species amplification in sixteen teleost fish species: seven sparid species and nine other species from different families. Approximately 78% of the primer pairs gave PCR products of expected size in gilthead sea bream, and as expected, the rate of successful amplification of sea bream EST-SSRs was higher in sparids, lower in other perciforms and even lower in species of the Clupeiform and Gadiform orders. We finally determined the polymorphism and the heterozygosity of 63 markers in a wild gilthead sea bream population; fifty-eight loci were found to be polymorphic with the expected heterozygosity and the number of alleles ranging from 0.089 to 0.946 and from 2 to 27, respectively. These tools and markers are expected to enhance the available genetic linkage map in gilthead sea bream, to assist comparative mapping and genome analyses for this species and further with other model fish species and finally to help advance genetic analysis for cultivated and wild populations and accelerate breeding programs.     

Exploiting EST databases for the development and characterization of EST-SSRs in the Pacific oyster (Crassostrea gigas)

A total of 147 microsatellite-containing expressed sequence tags (ESTs) (3.63%) were detected from 4053 ESTs of the Pacific oyster (Crassostrea gigas) in GenBank. The average density of simple sequence repeats (SSRs) was 1 per 8.25 kb of EST after redundancy elimination. Dinucleotide repeat motifs appeared to be the most abundant type. Sixteen new polymorphic EST-SSRs were developed. The number of alleles per locus varied from 3 to 12, with an average of 5.9 alleles per locus. Marker transferability was tested on 2 other Crassostrea species, and 14 loci gave successful amplifications in both species. Twenty EST-SSRs were tested on 3 families of C. gigas for examination of inheritance mode of EST-SSRs. Thirty-five tests of segregation ratios revealed 5 significant departures from expected Mendelian ratios, 4 of which confirmed Mendelian expectations when accounting for the presence of null alleles. Null alleles were detected at 3 loci (15.0%) of the 20 loci, and the frequency of null alleles at EST-SSRs was lower than that in genomic SSRs in C. gigas. The results obtained in this study suggest that C. gigas EST-SSRs will complement the currently available genomic SSR markers and may be useful for comparative mapping, marker-assisted selection, and evolutionary studies.     

An in silico mining for simple sequence repeats from expressed sequence tags of zebrafish, medaka, Fundulus, and Xiphophorus

Teleost fish genome projects involving model species are resulting in a rapid accumulation of genomic and expressed DNA sequences in public databases. The expressed sequence tags (ESTs) collected in the databases can be mined for the analysis of both structural and functional genomics. In this study, we in silico analyzed 49,430 unigenes representing a total of 692,654 ESTs from four model fish for their potential use in developing simple sequence repeats (SSRs), or microsatellites. After bioinformatical mining, a total of 3,018 EST derived SSRs (EST-SSRs) were identified for 2,335 SSR containing ESTs (SSR-ESTs). The frequency of identified SSR-ESTs ranged from 1.5% for Xiphophorus to 7.3% for zebrafish. The dinucleotide repeat motif is the most abundant SSR, accounting for 47%, 52%, 64%, and 78% for medaka, Fundulus, zebrafish, and Xiphophorus, respectively. Simulation analysis suggests that a majority of these EST-SSRs have sufficient flanking sequences for polymerase chain reaction (PCR) primer design. Comparative DNA sequence analyses of SSR-ESTs identified several cross-species SSRs and sequences that may be used as cross-reference genes in comparative studies. For example, the flanking sequences of one SSR (CTG)n within the pituitary tumor-transforming gene (PTTG) 1 interacting protein (PTTGIP), showed conservation spanning the medaka, Fundulus, human, and mouse genomes. This study provides a large body of information on EST-SSRs that can be useful for the development of polymorphic markers, gene mapping, and comparative genome analysis. Functional analysis of these SSR-ESTs may reveal their role in metabolism and gene evolution of these model species.     

Development of expressed sequence tag derived-simple sequence repeats in the genus Lilium

Although lily is the second largest flower crop in cutting flower commodity, only six simple sequence repeats SSRs have been reported. Thus, we developed expressed sequence tag derived-SSRs (EST-SSRs) for the Lilium genus. Among 2,235 unique ESTs, 754 ESTs contained SSR motifs, among which 165 ESTs were amenable to primer design. Among these 165 EST-SSRs, 131 EST-SSRs showed amplification in at least one Lilium species, and 76 EST-SSRs showed amplification in at least nine species. Of the 76 EST-SSRs, 47 showed amplification in all Lilium species analyzed. Using 10 breeding lines, we selected 21 EST-SSRs that had the highest number of alleles and polymorphism information content. The polymorphism information content values of these selected EST-SSRs ranged from 0.49 to 0.94 with an average of 0.76, which are higher than other plant species. The phylogenetic dendrogram derived from the amplification profiles of the 21 high polymorphic EST-SSRs was congruent with the genetic background of the 84 selected lily accessions and hybrids, which are available in commerce. Thus, the developed EST-SSRs will be very useful in germplasm management, genetic diversity analysis, cultivar finger printing, and molecular breeding in the lily.     

Development of EST-SSR Markers by Data Mining in Three Species of Shrimp: Litopenaeus vannamei, Litopenaeus stylirostris, and Trachypenaeus birdy

We report on the data mining of publicly available Litopenaeus vannamei expressed sequence tags (ESTs) to generate simple sequence repeat (SSRs) markers and on their transferability between related Penaeid shrimp species. Repeat motifs were found in 3.8% of the evaluated ESTs at a frequency of one repeat every 7.8 kb of sequence data. A total of 206 primer pairs were designed, and 112 loci were amplified with the highest success in L. vannamei. A high percentage (69%) of EST-SSRs were transferable within the genus Litopenaeus. More than half of the amplified products were polymorphic in a small testing panel of L. vannamei. Evaluation of those primers in a larger testing panel showed that 72% of the markers fit Hardy-Weinberg equilibrium, which shows their utility for population genetic analysis. Additionally, a set of 26 of the EST-SSRs were evaluated for Mendelian segregation. A high percentage of monomorphic markers (46%) proved to be polymorphic by singles-stranded conformational polymorphism analysis. Because of the high number of ESTs available in public databases, a data mining approach similar to the one outlined here might yield high numbers of SSR markers in many animal taxa.     

小麦EST-SSRs的通用性研究

小麦EST—SSRs是一种从小麦ESTs序列中开发的新型SSR标记。根据小麦ESTs设计的597对EST—SSRs引物在小麦、玉米、水稻和大豆中的有效扩增比率分别为80％、65．8％、62．1％和58．1％。其中,255对EST—SSRs引物(42．7％)在四种作物中都能获得预期产物。序列相似性比较发现,255个通用EST—SSRs标记所对应的基因80％功能已知,其中30％与蛋白质定位相关,其余70％参与细胞代谢、转录调控、细胞发育、防御体系及其他所涉及的生理生化过程。由于不同物种间基因保守性,本研究发掘的EST—SSRs引物对于小麦功能基因组和比较基因组研究都有重要意义。     

Robust simple sequence repeat markers for spruce (Picea spp.) from expressed sequence tags

Traditionally, simple sequence repeat (SSR) markers have been developed from libraries of genomic DNA. However, the large, repetitive nature of conifer genomes makes development of robust, single-copy SSR markers from genomic DNA difficult. Expressed sequence tags (ESTs), or sequences of messenger RNA, offer the opportunity to exploit single, low-copy, conserved sequence motifs for SSR development. From a 20,275-unigene spruce EST set, we identified 44 candidate EST-SSR markers. Of these, 25 amplified and were polymorphic in white, Sitka, and black spruce; 20 amplified in all 23 spruce species tested; the remaining five amplified in all except one species. In addition, 101 previously described spruce SSRs (mostly developed from genomic DNA), were tested. Of these, 17 amplified across white, Sitka, and black spruce. The 25 EST-SSRs had approximately 9% less heterozygosity than the 17 genomic-derived SSRs (mean H=0.65 vs 0.72), but appeared to have less null alleles, as evidenced by much lower apparent inbreeding (mean F=0.046 vs 0.126). These robust SSRs are of particular use in comparative studies, and as the EST-SSRs are within the expressed portion of the genome, they are more likely to be associated with a particular gene of interest, improving their utility for quantitative trait loci mapping and allowing detection of selective sweeps at specific genes.     

EST databases as a source for molecular markers: lessons from Helianthus

Expressed sequence tag (EST) databases represent a potentially valuable resource for the development of molecular markers for use in evolutionary studies. Because EST-derived markers come from transcribed regions of the genome, they are likely to be conserved across a broader taxonomic range than are other sorts of markers. This paper describes a case study in which the publicly available cultivated sunflower (Helianthus annuus) EST database was used to develop simple sequence repeat (SSR) markers for use in the genetic analysis of a rare sunflower species, Helianthus verticillatus, as well as the more widespread Helianthus angustifolius. EST-derived SSRs were found to be more than 3 times as transferable across species as compared with anonymous SSRs (73% vs. 21%, respectively). Moreover, EST-SSRs whose primers were located within protein-coding sequence were more readily transferable than those derived from untranslated regions, and the former loci were no less variable than the latter. The utility of existing EST databases as a means for facilitating population genetic analyses in plants was further explored by cross-referencing publicly available EST resources against available lists of rare or invasive flowering plant taxa. This survey revealed that more than one-third of all plant-derived EST collections of sufficient size could conceivably serve as a source of EST-SSRs for the analysis of rare, endangered, or invasive plant species worldwide.     

中国板栗EST-SSR信息分析及其通用性

从已公布的壳斗科(Fagaceae)基因组数据库中获得48501条中国板栗(Castanea mollissima)EST序列,其中8479条EST含有12116个长度大于10 bp的SSR,EST-SSR的出现频率为24.98%。二核苷酸重复和三核苷酸重复是中国板栗EST-SSR最主要的重复类型,分别占总数的38.05%和42.20%。对中国板栗EST-SSR标记在茅栗(C.seguinii)和锥栗(C.henryi)上的通用性检测表明,677对中国板栗EST-SSR引物皆能扩增,证实这些引物在栗属中国特有种间具有很高的通用性。115个位点的多态性分析表明,中国板栗、茅栗和锥栗分别有106、108和101个位点表现出多态性。在575个等位基因中,有260个(45.22%)是3物种的共有等位基因,各物种都有各自特有的等位基因。这项研究为中国栗属植物EST-SSR标记的建立及应用奠定了基础。     

In silico mining for simple sequence repeat loci in a pineapple expressed sequence tag database and cross-species amplification of EST-SSR markers across Bromeliaceae

A collection of 5,659 expressed sequence tags (ESTs) from pineapple [Ananas comosus (L.) Merr.] was screened for simple sequence repeats (EST-SSRs) with motif lengths between 1 and 6 bp. Lower thresholds of 15, 7 and 5 repeat units were used to define microsatellites of the mono-, di-, and tri- to hexanucleotide repeat type, respectively. Based on these criteria, 696 SSRs were identified among 3,389 EST unigenes, together representing 2,840 kb. This corresponds to an average density of one SSR every 4.1 kb of non-redundant EST sequences. Dinucleotide repeats were most abundant (38.4% of all SSRs) followed by trinucleotide repeats (38.1%). Flanking primer pairs were designed for 537 EST-SSR loci, and 49 of these were screened for their functionality in 12 accessions of A. comosus, 14 accessions of 5 additional Ananas species and 1 species of Pseudananas. Distinct PCR products of the expected size range were obtained with 36 primer pairs. Eighteen loci analyzed in more detail were all polymorphic in pineapple, and primer pairs flanking these loci also generated PCR products from a wide range of genera and species from six subfamilies of the Bromeliaceae. The potential to reveal polymorphism in a heterologous target species was demonstrated in Deuterocohnia brevifolia (subfamily Pitcairnioideae).     

First Microsatellite Markers Developed from Cupuassu ESTs: Application in Diversity Analysis and Cross-Species Transferability to Cacao

The cupuassu tree (Theobroma grandiflorum) (Willd. ex Spreng.) Schum. is a fruitful species from the Amazon with great economical potential, due to the multiple uses of its fruit´s pulp and seeds in the food and cosmetic industries, including the production of cupulate, an alternative to chocolate. In order to support the cupuassu breeding program and to select plants presenting both pulp/seed quality and fungal disease resistance, SSRs from Next Generation Sequencing ESTs were obtained and used in diversity analysis. From 8,330 ESTs, 1,517 contained one or more SSRs (1,899 SSRs identified). The most abundant motifs identified in the EST-SSRs were hepta- and trinucleotides, and they were found with a minimum and maximum of 2 and 19 repeats, respectively. From the 1,517 ESTs containing SSRs, 70 ESTs were selected based on their functional annotation, focusing on pulp and seed quality, as well as resistance to pathogens. The 70 ESTs selected contained 77 SSRs, and among which, 11 were polymorphic in cupuassu genotypes. These EST-SSRs were able to discriminate the cupuassu genotype in relation to resistance/susceptibility to witches’ broom disease, as well as to pulp quality (SST/ATT values). Finally, we showed that these markers were transferable to cacao genotypes, and that genome availability might be used as a predictive tool for polymorphism detection and primer design useful for both Theobroma species. To our knowledge, this is the first report involving EST-SSRs from cupuassu and is also a pioneer in the analysis of marker transferability from cupuassu to cacao. Moreover, these markers might contribute to develop or saturate the cupuassu and cacao genetic maps, respectively.     

Mining for SSRs and FDMs from expressed sequence tags of Camellia sinensis

Simple Sequence Repeats (SSRs) developed from Expressed Sequence Tags (ESTs), known as EST-SSRs are most widely used and potentially valuable source of gene based markers for their high levels of crosstaxon portability, rapid and less expensive development. The EST sequence information in the publicly available databases is increasing in a faster rate. The emerging computational approach provides a better alternative process of development of SSR markers from the ESTs than the conventional methods. In the present study, 12,851 EST sequences of Camellia sinensis, downloaded from National Center for Biotechnology Information (NCBI) were mined for the development of Microsatellites. 6148 (4779 singletons and 1369 contigs) non redundant EST sequences were found after preprocessing and assembly of these sequences using various computational tools. Out of total 3822.68 kb sequence examined, 1636 (26.61%) EST sequences containing 2371 SSRs were detected with a density of 1 SSR/1.61 kb leading to development of 245 primer pairs. These mined EST-SSR markers will help further in the study of variability, mapping, evolutionary relationship in Camellia sinensis. In addition, these developed SSRs can also be applied for various studies across species.     

Development, characterization, and comparative analysis of polymorphism at common bean SSR loci isolated from genic and genomic sources.

Microsatellites or SSRs (single sequence repeats) have been used to construct and integrate genetic maps in crop species, including Phaseolus vulgaris. In the present study, 3 cDNA libraries generated by the Bean EST project (http://lgm.esalq.usp.br/BEST/), comprising a unigene collection of 3126 sequences and a genomic microsatellite-enriched library, were analyzed for the presence of SSRs. A total of 219 expressed sequence tags (ESTs) were found to carry 240 SSRs (named EST-SSR), whereas 714 genomic sequences contained 471 SSRs (named genomic-SSR). A subset of 80 SSRs, 40 EST-SSRs, and 40 genomic-SSRs were evaluated for molecular polymorphism in 23 genotypes of cultivated beans from the Mesoamerican and Andean genetic pools, including Brazilian cultivars and 2 related species. Of the common bean genotypes, 31 EST-SSR loci were polymorphic, yielding 2-12 alleles as compared with 26 polymorphic genomic-SSRs, accounting for 2-7 alleles. Cluster analysis from data using both genic and genomic-SSR revealed a clear separation between Andean and Mesoamerican beans. The usefulness of these loci for distinguishing bean genotypes and genetic mapping is discussed.     

Comparative population genetic analysis of bocaccio rockfish Sebastes paucispinis using anonymous and gene-associated simple sequence repeat loci

Comparative population genetic analyses of traditional and emergent molecular markers aid in determining appropriate use of new technologies. The bocaccio rockfish Sebastes paucispinis is a high gene-flow marine species off the west coast of North America that experienced strong population decline over the past 3 decades. We used 18 anonymous and 13 gene-associated simple sequence repeat (SSR) loci (expressed sequence tag [EST]-SSRs) to characterize range-wide population structure with temporal replicates. No F(ST)-outliers were detected using the LOSITAN program, suggesting that neither balancing nor divergent selection affected the loci surveyed. Consistent hierarchical structuring of populations by geography or year class was not detected regardless of marker class. The EST-SSRs were less variable than the anonymous SSRs, but no correlation between F(ST) and variation or marker class was observed. General linear model analysis showed that low EST-SSR variation was attributable to low mean repeat number. Comparative genomic analysis with Gasterosteus aculeatus, Takifugu rubripes, and Oryzias latipes showed consistently lower repeat number in EST-SSRs than SSR loci that were not in ESTs. Purifying selection likely imposed functional constraints on EST-SSRs resulting in low repeat numbers that affected diversity estimates but did not affect the observed pattern of population structure.     

Development of 16 polymorphic simple sequence repeat markers for Lycoris longituba from expressed sequence tags

Lycoris longituba is a tulip-like ornamental plant in China. We report on the data mining of L. longituba expressed sequence tags (ESTs) to generate simple sequence repeat (SSRs) markers. Eighteen EST-SSRs were isolated and validated for 32 individuals. These markers will be valuable for studying the genetic structure and diversity of populations for L. longituba.     

Mining and survey of simple sequence repeats in expressed sequence tags of dicotyledonous species.

Simple sequence repeat (SSR) markers are widely used in many plant and animal genomes due to their abundance, hypervariability, and suitability for high-throughput analysis. Development of SSR markers using molecular methods is time consuming, laborious, and expensive. Use of computational approaches to mine ever-increasing sequences such as expressed sequence tags (ESTs) in public databases permits rapid and economical discovery of SSRs. Most of such efforts to date focused on mining SSRs from monocotyledonous ESTs. In this study, we have computationally mined and examined the abundance of SSRs in more than 1.54 million ESTs belonging to 55 dicotyledonous species. The frequency of ESTs containing SSRs among species ranged from 2.65% to 16.82%. Dinucleotide repeats were found to be the most abundant followed by tri- or mono-nucleotide repeats. The motifs A/T, AG/GA/CT/TC, and AAG/AGA/GAA/CTT/TTC/TCT were the predominant mono-, di-, and tri-nucleotide SSRs, respectively. Most of the mononucleotide SSRs contained 15-25 repeats, whereas the majority of the di- and tri-nucleotide SSRs contained 5-10 repeats. The comprehensive SSR survey data presented here demonstrates the potential of in silico mining of ESTs for rapid development of SSR markers for genetic analysis and applications in dicotyledonous crops.     

Development of 14 EST-SSRs for Betula maximowicziana and their applicability to related species

Simple sequence repeats (SSRs) markers were developed for Betula maximowicziana using 2698 expressed sequence tags (ESTs) from the NCBI database. Out of 112 designed primer pairs, 54 showed clear PCR amplification and 14 of these revealed polymorphism in eight individuals sampled across the species’ range. The number of alleles detected and the expected heterozygosity ranged from 1 to 3 and 0.000 to 0.570, respectively, when these 14 loci were examined in 49 individuals from a single population. In the cross species transferability test, eight of the 14 loci were also polymorphic in all four of the diploid, tetraploid and hexaploid Betula species examined. These results showed high transferability of the developed EST-SSRs and that these markers are likely to be useful in studies of the population genetics of species in the genus Betula.     

Study of simple sequence repeat (SSR) markers from wheat expressed sequence tags (ESTs)

The increasing availability of expressed sequence tags (ESTs) in wheat (Triticum aestivum) and related cereals provides a valuable resource of non-anonymous DNA molecular markers. We examined 170,746 wheat ESTs from the public (International Triticeae EST Cooperative) and Génoplante databases, previously clustered in contigs, for the presence of di- to hexanucleotide simple sequence repeats (SSRs). Analysis of 46,510 contigs identified 3,530 SSRs, which represented 7.5% of the total number of contigs. Only 74% of the sequences allowed primer pairs to be designed, 70% led to an amplification product, mainly of a high quality (68%), and 53% exhibited polymorphism for at least one cultivar among the eight tested. Even though dinucleotide SSRs were less represented than trinucleotide SSRs (15.5% versus 66.5%, respectively), the former showed a much higher polymorphism level (83% versus 46%). The effect of the number and type of repeats is also discussed. The development of new EST-SSRs markers will have important implications for the genetic analysis and exploitation of the genetic resources of wheat and related species and will provide a more direct estimate of functional diversity.     

Genic Microsatellite Markers in Brassica rapa: Development, Characterization, Mapping, and Their Utility in Other Cultivated and Wild Brassica Relatives

Genic microsatellite markers, also known as functional markers, are preferred over anonymous markers as they reveal the variation in transcribed genes among individuals. In this study, we developed a total of 707 expressed sequence tag-derived simple sequence repeat markers (EST-SSRs) and used for development of a high-density integrated map using four individual mapping populations of B. rapa. This map contains a total of 1426 markers, consisting of 306 EST-SSRs, 153 intron polymorphic markers, 395 bacterial artificial chromosome-derived SSRs (BAC-SSRs), and 572 public SSRs and other markers covering a total distance of 1245.9 cM of the B. rapa genome. Analysis of allelic diversity in 24 B. rapa germplasm using 234 mapped EST-SSR markers showed amplification of 2 alleles by majority of EST-SSRs, although amplification of alleles ranging from 2 to 8 was found. Transferability analysis of 167 EST-SSRs in 35 species belonging to cultivated and wild brassica relatives showed 42.51% (Sysimprium leteum) to 100% (B. carinata, B. juncea, and B. napus) amplification. Our newly developed EST-SSRs and high-density linkage map based on highly transferable genic markers would facilitate the molecular mapping of quantitative trait loci and the positional cloning of specific genes, in addition to marker-assisted selection and comparative genomic studies of B. rapa with other related species.