Next Article in Journal
A Frameshift Mutation in KIT is Associated with  White Spotting in the Arabian Camel
Previous Article in Journal
Alternative Splicing Profile and Sex-Preferential Gene Expression in the Female and Male Pacific Abalone Haliotis discus hannai
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Development of Gene‐Based SSR Markers in Winged  Bean (Psophocarpus tetragonolobus (L.) DC.) for  Diversity Assessment

1
Biotechnology Research Centre, School of Biosciences, Faculty of Science, University of Nottingham Malaysia Campus, Jalan Broga, 43500 Semenyih, Selangor Darul Ehsan, Malaysia
2
Crops For the Future, Jalan Broga, 43500 Semenyih, Selangor Darul Ehsan, Malaysia
3
Deep Seq, Faculty of Medicine and Health Sciences, Queen’s Medical Centre, University of Nottingham, Nottingham NG7 2UH, UK
4
Department of Export Agriculture, Faculty of Agricultural Sciences, Sabaragamuwa University of Sri Lanka, Belihuloya 70140, Sri Lanka
5
School of Biosciences, Faculty of Science, University of Nottingham Sutton Bonington Campus, Sutton Bonington, Leicestershire LE12 5RD, UK
*
Author to whom correspondence should be addressed.
Submission received: 20 December 2016 / Accepted: 3 March 2017 / Published: 9 March 2017
(This article belongs to the Section Plant Genetics and Genomics)

Abstract

:
Winged bean (Psophocarpus tetragonolobus) is an herbaceous multipurpose legume grown in hot and humid countries as a pulse, vegetable (leaves and pods), or root tuber crop depending on local consumption preferences. In addition to its different nutrient-rich edible parts which could contribute to food and nutritional security, it is an efficient nitrogen fixer as a component of sustainable agricultural systems. Generating genetic resources and improved lines would help to accelerate the breeding improvement of this crop, as the lack of improved cultivars adapted to specific environments has been one of the limitations preventing wider use. A transcriptomic de novo assembly was constructed from four tissues: leaf, root, pod, and reproductive tissues from Malaysian accessions, comprising of 198,554 contigs with a N50 of 1462 bp. Of these, 138,958 (70.0%) could be annotated. Among 9682 genic simple sequence repeat (SSR) motifs identified (excluding monomer repeats), trinucleotide-repeats were the most abundant (4855), followed by di-nucleotide (4500) repeats. A total of 18 SSR markers targeting di- and tri-nucleotide repeats have been validated as polymorphic markers based on an initial assessment of nine genotypes originated from five countries. A cluster analysis revealed provisional clusters among this limited, yet diverse selection of germplasm. The developed assembly and validated genic SSRs in this study provide a foundation for a better understanding of the plant breeding system for the genetic improvement of winged bean.

1. Introduction

Winged bean (Psophocarpus tetragonolobus (L.) DC.) (2n = 2x = 18) is a tropical perennial vine species, classified in the family of Fabaceae and subfamily of Papilionoideae, that is cultivated mainly at a subsistence scale in hot and humid countries across India, Southeast Asia, and the Western Pacific islands, with a presence in a number of African countries as well [1,2,3,4,5]. It is grown for its green pods, tuberous roots, and mature seeds, all of which have received attention for their nutritional content in the past, as comprehensively described in ‘The Winged Bean—A high-protein crop for the tropics’ from the National Academy of Science in 1981 [1]. Initial interest was drawn to high crude protein levels in seeds, which are comparable to soybean [6,7,8]. Its vining nature and nitrogen fixation activity have seen it used as a cover crop and also incorporated into rotation or intercropping systems [9,10,11]. As such, winged bean could be a good candidate for diversifying diets to improve nutritional security, based on complex and more sustainable agricultural systems [12]. Despite its potential, winged bean has received limited research investment for developing molecular tools that can support breeding programmes, until recently. Recent reports include the development of inter-Simple Sequence Repeats (iSSRs) and Randomly Amplified Polymorphic DNA (RAPD) markers for genetic diversity and for clonal fidelity analyses and two small transcriptomic assemblies derived from a mix of leaf, bud, and shoot of Sri Lankan accessions and leaf tissue from a Nigerian accession, respectively [13,14,15,16,17]. Given that winged bean is believed to be largely self-pollinated, heterozygosity would be expected to be low, although a formal assessment is needed and the species does produce large flowers, suggesting a contribution from insect pollination, as recorded by Erskine [18]. Thus, molecular breeding will facilitate utilisation of genetic resources in winged bean breeding, especially among accessions, through combining beneficial traits. Molecular markers that are tightly linked to important agronomic traits are a precondition for undertaking molecular breeding in plants. The genetic basis of traits in winged bean remains largely unexplored, and to date there has not been any genetic linkage map reported for this crop, although controlled crosses have been reported [19,20,21,22].
In this study, we generated RNA-seq data from four tissues (leaf, root, reproductive tissues, and pod) of six locally grown accessions, followed by the identification of SSR-containing sequences and validation of a subset of genic SSR markers. To our knowledge, this is the first application of within-species genic SSR markers in winged bean accessions. The data will help to begin the development of comprehensive genetic information and tools to facilitate future breeding programmes, as well as allow the levels of natural inbreeding to be determined, to allow appropriate breeding schemes to be devised. The transcriptome will allow us to gain a better understanding of the phylogenetic relationships between winged bean and other leguminous and model plants.

