EST derived SSR markers for comparative mapping in wheat and rice

Genome-wide identification of microsatellites for mapping, genetic diversity and cross-transferability in wheat (Triticum spp)

Wheat (Triticum aestivum L.) is a crucial global staple crop, and is consistently being improved to enhance yield, disease resistance, and quality traits. However, the development of molecular markers is a challenging task due to its hexaploid genome. Molecular marker system such as simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) are helpful for breeding, but SNP has limitations due to its development cost and its conversion to breeder markers. The study proposed an in-silico approach, by utilizing the low-cost transcriptome sequencing of two parental lines, 'TAC 75' and 'WH 1105', to identify polymorphic SSRs for mapping in a recombinant inbred line (RIL) population. This study introduces a new approach to bridge wheat genetics intricacies and next-generation sequencing potential. It presents a comprehensive genome-wide SSR distribution using IWGSC CS RefSeq v2.1 genome assembly and to identify 189 polymorphic loci through in-silico strategy. Of these, 54.76% showed polymorphism between parents, surpassing the traditional low polymorphic success rate. A RIL population screening validated these markers, demonstrating the fitness of identified markers through chi-square tests. The designed SSRs were also validated for genetic diversity analysis in a subset of 37 Indian wheat genotypes and cross-transferability in the wild/relative wheat species. In diversity analysis, a subset of 38 markers revealed 95 alleles (2.5 allele/locus), indicating substantial genetic variation. Population structure analysis unveiled three distinct groups, supported by phylogenetic and PCoA analyses. Further the polymorphic SSRs were also analyzed for SSR-gene association using gene ontology analysis. By utilizing the developing seed transcriptome data within parental lines, the study has enhanced the polymorphic SSR identification precision and facilitated in the RIL population. The undertaken study pioneers the use of transcriptome sequencing and genetic mapping to overcome challenges posed by the intricate wheat genome. This approach offers a cost-effective, less labour-intensive alternative to conventional methods, providing a platform for advancing wheat breeding research.

Occurrence of simple sequence repeats in cDNA sequences of safflower (Carthamus tinctorius) reveals the importance of SSR-containing genes for cell biology and dynamic response to environmental cues

Safflower (Carthamus tinctorius) is a diploid crop plant belonging to the family Asteraceae and is well known as one of important oilseed crops due to edible oil containing unsaturated fatty acids. In recent years it is gaining increased attention for food, pharmaceutical and industrial uses, and hence the updating its breeding methods is necessary. Genic simple sequence repeats (SSRs) in addition of being desire molecular markers, are supposed to influence gene function and the respective phenotype. This study aimed to identify SSRs in cDNA sequences and further analysis of the functional features of the SSR-containing genes to elucidate their role in biological and cellular processes. We identified 1,841 SSR regions in 1,667 cDNA sequences. Among all types of repeats, trinucleotide repeats were the most abundant (35.7%), followed by hexanucleotide (29.6%) and dinucleotide repeats (22.0%). Thirty five SSR primer pairs were validated by PCR reaction, detected a high rate of polymorphism (>57%) among safflower accessions, physically mapped on safflower genome and could clearly discriminate the cultivated accessions from wild relatives. The cDNA-derived SSR markers are suitable for evaluation of genetic diversity, linkage and association mapping studies and genome-based breeding programmes. Occurrence of SSR repeats in biologically-important classes of proteins such as kinases, transferases and transcription factors was inferred from functional analyses, which along with variability of their repeat copies, can endow the cell and whole organism the flexibility of facing with continuously changing environment, and indicate a structure-based evolution mechanism of the genome which acts as an up-to-dating tool for the cell and whole origanism, which is realized in GO terms such as involvement of most SSR-containing genes in biological, cellular and metabolic processes, especially in response to stimulus, response to stress, interaction to other organisms and defense responses.

Analysis of genetic diversity and relationships of Perilla frutescens using novel EST-SSR markers derived from transcriptome between wild-type and mutant Perilla

BACKGROUND

Perilla frutescens (Lamiaceae) is distributed in East Asia and is classified into var. frutescens and crispa. P. frutescens is multipurpose crop for human health because of a variety of secondary metabolites such as phenolic compound and essential oil. However, a lack of genetic information has hindered the development and utilization of Perilla genotypes.

METHODS AND RESULTS

This study was performed to develop expressed sequence tag-simple sequence repeat (EST-SSR) markers from P. frutescens var. crispa (wild type) and Antisperill (a mutant cultivar) and used them to assess the genetic diversity of, and relationships among, 94 P. frutescens genotypes. We obtained 65 Gb of sequence data comprising 632,970 transcripts by de novo RNA-sequencing. Of the 14,780 common SSRs, 102 polymorphic EST-SSRs were selected using in silico polymerase chain reaction (PCR). Overall, successful amplification from 58 EST-SSRs markers revealed remarkable genetic diversity and relationships among 94 P. frutescens genotypes. In total, 268 alleles were identified, with an average of 4.62 alleles per locus (range 2-11 alleles/locus). The average polymorphism information content (PIC) value was 0.50 (range 0.04-0.86). In phylogenetic and population structure analyses, the genotypes formed two major groups: Group I (var. crispa) and Group II (var. frutescens).

CONCLUSION

This results suggest that 58 novel EST-SSR markers derived from wild-type cultivar (var. crispa) and its mutant cultivar (Antisperill) have potential uses for population genetics and recombinant inbred line mapping analyses, which will provide comprehensive insights into the genetic diversity and relationship of P. frutescens.

Genome survey of Chinese fir (Cunninghamia lanceolata): Identification of genomic SSRs and demonstration of their utility in genetic diversity analysis

Chinese fir (Cunninghamia lanceolata) is an important coniferous species that accounts for 20–30% of the total commercial timber production in China. Though traditional breeding of Chinese fir has achieved remarkable success, molecular-assisted breeding has made little progress due to limited availability of genomic information. In this study, a survey of Chinese fir genome was performed using the Illumina HiSeq Xten sequencing platform. K-mer analysis indicated that Chinese fir has a large genome of approximately 11.6 Gb with 74.89% repetitive elements and is highly heterozygous. Meanwhile, its genome size was estimated to be 13.2 Gb using flow cytometry. A total of 778.02 Gb clean reads were assembled into 10,982,272 scaffolds with an N50 of 1.57 kb. In total, 362,193 SSR loci were detected with a frequency of 13.18 kb. Dinucleotide repeats were the most abundant (up to 73.6% of the total SSRs), followed by trinucleotide and tetranucleotide repeats. Forty-six polymorphic pairs were developed, and 298 alleles were successfully amplified from 199 Chinese fir clones. The average PIC value was 0.53, indicating that the identified genomic SSR (gSSR) markers have a high degree of polymorphism. In addition, these breeding resources were divided into three groups, and a limited gene flow existed among these inferred groups.

Antimicrobial and antioxidant activities of hydrocarbon and oxygenated monoterpenes against some foodborne pathogens through in vitro and in silico studies

The present work describes the antimicrobial action of 25 monoterpenes (six hydrocarbons, five ketones, two aldehydes, six alcohols and six acetate analogues) against Gram-negative Escherichia coli and Gram-positive Staphylococcus aureus and antifungal activity against Aspergillus flavus. The antibacterial activity was evaluated by broth microdilution technique as a minimum inhibitory concentration (MIC) and the antifungal activity was performed by mycelia radial growth technique as the effective concentration causing 50% inhibition of the mycelial growth (EC₅₀). The results showed that thymol and α-terpineol were the most potent against E. coli (MIC = 45 and 55 mg/L, respectively) and S. aureus (MIC = 135 and 225 mg/L, respectively). The results also showed that thymol displayed the maximum antifungal action against A. flavus with EC₅₀ 20 mg/L. Furthermore, the antioxidant activity was determined using N,N-dimethyl-1,4-phenylenediamine (DMPD) and the results showed that geraniol were the most potent compound (IC₅₀ = 19 mg/L). Molecular docking studies indicated that the compounds displayed different binding interactions with the amino acid residues at the catalytic sites of N5-carboxyaminoimidazole synthetase and oxysterol binding protein Osh4 enzymes. Non-covalent interactions including van der Waals, hydrogen bonding as well as hydrophobic were observed between the compounds and the enzymes. A significant relationship was found between the docking score and the biological activity of the tested monoterpenes compared to the ceftriaxone and carbendazim as standard bactericide and fungicide, respectively. In silico ADMET properties were also performed and displayed potential for the development of promising antimicrobial agents. For these reasons, these compounds may be considered as potential ecofriendly alternatives in food preservation to delay or prevent the microbial infection and prolong the shelf life of food products.

Identification, analysis and development of salt responsive candidate gene based SSR markers in wheat

BACKGROUND

Salinity severely limits wheat production in many parts of the world. Development of salt tolerant varieties represents the most practical option for enhancing wheat production from these areas. Application of marker assisted selection may assist in fast tracking development of salt tolerant wheat varieties. However, SSR markers available in the public domain are not specifically targeted to functional regions of wheat genome, therefore large numbers of these need to be analysed for identification of markers associated with traits of interest. With the availability of a fully annotated wheat genome assembly, it is possible to develop SSR markers specifically targeted to genic regions. We performed extensive analysis to identify candidate gene based SSRs and assessed their utility in characterizing molecular diversity in a panel of wheat genotypes.

RESULTS

