Genetic diversity and relationships in mulberry (genus Morus) as revealed by RAPD and ISSR marker assays
© Awasthi et al; licensee BioMed Central Ltd. 2004
Received: 11 August 2003
Accepted: 10 January 2004
Published: 10 January 2004
The genus Morus, known as mulberry, is a dioecious and cross-pollinating plant that is the sole food for the domesticated silkworm, Bombyx mori. Traditional methods using morphological traits for classification are largely unsuccessful in establishing the diversity and relationships among different mulberry species because of environmental influence on traits of interest. As a more robust alternative, PCR based marker assays including RAPD and ISSR were employed to study the genetic diversity and interrelationships among twelve domesticated and three wild mulberry species.
RAPD analysis using 19 random primers generated 128 discrete markers ranging from 500–3000 bp in size. One-hundred-nineteen of these were polymorphic (92%), with an average of 6.26 markers per primer. Among these were a few putative species-specific amplification products which could be useful for germplasm classification and introgression studies. The ISSR analysis employed six anchored primers, 4 of which generated 93 polymorphic markers with an average of 23.25 markers per primer. Cluster analysis of RAPD and ISSR data using the WINBOOT package to calculate the Dice coefficient resulted into two clusters, one comprising polyploid wild species and the other with domesticated (mostly diploid) species.
These results suggest that RAPD and ISSR markers are useful for mulberry genetic diversity analysis and germplasm characterization, and that putative species-specific markers may be obtained which can be converted to SCARs after further studies.
Mulberry (genus Morus) is an economically important plant used for sericulture, as it is the sole food plant for the domesticated silkworm, Bombyx mori. The genus Morus, which is widely distributed in Asia, Europe, North and South America, and Africa, is cultivated extensively in East, Central and South Asia for silk production. A few species of mulberry are also valued for their edible fruit (M. alba, M. indica and M. laevigata), timber (M. laevigata and M. serrata). Whereas it has been widely believed that mulberry species originated on the low slopes of the Himalayas bordering China and India, the study of Hou suggests a multicentered origin . Since the classification of the genus Morus is mainly based on morphological characteristics, considerable differences exist among systematists as to the number of species that exist in this genus [2–7]. So far, more than 150 species of mulberry have been cited in the Index Kewensis, but a majority of them have been treated either as synonyms or as varieties rather than species, and some have been transferred to allied genera. A study carried out by Koidzumi in 1917  recognised 24 species and one variety under the genus Morus based on the style length in female flowers and the nature of the stigma in male flowers. In contrast, more than 60 years later by analyzing the electrophoretic patterns of seven enzymes and sap proteins in 131 varieties of three mulberry species, M. bombycis, M. alba, and M. latifolia Hirano categorised them into seven varietal groups, and affinities among them . Because of environmental influence, phenotypic traits in many cases fail to serve as unambiguous markers for systematics and diversity analysis . Moreover, most of the putative mulberry species are dioecious and can cross-pollinate among themselves to produce fertile hybrids, suggesting that they have relatively close genetic relationships. Such a high degree of cross-species reproductive success is not encountered often in nature, and has thus created considerable doubt with regard to the species status of mulberry.
Molecular markers successfully developed during the last two decades have largely overcome the problems that are associated with phenotype-based classification. Initially, isozymes [10–12] and Restriction Fragment Length Polymorphisms (RFLPs) [13–17] served as reliable markers for genetic analyses in plants. But PCR based techniques developed in recent years such as Random Amplified Polymorphic DNA (RAPDs) [18, 19], Inter Simple Sequence Repeats (ISSR) , Amplified Fragment Length Polymorphisms (AFLPs) , and Simple Sequence Repeats (SSRs) , also called microsatellites, provide DNA markers that are dispersed throughout plant genomes  and are easier to reproduce and analyse. High levels of polymorphism and their co-dominant nature have made SSRs ideal markers for studying genetic diversity in plants [24–26]. However, the time and cost of identifying SSR motifs and designing primers for regions flanking SSRs have restricted the widespread use of microsatellites in plants [27, 28]. ISSR markers, which show dominant inheritance, use SSR repeat-anchored primers and are being used as an alternate tool in diversity studies. ISSR markers are useful in detecting genetic polymorphisms among accessions by generating a large number of markers that target multiple microsatellite loci distributed across the genome. Further, they are simpler to use than the SSR technique as prior knowledge of the target sequences flanking the repeat regions is not required [20, 29–31].
So far only a few attempts have been made to characterise the genetic diversity in mulberry by using molecular markers. These include AFLP based marker analysis , RAPD and DAMD profiles , ISSR based analysis , and genetic polymorphism estimation in mulberry hybrids using RAPDs . In the current study we report the use of RAPD and ISSR markers for assessing the genetic diversity and relationships among 12 cultivated and 3 wild mulberry species.
Morus species used in the present study.
Species (Var. Acc. No)
Ploidy and chromosome number
Country of origin
M. alba Linn. (Var. S-1 Acc. No. ME-0065)
2n = 2x = 28
Style absent or indistinct, stigma papillose. Leaf ovate or wide ovate, lobed and unlobed, surface smooth without wrinkles, apex acute or acuminate, margin usually blunt serrate or dentate, few serrulate or denticulate, base indented or shallowly cordate
M. tiliaefolia Makino.* (Acc No. ME-0126)
2n = 6x = 84
Style absent or indistinct, stigma densely pubescent. Leaf tomentose on the under surface.
M. ihou (SER) (Acc No. ME-0107)
2n = 2x = 28
Style indistinct, stigma papillose. Leaf entire, comparatively large with surface slightly pubescent or smooth, glossy and shrivelled. Leaf margin serrate, base lobate.
M. rotundiloba Koidz (Acc No. ME-0095)
2n = 2x = 28
Style long, stigma lancecolate, minutely papillose. Leaf margin dentato-serrate and denticels are nearly equal. Leaf is usually trilobate and apex of the lateral lobes is rotundate.
M. bombycis Koidz. (Acc No. ME-0028)
2n = 3x = 42
Female flower with distinct long style, stigma papillose or pilose. Leaf mostly dentato-serrate, crenato-dentate rarely argutidentato-serrate margin. Bud scales glabrous.
M. rubra Linn. (Acc No. ME-0110)
2n = 2x = 28
Style absent or indistinct, stigma papillose, leaf ovate, margin crenate, apex acuminate, base cordate. Branches slightly drooping.
M. nigra Linn. (Acc No. ME-0008)
2n = 2x = 28
Style absent, stigma pubiscent. Leaf generally ovate or cordate, apex caudate, margin serrate having erose serrae, base truncate or shallowly cordate.
M. multicaulis perr. (Acc No. ME-0006)
2n = 2x = 28
Style absent or indistinct, stigma papillose with protuberances inside. Leaf cordate generally unlobed, surface has wrinkled appearance, margin serrate, apex sharp obtuse, base cordate.
M. indica Linn. (Acc No. MI-0068)
2n = 2x = 28
Style long, stigma hairy connate below. Leaf ovate margin sharply serrate often deeply lobed, scaberulous, apex caudate to acuminate. Male sepals elliptic and hairy. Female spikes short ovoid, sepals obovate-outer keeled and inner flat.
M. latifolia Poir. (Acc No. ME-0066)
2n = 2x = 28
Style absent or indistinct, stigma papilose. Leaf broad and thick, margin serrate with serrae being acute, apex short acute, base deeply cordate.
M. sinensis Hort. (Var. MR-2 Acc No. MI-0025)
2n = 2x = 28
Style short or indistinct, stigma short with blunt hair. Leaf narrow ovate smooth, entire, margin serrate, apex acuminate, base truncate.
M. laevigata Wall. * (Acc. No. MI-0164)
2n = 3x = 42
Style very short, stigma papillose inside. Flower in spikes of 4–16 cm length. Leaf ovate or wide ovate, comparatively larger, no hairs on ventral part, margin serrulate, apex long cuspidate, base shallow cordate to truncate.
M. australis Poir (Acc No. ME-0001)
2n = 2x = 28
Style long, stigma pilose, shoot errect. Under surface of leaf light green and slightly hairy, leaf ovate, ovato-oblong or ovato-elliptical, apex caudate, margin serrate with short prickles.
M. serrata Roxb. * (Acc. No. ME-0081)
2n = 4x = 56
Style very short or absent, stigma densely pubescent. Leaf rough usually with hairs, margin typically serrate, apex cuspidate or caudate, base cordate.
M. cathayana Hemsi. (Acc No. ME-0018)
2n = 2x = 28
Style absent, stigma pubescent. Leaf cordate and ventral side covered with thick hair, margin sparsely serrate, base cordate or truncate.
List of the ISSR primers used in the present study.
(GT)8 (A/G) Y
Distance matrix values based on RAPD data between 15 mulberry species.
Putative species specific RAPD markers.
Band size (bp)
Distance matrix values based on ISSR data between 15 mulberry species.
Mulberry is a perennial, heterogeneous outbreeding tree, the leaves of which are the exclusive food of the silk secreting insect, B. mori. Outbreeding in cultivated mulberry species is common and inter-species hybridization is often observed. Because of phenotypic plasticity, the occurrence of interspecific hybridization, mutation, and the absence of an unambiguous set of criteria for designating a true species, there is confusion in the systematic classification of mulberry. The present work evaluates the genetic diversity and relationships among fifteen mulberry species using RAPD and ISSR markers.
The wide variation in genetic distance among the different mulberry species revealed by both RAPD and ISSR techniques reflected a high level of polymorphism at the DNA level. Earlier studies by Sharma et al. using the AFLP technique  also showed a large genetic variation among different Morus genotypes. Such a high level of polymorphism, especially in cultivated mulberry, reflects the outcrossing nature of the species. Cluster analysis of RAPD and ISSR data using UPGMA revealed that the three wild species, namely, M. laevigata, M. serrata, and M. tiliaefolia, are genetically distant from the domesticated species studied here. The distinctness of the wild species as revealed by both RAPD and ISSR in the present study and that of Sharma et al.  can be attributed to their geographical isolation, which is in strong contrast to the outbreeding and high heterozygosity which have accompanied the long history of cultivation of domesticated mulberry species.
The ISSR profiles of mulberry generated by CA/TG repeat anchored primers showed that these repeats are abundant in the Morus genome. Vijayan and Chatterjee  also observed amplification of AC rich repeat based ISSR primers. However, this type of repeat is less abundant in plant genomes compared to other kinds of repeats , and the abundance of other classes of microsatellite repeats in the mulberry genome awaits to be demonstrated. The RAPD analysis revealed a close relationship between M. cathayana and M. rotundiloba, whereas these two species were found to be distinct in ISSR analysis. Such variation between RAPD and ISSR may be due to the fact that the PCR amplified profiles in the two marker assays originated from different repetitive and non-repetitive regions of the genomes, and the possibility that many co-migrating bands may be non-homologous, producing a background noise that could influence the results [37, 38]. Similar to the RAPD analysis, the ISSR results showed no correlation to ploidy status or to the number of amplified products. For example, in the triploid species, M. laevigata and M. bombycis, 23 and 35 polymorphic products were observed, respectively, whereas the tetraploid species M. serrata revealed 20 ISSR products, and in the hexaploid, M. tiliaefolia, 25 products were amplified. These results are consistent with the earlier studies of RAPD patterns in Chrysanthemum by Wolf & Rijn , which show that the ploidy level of a plant does not appear to influence the number of fragments amplified per primer.
Classification of Morus based on phenotypic variations or isozyme patterns should be reconsidered in the context of molecular analyses by RAPD and ISSR as well as that of AFLP . Hirano  studied 131 varieties of cultivated mulberry morphologically classified as belonging to three species, viz., M. bombycis, M. alba, and M. latifolia. Examining isozymes and several sap proteins, he found no significant difference among these three species, prompting him to conclude that all three species are the same. However, our study as well as the AFLP study unambiguously place M. alba as a separate species distinct from the other two, which are likely to be independent species. Thus, RAPD and ISSR based molecular markers were able to distinguish differences between the species which were indistinguishable by isozyme based markers. Similarly, Koidzumi [3, 4] considered that M. lhou, M. multicaulis, and M. latifolia are similar and belong to a single species. However, their wide separation in the present study contradicts this observation. Additional phylogenetic studies using chloroplast or mitochondrial gene sequences or appropriate nuclear gene sequences can help to evaluate the systematic positions of these species.
The present study reveals that PCR based fingerprinting techniques, RAPD and ISSR, are informative for estimating the extent of genetic diversity as well as to determine the pattern of genetic relationships between different species of Morus, with polymorphism levels sufficient to establish informative fingerprints with relatively fewer primer sets. The information obtained from the present study could be of practical use for mapping the mulberry genome as well as for classical breeding. The genetic similarity of wild mulberry with other cultivated species is low as indicated by the separation of these two groups in both RAPD and ISSR analyses. The informative primers identified in our studies will be useful in genetic analysis of mulberry accessions in germplasm holdings. The putative species-specific bands can be used as probes to ascertain whether they are in low or high copy numbers in the mulberry genome, and such specific bands may be used for genotype characterization and grouping germplasm accessions. Further, putative species-specific RAPD markers could be converted to sequence characterized amplification regions (SCARs) after sequencing and designing primer pairs to develop robust species specific markers. The study also provides a basis for mulberry breeders to make informed choices on selection of parental material based on genetic diversity to help overcome some of the problems usually associated with a tree crop improvement program.
The wild and cultivated mulberry species were obtained from the Central Sericultural Germplasm Research Center, Hosur, India. The species used in the present study along with their key morphological characters and country of origin are listed in Table 1.
DNA was extracted from fresh young leaves of five individuals of a species using a Nucleon Phytopure System (Amersham Pharmacia Biotech, UK) according to the manufacturer's instruction and pooled. The genomic DNA was quantified on 0.8% agarose gels and diluted to uniform concentration (10 ng/μl) for RAPD and ISSR analysis.
PCR reactions were performed according to the protocol of Williams et al. . Briefly, PCR amplifications were carried out in an MJ Research Thermal Cycler PTC-200, in a reaction volume of 20 μl containing 20 mM Tris-HCl (pH 8.4), 50 mM KCl, 2 mM MgCl2, 0.2 μM primer, 0.1 mM each of dATP, dTTP, dCTP and dGTP, 0.5 U of Taq DNA polymerase (Gibco BRL, USA) and 20 ng of template DNA. The 19 random primers used in the current study were obtained from Operon Technologies Inc. (Alameda, CA, USA); these included OPW-1-4, OPY-1-11, OPY-13-14, and OPY-16-17. Amplification reactions were carried out using the following cycle profile: Initial denaturation at 93°C for 2 min followed by 45 cycles at 93°C for 1 min, 35°C for 1 min, 72°C for 2 min and a final 7 min extension at 72°C. PCR products were electrophoresed on a 1.5% agarose gel according to Sambrook et al.  in 1X TBE buffer (89 mM Tris-borate, 2 mM EDTA, pH 8.0) and stained with ethidium bromide. The gel image was recorded using a Gel Documentation System (UVP, UK).
The ISSR primers were synthesized on an Applied Biosystem DNA synthesizer based on core repeats , anchored either at the 5' or 3' end (Table 2). Amplification reactions were carried out in 20 μl containing 10 mM Tris-HCl (pH 8.4), 50 mM KCl, 0.1% Triton X-100, 2 mM MgCl2, 10 μM of each primer, 0.1 mM each of dCTP, dGTP, dTTP and 0.075 mM of cold dATP with 2 μCi of α-32P dATP, (6000 Ci/mmol from BRIT, JONAKI, Hyderabad, India), 1 U of Taq DNA polymerase (Gibco BRL, USA) and 20 ng of template DNA. PCR reactions were performed using an MJ Research Thermal Cycler PTC-200 with following amplification conditions: Initial denaturation of 2 min at 94°C, followed by 30 cycles of 94°C for 30s, annealing at 52°C for 45s, extension at 72°C for 2 min and a final extension at 72°C for 7 min. PCR reactions were terminated by adding 13.2 μl of stop solution (95% formamide, 20 mM EDTA, 0.5% bromophenol blue and xylene cyanol). 4 μl of sample DNA was denatured at 75°C for 2 min, chilled on ice and then run on a non-denaturing sequencing gel containing 6% polyacrylamide, 3 M urea and 1X TBE at 900 V for 20 hr. The gels were dried and exposed for 2–10 hr on X-ray film (Kodak-Biomax) for autoradiography.
DNA banding patterns generated by RAPD and ISSR were scored for the presence (1) or for absence (0) of each amplified band. All RAPD and ISSR assays were repeated twice and only the reproducible bands were scored. For considering a marker as polymorphic, the absence of an amplified product in at least one species was used as a criterion. For genetic distance analysis, WINDIST software of the WIN BOOT package was used with the NTSYS format. Cluster analysis was based on similarity matrices using the unweighted pair group method analysis (UPGMA) program in the WIN BOOT software package . The Dice coefficient was used for dendrogram construction with a sample number of 100.
The authors thank the Central Silk Board (Ministry of Textiles – Govt. of India), Bangalore and the Department of Biotechnology, Govt. of India, New Delhi for financial assistance (to JN) for the study. We also thank Prof. Marian Goldsmith for critical reading of the manuscript and useful suggestions.
- Hou YJ: Mulberry breeding. Sericulture Department, Zhejiang Agriculture University, Hangzhou, China. 1994, 4-
- Linnaeus C: Species Plantarum. 1753, Stockholm, 2: 986-
- Koidzumi G: Taxonomical discussion on Morus plants. Bull Imp Sericult Exp Stat. 1917, 3: 1-62.
- Koidzumi G: Synopsis specierum generis Mori. Bull Imp Sericult Exp Stat. 1923, 11: 1-50.
- Hotta T: Taxonomical studies on the Morus plants and their distributions in Japan and its viscinties. Japanese Society for Promotion of Science, Ueno Park, Tokyo. 1958, 1-161.
- Katsumata T: Mulberry species in west Jawa and their peculiarities. J Seric Sci Jpn. 1972, 42: 213-23.
- Airy Shaw HK: A dictionary of flowering plants and ferns. Edited by: Willis JC, Airy Shaw HK. 1973, Cambridge University Press, London, 761-8th
- Hirano H: Thremmetological studies of protein variation in mulberry. Bull Sericul Exp Sta. 1980, 28: 67-186.
- Wang ZY, Tanksley SD: Restriction fragment length polymorphism in Oryza sativa L. Genome. 1989, 32: 1113-8.View Article
- Jana S, Pietrzak LN: Comparative assessment of genetic diversity in wild and primitive cultivated barley in a centre of diversity. Genetics. 1988, 119: 981-90.PubMed CentralPubMed
- Chengyin L, Weihua L, Mingjum L: Relationship between the evolutionary relatives and the variation of esterase isozymes in tea plant. J Tea Sci. 1992, 12: 15-20.
- Jelinski DE, Cheliak WM: Genetic diversity and spatial subdivision of Populus tremuloides (Salicaceae) in a heterogeneous landscape. Am J Bot. 1992, 79: 728-36.View Article
- Botstein D, White RL, Skolnik M, Davis RW: Construction of a linkage map in man using restriction fragment length polymorphism. Am J Hum Genet. 1980, 32: 314-31.PubMed CentralPubMed
- Tanksley SD, Young ND, Paterson AH, Bonierbale MW: RFLP mapping in plant breeding: New tools for an old science. Biotechnology. 1989, 7: 257-64.View Article
- Miller JC, Tanksley SD: Effects of restriction enzymes, probe source, and probe length on detecting restriction fragment length polymorphism in tomato. Theor Appl Genet. 1990, 80: 385-9.PubMed
- Wang ZY, Second G, Tanksley SD: Polymorphism and phylogenetic relationships among species in the genus Oryza as determined by analysis of nuclear RFLPs. Theor Appl Genet. 1992, 83: 565-81.View ArticlePubMed
- Beckmann JS, Soller M: Restriction fragment length polymorphism and genetic improvement of agricultural species. Euphytica. 1986, 35: 111-24.View Article
- Williams JG, Kubelik AR, Livak J, Rafalski JA, Tingey SV: DNA polymorphism amplified by arbitrary primers are useful as genetic markers. Nucleic Acids Res. 1990, 18: 6531-5.PubMed CentralView ArticlePubMed
- Welsh J, McClelland M: Fingerprinting genomes using PCR with arbitrary primers. Nucleic Acids Res. 1990, 18: 7213-8.PubMed CentralView ArticlePubMed
- Zietkiewicz E, Rafalski A, Labuda D: Genome fingerprinting by simple sequence repeat (SSR-) anchored polymerase chain reaction amplification. Genomics. 1994, 20: 176-83. 10.1006/geno.1994.1151.View ArticlePubMed
- Vos P, Hogers R, Bleeker M, Reijans M, van de Lee T, Hornes M, et al: AFLP: a new technique for DNA Fingerprinting. Nucleic Acids Res. 1995, 23: 4407-14.PubMed CentralView ArticlePubMed
- Weber JL, May PE: Abundant class of human DNA polymorphisms which can be typed using the polymerase chain reaction. Am J Hum Genet. 1989, 44: 388-96.PubMed CentralPubMed
- Wu KS, Tanksley SD: Abundance, polymorphism and genetic mapping of microsatellites in rice. Mol Gen Genet. 1993, 241: 225-35.View ArticlePubMed
- Akkaya MS, Bhagwat AA, Cregan PB: Length polymorphism of simple sequence repeat DNA in soybean. Genetics. 1992, 132: 1131-9.PubMed CentralPubMed
- Morgante M, Olivieri AM: PCR-amplified microsatellites as markers in plant genetics. Plant J. 1993, 3: 175-82. 10.1046/j.1365-313X.1993.t01-9-00999.x.View ArticlePubMed
- Plaschke J, Ganal MW, Roder MS: Detection of genetic diversity in closely related bread wheat using microsatellite markers. Theor Appl Genet. 1995, 91: 1001-7.PubMed
- Beckmann JS, Soller M: Toward a unified approach to genetic mapping of eukaryotes based on sequence tagged microsatellite sites. Biotechnology. 1990, 8: 930-2.View ArticlePubMed
- Roder MS, Plaschke J, Konig SU, Boner A, Sorrells ME, Tanksley SD, Ganal MW: Abundance, variability and chromosomal location of microsatellites in wheat. Mol Gen Genet. 1995, 246: 327-33.View ArticlePubMed
- Tsumura Y, Ohba K, Strauss SH: Diversity and inheritance of inter-simple sequence repeat polymorphism in Douglas-fir (Pseudotsuga menziesii) and Sugi (Cryptomeria japonica). Theor Appl Genet. 1996, 92: 40-5.View ArticlePubMed
- Nagaoka T, Ogihara Y: Applicability of inter-simple sequence repeat polymorphisms in wheat for use as DNA markers in comparison to RFLP and RAPD markers. Theor Appl Genet. 1997, 94: 597-602. 10.1007/s001220050456.View Article
- Nagaraju J, Kathirvel M, Kumar RR, Siddiq EA, Hasnain SE: Genetic analysis of traditional and evolved Basmati and non-Basmati rice varieties by using fluorescence-based ISSR-PCR and SSR markers. Proc Natl Acad Sci U S A. 2002, 99: 5836-41. 10.1073/pnas.042099099.PubMed CentralView ArticlePubMed
- Sharma A, Sharma R, Machii H: Assessment of genetic diversity in a Morus germplasm collection using fluorescence-based AFLP markers. Theor Appl Genet. 2000, 101: 1049-55. 10.1007/s001220051579.View Article
- Bhattacharya E, Ranade SA: Molecular distinction amongst varieties of Mulberry using RAPD and DAMD profiles. BMC Plant Biol. 2001, 1: 3-10.1186/1471-2229-1-3.PubMed CentralView ArticlePubMed
- Vijayan K, Chatterjee SN: ISSR profiling of Indian cultivars of mulberry (Morus spp.) and its relevance to breeding programs. Euphytica. 2003, 131: 53-63. 10.1023/A:1023098908110.View Article
- Lou CF, Zhang YZ, Zhou JM: Polymorphisms of genomic DNA in parents and their resulting hybrids in mulberry Morus. Sericologia. 1998, 38: 437-49.
- Morgante M, Hanafey M, Powell W: Microsatellites are preferentially associated with nonrepetitive DNA in plant genomes. Nat Genet. 2002, 30: 194-200. 10.1038/ng822.View ArticlePubMed
- Thormann CE, Ferreira ME, Camargo LEA, Tivang JG, Osborn TC: Comparison of RFLP and RAPD Markers to Estimating Genetic Relationships Within and Among Cruciferous Species. Theor Appl Genet. 1994, 88: 973-80.View ArticlePubMed
- Rieseberg LH: Homology among RAPD fragments in interspecific comparisons. Mol Ecol. 1996, 5: 99-105.View Article
- Wolff K, Peters-Van Rijn J: Rapid detection of genetic variability in Chrysanthemum (Dendranthema Grandiflora Tzvelev) using random primers. Heredity. 1993, 71: 335-41.View ArticlePubMed
- Sambrook J, Fritsch EF, Maniatis T: Molecular cloning A laboratory manual. 1989, New York: Cold Spring Harbor Laboratory, Cold Spring Harbor, 2
- Yap IV, Nelson RJ: WINBOOT: A program for performing bootstrap analysis of binary data to determine the confidence limits of UPGMA-based dendrograms. 1996, IRRI Discussion paper Series No. 14, International Rice Research Institute, Manila, Philippines
This article is published under license to BioMed Central Ltd. This is an Open Access article: verbatim copying and redistribution of this article are permitted in all media for any purpose, provided this notice is preserved along with the article's original URL.