Development of Gene-Based SSR Markers in Rice Bean (Vigna umbellata L.) Based on Transcriptome Data

Rice bean (Vigna umbellata (Thunb.) Ohwi & Ohashi) is a warm season annual legume mainly grown in East Asia. Only scarce genomic resources are currently available for this legume crop species and no simple sequence repeat (SSR) markers have been specifically developed for rice bean yet. In this study, approximately 26 million high quality cDNA sequence reads were obtained from rice bean using Illumina paired-end sequencing technology and assembled into 71,929 unigenes with an average length of 986 bp. Of these unigenes, 38,840 (33.2%) showed significant similarity to proteins in the NCBI non-redundant protein and nucleotide sequence databases. Furthermore, 30,170 (76.3%) could be classified into gene ontology categories, 25,451 (64.4%) into Swiss-Prot categories and 21,982 (55.6%) into KOG database categories (E-value < 1.0E-5). A total of 9,301 (23.5%) were mapped onto 118 pathways using the Kyoto Encyclopedia of Genes and Genome (KEGG) pathway database. A total of 3,011 genic SSRs were identified as potential molecular markers. AG/CT (30.3%), AAG/CTT (8.1%) and AGAA/TTCT (20.0%) are the three main repeat motifs. A total of 300 SSR loci were randomly selected for validation by using PCR amplification. Of these loci, 23 primer pairs were polymorphic among 32 rice bean accessions. A UPGMA dendrogram revealed three major clusters among 32 rice bean accessions. The large number of SSR-containing sequences and genic SSRs in this study will be valuable for the construction of high-resolution genetic linkage maps, association or comparative mapping and genetic analyses of various Vigna species.     

Development of genomic simple sequence repeat markers in opium poppy by next-generation sequencing

Opium poppy (Papaver somniferum L.) is an important pharmaceutical crop with very few genetic marker resources. To expand these resources, we sequenced genomic DNA using pyrosequencing technology and examined the DNA sequences for simple sequence repeats (SSRs). A total of 1,244,412 sequence reads were obtained covering 474 Mb. Approximately half of the reads (52 %) were assembled into 166,724 contigs representing 105 Mb of the opium poppy genome. A total of 23,283 non-redundant SSRs were identified in 18,944 contigs (11.3 % of total contigs). Trinucleotide and tetranucleotide repeats were the most abundant SSR repeats, accounting for 49.0 and 27.9 % of all SSRs, respectively. The AAG/TTC repeat was the most abundant trinucleotide repeat, representing 19.7 % of trinucleotide repeats. Other SSR repeat types were AT-rich. A total of 23,126 primer pairs (98.7 % of total SSRs) were designed to amplify SSRs. Fifty-three genomic SSR markers were tested in 37 opium poppy accessions and seven Papaver species for determination of polymorphism and transferability. Intraspecific polymorphism information content (PIC) values of the genomic SSR markers were intermediate, with an average 0.17, while the interspecific average PIC value was slightly higher, 0.19. All markers showed at least 88 % transferability among related species. This study increases sequence coverage of the opium poppy genome by sevenfold and the number of opium poppy-specific SSR markers by sixfold. This is the first report of the development of genomic SSR markers in opium poppy, and the genomic SSR markers developed in this study will be useful in diversity, identification, mapping and breeding studies in opium poppy.     

Identification,characterization, validation and cross-species amplification of genic-SSRs in Indian Mustard (<Emphasis Type="Italic">Brassica juncea</Emphasis>)

Brassica juncea is an economically important oilseed crop worldwide. It has limited genomic resources at present. We generated 47,962,057 expressed sequence reads which were assembled into 45,280 unigenes. A total of 4108 SSR loci (≥10 bp) were identified in these unigenes. Trinucleotide was the most frequent repeat unit (59.91 %) followed by di- (38.66 %), tetra - (0.71 %), hexa - (0.49 %) and pentanucleotide repeats (0.24 %). Primers were designed for 2863 SSR loci among which 460 were selected for primer synthesis. A total of 339 loci amplified successfully of which 134 (39.5 %) exhibited polymorphism among six B. juncea genotypes with PIC values ranging from 0.18 to 0.81. Further, 25 polymorphic SSRs were used for analysis of genetic variability in 25 genotypes of Brassicas and their wild relatives. Two to five alleles with PIC values 0.22–0.66 were detected at these loci. The dendrogram grouped the genotypes according to their known pedigree/systematic position.     

Novel and Stress Relevant EST Derived SSR Markers Developed and Validated in Peanut

