In 2003, Kandelia obovata was identified as a new mangrove species differentiated from Kandelia candel. However, little is known about their chloroplast (cp) genome differences and their possible ecological significance. In this study, 25 whole cp genomes, with seven samples of K. candel from Malaysia, Thailand, and Bangladesh and 18 samples of K. obovata from China, were sequenced for comparison. The cp genomes of both species encoded 128 genes, namely 83 protein-coding genes, 37 tRNA genes, and eight rRNA genes, but the cp genome size of K. obovata was ~2 kb larger than that of K. candle due to the presence of more and longer repeat sequences. Of these, tandem repeats and simple sequence repeats exhibited great differences. Principal component analysis based on indels, and phylogenetic tree analyses constructed with homologous protein genes from the single-copy genes, as well as 38 homologous pair genes among 13 mangrove species, gave strong support to the separation of the two species within the Kandelia genus. Homologous genes ndhD and atpA showed intraspecific consistency and interspecific differences. Molecular dynamics simulations of their corresponding proteins, NAD(P)H dehydrogenase chain 4 (NDH-D) and ATP synthase subunit alpha (ATP-A), predicted them to be significantly different in the functions of photosynthetic electron transport and ATP generation in the two species. These results suggest that the energy requirement was a pivotal factor in their adaptation to differential environments geographically separated by the South China Sea. Our results also provide clues for future research on their physiological and molecular adaptation mechanisms to light and temperature.

Molluscs are an important component of the mangrove ecosystem, and the vertical distributions of molluscan species in this ecosystem are primarily dictated by tidal inundation. Thus, sea-level rise (SLR) may have profound effects on mangrove mollusc communities. Here, we used two dynamic empirical models based on measurements of surface elevation change, sediment accretion and zonation patterns of molluscs to predict changes in molluscan spatial distributions in response to different sea-level rise rates in the mangrove forests of Zhenzhu Bay (Guangxi, China). The change in surface elevation was 4.76–9.61 mm a−1 during the study period (2016–2020), and the magnitude of surface-elevation change decreased exponentially as original surface elevation increased. Based on our model results, we predicted that mangrove molluscs might successfully adapt to a low rate of SLR (marker-horizon model: 2–4.57 mm a−1; plate model: 2–5.20 mm a−1) by 2100, with molluscs moving seaward and those in the lower intertidal zones expanding into newly available zones. However, as SLR rate increased (marker-horizon model: 4.57–8.14 mm a−1; plate model: 5.20–6.88 mm a−1), our models predicted that surface elevations would decrease beginning in the high intertidal zones and gradually spreading to the low intertidal zones. Finally, at high rates of SLR (marker-horizon model: 8.14–16.00 mm a−1; plate model: 6.88–16.00 mm a−1), surface elevations were predicted to decrease across the elevation gradient, with molluscs moving landward and species in higher intertidal zones would be blocked by landward barriers. Tidal inundation and the consequent increase in interspecific competition and predation pressure were predicted to threaten the survival of many molluscan groups in higher intertidal zones, especially species at the landward edge of the mangroves. Thus, future efforts to conserve mangrove floral and faunal diversity should prioritize species restricted to landward mangrove areas.