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).     

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.     

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.     

Exploitation of a turbot (Scophthalmus maximus L.) immune-related expressed sequence tag (EST) database for microsatellite screening and validation

In this study, we identified and characterized 160 microsatellite loci from an expressed sequence tag (EST) database generated from immune-related organs of turbot (Scophthalmus maximus). A final set of 83 new polymorphic microsatellites were validated after the analysis of 40 individuals of Atlantic origin including both wild and farmed individuals. The allele number and the expected heterozygosity ranged from 2 to 18 and from 0.021 to 0.951, respectively. Evidences of null alleles at moderate-high frequencies were detected at six loci using population data. None of the analysed loci showed deviations from Mendelian segregation after the analysis of five full-sib families including approximately 92 individuals/family. The markers are used to consolidate the turbot genetic map, and because they are mostly EST-derived, they will be very useful for comparative genomic studies within flatfishes and with model fish species. Using an in silico approach, we detected significant homologies of microsatellite sequences with the EST databases of the flatfish species with highest genomic resources (Senegalese sole, Atlantic halibut, bastard halibut) in 31% of these turbot markers. The conservation of these microsatellites within Pleuronectiformes will pave the way for anchoring genetic maps of different species and identifying genomic regions related to productive traits.     

Identification and characterization of SSRs from soybean (Glycine max) ESTs

Microsatellites, or simple sequence repeats (SSRs), are very useful molecular markers for a number of plant species. We used a new publicly available module (TROLL) to extract microsatellites from the public database of soybean expressed sequence tag (EST) sequences. A total of 12,833 sequences containing di- to penta-type SSRs were identified from 200,516 non-redundant soybean ESTs. On average, one SSR was found per 7.25?kb of EST sequences, with the tri-nucleotide motifs being the most abundant. Primer sequences flanking the SSR motifs were successfully designed for 9,638 soybean ESTs using the software primer3.0 and only 59 pairs of them were found in earlier studies. We synthesized 124 pairs of the primers to determine the polymorphism and heterozygosity among eight genotypes of soybean cultivars, which represented a wide range of the cultivated soybean cultivars. PCR amplification products with anticipated SSRs were obtained with 81 pairs of primers; 36 PCR products appeared to be homozygous and the remaining 45 PCR products appeared to be heterozygous and displayed polymorphism among the eight cultivars. We further analysed the EST sequences containing 45 polymorphic EST-SSR markers using the programs BLASTN and BLASTX. Sequence alignment showed that 29 ESTs have homologous sequences and 15 ESTs could be classified into a Uni-gene cluster with comparatively convincing protein products. Among these 15 ESTs belonging to a Uni-gene cluster, 9 SSRs were located in 3'-UTR, 4 SSRs were located in the intron region and 2 SSRs were located in the CDS region. None of these SSRs was located in the 5'-UTR. These novel SSRs identified in the ESTs of soybean provide useful information for gene mapping and cloning in future studies.     

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.     

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.     

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.     

Development and characterization of 60 novel EST-SSR markers in half-smooth tongue sole Cynoglossus semilaevis

Sixty novel simple sequence repeat (SSR) markers were developed from expressed sequence tags (EST) of half-smooth tongue sole Cynoglossus semilaevis exploited in the laboratory. The number of alleles, observed and expected heterozygosity per locus ranged from two to 16, from 0·0833 to 1·0000 and from 0·0816 to 0·913, respectively. Of these SSRs, 20 had significant homology to known genes by BLASTx (basic local alignment search tool x) search. For cross-species amplification, there are 53 positive amplifications in Japanese flounder Paralichthys olivaceus with 12 polymorphic loci and 51 positive amplifications in Senegalese sole Solea senegalensis with 11 polymorphic loci. These new EST-SSR markers will be useful for genetic studies and genome mapping of C. semilaevis and its closely related fishes.     

Utility of EST-derived SSRs as population genetics markers in a beetle