2. Material and Methods

2.1. Plant Material, RNA Extraction, Complementary DNA (cDNA) Library Construction, and Sequencing

A total of six locally grown winged bean accessions (two derived from Malaysian Agricultural Research and Development Institute (MARDI) and four from local planters) were grown from August to December 2012 at Lady Bird Farm, Broga, Semenyih, Malaysia (Latitude: 2.9394 N; Longitude: 101.8971 E; Altitude: 45m asl). RNA was extracted separately from leaf, root, pod, and reproductive tissue (comprising of bud and flower) by pooling the respective tissues from all the six accessions. Extraction was performed from different tissue groups separately using TRIzol Reagent (Thermo Fisher Scientific, Waltham, MA, USA) followed by another round of purification using RNeasy MinElute Cleanup Kit (Qiagen, Hilden, Germany) before library preparation.
Total RNA was measured using the Qubit RNA BR assay kit (Thermo Fisher Scientific, Waltham, MA, USA). A total of 5 µg of RNA was used for enrichment of mRNA using NEBNext Poly(A) mRNA Magnetic Isolation Module (New England Biolabs, Beverly, MA, USA). RNA fragmentation was done using NEBNext Magnesium RNA Fragmentation Module (New England Biolabs, Beverly, MA, USA). Illumina stranded whole transcriptome sequencing libraries were prepared through a dUTP approach using NEBNext Ultra Directional RNA Library Prep Kit for Illumina (New Engladn Biolabs, Beverly, MA, USA). Libraries were gel purified using 2% E-Gel SizeSelect (Thermo Fisher Scientific, Waltham, MA, USA) and quality control was performed using bioanalyser HS kit (Agilent biotechnologies, Palo Alto, CA, USA). Quantification was done using qPCR (quantitative polymerase chain reaction) (Kapa Biosystems, Woburn, MA, USA). Equimolar amounts of barcoded libraries were mixed and subjected to 250 bp paired-end run using MiSeq V2 chemistry (Illumina, San Diego, CA, USA) on Illumina MiSeq sequencing platform according to manufacturer’s instruction.
For simple sequence repeat (SSR) marker development, a total of nine plants from five accessions—each one representing a geographical origin—from the International Institute of Tropical Agriculture (IITA) genebank and MARDI, were used (as listed in Table 1 below). Two individuals were used per accession, which were collected after a cycle of single seed descent purification from January to June 2013 at the Lady Bird Farm, except for the Malaysian line.

2.2. De Novo Transcriptome Assembly and Microsatellite Identification

Adaptors and low quality reads (below Q20) were trimmed using Scythe and Sickle [23], respectively, using default settings. Trimmed reads from all tissues were pooled to assemble a combined de novo assembly with Trinity version 2.2.0 pipeline using the strand specificity option. Trinity assembled transcripts were annotated with Trinotate software suite version 1.1 [24], with a blast e-value threshold of 1 × 10−5 from NCBI-BLAST [25], HMMER/PFAM [26], SignalP [27], EMBL eggNOG [28], and Gene Ontology (GO) [29] databases. The data is deposited in NCBI Sequence Read Archive (BioProject ID PRJNA374598) under the accession number of SRP099538; SRR5252646 (root), SRR5252647 (reproductive tissue), SRR5252648 (pod), and SRR5252649 (leaf).
This was followed by the identification of microsatellites using MIcroSAtellite (MISA) Perl script program, based on a minimum number of repeats of six for di-, five for tri-, tetra-, penta- and hexa-nucleotide repeat motif (monomer repeats were excluded), whilst the maximum number of bases interrupting two SSRs in a compound microsatellite was 100 [30].