Our analysis revealed, 161 SSR motifs in 94 salt tolerance candidate genes of wheat. These SSR motifs were nearly equally distributed on the three wheat sub-genomes; 29.8% in A, 35.7% in B and 34.4% in D sub-genome. The maximum number of SSR motifs was present in exons (31.1%) followed by promoters (29.8%), 5'UTRs (21.1%), introns (14.3%) and 3'UTRs (3.7%). Out of the 65 candidate gene based SSR markers selected for validation, 30 were found polymorphic based on initial screening and employed for characterizing genetic diversity in a panel of wheat genotypes including salt tolerant and susceptible lines. These markers generated an average of 2.83 alleles/locus. Phylogenetic analysis revealed four clusters. Salt susceptible genotypes were mainly represented in clusters I and III, whereas high and moderate salt tolerant genotypes were distributed in the remaining two clusters. Population structure analysis revealed two sub-populations, sub-population 1 contained the majority of salt tolerant whereas sub-population 2 contained majority of susceptible genotypes. Moreover, we observed reasonably higher transferability of SSR markers to related wheat species.

CONCLUSION

We have developed salt responsive gene based SSRs in wheat for the first time. These were highly useful in unravelling functional diversity among wheat genotypes with varying responses to salt stress. The identified gene based SSR markers will be valuable genomic resources for genetic/association mapping of salinity tolerance in wheat.

Characterization and development of EST-SSR markers to study the genetic diversity and populations analysis of Jerusalem artichoke (Helianthus tuberosus L.)

In recent years, Jerusalem artichoke has received widespread attention as a novel source of sugar, biofuel, and animal feed. Currently, only few gDNA-SSRs derived from sunflower were verified in the Jerusalem artichoke; therefore, it is particularly important to develop SSR primer markers that belonged to Jerusalem artichoke resources. Using EST data to develop EST-SSR markers is simple and effective. In order to understand the general characteristics of SSR markers in Jerusalem artichoke EST sequences and accelerate the use of SSR markers in Jerusalem artichoke research. This study used 40,370 sequenced unigene fragments and MISA software to identify SSR loci. The 48 pairs of EST-SSR primers assessed for the identification of 45 varieties of Jerusalem artichoke. Cluster, genetic diversity parameters and AMOVA analysis was conducted using the genetic similarity coefficient, revealing genetic differences between 48 genetic material. A total of 1204 SSR loci were identified with 13 different types of repeats, distributed among 1020 EST sequences, of which trinucleotide repeats were the most common, accounting for 38.21% of the total SSR loci. Among the 44 repeat motifs, AG/CT, AAG/CTT, and ATC/ATG motifs had the highest frequencies, accounting for 22.45, 14.71, and 7.84% of all motifs, respectively. From these sequences, 48 pairs of EST-SSR primers were designed, and 22 primer pairs for loci with high polymorphism were selected to analyze the genetic diversity of 45 Jerusalem artichoke germplasm sources. The results indicated that the variation range of the effective number of alleles for 22 primers ranged between 1.7502 and 4.5660. The Shannon's information index ranged between 0.6200 and 1.6423. The variation range of PIC ranged between 0.3121 and 0.6662 with an average of 0.5184. Cluster analysis was conducted using the genetic similarity coefficient, revealing significant genetic differences between Asian and European genetic material. Cluster analysis revealed a relationship between the genotypes and geographic origins of the Jerusalem artichoke. The results of AMOVA as well as the genetic identity and genetic distance in the Jerusalem artichoke population showed that there presented certain genetic heterogeneity in Jerusalem artichoke genetic structure of 45 samples from seven different geographic populations. The Jerusalem artichoke EST-SSR marker system established in this study provides an effective molecular marker system for future research focused on Jerusalem artichoke genetic diversity and the breeding of new varieties.

Systemic Identification of Hevea brasiliensis EST-SSR Markers and Primer Screening

This research aimed to systematically identify and preliminarily validate the Hevea brasiliensis expressed sequence tag (EST) information using Simple Sequence Repeat (SSR) and provide evidence for further development of SSR molecular marker. The definition of general SSR features of Hevea EST splicing sequences and development of SSR primers founded the basis of diversity analysis and variety identification for Hevea tree resource. 1134 SSR loci were identified in the EST splicing sequence and distributed in 840 Unigene. The occurrence rate of SSR loci was 23.9%, and the average distribution distance of EST-SSR was 2.59 kb. The major repeat type was mononucleotide repeat motif, which accounted for 38.89%, while the corresponding value was 36.95% for dinucleotide repeat motif and 18.17% for trinucleotide repeat motif; the proportion of other motifs was only 5.99%. The superior repeat motifs for mononucleotide, dinucleotide, and trinucleotide were A/T, AG/CT, and AAG/CTT, respectively. 739 pair of primers were designed for 1134 SSR loci. PCR amplification was performed on Hevea Reyan5-11, Reyan87-6-47, and PR107, and 180 pairs of primers were selected which were able to amplify polymorphism bands.

Expressed Sequence Tags Analysis and Design of Simple Sequence Repeats Markers from a Full-Length cDNA Library in Perilla frutescens (L.)

Perilla frutescens is valuable as a medicinal plant as well as a natural medicine and functional food. However, comparative genomics analyses of P. frutescens are limited due to a lack of gene annotations and characterization. A full-length cDNA library from P. frutescens leaves was constructed to identify functional gene clusters and probable EST-SSR markers via analysis of 1,056 expressed sequence tags. Unigene assembly was performed using basic local alignment search tool (BLAST) homology searches and annotated Gene Ontology (GO). A total of 18 simple sequence repeats (SSRs) were designed as primer pairs. This study is the first to report comparative genomics and EST-SSR markers from P. frutescens will help gene discovery and provide an important source for functional genomics and molecular genetic research in this interesting medicinal plant.

CNMS: The preferred genic markers for comparative genomic, molecular phylogenetic, functional genetic diversity and differential gene regulatory expression analyses in chickpea

The intra/inter-genomic comparative mapping-based phylogenetic footprinting identified 5 paralogous and 656 orthologous genome-wide CNMS markers in the upstream sequences of chickpea genes. These CNMS markers revealed a high-degree of gene-based syntenic relationship between chickpea and Medicago genomes while minimum between chickpea and Vitis genomes. The time of divergence and duplication estimated using CNMS markers highlight the expected phylogenetic relationships between chickpea and six dicot (legume) species as well as occurrence of ancient genome (approximately 53 Mya) with small-scale recent segmental (approximately 10 Mya) duplication events in chickpea. A wider level of functional molecular diversity (14 to 88 percent) and admixed population genetic structure was detected among desi, kabuli and wild genotypes by genic CNMS markers at a genome-wide scale suggesting their utility in large-scale genetic analysis in chickpea. The subfunctionalization at the cis-regulatory element region and TFBS (transcription factor binding site) motif levels in the upstream sequences of CNMS marker-associated orthologous genes than the paralogues was predominant. Functional constraint might have considerable effect on these CNMScontaining regulatory elements controlling consistent orthologous gene expression in dicots. A rapid subfunctionalization based on diverge differential expression of paralogous CNMS marker-associated genes particularly those that underwent recent small-scale segmental duplication events in chickpea was apparent. The differential regulation of expression and subfunctionalization potential of Ultra CNMS marker-associated genes suggest their utility in deciphering the complex gene regulatory function as well as identification and targeted mapping of potential genes/QTLs governing vital agronomic traits in chickpea. The gene-based CNMS markers with desirable inherent genetic attributes like higher degree of comparative genome mapping, functional genetic diversity and differential gene regulatory expression potential can significantly propel the genomics-assisted chickpea crop improvement.

Molecular breeding in Brassica for salt tolerance: importance of microsatellite (SSR) markers for molecular breeding in Brassica

Salinity is one of the important abiotic factors for any crop management in irrigated as well as rainfed areas, which leads to poor harvests. This yield reduction in salt affected soils can be overcome by improving salt tolerance in crops or by soil reclamation. Salty soils can be reclaimed by leaching the salt or by cultivation of salt tolerance crops. Salt tolerance is a quantitative trait controlled by several genes. Poor knowledge about mechanism of its inheritance makes slow progress in its introgression into target crops. Brassica is known to be a good reclamation crop. Inter and intra specific variation within Brassica species shows potential of molecular breeding to raise salinity tolerant genotypes. Among the various molecular markers, SSR markers are getting high attention, since they are randomly sparsed, highly variable and show co-dominant inheritance. Furthermore, as sequencing techniques are improving and softwares to find SSR markers are being developed, SSR markers technology is also evolving rapidly. Comparative SSR marker studies targeting Arabidopsis thaliana and Brassica species which lie in the same family will further aid in studying the salt tolerance related QTLs and subsequent identification of the "candidate genes" and finding out the origin of important QTLs. Although, there are a few reports on molecular breeding for improving salt tolerance using molecular markers in Brassica species, usage of SSR markers has a big potential to improve salt tolerance in Brassica crops. In order to obtain best harvests, role of SSR marker driven breeding approaches play important role and it has been discussed in this review especially for the introgression of salt tolerance traits in crops.

Identification and analysis of novel salt responsive candidate gene based SSRs (cgSSRs) from rice (Oryza sativa L.)

BACKGROUND

