We ran spatial auto-correlation analysis in SPAGEDI (Hardy and Vekemans, 2002 (
link)) to test for the presence of fine-scale genetic structure at the five populations with more than 20 individuals (WX1, WX2, WX3, NPA2 and PXH1). Pairwise kinship coefficients were calculated between all individuals (
Fij) within each population (Loiselle et al., 1995 (
link); Kalisz et al., 2001 (
link)), mean
Fij was derived for each distance interval,
d, and this was plotted against distance in metres. The software requires that the number of pairwise comparisons is kept constant across all distance intervals. Mean
Fij(
d) estimates were calculated for intervals defined as 0–20 m (
d = 5 m), 21–50 m (
d = 10 m), 51–100 m (
d = 50 m or end-point) and 101–600 m (
d = 100 m or end-point), and 95% confidence intervals (CI) associated with the null hypothesis of no genetic structure [
Fij(
d) = 0] were constructed using 1,000 random permutations. Significant positive or negative structure was inferred if the CIs did not overlap.
We then regressed the slope
bLF(d) [linear regression of
Fij (
d) on ln (
d)] to test whether there was significant deviation from the null hypothesis of no genetic structure [
bLF(d) = 0]. To compare overall intensity of fine-scale genetic structure among populations, we also calculated the
Sp statistic (Vekemans and Hardy, 2004 (
link)), given by
Sp = -
bLF(d)/[1-F (
d1)], where F(
d1) is the average kinship coefficient between individuals of the first distance class (i.e., 0–20 m,
d = 5 m),
Fij.
Finally, we estimated the relative contribution of pollen (σ
p) and seed (σ
s) dispersal to total gene flow, σ (Heuertz, 2010 (
link)). Using the average
Fij(
d) for all samples from each population, we regressed the residuals [
f(
d):
Fij(
d) -
Fij(
d)
exp] on ln(
d) by a polynomial regression of the third power:
f(
d) = a + b ln(
d) + c [ln(
d)]
2 + d [ln(
d)]
3, where
Fij(
d)
exp is the dependent variable of the linear regression equation at independent variable ln(
d). The curvature of
f(
d) is given by the second derivative,
k = 2c + 6d*ln (
d1), where
d1 is the average distance of the first distance class. A concave curve at short distances or
k >0 suggests more restricted seed dispersal than pollen dispersal (σ
s ≪ σ
p), whereas a convex shape or
k <0 suggests more restricted pollen dispersal or no particular restriction in seed dispersal (σ
s ≥ σ
p) (Vekemans and Hardy, 2004 (
link)). Statistics were calculated in SPAGEDI (Hardy and Vekemans, 2002 (
link)) and SPSS 22.0 (IBM Corp., New York, USA).
Zhang Z., Li J., Suddee S., Bouamanivong S., Averyanov L.V, & Gale S.W. (2023). Exploring island syndromes: Variable matrix permeability in Phalaenopsis pulcherrima (Orchidaceae), a specialist lithophyte of tropical Asian inselbergs. Frontiers in Plant Science, 14, 1097113.