Wen-Xun Lu
Bing-Huang Zhang1†- 1College of the Environment and Ecology, Xiamen University, Xiamen, China
- 2School of Life Sciences, Peking University, Beijing, China
Temperature is one of the climatic factors that shape the geographic distribution of plant populations. Mangroves are temperature-sensitive plants, and their distributions are severely limited by low temperatures. It is unknown, however, to what extent temperature contributes to their population differentiation and evolution. Kandelia obovata (Rhizophoraceae) is a mangrove species with high cold tolerance in the Northern Hemisphere. We investigated the phenotypic responses of an artificial population of K. obovata, with plants transplanted from different source populations, to extremely low temperatures during winter of 2015–2016 in Yueqing County (28°20′N), Zhejiang Province of China. Using two binary traits, “with/without leaves alive on the branches” and “with/without alive buds on the tips of branches,” we classified plants in this artificial population into strong, moderate and poor cold resistance groups. We further assessed the genetic diversity, structure and differentiation of these three groups, as well as five natural populations along a latitudinal gradient using ten nuclear and six plastid microsatellite markers. Microsatellite data revealed genetic differentiation among the natural populations along the latitudinal gradient. Molecular data indicated that the cold tolerance of three groups in the artificial population was associated with their geographic origins, and that the most cold-tolerant group came from the northernmost natural population. Our study thus indicates that natural populations of K. obovata may have evolved divergent capacity of cold tolerance.
Introduction
Climate is one of the primary factors that shape the natural distribution, genetic and geographic structures of plant populations (MacArthur, 1984; Woodward, 1987; Etterson, 2004; Lomolino et al., 2010). Available evidence has shown that habitat preference, with a genetic basis, can be maintained by natural selection, which may eventually generate a series of genetic lineages corresponding to different environments (Lenormand, 2002; Kawecki and Ebert, 2004; Sork, 2016). Temperature, for example, can induce altitudinal and latitudinal differentiation of plant populations (Saxe et al., 2001; Jump et al., 2006; Halbritter et al., 2018). In addition, there are several case studies, such as on Pseudotsuga menziesii (Pinaceae) and Cicer arietinum (Fabaceae), that investigated the relationships of population genetic structure with phenotypic traits and with cold tolerance (Eckert et al., 2009; Sani et al., 2018).
Mangroves are temperature-sensitive plants. Both the annual coldest air temperature and the winter sea surface temperature can directly or indirectly affect mangroves (Chapman, 1977; Duke et al., 1998; Lovelock et al., 2016). Therefore, their distributions appear to be strongly limited by temperature in spite of some outliers, such as Avicennia marina var. australasica (Acanthaceae) in New Zealand (38°45′S), Avicennia germinans (Acanthaceae) in Bermuda (32°20′N), and Kandelia obovata in Japan (31°N), which may move from mangrove-rich regions to their latitudinal limits by the aids of warm ocean currents (Tomlinson, 2016). According to their adaptability to temperature, mangroves can be grouped into three thermal categories: thermophilic stenotopic species, thermophilic eurytopic species, and winter-resistant eurytopic species (Li and Lee, 1997; Hu et al., 2012). Since the 1980s, a series of studies have investigated the physiology, ecology, and genetics of cold resistance of mangroves (Lin et al., 1994; Yang and Lin, 1997, 1998; Yang et al., 1999, 2003; Wang et al., 2011; Chen et al., 2012, 2017; Zheng et al., 2012; Peng et al., 2015; Su et al., 2019). These studies revealed that morphological and physiological traits, such as tree defoliation, leaf photosynthesis rate, soluble sugar content, and leaf conductivity are indicators of the cold resistance of mangroves. Yang and Lin (1997) measured the cold-resistant ability of leaves of two mangrove species, K. obovata and Aegiceras corniculatum (Myrsinaceae), and found that their cold-resistant ability increased with decreasing air temperature. In a previous study, where the cold-hardiness of leaves of 16 mangrove species was measured by conductivity method, Yang and Lin (1998) found a latitudinal pattern of population differentiation, that is, a positive relationship between latitude and cold tolerance. Wu et al. (2018) examined the relationship between temperature and the distribution of mangrove species, and found that the species composition of a mangrove assemblage was significantly correlated with winter temperature. Yet, it is unknown to what extent temperature contributes to the population structure and evolution of mangrove species. We currently lack evidence that populations of a single mangrove species may differentiate under temperature control.
Kandelia obovata is a widely distributed mangrove species, and its distribution range contains the Ganges Delta, Burma, Southeast Asia, the Ryukyu Islands, and even southern Japan (Tomlinson, 2016). As a monoecious mangrove species, K. obovata is self-compatible, and depends on small insects for pollination (Tomlinson, 2016). It is regarded as one of the mangrove species with the highest cold tolerance in the Northern Hemisphere. Its natural northern range limit in China is Fuding City of Fujian Province (27°26′N) with an average winter temperature of 9.9°C (Lin, 1999; Wang et al., 2011). In 1957, K. obovata was successfully transplanted from Hainan Province to the Yueqing County (28°20′N) of Wenzhou City in Zhejiang Province for the first time (Mu et al., 2005). Subsequently, propagules of K. obovata from other areas, such as Fuding County in Fujian Province and Cangnan County in Zhejiang Province, were gradually introduced to Yueqing (Lin, 1999; Mu et al., 2005; Chen et al., 2019). Nowadays, K. obovata plants introduced from different areas and at different periods have survived in Yueqing, which forms a mixed artificial population (Zheng et al., 2016). In a study on the population genetics of K. obovata along the southeastern coast of China, Ruan et al. (2013) found significant differentiation between most population pairs using AFLP markers. Other than that, Yang and Lin (1998) demonstrated that the cold-hardiness of leaves in different K. obovata populations increased with increasing latitude. However, it is unknown whether the difference in cold resistance of K. obovata populations is associated with their genetic differentiation. From 2008 to 2016, several extreme cold climatic events occurred in South China, having a pervasive influence on mangroves (Chen et al., 2017; Su et al., 2019). Especially, a chilling event in the winter of 2015–2016 caused differing degrees of defoliation and bud withering in accessions of K. obovata that were introduced from different natural populations to Yueqing. This phenomenon suggests different adaptability of source populations of K. obovata to extremely low temperatures.
To understand whether the observed divergence of cold resistance has genetic underpinnings, we investigated the aforementioned artificially mixed population of K. obovata in a “common garden” condition in comparison with five natural populations along a latitudinal gradient. We analyzed phenotypic responses of K. obovata to coldness at individual levels in Yueqing. Besides, we analyzed the population genetic diversity, structure and differentiation of Yueqing as well as five natural populations using ten nuclear and six plastid microsatellite markers. Furthermore, we compared the cold tolerance with the genetic data in Yueqing, and traced the geographic origin of individuals in this artificial population. Finally, we test the hypothesis that cold tolerance of K. obovata may exhibit a latitudinal gradient.
Materials and Methods
Study Sites and Timing
According to the available literature, Fuding (FD) was one of the known exact source populations for the Yueqing (YQ) artificial population (Chen et al., 2019). Populations from Hainan Island and the south of Fujian Province were also introduced to YQ (Li et al., 2001; Mu et al., 2005). Therefore, our experiments focused on the YQ artificially “mixed” population; in parallel, five populations along the southeast coast of China, including one exact source population (FD), two relatively accurate sources (Haikou, HK and Yunxiao, YX) and two possible provenances (Shenzhen, SZ and Fangchenggang, FCG) of YQ, were studied for comparison (Figure 1A and Table 1).
Figure 1. (A) Map of the study areas in China with details of the localities of six K. obovata populations and their average temperatures in January (from China Meteorological Data Sharing Service System). (B) K. obovata grown in YQ, from top to bottom, before the chilling event, 5 months after the chilling event, and 10 months after the event, respectively. The red boxes in the pictures indicate the same individual.
Table 1. Geographic information and meteorological characteristics of the areas where the studied populations are located.
The chilling event that prompted this study occurred in YQ during the winter of 2015–2016 when the temperature fell to −5°C1. On June 17, 2016, we conducted the first survey and collected materials in YQ. Over a growing season approximately 10 months after the chilling event, our second “recovery” survey and materials sampling started on November 12, 2016 (Figure 1B). The investigation and sampling of five natural populations were conducted between August 20 and November 21, 2016.
Sample Collection and Analysis of Morphological Traits
Almost all surviving adult individuals in YQ were sampled (126 samples for the first and 34 samples for the second investigation). From the five natural populations, 117 individuals were sampled with a distance of more than 10 m between each pair of individuals. In total, 277 individuals were sampled from six sites. Leaves were collected from each selected individual and then desiccated with silica gel, and then preserved at room temperature until DNA extraction. For the measurement of morphological traits, three mature leaves of each individual were used to measure leaf length (LL) and leaf width (LW), except for some individuals of the YQ population with a defoliation rate over 90%. In total, we obtained measurements of six morphological traits, e.g., plant height (PH), stem basal diameter (SBD), diameter at breast height (DBH), crown diameter (CD), leaf length (LL), and leaf width (LW; Supplementary Table 1; five traits for YQ and FD as their plant height was lower than 1.3 m rendering DBH missing data). Since the YQ population had experienced the chilling event, its individuals were recorded for several traits that may be directly related to cold stress, such as growth form (GF), defoliation percentage (DP), chilling damage to branches (CDB), and germination position (GP; Supplementary Table 1). All of the morphological traits measurement data used for subsequent analyses is listed in Supplementary Data Sheet 1. The average value, standard error (SE), and coefficient of variation (CV) of the above six traits and these cold-stress related traits were analyzed using Excel 2019 (Microsoft Inc., Redmond, WA, United States) and SPSS 22.0 software (SPSS Inc., Chicago, IL, United States). In addition, five quantitative variables (PH, SBD, CD, CDB, and DP) were standardized by subtracting the mean value, which was divided by the standard deviation to remove scale effects. A principal component analysis (PCA) was performed on those five quantitative traits using R version 4.0.02 to depict the differences among different individuals in YQ.
DNA Extraction and Screening of Microsatellite Fragments
Total genomic DNA was extracted from dried leaf samples using the Plant Genomic DNA Kit (Tiangen) following the standard protocol provided by the manufacturer. DNA extractions were checked for concentration and purity using NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific), and subsequently diluted to a concentration of 5–20 ng/μL and stored at −20°C. For PCR amplification, 16 microsatellite markers (ten nSSR markers and six cpSSR markers) developed by previous studies were used (Supplementary Table 2; Sugaya et al., 2002; Islam et al., 2006, 2008), and PCR reactions were performed following the published protocol (Sugaya et al., 2002; Islam et al., 2006, 2008). Forward primers for PCR amplifications were fluorescence-labeled with HEX (green) or ROX (red). Amplification products were screened using capillary electrophoresis (CE) [Sangon Biotech (Shanghai) Co., Ltd.]. For samples that exhibited ambiguous bands, amplifications were repeated twice with sterilized water used as the negative control to exclude systemic error.
Nuclear Simple Sequence Repeats Genetic Diversity and Differentiation
Although 160 samples were collected in YQ, there were 12 samples with missing data, and therefore a total of 148 individuals in YQ and 114 individuals in five natural populations (three samples with missing data were excluded) were included in the subsequent analysis. The simple sequence repeats (SSRs) data (both nSSR loci and plastid SSR loci) used for subsequent analyses is listed in Supplementary Data Sheet 2. For both YQ artificial population and five natural populations, GENEPOP 3.4 (Raymond and Rousset, 1995) was used to test for null alleles, the Hardy--Weinberg equilibrium (HWE) and for linkage disequilibrium (LD). The nSSR loci under selection were detected with BayeScan 2.03. Both neutral loci and all loci were used to calculate genetic diversity, but only neutral loci were used to infer population structure. Genetic diversity parameters including the number of alleles (Na), number of effective alleles (Ne), Shannon diversity index (I), observed heterozygosity (Ho), expected heterozygosity (He), and unbiased expected heterozygosity (uHe) were calculated using GenAlEx 6.5 (Peakall and Smouse, 2012). Population pairwise FST was computed using Genetix 4.05 with 1,000 permutations to calculate the significance of FST (Belkhir et al., 1996-2004). The analysis of molecular variance (AMOVA) was performed to assess the distribution of genetic variance among different populations using Φ-statistics in GenAlEx 6.5 (Peakall and Smouse, 2012). Two AMOVA’s were conducted, one is to group individuals by sampling locality, and the other to group individuals by different cold resistance in YQ population. The potential structure of the populations (YQ and five natural populations) was inferred from a Bayesian approach using STRUCTURE 2.3.3 based on the admixture model with correlated allele frequencies between populations (Pritchard et al., 2000). For each simulated value of K (the number of clusters, varying from 1 to 10), 20 independent runs were performed with the Markov chain Monte Carlo (MCMC) run for 1,000,000 replications and a burn-in period of 100,000 interactions. The most likely number of K was determined following the method described by Evanno et al. (2005).
Plastid Simple Sequence Repeats Genetic Diversity and Geographic Patterns
GenAlEx 6.5 was used to calculate genetic diversity parameters, the number of alleles (Na), number of effective alleles (Ne), Shannon diversity index (I), haplotype diversity (h), and unbiased diversity (uh) (Peakall and Smouse, 2012). AMOVA was performed to assess the distribution of genetic variance among different cold resistance groups in YQ using GenAlEx 6.5. Plastid haplotypes were analyzed using Arlequin 3.5 (Excoffier and Lischer, 2010), and the haplotype network was constructed with NETWORK 10.14 using the Median Joining method. BAPS 6 was used with a spatial genetic mixture model to infer the geographic patterns of cpSSRs among different populations (Corander et al., 2008a,b). We ran the model for ten replicates with K ranging between 1 and 20 using the option of “clustering with linked loci” and chose the K with the highest log-likelihood. The presence of phylogeographic structure was determined using the ratio of NST (considering the mutational distances between haplotypes) and GST (depending on the frequencies of haplotypes) with 1,000 random permutations in Permut cpSSR 2.0 (Pons and Petit, 1996).
Results
Phenotypic Responses to Chilling Event
Based on the phenotypic response of individuals to chilling event in the YQ population during the winter of 2015–2016, we could classify them into three categories of cold resistance: the group with strong resistance to cold (SRC), the group with moderate resistance to cold (MRC) and the group with poor resistance to cold (PRC) (Figure 2). The results of PCA on five quantitative morphological traits also showed differentiation between the above three categories (Supplementary Figure 1). Further, we analyzed the traits of three categories in the YQ population in response to low temperature. The chilling damage to branches (estimated as the proportion of completely damaged branches) of the SRC group was 40%, which was much lower than the 89.7 and 100.0% estimated for MRC and PRC groups (Table 2). The defoliation percentage of SRC group was 57.1%, which was also much lower than the 96.1 and 100.0% estimated for the other two groups. However, certain “dead” individuals in the PRC group had new branches, and leaves grew out at the base of the trunk 10 months after the chilling event (Figure 2). Compared to the other two groups, the SRC group showed a relatively large variation of chilling damage to branches and defoliation percentage. Among all analyzed traits, the coefficient of variation was largest in plant height and stem basal diameter for all three cold-resistance groups (Table 2).
Figure 2. Criteria for classifying individuals into strong, moderate, and poor resistance to cold. Images of three groups of plants were taken after the chilling event in the winter of 2015–2016. The red arrow points to the newly grown branches of the “dead” individual during the second survey.
Table 2. Phenotypic traits of three groups with different cold resistance in the YQ population after the chilling event during the winter of 2015–2016.
Recent studies have shown that leaf size was often correlated with mean annual temperature (Moles et al., 2014). To test this, we measured leaf sizes of individuals from mixed population YQ and from the five natural populations. Regarding the leaf trait analysis, certain individuals in YQ did not have enough leaves for measurement. Eventually, 238 leaves were available for analysis (77 leaves from SRC, 59 leaves from MRC, and 102 leaves from PRC). We found that SRC had the smallest leaves and PRC had the largest (Table 2). In parallel, we measured all six phenotypic traits in five natural populations and similarly found a gradual decrease in leaf size along latitude (Supplementary Table 3).
Genetic Diversity, Variation and Structure Revealed by Nuclear Simple Sequence Repeats
Of the ten loci, the mean frequency of null alleles was 10%, ranging from 2.2 to 11.9%. Across the five natural populations, no signs of LD were found, while 10 out of 45 pairwise comparisons showed significant LD in YQ population (p < 0.01). Four nSSR loci (Kaca01, Kaca04, Kaca12, and Kcan004) deviated significantly from HWE (p < 0.01). One locus (Kcan011) was identified by BayeScan as an outlier and the remaining nine loci were considered neutral (Supplementary Table 4). In total, 176 alleles of 10 nuclear microsatellite loci were detected in six populations, and the number of alleles per locus (Na) ranged from 8 to 35, with a mean of 17.6 per locus. The genetic diversity parameters for each nSSR locus were shown in Supplementary Table 5.
Genetic diversity parameters of six populations are shown in Table 3. In general, genetic diversity of natural populations decreased with increasing latitude, except for population HK. As an artificial mixed population, various genetic diversity parameters of the artificial YQ population, including Ne, Ho, He, and I, were generally within the range of natural populations (Table 3). When Kcan011 was excluded, genetic diversity decreased slightly based on 9 neutral nSSR loci, but the pattern of genetic diversity among above populations did not change (Supplementary Table 6). In addition to the geographic grouping, three groups of K. obovata with different cold resistance in YQ population were analyzed based on 9 neutral nSSR loci and all nSSR loci separately. We found different levels of polymorphism at nuclear SSR loci: the MRC group had the highest genetic diversity, followed by the PRC group, and the SRC group displayed the lowest genetic diversity (Supplementary Tables 7, 8).
Table 3. Genetic diversities within each of the six populations of K. obovata based on ten nSSR loci and six cpSSR loci.
Genetic differentiation among six populations assessed by pairwise FST is listed in Table 4. The pairwise measure of differentiation in all populations was significant, but the degree of genetic differentiation among YX, SZ, and FCG populations was relatively low. AMOVA based on nuclear SSRs among five natural populations (Table 5) revealed that the genetic variation mainly existed within populations (72%), only 28% occurred among different populations (p < 0.001). Similarly, AMOVA based on nuclear microsatellites in YQ revealed that major genetic differences existed within different cold resistance groups (87%), and only 13% occurred among different groups (p < 0.001) (Supplementary Table 9).
Table 4. Population pairwise FST of the six populations of K. obovata.
Table 5. Summary of the AMOVA results for five natural populations of K. obovata.
The Bayesian clustering analysis for six populations using nine neutral loci indicated that K = 3 was the optimal parameter of their genetic structure (Figure 3A). These three groups were distributed among the six populations as follows: one group contained all individuals in FD and some individuals in YQ (there mainly, but not exclusively, in the SRC group); one group consisted of all individuals of populations HK and FCG, almost all individuals in SZ (with some admixture), and several admixed individuals YX; the last group was composed of few individuals in YX and most individuals in YQ (nearly all MRC and most PRC individuals). Consequently, there was no clear association between phenotypic and genetic variation, as phenotypic groups of different cold resistance in YQ did not correlate with genetic groups.
Figure 3. Estimated genetic structures of six populations based on nine nSSR loci (A) and on six cpSSR loci (B). (C) Plastid haplotype network of YQ and five natural populations, each circle represents a haplotype, and the color within the circle represents the populations or groups that share that haplotype. Black vertical lines between haplotype 32 and 33 represent mutation steps. See Table 1 for population abbreviations.
Haplotype Diversity and Geographic Patterns Inferred From Plastid Simple Sequence Repeats
The genetic polymorphisms at the six plastid microsatellite loci were relatively low (Supplementary Table 10). Only 24 alleles were detected across six cpSSR loci in all populations with Na varying from 2 to 8, the genetic diversity h ranging from 0.008 to 0.719, and Shannon’s index I ranging from 0.025 to 1.516. Regarding plastid microsatellite polymorphisms based on collection sites, genetic diversity decreases with increasing latitude in natural populations, and these genetic diversity parameters of YQ population were also within the range of the natural populations (Table 3). Interestingly, genetic diversity similarly decreases with increasing cold resistance in YQ population. The highest plastid polymorphism was observed in PRC (Na = 2.000, Ne = 1.577, h = 0.219, I = 0.384) and the lowest was in SRC (Na = 1.500, Ne = 1.201, h = 0.122, I = 0.200) (Supplementary Table 7). AMOVA based on cpSSR in YQ revealed that the plastid microsatellites variation was 63% among different cold resistance groups and 37% within groups (p < 0.001), consistent with the haplotype sharing between different groups (Supplementary Table 11).
A total of 15 plastid microsatellite haplotypes were found across three phenotypic groups of the artificial population (Figure 3C and Supplementary Table 12). Only one haplotype (H29) was shared between SRC and MRC, and the most frequent haplotypes were different among the three groups, indicating different origins of these three groups. For example, the most frequent haplotype of SRC was H35 (87%), and that of MRC was H29 (59%). In total, 36 haplotypes were found in six populations, of which six were shared by the natural populations and the artificial population in YQ (Figure 3C and Supplementary Table 12). Haplotypes H2 and H4 found in the southernmost natural population HK were shared with PRC. Haplotypes H22, H25, and H29 found in the medium latitude populations, including YX, SZ, and FCG, were also found in MRC in YQ. The most frequent haplotype H35 in SRC was primarily observed in the northernmost natural population of FD.
The population structure analysis based on cpSSRs identified 12 different clusters through BAPS software (K = 12 was the optimal parameter, log-likelihood values = −352.618) (Figure 3B). In detail, in the YQ population, groups with different cold resistance each occupied different clusters, and there were almost no shared clusters between different groups. In natural populations, FD and HK occupied their own unique clusters, while there was cluster sharing among YX, SZ, and FCG. Interestingly, SRC in YQ population shared some clusters with population FD, MRC shared the same clusters with three populations (YX, SZ, and FCG), and PRC shared some clusters with population HK.
The plastid haplotype network using the Median Joining methods indicated that the 36 haplotypes were grouped into two major clusters (Figure 3C), separated by sixteen mutational steps. Cluster I consisted of 19 haplotypes, and cluster II consisted of 17 haplotypes. Interestingly, MRC and PRC appeared in both clusters, whereas SRC only contained haplotypes belonging to cluster I. In addition, FD and HK populations hardly shared haplotypes with other natural populations, indicating significant phylogeographic structure across the five populations of K. obovata. The presence of phylogeographic structure was also confirmed by RST (0.308) being significantly higher than GST (0.217) (P < 0.001).
Discussion
Low temperature has been regarded as one of the primary factors affecting the latitudinal range limits of mangrove species (Lugo and Patterson-Zucca, 1977; Duke et al., 1998; Stuart et al., 2007; Ross et al., 2009; Tomlinson, 2016; Osland et al., 2017; Wu et al., 2018). Yang and Lin (1998) showed an increase of cold resistance within K. obovata along latitude. The chilling event that occurred during the winter of 2015–2016 in Yueqing provided us with an ideal opportunity to obtain data of phenotypic response of K. obovata to low temperature and to compare the cold-resistant capability between populations of different geographic origins.
Among several chilling injuries, defoliation is a typical phenotypic response found in mangroves (Ellis et al., 2006; Tomlinson, 2016). A previous study reported that plants of K. obovata in the Shacheng Bay (27°17′N, 120°18′E) displayed leaf scorch, defoliation, and bud withering after frost (Wang et al., 2011). However, another study showed no defoliation in K. obovata populations during the 2008 chilling event (Chen et al., 2017). In this study, we found that individuals of K. obovata in the artificially mixed population YQ could be divided into three different phenotypic groups according to their appearance of defoliation and bud withering (Figure 2). Furthermore, the reduced leaf length and width in the YQ population was associated with an increase of cold resistance, and the same trend was observed among natural populations along the latitudinal gradient, which may indicate that smaller leaf size was more adapted to low temperatures (Rosbakh et al., 2015). Nevertheless, the standard deviation and CV of phenotypic traits within each cold resistance type were relatively large, probably due to mixed individuals introduced from multiple natural populations since the 1950s.
Climate change can drive local adaptation of plants, especially trees (Petit and Hampe, 2006). The adaptive diversification of plant populations can be strongly shaped by historical and contemporary climatic factors (Turesson, 1925; Linhart and Grant, 1996; Etterson, 2004; Jordan et al., 2017; Cuervo-Alarcon et al., 2018). As an important selective agent, the temperature may affect the population adaptation of a species over altitude and latitude gradients (Saxe et al., 2001; Kitamura et al., 2020). Mangroves are temperature-sensitive plants, and thus low temperature may exert strong selective forces on natural populations. Temperature is likely a critical factor in determining mangrove community assemblage and species diversity (Yang and Lin, 1998; Wang et al., 2011; Ruan et al., 2013; Wu et al., 2018).
Despite that microsatellites are usually regarded as neutral genetic markers, they may also reflect population genetic composition shaped by natural selection (probably due to being linked to the selection-targeting loci), and thus, can be used for detecting ecological adaptation (Thibert-Plante and Hendry, 2010; Jay et al., 2012). Especially plastid genomes usually lack DNA recombination, so that neutral markers may be more likely linked to selection-targeting loci. Previous case studies, such as on oak, coast redwood and runner bean, demonstrated that plastid neutral makers clearly revealed genetic structure related to adaptation (Rodriguez et al., 2013; Yan et al., 2019; Breidenbach et al., 2020; Danusevičius et al., 2021). In the present study, the plastid microsatellites revealed genetic differentiation among three groups of different cold tolerance in YQ, the artificial population, and a significant geographical structure of haplotypes among natural populations along the latitudinal gradient. We found that the population at high latitude (FD) mainly shared haplotypes with the SRC group in YQ. Similarly, the population at low latitudes (HK) mainly shared haplotypes with the PRC phenotype group in YQ (Figure 3C and Supplementary Table 12). These results suggest that the cold resistance difference between SRC, MRC and PRC groups in the artificial populations may be explained by their origin from northern, mid-latitude and southern populations in nature. If plants in the artificial population represent an unbiased sample from the corresponding natural populations, our study indicates these natural populations are differentiated in cold resistance, which is likely the consequence of local adaption to different temperatures across latitudes.
Nevertheless, by using nuclear microsatellites we failed to detect a clear correlation between phenotypic and genetic patterns as revealed by the plastid markers. One reason could be that, concerning neutral genetic markers, plastid markers in this case may be more effective than nuclear markers in providing information of ecological adaptation due to the lack of DNA recombination in the plastome and thus stronger linkage to the selection-targeted loci. This assumes that at least some of the relevant loci are in the plastid genome, which requires further studies. A second reason may be the extent of gene flow. Previous studies reported that pollen dispersal distances are often long in tropical tree species, which depend on insects for pollination (Ward et al., 2005; Geng et al., 2008). In theory, natural selection is expected to decrease genetic diversity within a population while increasing the differentiations among populations. In contrast, strong gene flow may reduce divergence among populations while increasing variation within individual populations. Therefore, nuclear gene flow, including contributions from pollen flow over larger distances, may have blurred the imprint of natural selection on the population differentiation. Finally, another possible explanation is that we failed to collect the precise source populations of all individuals in the YQ population.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author Contributions
W-XL and S-CY conceived the idea, designed the experiments, and collected the plant materials. W-XL and B-HZ conducted the molecular work. W-XL performed the analyses and drafted the manuscript. W-XL, S-CY and Y-YZ revised the manuscript. All authors contributed to the study and approved the submitted version.
Funding
This study was funded by the National Key Research and Development Program of China (2017YFC0506103), Ocean Public Fund Research Project Grants of State Oceanic Administration (No. 201305021-2), and National Natural Science Foundation of China (No. 30972334).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors wish to thank Dr. Yan-Ping Guo (Beijing Normal University) for the helpful feedback on the manuscript. Additionally, the authors thank all reviewers for their comments and suggestions.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.695746/full#supplementary-material
Footnotes
- ^ http://www.tianqihoubao.com/lishi/leqing.html
- ^ https://www.r-project.org/
- ^ http://cmpg.unibe.ch/software/BayeScan/download.html
- ^ https://www.fluxus-engineering.com
References
Belkhir, K., Borsa, P., Chikhi, L., Raufaste, N., and Bonhomme, F. (1996-2004). Genetix 4.05, Logiciel Sous Windows TM Pour la Génétique des Populations. Laboratoire Génome, Populations, Interactions, CNRS UMR 5000. Montpellier: Université de Montpellier II.
Breidenbach, N., Gailing, O., and Krutovsky, K. V. (2020). Genetic structure of coast redwood (Sequoia sempervirens [D. Don] Endl.) populations in and outside of the natural distribution range based on nuclear and chloroplast microsatellite markers. PLoS One 15:e0243556. doi: 10.1371/journal.pone.0243556
Chapman, V. J. (1977). “Introduction,” in Ecosystems of the World, ed. V. J. Chapman (Amsterdam: Elsevier), 1–29. doi: 10.1093/actrade/9780198837572.003.0005
Chen, L. Z., Du, X. N., Lu, L. M., Zheng, S. J., and Yang, S. C. (2012). Synergistic effects of low temperature in winter and ebb tide at night on Sonneratia apetala seedlings growth and key eco-physiological traits (in Chinese). Chin. J. Appl. Ecol. 23, 953–958.
Chen, L. Z., Wang, W. Q., Li, Q. Q., Zhang, Y. H., Yang, S. C., Osland, M. J., et al. (2017). Mangrove species’ responses to winter air temperature extremes in China. Ecosphere 8:e01865. doi: 10.1002/ecs2.1865
Chen, Q. X., Yang, S., Wang, J. W., Liu, X., Zheng, J., and Deng, R. J. (2019). Development history and discussion of mangrove forest in Zhejiang Province (in Chinese). J. Zhejiang Agric. Sci. 60, 1177–1181.
Corander, J., Marttinen, P., Sirén, J., and Tang, J. (2008a). Enhanced Bayesian modelling in BAPS software for learning genetic structures of populations. BMC Bioinformatics 9:539. doi: 10.1186/1471-2105-9-539
Corander, J., Sirén, J., and Arjas, E. (2008b). Bayesian spatial modeling of genetic population structure. Comput. Stat. 23, 111–129. doi: 10.1007/s00180-007-0072-x
Cuervo-Alarcon, L., Arend, M., Mueller, M., Sperisen, C., Finkeldey, R., and Krutovsky, K. V. (2018). Genetic variation and signatures of natural selection in populations of European beech (Fagus sylvatica L.) along precipitation gradients. Tree Genet. Genomes 14:84. doi: 10.1007/s11295-018-1297-2
Danusevičius, D., Baliuckas, V., Buchovska, J., and Kembrytë, R. (2021). Geographical structuring of Quercus robur (L.) chloroplast DNA haplotypes in Lithuania: recolonization, adaptation, or overexploitation effects? Forests 12:831. doi: 10.3390/f12070831
Duke, N. C., Ball, M. C., and Ellison, J. C. (1998). Factors influencing biodiversity and distributional gradients in mangroves. Glob. Ecol. Biogeogr. Lett. 7, 27–47. doi: 10.2307/2997695
Eckert, A. J., Bower, A. D., Wegrzyn, J. L., Pande, B., Jermstad, K. D., Krutovsky, K. V., et al. (2009). Association genetics of coastal Douglas fir (Pseudotsuga menziesii var. menziesii, Pinaceae). I. Cold-hardiness related traits. Genetics 182, 1289–1302. doi: 10.1534/genetics.109.102350
Ellis, W. L., Bowles, J. W., Erickson, A. A., Stafford, N., Bell, S. S., and Thomas, M. (2006). Alteration of the chemical composition of mangrove (Laguncularia racemosa) leaf litter fall by freeze damage. Estuar. Coast. Shelf Sci. 68, 363–371. doi: 10.1016/j.ecss.2006.02.017
Etterson, J. R. (2004). Evolutionary potential of Chamaecrista fasciculata in relation to climate change. 1. Clinal patterns of selection along an environmental gradient in the Great Plains. Evolution 58, 1446–1458. doi: 10.1111/j.0014-3820.2004.tb01726.x
Evanno, G., Regnaut, S., and Goudet, J. (2005). Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol. Ecol. 14, 2611–2620. doi: 10.1111/j.1365-294X.2005.02553.x
Excoffier, L., and Lischer, H. E. L. (2010). Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour. 10, 564–567. doi: 10.1111/j.1755-0998.2010.02847.x
Geng, Q., Lian, C., Goto, S., Tao, J., Kimura, M., Isalm, M. S., et al. (2008). Mating system, pollen and propagule dispersal, and spatial genetic structure in a high-density population of the mangrove tree Kandelia candel. Mol. Ecol. 17, 4724–4739. doi: 10.1111/j.1365-294X.2008.03948.x
Halbritter, A. H., Fior, S., Keller, I., Billeter, R., Edwards, P. J., Holderegger, R., et al. (2018). Trait differentiation and adaptation of plants along elevation gradients. J. Evol. Biol. 31, 784–800. doi: 10.1111/jeb.13262
Hu, H. Y., Chen, S. Y., Wang, W. Q., Dong, K. Z., and Lin, G. H. (2012). Current status of mangrove germplasm resources and key techniques for mangrove seedling propagation in China (in Chinese). Chin. J. Appl. Ecol. 23, 939–946.
Islam, M. S., Lian, C., Geng, Q., Kameyama, N., and Hogetsu, T. (2008). Chloroplast microsatellite markers for the mangrove tree species Bruguiera gymnorrhiza, Kandelia candel, and Rhizophora stylosa, and cross-amplification in other mangrove species. Conserv. Genet. 9, 989–993. doi: 10.1007/s10592-007-9424-x
Islam, M. S., Tao, J. M., Geng, Q. F., Lian, C. L., and Hogetsu, T. (2006). Isolation and characterization of eight compound microsatellite markers in a mangrove tree Kandelia candel (L.) Druce. Mol. Ecol. Notes 6, 1111–1113. doi: 10.1111/j.1471-8286.2006.01453.x
Jay, F., Manel, S., Alvarez, N., Durand, E. Y., Thuiller, W., Holderegger, R., et al. (2012). Forecasting changes in population genetic structure of alpine plants in response to global warming. Mol. Ecol. 21, 2354–2368. doi: 10.1111/j.1365-294X.2012.05541.x
Jordan, R., Hoffmann, A. A., Dillon, S. K., and Prober, S. M. (2017). Evidence of genomic adaptation to climate in Eucalyptus microcarpa: implications for adaptive potential to projected climate change. Mol. Ecol. 26, 6002–6020. doi: 10.1111/mec.14341
Jump, A. S., Hunt, J. M., Martinez-Izquierdo, J. A., and Penuelas, J. (2006). Natural selection and climate change: temperature-linked spatial and temporal trends in gene frequency in Fagus sylvatica. Mol. Ecol. 15, 3469–3480. doi: 10.1111/j.1365-294X.2006.03027.x
Kawecki, T. J., and Ebert, D. (2004). Conceptual issues in local adaptation. Ecol. Lett. 7, 1225–1241. doi: 10.1111/j.1461-0248.2004.00684.x
Kitamura, K., Uchiyama, K., Ueno, S., Ishizuka, W., Tsuyama, I., and Goto, S. (2020). Geographical gradients of genetic diversity and differentiation among the southernmost marginal populations of Abies sachalinensis revealed by EST-SSR polymorphism. Forests 11. doi: 10.3390/f11020233
Lenormand, T. (2002). Gene flow and the limits to natural selection. Trends Ecol. Evol. 17, 183–189. doi: 10.1016/S0169-5347(02)02497-7
Li, J. Q., Xu, H. F., Ye, Z. L., Gu, J. F., and Wang, R. Q. (2001). Introduction and afforestation technique of Kandelia candel to the north (in Chinese). J. Zhejiang For. Sci. Technol. 21, 51–53.
Li, M. S., and Lee, S. Y. (1997). Mangroves of China: a brief review. For. Ecol. Manag. 96, 241–259. doi: 10.1016/s0378-1127(97)00054-6
Lin, P., Shen, R. C., and Lu, C. Y. (1994). The characteristics of chilling-resistance on six mangrove plants (in Chinese). J. Xiamen Univ. 33, 249–252.
Linhart, Y. B., and Grant, M. C. (1996). Evolutionary significance of local genetic differentiation in plants (in Chinese). Annu. Rev. Ecol. Syst. 27, 237–277. doi: 10.1146/annurev.ecolsys.27.1.237
Lomolino, M., Riddle, B., Whittaker, R., and Brown, J. (2010). Biogeography. Sunderland, MA: Oxford University Press.
Lovelock, C. E., Krauss, K. W., Osland, M. J., Reef, R., and Ball, M. C. (2016). “The physiology of mangrove trees with changing climate,” in Tropical Tree Physiology, eds G. Goldstein and L. S. Santiago (Cham: Springer), 149–179.
Lugo, A. E., and Patterson-Zucca, C. (1977). The impact of low temperature stress on mangrove structure and growth. Trop. Ecol. 18, 149–161.
MacArthur, R. H. (1984). Geographical Ecology: Patterns in the Distribution of Species. Princeton, NJ: Princeton University Press.
Moles, A. T., Perkins, S. E., Laffan, S. W., Flores-Moreno, H., Awasthy, M., Tindall, M. L., et al. (2014). Which is a better predictor of plant traits: temperature or precipitation? J. Veg. Sci. 25, 1167–1180. doi: 10.1111/jvs.12190
Mu, A. Y., Liu, J. J., Yang, J. Q., Teng, Y. P., and Zhou, C. Q. (2005). Investigation of introduction mangroves of Kandelia candel Druce in the northest boundary of China (in Chinese). Prot. For. Sci. Technol. 1, 6–8.
Osland, M. J., Feher, L. C., Griffith, K. T., Cavanaugh, K. C., Enwright, N. M., Day, R. H., et al. (2017). Climatic controls on the global distribution, abundance, and species richness of mangrove forests. Ecol. Monogr. 87, 341–359. doi: 10.1002/ecm.1248
Peakall, R., and Smouse, P. E. (2012). GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research-an update. Bioinformatics 28, 2537–2539. doi: 10.1093/bioinformatics/bts460
Peng, Y. L., Wang, Y. S., Fei, J., Sun, C. C., and Cheng, H. (2015). Ecophysiological differences between three mangrove seedlings (Kandelia obovata, Aegiceras corniculatum, and Avicennia marina) exposed to chilling stress. Ecotoxicology 24, 1722–1732. doi: 10.1007/s10646-015-1488-7
Petit, R. J., and Hampe, A. (2006). Some evolutionary consequences of being a tree. Annu. Rev. Ecol. Evol. Syst. 37, 187–214. doi: 10.1146/annurev.ecolsys.37.091305.110215
Pons, O., and Petit, R. J. (1996). Measuring and testing genetic differentiation with ordered versus unordered alleles. Genetics 144, 1237–1245. doi: 10.1093/genetics/144.3.1237
Pritchard, J. K., Stephens, M., and Donnelly, P. (2000). Inference of population structure using multilocus genotype data. Genetics 155, 945–959. doi: 10.1093/genetics/155.2.945
Raymond, M., and Rousset, F. (1995). GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J. Hered. 86, 248–249. doi: 10.1093/oxfordjournals.jhered.a111573
Rodriguez, M., Rau, D., Angioi, S. A., Bellucci, E., Bitocchi, E., Nanni, L., et al. (2013). European Phaseolus coccineus L. landraces: population structure and adaptation, as revealed by cpSSRs and phenotypic analyses. PLoS One 8:e57337. doi: 10.1371/journal.pone.0057337
Rosbakh, S., Römermann, C., and Poschlod, P. (2015). Specific leaf area correlates with temperature: new evidence of trait variation at the population, species and community levels. Alp. Bot. 125, 79–86. doi: 10.1007/s00035-015-0150-6
Ross, M. S., Ruiz, P. L., Sah, J. P., and Hanan, E. J. (2009). Chilling damage in a changing climate in coastal landscapes of the subtropical zone: a case study from south Florida. Glob. Change Biol. 15, 1817–1832. doi: 10.1111/j.1365-2486.2009.01900.x
Ruan, Y., Huang, B. H., Lai, S. J., Wan, Y. T., Li, J. Q., Huang, S., et al. (2013). Population genetic structure, local adaptation, and conservation genetics of Kandelia obovata. Tree Genet. Genomes 9, 913–925. doi: 10.1007/s11295-013-0605-0
Sani, S. G. A. S., Chang, P. L., Zubair, A., Carrasquilla-Garcia, N., Cordeiro, M., Penmetsa, R. V., et al. (2018). Genetic diversity, population structure, and genetic correlation with climatic variation in chickpea (Cicer arietinum) landraces from Pakistan. Plant Genome 11:170067. doi: 10.3835/plantgenome2017.08.0067
Saxe, H., Cannell, M. G. R., Johnsen, B., Ryan, M. G., and Vourlitis, G. (2001). Tree and forest functioning in response to global warming. New Phytol. 149, 369–399. doi: 10.1046/j.1469-8137.2001.00057.x
Sork, V. L. (2016). Gene flow and natural selection shape spatial patterns of genes in tree populations: implications for evolutionary processes and applications. Evol. Appl. 9, 291–310. doi: 10.1111/eva.12316
Stuart, S. A., Choat, B., Martin, K. C., Holbrook, N. M., and Ball, M. C. (2007). The role of freezing in setting the latitudinal limits of mangrove forests. New Phytol. 173, 576–583. doi: 10.1111/j.1469-8137.2006.01938.x
Su, W., Ye, C., Zhang, Y., Hao, S., and Li, Q. Q. (2019). Identification of putative key genes for coastal environments and cold adaptation in mangrove Kandelia obovata through transcriptome analysis. Sci. Total Environ. 681, 191–201. doi: 10.1016/j.scitotenv.2019.05.127
Sugaya, T., Takeuchi, T., Yoshimaru, H., and Katsuta, M. (2002). Development and polymorphism of simple sequence repeat DNA markers for Kandelia candel (L.) Druce. Mol. Ecol. Notes 2, 65–66. doi: 10.1046/j.1471-8286.2002.00151.x
Thibert-Plante, X., and Hendry, A. P. (2010). When can ecological speciation be detected with neutral loci? Mol. Ecol. 19, 2301–2314. doi: 10.1111/j.1365-294X.2010.04641.x
Turesson, G. (1925). The plant species in relation to habitat and climate: contributions to the knowledge of genecological units. Hereditas 6, 147–236. doi: 10.1111/j.1601-5223.1925.tb03139.x
Wang, W., You, S., Wang, Y., Huang, L., and Wang, M. (2011). Influence of frost on nutrient resorption during leaf senescence in a mangrove at its latitudinal limit of distribution. Plant Soil 342, 105–115. doi: 10.1007/s11104-010-0672-z
Ward, M., Dick, C. W., Gribel, R., and Lowe, A. J. (2005). To self, or not to self…A review of outcrossing and pollen-mediated gene flow in neotropical trees. Heredity 95, 246–254. doi: 10.1038/sj.hdy.6800712
Wu, Y., Ricklefs, R. E., Huang, Z., Zan, Q., and Yu, S. (2018). Winter temperature structures mangrove species distributions and assemblage composition in China. Glob. Ecol. Biogeogr. 27, 1492–1506. doi: 10.1111/geb.12826
Yan, M., Liu, R., Li, Y., Hipp, A. L., Deng, M., and Xiong, Y. (2019). Ancient events and climate adaptive capacity shaped distinct chloroplast genetic structure in the oak lineages. BMC Evol. Biol. 19:202. doi: 10.1186/s12862-019-1523-z
Yang, S. C., Li, Y., and Lin, P. (2003). Change of leaf caloric value from Avicennia marina and Aegiceras corniculatum mangrove plants under cold stress. J. Oceanogr. Taiwan Strait 22, 46–52.
Yang, S. C., and Lin, P. (1997). Cold-resistance ability of two mangrove species Kandelia candel and Aegiceras corniculatum during their overwintering period (in Chinese). Chin. J. Appl. Ecol. 8, 561–565.
Yang, S. C., and Lin, P. (1998). Ecological studies on the resistance and adaptation to cold of some tidal mangrove species in China (in Chinese). Chin. J. Plant Ecol. 22, 60–67.
Yang, S. C., Lin, P., Li, Z. J., Wang, W. W., and Tsuneo, N. (1999). Effect of low night temperature on photosynthetic properties of mangrove seedlings. J. Xiamen Univ. 38, 617–621.
Zheng, C., Qiu, J., Liu, W., Huang, L., Chen, S., Huang, X., et al. (2012). Ecophysiological characteristics of higher-latitude transplanted mangrove Kandelia candel in strong tidal range area. Acta Ecol. Sin. 32, 4453–4461. doi: 10.5846/stxb201106190843
Keywords: adaptation, Kandelia obovata, low temperature, mangroves, genetic structure
Citation: Lu W-X, Zhang B-H, Zhang Y-Y and Yang S-C (2022) Differentiation of Cold Tolerance in an Artificial Population of a Mangrove Species, Kandelia obovata, Is Associated With Geographic Origins. Front. Plant Sci. 12:695746. doi: 10.3389/fpls.2021.695746
Received: 15 April 2021; Accepted: 31 December 2021;
Published: 03 February 2022.
Edited by:
Gerald Matthias Schneeweiss, University of Vienna, AustriaReviewed by:
Alejandra Vázquez-Lobo, Universidad Autónoma del Estado de Morelos, MexicoXiao-Yong Chen, East China Normal University, China
Lian Chunlan, The University of Tokyo, Japan
Copyright © 2022 Lu, Zhang, Zhang and Yang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sheng-Chang Yang, scyang@xmu.edu.cn
†These authors have contributed equally to this work and share first authorship