2.3. Microsatellite Markers Development and Scoring

Primer pairs were designed from sequences harbouring a minimum of 18-bases long microsatellites (i.e., minimum 9 and 6 repetitions for di- and tri-nucleotide motifs, respectively) (Table S1). Primer3 [31] and PrimerQuest (Integrated DNA Technologies, Coralville, Iowa, USA) were used for oligo design, with the latter used whenever the first was not able to design with standard parameters for the given sequence. Where possible, two pairs of primers were designed for a single target region, so as to have an alternative if the first pair failed to amplify the target region. DNA was extracted from the leaf of nine genotypes (Table 1) using a modified cetyltrimethylammonium bromide (CTAB) method [32]. In addition, an RNase digestion step was added and 1 volume of isopropanol was used instead of 2/3 volume. Primer screening and optimisation was carried out in a three primer system for fluorescent labelling [33] using an equimolar mixture of genotypes. Each 20 µL of PCR reaction consisted of 1× Buffer S, 200 µM dNTPs, 0.02 µM forward primer, 0.18 µM M13-dye labelled primer, 0.2 µM reverse primer, 20 ng of DNA, and 1 U of Taq DNA Polymerase (Vivantis, Subang Jaya, Selangor, Malaysia). The PCR was programmed for 3 min of initial denaturation at 94 °C, followed by 35 cycles of 1 min at 94 °C, 1 min at 60 °C, and 1 min at 72 °C, with a 10 min final elongation. Initial polymorphism evaluation of primers was performed across all genotypes using a long capillary fragment analyser (Advance Analytical Technologies – AATI, Ankeny, IA, USA) with default parameters, except for using 4 µL of samples, and 4 kV of separation voltage for 180 min per run. Data was then analysed with Advance Analytical PROSize 2.0 v1.3.1.1 (Advance Analytical Technologies – AATI, Ankeny, IA, USA) to identify potential SSR markers based on presence of polymorphic amplicons across genotypes.
PCR products of single genotypes from potentially polymorphic SSR markers were then separated using an ABI Genetic Analyser ABI3730XL using Peak Scanner v2 for scoring (Applied Biosystems, Foster City, CA, USA). Validated markers were subsequently characterised using Power Marker v3.5 [34] for major allelic frequency, alleles per marker, heterozygosity, and polymorphic information content (PIC).

2.4. Cluster Analysis

A hierarchical cluster analysis was performed with the Dice (also known as Nei and Li) similarity coefficient and unweighted pair-group method with arithmetic mean (UPGMA) algorithm in Genstat 18th Edition [35].

3. Results and Discussion

3.1. Transcriptome Assembly and In Silico Identification of Microsatellites

A total of four libraries generated 12.77 million reads of cleaned 250 bp paired ends (Table 2). The de novo assembly derived from all tissues produced a total of 198,554 contigs with an average size of 798 bp and an N50 of 1462 bp.
Out of 198,554 contigs, 138,958 (70.0%) could be annotated. Among them, 75,308 (54.2%), 69,172 (49.8%), 6499 (4.7%), 60,040 (43.2%), and 70,069 (50.4%) were found in NCBI-BLAST, HMMER/PFAM, SignalP, EMBL eggNOG, and GO databases, respectively, with no significant homology found from tmHMM on the prediction of transmembrane helices (Table S2). Figure 1 illustrates the abundance of transcripts classified based on gene ontology.
In this study, a total of 9682 putative SSR repeat motifs were identified from 8793 SSR containing sequences, which came from 4.4% of the total contig number in this assembly (Table S3). On average, there was one SSR locus for every 16.4 kbp of de novo assembly. After excluding mononucleotide motifs, trinucleotide repeats were the most abundant type (50.1%) (summarised in Table 3). This is consistent with Vatanparast et al.’s study [16], although hexamer motifs were not evaluated in this study. The most frequent dimer motifs were AG/GA/CT/TC type, followed by AT/TA, whereas for trimeric repeats, AAG/AGA/GAA/CTT/TCT/TTC were the most abundant (Figure 2 and Table 4). Both observations on the most common di- and tri-nucleotide repeat motif are in agreement with the winged bean transcriptome from Vatanparast et al. as well as with the soybean, medicago, and lotus Expressed Sequence Tag (EST)-SSR summarised by Jayashree et al. [16,36].

3.2 Development of SSR Markers and Cluster Analysis

