Search
2023 Volume 2
Article Contents
ARTICLE   Open Access    

Effects of habitat fragmentation on the coastal Vatica mangachapoi forest (Dipterocarpaceae) in Shimei Bay, Hainan Island, China

More Information
  • There is no significant difference in genetic diversity between the coastal and lowland V. mangachapoi populations.

    No fine-scale spatial genetic structure was detected in the coastal SM population.

    Strong gene flow between populations accounts for the homogenous genetic structure of the coastal and lowland populations.

    There is not enough time to accumulate differentiation between fragmented populations in the coastal forest.

  • Habitat fragmentation can cause isolation and decline of a formerly continuously distributed population, which leads to loss of genetic variation and increased risk of extinction. Vatica mangachapoi Blanco is a dominant tree species growing in the lowland rainforests of Hainan Island, China. Remarkably, this species dominates a coastal forest in Shimei Bay, Wanning City of Hainan Province (China). Due to logging, expansion of farmland and villages, and construction of tourism facilities, the coastal V. mangachapoi-dominated forest has become fragmented, threatening its future. To evaluate the effects of habitat fragmentation on this unique coastal forest, two V. mangachapoi populations (SM and RY) along the coast and one population in the lowland rainforest near the coast were selected, and their genetic diversity was assessed based on 12 SSR markers. In addition, the genetic structure of the three populations and gene flow among them, and the fine-scale spatial genetic structure (FSGS) of the SM population were also studied. The results show that the three V. mangachapoi populations had comparable levels of genetic variation, and differentiation among them is negligible ($ {F}_{\mathrm{s}\mathrm{t}} $ = 0.008 ~ 0.013). Model-based clustering, Principal co-ordinate analysis and the Neighbor-joining (NJ) methods consistently support a homogeneous genetic structure of the three populations, and strong gene flow was detected among them by MIGRATE analyses. Moreover, there is no significant FSGS in the SM population. A relatively short time since habitat fragmentation and gene flow mediated by seed dispersal might be the likely reasons for the high levels of genetic variation and an absence of genetic structure of the coastal V. mangachapoi populations. In conclusion, even though there are no significant effects of fragmentation on the coastal V. mangachapoi forest, strict protection is required to prevent further deforestation and fragmentation. Besides, saplings of V. mangachapoi should be planted in forest gaps to reconnect fragments of the coastal forest, which would be of benefit for the long-term survival of the tropical coastal V. mangachapoi-dominated forest.
    Graphical Abstract
  • 加载中
  • Supplemental Fig. S1 Neighbor-joining tree based on Nei's (1983) genetic distance for V. mangachapoi.
    Supplemental Table S1 Analysis of molecular variance (AMOVA) for the three V. mangachapoi populations.
  • [1]

    Ghazoul J. 2016. Dipterocarp Biology, Ecology, and Conservation. 1st Edition. UK: Oxford University Press. https://doi.org/10.1093/acprof:oso/9780199639656.001.0001

    [2]

    Ashton PS. 1988. Dipterocarp biology as a window to the understanding of tropical forest structure. Annual Review of Ecology and Systematics 19:347−70

    doi: 10.1146/annurev.es.19.110188.002023

    CrossRef   Google Scholar

    [3]

    Sodhi NS, Posa MRC, Lee TM, Bickford D, Koh LP, et al. 2010. The state and conservation of Southeast Asian biodiversity. Biodiversity and Conservation 19:317−28

    doi: 10.1007/s10531-009-9607-5

    CrossRef   Google Scholar

    [4]

    Edwards DP, Tobias JA, Sheil D, Meijaard E, Laurance WF. 2014. Maintaining ecosystem function and services in logged tropical forests. Trends in Ecology and Evolution 29:511−20

    doi: 10.1016/j.tree.2014.07.003

    CrossRef   Google Scholar

    [5]

    Hu YJ, Li YX. 1992. The Tropical of Rain Forest of Hainan Island. 1st Edition. Guangdong, China: Guangdong Higher Education Press

    [6]

    Hu YJ. 1983. The phytocoenological features and types of Dipterocap forest in Hainan island. Ecological Science 2:16−24

    Google Scholar

    [7]

    Xu H, Li YD, Luo SS, Chen DX, Lin MX. 2007. Review of Vatica mangachapoi, a national key protected plant on Hainan Island. Tropical Forestry 35:8−11

    Google Scholar

    [8]

    Liang SQ, Lin Y, Yang XB, Huang SM, Fu SX, et al. 1993. The Vatica hainanensis Forest of Li Ji, Wan Ning County, Hainan Island. Natural Science Journal of Hainan University 11:1−9

    Google Scholar

    [9]

    Yang XB, Hu RG. 2000. The floral components and soil properties of forest on the tropical sandy beach. Chinese Journal of Ecology 19:6−11

    Google Scholar

    [10]

    Hu RG, Liang SQ, Lin Y. 1995. A study on the soil properties of Vatica mangachapoi forest in Li Ji, Wan Ning County, Hainan Province. Natural Science Journal of Hainan University 13:203−10

    Google Scholar

    [11]

    Wu YH, Huang QM, Liang SQ, Lin Y. 1996. Disease and pest control of Vatica mangachapoi forest in Liji Town, Wanning City, Hainan Province. Tropical Forestry 24:47−51

    Google Scholar

    [12]

    Wang G, Zhao JM, Hao QY. 2012. Landscape fragmentation of coastal Vatica mangachapoi forest nature reserve in Shimei bay. Guangdong Agricultural Sciences 39:171−4

    Google Scholar

    [13]

    Kwak MM, Velterop O, Van Andel J. 1998. Pollen and gene flow in fragmented habitats. Applied Vegetation Science 1:37−54

    doi: 10.2307/1479084

    CrossRef   Google Scholar

    [14]

    Kolb A, Diekmann M. 2005. Effects of life-history traits on responses of plant species to forest fragmentation. Conservation Biology 19:929−38

    doi: 10.1111/j.1523-1739.2005.00065.x

    CrossRef   Google Scholar

    [15]

    Lowe AJ, Boshier D, Ward M, Bacles CFE, Navarro C. 2005. Genetic resource impacts of habitat loss and degradation; reconciling empirical evidence and predicted theory for neotropical trees. Heredity 95:255−73

    doi: 10.1038/sj.hdy.6800725

    CrossRef   Google Scholar

    [16]

    Ward M, Dick CW, Gribel R, Lowe AJ. 2005. To self, or not to self.. A review of outcrossing and pollen-mediated gene flow in neotropical trees. Heredity 95:246−54

    doi: 10.1038/sj.hdy.6800712

    CrossRef   Google Scholar

    [17]

    Vekemans X, Hardy OJ. 2004. New insights from fine-scale spatial genetic structure analyses in plant populations. Molecular Ecology 13:921−35

    doi: 10.1046/j.1365-294X.2004.02076.x

    CrossRef   Google Scholar

    [18]

    Sebbenn AM, Carvalho ACM, Freitas MLM, Moraes SMB, Gaino APSC, et al. 2011. Low levels of realized seed and pollen gene flow and strong spatial genetic structure in a small, isolated and fragmented population of the tropical tree Copaifera langsdorffii Desf. Heredity 106:134−45

    doi: 10.1038/hdy.2010.33

    CrossRef   Google Scholar

    [19]

    Suzuki E, Ashton PS. 1996. Sepal and nut size ratio of fruits of Asian Dipterocarpaceae and its implications for dispersal. Journal of Tropical Ecology 12:853−70

    doi: 10.1017/S0266467400010129

    CrossRef   Google Scholar

    [20]

    Smith JR, Bagchi R, Ellens J, Kettle CJ, Burslem DFRP, et al. 2015. Predicting dispersal of auto-gyrating fruit in tropical trees: a case study from the Dipterocarpaceae. Ecology and Evolution 5(9):1794−801

    doi: 10.1002/ece3.1469

    CrossRef   Google Scholar

    [21]

    Momose K, Yumoto T, Nagamitsu T, Kato M, Nagamasu H, et al. 1998. Pollination biology in a lowland dipterocarp forest in Sarawak, Malaysia. I. Characteristics of the plant-pollinator community in a lowland dipterocarp forest. American Journal of Botany 10:1477−501

    doi: 10.2307/2446404

    CrossRef   Google Scholar

    [22]

    Lee SL, Ng KKS, Ng CH, Tnah LH, Lee CT, et al. 2016. Spatial studies of Shorea parvifolia ssp parvifolia (Dipterocarpaceae) in a lowland and hill dipterocarp forest. Journal of Tropical Forest Science 28:309−17

    Google Scholar

    [23]

    Kettle CJ, Maycock CR, Ghazoul J, Hollingsworth PM, Khoo E, et al. 2011. Ecological implications of a flower size/number trade-off in tropical forest trees. PLoS One 6:e16111

    doi: 10.1371/journal.pone.0016111

    CrossRef   Google Scholar

    [24]

    Harata T, Nanami S, Yamakura T, Matsuyama S, Chong L, et al. 2012. Fine-scale spatial genetic structure of ten Dipterocarp tree species in a Bornean rain forest. Biotropica 44:586−94

    doi: 10.1111/j.1744-7429.2011.00836.x

    CrossRef   Google Scholar

    [25]

    Finger A, Kettle CJ, Kaiser-Bunbury CN, Valentin T, Mougal J, et al. 2012. Forest fragmentation genetics in a formerly widespread island endemic tree: Vateriopsis seychellarum (Dipterocarpaceae). Molecular Ecology 21:2369−82

    doi: 10.1111/j.1365-294X.2012.05543.x

    CrossRef   Google Scholar

    [26]

    Widiyatno, Indrioko S, Na’iem M, Purnomo S, Hosaka T, et al. 2017. Effects of logging rotation in a lowland dipterocarp forest on mating system and gene flow in Shorea parvifolia. Tree Genetics and Genomes 13:85

    doi: 10.1007/s11295-017-1167-3

    CrossRef   Google Scholar

    [27]

    Doyle JJ. 1990. Isolation of Plant DNA from Fresh Tissue. Focus 12:13−15

    Google Scholar

    [28]

    Guo JJ, Shang SB, Wang CS, Zhao ZG, Zeng J. 2017. Twenty microsatellite markers for the endangered Vatica mangachapoi (Dipterocarpaceae). Applications in Plant Sciences 5:1600134

    doi: 10.3732/apps.1600134

    CrossRef   Google Scholar

    [29]

    Peakall R, Smouse PE. 2006. GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Molecular Ecology Notes, 6:288−95

    doi: 10.1111/j.1471-8286.2005.01155.x

    CrossRef   Google Scholar

    [30]

    Nei M, Chesser RK. 1983. Estimation of fixation indices and gene diversities. Annals of Human Genetics 47:253−59

    doi: 10.1111/j.1469-1809.1983.tb00993.x

    CrossRef   Google Scholar

    [31]

    Liu K, Muse SV. 2005. PowerMarker: an integrated analysis environment for genetic marker analysis. Bioinformatics 21:2128−29

    doi: 10.1093/bioinformatics/bti282

    CrossRef   Google Scholar

    [32]

    Excoffier L, Lischer HEL. 2010. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Molecular Ecology Resources 10:564−7

    doi: 10.1111/j.1755-0998.2010.02847.x

    CrossRef   Google Scholar

    [33]

    Piry S, Luikart G, Cornuet JM. 1999. BOTTLENECK: a computer program for detecting recent reductions in the effective size using allele frequency data. Journal of Heredity 90:502−3

    doi: 10.1093/jhered/90.4.502

    CrossRef   Google Scholar

    [34]

    Falush D, Stephens M, Pritchard JK. 2003. Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics 164:1567−87

    doi: 10.1093/genetics/164.4.1567

    CrossRef   Google Scholar

    [35]

    Earl DA, VonHoldt BM. 2012. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conservation Genetics Resources 4:359−61

    doi: 10.1007/s12686-011-9548-7

    CrossRef   Google Scholar

    [36]

    Jakobsson M, Rosenberg NA. 2007. CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics 23:1801−6

    doi: 10.1093/bioinformatics/btm233

    CrossRef   Google Scholar

    [37]

    Rosenberg NA. 2004. DISTRUCT: a program for the graphical display of population structure. Molecular Ecology Notes 4:137−38

    doi: 10.1046/j.1471-8286.2003.00566.x

    CrossRef   Google Scholar

    [38]

    Tamura K, Stecher G, Kumar S. 2021. MEGA11: Molecular Evolutionary Genetics Analysis Version 11. Molecular Biology and Evolution 38:3022−27

    doi: 10.1093/molbev/msab120

    CrossRef   Google Scholar

    [39]

    Beerli P. 2006. Comparison of Bayesian and maximum-likelihood inference of population genetic parameters. Bioinformatics 22:341−5

    doi: 10.1093/bioinformatics/bti803

    CrossRef   Google Scholar

    [40]

    Hardy OJ, Vekemans X. 2002. SPAGeDi: a versatile computer program to analyse spatial genetic structure at the individual or population levels. Molecular Ecology Notes 2:618−20

    doi: 10.1046/j.1471-8286.2002.00305.x

    CrossRef   Google Scholar

    [41]

    Loiselle BA, Sork VL, Nason J, Graham C. 1995. Spatial genetic structure of a tropical understory shrub, Psychotria officinalis (Rubiaceae). American Journal of Botany 82:1420−25

    doi: 10.1002/j.1537-2197.1995.tb12679.x

    CrossRef   Google Scholar

    [42]

    Doligez A, Joly HI. 1997. Genetic diversity and spatial structure within a natural stand of a tropical forest tree species, Carapa procera (Meliaceae), in French Guiana. Heredity 79:72−82

    doi: 10.1038/hdy.1997.124

    CrossRef   Google Scholar

    [43]

    Cavers S, Degen B, Caron H, Lemes MR, Margis R, et al. 2005. Optimal sampling strategy for estimation of spatial genetic structure in tree populations. Heredity 95:281−89

    doi: 10.1038/sj.hdy.6800709

    CrossRef   Google Scholar

    [44]

    Chen W, Qi H, Fu R, Li D, Jiang D. 2021. Survey and population dynamics analysis of Vatica mangachapoi germplasm resource in Hainan Island. Molecular Plant Breeding 14:4846−54

    Google Scholar

    [45]

    Hamrick JL, Godt MJW, Sherman-Broyles SL. 1992. Factors influencing levels of genetic diversity in woody plant species. New Forest 6:95−124

    doi: 10.1007/BF00120641

    CrossRef   Google Scholar

    [46]

    Ony MA, Nowicki M, Boggess SL, Klingeman WE, Zobel JM, et al. 2020. Habitat fragmentation influences genetic diversity and differentiation: Fine-scale population structure of Cercis canadensis (eastern redbud). Ecology and Evolution 10:3655−70

    doi: 10.1002/ece3.6141

    CrossRef   Google Scholar

    [47]

    Bacles CFE, Burczyk J, Lowe AJ, Ennos RA. 2005. Historical and contemporary mating patterns in remnant populations of the forest tree Fraxinus excelsior L. Evolution 59:979−90

    doi: 10.1111/j.0014-3820.2005.tb01037.x

    CrossRef   Google Scholar

    [48]

    Kikuchi S, Shibata M, Tanaka H. 2015. Effects of forest fragmentation on the mating system of a cool-temperate heterodichogamous tree Acer mono. Global Ecology and Conservation 3:789−801

    doi: 10.1016/j.gecco.2015.04.005

    CrossRef   Google Scholar

    [49]

    Tito de Morais C, Ghazoul J, Maycock CR, Bagchi R, Burslem DFRP, et al. 2015. Understanding local patterns of genetic diversity in dipterocarps using a multi-site, multi-species approach: Implications for forest management and restoration. Forest Ecology and Management 356:153−65

    doi: 10.1016/j.foreco.2015.07.023

    CrossRef   Google Scholar

    [50]

    Smith JR, Ghazoul J, Burslem DFRP, Itoh A, Khoo E, et al. 2018. Are patterns of fine-scale spatial genetic structure consistent between sites within tropical tree species? PLoS One 13(3):e0193501

    doi: 10.1371/journal.pone.0193501

    CrossRef   Google Scholar

    [51]

    Lü B, Wang N, Liu S, Hao Q. 2015. Reproductive phenological characteristics of hainan coastal Vatica mangachapoi forest. Acta Ecologica Sinica 35:416−23

    doi: 10.5846/stxb201304010572

    CrossRef   Google Scholar

    [52]

    Tackenberg O. 2003. Modeling long-distance dispersal of plant diaspores by wind. Ecological Monographs 73:173−89

    doi: 10.1890/0012-9615(2003)073[0173:MLDOPD]2.0.CO;2

    CrossRef   Google Scholar

    [53]

    Maurer KD, Bohrer G, Medvigy D, Wright SJ. 2013. The timing of abscission affects dispersal distance in a wind-dispersed tropical tree. Functional Ecology 27:208−18

    doi: 10.1111/1365-2435.12028

    CrossRef   Google Scholar

    [54]

    Kettle CJ, Hollingsworth PM, Burslem DFRP, Maycock CR, Khoo E, et al. 2011. Determinants of fine-scale spatial genetic structure in three co-occurring rain forest canopy trees in Borneo. Perspectives in Plant Ecology, Evolution and Systematics 13:47−56

    doi: 10.1016/j.ppees.2010.11.002

    CrossRef   Google Scholar

    [55]

    Wang B. 1987. Approach to the horizontal zonation of monsoon forests. Acta Phytoecologica et Geobotanica Sinica 11:154−58

    Google Scholar

    [56]

    Liang S, Lin Y, Yang X, Huang S, Fu S, et al. 1994. The Vatica mangachapoi Forest of Li Ji, Wan Ning County, Hainan Island. Natural Science Journal of Hainan University 12:14−19

    Google Scholar

  • Cite this article

    Duan J, Wang H, Tang L. 2023. Effects of habitat fragmentation on the coastal Vatica mangachapoi forest (Dipterocarpaceae) in Shimei Bay, Hainan Island, China. Tropical Plants 2:8 doi: 10.48130/TP-2023-0008
    Duan J, Wang H, Tang L. 2023. Effects of habitat fragmentation on the coastal Vatica mangachapoi forest (Dipterocarpaceae) in Shimei Bay, Hainan Island, China. Tropical Plants 2:8 doi: 10.48130/TP-2023-0008

