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.     

EST-SSRs characterization and in-silico alignments with linkage map SSR loci in grape (Vitis L.) genome

11,581 grape (Vitis L.) EST-SSRs were produced and characterized from a total of 381,609 grape ESTs. Among the EST-SSRs, the tri repeat (5,560, 45.4%) represented the most abundant class of microsatellites in grape EST. Most of grape EST-SSR motifs fall within 18-24 bps in length. The EST-SSRs tri-repeats occurred a higher percentage in 5??-end (59.3%) than in 3??-end (48.3%). And EST-SSR tri-repeats had abundant codon repeats for putative amino acid runs as Proline, Arginine in grape ESTs. To better utilizing these markers, 142 of newly developed and validated EST SSR loci as well as 223 linkage map SSR loci were in silico aligned and mapped in grape genome. The orders of these SSR loci in the chromosomal physical locations and in the linkage groups were compared, and about twenty linkage map loci positions were switched or rearranged in grape genome. The EST-SSR markers extended the linkage map in grape genome. The method of in silico mapping reported in this study provided an initial collection for grape mapping resources. This approach offers great opportunities to understand the genetic variations in nucleotide sequences differences in physical map, and genetic recombination in linkage maps, as well as benefits for markers enrichment in a specific grape genome region for fine mapping or QTL mapping.     

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.     

Analysis of expressed sequence tags from grapevine flower and fruit and development of simple sequence repeat markers

A total of 6,230 EST sequences were produced from 7,561 clones in a cDNA library generated from grapevine (Vitis vinifera cv. ‘Summer Black’) flower and fruit tissues in this study. After cluster and assembly analysis of the datasets, 3,582 unigenes (GenBank accession numbers GW836604–GW840185) were established, among which 381 were new grapevine EST sequences. Out of the 381 new ESTs, 289 could be mapped on the 19 grapevine chromosomes. 913 unique ESTs with known or putative functions were assigned to 11 putative cellular roles. 540 potentially workable grapevine EST-SSRs were developed from 3,582 unigenes and about 42.6% of these unigenes were identified as true-to-type SSR loci and could amplify polymorphic bands from 22 individual plants of V. vinifera L, indicating that grapevine EST datasets are a valuable source for the development of functional simple sequence repeat (SSR) markers.     

茶树EST-SSRs分布特征及引物开发

为了在茶树中开发EST-SSRs功能性标记,利用生物信息学方法对NCBI网上公开的3288奈茶树（Camellia subebsus）ESTs序列进行EST-SSRs特征分析。剔除冗余序列,得到非冗余序列2083条。在非冗余序列中发现含不同重复基元SSRs的EST序列有385条,共486个EST-SSRs,平均相隔2．10kb出现1个SSR。在2～6bp的重复基元中,二核苷酸重复基元的SSRs出现频率最高（51．97％）,其次是三核苷酸（19．55％）。对所有的重复基元类型进行统计分析发现,所占比例最高的是AG／CT（47．74％）,其次分别是AT／TA（4．73％）和AAG／CTT（4．73％）。利用Prime5软件,设计了206对EST-SSRs引物,随机选用72对引物进行SSR扩增,发现31对引物可以扩增出条带,其中29对引物具有多态性,多态性比率为93．5％。这些EST-SSRs将有助于茶树基因组学方面的研究。     

Evidence for abundant transcription of non-coding regions in the Saccharomyces cerevisiae genome

Background

The ChickRH6 whole chicken genome radiation hybrid (RH) panel recently produced has already been used to build radiation hybrid maps for several chromosomes, generating comparative maps with the human and mouse genomes and suggesting improvements to the chicken draft sequence assembly. Here we present the construction of a RH map of chicken chromosome 2. Markers from the genetic map were used for alignment to the existing GGA2 (Gallus gallus chromosome 2) linkage group and EST were used to provide valuable comparative mapping information. Finally, all markers from the RH map were localised on the chicken draft sequence assembly to check for eventual discordances.

Results

Eighty eight microsatellite markers, 10 genes and 219 EST were selected from the genetic map or on the basis of available comparative mapping information. Out of these 317 markers, 270 gave reliable amplifications on the radiation hybrid panel and 198 were effectively assigned to GGA2. The final RH map is 2794 cR₆₀₀₀ long and is composed of 86 framework markers distributed in 5 groups. Conservation of synteny was found between GGA2 and eight human chromosomes, with segments of conserved gene order of varying lengths.

Conclusion