A total of 56 (targeting 42 dimer-repeat regions) and 78 (targeting 53 trinucleotide SSR) primer pairs were designed. Subsequently, 20 dinucleotide SSR primers and 26 trinucleotide SSR primers gave good amplification products at the expected size. After polymorphism evaluation using all genotypes in this study, 18 validated SSR markers (8 for di-nucleotide and 10 for tri-nucleotide repeated motifs; Table S1) were scored and are summarised in Table 5. The low validation rate of polymorphic markers is likely to be partly due to the limited number of accessions screened, and should increase with more accessions covering a broader range of geographical origins. Residual heterozygosity could still be observed within each accession (shaded values in Table 5), even where a cycle of line purification in a controlled environment has been carried out, indicating that further cycles are needed to obtain homozygous lines, in particular for Tpt10, Tpt53, and M3. This data, along with the winged bean large flower size, also suggest that such a purification process may need to be carried out under an insect-proof enclosed environment. Using these markers, an average of 2.5 and 2.4 alleles per locus for di- and tri-nucleotide SSRs, respectively, was observed (Table 6). Individual PIC values varied from 0.16 to 0.67, which is comparable to recent legume studies in pigeonpea [37], mungbean [38], and common bean [39], although lower than in cowpea [40] and bambara groundnut [41].
The cluster analysis from the SSR scores (Figure 3) showed a few clusters with the accessions originating from Papua New Guinea closely related to the Sri Lankan accession, but sharing the least similarity with the Malaysian and Indonesian materials, comparatively. To our knowledge, the genetic relationship between germplasm from Bangladesh and Malaysia are here investigated for the first time with molecular markers, and place the Bangladesh origin closer to the Sri Lankan and Papua New Guinean germplasm. Although the number of accessions used in this study is limited, they cover a reasonable range of germplasm from different origins.

4. Conclusions

A set of validated functional winged bean genic-SSR markers is reported here for the first time, to our best knowledge. The reported residual heterozygosity across screened genotypes has suggested that further investigation needs to be carried out on the rate of natural outcrossing in winged bean, in order to understand how genetic materials should be maintained, improved, and introduced into breeding programmes. The cluster analysis provides an initial insight into the potential for these markers to be used on a larger number of winged bean accessions, to carry out a more comprehensive diversity analysis with the evaluation of germplasm from genebanks and from commonly cultivated lines. Finally, this set of 18 microsatellite markers could also be used to contribute to genetic linkage maps in winged bean, with the integration of single nucleotide polymorphisms (SNPs) markers for higher density. Such a map would be the first backbone for linkage analysis and the genetic dissection of traits with agronomic importance in this legume.

Supplementary Materials

The following are available online at www.mdpi.com/2073-4425/8/3/100/s1, Table S1: the sequences of forward and reverse SSR-primers used in this study; Table S2: Functional annotation of assembled transcripts; Table S3: Identified microsatellites.

Acknowledgments

This work was financially supported by Kouk Foundation Berhad, Crops for the Future (CFF), and the University of Nottingham Malaysia Campus (UNMC) through the CFF-UNMC Doctoral Training Partnership scheme. The authors would like to thank Victoria Wright from Deep Seq. Seeds were contributed by the International Institute of Tropical Agriculture (IITA) and Malaysian Agricultural Research and Development Institute (MARDI). This paper is dedicated to the memory of Quin Nee Wong, who passed away in a tragic car accident on 12 October 2014.

Author Contributions

