Abstract
Random amplified polymorphic DNA markers (RAPD) were used to assess the genetic relatedness among accessions of four species of Desmodium and allied genera (Dendrolobium triangulare, Desmodium gangeticum, Desmodium heterocarpon ssp. heterocarpon, and Tadehagi triquetrum) originating from Northeast Vietnam. Since information on the genetic diversity of these species is deficient, the creation of baseline data is an important means for the development of more sustainable and cost-efficient conservation approaches which eventually result in more comprehensive ex situ germplasm collections. The species analyzed are native to tropical and subtropical Asia, Australia, and Oceania and possess a potential as forage and/or medicinal plants. Moderate levels of inter-accession diversity represented by 37.5% and 33.3% of polymorphic fragments (P%) and average Jaccard's similarity coefficients (JSCs) of 0.60 and 0.64 were found in D. heterocarpon and T. triquetrum, respectively, while moderate to high levels were detected in D. triangulare (P% = 52.9 and JSC = 0.61) and D. gangeticum (P% = 34.5 and JSC = 0.49). Mantel tests failed to reveal a correlation between geographic and genetic distances. Based on the results of this study, baseline data for further marker-assisted research are generated and future collecting and ex situ conservation strategies for the species studied are discussed.
Introduction
Within the Leguminosae, the tribe Desmodieae (Benth.) Hutch. is recognized for its substantial contribution to forage production in the tropics [43]. It comprises several species of agronomic interest such as Dendrolobium triangulare (Retz.) Schindl. (syn. Desmodium triangulare [Retz.] Merr.), Desmodium gangeticum (L.) DC., Desmodium heterocarpon (L.) DC. ssp. heterocarpon, and Tadehagi triquetrum (L.) H. Ohashi (syn. Desmodium triquetrum (L.) DC.) (Fig. 1). Due to their ability to adapt to low-fertility soils and to tolerate water drought, these species represent valuable plant genetic resources for forage production, soil protection and improvement in marginal smallholder farming systems of subhumid and humid tropical regions [10, 20, 41]. Moreover, the four legume species included in this study have a long history in traditional Asian health care [30, 33, 46].

Images of four species of Desmodium and allied genera analyzed with RAPDs: (a) Dendrolobium triangulare (Retz.) Schindl., (b) Desmodium gangeticum (L.) DC., (c) Desmodium heterocarpon (L.) DC. ssp. heterocarpon, (d) Tadehagi triquetrum (L.) H. Ohashi. Photos (a, d) by R. Schultze-Kraft, (b) by A. Ciprián (CIAT), and (c) by B. Heider.
With the exception of D. triangulare, a woody, perennial shrub, all analyzed species are woody herbs or subshrubs that mainly occur in fallow areas, thickets, forest borders, or along roadsides [4, 29]. The genus Desmodium and its allied genera, members of the subtribe Desmodiinae, are mainly native to the tropics and subtropics. Ohashi [28] records two centers of geographic distribution and differentiation: Southeast Asia and Mexico. According to Williams [47] Southeast Asia is a center of legume diversity. There, systematic exploration of native legume genetic resources started only about two decades ago when a series of collecting missions under the auspices of CIAT (International Center for Tropical Agriculture) was carried out in Malaysia, Indonesia, tropical China, Thailand, Papua New Guinea, and Vietnam. However, the mountainous and highly diverse north of Vietnam was not visited during previous collecting missions [37, 38]. At the same time, plant genetic resources of the North Vietnamese mountains are increasingly threatened due to ecosystem transformation and growing demand of natural resources in the wake of high population pressure [45]. For these reasons the area was chosen as the collection target in this study.
Molecular marker techniques have become increasingly important in conservation biology. Information about the underlying genetic diversity of species is essential for sustainable use and efficient conservation efforts of plant genetic resources [21, 40]. Unfortunately, such information is still inadequate concerning species of Desmodium and allied genera, especially at the DNA level. The traditional methodological approach to analyzing genetic diversity is based on morphological and/or agronomic traits. However, results of such studies might be biased by confounding environmental factors that influence the phenotype. Apart from D. heterocarpon, all Desmodium species so far analyzed were merely subjected to taxonomic, agronomic or, in the case of D. gangeticum, chemical studies [9, 22, 28, 41]. Isozyme electrophoresis allowed distinguishing D. heterocarpon ssp. ovalifolium from D. heterocarpon ssp. heterocarpon [27] but failed to discriminate the former at genotype level [17]. In a recent study, genetic diversity of a germplasm collection of D. heterocarpon ssp. ovalifolium held by CIAT was analyzed and a tentative core collection built up on the basis of random amplified polymorphic DNA (RAPD) marker data [6].
RAPDs provide a simple, fast, and cost efficient methodology to primarily screen the level of genetic variation of a great number of yet unknown genotypes [48]. Even though RAPDs are frequently criticized for unreliable reproducibility [15], dominance [25], and comigration [35], they have proven to be a useful tool in studies analyzing genetic variation in a substantial number of tropical legume species (e.g., [8], [16], [23]).
For these reasons, RAPD markers were used in the present study to determine the genetic relatedness among accessions of four Desmodium and allied genera collected in Northeast Vietnam. The RAPD methodology applied will also be useful for screening duplication and establishing core collections in gene banks.
Methods