Microsatellite, or simple sequence repeat (SSR), loci can be identified by mining expressed sequence tag (EST) databases, and where these are available, marker development time and expense can be decreased considerably over conventional strategies of probing the entire genome. However, it is unclear whether they provide information on population structure similar to that generated by anonymous genomic SSRs. We performed comparative population genetic analyses between EST-derived SSRs (EST-SSRs) and anonymous SSRs developed from genomic DNA for the same set of populations of the insect Diabrotica virgifera, a beetle in the family Chrysomelidae. Compared with noncoding, nontranscribed regions, EST-SSRs were generally less polymorphic but had reduced occurrence of null alleles and greater cross-species amplification. Neutrality tests suggested the loci were not under positive selection. Across all populations and all loci, the genomic and EST-SSRs performed similarly in estimating genetic diversity, F(IS), F(ST), population assignment and exclusion tests, and detection of distinct populations. These findings, therefore, indicate that the EST-SSRs examined can be used with confidence in future genetic studies of Diabrotica populations and suggest that EST libraries can be added as a valuable source of markers for population genetics studies in insects and other animals.     

Development of SSR markers for the phylogenetic analysis of almond trees from China and the Mediterranean region.

Expressed sequence tag (EST) derived simple sequence repeats (SSRs, microsatellites) were screened and identified from 3863 almond and 10 185 peach EST sequences, and the spectra of SSRs in the non-redundant EST sequences were investigated after sequence assembly. One hundred seventy-eight (12.07%) almond SSRs and 497 (9.97%) peach SSRs were detected. The EST-SSR occurs every 4.97 kb in almond ESTs and 6.57 kb in peach, and SSRs with di- and trinucleotide repeat motifs are the most abundant in both almond and peach ESTs. Twenty one EST-SSRs were thereafter, developed and used together with 7 genomic SSRs, to study the genetic relationship among 36 almond (P. communis Fritsch.) cultivars from China and the Mediterranean area, as well as 8 accessions of other related species from the genus Prunus. Both EST-derived and genomic SSR markers showed high cross-species transferability in the genus. Out of the 112 polymorphic alleles detected in the 36 cultivated almonds, 28 are specific to Chinese cultivars and 25 to the others. The 44 accessions were clustered into 4 groups in the phylogenetic tree and the 36 almond cultivars formed two distinct subgroups, one containing only Chinese cultivars and one of unknown origin and the other only those originating from the Mediterranean area, indicating that Chinese almond cultivars have a distinct evolutionary history from the Mediterranean almond. Our preliminary results indicated that common almond was more closely related to peach (P. persica (L.) Batsch.) than to the four wild species of almond, (P. mongolica Maxim., P. ledebouriana Schleche, P. tangutica Batal., and P. triloba Lindl.). The implications of these SSR markers for evolutionary analysis and molecular mapping of Prunus species are discussed.     

Eucalyptus microsatellites mined in silico: survey and evaluation

Eucalyptus is an important short rotation pulpy woody plant, grown widely in the tropics. Recently, many genomic programmes are underway leading to the accumulation of voluminous genomic and expressed sequence tag sequences in public databases. These sequences can be utilized for analysis of simple sequence repeats (SSRs) and single nucleotide polymorphism (SNPs) available in the transcribed genes. In this study, in silico analysis of 15,285 sequences representing partial and full-length mRNA from Eucalyptus species for their use in developing SSRs or microsatellites were carried out. A total of 875 EST-SSRs were identified from 772 SSR containing ESTs. Motif size of 6 for dinucleotide and 5 for trinucleotide, tetranucleotide, and pentanucleotides were considered in locating the microsatellites. The average frequency of identified SSRs was 12.9%. The dinucleotide repeats were the most abundant among the dinucleotide, trinucleotide and tetranucleotide motifs and accounted for 50.9% of the Eucalyptus genome. Primer designing analysis showed that 571 sequences with SSRs had sufficient flanking regions for polymerase chain reaction (PCR) primer synthesis. Evaluation of the usefulness of the SSRs showed that EST-derived SSRs can generate polymorphic markers as all the primers showed allelic diversity among the 16 provenances of E. tereticornis.     

Mining and charactering microsatellites in Cucumis melo expressed sequence tags from sequence database

Simple sequence repeats (SSRs) derived from expressed sequence tags (ESTs) are valuable markers because they represent transcribed regions and often have putative functions. We mined and characterized microsatellites in melon ESTs. Three hundred and eighty‐three SSR loci were identified in 309 of 3188 unigenes assembled by 5747 EST and mRNA sequences in GenBank with occurring frequency of 1/4.7 kb. Twenty‐two polymorphic EST‐SSR markers were developed with the mean allele number of 2.9 per locus and mean expected heterozygosity of 0.442. Amplification products were also detected by 15 pairs of primer in Cucumis sativus. Those informative EST‐SSR markers can be used in melon genetic improvement projects.     

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.     