Q.N.W. and A.S.T. carried out planting and sampling work and, together with W.K.H., SSR marker development and validation. A.S.T. and W.K.H. wrote the manuscript. S.M. performed library preparation for RNA-seq and M.B. performed bioinformatics analysis. A.K. contributed Sri Lanka genetic materials. F.M. helped to draft the manuscript and contributed to experimental design. S.Ms. conceived of the study, participated in its design and coordination, and revised the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. NAS. The Winged Bean: High-Protein Crop for the Humid Tropics, 2nd ed.; National Academy Press: Washington, DC, USA, 1981. [Google Scholar]
  2. Klu, G.Y.P. Induced mutations for accelerated domestication-a case study of winged bean. West Afr. J. Appl. Ecol. 2000, 1, 47–52. [Google Scholar]
  3. Khan, T.N. Papua New Guinea: A centre of genetic diversity in winged bean (Psophocarpus tetragonologus (L.) Dc.). Euphytica 1976, 25, 693–705. [Google Scholar] [CrossRef]
  4. Harder, D.K. Chromosome Counts in Psophocarpus. Kew Bull. 1992, 47, 529–534. [Google Scholar] [CrossRef]
  5. Harder, D.K.; Smartt, J. Further evidence on the origin of the cultivated winged bean, Psophocarpus tetragonolobus (L.) DC. Econ. Bot. 1992, 46, 187–191. [Google Scholar] [CrossRef]
  6. Prakash, D.; Misra, P.N.; Misra, P.S. Amino acid profile of winged bean (Psophocarpus tetragonolobus (L.) DC.): A rich source of vegetable protein. Plant Foods Hum. Nutr. 1987, 37, 261–264. [Google Scholar] [CrossRef] [PubMed]
  7. Kadam, S.S.; Salunkhe, D.K. Winged bean in human nutrition. Crit. Rev. Food Sci. Nutr. 1984, 21, 1–40. [Google Scholar] [CrossRef] [PubMed]
  8. Okezie, B.O.; Martin, F.W. Chemical composition of dry seeds and fresh leaves of winged bean varieties grown in the U.S. and Puerto Rico. J. Food Sci. 1980, 45, 1045–1051. [Google Scholar] [CrossRef]
  9. Anugroho, F.; Kitou, M.; Kinjo, K.; Kobashigawa, N. Growth and nutrient accumulation of winged bean and velvet bean as cover crops in a subtropical region. Plant Prod. Sci. 2010, 13, 360–366. [Google Scholar] [CrossRef]
  10. Banerjee, A.; Bagchi, D.K.; Si, L.K. Studies on the potential of winged bean as a multipurpose legume cover crop in tropical regions. Exp. Agric. 2008, 20, 297–301. [Google Scholar] [CrossRef]
  11. Hikam, S.; MacKown, C.T.; Poneleit, C.G.; Hildebrand, D.F. Growth and N accumulation in maize and winged bean as affected by N level and intercropping. Ann. Bot. 1991, 68, 17–22. [Google Scholar] [CrossRef]
  12. FAO. Coping with Climate Change—The Roles of Genetic Resources for Food and Agriculture; FAO: Rome, Italy, 2015. [Google Scholar]
  13. Mohanty, C.S.; Verma, S.; Singh, V.; Khan, S.; Gaur, P.; Gupta, P.; Nizar, M.A.; Dikshit, N.; Pattanayak, R.; Shukla, A.; et al. Characterization of winged bean (Psophocarpus tetragonolobus (L.) DC.) based on molecular, chemical and physiological parameters. Am. J. Mol. Biol. 2013, 3, 187–197. [Google Scholar] [CrossRef]
  14. Koshy, E.P.P.; Alex, B.K.K.; John, P. Clonal fidelity studies on regenerants of Psophocarpus tetragonolobus (L.) DC. using RAPD markers. Bioscan 2013, 8, 763–766. [Google Scholar]
  15. Chen, D.; Yi, X.; Yang, H.; Zhou, H.; Yu, Y.; Tian, Y.; Lu, X. Genetic diversity evaluation of winged bean (Psophocarpus tetragonolobus (L.) DC.) using inter-simple sequence repeat (ISSR). Genet. Resour. Crop Evol. 2015, 62, 823–828. [Google Scholar] [CrossRef]
  16. Vatanparast, M.; Shetty, P.; Chopra, R.; Doyle, J.J.; Sathyanarayana, N.; Egan, A.N. Transcriptome sequencing and marker development in winged bean (Psophocarpus tetragonolobus; Leguminosae). Sci. Rep. 2016, 6, 29070. [Google Scholar] [CrossRef] [PubMed]
  17. Chapman, M.A. Transcriptome sequencing and marker development for four underutilized legumes. Appl. Plant Sci. 2015, 3, 1400111. [Google Scholar] [CrossRef] [PubMed]
  18. Erskine, W. Measurements of the cross-pollination of winged bean in Papua New Guinea. SABRAO J. Breed. Genet. 1980, 12, 11–14. [Google Scholar]
  19. Eagleton, G.E. Evaluation of Genetic Resources in the Winged Bean (Psophocarpus tetragonolobus (L.) DC.) and Their Utilisation in the Development of Cultivars for Higher Latitudes. Ph.D. Dissertation, University of Western Australia, Perth, Australia, 1983. [Google Scholar]
  20. De Silva, H.N.; Omran, A. Diallel analysis of yield and yield components of winged bean (Psophocarpus tetragonolobus (L.) D.C.). J. Agric. Sci. 1986, 106, 485–490. [Google Scholar] [CrossRef]
  21. Erskine, W.; Khan, T.N. Inheritance of pigmentation and pod shape in winged bean. Euphytica 1977, 26, 829–831. [Google Scholar] [CrossRef]
  22. Erskine, B.Y.W. Heritability and combining ability of vegetative and phenological characters of winged beans (Psophocarpus tetragonolobus (L.) DC.). J. Agric. Sci. 1981, 96, 503–508. [Google Scholar] [CrossRef]
  23. Joshi, N.A.; Fass, J.N. Sickle: A Sliding-Window, Adaptive, Quality-Based Trimming Tool for FastQ Files, version 1.33; Software; 2011. Available online: https://github.com/najoshi/sickle (accessed on 13 March 2013).
  24. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.; et al. Full-length transcriptome assembly from RNA-seq data without a reference genome. Nat. Biotechnol. 2011, 29, 644–652. [Google Scholar] [CrossRef] [PubMed]
  25. Camacho, C.; Coulouris, G.; Avagyan, V.; Ma, N.; Papadopoulos, J.; Bealer, K.; Madden, T.L. BLAST+: Architecture and applications. BMC Bioinform. 2009, 10, 421. [Google Scholar] [CrossRef] [PubMed]
  26. Finn, R.D.; Bateman, A.; Clements, J.; Coggill, P.; Eberhardt, R.Y.; Eddy, S.R.; Heger, A.; Hetherington, K.; Holm, L.; Mistry, J.; et al. Pfam: The protein families database. Nucleic Acids Res. 2014, 42, D222–D230. [Google Scholar] [CrossRef] [PubMed]
  27. Petersen, T.N.; Brunak, S.; von Heijne, G.; Nielsen, H. SignalP 4.0: Discriminating signal peptides from transmembrane regions. Nat. Methods 2011, 8, 785–786. [Google Scholar] [CrossRef] [PubMed]
  28. Powell, S.; Forslund, K.; Szklarczyk, D.; Trachana, K.; Roth, A.; Huerta-Cepas, J.; Gabaldón, T.; Rattei, T.; Creevey, C.; Kuhn, M.; et al. eggNOG v4.0: Nested orthology inference across 3686 organisms. Nucleic Acids Res. 2014, 42, D231–D239. [Google Scholar] [CrossRef] [PubMed]
  29. Gene Ontology Consortium. The Gene Ontology (GO) database and informatics resource. Nucleic Acids Res. 2004, 32 (Suppl. 1), D258–D261. [Google Scholar]
  30. Thiel, T.; Michalek, W.; Varshney, R.K.; Graner, A. Exploiting EST databases for the development and characterization of gene-derived SSR-markers in barley (Hordeum vulgare L.). Theor. Appl. Genet. 2003, 106, 411–422. [Google Scholar] [CrossRef] [PubMed]
  31. Koressaar, T.; Remm, M. Enhancements and modifications of primer design program Primer3. Bioinformatics 2007, 23, 1289–1291. [Google Scholar] [CrossRef] [PubMed]
  32. Doyle, J. DNA Protocols for Plants. In Molecular Techniques in Taxonomy; Hewitt, G.M., Johnston, A.W.B., Young, J.P.W., Eds.; Springer: Berlin/Heidelberg, Germany, 1991; pp. 283–293. [Google Scholar]
  33. Schuelke, M. An economic method for the fluorescent labeling of PCR fragments. Nat. Biotechnol. 2000, 18, 233–234. [Google Scholar] [CrossRef] [PubMed]
  34. Liu, K.; Muse, S.V. PowerMarker: An integrated analysis environment for genetic marker analysis. Bioinformatics 2005, 21, 2128–2129. [Google Scholar] [CrossRef] [PubMed]
  35. VSN International. GenStat for Windows, 18th ed.; VSN International: Hemel Hempstead, UK, 2015. [Google Scholar]
  36. Jayashree, B.; Punna, R.; Prasad, P.; Bantte, K.; Hash, C.T.; Chandra, S.; Hoisington, D.A.; Varshney, R.K. A database of simple sequence repeats from cereal and legume expressed sequence tags mined in silico: Survey and evaluation. In Silico Biol. 2006, 6, 607–620. [Google Scholar] [PubMed]
  37. Dutta, S.; Kumawat, G.; Singh, B.P.; Gupta, D.K.; Singh, S.; Dogra, V.; Gaikwad, K.; Sharma, T.R.; Raje, R.S.; Bandhopadhya, T.K.; et al. Development of genic-SSR markers by deep transcriptome sequencing in pigeonpea [Cajanus cajan (L.) Millspaugh]. BMC Plant Biol. 2011, 11, 17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Gupta, S.K.; Bansal, R.; Gopalakrishna, T. Development and characterization of genic SSR markers for mungbean (Vigna radiata (L.) Wilczek). Euphytica 2014, 195, 245–258. [Google Scholar] [CrossRef]
  39. Blair, M.W.; Hurtado, N.; Chavarro, C.M.; Muñoz-Torres, M.C.; Giraldo, M.C.; Pedraza, F.; Tomkins, J.; Wing, R.; Varshney, R.; Graner, A.; et al. Gene-based SSR markers for common bean (Phaseolus vulgaris L.) derived from root and leaf tissue ESTs: An integration of the BMc series. BMC Plant Biol. 2011, 11, 50. [Google Scholar] [CrossRef] [PubMed]
  40. Gupta, S.K.; Gopalakrishna, T. Development of unigene-derived SSR markers in cowpea (Vigna unguiculata) and their transferability to other Vigna species. Genome 2010, 53, 508–523. [Google Scholar] [PubMed]
  41. Molosiwa, O.O.; Aliyu, S.; Stadler, F.; Mayes, K.; Massawe, F.; Kilian, A.; Mayes, S. SSR marker development, genetic diversity and population structure analysis of Bambara groundnut [Vigna subterranean (L.) Verdc.] landraces. Genet. Resour. Crop Evol. 2015, 62, 1225–1243. [Google Scholar] [CrossRef]