We obtained a radiation hybrid map of chicken chromosome 2. Comparison to the human genome indicated that most of the 8 groups of conserved synteny studied underwent internal rearrangements. The alignment of our RH map to the first draft of the chicken genome sequence assembly revealed a good agreement between both sets of data, indicative of a low error rate.     

Exploiting rice–sorghum synteny for targeted development of EST-SSRs to enrich the sorghum genetic linkage map

The sequencing and detailed comparative functional analysis of genomes of a number of select botanical models open new doors into comparative genomics among the angiosperms, with potential benefits for improvement of many orphan crops that feed large populations. In this study, a set of simple sequence repeat (SSR) markers was developed by mining the expressed sequence tag (EST) database of sorghum. Among the SSR-containing sequences, only those sharing considerable homology with rice genomic sequences across the lengths of the 12 rice chromosomes were selected. Thus, 600 SSR-containing sorghum EST sequences (50 homologous sequences on each of the 12 rice chromosomes) were selected, with the intention of providing coverage for corresponding homologous regions of the sorghum genome. Primer pairs were designed and polymorphism detection ability was assessed using parental pairs of two existing sorghum mapping populations. About 28% of these new markers detected polymorphism in this 4-entry panel. A subset of 55 polymorphic EST-derived SSR markers were mapped onto the existing skeleton map of a recombinant inbred population derived from cross N13 × E 36-1, which is segregating for Striga resistance and the stay-green component of terminal drought tolerance. These new EST-derived SSR markers mapped across all 10 sorghum linkage groups, mostly to regions expected based on prior knowledge of rice–sorghum synteny. The ESTs from which these markers were derived were then mapped in silico onto the aligned sorghum genome sequence, and 88% of the best hits corresponded to linkage-based positions. This study demonstrates the utility of comparative genomic information in targeted development of markers to fill gaps in linkage maps of related crop species for which sufficient genomic tools are not available.     

Towards an efficient computational mining approach to identify EST-SSR markers

Microsatellites are the markers of choice due to their high abundance reproducibility, degree of polymorphism and co-dominant nature. These are mainly used for studying the genetic variability in different species and Marker assisted selection. Expressed Sequence Tags (ESTs) serve as the main resource for Simple Sequence Repeats (SSRs). The computational approach for detecting SSRs and developing SSR markers from EST-SSRs is preferred over the conventional methods as it reduces time and cost to a great extent. The available EST sequence databases, various web interfaces and standalone tools provide the platform for an easy analysis of the EST sequences leading to the development of potential EST-SSR Markers. This paper is an overview of in silico approach to develop SSR Markers from the EST sequence using some of the most efficient tools that are available freely for academic purpose.     

A genome-wide set of 106 microsatellite markers for the blue tit (Cyanistes caeruleus)

We have characterized a set of 106 microsatellite markers in 26-127 individual blue tits (Cyanistes caeruleus), and assigned their location on the zebra finch (Taeniopygia guttata) and on the chicken (Gallus gallus) genome on the basis of sequence homology. Thirty-one markers are newly designed from zebra finch EST (expressed sequence tags) sequences, 22 markers were developed by others from EST sequences using different methods and the remaining 53 loci were previously designed or modified passerine markers. The 106 microsatellite markers are distributed over 26 and 24 chromosomes in the zebra finch and in the chicken genome respectively and the number of alleles varies between 2 and 49. Eight loci deviate significantly from Hardy-Weinberg equilibrium and show a high frequency of null alleles, and three pairs of markers located in the same chromosome appear to be in linkage disequilibrium. With the exception of these few loci, the polymorphic microsatellite markers presented here provide a useful genome-wide resource for population and evolutionary genetic studies of the blue tit, in addition to their potential utility in other passerine birds.     

Development of simple sequence repeat markers and diversity analysis in alfalfa (Medicago sativa L.)

Efficient and robust molecular markers are essential for molecular breeding in plant. Compared to dominant and bi-allelic markers, multiple alleles of simple sequence repeat (SSR) markers are particularly informative and superior in genetic linkage map and QTL mapping in autotetraploid species like alfalfa. The objective of this study was to enrich SSR markers directly from alfalfa expressed sequence tags (ESTs). A total of 12,371 alfalfa ESTs were retrieved from the National Center for Biotechnology Information. Total 774 SSR-containing ESTs were identified from 716 ESTs. On average, one SSR was found per 7.7 kb of EST sequences. Tri-nucleotide repeats (48.8 %) was the most abundant motif type, followed by di—(26.1 %), tetra—(11.5 %), penta—(9.7 %), and hexanucleotide (3.9 %). One hundred EST–SSR primer pairs were successfully designed and 29 exhibited polymorphism among 28 alfalfa accessions. The allele number per marker ranged from two to 21 with an average of 6.8. The PIC values ranged from 0.195 to 0.896 with an average of 0.608, indicating a high level of polymorphism of the EST–SSR markers. Based on the 29 EST–SSR markers, assessment of genetic diversity was conducted and found that Medicago sativa ssp. sativa was clearly different from the other subspecies. The high transferability of those EST–SSR markers was also found for relative species.     