Collecting sites of Desmodium and allied genera in Bac Kan province, Northeast Vietnam.
Ten plants of each accession were propagated from scarified seeds and maintained under greenhouse conditions. Fresh leaf material was sampled from six-month-old plantlets. In the case of T. triquetrum and D. heterocarpon, genetic analyses were performed on a maximum number of 10 individual plants per accession in order to obtain sufficient leaf material for DNA extraction. Growth of D. gangeticum and D. triangulare was more vigorous, yet few individuals failed to survive. In order to standardize methodology, the number of individuals per bulk was limited to five plants per accession forming the a bulk. The surplus was placed in the b bulk. DNA of a and b bulks was separately isolated and analyzed. The division into an a and b bulk was also used to approximate the ideal number of plants per bulk.
Polymerase chain reaction (PCR) amplifications were performed in reaction mixtures containing 10 ng template DNA, 0.2 μM primer, 100 μM dNTP's, 10 mM PCR buffer, 1.5 mM MgCl2, and 1 unit Taq-polymerase in a total volume of 10 μl. Amplifications were carried out in a thermo-cycler (T Gradient, Biometra) with the following program: initial denaturation step of 3 min at 94 °C followed by 34 cycles of 1 min at 94 °C, 1 min at 36 °C, 2 min at 72 °C, and a final extension step of 5 min at 72 °C. DNA fragments were resolved by electrophoresis in 1.5% agarose gels run at 4 V cm1. Gels were stained with ethidium bromide, visualized under UV-light, and photographed by a CCD video camera (Raytest). To ensure reproducibility of banding patterns, each amplification reaction was repeated twice. Four primer kits (L, M, N, and U; Roth), each comprising 20 decamer primers, were screened for each of the four analyzed species. Only those primers were selected that showed well-defined and reproducible banding patterns (Table 1).
Oligonucleotide primers employed in RAPD analysis, their sequence, number of bands obtained per primer, size range of bands, and specific size of polymorphic (polym.) fragments.
Results
Dendrolobium triangulare (Retz.) Schindl.

Examples of RAPD profiles among accessions of Dendrolobium triangulare, Desmodium gangeticum, D. heterocarpon, and Tadehagi triquetrum obtained with the respective primers. Numbered lanes display the RAPD profile of different accessions; lane M: molecular weight marker (100 bp DNA-ladder, Roche Diagnostics); lane C: negative control (H2Odest. instead of template DNA). *Polymorphism at 1200 bp not shown on this gel.

UPGMA dendrogram illustrating genetic relationships among 35 accessions of Dendrolobium triangulare. The numerical scale of JSCs indicates increasing similarity.
Desmodium gangeticum (L.) DC.

UPGMA dendrogram illustrating genetic relationships among 37 accessions of Desmodium gangeticum. The numerical scale of JSCs indicates increasing similarity.
Desmodium heterocarpon (L.) DC. ssp. heterocarpon

UPGMA dendrogram illustrating genetic relationships among 66 accessions of Desmodium heterocarpon. The numerical scale of JSCs indicates increasing similarity.
Tadehagi triquetrum (L.) H. Ohashi