Figures(6)  /  Tables(3)

Article Metrics

Article views(3379) PDF downloads(401)

Other Articles By Authors

ARTICLE   Open Access    

Effects of habitat fragmentation on the coastal Vatica mangachapoi forest (Dipterocarpaceae) in Shimei Bay, Hainan Island, China

Tropical Plants  2 Article number: 8  (2023)  |  Cite this article

Abstract: Habitat fragmentation can cause isolation and decline of a formerly continuously distributed population, which leads to loss of genetic variation and increased risk of extinction. Vatica mangachapoi Blanco is a dominant tree species growing in the lowland rainforests of Hainan Island, China. Remarkably, this species dominates a coastal forest in Shimei Bay, Wanning City of Hainan Province (China). Due to logging, expansion of farmland and villages, and construction of tourism facilities, the coastal V. mangachapoi-dominated forest has become fragmented, threatening its future. To evaluate the effects of habitat fragmentation on this unique coastal forest, two V. mangachapoi populations (SM and RY) along the coast and one population in the lowland rainforest near the coast were selected, and their genetic diversity was assessed based on 12 SSR markers. In addition, the genetic structure of the three populations and gene flow among them, and the fine-scale spatial genetic structure (FSGS) of the SM population were also studied. The results show that the three V. mangachapoi populations had comparable levels of genetic variation, and differentiation among them is negligible ($ {F}_{\mathrm{s}\mathrm{t}} $ = 0.008 ~ 0.013). Model-based clustering, Principal co-ordinate analysis and the Neighbor-joining (NJ) methods consistently support a homogeneous genetic structure of the three populations, and strong gene flow was detected among them by MIGRATE analyses. Moreover, there is no significant FSGS in the SM population. A relatively short time since habitat fragmentation and gene flow mediated by seed dispersal might be the likely reasons for the high levels of genetic variation and an absence of genetic structure of the coastal V. mangachapoi populations. In conclusion, even though there are no significant effects of fragmentation on the coastal V. mangachapoi forest, strict protection is required to prevent further deforestation and fragmentation. Besides, saplings of V. mangachapoi should be planted in forest gaps to reconnect fragments of the coastal forest, which would be of benefit for the long-term survival of the tropical coastal V. mangachapoi-dominated forest.

    • Trees from family Dipterocarpaceae serve an important ecosystem function in the rainforest community of Asian tropical forests, where 20%−50% of the canopy layer belong to this family[1]. With high-quality timber that has a high economic value, dipterocarp forests also form a major pillar of the global tropical timber trade[2]. Due to long-term over-harvesting and land use change, tropical rainforests have become severely fragmented, and a large number of dipterocarps are today listed as endangered species and are at risk of extinction[3]. As a result, ecosystem services of Asian tropical rainforests in which dipterocarps are the dominant species have been seriously impaired[4]. Therefore, conservation genetics studies focusing on Dipterocarpaceae are urgently needed[1].

      The tropical forests in Hainan Island, China are located at the northern edge of tropical Asia and are distinct from the typical tropical rainforests of Southeast Asia in terms of species composition, community structure and appearance as a consequence of the influence of the Asian monsoon[5]. Only three species of Dipterocarpaceae, Hopea hainanensis Merrill & Chun, H. reticulata Tardieu and Vatica mangachapoi Blanco can be found on this island. Although the species diversity of Dipterocarpaceae in Hainan Island has been greatly reduced as compared to that in Southeast Asia, the three species, especially V. mangachapoi (Fig. 1), play a key role in community assembly and ecosystem functioning of the lowland rainforests of this island[57]. It is remarkable that V. mangachapoi has developed into a continuously distributed coastal forest growing on sand substrate with 25 kilometers-long and 400 to 500 meters-wide at Shimei Bay, Wanning City (China), which is estimated to be at least 4000 years old[8, 9]. A study showed that soil moisture and organic matters of the sand substrate are much lower than those of normal tropical soils[10]. The formation of the coastal V. mangachapoi-dominated forest on barren and harsh sandy beach is unique and rare in itself, which could serve as an example to study the underlying physiological and genetical adaptation of V. mangachapoi to arid and poor substrate. In recent years, due to such factors as coastal development and village expansion, the area of the coastal V. mangachapoi-dominated forest in Shimei Bay has been reduced, and the formerly intact population has fragmented into several isolated patches. Coupled with the presence of forest gaps and fungal disease caused by human interference, the survival of the coastal V. mangachapoi is seriously threatened[9, 11, 12]. Conservation management is thus needed to protect this unique coastal forest dominated by V. mangachapoi.

      Figure 1. 

      Morphology of Vatica mangachapoi. (a) An individual tree, (b) flowers, (c) fruits.

      Habitat fragmentation can cause an intact population into small, isolated patches, with reduced gene flow between patches and increased genetic drift and inbreeding within patches[1316]. If seed dispersal is limited and restricted within patches, genotypes are likely to be spatially clustered, producing a strong fine-scale spatial genetic structure (FSGS)[17, 18]. Vatica magachapoi has winged fruits, which may promote seed dispersal by wind. Studies indeed showed that the dispersal distance of winged fruits are generally further than non-winged fruits in dipterocarps[19, 20]. On the other hand, the strength of FSGS was also affected by pollen dispersal. Dipterocarps that are pollinated by large insects, such as bees, can achieve longer distance of pollen flow than those pollinated by small insects, such as thrips, because large pollinators can move further than small ones[21, 22]. The limited seed and pollen dispersal were confirmed as the main reason for significant FSGS of most dipterocarps in fragmented habitats[2326]. Is there a significant FSGS in the coastal forest of V. mangachapoi? Is genetic diversity lower in the coastal V. mangachapoi populations than the undisturbed rainforest populations nearby? Does significant genetic differentiation occur between patches of the coastal V. mangachapoi forest? These questions are important for the conservation and management of the unique coastal dipterocarp forest but remain to be resolved.

      To answer the above questions, two coastal populations of V. mangachapoi (SM and RY) were sampled, and one population in the lowland rainforest near the coast (TT) were further collected for comparison. Genetic diversity, structure and population differentiation were assessed for the three V. mangachapoi populations using 12 SSR markers. Gene flow was estimated to test whether habitat fragmentation interrupted genetic exchange between them or not. Finally, FSGS were analyzed using the SM population to show whether significant spatial genetic structure has occurred within patches. Answering these questions could shed light on the conservation and continued survival of this unique coastal V. mangachapoi-dominated forest.

    • The study area is located in the Provincial Nature Reserve of V. mangachapoi, Wanning City, Hainan Province (China). Two V. mangachapoi populations (SM and RY) in the coastal forest separated by villages, roads and human facilities, and one population (TT) in the lowland rainforest near the coast were selected (Table 1, Fig. 2). In total, 188 V. mangachapoi trees individually spaced out more than 25 m apart and with DBH > 5 cm were sampled. Mature leaves lacking disease spots were selected, dried by silica gel and then were stored in a −20 °C refrigerator. The voucher specimens of V. mangachapoi were kept in Hainan University (Hainan, China).

      Table 1.  Genetic diversity indices of the three V. mangachapoi populations based on 12 SSR markers.

      PopulationLocationNNaNeHoHeFis
      SM110.26691° E, 18.66671° N9183.6470.5470.6900.207
      RY110.17952° E, 18.59768° N3973.7040.6050.7000.142
      TT110.24941° E, 18.67744° N587.53.6940.5660.6920.167
      Average7.53.6820.5720.6940.172

      Figure 2. 

      Geographic distribution of the coastal V. mangachapoi-dominated forest and the locations of the three sampled V. mangachapoi populations. Gene flow among them was estimated, with the width of lines being proportional to the intensity of gene flow.

    • A modified CTAB method[27] was used to extract the genomic DNA of V. mangachapoi. Twelve pairs of polymorphic SSR primers developed by Guo et al.[28] were used in this study. PCR amplification were performed in a total volume of 10 μL, containing 1.0 μL of genomic DNA (around 50 ng), 5.0 μL of Taq PCR Master Mix (GeneTech), 0.5 μL of forward and reverse primers, and 3.0 μL of ddH2O. Amplification was carried out as follows: pre-denaturation at 95 °C for 5 min, followed by 35 cycles of denaturation at 95 °C for 20 s, annealing at 52−62 °C for 15 s, extension at 72 °C for 30 s, and finally extension at 72 °C for 7 min. PCR products were separated by capillary electrophoresis using ABI3730xl (Applied Biosystem) and SSR genotypes were analyzed by the GeneMarker software.

    • Population genetic parameters, including number of alleles (Na), effective number of alleles (Ne), observed heterozygosity (Ho), expected heterozygosity (He), inbreeding coefficient (Fis) and genetic differentiation among populations ($ {F}_{\mathrm{s}\mathrm{t}} $) were estimated using GenAlex 6.51[29]. Polymorphism information content (PIC) and Nei & Chessers[30] genetic distances were calculated by PowerMarker 3.25[31]. Analysis of molecular variance (AMOVA) of the V. mangachapoi population was performed using Arlequin 3.5[32]. Potential population bottleneck was examined by Wilcoxon sign-rank test and model shift using the Stepwise Mutation Model and Two-phased Mutation Model implemented in Bottleneck 1.3.2[33].

      Bayesian clustering analysis was performed using Structure 2.3.4[34]. The values of K were set from 1 to 10, and for each value of K, 10 independent replicates were run with 100,000 burn-in iterations followed by 200,000 MCMC (Markov chain Monte Carlo) iterations. The best K was determined according to the delta K of STRUCTURE Harvester[35]. The results of the 10 replicates were combined by the Greedy algorithm implemented in Clumpp1.1.2[36], and the result of individual clustering were drawn using Distruct1.1[37]. NJ trees were constructed using MEGA 11.0[38] based on Nei & Chessers[30] genetic distances. Principal co-ordinate analysis (PCoA) was performed with GenAlex 6.51.

      The effective population size (θ) of the three V. mangachapoi populations and the migration rate (M) between them were calculated using MIGRATE[39], a software based on the coalescent theory and Bayesian inference to estimate values of parameters of a user-specified population model. MIGRATE analyses were run under a Brownian motion model, and four heat chains with different temperatures of 1.0, 1.5, 3.0 and 1.0 × 105 were simulated. Gene flow was calculated according to the equation Nm = θ*M/x. For SSR markers, x was set to 4. Three independent replications were run to ensure the convergence of the Markov chain Monte Carlo methods implemented in MIGRATE.

      The fine-scale spatial genetic structure (FSGS) within-population was assessed using SPAGeDi 1.5[40]. The kinship coefficients (Fij, kinship coefficients) between any two individuals were calculated and regressed against the natural logarithm of the spatial distance to obtain the regression slope bF[41]. We divided the distance between any pair of individuals sampled from the SM population into 10 distance classes (35, 50, 75, 100, 150, 300, 500, 700, 850, 1,100 m), with at least 30 pairs of individuals per distance class[42, 43].

      The 95% confidence intervals of Fij were calculated from 9,999 permutations of spatial distance among pairs of adults for 10 distance classes. If the Fij was higher than the upper bound of the 95% confidence interval, there is significant spatial genetic structure in population and a high level of genetic similarity among individuals; if the Fij fell within the 95% confidence interval, there is no spatial genetic structure in population and individuals were considered to be spatially randomly distributed; if the Fij was less than the lower bound of the 95% confidence interval, individuals were considered to be uniformly distributed in space without spatial genetic structure. The value of the Sp statistic reflects the strength of FSGS and is defined as Sp = -bF / (1-F(1) ), where bF is the regression slope, and F(1) is the mean pairwise kinship coefficient of the first distance class.

    • There is no significant difference in the level of genetic diversity between the coastal (SM and RY) and the rainforest (TT) populations (Table 1). The inbreeding coefficient was greater than 0, indicating inbreeding and an excess of homozygotes in the three V. mangachapoi populations. Totally 90 alleles were detected from the 12 SSR loci, and the number of alleles at a single locus ranged from 4.667 to 11.000, with an average of 7.500 alleles per locus. The observed and expected heterozygosity ranged from 0.303 to 0.769 and from 0.415 to 0.808, respectively. The primer sequences, range of allele sizes and genetic diversity indices of the 12 SSR loci are shown in Table 2.

      Table 2.  Primer sequences, allele size and genetic diversity indices of the 12 SSR markers.

      LociPrimer sequences (5’-3’)Repeat
      motif
      Allele sizeGenAlexPowerMarker
      NaNeHoHePIC
      VM1F:GAACCCTTATTGGCCTGCCTAC(AT)11166−1847.3334.2310.7400.7630.7430
      R:GGGACCAAATGACTTGAGTAATCT
      VM2F:ACCCTAACAATTCTCTTTGTTTCCT(TAA)11152−1959.6674.1200.5130.7550.7364
      R:CCCCAATCTCAGTAAGGACTCA
      VM3F:CTTGTGTCGAGCATGCATGTAT(AT)11175−1918.3334.8570.7610.7930.7659
      R:TGCTGGCCTTTTATGTTAGGGT
      VM4F:ATAGCAGGCACTTCGGAAGTAC(TA)8261−2778.6674.6130.3700.7810.7533
      R:CCTGAGAAACAAAGCAACGCAT
      VM5F:GCACTAGCACTAGCACTAGCTT(CT)11218−2264.6672.9080.6290.6510.6026
      R:GGCTTTTCCAATTTCCATGGCT
      VM6F:AGTTAAGGGACCAAATTTAGCGT(TA)7259−2695.0002.7940.5930.6360.5902
      R:GTGTTTGTCAACTGGGCTTCAA
      VM7F:CCCATGTGCTAGGCTAATGCTA(AT)6229−2395.0002.3940.3030.5820.5409
      R:AAATCAGCATGAAACTTCTCCATT
      VM8F:CACCACCACAGGCTTGAGTATA(TA)7168−1825.6671.7220.3740.4150.4044
      R:GAAGGCCAACTAATCAAGCTGC
      VM9F:TCATTTCTGTCTCACTCGACCC(TTC)10148−1685.6673.0100.6390.6660.6097
      R:TCATCGACGAATCACTGTTCGA
      VM10F:ACGGATAAGTTAACGGACTAGACA(TA)10215−2279.3334.7130.5680.7760.7997
      R:AGATTTTCCCCCAGTCATCGAC
      VM11F:GCTGGCACTTAGGATGCCTTAA(ATT)11138−15011.0003.5640.6100.7020.6657
      R:AGCAACCAATTAGCTCAAATCAA
      VM12F:GGGCAGCCTCGTAAATCAATTAC(ATT)13225−2499.6675.2530.7690.8080.7958
      R:ATTACCTGGCACAACCTTAGCC

      Genetic differentiation was weak among the three V. mangachapoi populations (Fst = 0.008~0.013). The result of AMOVA showed that 99% of genetic variation was partitioned within population, in line with little divergence among populations (Supplemental Table S1).

      The Wilcoxon sign rank test found that the p-values were not significant under either the S.M.M or the T.P.M model for the three V. mangachapoi populations, and their allele frequency distributions were generally L-shaped (Fig. 3), indicating that the three populations have not experienced genetic bottlenecks recently.

      Figure 3. 

      Allele frequency distribution of the three V. mangachapoi populations.

    • STRUCTURE analysis found that delta K was maximized at K = 2, indicating two genetic clusters of the studied V. mangachapoi populations. The distribution of the two clusters did not differ significantly between the three populations, and this is also true for K = 3 or 5 (Fig. 4). Consistent with the results of STRUCTURE analyses, NJ tree (Supplemental Fig. S1) and PCoA analysis (Fig. 5) also suggest a homogeneous genetic structure of the three V. mangachapoi populations.

      Figure 4. 

      Results of STRUCTURE analysis. (a) Best K determined using the delta K method. (b) Log probabilities and delta K values for K from two to ten. (c) The results of individual assignment at K = 2, 3 and 5. Each vertical bar represents an individual, and the proportion of the colors corresponds to the posterior probability of genetic clusters assigned to each individual.

      Figure 5. 

      Principal co-ordinate analysis (PCoA) based on Nei & Chessers[30] genetic distance among individual samples of V. mangachapoi.

    • The effective population sizes of the three V. mangachapoi populations estimated by MIGRATE were similar, but the intensity of gene flow varied among them (Fig. 2, Table 3). The gene flows from RY to the other two populations were less than their reverse gene flows, however, the gene flows from SM to the other two populations were greater than their reverse gene flows. These results suggested that gene flows between the V. mangachapoi populations were asymmetric.

      Table 3.  Mutation-scaled migration rate, effective population size and gene flow estimated by program MIGRATE.

      Direction of
      gene flow
      Migration
      rate (M)
      Effective population
      size (θ)
      Gene flow (Nm)
      SM→RY129.615θSM = 0.097903.172327
      SM→TT179.8394.401560
      RY→SM63.241θRY = 0.096861.531380
      RY→TT117.4902.845020
      TT→SM115.412θTT = 0.097462.812013
      TT→RY134.0623.266421
    • No spatial genetic structure was detected at any of the 10 distance classes in the SM population (Fig. 6). The values of Fij were less than zero over multiple distance classes, indicating that individual trees of V. mangachapoi were spatially uniformly distributed. Based on the mean affinity (F(1)) for the first distance class (0.0151) and the regression slope bF (−0.004605), the strength of FSGS (Sp) was derived as 0.004675 for the V. mangachapoi population in Shimei Bay.

      Figure 6. 

      Fine-scale genetic structure of V. mangachapoi in Shimei Bay. The solid line represents the mean Kinship coefficient F (Loiselle et al.[41]), and the dashed lines represent the 95% confidence intervals of the mean Kinship coefficient F.

    • Due to rapid coastal development and village expansion, the coastal population of V. mangachapoi declined and fragmented into several small and discontinuous patches[12] (Fig. 2). However, there was no significant difference in genetic diversity between the coastal (SM and RY) and the rainforest (TT) populations (Table 1), and the effective sizes of the three populations are quite close (Table 3). Besides, no recent bottleneck could be detected in these populations (Fig. 3). The results indicated that the sizes of the two coastal populations are still large to maintain a comparative level of genetic diversity relative to that of the rainforest population nearby[44]. In addition, the three V. mangachapoi populations were demonstrated to share a homogenous genetic structure and there is little differentiation between them (Figs 4 & 5). In summary, patterns of SSR variation observed in V. mangachapoi suggested either genetic connection through gene flow or not enough time to accumulate divergence after fragmentation of the coastal forest dominated by V. mangachapoi.

      Frequent gene flow can prevent rapid loss of genetic variation and differentiation between patches in a fragmented population[45, 46]. Vateriopsis seychellarum is endemic in the Seychelles, after a long period of logging, only a few hundred adult trees remained. Nevertheless, a relatively high level of genetic diversity was found in this species. Long-distance gene flow between isolated patches of Va. seychellarum was considered as the main reason to maintain genetic variation in this species[25]. If pollinators can travel across the gaps created by fragmentation, pollen-mediated gene flow could be maintained between patches, as a result, genetic drift happened within patches would be mitigated in the short term[26, 47, 48]. In this study, frequent gene flow (Nm > 1) was detected between the three V. mangachapoi populations. Besides, the time of fragmentation of the coastal V. mangachapoi forest is relatively short comparing with the generation time of this species. Differentiation between patches is probably impeded by gene flow or there is not enough time to accumulate significant divergence among them, and genetic variation could be largely maintained within the coastal V. mangachapoi forest.

    • No significant FSGS was detected in the SM population of V. mangachapoi, which may result from long distance dispersal of seeds and/or pollens within population[24, 49]. Two of the five sepals of V. mangachapoi flowers keep growing and develop into functional wings of the fruits that would promote seed dispersal by wind[19, 50]. Hainan Island will experience frequent Pacific typhoon from June to October each year, and the time of fruit ripening of V. mangachapoi (mid to late July-August) coincides well with the activities of the Pacific typhoon[51]. Moreover, Wanning City is one of the locations where typhoons frequently make landfall on Hainan Island. The coastal V. mangachapoi forest at Shimei Bay would be highly likely to encounter a Pacific typhoon during its fruit ripening period. Strong convection currents can carry winged fruits of V. mangachapoi into the upper air and achieve a long-distance horizontal dispersal with hundreds of meters, which may account for the lack of significant FSGS in the SM population[52, 53].

      Pollen-mediated gene flow can also influence the strength of FSGS, and restricted gene flow generally results in significant FSGS[22, 24, 54]. Kettle et al. studied the FSGS of three dipterocarp species from the tropical rainforests at Borneo[54]. No clear signal of FSGS was detected in Dipterocarpus grandiflorus, a species pollinated by large pollinators with strong mobility, which may facilitate long-distance pollen flow. On the contrary, S. xanthophylla and Parashorea tomentella had significant FSGS, probably because they were pollinated by small pollinators and consequently short distances of pollen flow. Lee et al. studied the FSGS of S. parvifolia populations in different habitats[22]. They found that populations in montane rainforests, which were mainly pollinated by large pollinators, had no FSGS, whereas populations in lowland rainforests, which were pollinated by small pollinators, had significant FSGS. Lee et al. suggested that it is the restricted pollen flow that leads to strong FSGS in the lowland dipterocarp rainforests[22]. However, as the distance of pollen-mediated gene flow of V. mangachapoi is unclear at present, the relative contribution of seed and pollen dispersal to the random distribution of genotypes in space deserve further studies.

    • In this study, we demonstrated that there was no significant difference in genetic diversity among the two coastal V. mangachapoi populations and one rainforest population near the coast. Moreover, little differentiation and frequent gene flow were detected between the three populations. Even though the coastal populations maintain a relatively high level of genetic diversity, the extremely simple community structure and poor species richness of the coastal V. mangachapoi-dominated forests indicate that this unique community is much more fragile than other dipterocarp communities[55, 56]. In addition, the coastal V. mangachapoi-dominated forest is likely to degrade into a sandy scrub community or bare sandy beach in the near future if intense anthropogenic disturbance persists. Therefore, further logging and invasion of the coastal forest must be strictly prohibited, and saplings of V. mangachapoi should be planted in forest gaps to promote the restoration of the coastal V. mangachapoi-dominated forest landscape.

    • The fragmented coastal V. mangachapoi-dominated forest in Shimei Bay have not yet exhibited significant genetic differentiation and diversity loss. The winged fruits of V. mangachapoi may promote seed dispersal and maintain gene flow between populations, which could mitigate genetic drift and lead to random distribution of genotypes within population. In addition, comparing with the generation time of V. mangachapoi, the time of fragmentation of the coastal V. mangachapoi forest is relatively short, so there is not enough time to accumulate differentiation for populations from the fragmented forest. Based on the above findings, we suggest to strengthen the protection of the coastal V. mangachapoi forest to prevent further deforestation. Besides, saplings of V. mangachapoi should be planted to connect isolated populations and facilitate the restoration of the unique coastal V. mangachapoi forest.

      • We thank W-Q Xiang, S-Q Zhu and Y Cai of Hainan University for their help on data analysis and collection of samples. We also thank the Provincial Nature Reserve of Vatica mangachapoi at Wanning City for their help on field research. This study was supported by the National Natural Science Foundation of China (Grant No. 32060236) and the Key R&D Projects of Hainan Province (Grant No. ZDYF2022XDNY260).

      • The authors declare that they have no conflict of interest.

      • Received 6 May 2023; Accepted 14 June 2023; Published online 5 July 2023

      • There is no significant difference in genetic diversity between the coastal and lowland V. mangachapoi populations.

        No fine-scale spatial genetic structure was detected in the coastal SM population.

        Strong gene flow between populations accounts for the homogenous genetic structure of the coastal and lowland populations.

        There is not enough time to accumulate differentiation between fragmented populations in the coastal forest.

      • Supplemental Fig. S1 Neighbor-joining tree based on Nei's (1983) genetic distance for V. mangachapoi.
      • Supplemental Table S1 Analysis of molecular variance (AMOVA) for the three V. mangachapoi populations.
      • Copyright: © 2023 by the author(s). Published by Maximum Academic Press on behalf of Hainan University. This article is an open access article distributed under Creative Commons Attribution License (CC BY 4.0), visit https://creativecommons.org/licenses/by/4.0/.
    Figure (6)  Table (3) References (56)
  • About this article
    Cite this article
    Duan J, Wang H, Tang L. 2023. Effects of habitat fragmentation on the coastal Vatica mangachapoi forest (Dipterocarpaceae) in Shimei Bay, Hainan Island, China. Tropical Plants 2:8 doi: 10.48130/TP-2023-0008
    Duan J, Wang H, Tang L. 2023. Effects of habitat fragmentation on the coastal Vatica mangachapoi forest (Dipterocarpaceae) in Shimei Bay, Hainan Island, China. Tropical Plants 2:8 doi: 10.48130/TP-2023-0008

Catalog

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return