Chromosome bin map of expressed sequence tags in homoeologous group 1 of hexaploid wheat and homoeology with rice and Arabidopsis

A total of 944 expressed sequence tags (ESTs) generated 2212 EST loci mapped to homoeologous group 1 chromosomes in hexaploid wheat (Triticum aestivum L.). EST deletion maps and the consensus map of group 1 chromosomes were constructed to show EST distribution. EST loci were unevenly distributed among chromosomes 1A, 1B, and 1D with 660, 826, and 726, respectively. The number of EST loci was greater on the long arms than on the short arms for all three chromosomes. The distribution of ESTs along chromosome arms was nonrandom with EST clusters occurring in the distal regions of short arms and middle regions of long arms. Duplications of group 1 ESTs in other homoeologous groups occurred at a rate of 35.5%. Seventy-five percent of wheat chromosome 1 ESTs had significant matches with rice sequences (E < or = e(-10)), where large regions of conservation occurred between wheat consensus chromosome 1 and rice chromosome 5 and between the proximal portion of the long arm of wheat consensus chromosome 1 and rice chromosome 10. Only 9.5% of group 1 ESTs showed significant matches to Arabidopsis genome sequences. The results presented are useful for gene mapping and evolutionary and comparative genomics of grasses.     

Development of microsatellite markers in autopolyploid sugarcane and comparative analysis of conserved microsatellites in sorghum and sugarcane

Sugarcane has become an increasingly important first-generation biofuel crop in tropical and subtropical regions. It has a large, complex, polyploid genome that has hindered the progress of genomic research and marker-assisted selection. Genetic mapping and ultimately genome sequence assembly require a large number of DNA markers. Simple sequence repeats (SSRs) are widely used in genetic mapping because of their abundance, high rates of polymorphism, and ease of use. The objectives of this study were to develop SSR markers for construction of a saturated genetic map and to characterize the frequency and distribution of SSRs in a polyploid genome. SSR markers were mined from expressed sequence tag (EST), reduced representation library genomic sequences, and bacterial artificial chromosome (BAC) sequences. A total of 5,675 SSR markers were surveyed in a segregating population. The overall successful amplification and polymorphic rates were 87.9 and 16.4%, respectively. The trinucleotide repeat motifs were most abundant, with tri- and hexanucleotide motifs being the most abundant for the ESTs. BAC and genomic SSRs were mostly AT-rich while the ESTs were relatively GC-rich due to codon bias. These markers were also aligned to the sorghum genome, resulting in 1,203 markers mapped in the sorghum genome. This set of SSRs conserved in sugarcane and sorghum would be the most informative for mapping quantitative trait loci in sugarcane and for comparative genomic analyses. This large collection of SSR markers is a valuable resource for sugarcane genomic research and crop improvement.     

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-SSRs的信息分析及分子标记筛选

昆虫EST资源库的扩充为开发新的分子标记提供了宝贵的资源。本研究对NCBI的EST数据库中来源于麦红吸浆虫Sitodiplosis mosellana唾腺的1 217条EST序列进行了unigene组装、 SSR信息分析和EST-SSR分子标记筛选。结果表明： 在1 047个unigenes中共找到141个SSR位点, 分布于106个（10.12%）unigenes中, 平均每3.49 kb出现一个SSR位点。在1~6碱基重复基元中, 1~3碱基是主要重复类型, 占总SSR的97.16%以上。A/T（31.21%）, AC/GT（15.60%）和AAC/GTT（9.22%）分别是单、 双和三碱基中占优势的重复基元类型。利用Primer Premier 5.0软件对查找的EST SSRs进行引物设计, 并以麦红吸浆虫基因组DNA为模板, 对从中选出的26对SSR引物进行多态性检测。结果有20对（76.92%）引物能扩增出清晰的目的条带, 并且其中9对（45%）引物表现出多态性。多态性分析结果表明, 从9对EST-SSR引物中, 共检测到51个等位基因, 平均每个位点含有等位基因数为5.67, 平均期望杂合度为0.65, 平均多态信息含量为0.60。本研究能够为今后麦红吸浆虫的种群遗传结构与遗传多样性研究提供帮助。     

