- Research article
- Open Access
Determination of the genetic structure of remnant Morus boninensis Koidz. trees to establish a conservation program on the Bonin Islands, Japan
© Tani et al; licensee BioMed Central Ltd. 2006
- Received: 29 May 2006
- Accepted: 11 October 2006
- Published: 11 October 2006
Morus boninensis, is an endemic plant of the Bonin (Ogasawara) Islands of Japan and is categorized as "critically endangered" in the Japanese red data book. However, little information is available about its ecological, evolutionary and genetic status, despite the urgent need for guidelines for the conservation of the species. Therefore, we adopted Moritz's MU concept, based on the species' current genetic structure, to define management units and to select mother tree candidates for seed orchards.
Nearly all individuals of the species were genotyped on the basis of seven microsatellite markers. Genetic diversity levels in putative natural populations were higher than in putative man-made populations with the exception of those on Otouto-jima Island. This is because a limited number of maternal trees are likely to have been used for seed collection to establish the man-made populations. A model-based clustering analysis clearly distinguished individuals into nine clusters, with a large difference in genetic composition between the population on Otouto-jima Island, the putative natural populations and the putative man-made populations. The Otouto-jima population appeared to be genetically differentiated from the others; a finding that was also supported by pairwise FST and RST analysis. Although multiple clusters were detected in the putative man-made populations, the pattern of genetic diversity was monotonous in comparison to the natural populations.
The genotyping by microsatellite markers revealed strong genetic structures. Typically, artificial propagation of this species has ignored the genetic structure, relying only on seeds from Otouto-jima for replanting on other islands, because of a problem with inter-specific hybridization on Chichi-jima and Haha-jima Islands. However, this study demonstrates that we should be taking into consideration the genetic structure of the species when designing a propagation program for the conservation of this species.
- Seed Orchard
- Genetic Composition
- Maternal Tree
- Operational Population
- Pairwise Kinship
Morus boninensis, a plant native to the Bonin Islands (typical oceanic islands, located 1,000 km south of Tokyo, Japan), is only endemic to Otouto-jima, Chichi-jima and Haha-jima Islands; it is categorized as "critically endangered" in the Japanese Red Data Book . This species is a typical case in which there is little information about the species, although guidelines are urgently needed to aid in its conservation. There are fewer than about 170 remaining trees and natural regeneration does not seem to be occurring at present (Yoshimaru et al. unpublished data). The reason for the degradation of the species was intensive logging during the last quarter of the 19th century and the start of the 20th century (details described in ). Although Morus boninensis used to be one of the main species constituting the canopy in the moist tall forest on the Bonin Islands, some invasive trees, mainly Bischofia javanica, have replaced it in recent years [3–5]. In our field observations, seedling recruitment has not been observed since 1995. Yoshimaru et al. (unpublished data) estimated that the mortality rate of the mature individuals is between 0.56% and 3.56% per year in each population. Furthermore, hybridization with the introduced species, M. acidosa, has been observed and has been confirmed by molecular marker analysis . To promote the propagation of the next generation, selection of mother trees should be considered to maximize evolutionary success based on the concept of the Evolutionary Significant Unit (ESU, ). To achieve this, it is best practice to define ESUs based on genetic as well as ecological information. However, there is no ecological information about the species. Furthermore, the Bonin Islands are a typical example of the changing balance in Japan between bio-diversity and single-minded development, between the desire to conserve native species and the desire to satisfy human desires, and between the modesty and creativity of local peoples and the arrogance and insensitivity inherent in massive public works funding. Therefore, it is urgent that guidelines for conducting ex situ conservation and promoting the propagation of individuals for the next generation are put in place. One proposal by Moritz  was that the population ESU should be defined by the reciprocal monophyletic relationship based on mtDNA alleles and significant divergence of allele frequencies at nuclear loci (Moritz's Management unit, MU). Although Crandall et al  identified several conceptual and practical problems with the effectiveness of the use of a historical population structure, as defined by molecular genetic techniques, the concept has been adopted in various applied studies of animals to define conservation units based on ESUs [9–11]. Because of the pressing nature of our work, we have adopted Moritz's MU concept to define management units and aid in the selection of mother tree candidates for the seed orchards. This is based on the current genetic structure, since only genetic information is available at present. In this paper, we present a description of the current genetic structure of the species, genetic differentiation between populations and kinship within clustered individuals based on microsatellite markers. These data can be used to establish a conservation program for the species.
Genetic variation within the operational populations
Information of sampling sites and trees, such as background of population, location, tree size and vegetation
Background of the population
Latitude and Longitude
Tree size (DBH)
Various type of vegetation
Introduced species dominating
Average genetic variation of each locus for five operational populations of Morus boninensis. Standard deviation of the mean is in parenthesis. Clusters 2, 6 and 7 are groups assigned by the model-based clustering analysis with a 90% confidence interval
Fixation indices for five operational populations of Morus boninensis, with and without double reduction. Clusters 2, 6 and 7 are groups assigned by the model-based clustering analysis with a 90% confidence interval
F (Chromosome, Ce)
F (Chromatid, Cd)
Genetic differentiation among the operational populations
The global F and R statistics for seven microsatellite loci
Individual based analysis using the model-based clustering method
Genetic diversity and kinship structure within the clusters
Genetic diversity within populations
The model-based clustering method demonstrated that the genetic structure of the M. boninensis remnant trees did not correspond perfectly to the operational populations based on geographical distribution (Fig. 4). Although the OTJ operational population did coincide with the model based clustering, other operational populations contained multiple clusters. However, in determining operational populations to be used in species conservation, there are two considerations of practical importance: 1) the convenience of conducting controlled crossing using trees in close physical proximity, whilst knowing their individual genetic relatedness; and 2) the ability to sample seed that may have been open pollinated from nearby paternal trees. Therefore, it is important to understand the genetic composition at two levels: the operational population and the individual level. The genetic diversity levels in OTJ and CCJ were lower than in the other populations (Table 2). Although the OTJ population contains the secondly largest number of trees of all the operational populations, the genetic composition of each tree is very similar, as demonstrated by the model-based clustering (Fig. 4). This might be the result of a founding event, with a small number of maternal trees having produced the seed to establish the extant OTJ population. However, the low genetic diversity in the CCJ population is different. Four different genetic groups were identified in the CCJ population by model-based clustering analysis. Out of the four genetic groups, some of the trees comprising three of them (blue/red, yellow and green in Fig 4) are located in the northern part of Chichi-jima Island. This is an area inhabited by immigrant people. Trees in these same genetic groups were also present in KWK on Haha-jima Island. Therefore, these trees in CCJ were probably planted by people using seed sources from Haha-jima Island. In contrast, the vegetation of the south of Chichi-jima is, to date, well preserved although there was some disturbance during World War II. These trees have retained their endemic genetic composition (purple in Fig. 4) on Chichi-jima Island, which probably represents the native genetic composition of the species here. We only found nine remnant trees classified with the purple genetic composition in the CCJ population. Although the expected heterozygosity level of the nine remnant trees was smaller than that of CCJ (0.628 vs. 0.674), their allelic richness exceeded that of CCJ (2.476 vs. 2.429). This might be because a limited number of maternal trees were used as the seed source for planting the trees in the northern part of Chichi-jima Island, or because of bi-parental inbreeding of parent trees of this population.
On Haha-jima Island, we observed a similar pattern of genetic diversity as for the CCJ population. Trees on Haha-jima Island were divided into four operational populations. One of the current authors (Y. Hoshi) remembered that the species had been grown in plantations on this Island. He recollected that part of KWK and LPS were planted by humans. Although the origin of the IGM population was not clear, the small, relatively uniform trees in IGM are distributed in quite a small area. This might indicate that, whilst these trees were not planted, they are the result of simultaneous regeneration from a few maternal trees, such as those of an original plantation. The amounts of within population genetic diversity of plantation trees (LPS) and IGM were lower than those of SKM and KWK. According to the model-based clustering, we identified a similar pattern of genetic composition between the planted populations (part of KWK and LPS) and IGM – the clusters represented by mostly green and pink colours (Fig. 4). This might also indicate an origin from a limited number of maternal trees or bi-parental inbreeding. In contrast, the rest of the trees in the KWK and SKM populations exhibited a complicated genetic composition, and the SKM population contained the highest amount of genetic diversity of all the operational populations, despite having the smallest sample size. In terms of genetic diversity and size of tree (data not shown), most trees in the SKM population and part of the KWK population are likely to be natural remnants.
The evidence of population bottlenecks and inbreeding
The individuals in the SKM population exhibited many components of different clusters (Fig. 4). As a result of its genetic characteristics, size of trees and vegetation pattern on the Islands, we deduced that this population is the closest to being natural. When we compared the genetic diversity level between the OTJ (almost equivalent to cluster 6) and SKM populations, the number of alleles, as well as the observed and expected heterozygosity was lower in the former. According to one theoretical study, the amount of reduction in heterozygosity depends not only on the "bottleneck" size but also on the rate of population growth after passing through the bottleneck, while the loss of alleles largely depends on the size of the "bottleneck" . If the SKM population produced many generations after the bottleneck event, the heterozygosity level would be reduced. However, we have assumed that the SKM population produced few generations after population decline, because M. boninensis is known to be an extremely long-lived woody species . Although the SKM population has maintained a high level of genetic diversity to date, the decline in genetic diversity in OTJ might be the result of bottleneck events, such as much severe logging, long-term population decline, or catastrophic population destruction.
When we estimated the pairwise kinships between individuals within a cluster and compared their average between clusters, the difference in average pairwise kinships between cluster 6 and all individuals was significant (Fig. 5). This meant that bi-parental inbreeding may have occurred during the bottleneck event in the OTJ population. However, the level of inbreeding was, supposedly, not strong, although expression of strong inbreeding depression for the species may have been masked by mass mortality of inbred individuals at the seedling stage. This speculation is also supported by the non-significant and negative fixation indices in cluster 6 and OTJ (Table 3). In contrast, potentially man-made populations, such as parts of KWK and LPS, and IGM, contained multiple clusters. However, most individuals within these populations were classified into two main clusters (green and pink), supported by the 90% confidence interval in the model-based clustering analysis (Fig. 4). Furthermore, we detected higher pairwise kinship values within clusters (green and pink) than those for all pairs of individuals (Fig. 5). This could be the consequence of a limited number of maternal trees producing the seeds from which the man-made populations were propagated. In contrast, although the mean of the pairwise kinship of SKM was somewhat high value, the mean of pairwise kinship of the putative natural populations and mixed background population, SKM and CCJ, is not significantly different from zero. This suggests that these populations might not have experienced intensive bi-parental inbreeding, as occurred in the putative man-made populations. If we collect seed from a few mother trees for the propagation of the next generation, even open pollination will potentially cause bi-parental inbreeding and a decline in the genetic diversity of the new generation. Therefore, it is very important to construct seed orchards using clones of trees from natural populations and/or to perform controlled crossing between them. Furthermore, some pairs in SKM population represented high level of pairwise kinships between them, which implied that we must consider selection of parental trees to construct seed orchards and combinations of parental trees for control crossing even if parental trees are only selected from the natural populations
Establishment of conservation units and a propagation program
Ex situ conservation would seem to be necessary to conserve M. boninensis, because it is difficult to obtain seeds free of hybridization with M. accidosa. The exception is the OTJ population . In addition, seedling establishment is extremely rare because of the competition with introduced alien species, such as Bischofia javanica . Therefore, selection of pure M. boninensis seedlings, controlled crossing and ex situ conservation are all necessary to ensure propagation of the species. To achieve this, we must consider Moritz's management units in order to select appropriate mother trees and we must undertake controlled crossing to obtain seed sources for propagation. In terms of genetic diversity, there was a large difference between the putative natural populations, with the exception of OTJ and IGM (south of CCJ, a part of KWK, most of SKM) and putative plantation trees. According to the model-based clustering analysis, planted individuals are classified into two main clusters, green and pink. This genetic differentiation, which is supported by the model-based clustering, should be used in selecting the MUs. We were able to identify at least seven main MUs for M. boninensis, which were OTJ, natural trees in CCJ, natural trees in KWK, SKM, man-made trees in CCJ and KWK and two major clusters including individuals in IGM and LPS (Green and Pink). Natural populations retained many elements of the different clusters. On the other hand, most of the plantation trees were simply classified into a single cluster. Hence, natural individuals were prime candidates for use as maternal trees in controlled crossing and ex situ conservation by grafting. According to analyses of pairwise FST and RST and the model-based clustering, we should avoid controlled crossing and instead should establish new seed orchards by open pollination of individuals from different clusters. To date, the OTJ population has been the only one to produce pure M. boninensis seeds. These seeds, therefore, have been the only ones used for propagation. There is a hybridization problem with the introduced species in the KWK and SKM populations. Because the MU of OTJ is definitely different from the other natural populations, we need to independently conduct controlled crosses or ex situ conservation (and subsequently develop seed orchards) within each cluster. Only then will propagated seedlings be available to use for propagation of the species in the forest according to the MUs identified in this study.
Genotyping using microsatellite markers revealed that the pattern of genetic variation was different between the OTJ population, the north CCJ putative man-made population, the south CCJ putative natural population, the natural population in Haha-jima (most of SKM and a part of KWK), IGM and the man-made population in Haha-jima (a part of KWK and LPS). These differences should guide the selection of the MUs. The putative natural populations exhibited higher genetic diversity than those of the man-made populations, IGM and OTJ. These are, therefore, important genetic resources for the propagation of a new generation. To date, propagation of the species has ignored its genetic structure by relying only on seeds from the OTJ population, even for planting on the other islands. This has been because of the inter-specific hybridization problem on Chichi-jima and Haha-jima Islands. However, this study has demonstrated the high level of genetic differentiation within the species. This needs to be considered with respect to any propagation program aimed at conserving the evolutionary range of the species.
Field survey and collection of samples
Between 1998 and 2002, we conducted a field survey to find the unknown remnant trees of M. boninensis. Summary of field survey concerning with background of each population, location, tree size and vegetation was shown in Table 1. Identification of true M. boninensis or hybrids between M. boninensis and M. acidosa was conducted, following sampling, using the three SCAR markers . This allowed us to exclude hybrids from subsequent population genetic analyses. To date, we have found and collected petioles from 35 trees on Otouto-jima Island (OTJ), 20 trees on Chichi-jima Island (CCJ) and 109 trees on Haha-jima Island (Yoshimaru et al. unpublished data). Four areas were recognized in terms of their geography and the distributional density of the remnant trees on Haha-jima Island. Therefore, we assumed that these four geographical areas would be operational populations. We gave them the names: Iguma-wan (IGM), Kuwanoki-yama (KWK), Sekimon (SKM) and lower plateau of Sekimon (LPS). We were able to recognize some sub-groups of trees within two operational populations: KWK and SKM. KWK contained Uchu-sawa (UCU) and Kuwanoki-yama test field (KWT). SKM contained the upper plateau of Sekimon (UPS), Kiri-hama (KRH) and Naga-hama (NGH). However, we did not deal with these tree clusters as operational populations because they contained relatively few individuals. Thus, 23, 23, 18 and 45 remnant trees were identified and sampled in the IGM, KWK, SKM and LPS operational populations, respectively. The petioles from these trees were stored at -20°C until the DNA was extracted.
DNA preparation and genotyping of microsatellite loci
DNA was extracted from 100 mg of the petioles of each of the sampled remnant trees using a DNeasy plant-mini Kit (QIAGEN, Hilden, Germany) and following the manufacturer's instructions. The genomic DNA concentration was adjusted to 5 ng/μL by dilution after measuring the concentration of extracted DNA using a spectrophotometer (Amersham Biosciences, Little Chalfont, UK). For all trees, the genotypes of seven microsatellite markers were assigned according to the procedure described by Tani et al. . Because M. boninensis is a putative autotetraploid species, we determined the allele copy number of partial heterozygotes based on dosage of electrophoretogram peaks.
Operational populations based analyses
We based the composition of the six operational populations on geographical distribution and individual tree aggregations. For these we calculated population genetic statistics as follows: the number of alleles (A); the average number of alleles per individual (Ai); the average number of four allele genotypes (G); and observed heterozygosity (Ho). These were obtained for each locus using the program AUTOTET . Two types of expected heterozygosities (He) and fixation indices (Fis) were calculated under two assumptions: 1) random mating and random chromosome segregation (RceS); and 2) random mating and some level of chromatid segregation (RcdS), for which we selected the maximum theoretical double reduction rate, α = 1/7 [17–19]. Using SPAGeDi software, F-statistics and R-statistics, based on the infinite allele model and the stepwise mutation model, respectively, were estimated in order to evaluate genetic diversity between the operational populations [20–22]. Random resampling of individuals was parmuteted at 20,000 times to obtain a confidence interval for the estimators of population differentiation and the inbreeding coefficient. The genetic relationships between the operational populations were elucidated using neighbor-joining trees based on pairwise FST and RST, with the aid of two software packages, Mega ver. 2.1 and SPAGeDi [20, 23, 24].
Individual based analyses
To identify any unknown population structure, we used a model-based clustering method implemented by the program Structure ver. 2.0. This estimated the number (K) of clusters into which the sample data (X) were fitted with posterior probability Pr(X|K), using a model with admixture and uncorrelated allele frequency . We conducted 106 iterations following a burn-in period of at least 30,000 iterations. The genotypes of the Mos0157-2 locus were omitted from the data set because Mos0157-2 was constructed from same sequence as Mos0157-1, and hence these loci were closely linked to each other. K provides only a rough guide for determining which models are consistent with the data. Therefore, we examined various K values, from 1 to 20, within the simulation, and searched for a suitable K value to maximise the posterior probability Pr(X|K). The Structure software estimates the proportion of ancestry (Q value) from each of the K clusters for each individual and the 90% probability intervals. Assignment of individuals into inferred populations was conducted using the Q values and their probability intervals. Only individuals supported by 90% probability intervals were assigned to each cluster. A kinship coefficient between individuals within a cluster was estimated using SPAGeDi software [13, 20]. The amount of kinship within each cluster was estimated by comparing the probability of two individuals in a cluster having identical genes to the probability that two individuals chosen at random from all the samples are identical.
The authors would like to thank N. Savory Jr. and N. Tanaka for information about remnant M. boninensis trees, I. Okochi for management of the project, C. Yamanaka and Y. Taguchi for laboratory assistance, H. Iwata for advice on data analysis and Y. Tsumura, S. L. Lee and two anonymous reviewers for valuable comments on the manuscript. This research was partly supported by a grant to finance a 'Restoration study of the Ogasawara forest ecosystem due to exclusion of the effects of naturalized organisms' from the Japanese Ministry of the Environment.
- Environment_Agency_of_Japan: Threatened Wildlife of Japan - Red Data Book 2nd ed.-. 2000, Tokyo , Japan Wildlife Research Center, 8, Vascular Plants:Google Scholar
- Tani N, Kawahara T, Yoshimaru H, Hoshi Y: Development of SCAR markers distinguishing pure seedlings of the endangered species Morus boninensis from M. boninensis x M. acidosa hybrids for conservation in Bonin (Ogasawara) Islands. Conserv Genet. 2003, 4 (5): 605-612. 10.1023/A:1025655331429.View ArticleGoogle Scholar
- Shimizu Y: Vegetation of Mt. Kuwanoki in the Bonin (Ogasawara) Islands with reference to the invasion of an introduced tree species Bischofia javanica. Instutute for Applied Geography, Komazawa University, Tokyo, Regional Views. 1988, 1: 31-46.Google Scholar
- Tanimoto T, Toyota T: Survivorship and growth of Akagi (Bischofia javanica BI) seedlings under the forest canopy and different temperature conditions. Bulletin of the Forestry and Forest Products Research Institute. 1996, 370: 1-19.Google Scholar
- Yamashita N, Tanaka N, Hoshi Y, Kushima H, Kamo K: Seed and seedling demography of invasive and native trees of subtropical Pacific islands. J Veg Sci. 2003, 14 (1): 15-24. 10.1658/1100-9233(2003)014[0015:SASDOI]2.0.CO;2.View ArticleGoogle Scholar
- Crandall KA, Bininda-Emonds ORP, Mace GM, Wayne RK: Considering evolutionary processes in conservation biology. Trends Ecol Evol. 2000, 15 (7): 290-295. 10.1016/S0169-5347(00)01876-0.View ArticlePubMedGoogle Scholar
- Guo N, McCormack G: Coming to terms with nature: development dilemmas on the Ogasawara Islands. Japan Forum. 2001, 13 (2): 177-193. 10.1080/09555800120081385.View ArticleGoogle Scholar
- Moritz C: Defining Evolutionarily-Significant-Units for Conservation. Trends Ecol Evol. 1994, 9 (10): 373-375. 10.1016/0169-5347(94)90057-4.View ArticlePubMedGoogle Scholar
- Douglas MR, Brunner PC: Biodiversity of Central Alpine Coregonus (Salmoniformes): Impact of one-hundred years of management. Ecol Appl. 2002, 12 (1): 154-172.View ArticleGoogle Scholar
- Hedrick PW, Parker KM, Lee RN: Using microsatellite and MHC variation to identify species, ESUs, and MUs in the endangered Sonoran topminnow. Mol Ecol. 2001, 10 (6): 1399-1412. 10.1046/j.1365-294X.2001.01289.x.View ArticlePubMedGoogle Scholar
- Holland BS, Hadfield MG: Islands within an island: phylogeography and conservation genetics of the endangered Hawaiian tree snail Achatinella mustelina. Mol Ecol. 2002, 11 (3): 365-375. 10.1046/j.1365-294X.2002.01464.x.View ArticlePubMedGoogle Scholar
- Yoshida K, Oka S: Ecological characteristics of Morus boninensis reconstructed from its remaining stumps in the Sekimon Region of Haha-jima Island, Ogasawara (Bonin) Islands, northwestern Pacific. Jpn J Histor Bot. 2000, 9 (1): 21-28.Google Scholar
- Loiselle BA, Sork VL, Nason J, Graham C: Spatial Genetic-Structure of a Tropical Understory Shrub, Psychotria Officinalis (Rubiaceae). Am J Bot. 1995, 82 (11): 1420-1425. 10.2307/2445869.View ArticleGoogle Scholar
- Nei M, Maruyama T, Chakraborty R: Bottleneck Effect and Genetic-Variability in Populations. Evolution. 1975, 29 (1): 1-10. 10.2307/2407137.View ArticleGoogle Scholar
- Tani N, Kawahara T, Yoshimaru H: Development and diversity of microsatellite markers for endangered species, Morus boninensis Koidz., to establish conservation program. Mol Ecol Notes. 2005, 5 (2): 398-400. 10.1111/j.1471-8286.2005.00940.x.View ArticleGoogle Scholar
- Thrall PH, Young A: AUTOTET: A program for analysis of autotetraploid genotypic data. J Hered. 2000, 91 (4): 348-349. 10.1093/jhered/91.4.348.View ArticlePubMedGoogle Scholar
- Bever JD, Felber F: The theoretical population genetics of autopolyploidy. Oxford surveys in evolutionary biology. Edited by: Antonovics J, Futuyma D. 1992, New York , Oxford University Press, 185-217.Google Scholar
- Geiringer H: Chromatid Segregation of Tetraploids and Hexaploids. Genetics. 1949, 34 (6): 665-684.PubMed CentralPubMedGoogle Scholar
- Haldane JBS: Theoretical genetics of autopolyploids. J Genet. 1930, 22: 359-372.View ArticleGoogle Scholar
- Hardy OJ, Vekemans X: SPAGEDi: a versatile computer program to analyse spatial genetic structure at the individual or population levels. Mol Ecol Notes. 2002, 2 (4): 618-620. 10.1046/j.1471-8286.2002.00305.x.View ArticleGoogle Scholar
- Rousset F: Equilibrium values of measures of population subdivision for stepwise mutation processes. Genetics. 1996, 142 (4): 1357-1362.PubMed CentralPubMedGoogle Scholar
- Slatkin M: A Measure of Population Subdivision Based on Microsatellite Allele Frequencies. Genetics. 1995, 139 (1): 457-462.PubMed CentralPubMedGoogle Scholar
- Kumar S, Tamura K, Jakobsen IB, Nei M: MEGA2: molecular evolutionary genetics analysis software. Bioinformatics. 2001, 17 (12): 1244-1245. 10.1093/bioinformatics/17.12.1244.View ArticlePubMedGoogle Scholar
- Saitou N, Nei M: The Neighbor-Joining Method - a New Method for Reconstructing Phylogenetic Trees. Mol Biol Evol. 1987, 4 (4): 406-425.PubMedGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155 (2): 945-959.PubMed CentralPubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.