小麦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引物对于小麦功能基因组和比较基因组研究都有重要意义。     

Application of functional genomic information to develop efficient EST-SSRs for the chicken (Gallus gallus)

Many years of domestication and breeding have given rise to the wide range of chicken breeds that exist today; however, an increasing number of local chicken breeds are under threat of extinction. A comprehensive characterization of chicken markers (especially type I markers) is needed to monitor and conserve genetic diversity in this species. The explosion of genomics and functional genomics information in recent years has opened new possibilities for the generation of molecular markers. We analyzed a large number of expressed sequence tags (ESTs) to test the possibility of using EST-derived microsatellite markers for investigating the Gallus gallus genome. Chromosomal locations for the majority of these SSRs were predicted. Of the 31,576 unigenes assembled from the 544,150 redundant EST sequences, 1757 SSR markers were discovered on 1544 ESTs, using the SSRLocator software, with an average density of 28.7 kb per SSR. The dimer motifs were the most frequent (46.38%), followed by trimeric (38.58%), tetrameric (10.19%), pentameric (4.5%), and hexameric (<1%) markers. Different from the case for cattle and sheep, AT/TA was the most abundant dimeric repeat, accounting for 41.71% of all dimeric repeats in the chicken ESTs. The EST-SSR distribution was not uniform among the chromosomes; the majority of the EST-SSRs were located on chromosomes GGA2 and GGA10. We found that most of the EST-SSRs are involved in positive regulation of cellular and metabolic processes. This is the first time that EST sequences have been mined to find chicken microsatellites. On average, 3.8% of the G. gallus UniGene sequences could be exploited for development of EST-SSRs, indicating a good source for molecular markers as well as for functional genome analysis.     

Bin mapping of genomic and EST-derived SSRs in melon (Cucumis melo L.)

We report the development of 158 primer pairs flanking SSR motifs in genomic (gSSR) and EST (EST-SSR) melon sequences, all yielding polymorphic bands in melon germplasm, except one that was polymorphic only in Cucurbita species. A similar polymorphism level was found among EST-SSRs and gSSRs, between dimeric and trimeric EST-SSRs, and between EST-SSRs placed in the open reading frame or any of the 5′- or 3′-untranslated regions. Correlation between SSR length and polymorphism was only found for dinucleotide EST-SSRs located within the untranslated regions, but not for trinucleotide EST-SSRs. Transferability of EST-SSRs to Cucurbita species was assayed and 12.7% of the primer pairs amplified at least in one species, although only 5.4% were polymorphic. A set of 14 double haploid lines from the cross between the cultivar “Piel de Sapo” and the accession PI161375 were selected for the bin mapping approach in melon. One hundred and twenty-one SSR markers were newly mapped. The position of 46 SSR loci was also verified by genotyping the complete population. A final bin-map was constructed including 80 RFLPs, 212 SSRs, 3 SNPs and the Nsv locus, distributed in 122 bins with an average bin length of 10.2 cM and a maximum bin length of 33 cM. Map density was 4.2 cM/marker or 5.9 cM/SSR. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Assessment of genetic diversity and population structure of Chinese wild almond, Amygdalus nana, using EST- and genomic SSRs

Amygdalus nana L., commonly known as wild almond, is an endangered wild relative of cultivated almond, which has great potential in almond crop breeding. In this study, we used microsatellite (SSR) loci derived from both expressed sequence tag (EST) and anonymous genomic sequence to explore the genetic diversity and population structure of A. nana in Xinjiang of China. Seven natural populations were collected across the whole distribution of A. nana in China, including populations from both inside (four populations) and outside (three populations) the established protected areas. A total of 22 and 19 alleles were detected from the seven pairs of EST and genomic SSR loci, respectively. Generally, the genomic SSRs showed lower levels of variation than EST-SSRs, which may partially due to the higher cross-species transferability in EST-SSRs than in genomic SSRs. The population-level genetic diversity (A = 1.84, P = 50.00%, H_o = 0.3491, H_E = 0.2271) was lower than cultivated almond and several wild fruit species with similar breeding system. Most of the genetic variation (82.16%) was partitioned within populations. In particular, the population collected from Tacheng County (outside the protected areas) had the highest levels of genetic diversity and had significantly different genetic constitution from other populations.