Majority of the Asian people depend on rice for nutritional energy. Rice cultivation and yield are severely affected by soil salinity stress worldwide. Marker assisted breeding is a rapid and efficient way to develop improved variety for salinity stress tolerance. Genomic microsatellite markers are an elite group of markers, but there is possible uncertainty of linkage with the important genes. In contrast, there are better possibilities of linkage detection with important genes if SSRs are developed from candidate genes. To the best of our knowledge, there is no such report on SSR markers development from candidate gene sequences in rice. So the present study was aimed to identify and analyse SSRs from salt responsive candidate genes of rice.

RESULTS

In the present study, based on the comprehensive literature survey, we selected 220 different salt responsive genes of rice. Out of them, 106 genes were found to contain 180 microsatellite loci with, tri-nucleotide motifs (56%) being most abundant, followed by di-(41%) and tetra nucleotide (2.8%) motifs. Maximum loci were found in the coding sequences (37.2%), followed by in 5'UTR (26%), intron (21.6%) and 3'UTR (15%). For validation, 19 primer sets were evaluated to detect polymorphism in diversity analysis among the two panels consisting of 17 salt tolerant and 17 susceptible rice genotypes. Except one, all primer sets exhibited polymorphic nature with an average of 21.8 alleles/primer and with a mean PIC value of 0.28. Calculated genetic similarity among genotypes was ranged from 19%-89%. The generated dendrogram showed 3 clusters of which one contained entire 17 susceptible genotypes and another two clusters contained all tolerant genotypes.

CONCLUSION

The present study represents the potential of salt responsive candidate gene based SSR (cgSSR) markers to be utilized as novel and remarkable candidate for diversity analysis among rice genotypes differing in salinity response.

A genetic linkage map of sole (Solea solea): a tool for evolutionary and comparative analyses of exploited (flat)fishes

Linkage maps based on markers derived from genes are essential evolutionary tools for commercial marine fish to help identify genomic regions associated with complex traits and subject to selective forces at play during exploitation or selective breeding. Additionally, they allow the use of genomic information from other related species for which more detailed information is available. Sole (solea solea L.) is a commercially important flatfish species in the North Sea, subject to overexploitation and showing evidence of fisheries-induced evolutionary changes in growth- and maturation-related traits. Sole would definitely benefit from a linkage map to better understand how evolution has shaped its genome structure. This study presents a linkage map of sole based on 423 single nucleotide polymorphisms derived from expressed sequence tags and 8 neutral microsatellite markers. The total map length is 1233.8 cM and consists of 38 linkage groups with a size varying between 0 to 92.1 cM. Being derived from expressed sequence tags allowed us to align the map with the genome of four model fish species, namely medaka (Oryzias latipes), Nile tilapia (Oreochromis niloticus), three-spined stickleback (Gasterosteus aculeatus) and green spotted pufferfish (Tetraodon nigroviridis). This comparison revealed multiple conserved syntenic regions with all four species, and suggested that the linkage groups represent 21 putative sole chromosomes. The map was also compared to the linkage map of turbot (Scophthalmus maximus), another commercially important flatfish species and closely related to sole. For all putative sole chromosomes (except one) a turbot homolog was detected, confirming the even higher degree of synteny between these two flatfish species.

Characterization of EST-SSR markers in durum wheat EST library and functional analysis of SSR-containing EST fragments

The goal of this study is to identify characterization of expressed sequence tag (EST)-simple sequence repeats (SSR) markers from EST library of durum wheat and functional analysis of SSR-containing EST sequences for application in comparative genomics and breeding. 19,141 sequences were analyzed among which 18,937 ESTs were selected. Consistent with MISA results, 313 EST-SSRs were yielded. The final EST-SSRs were compared to the GenBank non-redundant database using BLASTX and classified based on these functions. Results indicated that the perfect EST-SSRs are the most frequent. The TTG/CTG imperfect EST-SSR had gamma-gliadin putative function that can be appropriate for durum wheat. Also, the mononucleotides and trinucleotides were the most frequent. Findings suggested that the identified EST-SSRs could be categorized into 83 types. Motifs TTG in trinucleotides and TC in dinucleotides had the highest frequency. TTG is the new motif in durum wheat identified in this study. We identified new EST-SSRs with more than trinucleotide and detected motifs that have potential to code amino acids. Arginine was the most frequent amino acid. Enzymes had the highest frequency among predicted functions. EST-SSRs have been identified in this study can be used for developing ESS-SSR-based detection tool for durum wheat in future studies and will be a useful resource for molecular breeding, genetics, genomics, and environmental stress studies. Motifs coding amino acids could be used as a new source of functional markers and biological study. In addition to, designed new PCR primer pairs are new resources for to identify useful alleles in transcription factors, storage proteins, and enzymes which incorporated them again into the cultivated material.

De novo assembly and annotation of Salvia splendens transcriptome using the Illumina platform

BACKGROUND

As an important perennial herbaceous flower, Salvia splendens possesses high ornamental value. Understanding its branching processes may help scientists select the best plant type. Although Salvia splendens is a frequently-used horticultural flower, only limited transcriptomic or genomic research is available in public databases. In the present study, we, for the first time, constructed a comprehensive dataset for Salvia splendens through de novo high-throughput transcriptome sequencing.

METHODOLOGY/PRINCIPAL FINDINGS

We performed de novo transcriptome sequencing on two different branching type plants (Strain 35 and Cailinghong) using the Illumina paired-end sequencing technology. For Strain 35, a total of 16,488,829 reads were generated and assembled into 38,498 unigenes, with a mean length of approximately 779 bp. For Cailinghong, 16,464,713 reads were generated and assembled into 34,302 unigenes, with a mean length of approximately 812 bp. Moreover, a total of 49,310 unigenes for Salvia splendens were identified, among them 33,925 (68.80%) were annotated in the non-redundant NCBI database, 25,371 (51.45%) were annotated in the Swiss-Prot database, while 24,888 (50.47%) and 9,896 (20.07%) unigenes were assigned to gene ontology categories and clusters of orthologous groups, respectively. Using the Kyoto Encyclopedia of Genes and Genomes pathway database, we identified 134 differently expressed unigenes between Strain 35 and Cailinghong, and then these unigenes were mapped to 79 pathways. In addition, we detected 2,453 simple sequence repeats (SSRs).

CONCLUSIONS

We obtained a comprehensive transcriptomic information from this work and provided a valuable resource of transcript sequences of Salvia splendens in public databases. Moreover, some candidate genes potentially involved in branching were identified. Furthermore, numerous obtained SSRs might contribute to marker-assisted selection. These data could be further utilized in functional genomics studies on Salvia splendens.

Transcriptome sequencing and De Novo analysis of Youngia japonica using the illumina platform

Youngia japonica, a weed species distributed worldwide, has been widely used in traditional Chinese medicine. It is an ideal plant for studying the evolution of Asteraceae plants because of its short life history and abundant source. However, little is known about its evolution and genetic diversity. In this study, de novo transcriptome sequencing was conducted for the first time for the comprehensive analysis of the genetic diversity of Y. japonica. The Y. japonica transcriptome was sequenced using Illumina paired-end sequencing technology. We produced 21,847,909 high-quality reads for Y. japonica and assembled them into contigs. A total of 51,850 unigenes were identified, among which 46,087 were annotated in the NCBI non-redundant protein database and 41,752 were annotated in the Swiss-Prot database. We mapped 9,125 unigenes onto 163 pathways using the Kyoto Encyclopedia of Genes and Genomes Pathway database. In addition, 3,648 simple sequence repeats (SSRs) were detected. Our data provide the most comprehensive transcriptome resource currently available for Y. japonica. C4 photosynthesis unigenes were found in the biological process of Y. japonica. There were 5596 unigenes related to defense response and 1344 ungienes related to signal transduction mechanisms (10.95%). These data provide insights into the genetic diversity of Y. japonica. Numerous SSRs contributed to the development of novel markers. These data may serve as a new valuable resource for genomic studies on Youngia and, more generally, Cichoraceae.

Isolation and characterization of a wheat--Psathyrostachys huashanica 'Keng' 3Ns disomic addition line with resistance to stripe rust

We isolated a wheat germplasm line, 22-2, which was derived from common wheat (Triticum aestivum '7182') and Psathyrostachys huashanica 'Keng' (2n = 2x = 14, NsNs). Genomic composition and homoeologous relationships of 22-2 was analyzed using cytology, genomic in situ hybridization (GISH), EST-SSR, and EST-STS to characterize the alien chromatin in the transfer line. The cytological investigations showed that the chromosome number and configuration were 2n = 44 = 22 II. Mitotic and meiotic GISH using P. huashanica genomic DNA as the probe indicated that 22-2 contained a pair of P. huashanica chromosomes. The genomic affinities of the introduced P. huashanica chromosomes were determined by EST-SSR and EST-STS using multiple-loci markers from seven wheat homoeologous groups between the parents and addition line. One EST-SSR and 17 EST-STS markers, which were located on the homoeologous group 3 chromosomes of wheat, amplified polymorphic bands in 22-2 that were unique to P. huashanica. Thus, these markers suggested that the introduced Ns chromosome pair belonged to homoeologous group 3, so we designated 22-2 as a 3Ns disomic addition line. Based on disease reaction to mixed races (CYR31, CYR32, and Shuiyuan14) of stripe rust in the adult stages, 22-2 was found to have high resistance to stripe rust, which was possibly derived from its P. huashanica parent. Consequently, the new disomic addition line 22-2 could be a valuable donor source for wheat improvement depending on the excellent agronomic traits, especially, the introduction of novel disease resistance genes into wheat during breeding programs.