With the aim to increase the number of functional markers in resource poor crop like cultivated peanut (Arachis hypogaea), large numbers of available expressed sequence tags (ESTs) in the public databases, were employed for the development of novel EST derived simple sequence repeat (SSR) markers. From 16424 unigenes, 2784 (16.95%) SSRs containing unigenes having 3373 SSR motifs were identified. Of these, 2027 (72.81%) sequences were annotated and 4124 gene ontology terms were assigned. Among different SSR motif-classes, tri-nucleotide repeats (33.86%) were the most abundant followed by di-nucleotide repeats (27.51%) while AG/CT (20.7%) and AAG/CTT (13.25%) were the most abundant repeat-motifs. A total of 2456 EST-SSR novel primer pairs were designed, of which 366 unigenes having relevance to various stresses and other functions, were PCR validated using a set of 11 diverse peanut genotypes. Of these, 340 (92.62%) primer pairs yielded clear and scorable PCR products and 39 (10.66%) primer pairs exhibited polymorphisms. Overall, the number of alleles per marker ranged from 1-12 with an average of 3.77 and the PIC ranged from 0.028 to 0.375 with an average of 0.325. The identified EST-SSRs not only enriched the existing molecular markers kitty, but would also facilitate the targeted research in marker-trait association for various stresses, inter-specific studies and genetic diversity analysis in peanut.     

Development of genomic simple sequence repeat markers in faba bean by next-generation sequencing

Faba bean (Vicia faba L.) is an important food legume crop with a huge genome. Development of genetic markers for faba bean is important to study diversity and for molecular breeding. In this study, we used Next Generation Sequencing (NGS) technology for the development of genomic simple sequence repeat (SSR) markers. A total of 14,027,500 sequence reads were obtained comprising 4,208 Mb. From these reads, 56,063 contigs were assembled (16,367 Mb) and 2138 SSRs were identified. Mono and dinucleotides were the most abundant, accounting for 57.5 % and 20.9 % of all SSR repeats, respectively. A total of 430 primer pairs were designed from contigs larger than 350 nucleotides and 50 primers pairs were tested for validation of SSR locus amplification. Nearly all (96 %) of the markers were found to produce clear amplicons and to be reproducible. Thirty-nine SSR markers were then applied to 46 faba bean accessions from worldwide origins, resulting in 161 alleles with 87.5 % polymorphism, and an average of 4.1 alleles per marker. Gene diversity (GD) of the markers ranged from 0 to 0.48 with an average of 0.27. Testing of the markers showed that they were useful in determining genetic relationships and population structure in faba bean accessions.     

Developing and characterising Ricinus communis SSR markers by data mining of whole-genome sequences

Ricinus communis is a versatile industrial oil crop that is cultivated worldwide. Genetic improvement and marker-assisted breeding of castor bean have been slowed owing to the lack of abundant and efficient molecular markers. As co-dominant markers, simple sequence repeats (SSRs) are useful for genetic evaluation and molecular breeding. The recently released whole-genome sequence of castor bean provides useful genomic resources for developing markers on a genome-wide scale. In the present study, the distribution and frequency of microsatellites in the castor bean genome were characterised and numerous SSR markers were developed using genomic data mining. In total, 18,647 SSR loci at a density of one SSR per 18.89 Kb in the castor bean genome sequence (representing approximately 352.27 Mb) were identified. Dinucleotide repeats were the most frequently observed microsatellites, although the AAT repeat motif was also prevalent. Using six cultivars as screening samples, 670 polymorphic SSR markers from 1,435 primer pairs (46.7 %) were developed. Trinucleotide motif loci contained a higher proportion of polymorphisms (48.5 %) than dinucleotide motif loci (39.2 %). The polymorphism level in the SSR loci was positively correlated with the increasing number of repeat units in the microsatellites. The phylogenetic relationship among 32 varieties was evaluated using the developed SSR markers. Cultivars developed at the same institute clustered together, suggesting that these cultivars have a narrow genetic background. The large number of SSR markers developed in this study will be useful for genetic mapping and for breeding improved castor-oil plants. These markers will also facilitate genetic and genomic studies of Euphorbiaceae.     

Development of EST-SSR markers in castor bean (Ricinus communis L.) and their utilization for genetic purity testing of hybrids

Expressed sequence tag (EST) databases offer opportunity for the rapid development of simple sequence repeat (SSR) markers in crops. Sequence assembly and clustering of 57?895 ESTs of castor bean resulted in the identification of 10?960 unigenes (6459 singletons and 4501 contigs) having 7429 SSRs. On an average, the unigenes contained 1 SSR for every 1.23?kb of unigene sequence. The identified SSRs mostly consisted of dinucleotide (62.4%) and trinucleotide (33.5%) repeats. The AG class was the most common among the dinucleotide motifs (68.9%), whereas the AAG class (25.9%) was predominant among the trinucleotide motifs. A total of 611 primer pairs were designed for the SSRs, having repeat length more than or equal to 20 nucleotides, of which a set of 130 markers were tested and 92 of these yielding robust amplicons were analyzed for their utility in genetic purity assessment of castor bean hybrids. Nine markers were able to detect polymorphism between the parental lines of nine commercial castor bean hybrids (DCH-32, DCH-177, DCH-519, GCH-2, GCH-4, GCH-5, GCH-6, GCH-7, and RHC-1), and their utility in genetic purity testing was demonstrated. These novel EST-SSR markers would be a valuable addition to the growing molecular marker resources that could be used in genetic improvement programmes of castor bean.     