褐飞虱EST资源的微卫星信息分析

表达序列标签（expressed sequence tags,ESTs）是开发微卫星标记的一个重要的资源。褐飞虱Nilaparvata lugens （Stål） EST序列的公布为开发EST-SSRs提供了宝贵的数据资源,本研究利用生物信息学对NCBI公共数据库中的37 398条褐飞虱ESTs序列进行EST-SSRs特征分析,得到全长为7 619 324 kb的无冗余EST 9 852条。按照3个不同的查找标准在这些序列中搜索SSR。查找结果显示：褐飞虱EST-SSRs主要重复基元以1～3碱基为主,占总EST-SSR的95％以上。在单碱基重复基元中,A/T是占优势的重复基元,在二相重复类型中,AG/CT重复基元出现的频率最多,而AAG/CTT是三相重复中占绝对优势的重复基元。在褐飞虱EST-SSRs中未查找到GC重复基元。以100 bp为参照,在3种查找标准下含有SSR的EST序列中两端侧翼序列均≥100 bp的序列分别为738,89和42个。通过分析褐飞虱EST-SSRs标记可以为褐飞虱和近缘种的SSR标记的开发提供信息,同时通过分析褐飞虱EST-SSRs的分布频率和分布特征可以为昆虫EST-SSRs的研究提供借鉴和参考。     

Development,characterization, validation,and mapping of SSRs derived from <Emphasis Type="Italic">Theobroma cacao</Emphasis> L.<Emphasis Type="Italic">–Moniliophthora perniciosa</Emphasis> interaction ESTs

In this study, we report results of the detection and analysis of SSR markers derived of cacao–Moniliophthora perniciosa expressed sequence tags (ESTs) in relation to cacao resistance to witches’ broom disease (WBD), and we compare the polymorphism of those ESTs (EST-simple sequence repeat (SSR)) with classical neutral SSR markers. A total of 3,487 ESTs was used in this investigation. SSRs were identified in 430 sequences: 277 from the resistant genotype TSH 1188 and 153 from the susceptible one Catongo, totalizing 505 EST-SSRs with three types of motifs: dinucleotides (72.1%), trinucleotides (27.3%), and tetranucleotides (0.6%). EST-SSRs were classified into 16 main categories; most of the EST-SSRs belonged to “Unknown function” and “No homology” categories (45.82%). A high frequency of SSRs was found in the 5’UTR and in the ORF (about 27%) and a low frequency was observed in the 3’UTR (about 8%). Forty-nine EST-SSR primers were designed and evaluated in 21 cacao accessions, 12 revealed polymorphism, having 47 alleles in total, with an average of 3.92 alleles per locus. On the other hand, the 11 genomic SSR markers revealed a total of 47 alleles, with an average of 5.22 alleles per locus. The association of EST-SSR with the genomic SSR enhanced the analysis of genetic distance among the genotypes. Among the 12 polymorphic EST-SSR markers, two were mapped on the F₂ Sca 6 × ICS 1 population reference for WBD resistance.     

Mining and characterization of EST-SSR markers for <Emphasis Type="Italic">Zingiber officinale</Emphasis> Roscoe with transferability to other species of Zingiberaceae

Zingiber officinale is a model spice herb, well known for its medicinal value. It is primarily a vegetatively propagated commercial crop. However, considerable diversity in its morphology, fiber content and chemoprofiles has been reported. The present study explores the utility of EST-derived markers in studying genetic diversity in different accessions of Z. officinale and their cross transferability within the Zingiberaceae family. A total of 38,115 ESTs sequences were assembled to generate 7850 contigs and 10,762 singletons. SSRs were searched in the unigenes and 515 SSR-containing ESTs were identified with a frequency of 1 SSR per 25.21 kb of the genome. These ESTs were also annotated using BLAST2GO. Primers were designed for 349 EST-SSRs and 25 primer pairs were randomly picked for EST SSR study. Out of these, 16 primer pairs could be optimized for amplification in different accessions of Z. officinale as well as other species belonging to Zingiberaceae. GES454, GES466, GES480 and GES486 markers were found to exhibit 100% cross-transferability among different members of Zingiberaceae.