Next generation characterisation of cereal genomes for marker discovery

Cereal crops form the bulk of the world’s food sources, and thus their importance cannot be understated. Crop breeding programs increasingly rely on high-resolution molecular genetic markers to accelerate the breeding process. The development of these markers is hampered by the complexity of some of the major cereal crop genomes, as well as the time and cost required. In this review, we address current and future methods available for the characterisation of cereal genomes, with an emphasis on faster and more cost effective approaches for genome sequencing and the development of markers for trait association and marker assisted selection (MAS) in crop breeding programs.

Development and characterization of a Psathyrostachys huashanica Keng 7Ns chromosome addition line with leaf rust resistance

The aim of this study was to characterize a Triticum aestivum-Psathyrostachys huashanica Keng (2n = 2x = 14, NsNs) disomic addition line 2-1-6-3. Individual line 2-1-6-3 plants were analyzed using cytological, genomic in situ hybridization (GISH), EST-SSR, and EST-STS techniques. The alien addition line 2-1-6-3 was shown to have two P. huashanica chromosomes, with a meiotic configuration of 2n = 44 = 22 II. We tested 55 EST-SSR and 336 EST-STS primer pairs that mapped onto seven different wheat chromosomes using DNA from parents and the P. huashanica addition line. One EST-SSR and nine EST-STS primer pairs indicated that the additional chromosome of P. huashanica belonged to homoeologous group 7, the diagnostic fragments of five EST-STS markers (BE404955, BE591127, BE637663, BF482781 and CD452422) were cloned, sequenced and compared. The results showed that the amplified polymorphic bands of P. huashanica and disomic addition line 2-1-6-3 shared 100% sequence identity, which was designated as the 7Ns disomic addition line. Disomic addition line 2-1-6-3 was evaluated to test the leaf rust resistance of adult stages in the field. We found that one pair of the 7Ns genome chromosomes carried new leaf rust resistance gene(s). Moreover, wheat line 2-1-6-3 had a superior numbers of florets and grains per spike, which were associated with the introgression of the paired P. huashanica chromosomes. These high levels of disease resistance and stable, excellent agronomic traits suggest that this line could be utilized as a novel donor in wheat breeding programs.

Genetic basis for spontaneous hybrid genome doubling during allopolyploid speciation of common wheat shown by natural variation analyses of the paternal species

The complex process of allopolyploid speciation includes various mechanisms ranging from species crosses and hybrid genome doubling to genome alterations and the establishment of new allopolyploids as persisting natural entities. Currently, little is known about the genetic mechanisms that underlie hybrid genome doubling, despite the fact that natural allopolyploid formation is highly dependent on this phenomenon. We examined the genetic basis for the spontaneous genome doubling of triploid F1 hybrids between the direct ancestors of allohexaploid common wheat (Triticum aestivum L., AABBDD genome), namely Triticumturgidum L. (AABB genome) and Aegilopstauschii Coss. (DD genome). An Ae. tauschii intraspecific lineage that is closely related to the D genome of common wheat was identified by population-based analysis. Two representative accessions, one that produces a high-genome-doubling-frequency hybrid when crossed with a T. turgidum cultivar and the other that produces a low-genome-doubling-frequency hybrid with the same cultivar, were chosen from that lineage for further analyses. A series of investigations including fertility analysis, immunostaining, and quantitative trait locus (QTL) analysis showed that (1) production of functional unreduced gametes through nonreductional meiosis is an early step key to successful hybrid genome doubling, (2) first division restitution is one of the cytological mechanisms that cause meiotic nonreduction during the production of functional male unreduced gametes, and (3) six QTLs in the Ae. tauschii genome, most of which likely regulate nonreductional meiosis and its subsequent gamete production processes, are involved in hybrid genome doubling. Interlineage comparisons of Ae. tauschii's ability to cause hybrid genome doubling suggested an evolutionary model for the natural variation pattern of the trait in which non-deleterious mutations in six QTLs may have important roles. The findings of this study demonstrated that the genetic mechanisms for hybrid genome doubling could be studied based on the intrinsic natural variation that exists in the parental species.

Mapping QTL for grain yield and other agronomic traits in post-rainy sorghum [Sorghum bicolor (L.) Moench]

Sorghum, a cereal of economic importance ensures food and fodder security for millions of rural families in the semi-arid tropics. The objective of the present study was to identify and validate quantitative trait loci (QTL) for grain yield and other agronomic traits using replicated phenotypic data sets from three post-rainy dry sorghum crop seasons involving a mapping population with 245 F9 recombinant inbred lines derived from a cross of M35-1 × B35. A genetic linkage map was constructed with 237 markers consisting of 174 genomic, 60 genic and 3 morphological markers. The QTL analysis for 11 traits following composite interval mapping identified 91 QTL with 5-12 QTL for each trait. QTL detected in the population individually explained phenotypic variation between 2.5 and 30.3 % for a given trait and six major genomic regions with QTL effect on multiple traits were identified. Stable QTL across seasons were identified. Of the 60 genic markers mapped, 21 were found at QTL peak or tightly linked with QTL. A gene-based marker XnhsbSFCILP67 (Sb03g028240) on SBI-03, encoding indole-3-acetic acid-amido synthetase GH3.5, was found to be involved in QTL for seven traits. The QTL-linked markers identified for 11 agronomic traits may assist in fine mapping, map-based gene isolation and also for improving post-rainy sorghum through marker-assisted breeding.

Development of microsatellite markers by transcriptome sequencing in two species of Amorphophallus (Araceae)

Background

Amorphophallus is a genus of perennial plants widely distributed in the tropics or subtropics of West Africa and South Asia. Its corms contain a high level of water-soluble glucomannan; therefore, it has long been used as a medicinal herb and food source. Genetic studies of Amorphophallus have been hindered by a lack of genetic markers. A large number of molecular markers are required for genetic diversity study and improving disease resistance in Amorphophallus. Here, we report large scale of transcriptome sequencing of two species: Amorphophallus konjac and Amorphophallus bulbifer using deep sequencing technology, and microsatellite (SSR) markers were identified based on these transcriptome sequences.

Results

cDNAs of A. konjac and A. bulbifer were sequenced using Illumina HiSeq™ 2000 sequencing technology. A total of 135,822 non-redundant unigenes were assembled from about 9.66 gigabases, and 19,596 SSRs were identified in 16,027 non-redundant unigenes. Di-nucleotide SSRs were the most abundant motif (61.6%), followed by tri- (30.3%), tetra- (5.6%), penta- (1.5%), and hexa-nucleotides (1%) repeats. The top di- and tri-nucleotide repeat motifs included AG/CT (45.2%) and AGG/CCT (7.1%), respectively. A total of 10,754 primer pairs were designed for marker development. Of these, 320 primers were synthesized and used for validation of amplification and assessment of polymorphisms in 25 individual plants. The total of 275 primer pairs yielded PCR amplification products, of which 205 were polymorphic. The number of alleles ranged from 2 to 14 and the polymorphism information content valued ranged from 0.10 to 0.90. Genetic diversity analysis was done using 177 highly polymorphic SSR markers. A phenogram based on Jaccard’s similarity coefficients was constructed, which showed a distinct cluster of 25 Amorphophallus individuals.

Conclusion

A total of 10,754 SSR markers have been identified in Amorphophallus using transcriptome sequencing. One hundred and seventy-seven polymorphic markers were successfully validated in 25 individuals. The large number of genetic markers developed in the present study should contribute greatly to research into genetic diversity and germplasm characterization in Amorphophallus.

Development of eSSR-Markers in Setaria italica and Their Applicability in Studying Genetic Diversity, Cross-Transferability and Comparative Mapping in Millet and Non-Millet Species

Foxtail millet (Setariaitalica L.) is a tractable experimental model crop for studying functional genomics of millets and bioenergy grasses. But the limited availability of genomic resources, particularly expressed sequence-based genic markers is significantly impeding its genetic improvement. Considering this, we attempted to develop EST-derived-SSR (eSSR) markers and utilize them in germplasm characterization, cross-genera transferability and in silico comparative mapping. From 66,027 foxtail millet EST sequences 24,828 non-redundant ESTs were deduced, representing ~16 Mb, which revealed 534 (~2%) eSSRs in 495 SSR containing ESTs at a frequency of 1/30 kb. A total of 447 pp were successfully designed, of which 327 were mapped physically onto nine chromosomes. About 106 selected primer pairs representing the foxtail millet genome showed high-level of cross-genera amplification at an average of ~88% in eight millets and four non-millet species. Broad range of genetic diversity (0.02–0.65) obtained in constructed phylogenetic tree using 40 eSSR markers demonstrated its utility in germplasm characterizations and phylogenetics. Comparative mapping of physically mapped eSSR markers showed considerable proportion of sequence-based orthology and syntenic relationship between foxtail millet chromosomes and sorghum (~68%), maize (~61%) and rice (~42%) chromosomes. Synteny analysis of eSSRs of foxtail millet, rice, maize and sorghum suggested the nested chromosome fusion frequently observed in grass genomes. Thus, for the first time we had generated large-scale eSSR markers in foxtail millet and demonstrated their utility in germplasm characterization, transferability, phylogenetics and comparative mapping studies in millets and bioenergy grass species.