Figure 1. Distribution of first level gene ontology classification of the de novo assembly.
Figure 1. Distribution of first level gene ontology classification of the de novo assembly.
Genes 08 00100 g001
Figure 2. The number distribution of different microsatellite motif types identified. SSR, simple sequence repeat.
Figure 2. The number distribution of different microsatellite motif types identified. SSR, simple sequence repeat.
Genes 08 00100 g002
Figure 3. A dendrogram of the genetic relationship between genotypes from Papua New Guinea (Tpt10), Sri Lanka (SLS319), Bangladesh (Tpt53), Indonesia (Tpt17), and Malaysia (M3).
Figure 3. A dendrogram of the genetic relationship between genotypes from Papua New Guinea (Tpt10), Sri Lanka (SLS319), Bangladesh (Tpt53), Indonesia (Tpt17), and Malaysia (M3).
Genes 08 00100 g003
Table 1. Winged bean accessions (two individuals per origin, except Malaysian line) used and their origins.
Table 1. Winged bean accessions (two individuals per origin, except Malaysian line) used and their origins.
IndividualsOrigin
Tpt53-9-8
Tpt53-9-10
Bangladesh
Tpt17-6-3
Tpt17-6-8
Indonesia
M3-3Malaysia
Tpt10-7-5
Tpt10-7-7
Papua New Guinea
SLS319-10-3
SLS319-10-4
Sri Lanka
Table 2. Summary statistics for the de novo assembled transcriptome.
Table 2. Summary statistics for the de novo assembled transcriptome.
TissueLeafPodReproductive TissueRoot
Number of raw read/base (bp)3,150,356/1,544,004,8223,973,092/1,868,456,6803,544,968/1,719,303,6323,873,893/1,859,527,511
Numbers of trimmed read/base (bp)3,113,502/1,438,258,1803,157,832/1,301,766,1133,199,527/1,431,024,1413,303,324/1,461,908,151
Number of contigs/base (bp)198,554/158,382,439
Average contig size (bp)798
N501462
Table 3. In silico identification of microsatellites from the mixed tissue assembly. SSR, simple sequence repeat.
Table 3. In silico identification of microsatellites from the mixed tissue assembly. SSR, simple sequence repeat.
Total number of sequences examined198,554
Total size of examined sequences (bp)158,382,439
Total number of identified SSRs9682
Number of SSR containing sequences8793
Number of sequences containing more than one SSR780
Number of SSRs present in compound formation352
Number of dimer-repeat4500
Number of trimer-repeat4855
Number of tetramer-repeat279
Number of pentamer-repeat48
Table 4. Frequency distribution of di- and tri-nucleotide motif repeat in this de novo assembly.
Table 4. Frequency distribution of di- and tri-nucleotide motif repeat in this de novo assembly.
Number of Repeat MotifTotal%
Di-nucleotide5678910>10
AC/GT/CA/TG-2561386337121051611.5
AG/CT/GA/TC-995543330391434189288264.0
AT/TA-407201116711310768106323.6
CG/GC-3810000390.9
Total-1696 (37.7%)883 (19.6%)560 (12.4%)541 (12.0%)553 (12.3%)267 (5.9%)4500
Tri-nucleotide
AAC/ACA/CAA/GTT/TGT/TTG27913150170004779.8
AAG/AGA/GAA/CTT/TCT/TTC61240535111000137928.4
AAT/ATA/TAA/TTA/TAT/ATT3051451121500057711.9
ACC/CAC/CCA/GGT/GTG/TGG3076255100004348.9
ACG/CGA/GAC/CGT/GTC/TCG90661270001753.6
ACT/CTA/TAC/AGT/TAG/GTA36833000501.0
AGC/CAG/GCA/TGC/CTG/GCT2711053890004238.7
AGG/GGA/GAG/TCC/CTC/CCT24711575110004489.2
ATC/CAT/TCA/GAT/ATG/TGA3118324340004529.3
CCG/CGC/GCC/GGC/GCG/CGG2471305580004409.1
Total2705 (55.7%)1250 (25.7%)775 (16.0%)125 (2.6%)0004855
Table 5. Scores of 18 SSR markers from nine winged bean individuals.
Table 5. Scores of 18 SSR markers from nine winged bean individuals.
MarkerPapua New GuineaIndonesiaBangladeshSri LankaMalaysia
Tpt10-7-5Tpt10-7-7Tpt17-6-3Tpt17-6-8Tpt53-9-8Tpt53-9-10SLS319-10-3SLS319-10-4M3-3
P27.2205199/205199205205205205205205
P43.2199199195195197199199199195
Pt1.1335335339339339335/339335335339
Pt10226/228228226226228228228228228
Pt14358358352352350350358358354
Pt24219217/219217217219219219219217
Pt7.2426/432426426426428428426426426/428
WB17198198198198198194/198198198198
Pt53315309315315309/315309/315312312315
Pt58255/261255/261261261261261261261261
Pt65.1273273267267267267267267267/273
Pt67.1293293296296293293296296293/296
Pt68.1226226229229226223/226223223/226226/235
Pt76.1203203203203209209209209209
Pt78.1306/309306/309306306309309306306309
Pt85.1276/279276/279276276276276276276279
Pt93.1266266272272266/272272266266276
Pt99.2189/195189/195195195189189189189195
Table 6. A summary of data analysis of 18 SSR markers. PIC, polymorphic information content.
Table 6. A summary of data analysis of 18 SSR markers. PIC, polymorphic information content.
MarkerSSR MotifMajor Allele FrequencyNo. of AllelesHeterozygosityPIC
P27.2TA0.8320.110.24
P43.2TA0.56300.49
Pt1.1CT0.520.110.38
Pt10TC0.7220.110.32
Pt14TG0.44400.64
Pt24GT0.6120.110.36
Pt7.2TC0.6730.220.4
WB17GA0.9420.110.1
Average dimer SSR markers0.662.50.10.37
Pt53CGC0.5630.220.53
Pt58TAG0.8920.220.18
Pt65.1CAG0.7220.110.32
Pt67.1AGA0.520.110.38
Pt68.1AAC0.540.330.59
Pt76.1CGC0.56200.37
Pt78.1AAC0.5620.220.37
Pt85.1GCG0.7820.220.29
Pt93.1TGT0.530.110.5
Pt99.2TTC0.5620.220.37
Average trimer SSR marker0.612.40.180.39

