The spatial distribution of geographical distances at five natural populations of
Camellia japonica
in Busan, Korea was studied. The four plots (Mollundae, Gadeok-do, Du-do, and Jwiseum) of
C. japonica
were uniformly distributed in the forest community and only one plot (Amnam-dong) was aggregately distributed in the forest community. Morisita index is related to the patchiness index showed that the plot 20m×50m had an overly steep slope when the area was larger than 20m×20m, which indicated that the degree of aggregation increased significantly with increasing quadrat sizes, while the patchiness indices did not change from the plot 5m×10m to 10m×10m. The spatial structure was quantified by Moran’s
I
, a coefficient of spatial autocorrelation. Ten of the significant values (76.9%) were positive, indicating similarity among individuals in the first 4 distance classes (80 m), i.e., pairs of individuals with dissimilarity characteristics can separate by more than 100 m.
Introduction
In recent decades, there was much increase in the statistical tools used in spatial ecology
[13]
. Botanists, ecologists, geographers and plant evolutionary biologists have long recognized that plants are not distributed at random within communities but are rather clustered in distinct patches
[8
,
20]
. Environmental heterogeneity is usually cited as playing a critical role in determining the spatial structure, but colonization patterns and stochastic events affecting establishment and mortality are also important
[18]
. More recently, plant evolutionary biologists have demonstrated that genetic variations in plant populations are also nonrandomly distributed
[6]
. This nonrandom distribution of genetic variation is often referred to as the genetic structure of a population
[9]
. The genetic structure is an integral part of the process of population genetics
[5]
. Population structure interacts with a number of factors: microenvironmental heterogeneity, mortality due to stochastic events
[19]
, and mating systems that feature limited dispersal of seed or pollen
[5]
.
Of the several methods of describing the spatial distribution of plant community the simplest way is percentage distribution of individuals over the geographical areas. Another methodology usually adopted is to list the geographical areas of a given class into rank order which enables comparison of ranking from individual to individual. In this report, the several statistical tools of percentage distribution and population structure of the geographical areas are used to study the spatial distribution of
Camellia japonica
in Busan.
In theory, genetic differentiation over short distances may occur either as a result of spatially variable selection or localized genetic drift, provided that gene flow is sufficiently restricted
[14]
. Indirect evidence for genetic correlations between neighboring plants has been obtained from data on mating systems
[5]
. Localized seed and pollen dispersals produce family clusters within these populations
[7]
. Several studies have revealed decreased seed set and seed survivorship from mating between genetically similar near-neighbors, which has been interpreted as inbreeding depression
[10]
.
In the wild,
C. japonica
is found in mainland China (Shandong, east Zhejiang), Taiwan, southern Korea and southern Japan.
C. japonica
is sometimes called the rose of winter, it belongs to the Theaceae family. An edible oil is obtained from the seeds of this species. The leaves are a tea substitute. The dried flowers used as a vegetable or mixed with gelatinous-rice to make a Japanese food called ‘mochi’.
Saha-gu and Gangseo-gu locate in south of the Korean, and at the western part of the North Pacific Ocean. A sample of a large (more than 300 individuals) natural population of wild species
C. japonica
collected at Saha-gu and Gangseogu in Korea was used in this study.
The purpose of this paper was to describe a statistical analysis for detecting a species association, which is valid even when the assumption of within- species spatial randomness is violated. The purpose of this study is to find if there a spatial structure within four populations of
C. japonica
and 2) if so, what is the spatial pattern and if it is the same for all populations?
Materials and Methods
- Study area
We conducted the spatial analysis in the communities of
Camellia japonica
at Saha-gu and Gangseo-gu in Busan-si (
Fig. 1
). This area is on the southern margin of Busan. It has a temperate climate with a little hot and long summer. In this region the mean annual temperature is 14.7℃ with the maximum temperature being 29.4℃ in August and the minimum −0.6℃ in January. Mean annual precipitation is about 1519.1 mm with most rain falling period between June and August.
The five studied populations of Camellia japonica in Busan, Korea. ANM: Amnam-dong, DUD: Du=do, JWI: Jwiseum, GAD: Gadeok-do
- Sampling procedure
We established nine plots with an area of 20 m × 160 m each around four populations at Saha-gu and one population at Gangseo-gu in Busan, 2014. We randomly located quadrates in each plot which we established populations. The quadrat sizes were 5 m × 5 m, 5 m × 10 m, 10 m × 10 m, 10 m × 20 m, 20 m × 20 m, and 20 m × 50 m. We mapped all plants to estimate population density.
- Index calculation and data analysis
The spatial pattern of
C. japonica
was analyzed according to the Neatest Neighbor Rule
[3
,
11]
with Microsoft Excel 2010.
Average viewing distance (
rA
) was calculated as follows:
Where
ri
is the distance from the individual to its nearest neighbor.
N
is the total number of individuals within the quadrat.
The expectation value of mean distance of individuals within a quadrat (
rB
) was calculated as follows:
Where
D
is population density and
D
is the number of individuals per plot size.
R = rA / rB
When
R
>1, it is a uniform distribution,
R
=1, it is a random distribution,
R
<1, it is an aggregated distribution.
The significance index of the deviation of
R
that departs from the number of “1” is calculated from the following formula
[11]
.
When
CR
>1.96, the level of the significance index of the deviation of R is 5%, and when
CR
>2.58, the level is 1%.
We calculated the degree of population aggregation under different sizes of plots by dispersion indices: index of clumping or the index of dispersion (
C
), aggregation index (
CI
), mean crowding (
M
*), patchiness index (
PAI
), negative binominal distribution index K, Ca indicators (Ca is the name of one index)
[12
] and Morisita index (IM) were calculated with Microsoft Excel 2010. The formulae are as follows:
Index of dispersion: C = S2/m
Aggregation index
Mean crowding
= m+CI=m+C−1−1
Patchiness index
Aggregation intensity
Ca indicators Ca = 1/k
IM =
Where
S
2
is variance and
m
is mean density of
C. japonica.
When
C, M*, PAI
>1, it means aggregately distributed, when
C, M*, PAI
<1, it means uniformly distributed, when
CI, PA
, Ca >0, it means aggregately distributed, and when
CI, PA
, Ca <0 it means uniformly distributed.
We used the mean aggregation number to find the reason for the aggregation of
C. japonica
[1]
.
δ = mr/2k
Where r is the value of chi-square when the degree of freedom is 2k and k is the aggregation intensity.
- Spatial structure
Numerical simulations of previous analyses were performed to investigate the significant differences at various distance scales, i.e., 10.0, 20.0 m, and so on. However, no significant population structure was found within the 20.0 m distance classes by means of Moran's
I
, and a significant population structure was revealed beyond 20.0-m. Thus, the distance classes are 0-20.0 m (class I), 20.0-40.0 m (class II), 40.0-60.0 m (class III), 60.0-80.0 m (class IV), 80.0-100.0 m (class V), 100.0-120.0 m (class VI), 120.0-140.0 m (class VII), and 140.0-160.0 m (class VIII). The codes of classes are the same as in the distance classes and are listed
Table 1
.
Spatial patterns ofCamelliajaponica individuals at five populations in Busan
R and CR were shown in text.
The spatial structure was quantified by Moran's
I
, a coefficient of spatial autocorrelation (SA)
[15
,
16]
. As applied in this study, Moran's
I
quantifies the similarity of pairs of spatially adjacent individuals relative to the population sample as a whole. The value of
I
ranges between +1(completely positive autocorrelation, i.e., paired individuals have identical values) and −1(completely negative autocorrelation). Each plant was assigned a value depending on the presence or absence of a specific individual. If the
i
th plant was a homozygote for the individual of interest, the assigned pi value was 1, while if the individual was absent, the value 0 was assigned.
Pairs of sampled individuals were classified according to the Euclidian distance, dij, so that class k included dij satisfying k − 1 < dij < k + 1, where k ranges from 1 to 7. The interval for each distance class was 20 m. Moran's
I
statistic for class k was calculated as follows:
where Zi is
pi
− p (p is the average of
pi
); Wij is 1 if the distance between the
i
th and
j
th plants is classified into class k; otherwise, Wij is 0; n is the number of all samples and S is the sum of
in class k. Under the randomization hypothesis,
I
(k) has the expected value u1 = −1/(n − 1) for all k. Its variance, u2, has been given, for example, in Sokal and Oden (1978a). Thus, if an individual is randomly distributed for class k, the normalized
I
(k) for the standard normal deviation (SND) for the plant genotype, g(k) = {
I
(k) − u1}/u2
1/2
, asymptotically has a standard normal distribution
[3]
. Hence, SND g(k) values exceeding 1.96, 2.58, and 3.27 are significant at the probability levels of 0.05, 0.01, and 0.001, respectively.
Results
- The spatial pattern of individuals
Population densities (
D
) varied from 0.097 to 0.190, with a mean of 0.130 (
Table 1
). The values of spatial distance (the rete of observed distance-to-expected distance) among the nearest individuals were higher than 1 but Amnam-dong was lower than 1. The four plots (Mollundae, Gadeok-do, Du-do, and Jwiseum) of
C. japonica
were uniformly distributed in the forest community and only one plot (Amnam-dong) was aggregately distributed in the forest community (
Table 1
).
- The degree of population aggregation
Dispersion index (
C
) were higher than 1 except for three quadrats (5 m × 5 m, 5 m × 10 m, and 10 m × 10 m) of Amnam-dong (
Table 2
). As the sizes of quadrat were greater, the values of
C. japonica
were high. Thus aggregation indices were positive except for plots which indicate a clumped distribution. The values of
PI
and
Ca
were shown greater than zero. The values of
PAI
except three quadrats of Anman-dong were greater than 1 (
Table 2
). Thus, the most individuals of
C. japonica
were clustered and the distribution pattern of the
C. japonica
was quadrat-sampling dependent. When the sampling quadrat in Amnam-dong was smaller than 10 m × 10 m,
C
. japonica were aggregately distributed, and when the sampling quadrat was greater than 10 m x 10 m, the aggregation index showed the trend of being uniformly distributed for
C. japonica
. The mean crowding (
M
*) and aggregation intensity (
Ca
) indicator were higher with a big quadrat.
Changes in gathering strength ofCamellia japonicaat different sampling quadrat sizes
Aggregation indices were shown in text.
Morisita index (
IM
) related to the patchiness index (
PAI
) showed that the plot 20 m × 50 m had an overly steep slope when the area was larger than 20 m × 20 m, which indicated that the degree of aggregation increased significantly with increasing quadrat sizes, while the patchiness indices did not change from the plot 5 m × 10 m to 10 m × 10 m (
Fig. 2
).
The curves of patchiness in four populations of Camellia japonica using values of Green index.
The mean aggregation number analysis showed that the reasons for aggregation of
C. japonica
differed in quadrats with different plot sizes (
Table 3
). The cluster at 5 m × 5 m quadrat was determined by environmental factors. When the sizes were greater than 10 m × 10 m quadrat, the clusters were determined by both species characteristics and environmental factors.
Spatial autocorrelation coefficients (Moran'sI) among five populations ofCamellia japonicafor seven distance classes
* p<0.05, ** p<0.01, *** p<0.001.
The changes in the mean aggregation numbers for five populations.
- Analysis of spatial autocorrelation
The spatial autocoefficient, Moran's
I
is presented in
Table 3
. Separate counts for each type of joined individuals and for each distance class of separation were tested for significant deviation from random expectations by calculating the SND. Moran's
I
of
C. japonica
significantly differed from the expected value in only 13 of 35 cases (37.1%). Thirteen of these values (37.1%) were negative, indicating a partial dissimilarity among pairs of individuals in the seven distance classes. Ten of the significant values (76.9%) were positive, indicating similarity among individuals in the first four distance classes, i.e., pairs of individuals can separate by more than 100 m. Namely, significant aggregations were partially observed within IV classes. As a matter of course, the negative SND values at classes IV, V, and VI. Thus, dissimilarity among pairs of individuals could be found by more than 100 m.
The comparison of Moran’s
I
values to a logistic regression indicated that a highly significant percentage of individual dispersion in
C. cammelia
populations of the Saha-gu and Gangseo-gu could be explained by isolation by distance.
Discussion
When the value of δ is less than 2, the aggregation is mainly caused by the environmental factors
[11]
. When δ is higher than 2, the aggregation is mainly caused by both species characteristics and environmental factors
[11]
. We recognized that the important environmental factors might be considered competition, growth rate, little decomposition, light, and below-ground resources. The characteristics of the
C. japonica
included primarily their life history, artificial disturbance, and population density. Life history theory seeks to understand the variation in traits such as growth rate, number and size of offsprings and life span observed in nature, and to explain them as evolutionary adaptations to environmental conditions
[17]
. The cluster was determined by environmental factors when the sampling quadrat was smaller than 10 m × 10 m (<2). Artificial disturbance such as constitutional roads is an important environmental factor affecting
C. japonica
in Amnam-dong. At the plots which had fewer
C. japonica
, the cluster was mainly determined by
C. japonica
themselves; the mean value of the aggregation index changed irregularly with the variation in plot sizes.
A significant positive value of Moran's
I
indicated that pairs of individuals separated by distances that fell within distance class IV had similar individuals, whereas a significant negative value indicated that they had dissimilar individuals. The overall significance of individual correlograms was tested using Bonferroni's criteria. The results revealed that patchiness similarity was shared among individuals within up to a scale of an 80 m distance. Thus it was looked for the presence of dispersion correlations between neighbors at this scale.
The results from this study are consistent with the supposition that a plant population is subdivided into local demes, or neighborhoods of related individuals
[4
,
7]
. Previous reports on the local distribution of genetic variability suggested that microenvironmental selection and limited gene flow are the main factors causing substructuring of alleles within a population
[5]
.
The community of
C. japonica
at Gadeok-do is the Busan Natural Monument No. 38. There are some communities of the Korean Natural Monuments of
C. japonica
: No. 66 (Daecheon-do, Ungjin-gun, Gyeonggi-do). No. 161 (Baekryeonsa, Gangjin-gun, Jeon-nam), No. 169 (Marhang-ri, Seocheon-gun, Chung-nam), No. 184 (Seonsa, Gochang-gun, Jeon-buk), No. 489 (Okrhongsa, Gwangyang-ci, Jeon-nam).
C. japonica
is one of very important resource of East Asia. The camellia in Europe was brought details of over 30 varieties back from Asia. Camellias were introduced into Europe during the 18th century and had already been cultivated in the Orient for thousands of years.
In conclusion,
C. japonica
populations within the Saha-gu and Gangseo-gu were observed a strong spatial structure. Neighboring patches of
C. japonica
were predominantly 80 to 100 m apart on average. The present study demonstrates that a spatial structure of
C. japonica
in the Saha-gu and Gangseo-gu populations could be explained by isolation by distance, limited gene flow, and topography. The results of this study were used as systematic conservation planning which is an effective way to seek and identify efficient and effective types of reserve design to capture or sustain the highest priority biodiversity values and to work with communities in support of local ecosystems.
Acknowledgements
This work was supported by Dong-eui University Grant (2014AA474).
Arbous A. G.
,
Kerrich J. E
1951
Accident statistics and the concept of accident proneness
Biometrics
http://dx.doi.org/10.2307/3001656
7
340 -
342
DOI : 10.2307/3001656
Clark P. J.
,
Evans F. C
1954
Distance to nearest neighbor as a measure of spatial relationships in populations
Ecology
http://dx.doi.org/10.2307/1931034
35
445 -
453
DOI : 10.2307/1931034
Cliff A. D.
,
Ord J. K.
1971
Spatial autocorrelation
Pion
London
Epperson B. K.
,
Chung M. G.
,
Telewski F. W.
2003
Spatial pattern variation in a contact zone of Pinus ponderosa and P. arizonica (Pinaceae)
Am J Bot
http://dx.doi.org/10.3732/ajb.90.1.25
90
25 -
31
DOI : 10.3732/ajb.90.1.25
Evans K. L.
,
Duncan R. P.
,
Blackburn T. M.
,
Crick H. Q. P
2005
Investigating geographic variation in clutch size using a natural experiment
Functional Ecol
http://dx.doi.org/10.1111/j.1365-2435.2005.01016.x
19
616 -
624
DOI : 10.1111/j.1365-2435.2005.01016.x
Garnier L. K. M.
,
Durand J.
,
Dajoz I.
2002
Limited seed dispersal and microspatial population structure of an agamospermous grass of west African savannahs, Hyparrhenia diplandra (Poaceae)
Am J Bot
http://dx.doi.org/10.3732/ajb.89.11.1785
89
1785 -
1791
DOI : 10.3732/ajb.89.11.1785
Hamrick J. L.
,
Godt M. J. W.
,
Sherman-Broyles S. L
1992
Factors influencing levels of genetic diversity in woody plant species
New Forests
http://dx.doi.org/10.1007/BF00120641
6
95 -
124
DOI : 10.1007/BF00120641
Hayakawa T.
,
Tomaru N.
,
Yamamoto S
2004
Stem distribution and clonal structure of Chamaecyparis pisifera growing in an old-growth beech-conifer forest
Ecol Res
http://dx.doi.org/10.1111/j.1440-1703.2004.00652.x
19
411 -
420
DOI : 10.1111/j.1440-1703.2004.00652.x
Ishida T. A.
,
Kimura M. T.
2003
Assessment of within-population genetic structure in Quercus crispula and Q. dentate by amplified fragment length polymorphism analysis
Ecol Res
http://dx.doi.org/10.1046/j.1440-1703.2003.00583.x
18
619 -
623
DOI : 10.1046/j.1440-1703.2003.00583.x
Lian X.
,
Jiang Z.
,
Ping X.
,
Tang S.
,
Bi J.
,
Li C
2012
Spatial distribution pattern of the steppe toad-headed lizard (Phrynocephalus frontalis) and its influencing factors
Asian Herpet Res
http://dx.doi.org/10.3724/SP.J.1245.2012.00046
3
46 -
51
DOI : 10.3724/SP.J.1245.2012.00046
Lloyd M
1967
Mean crowding
J Anim Ecol
http://dx.doi.org/10.2307/3012
36
1 -
30
DOI : 10.2307/3012
Reiczigel J.
,
Lang Z.
,
Rozsa L.
,
Tothmeresz B
2005
Properties of crowding indices and statistical tools to analyze parasite crowding data
J Parasitol
http://dx.doi.org/10.1645/GE-281R1
91
245 -
252
DOI : 10.1645/GE-281R1
Schoen D. J.
,
Latta R. G.
1989
Spatial autocorrelation of genotypes in populations of Impatiens pallida and Impatiens capensis
Heredity
http://dx.doi.org/10.1038/hdy.1989.90
63
181 -
189
DOI : 10.1038/hdy.1989.90
Sokal R. R.
,
Oden N. L.
1978
Spatial autocorrelation in biology 1. methodology
Biol J Linn Soc
http://dx.doi.org/10.1111/j.1095-8312.1978.tb00013.x
10
199 -
228
DOI : 10.1111/j.1095-8312.1978.tb00013.x
Sokal R. R.
,
Oden N. L.
1978
Spatial autocorrelation in biology 2. some biological implications and four applications of evolutionary and ecological interest
Biol Linn Soc
http://dx.doi.org/10.1111/j.1095-8312.1978.tb00014.x
10
229 -
249
DOI : 10.1111/j.1095-8312.1978.tb00014.x
Souza A. F.
,
Martins F. R.
2004
Microsite specialzation and spatial distribution of Geonoma brevispatha, a clonal palm in south-eastern Brazil
Ecol Res
http://dx.doi.org/10.1111/j.1440-1703.2004.00670.x
19
521 -
532
DOI : 10.1111/j.1440-1703.2004.00670.x
Stearns S.
1992
The Evolution of Life Histories
Oxford Univ. Press
UK
Walter R.
,
Epperson B. K
2004
Microsatatellite analysis of structure among seedling in populations of Pinus strobus (Pinaceae)
Am J Bot
http://dx.doi.org/10.3732/ajb.91.4.549
91
549 -
557
DOI : 10.3732/ajb.91.4.549
Woodward F. I.
1987
Climate and plant distribution
Cambridge Univ. Press
Cambridge, UK
Xia B.
,
Abbott D
1987
Edible seaweeds of China and their place in the Chinese diet
Econ Bot
http://dx.doi.org/10.1007/BF02859049
41
341 -
353
DOI : 10.1007/BF02859049