Exploiting the transcriptome of Euphrates Poplar, Populus euphratica (Salicaceae) to develop and characterize new EST-SSR markers and construct an EST-SSR database

BACKGROUND

Microsatellite markers or Simple Sequence Repeats (SSRs) are the most popular markers in population/conservation genetics. However, the development of novel microsatellite markers has been impeded by high costs, a lack of available sequence data and technical difficulties. New species-specific microsatellite markers were required to investigate the evolutionary history of the Euphratica tree, Populus euphratica, the only tree species found in the desert regions of Western China and adjacent Central Asian countries.

METHODOLOGY/PRINCIPAL FINDINGS

A total of 94,090 non-redundant Expressed Sequence Tags (ESTs) from P. euphratica comprising around 63 Mb of sequence data were searched for SSRs. 4,202 SSRs were found in 3,839 ESTs, with 311 ESTs containing multiple SSRs. The most common motif types were trinucleotides (37%) and hexanucleotides (33%) repeats. We developed primer pairs for all of the identified EST-SSRs (eSSRs) and selected 673 of these pairs at random for further validation. 575 pairs (85%) gave successful amplification, of which, 464 (80.7%) were polymorphic in six to 24 individuals from natural populations across Northern China. We also tested the transferability of the polymorphic eSSRs to nine other Populus species. In addition, to facilitate the use of these new eSSR markers by other researchers, we mapped them onto Populus trichocarpa scaffolds in silico and compiled our data into a web-based database (http://202.205.131.253:8080/poplar/resources/static_page/index.html).

CONCLUSIONS

The large set of validated eSSRs identified in this work will have many potential applications in studies on P. euphratica and other poplar species, in fields such as population genetics, comparative genomics, linkage mapping, QTL, and marker-assisted breeding. Their use will be facilitated by their incorporation into a user-friendly web-based database.

Pearl millet [Pennisetum glaucum (L.) R. Br.] consensus linkage map constructed using four RIL mapping populations and newly developed EST-SSRs

BACKGROUND

Pearl millet [Pennisetum glaucum (L.) R. Br.] is a widely cultivated drought- and high-temperature tolerant C4 cereal grown under dryland, rainfed and irrigated conditions in drought-prone regions of the tropics and sub-tropics of Africa, South Asia and the Americas. It is considered an orphan crop with relatively few genomic and genetic resources. This study was undertaken to increase the EST-based microsatellite marker and genetic resources for this crop to facilitate marker-assisted breeding.

RESULTS

Newly developed EST-SSR markers (99), along with previously mapped EST-SSR (17), genomic SSR (53) and STS (2) markers, were used to construct linkage maps of four F7 recombinant inbred populations (RIP) based on crosses ICMB 841-P3 × 863B-P2 (RIP A), H 77/833-2 × PRLT 2/89-33 (RIP B), 81B-P6 × ICMP 451-P8 (RIP C) and PT 732B-P2 × P1449-2-P1 (RIP D). Mapped loci numbers were greatest for RIP A (104), followed by RIP B (78), RIP C (64) and RIP D (59). Total map lengths (Haldane) were 615 cM, 690 cM, 428 cM and 276 cM, respectively. A total of 176 loci detected by 171 primer pairs were mapped among the four crosses. A consensus map of 174 loci (899 cM) detected by 169 primer pairs was constructed using MergeMap to integrate the individual linkage maps. Locus order in the consensus map was well conserved for nearly all linkage groups. Eighty-nine EST-SSR marker loci from this consensus map had significant BLAST hits (top hits with e-value ≤ 1E-10) on the genome sequences of rice, foxtail millet, sorghum, maize and Brachypodium with 35, 88, 58, 48 and 38 loci, respectively.

CONCLUSION

The consensus map developed in the present study contains the largest set of mapped SSRs reported to date for pearl millet, and represents a major consolidation of existing pearl millet genetic mapping information. This study increased numbers of mapped pearl millet SSR markers by >50%, filling important gaps in previously published SSR-based linkage maps for this species and will greatly facilitate SSR-based QTL mapping and applied marker-assisted selection programs.

Limitations, progress and prospects of application of biotechnological tools in improvement of bamboo-a plant with extraordinary qualities

Bamboos (family Poaceae) are the most beautiful and useful plants on the Earth, mainly found in the tropical and sub-tropical regions of the world. Bamboos are fast growing and early maturing, but lack of proper management of bamboo resources is leading to rapid reduction of the existing bamboosetum. Bamboo propagation through seeds is limited due to long flowering cycle of upto 120 years, seed sterility and short seed viability. Infrequent and unpredictable flowering events coupled with peculiar monocarpic behaviour i.e. flowering once before culm death, and extensive genome polyploidization are additional challenges for this woody group. Similarly, vegetative propagation by cuttings, offsets and rhizomes are also inadequate to cope up with the demand of planting stock due to large propagule size, limited availability, seasonal dependence, low multiplication rate and rooting percentage. Therefore, attempts have been made to propagate bamboos through in vitro techniques. In vitro flowering has also been achieved successfully in some bamboo species. Classification systems proposed to date need further support, as taxonomic delineation at lower levels is still lacking sufficient resolution. Tremendous advancement in molecular markers holds the promise to address the needs of bamboo taxonomy (systematics and identification) and diversity studies. Successful application of molecular marker techniques has been achieved in several bamboo species although, more studies are required to understand the population structure and genetic diversity of bamboos in a better way. In addition, some efforts have also been made to clone important genes from bamboos and also for genetic transformation using Agrobacterium and particle bombardment methods. An overview of the recent developments made in improvement of bamboos through in vitro propagation, molecular marker technologies, cloning, and transformation and transgenics has been presented. The future potential of improvement of bamboos using modern biotechnological tools has also been discussed.

Transcriptome sequencing and de novo analysis for Ma bamboo (Dendrocalamus latiflorus Munro) using the Illumina platform

BACKGROUND

Bamboo occupies an important phylogenetic node in the grass family with remarkable sizes, woodiness and a striking life history. However, limited genetic research has focused on bamboo partially because of the lack of genomic resources. The advent of high-throughput sequencing technologies enables generation of genomic resources in a short time and at a minimal cost, and therefore provides a turning point for bamboo research. In the present study, we performed de novo transcriptome sequencing for the first time to produce a comprehensive dataset for the Ma bamboo (Dendrocalamus latiflorus Munro).

RESULTS

The Ma bamboo transcriptome was sequenced using the Illumina paired-end sequencing technology. We produced 15,138,726 reads and assembled them into 103,354 scaffolds. A total of 68,229 unigenes were identified, among which 46,087 were annotated in the NCBI non-redundant protein database and 28,165 were annotated in the Swiss-Prot database. Of these annotated unigenes, 11,921 and 10,147 unigenes were assigned to gene ontology categories and clusters of orthologous groups, respectively. We could map 45,649 unigenes onto 292 pathways using the Kyoto Encyclopedia of Genes and Genomes Pathway database. The annotated unigenes were compared against Moso bamboo, rice and millet. Unigenes that did not match any of those three sequence datasets are considered to be Ma bamboo unique. We predicted 105 unigenes encoding eight key enzymes involved in lignin biosynthesis. In addition, 621 simple sequence repeats (SSRs) were detected.

CONCLUSION

Our data provide the most comprehensive transcriptomic resource currently available for D. latiflorus Munro. Candidate genes potentially involved in growth and development were identified, and those predicted to be unique to Ma bamboo are expected to give a better insight on Ma bamboo gene diversity. Numerous SSRs characterized contributed to marker development. These data constitute a new valuable resource for genomic studies on D. latiflorus Munro and, more generally, bamboo.

Development and characterization of EST-SSR markers from Scapharca broughtonii and their transferability in Scapharca subcrenata and Tegillarca granosa

Twenty-five novel EST-derived simple sequence repeat (EST-SSR) markers were developed in the ark shell Scapharca broughtonii. Polymorphisms of these EST-SSR markers were evaluated in 48 wild individuals collected from Shidao, Shandong Province, China. A total of 202 alleles were detected at 25 loci. The numbers of alleles per locus ranged from 4 to 14, with an average of 8.08. The observed and expected heterozygosities varied from 0.2917 to 1.000 and from 0.3570 to 0.9002, respectively. After sequential Bonferroni correction for multiple tests, only one locus was found to deviate from Hardy-Weinberg equilibrium. Twenty-five EST-SSR markers showed a high rate of across-species transferability (100%) in Scapharca subcrenata and a low rate of across-genus transferability (20%) in Tegillarca granosa. These EST-SSRs will be helpful for QTL mapping, molecular breeding and investigation of population genetic diversity in ark shell S. broughtonii and other Scapharca species.

Yellow lupin (Lupinus luteus L.) transcriptome sequencing: molecular marker development and comparative studies

Background

Yellow lupin (Lupinus luteus L.) is a minor legume crop characterized by its high seed protein content. Although grown in several temperate countries, its orphan condition has limited the generation of genomic tools to aid breeding efforts to improve yield and nutritional quality. In this study, we report the construction of 454-expresed sequence tag (EST) libraries, carried out comparative studies between L. luteus and model legume species, developed a comprehensive set of EST-simple sequence repeat (SSR) markers, and validated their utility on diversity studies and transferability to related species.

Results

Two runs of 454 pyrosequencing yielded 205 Mb and 530 Mb of sequence data for L1 (young leaves, buds and flowers) and L2 (immature seeds) EST- libraries. A combined assembly (L1L2) yielded 71,655 contigs with an average contig length of 632 nucleotides. L1L2 contigs were clustered into 55,309 isotigs. 38,200 isotigs translated into proteins and 8,741 of them were full length. Around 57% of L. luteus sequences had significant similarity with at least one sequence of Medicago, Lotus, Arabidopsis, or Glycine, and 40.17% showed positive matches with all of these species. L. luteus isotigs were also screened for the presence of SSR sequences. A total of 2,572 isotigs contained at least one EST-SSR, with a frequency of one SSR per 17.75 kbp. Empirical evaluation of the EST-SSR candidate markers resulted in 222 polymorphic EST-SSRs. Two hundred and fifty four (65.7%) and 113 (30%) SSR primer pairs were able to amplify fragments from L. hispanicus and L. mutabilis DNA, respectively. Fifty polymorphic EST-SSRs were used to genotype a sample of 64 L. luteus accessions. Neighbor-joining distance analysis detected the existence of several clusters among L. luteus accessions, strongly suggesting the existence of population subdivisions. However, no clear clustering patterns followed the accession’s origin.

Conclusion

L. luteus deep transcriptome sequencing will facilitate the further development of genomic tools and lupin germplasm. Massive sequencing of cDNA libraries will continue to produce raw materials for gene discovery, identification of polymorphisms (SNPs, EST-SSRs, INDELs, etc.) for marker development, anchoring sequences for genome comparisons and putative gene candidates for QTL detection.

Development and validation of genic-SSR markers in sesame by RNA-seq

Background

Sesame (Sesamum indicum L.) is one of the most important oil crops; however, a lack of useful molecular markers hinders current genetic research. We performed transcriptome sequencing of samples from different sesame growth and developmental stages, and mining of genic-SSR markers to identify valuable markers for sesame molecular genetics research.

Results

In this study, 75 bp and 100 bp paired-end RNA-seq was used to sequence 24 cDNA libraries, and 42,566 uni-transcripts were assembled from more than 260 million filtered reads. The total length of uni-transcript sequences was 47.99 Mb, and 7,324 SSRs (SSRs ≥15 bp) and 4,440 SSRs (SSRs ≥18 bp) were identified. On average, there was one genic-SSR per 6.55 kb (SSRs ≥15 bp) or 10.81 kb (SSRs ≥18 bp). Among perfect SSRs (≥18 bp), di-nucleotide motifs (48.01%) were the most abundant, followed by tri- (20.96%), hexa- (25.37%), penta- (2.97%), tetra- (2.12%), and mono-nucleotides (0.57%). The top four motif repeats were (AG/CT)n [1,268 (34.51%)], (CA/TG)n [281 (7.65%)], (AT/AT)n [215 (5.85%)], and (GAA/TTC)n [131 (3.57%)]. A total of 2,164 SSR primer pairs were identified in the 4,440 SSR-containing sequences (≥18 bp), and 300 SSR primer pairs were randomly chosen for validation. These SSR markers were amplified and validated in 25 sesame accessions (24 cultivated accessions, one wild species). 276 (92.0%) primer pairs yielded PCR amplification products in 24 cultivars. Thirty two primer pairs (11.59%) exhibited polymorphisms. Moreover, 203 primer pairs (67.67%) yielded PCR amplicons in the wild accession and 167 (60.51%) were polymorphic between species. A UPGMA dendrogram based on genetic similarity coefficients showed that the correlation between genotype and geographical source was low and that the genetic basis of sesame in China is narrow, as previously reported. The 32 polymorphic primer pairs were validated using an F₂ mapping population; 18 primer pairs exhibited polymorphisms between the parents, and 14 genic-SSRs could be integrated into 9 main linkage groups.

Conclusions

2,164 genic-SSR markers have been developed in sesame using transcriptome sequencing. 276 of 300 validated primer pairs successfully yielded PCR amplicons in 24 cultivated sesame accessions. These markers increase current SSR marker resources and will greatly benefit genetic diversity, qualitative and quantitative trait mapping and marker-assisted selection studies in sesame.

Large-scale development of expressed sequence tag-derived simple sequence repeat markers and diversity analysis in Arachis spp

Large-scale development of expressed sequence tag simple sequence repeat (EST-SSR) markers was performed in peanut (Arachis hypogaea L.) to obtain more informative genetic markers. A total of 10,102 potential non-redundant EST sequences, including 3,445 contigs and 6,657 singletons, were generated from cDNA libraries of the gynophore, roots, leaves and seedlings. A total of 3,187 primer pairs were designed on flanking regions of SSRs, some of which allowed one and two base mismatches. Among the 3,187 markers generated, 2,540 (80%) were trinucleotide repeats, 302 (9%) were dinucleotide repeats, and 345 (11%) were tetranucleotide repeats. Pre-polymorphic analyses of 24 Arachis accessions were performed using 10% polyacrylamide gels. A total of 1,571 EST-SSR markers showing clear polymorphisms were selected for further polymorphic analysis with a Fluoro-fragment Analyzer. The 16 Arachis accessions examined included cultivated peanut varieties as well as diploid species with the A or B genome. Altogether 1,281 (81.5%) of the 1,571 markers were polymorphic among the 16 accessions, and 366 (23.3%) were polymorphic among the 12 cultivated varieties. Diversity analysis was performed and the genotypes of all 16 Arachis accessions showed similarity coefficients ranging from 0.37 to 0.97. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s11032-011-9604-8) contains supplementary material, which is available to authorized users.