Share and Cite

MDPI and ACS Style

Wong, Q.N.; Tanzi, A.S.; Ho, W.K.; Malla, S.; Blythe, M.; Karunaratne, A.; Massawe, F.; Mayes, S. Development of Gene‐Based SSR Markers in Winged  Bean (Psophocarpus tetragonolobus (L.) DC.) for  Diversity Assessment. Genes 2017, 8, 100. https://0-doi-org.brum.beds.ac.uk/10.3390/genes8030100

AMA Style

Wong QN, Tanzi AS, Ho WK, Malla S, Blythe M, Karunaratne A, Massawe F, Mayes S. Development of Gene‐Based SSR Markers in Winged  Bean (Psophocarpus tetragonolobus (L.) DC.) for  Diversity Assessment. Genes. 2017; 8(3):100. https://0-doi-org.brum.beds.ac.uk/10.3390/genes8030100

Chicago/Turabian Style

Wong, Quin Nee, Alberto Stefano Tanzi, Wai Kuan Ho, Sunir Malla, Martin Blythe, Asha Karunaratne, Festo Massawe, and Sean Mayes. 2017. "Development of Gene‐Based SSR Markers in Winged  Bean (Psophocarpus tetragonolobus (L.) DC.) for  Diversity Assessment" Genes 8, no. 3: 100. https://0-doi-org.brum.beds.ac.uk/10.3390/genes8030100

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop