Many species showed evidence of genetic structuring, particularly when there are physical barriers to dispersal [1–3] and mating constraints [4,5]. Physical barriers may exist between regions as the presence of natural barriers such as mountains, valleys, cliffs or oceans that lead to the restriction of gene flow [1–3] or may also be caused by habitat fragmentation due to anthropogenic disturbances [5,6]. In such a case, genetic differences occurred as a result of different patterns of allele frequency distribution between populations . In other cases, genetic differentiation may exist as a result of mating barriers due to the differences on the floral structures and biology [5,7–9] among individuals within a population, that in turn lead to a mating incompatibility [4,5]. In such cases, reproductive failure could more be attributed to the mating barrier of pollination due to the differences on stylus length , structures of the corolla, filament length, stigma-anthers position , pistil dysfunction, male sterility  and / or asynchrony on sexual organs maturity .
Distributed naturally along China, India, Indonesia and The Philippines , an economic-important species Santalum album L. (Santalaceae)—formerly called sandalwood—has been categorized as threatened species due to its significant degradation and / or habitat loss caused by demand on its wood and oil . Its heartwood, containing 1.5 % to 5 % of β-santalol—a strong-specific fragrance of oil—is widely used for wood carving, religious and medicinal purposes. Its oil has been used as materials for cosmetics, prime sources of perfumes and aroma therapy, and was presumed to contain anti-melanoma compounds [4,11]. Recently, this species found to be extinct in the wild in most of its native in Eastern parts of Indonesia .
Despite a significant degradation in its origin in the South-eastern islands of Indonesia, new landraces of sandalwood emerged in Gunung Sewu Geopark, a 1 300 km mountainous limestone zones in the central part of Java island, Indonesia. The oldest sandals' herbarium specimen (dated by the year 1853) in Java island was found in Imogiri District, and another specimen (dated by the year 1960) was collected from Nglipar District; both were part of Gunung Sewu Geopark. Previously, geographical events which were started from about 1.8 10 yr ago—involved tectonic movements, volcanic activities, and seawater erosion—has derived this area into various landscape structures differed in altitude, elevation, soils and microclimate conditions [13,14]. Considering a rapid degradation of sandalwood population in its origin, the occurrence of new landraces in Gunung Sewu may provide a promising source for any of reintroduction and rehabilitation efforts.
Many studies showed clear evidence that differences in floral structures may affect the mating systems [15,16], which in turn resulted in different pollination success and reproductive outputs [5–9]. Furthermore, the differences on floral sexual organs may affect pollen-pistil interaction which in turn resulting in interspecific-mating incompatibility, as reported for S. album in India , S. album, S. lanceolatum and S. spicatum in Western Australia  and S. lanceolatum in Victoria Australia . Sandalwood in Gunung Sewu occurred in various types of landscapes; some of them were bordered by natural barriers such as dense woody forests and rocky limestone walls. Beside of these physical barriers, sandalwood in Gunung Sewu tended to experience mating constraints since some of the population failed to produce mature fruits. Each population in Gunung Sewu consisted of at least three sandal variants that were distinguished by their floral and leaves structures. The preliminary study divided these variants into three groups (YBF, refers to “yellow big flower”; RBF, “big red flower”; and RSF, “red small flower,” respectively), which were differed in floral structures and longevity. These differences were considered to be under genetic controls, while the variation among sites was affected more by environmental differences. In order to arrange the genetically-based conservation strategy, this previous finding recommended further investigation of genetic differences and crossing abilities among variants. In this study, the study combined field observations, hand-pollination with biochemical marker-gene analysis to compare genetic diversity, mating systems and cross ability of three sandals variants among four populations representing geographical zones in Gunung Sewu Geopark, Indonesia.
2. Materials and Methods
Gunung Sewu consisted of more than ten sandalwood populations in the form of both planted and naturally regenerated stands. However, study sites only compared four populations which were representing distinctively different population structures: one of population (the basin of Bleberan) in the Middle Zone, two (the highland of Nglanggeran and the lowland cave of Bejiharjo) in the Northern Zone, and one (the karst area Petir) in the Southern Zone, respectively. Each of population is separated by 25 km to 40 km. These sites are at different altitude, experienced different climatic regimes and having clear ecological differences (Table 1; Fig. 1).
Sandalwood is a long-lived, yearly ﬂowering perennial of semi-arid and tropical region throughout south-eastern and middle parts of Indonesia. It is a shrub attaining a height of about 12 m and a girth of 10 cm to 30 cm [4,11], mainly vegetatively reproduced by root suckers [5,6], particularly under marginal condition. Generally, sandalwood in Indonesia flowered twice a year: at the beginning of dry season in May to September, and at the beginning of rainy season in November to March; with 4 mo to 5 mo flowering period [20,21]. Flowers pollinated by insects belong to hymenopterans, lepidopterans and dipterans [16,25]. The previous study showed evidence of dichogamy , highly outcrossing rate and self-incompatibility [3,5,11,12]. However, partially self-compatibility has been reported for S. accuminatum  and S. album [21,23], particularly under isolated condition. Santalum has strong ability to produce root suckers; the vegetatively propagated sprouting emerged from the roots . Clonality event occurred when most of the off-springs in the population were derived from a genetically identical individual [5,6]. Previous researches documented pollen limitation due to less pollen production [16,18] and male sterility . Its reproductive success was also very low, ranging from 0.03 % in natural population to less than 20 % in ex-situ plantation , however, cross-pollination tended to enhance seed set .
The preliminary study divided sandalwood in Gunung Sewu into three variants (YBF, refers to “big yellow flower”; RBF, “big red flower”; and RSF, “red small flower,” respectively) differed by floral structures. RSF and RBF dominated by red and maroon colors, while YBF is more yellowish to orange. RBF and YBF possessed bigger perigonium, longer sexual organs but with shorter longevity, and similar/lower position of the stylus to stamens. RSF flowers are the smaller, similar/higher position of the stylus and longer longevity. RSF produced more flowers per inflorescence. Sandalwood flowered twice a year in all of sites and variants, however, the onset and duration varied. YBF flowered earliest while RBF was the latest. RSF possessed the longest flowering period. Flowering differences among variants were considered to be under genetic controls, while the variation among sites was affected more by environmental differences. Flowering varied among sites due to the altitude, soils and climatic differences. Sandalwood in lower altitude, drier and warmer sites flowered earlier and shorter.
Population structure measurements
Measurements in early 2016 were made to determine the population width and size [5,23]. Populations were defined as spatially discrete clusters of plants, separated from others by at least 500 m . A 10 m 10 m grid was made on each of 1 ha of width, where the number of individuals was counted . Measurements in effective population size as the number of flowering individuals divided by the total number of adult individuals . Population density conformed to the number of single individuals per ha . To measure the extent of clonality, it will be counted the proportion of vegetatively vs. reproductively propagated individuals within each of populations. Visibly separated stands  which have different vegetative and sexual organs morphology  were regarded as an individual, whereas stem separated by less than 10 cm  and/or shared the same root systems  were regarded as part of the same individual. Individual maturity determined by the presence of flowering [5,23].
Hand-pollination treatments and reproductive parameters measurements
The degree of cross- and self-compatibility was assessed over the 2016's flowering period by performing hand pollination experiments: i) hand self-pollination; ii) hand cross-pollination; and iii) natural open pollination . To support cross-pollination, filaments were emasculated at anthesis to remove immature anthers. Receptive stigma was pollinated by a selected pollen source (male parent) by applying the mature anther to the stigma in 3 times replication to ensure the pollen transfers. Each inﬂorescence was then tagged and isolated using fine fabrics. To perform self-pollination, single inflorescence was isolated using fine fabrics prior to anthesis. Hand pollination was conducted in four populations (total N for self hand-pollination = 279 single flowers in Nglanggeran, 551 in Bejiharjo, 301 in Bleberan, and 303 in Petir; N for open pollination = 257 single flowers in Nglanggeran, 405 in Bejiharjo, 214 in Bleberan, and 209 in Petir; N for cross hand-pollination = 300 in each population). These four sites represented each of the geographical zones in Gunung Sewu. To measure reproductive success of the seed sets, the flowers and mature seeds were counted and the value of Reproductive Success (RS)—correspond to the ability of reproductive organs to form mature seeds—was measured following the formula, RS = (Fruit/Flower) (Seed/Ovule).
Genetic diversity was measured spatially across four populations along geographical gradients in Gunung Sewu. To measure genetic diversity of parents, at the end of rainy season in June 2016, juvenile leaves were sampled from randomly chosen individuals. In order to avoid the re-sampling of the same siblings, particularly at the strongly clonalized populations, it supposed to be ensured that the sampled trees were propagated from sexual reproduction. The flowering trees were of the first priority to be sampled. Samples were wrapped, frozen in ice packs and taken to the laboratory for allozyme extraction and electrophoresis. Previous study gained three enzymes, shikimate dehydrogenase (E.C. 22.214.171.124.), esterase (E.C. 3.1.1.) and diaphorase (E.C. 126.96.36.199.) which observed to be polymorphic. Zymogram phenotypes that were interpretable were found for only six loci, shikimate dehydrogenase Skd-1, esterase Est-1, Est-2, and Est-3, and diaphorase Dia-1 and Dia-2. Electrophoretic procedures were conducted with vertical polyacrilamide gel electrophoresis procedures following David-Ornstein method . The leaves were homogenized in modified extraction buffer and centrifuged at 1570.78 rad/sec for 15 min at 4 C. The supernatant was loaded onto polyacrylamide vertical slab (Sigma Inc., USA) gels and electrophoresed at 4 C at 220 V and 200 mA current for about 3 h. After electrophoresis, the gels were stained using staining solution of each enzyme system, and the allozyme gels were genetically interpreted.
At each of the allozyme locus, the frequency of each allele and the genotype were counted. For each locus the number of heterozygote genotype were counted and expressed as percent observed heterozygosity (H). The observed heterozygosity was then pooled and averaged over all loci to determine the percent observed heterozygosity for a population. The expected heterozygosity (H) for each locus and over all loci for each population in Hardy–Weinberg equilibrium was counted following the formula: H = 1 – Σ p, where p refers to the i's allele frequency of a given population. Total population heterozygosity (H) was counted following the formula: H = 1 – Σ p, where p refers to the i's allele frequency of the populations in total. Fixation index, the deviation from expected frequencies under Hardy–Weinberg equilibrium, was measured following the formula: F = 1 – H/ H. To measure genetic differentiation, genetic diversity between populations (D) was calculated by reducing the value of total heterozygosity (H) to heterozygosity within population (H). The mean proportion of total gene diversity at polymorphic loci due to differences between populations (G) was measured by dividing D to H
Test of analysis of variance (ANOVA, LSD test) performed among parameters within populations and within floral variants. ANOVA was conducted to determine differences on the population structures, observed heterozygosity and fixation index among sites and among floral variants. ANOVA was also conducted in determining the differences of Reproductive Success resulted from different hand-pollination treatments, among floral variants, and among sites. Similar analyses were conducted to examine the differences between intra- and inter-specific outcrossed hand-pollination, and the differences among each of floral variants as the female parent. Multiple linear regressions based on population size and degree of clonality were carried out to examine the relationship between the population structures and the predictor variables observed heterozygosity (Ho) and inbreeding coefficient (FIS). The same method was applied to estimate the effect of both parental observed heterozygosity (Ho) and parental inbreeding coefficient (FIS) to the offspring heterozygosity and seedling recruitments. Some of the variables were subjected to logarithmic transformation to obtain normal data distributions. For the multiple regression analysis, a backward stepwise procedure was applied, with the final model including only variables with a significant (p 0.05) effect on the dependent variable. All statistical analyses were carried out with SPSS (Version 16.0, SPSS Inc.).
Allele frequency and distribution
Six loci showed considerable polymorphism (Est-1, Est-2, Est-3, Dia-1, Dia-2, and Skd-1) in most of populations, except Petir in which Dia-1 loci is monomorphic. These loci were also polymorphic in most of flower variants, except YBF in which many of loci are monomorphic. The YBF variant possessed monomorphic loci for Est-2 in Bleberan and Nglanggeran site, and for Dia-1 in Petir. In Bejiharjo, all of loci, except Dia-1, were monomorphic for YBF variant. However, the monomorphic loci of each variants were differed with sites. In Petir, most of monomorphic loci (Est-2, Est-3, and Dia-1) occurred in RBF. While in Nglanggeran, the Est-2 loci was monomorphic for YBF and RSF variants. Allele distributions differed with floral variants. The high frequency of “c” allele of Dia-1, which was considered rare at any of sites and variants, was occurred in YBF variant at all sites. In Skd-1, allele “a” dominated the RBF and RSF variants, while those of YBF was dominated by allele “c” (data not shown).
Heterozygosity and fixation index
Observed heterozygosity varied significantly with sites (Ho 0.021 to 0.327; F stat. 0.008, P value 8.016; respectively). A moderate level of diversity was maintained among sites (H 0.200; H 0.265; D= 6.5 %; G= 24.41 %, respectively). In this study, populations which experienced heavy exploitation and/or grew on marginal-rocky conditions were more clonalized. They existed as genotypically-identical parent trees, and recruited individuals mostly by root suckers. In general, the more clonalized populations exhibited much lower level of heterozygosity and tended to be more inbreeder. The F values, which represented the parental inbreeding coefficient, were also differed among sites in line with the differences on heterozygosity and degree of clonality. Of the four populations (each consisted of three variants) tested for conformity to Hardy–Weinberg equilibrium, only the most clonalized ones (Petir and Bejiharjo) showed a significant departure due to an excess of homozygotes, indicating a very high level of parental inbreeding. In general, the more clonalized populations exhibited much higher level of parental inbreeding, and lower level of heterozygosity. However, there was no significant differences on heterozygosity among floral variants (H 0.227 to 0.279; F stat. 0.787, P value 0.246; respectively). Some of diversity was maintained among variants (H 0.248; H 0.265; D= 1.7 %; G= 6.45 %, respectively) (Table 2).
Controlled pollination: The crossability and reproductive outputs
Hand-self pollination treatments produced very few reproductive outputs, or even no seeds were produced. Particularly in the population which possessed higher heterozygosity and tended to be more outcrosser, no signs of fertilization were observed in any of the self-pollinated ﬂowers. All ﬂowers were abscissed right after the end of anthesis phase. However, the highly clonalized populations which tended to be more inbreeders were able to produce flowers in hand-self pollination (Figure 2).
The highest level of Reproductive Success (RS) was achieved from interspecific cross-pollinated flowers for all variants. The interspecific crossing, a cross-pollination within the same variant, performed the highest level of RS. Particularly in the more genetically diverse populations, the seed sets from the interspecific cross-pollination were even higher than those resulting from the open pollination. However, a very low level of RS was observed in the intra-specific cross-pollination, a crossing between variants. The RSF variant even observed to be incompatible at any of cross-pollination treatments, resulting in 0 % mean of RS. Mature fruits were only observed in the reciprocal cross-pollination of YBF RBF variants.
ANOVA showed no differences in the Reproductive Success among both floral variants and populations. However, a significant Reproductive Success was observed among hand pollination treatments, in which the interspecific cross-pollination always resulted in the highest RS. The three of floral variants showed similar cross-ability while reciprocally-crossed to other variants (Table 3).
Genetic diversity and mating systems
Differences in genetic structures were observed among various types of population structures in this study. This finding showed evidence that the more clonalized populations exhibited a much lower level of heterozygosity. However, departure from prediction was reported for some small and isolated populations with the wider genetic base, which tended to have a higher level of genetic diversity and lower fixation index. It is considered that the parental genetic base played a role in maintaining outcrossing within such isolated populations.
Parental genetic base, which in this study is regarded to the heterozygosity of parent trees, is strongly determined by the proportion of vegetatively vs. reproductively propagated parents in population and their genetic composition. The degree of clonality strongly affected the mating systems and thus determined the level of cross- and self-mating within the population. Since sandalwood has strong ability to reproduce vegetatively by root suckers, thus the landscape features which promote vegetative propagation, such as the rocky and shallow soils, may contribute to the formation of a large portion of root suckers which were genotypically identic, and therefore enhance the level of inbreeding within the population.
Mating systems, which in this study is regarded to the parental inbreeding coefficient, were affected by the degree of parental clonality and heterozygosity. Mean F values indicated that only population undergoes strong clonality were deviated from Hardy–Weinberg equilibrium, showing a high level of inbreeding and an excess of homozygosity. As also reported in S. insulare, the fixation index F was positive and significantly different from zero for all the isolated islands . The authors proposed that the assumption of an increased selfing in the isolated populations appeared to be the main explanatory factor. An evolution of the mating system towards self-pollination could result from a lack of pollinators in some islands  and lack of outcross-mating partners in the heavily clonalized [5,6] or lower genetic bases'  populations.
In this study, the clonalized populations exhibited lower heterozygosity and tended to be more inbreeders. Studies on other sandals species also reported similar influences of fragmentation and clonality to the mating systems. RFLP-based study on S. spicatum in Western Australia showed that isolated populations are less diverse than bigger populations . The F value which was close to zero (F -0.026) indicated a random mating, but inbreeding tended to dominate mating in some of the isolated populations . Microsatellite-based study on S. insular in insular Pacific reported that significant excess of homozygosity, indicating a high level of inbreeding, was only observed within isolated islands (H 0.28; F 0.12; 58 % level of clonalities) and not in mainland (H 0.49; F -0.07), respectively . RAPD-based analysis on five isolated and clonalized populations of S. lanceolatum in south-eastern Victoria, Australia reported a very low heterozygosity (H 0.06) due to a very high inbreeding depression .
Higher gene diversity and outcrossing rate despite the small population size can be partly related to the mating system of the outbreeding plants . The increase of selfing rate, which lowering heterozygosity, may occur in the outbreeding taxa populations which subjected to the low parental genetic base [1,28], strongly clonalized [5,6] or isolated condition [24,29]. Clonality events may favour high levels of geitonogamy which increase the homozygosity. Moreover, a high level of inbreeding in naturally outbreeding taxa may favour population bottlenecks, genetic drift and inbreeding depression , such conditions which increase rare and missing alleles [1,2,20], reduce heterozygosity [6,8] and reduce reproductive fitness [18,22]. Landscape structures regarded to the rock outcrops, slopes, and natural barriers may also lead to the differences of genetic diversity and mating systems .
Since sandalwood is naturally an outbreeder, it is hypothesized that it will have less ability to cope with inbreeding depression. However, as a result of selection conformed to “reproductive assurance theory” [29,24], the small, isolated, clonalized and/or lower genetic bases' populations will be more inbreeder, and therefore will lack the negative effects of inbreeding depression. Both populations which were highly clonalized, the Petir and Bejiharjo, were confirming the “reproductive assurance theory,” in which selfing taxa occurred in isolated or clonalized population to enhance the ability for self-fertilization. These two clonalized populations were the only populations observed to produce high selfed-seeds, particularly at Petir. Both clonalized populations were also produced very few, even no outcrossed-seeds, in hand cross-pollination. However, a very low number of selfed-seeds produced in Bejiharjo were considered as the effect of the early-acting inbreeding depression. In contrast, the evidence of late-acting inbreeding depression was observed in the bigger population Petir which has been exposed to a clonality condition and therefore has been dominated by the self-mating system, for a relatively longer time period. In this bigger and denser but highly clonalized population, the effect of inbreeding depression will not be expressed until the reproductive outputs reached the seedling level. The success of reproductive at populations which have been subjected to inbreeding for a long time period could be seen as a result of possession of pre-adapted breeding systems, as reported for C. euphrasioides in the highland of Andes, Chile , S. spicatum in Western Australia  and S. australocaledonium in New Caledonia .
Crossability between sandalwood floral variants
The hand-pollination experiments in this study indicated that sandalwood is an allogamous and self-incompatible species. Naturally, sandalwood observed to be a highly outcrossed species [3,16,17] with partially  or totally  self-incompatible. Factors that may reduce or inhibit outcross-pollination processes, therefore, could have considerable effects on reproduction aspects that result in the failure of reproductive processes and thus reducing the reproductive outputs [8,15].
Intraspecific crossed-pollination in this study resulted in a very low seed sets. The RSF variant even observed to be incompatible at any of cross-pollination treatments. Mature fruits were only observed in the reciprocal cross-pollination of YBF RBF variants. Size of the corolla, length of filament and anthers position were positively correlated to the quality of pollen transferred. Hence the heterostylous and heteranthery flowers tended to be more outcrosser [8,9]. In this study, YBF and RBF possessed similar floral structures: bigger perigonium, longer sexual organs but with shorter longevity, and similar/lower position of the stylus to stamens. Whilst, RSF flowers are the smaller, similar/higher position of the stylus, and possessed longer longevity. In this study, crossing between variants having similar floral structures gained more seed sets. Other study found differences on the size of reproductive structures on twenty-five biotically pollinated plants of the Chaco Serrano Forest, Argentina, that was resulting in the differences of seed sets . The mating systems of narrow endemic Anthirrhinum microphyllum was also strongly related to the characters of sexual organs .
A barrier to intermate, which lead to a mating-incompatibility, may also occur as a result of the difference in gene composition . These genetic differences existed as a result of different allele distribution between populations ; a case that has also observed in this study in which there were different patterns of allele frequency distribution among sandalwood variants. The existence of genetic differentiation among variants, along with mating incompatibility between variants, has brought clear evidence that population has had subdivided and differentiated into small groups of the population according to the floral variants. This population differentiation and fragmentation may dramatically reduce population size, which along with the genetic differentiation, are the main factors causing inbreeding and restricted gene flow [1–3]. Within-population genetic differentiation, which increases inbreeding and reduces heterozygosity, was also observed with S. spicatum  and S. austrocaledonicum . These results on clear clustering of sandal populations suggested that the gene ﬂow and mating opportunities might have been reduced . Studies have also suggested that different pattern of allele frequency distribution , a low seed sets and/or restricted gene flow [1,4] can create such genetic clustering in natural populations. Hence, the extent of the heterozygote deficit could in part be attributed to the lack of sufficient mating partners, which in turn leading to the successive inter-mating between small groups of plants.
The sandalwood floral variants in this study were considered to be classified in the different gene pool according to their inter-mating ability. Some authors proposed classifying each plant and its related species by gene pools rather than by formal taxonomy . Primary gene pool (GP-1) is according to the members of the “probably in the same species" that can intermate freely in which crossing and gene transfer is easy, hybrids are generally fertile with good chromosome pairing, and gene segregation is approximately normal. In secondary gene pool (GP-2), the members are probably normally classified as different species to the primary gene pool. However, these species are closely related and can cross and produce at least some fertile hybrids though there are some reproductive barriers. The members of the tertiary gene pool (GP-3) are more distantly related to the primary gene pool. The primary and tertiary gene pools can be intermated, but gene transfer between them is almost impossible.
5. Conclusion and Recommendation
Observed heterozygosity varied significantly with sites (H 0.021 to 0.327; F stat. 0.008, P value 8.016; respectively). A moderate level of diversity was maintained among sites (H 0.200; H 0.265; D= 6.5 %; G= 24.41 %, respectively). The more clonalized populations exhibited a much lower level of heterozygosity and tended to be more inbreeder. The F values, which represented the parental inbreeding coefficient, were also differed among sites in line with the differences in heterozygosity and the degree of clonality. However, there were no significant differences in heterozygosity among floral variants (H 0.227 to 0.279; F stat. 0.787, P value 0.246; respectively). Some of diversity was maintained among variants (H 0.248; H 0.265; D= 1.7 %; G= 6.45 %, respectively).
No differences found on the Reproductive Success among both floral variants and populations. However, a significant Reproductive Success was observed among hand pollination treatments. The highest level of Reproductive Success (RS) was achieved from interspecific cross-pollinated flowers for all variants. However, a very low level of RS was observed in the intraspecific cross-pollination. The RSF variant even observed to be incompatible at any of cross-pollination treatments. Mature fruits were only observed in the reciprocal cross-pollination of YBF RBF variants. The population with higher heterozygosity and more outcrosser failed to produce seeds in hand-self pollination. However, the highly clonalized populations which tended to be more inbreeders were able to produce flowers in hand-self pollination.
Differences in floral structures, as well as the genetical differentiation, might have considerable effects on the success of reproductive processes. The hand-pollination experiments in this study clearly indicate that sandalwood is an outcrossing, allogamous and self-incompatible species. Hence, sexual reproductive failure was considered a result of inbreeding depression due to a very low genetic base and/or the strong clonality. Maintaining the evolutionary processes within each population should be carried out with a different strategy according to the genetic diversity, degree of clonality and the composition of floral variants in the population. Some of small, fragmented sandalwood plantations in Gunung Sewu are tending towards clonality; and the sub-structuring of the population was strengthened by the occurrence of three different variants that were failed to intermate each other. The fragmentation of these populations into small stands of almost identic genotypes, exhibiting little or no sexual reproduction, suggests a need for a broader conservation strategy. In order to maintain local adaptations that contribute to the adaptive diversity of the species, it is recommended to conduct enrichment planting and facilitate natural regeneration by sexual reproduction in each of floral variants. Restoring the population by replanting areas is using the same variants, but with the broader genetic base, and also recommended. Since then, variants were considered to be classified in different gene pool due to its difference on floral structures, allele distribution, and crossability; and supposed to be managed under a different management unit.
These results are part of long-term research on sandalwood genetic conservation and improvement; which is supported by the grant of Hibah Peningkatan Kapasitas Peneliti Dosen Muda year 2016 by the Universitas Gadjah Mada in the grant number. 1613/UN1-P.III/LT/DIT-LIT/2016. The laboratory work and analyses were done in the Genetics Laboratory of Silviculture Department, Universitas Gadjah Mada in Yogyakarta, Indonesia, and many thanks to Professor Mohammad Na'iem as the head of the laboratory. Many thanks go to Untung Maryanto for laboratory work, and Sunar (Bejiharjo), Andiyanto (Nglanggeran), Kustini (Bleberan), and Suparsono (Petir) for field assistances.