Genome-wide single nucleotide polymorphism discovery and validation in adzuki bean

Adzuki bean, also known as red bean (Vigna angularis), with 2n = 22 chromosomes, is an important legume crop in East Asian countries, including China, Japan, and Korea. For single nucleotide polymorphism (SNP) discovery, we used Vigna accessions, V. angularis IT213134 and its wild relative V. nakashimae IT178530, because of the lack of DNA sequence polymorphism in the cultivated species. Short read sequences of IT213134 and IT178530 of approximately 37 billion and 35 billion bp were produced using the Illumina HiSeq 2000 system to a sequencing depth of 61.5× and 57.7×, respectively. After de novo assembly was carried out with trimmed HiSeq reads from IT213134, 98,441 contigs of various sizes were produced with N50 of 13,755 bp. Using Burrows–Wheeler Aligner software, trimmed short reads of V. nakashimae IT178530 were successfully mapped to IT213134 contigs. All sequence variations at the whole-genome level were examined between the two Vigna species. Of the 1,565,699 SNPs, 59.4 % were transitions and 40.6 % were transversions. A total of 213,758 SNPs, consisting of 122,327 non-synonymous and 91,431 synonymous SNPs, were identified in coding sequences. For SNP validation, 96 SNPs in the genic region were chosen from among IT213134 contigs longer than 10 kb. Of these 96 SNPs, 88 were confirmed by Sanger sequencing of 10 adzuki bean genotypes from various geographic origins as well as IT213134 and its wild relative IT178530. These genome-wide SNP markers will enrich the existing Vigna resources and, specifically, could be of value for constructing a genetic map and evaluating the genetic diversity of adzuki bean.     

Discovery and experimental analysis of microsatellites in an oil woody plant Camellia chekiangoleosa

Oil camellia trees are important woody plants for the production of high-quality cooking oil. On the contrary to their economic importance, their genetic and genomic resources are very limited, which greatly hamper the genetic studies on oil camellia trees. Microsatellites or simple sequence repeats (SSRs) have great value in many aspects of genetic analyses due to their high polymorphism and codominant inheritance. In this study, we report the large-scale development and characterization of SSR markers derived from genomic sequences of Camellia chekiangoleosa by high-throughput pyrosequencing technology. A total of 1,091,393 genomic shotgun reads were generated using Roche 454 FLX sequencer, the average read length was 319 bp, and the total sequence throughput was 347.9 Mb. These sequences were assembled into 35,315 contigs with total length of 14.8 Mb and the N50 contig size of 770 bp. By analyzing with microsatellite (MISA), a total of 5,844 perfect microsatellites were detected from the assembled sequences. Among them, tetranucleotide repeats were found to be the most frequent microsatellites in the genome of C. chekiangoleosa, and all the dominant repeat motifs for different types of SSRs were detected to be rich in A/T. Experimental analysis with 900 SSR primer pairs revealed that 66 % of them succeeded in PCR amplification. Further investigation with 345 SSR primer pairs showed that a relatively high percentage of primers amplified polymorphic loci (31.9 %). Experimental data also revealed that, overall, long microsatellite repeats (>20 bp) were more variable than the short ones (<20 bp) in the genome of oil camellia tree.     

Development of genomic simple sequence repeat markers for linseed using next-generation sequencing technology

Linseed (Linum usitatissimum L.) is regarded as a cash crop of tomorrow because of the presence of nutraceutically important ??-linolenic acid (ALA) and lignan. However, only limited breeding progress has been made in this crop, mainly due to the lack of sufficient genetic and genomic resources. Among these, simple sequence repeats (SSR) are useful DNA markers for diversity analysis, genetic mapping and tagging traits because of their co-dominant and highly polymorphic nature. In order to develop SSR markers for linseed, we used three microsatellite isolation methods, viz., PCR Isolation of Microsatellite Arrays (PIMA), 5??-anchored PCR method, and Fast Isolation by AFLP of Sequences COntaining repeats (FIASCO). The amplified products from these methods were pooled and sequenced using the 454 GS-FLX platform. A total of 36,332 reads were obtained, which assembled into 2,183 contigs and 2,509 singlets. The contigs and the singlets contained 1,842 microsatellite motifs, with dinucleotide motifs as the most abundant repeat type (54%) followed by trinucleotide motifs (44%). Based on this, 290 SSR markers were designed, 52 of which were evaluated using a panel of 27 diverse linseed genotypes. Among the three enrichment methods, the 5??-anchored PCR method was most efficient for isolation of microsatellites, while FIASCO was most efficient for developing SSR markers. We show the utility of next-generation sequencing technology for efficiently discovering a large number of microsatellite markers in non-model plants.