Generation and analysis of expressed sequence tags from a normalized cDNA library of young leaf from Ma bamboo (Dendrocalamus latiflorus Munro)

Ma bamboo (Dendrocalamus latiflorus Munro) belongs to Dendrocalamus genus, Bambusease tribe, Bambusoideae subfamily, Poaceae family. It is a representative species of clumping bamboo, and a principal commercial species for various construction purposes using mature culms and for human consumption using young shoots. A normalized cDNA library was constructed from young leaves of Ma bamboo and 9,574 high-quality ESTs were generated, from which 5,317 unigenes including 1,502 contigs and 3,815 singletons were assembled. The unigenes were assigned into different gene ontology (GO) categories and summarized into 13 broad biologically functional groups according to similar functional characteristics or cellular roles by BLAST search against public databases. Eight hundred and ninety-one unigenes were assigned by KO identifiers and mapped to six KEGG biochemical pathways. The transcripts involved in biosynthesis of secondary metabolites such as cytochrome 450, flavonol synthase/flavanone 3-hydroxylase, and dihydroflavonol-4-reductase were well represented by 14 unigenes in the unigene set. The candidate genes involved in phytohormone metabolism, signal transduction and encoding cell wall-associated receptor kinases were also identified. Sixty-seven unigenes related to plant resistance (R) genes, including RPP genes, RGAs and RDL/RF genes, were discovered. These results will provide genome-wide knowledge about the molecular physiology of Ma bamboo young leaves and tools for advanced studies of molecular mechanism underlying leaf growth and development.

Development and characterization of microsatellite markers for the Cape gooseberry Physalis peruviana

Physalis peruviana, commonly known as Cape gooseberry, is an Andean Solanaceae fruit with high nutritional value and interesting medicinal properties. In the present study we report the development and characterization of microsatellite loci from a P. peruviana commercial Colombian genotype. We identified 932 imperfect and 201 perfect Simple Sequence Repeats (SSR) loci in untranslated regions (UTRs) and 304 imperfect and 83 perfect SSR loci in coding regions from the assembled Physalis peruviana leaf transcriptome. The UTR SSR loci were used for the development of 162 primers for amplification. The efficiency of these primers was tested via PCR in a panel of seven P. peruviana accessions including Colombia, Kenya and Ecuador ecotypes and one closely related species Physalis floridana. We obtained an amplification rate of 83% and a polymorphic rate of 22%. Here we report the first P. peruviana specific microsatellite set, a valuable tool for a wide variety of applications, including functional diversity, conservation and improvement of the species.

Microsatellites in different Potyvirus genomes: survey and analysis

Simple sequence repeats (SSRs) have been extensively used for various genetic and evolutionary studies in eukaryotic and prokaryotic organisms, while few relevant researches have been made in viruses. The Potyvirus is a fine system to study roles and evolution of SSRs in viruses. The densities, relative abundances, compositions and evolutionary inferences of SSRs in 45 different Potyvirus genomes have been analyzed in this study. Results showed that the densities and relative abundances of SSRs are similar in all those Potyvirus genomes. The number of SSRs decreases with an increase in the length of repeat unit. Dinucleotide repeats are the most abundant and followed by trinucleotide repeats, and the numbers of tetra-, penta- and hexanucleotide repeats are very small. Repeats of AC/CA, AG/GA and AAG/GAA predominate, whereas repeats of CG/GC, ATA and CAC are rare. The genome sizes of the Potyvirus species have little influence on the total number and relative abundance of SSRs. Our study suggested that the variety of SSRs may be related to the genome diversity of Potyvirus. Maybe Potyvirus and HIV genomes have the similar evolution mode and parallel evolution level.

Microsatellite isolation and marker development in carrot - genomic distribution, linkage mapping, genetic diversity analysis and marker transferability across Apiaceae

BACKGROUND

The Apiaceae family includes several vegetable and spice crop species among which carrot is the most economically important member, with ~21 million tons produced yearly worldwide. Despite its importance, molecular resources in this species are relatively underdeveloped. The availability of informative, polymorphic, and robust PCR-based markers, such as microsatellites (or SSRs), will facilitate genetics and breeding of carrot and other Apiaceae, including integration of linkage maps, tagging of phenotypic traits and assisting positional gene cloning. Thus, with the purpose of isolating carrot microsatellites, two different strategies were used; a hybridization-based library enrichment for SSRs, and bioinformatic mining of SSRs in BAC-end sequence and EST sequence databases. This work reports on the development of 300 carrot SSR markers and their characterization at various levels.

RESULTS

Evaluation of microsatellites isolated from both DNA sources in subsets of 7 carrot F2 mapping populations revealed that SSRs from the hybridization-based method were longer, had more repeat units and were more polymorphic than SSRs isolated by sequence search. Overall, 196 SSRs (65.1%) were polymorphic in at least one mapping population, and the percentage of polymophic SSRs across F2 populations ranged from 17.8 to 24.7. Polymorphic markers in one family were evaluated in the entire F2, allowing the genetic mapping of 55 SSRs (38 codominant) onto the carrot reference map. The SSR loci were distributed throughout all 9 carrot linkage groups (LGs), with 2 to 9 SSRs/LG. In addition, SSR evaluations in carrot-related taxa indicated that a significant fraction of the carrot SSRs transfer successfully across Apiaceae, with heterologous amplification success rate decreasing with the target-species evolutionary distance from carrot. SSR diversity evaluated in a collection of 65 D. carota accessions revealed a high level of polymorphism for these selected loci, with an average of 19 alleles/locus and 0.84 expected heterozygosity.

CONCLUSIONS

The addition of 55 SSRs to the carrot map, together with marker characterizations in six other mapping populations, will facilitate future comparative mapping studies and integration of carrot maps. The markers developed herein will be a valuable resource for assisting breeding, genetic, diversity, and genomic studies of carrot and other Apiaceae.

Transcriptome sequencing of lentil based on second-generation technology permits large-scale unigene assembly and SSR marker discovery

Background

Lentil (Lens culinaris Medik.) is a cool-season grain legume which provides a rich source of protein for human consumption. In terms of genomic resources, lentil is relatively underdeveloped, in comparison to other Fabaceae species, with limited available data. There is hence a significant need to enhance such resources in order to identify novel genes and alleles for molecular breeding to increase crop productivity and quality.

Results

Tissue-specific cDNA samples from six distinct lentil genotypes were sequenced using Roche 454 GS-FLX Titanium technology, generating c. 1.38 × 10⁶ expressed sequence tags (ESTs). De novo assembly generated a total of 15,354 contigs and 68,715 singletons. The complete unigene set was sequence-analysed against genome drafts of the model legume species Medicago truncatula and Arabidopsis thaliana to identify 12,639, and 7,476 unique matches, respectively. When compared to the genome of Glycine max, a total of 20,419 unique hits were observed corresponding to c. 31% of the known gene space. A total of 25,592 lentil unigenes were subsequently annoated from GenBank. Simple sequence repeat (SSR)-containing ESTs were identified from consensus sequences and a total of 2,393 primer pairs were designed. A subset of 192 EST-SSR markers was screened for validation across a panel 12 cultivated lentil genotypes and one wild relative species. A total of 166 primer pairs obtained successful amplification, of which 47.5% detected genetic polymorphism.

Conclusions

A substantial collection of ESTs has been developed from sequence analysis of lentil genotypes using second-generation technology, permitting unigene definition across a broad range of functional categories. As well as providing resources for functional genomics studies, the unigene set has permitted significant enhancement of the number of publicly-available molecular genetic markers as tools for improvement of this species.

Development of soybean aphid genomic SSR markers using next generation sequencing

Simple sequence repeats (SSRs) or microsatellites are very useful molecular markers, owing to their locus-specific codominant and multiallelic nature, high abundance in the genome, and high rates of transferability across species. The soybean aphid (Aphis glycines Matsumura) has become the most damaging insect pest of soybean (Glycine max (L.) Merr.) in North America, since it was first found in the Midwest of the United States in 2000. Biotypes of the soybean aphid capable of colonizing newly developed aphid-resistant soybean cultivars have been recently discovered. Genetic resources, including molecular markers, to study soybean aphids are severely lacking. Recently developed next generation sequencing platforms offer opportunities for high-throughput and inexpensive genome sequencing and rapid marker development. The objectives of this study were (i) to develop and characterize genomic SSR markers from soybean aphid genomic sequences generated by next generation sequencing technology and (ii) to evaluate the utility of the SSRs for genetic diversity or relationship analyses. In total 128 SSR primer pairs were designed from sequences generated by Illumina GAII from a reduced representation library of A. glycines. Nearly 94% (120) of the primer pairs amplified SSR alleles of expected size and 24 SSR loci were polymorphic among three aphid samples from three populations. The polymorphic SSRs were successfully used to differentiate among 24 soybean aphids from Ohio and South Dakota. Sequencing of PCR products of two SSR markers from four aphid samples revealed that the allelic polymorphism was due to variation in the SSR repeats among the aphids. These markers should be particularly useful for genetic differentiation among aphids collected from soybean fields at different localities and regions. These SSR markers provide the soybean aphid research community with the first set of PCR-based codominant markers developed from the genomic sequences of A. glycines.

SSR and EST-SSR-based genetic linkage map of cassava (Manihot esculenta Crantz)

Simple sequence repeat (SSR) markers provide a powerful tool for genetic linkage map construction that can be applied for identification of quantitative trait loci (QTL). In this study, a total of 640 new SSR markers were developed from an enriched genomic DNA library of the cassava variety 'Huay Bong 60' and 1,500 novel expressed sequence tag-simple sequence repeat (EST-SSR) loci were developed from the Genbank database. To construct a genetic linkage map of cassava, a 100 F(1) line mapping population was developed from the cross Huay Bong 60 by 'Hanatee'. Polymorphism screening between the parental lines revealed that 199 SSRs and 168 EST-SSRs were identified as novel polymorphic markers. Combining with previously developed SSRs, we report a linkage map consisted of 510 markers encompassing 1,420.3 cM, distributed on 23 linkage groups with a mean distance between markers of 4.54 cM. Comparison analysis of the SSR order on the cassava linkage map and the cassava genome sequences allowed us to locate 284 scaffolds on the genetic map. Although the number of linkage groups reported here revealed that this F(1) genetic linkage map is not yet a saturated map, it encompassed around 88% of the cassava genome indicating that the map was almost complete. Therefore, sufficient markers now exist to encompass most of the genomes and efficiently map traits in cassava.

Development of genic-SSR markers by deep transcriptome sequencing in pigeonpea [Cajanus cajan (L.) Millspaugh]

BACKGROUND

Pigeonpea [Cajanus cajan (L.) Millspaugh], one of the most important food legumes of semi-arid tropical and subtropical regions, has limited genomic resources, particularly expressed sequence based (genic) markers. We report a comprehensive set of validated genic simple sequence repeat (SSR) markers using deep transcriptome sequencing, and its application in genetic diversity analysis and mapping.

RESULTS

In this study, 43,324 transcriptome shotgun assembly unigene contigs were assembled from 1.696 million 454 GS-FLX sequence reads of separate pooled cDNA libraries prepared from leaf, root, stem and immature seed of two pigeonpea varieties, Asha and UPAS 120. A total of 3,771 genic-SSR loci, excluding homopolymeric and compound repeats, were identified; of which 2,877 PCR primer pairs were designed for marker development. Dinucleotide was the most common repeat motif with a frequency of 60.41%, followed by tri- (34.52%), hexa- (2.62%), tetra- (1.67%) and pentanucleotide (0.76%) repeat motifs. Primers were synthesized and tested for 772 of these loci with repeat lengths of ≥ 18 bp. Of these, 550 markers were validated for consistent amplification in eight diverse pigeonpea varieties; 71 were found to be polymorphic on agarose gel electrophoresis. Genetic diversity analysis was done on 22 pigeonpea varieties and eight wild species using 20 highly polymorphic genic-SSR markers. The number of alleles at these loci ranged from 4-10 and the polymorphism information content values ranged from 0.46 to 0.72. Neighbor-joining dendrogram showed distinct separation of the different groups of pigeonpea cultivars and wild species. Deep transcriptome sequencing of the two parental lines helped in silico identification of polymorphic genic-SSR loci to facilitate the rapid development of an intra-species reference genetic map, a subset of which was validated for expected allelic segregation in the reference mapping population.

CONCLUSION

We developed 550 validated genic-SSR markers in pigeonpea using deep transcriptome sequencing. From these, 20 highly polymorphic markers were used to evaluate the genetic relationship among species of the genus Cajanus. A comprehensive set of genic-SSR markers was developed as an important genomic resource for diversity analysis and genetic mapping in pigeonpea.

Mining of simple sequence repeats in the Genome of Gentianaceae

Simple sequence repeats (SSRs) or short tandem repeats are short repeat motifs that show high level of length polymorphism due to insertion or deletion mutations of one or more repeat types. Here, we present the detection and abundance of microsatellites or SSRs in nucleotide sequences of Gentianaceae family. A total of 545 SSRs were mined in 4698 nucleotide sequences downloaded from the National Center for Biotechnology Information (NCBI). Among the SSR sequences, the frequency of repeat type was about 429 -mono repeats, 99 -di repeats, 15 -tri repeats, and 2 --hexa repeats. Mononucleotide repeats were found to be abundant repeat types, about 78%, followed by dinucleotide repeats (18.16%) among the SSR sequences. An attempt was made to design primer pairs for 545 identified SSRs but these were found only for 169 sequences.