UPGMA dendrogram illustrating genetic relationships among 24 accessions of Tadehagi triquetrum. The numerical scale of JSCs indicates increasing similarity.
Discussion
RAPD analyses were useful to determine genetic relatedness among accessions of four species of Desmodium and allied genera collected in North Vietnam. Moderate levels of diversity existed among accessions of D. heterocarpon and T. triquetrum while moderate to high levels were found in D. gangeticum and D. triangulare.
The level of polymorphism observed among accessions of D. gangeticum (34.5%), D. heterocarpon (37.5%), and T. triquetrum (33.3%) was among the lowest compared to many studies of tropical herbaceous legume species. Liu [23, 24] reported 68.9% of polymorphism among Lablab purpureus and 86.6% among Stylosanthes scabra genotypes. RAPD analyses of Vigna angularis (49%), Pueraria montana var. lobata (54.3%), and P. phaseoloides (45.5%) yielded slightly higher levels of polymorphism [13, 49]. Comparisons of D. triangulare (52.9%) with RAPD marker studies of other tropical legume trees and shrubs resulted in a relatively high level of polymorphic loci for this species. While higher levels of polymorphism were reported for Cratylia argentea (91.8%) [1], and Flemingia macrophylla (95%) [12], lower levels were found in Argentinian Acacia species (19.1%) [3] and Prosopis species (28.3%) [16]. However, it should be recognised that the chosen primer selection approach may not have resulted in unbiased estimates of levels of polymorphism. Preliminary primer screens were run on a selection of 10 accessions per species and only those primers that generated polymorphic fragments were selected. In the case of uniform amplification, primers were discarded and not used for full screening of all accessions which might have interfered with data accuracy.
Average JSCs obtained in this study (D. heterocarpon, 0.60, T. triquetrum, 0.64, and D. gangeticum, 0.49) were comparable to previous RAPD marker studies of herbaceous legumes. Xu et al. [49] found JSCs of 0.73 in cultivated forms and 0.60 in wild and weedy populations of V. angularis, while Yee et al. [50] and Heider et al. [13] recorded JSCs of 0.52 in V. angularis and P. phaseoloides, respectively. Results of a recent molecular marker study showed a comparatively low JSC in Pueraria montana var. lobata (0.35) [13]. RAPD analyses of tropical tree and shrub legumes revealed average JSCs that were slightly higher in F. macrophylla (0.67) [12] than in D. triangulare (0.61). However, the fact that a large number of previous studies on genetic diversity of tropical legumes did not disclose the overall JSC, used different diversity indices, or employed different molecular marker techniques, represents an obstacle to direct comparisons.
Several authors experimented with the optimal number of individual plants that were pooled in order to reduce the sensitivity of RAPDs to DNA sequence variation and thus to ensure a representative estimation of genetic diversity in heterogeneous populations [7, 18, 19]. The recommended number of individuals per bulk ranged between three and 20. Based on these references, bulk sizes in the study of D. heterocarpon and T. triquetrum was adjusted to 10 individuals in each bulk. Due to low germination rates and low seedling vigor, the number of surviving greenhouse individuals restricted the bulk sizes of D. triangulare and D. gangeticum to five individual plants each (a bulk). Any additional plants formed b bulks that consisted of up to five individuals. The fact that 70.4% of a and b bulks of D. triangulare and 43.2% of D. gangeticum were not identical and therefore did not coincide in the same cluster, illustrates clearly that bulk sizes were too small to dilute rare alleles. In general, Desmodieae are assumed to be predominantly self-pollinated, but outcrossing rates, ranging from 0.01% to 10%, are also likely to occur as in other tropical legumes [26, 32]. Given the rudimentary knowledge available on breeding systems of the analyzed species, the degree of heterozygosity inherent in accessions of D. triangulare and D. gangeticum seems to be much higher than originally assumed. Seeds were sampled as accessions representing the natural population from which they were collected. The aim of collecting accessions is to capture the entire range of genetic variation, as much as possible, in a given population [2]. Thus individual seedlings were bulked for RAPD analysis to represent the genetic constitution of the original accession. The exclusive use of bulk samples, however, does not allow the use of statistical parameters of population genetics and impedes the assessment of the species' breeding system. There was no substantial information available on the reproduction mode of the four studied species. However, such information is decisive for the amount and partitioning of genetic variation in and among populations [11]. Therefore, in order to draw firm conclusions on intra-accession diversity, individual plants composing an accession should additionally be analyzed in future studies, a measure that might aid in improving data quality and thereby contribute to optimizing collection and conservation management.
Implications for conservation
Desmodium and allied genera play an important role in different research and development endeavors to address declining soil fertility, erosion, overgrazing, and dry-season forage constraints in marginal areas of the tropics and subtropics. Desmodium heterocarpon ssp. heterocarpon attracted research interest due to its ability to persist in highly productive grass/legume mixtures in humid subtropical and tropical regions, while tolerating heavy grazing, acid soils, drought, temporary water-logging, and moderate frost periods [10]. The carpon desmodium cultivar “Florida” of this species obtained economic significance as a forage legume in peninsular Florida [39]. Tadehagi triquetrum (L.) H. Ohashi is agronomically interesting because of its capacity to provide crude protein for livestock feeding during dry seasons and for its good adaptation to acid, low-fertility soils [42]. Dendrolobium triangulare (Retz.) Schindl. provides fuelwood and is suitable for soil protection and improvement. It has minor forage potential due to variable levels of condensed tannins [22]. Desmodium gangeticum (L.) DC. was found to be highly productive in shady environments such as under plantation crops [41]. Particularly, D. gangeticum is considered a plant resource of pharmaceutical importance which serves as a remedy for diarrhea, dysentery, and several eye diseases. These curative properties are due to a number of alkaloids, a pterocarpanoid (gangetin), flavons, and isoflavanoid glycosides as active compounds with antifungal, antileishmanial, antioxidant, antiradical, and anticholinesterase activities [e.g., 9].
The mountainous north of Vietnam experienced dramatic environmental and social changes during the last decades threatening the biological diversity of the region [45]. Desmodieae as well as other wild-growing plant species are subjected to genetic erosion because of habitat destruction and resource overexploitation as a consequence of increasing demographic pressure. Since the vast majority of tropical forage legume cultivars are unimproved selections from wild populations [34], the enormous yet untapped legume diversity of the North Vietnamese mountainous regions should be explored in order to avoid reliance on a narrow germplasm base. This includes further germplasm collecting as well as efficient ex situ conservation measures. Conservation of plant genetic diversity in gene banks requires balancing carefully between numerically extensive collections and the feasibility of maintaining and managing germplasm storage. Knowledge about genetic diversity in a target taxon is also indispensable when considering the establishment of a core collection or eliminating redundancy in collections [21, 40].
The heterogeneous clustering of a and b bulks in D. triangulare (Fig. 4) and D. gangeticum accessions (Fig. 5) imposes difficulties to the establishment of an ex situ conservation approach for these species. But generally, genotypes that are redundant may be removed from ex situ storage while outliers and particularly distinct accessions should be included. Clusters in which numerous extremely related accessions occur may be condensed to subsamples. However, such selection must be carried out with circumspection as RAPD diversity does not necessarily reflect diversity at quantitative trait loci, and important genotype characteristics might become lost [44]. Thus morphological plant characteristics should also be considered before discarding germplasm. In the case of D. triangulare and D. gangeticum, each cluster VII certainly encompasses the most distinct accessions (Figs. 4 and 5). Accessions that did not represent unique genotypes, i.e., those that formed clusters at a JSC of 1, are candidates for removal. The D. triangulare dendrogram contains several of these, particularly in clusters I and II. The first three clusters of D. gangeticum could also be reduced.
Within the D. heterocarpon (Fig. 6) collection accessions 491, 512, 544, and 743 were the most distinct accessions and therefore suggest themselves for ex situ conservation. Apart from these, a representative subsample of each cluster, eliminating the numerous duplicates, would suffice.
A comparison of passport data between the two main clusters of D. heterocarpon revealed that accessions forming the upper division were collected from humid and frequently shady habitats growing in heavy soils characterized by increased nitrogen, calcium, magnesium, and organic matter contents while accessions from the lower division grew in poorer soils and in more marginal, sun-prone, and therefore drier environments (data not shown). Thus it appears that the two main divisions consist of two different ecotypes.
High genetic differentiation was found among three of 24 T. triquetrum accessions (944, 945, and 918), while the remaining accessions harbored high similarity values among them. The outlier (918) of T. triquetrum exhibited three unique markers and as a consequence may contain unique genetic characteristics, which indicates that significant distinctiveness was captured here. However, neither taxonomic plant characteristics nor collecting site-related passport data could explain this distinctiveness. In general, data show that T. triquetrum represented the least diverse collection among the four species analyzed but also the smallest. Therefore, it seems advisable to conduct further, more extensive collecting missions at larger distances in order to collect more genetic diversity that would qualify for ex situ conservation. As for the samples collected in Bac Kan province, the outlier, accessions comprising cluster III, and a representative subsample of clusters I and II (Fig. 7) should be included if aiming at efficient ex situ conservation.
Future collecting missions should target populations at large distances covering altitudinal gradients in order to ensure that the ecological and geographic ranges of the respective species are included [2]. Outcrossing must be assumed in all four species studied. Hence, individuals within distant populations should be comprehensively collected instead of extracting a few individuals from many populations at short distances. A manageable and efficient ex situ conservation approach aims at capturing maximum genetic diversity within stored material, comprising outliers and genetically distinct germplasm, and representative subsamples of very similar or even redundant genotypes. Ideally such ex situ conservation would be based on both molecular and morphological data in order not to lose valuable adaptive variation.
Concluding remarks
Generated data showed that RAPD markers, despite their widely discussed disadvantages, are suitable tools to scan anonymous DNA of species with high secondary metabolite content such as Desmodium and allied genera. However, limited reliability of RAPD molecular markers should be improved by increasing the number of RAPD fragments. As no correlation between the pattern of genetic clustering and geographic distance or passport data was detected, further data analyses will be necessary in order to relate genetic diversity to greater geographic distance. Increasing emphasis is being placed on empirical data of genetic diversity to assess and manage genetic resources in a sustainable way. Therefore this study may contribute to provide germplasm conservationists with data in order to facilitate conservation efforts.
Footnotes
Acknowledgements
Research was financially supported by the Vater-und-Sohn-Eiselen-Stiftung, Ulm, Germany. The authors are also thankful to the Vietnam Agricultural Science Institute, Hanoi, Vietnam, especially to Prof. Dr. Tran Dinh Long and Mr. Ha Dinh Tuan for their advice und cooperation during germplasm collection.