Genetic structure and ecogeographical adaptation in wild barley (Hordeum chilense Roemer et Schultes) as revealed by microsatellite markers

BACKGROUND

Multi-allelic microsatellite markers have become the markers of choice for the determination of genetic structure in plants. Synteny across cereals has allowed the cross-species and cross-genera transferability of SSR markers, which constitute a valuable and cost-effective tool for the genetic analysis and marker-assisted introgression of wild related species. Hordeum chilense is one of the wild relatives with a high potential for cereal breeding, due to its high crossability (both interspecies and intergenera) and polymorphism for adaptation traits. In order to analyze the genetic structure and ecogeographical adaptation of this wild species, it is necessary to increase the number of polymorphic markers currently available for the species. In this work, the possibility of using syntenic wheat SSRs as a new source of markers for this purpose has been explored.

RESULTS

From the 98 wheat EST-SSR markers tested for transferability and polymorphism in the wild barley genome, 53 primer pairs (54.0%) gave cross-species transferability and 20 primer pairs (20.4%) showed polymorphism. The latter were used for further analysis in the H. chilense germplasm. The H. chilense-Triticum aestivum addition lines were used to test the chromosomal location of the new polymorphic microsatellite markers. The genetic structure and diversity was investigated in a collection of 94 H. chilense accessions, using a set of 49 SSR markers distributed across the seven chromosomes. Microsatellite markers showed a total of 351 alleles over all loci. The number of alleles per locus ranged from two to 27, with a mean of 7.2 alleles per locus and a mean Polymorphic Information Content (PIC) of 0.5.

CONCLUSIONS

According to the results, the germplasm can be divided into two groups, with morphological and ecophysiological characteristics being key determinants of the population structure. Geographic and ecological structuring was also revealed in the analyzed germplasm. A significant correlation between geographical and genetic distance was detected in the Central Chilean region for the first time in the species. In addition, significant ecological influence in genetic distance has been detected for one of the population structure groups (group II) in the Central Chilean region. Finally, the association of the SSR markers with ecogeographical variables was investigated and one marker was found significantly associated with precipitation. These findings have a potential application in cereal breeding.

EST-derived SSR markers used as anchor loci for the construction of a consensus linkage map in ryegrass (Lolium spp.)

BACKGROUND

Genetic markers and linkage mapping are basic prerequisites for marker-assisted selection and map-based cloning. In the case of the key grassland species Lolium spp., numerous mapping populations have been developed and characterised for various traits. Although some genetic linkage maps of these populations have been aligned with each other using publicly available DNA markers, the number of common markers among genetic maps is still low, limiting the ability to compare candidate gene and QTL locations across germplasm.

RESULTS

A set of 204 expressed sequence tag (EST)-derived simple sequence repeat (SSR) markers has been assigned to map positions using eight different ryegrass mapping populations. Marker properties of a subset of 64 EST-SSRs were assessed in six to eight individuals of each mapping population and revealed 83% of the markers to be polymorphic in at least one population and an average number of alleles of 4.88. EST-SSR markers polymorphic in multiple populations served as anchor markers and allowed the construction of the first comprehensive consensus map for ryegrass. The integrated map was complemented with 97 SSRs from previously published linkage maps and finally contained 284 EST-derived and genomic SSR markers. The total map length was 742 centiMorgan (cM), ranging for individual chromosomes from 70 cM of linkage group (LG) 6 to 171 cM of LG 2.

CONCLUSIONS

The consensus linkage map for ryegrass based on eight mapping populations and constructed using a large set of publicly available Lolium EST-SSRs mapped for the first time together with previously mapped SSR markers will allow for consolidating existing mapping and QTL information in ryegrass. Map and markers presented here will prove to be an asset in the development for both molecular breeding of ryegrass as well as comparative genetics and genomics within grass species.

Development of a EST dataset and characterization of EST-SSRs in a traditional Chinese medicinal plant, Epimedium sagittatum (Sieb. Et Zucc.) Maxim

Background

Epimedium sagittatum (Sieb. Et Zucc.) Maxim, a traditional Chinese medicinal plant species, has been used extensively as genuine medicinal materials. Certain Epimedium species are endangered due to commercial overexploition, while sustainable application studies, conservation genetics, systematics, and marker-assisted selection (MAS) of Epimedium is less-studied due to the lack of molecular markers. Here, we report a set of expressed sequence tags (ESTs) and simple sequence repeats (SSRs) identified in these ESTs for E. sagittatum.

Results

cDNAs of E. sagittatum are sequenced using 454 GS-FLX pyrosequencing technology. The raw reads are cleaned and assembled into a total of 76,459 consensus sequences comprising of 17,231 contigs and 59,228 singlets. About 38.5% (29,466) of the consensus sequences significantly match to the non-redundant protein database (E-value < 1e-10), 22,295 of which are further annotated using Gene Ontology (GO) terms. A total of 2,810 EST-SSRs is identified from the Epimedium EST dataset. Trinucleotide SSR is the dominant repeat type (55.2%) followed by dinucleotide (30.4%), tetranuleotide (7.3%), hexanucleotide (4.9%), and pentanucleotide (2.2%) SSR. The dominant repeat motif is AAG/CTT (23.6%) followed by AG/CT (19.3%), ACC/GGT (11.1%), AT/AT (7.5%), and AAC/GTT (5.9%). Thirty-two SSR-ESTs are randomly selected and primer pairs are synthesized for testing the transferability across 52 Epimedium species. Eighteen primer pairs (85.7%) could be successfully transferred to Epimedium species and sixteen of those show high genetic diversity with 0.35 of observed heterozygosity (Ho) and 0.65 of expected heterozygosity (He) and high number of alleles per locus (11.9).

Conclusion

A large EST dataset with a total of 76,459 consensus sequences is generated, aiming to provide sequence information for deciphering secondary metabolism, especially for flavonoid pathway in Epimedium. A total of 2,810 EST-SSRs is identified from EST dataset and ~1580 EST-SSR markers are transferable. E. sagittatum EST-SSR transferability to the major Epimedium germplasm is up to 85.7%. Therefore, this EST dataset and EST-SSRs will be a powerful resource for further studies such as taxonomy, molecular breeding, genetics, genomics, and secondary metabolism in Epimedium species.

Microsatellites in Brassica unigenes: relative abundance, marker design, and use in comparative physical mapping and genome analysis

Microsatellites present in the transcribed regions of the genome have the potential to reveal functional diversity. Unigene sequence databases are the sources of such genic microsatellites with unique flanking sequences and genomic locations even in complex polyploids. The present study was designed to assay the unigenes of Brassica napus and B. rapa for various microsatellite repeats, and to design markers and use them in comparative genome analysis and study of evolution. The average frequency of microsatellites in Brassica unigenes was one in every 7.25 kb of sequence, as compared with one in every 8.57 kb of sequence in Arabidopsis thaliana . Trinucleotide motifs coding for serine and the dinucleotide motif GA were most abundant. We designed 2374 and 347 unigene-based microsatellite (UGMS) markers including 541 and 58 class I types in B. napus and B. rapa, respectively, and evaluated their use across diverse species and genera. Most of these markers (93.3%) gave successful amplification of target microsatellite motifs, which was confirmed by sequencing. Interspecific polymorphism between B. napus and B. rapa detected in silico for the UGMS markers was 4.16 times higher in 5′ untranslated regions than in coding sequences. Physical anchoring of Brassica UGMS markers on the A. thaliana genome indicated their significance in studying the evolutionary history of A. thaliana genomic duplications in relation to speciation. Comparative physical mapping identified 85% of Brassica unigenes as single copy and gave clues for the presence of conserved primordial gene order. Complex chromosomal rearrangements such as inversions, tandem and segmental duplications, and insertions/deletions were evident between A. thaliana and B. rapa genomes. The results obtained have encouraging implications for the use of UGMS markers in comparative genome analysis and for understanding evolutionary complexities in the family Brassicaceae.

Mining of expressed sequence tag libraries of cacao for microsatellite markers using five computational tools

Expressed sequence tags (ESTs) provide researchers with a quick and inexpensive route for discovering new genes, data on gene expression and regulation, and also provide genic markers that help in constructing genome maps. Cacao is an important perennial crop of humid tropics. Cacao EST sequences, as available in the public domain, were downloaded and made into contigs. Microsatellites were located in these ESTs and contigs using five softwares (MISA, TRA, TROLL, SSRIT and SSR primer). MISA gave maximum coverage of SSRs in cacao ESTs and contigs, although TRA was able to detect higher order (5-mer) repeats. The frequency of SSRs was one per 26.9 kb in the known set of ESTs. One-third of the repeats in EST-contigs were found to be trimeric. A few rare repeats like 21-mer repeat were also located. A/T repeats were most abundant among the mononucleotide repeats and the AG/GA/TC/CT type was the most frequent among dimerics. Flanking primers were designed using Primer3 program and verified experimentally for PCR amplification. The results of the study are made available freely online database (http://riju.byethost31.com/cocoa/). Seven primer pairs amplified genomic DNA isolated from leaves were used to screen a representative set of 12 accessions of cacao.