Population genetic structure and variability in Lindera glauca (Lauraceae) indicates low levels of genetic diversity and skewed sex ratios in natural populations in mainland China

Lindera glauca (Lauraceae) is a tree of economic and ecological significance that reproduces sexually and asexually via apomictic seeds. It is widely distributed in the low-altitude montane forests of East Asia. Despite the potential implications of a mixed reproductive system in terms of genetic diversity, few studies have focused on this aspect. In this study, the genetic structure of wild populations of L. glauca was investigated via genetic analyses. Overall, 13 nuclear microsatellites (nSSRs) and five chloroplast microsatellites (cpSSRs) were used to genotype 300 individual plants, taken from 20 wild populations (a small sample size in some wild populations is due to the limitation of its specific reproduction, leading to certain limitations in the results of this study) and two cultivated populations ranging across nearly the entire natural distribution of mainland China. The populations exhibited low levels of genetic diversity (nSSR: AR = 1.75, Ho = 0.32, He = 0.36; cpSSR: Nb = 2.01, Hrs = 0.40), and no significant effect of isolation by distance between populations existed, regardless of marker type (nSSR: R2 = 0.0401, P = 0.068; cpSSR: R2 = 0.033, P = 0.091). Haplotype networks showed complex relationships among populations, and the H12 haplotype was predominant in most populations. Analyses of molecular variance obtained with nuclear markers (Fsc = 0.293, FST = 0.362) and chloroplast markers (Fsc = 0.299, FST = 0.312) were similar. The migration ratio of pollen flow versus seed flow in this study was negative (r = −1.149). Results suggest that weak barriers of dispersal between populations and/or the similarity of founders shared between neighbors and distant populations are indicative of the gene flow between populations more likely involving seeds. Wild L. glauca in mainland China was inferred to have highly skewed sex ratios with predominant females. In addition, some populations experienced a recent bottleneck effect, especially in Gujianshan, Chongqing, and southwest China (population GJS). It is suggested that few wild male individuals should be conserved in order to maintain overall genetic diversity in the wild populations of this species. These findings provide important information for the sustainable utilization and preservation of the overall genetic diversity of L. glauca.


INTRODUCTION
Plant populations respond to the changing environment and climate via phenotype shifts (Nicotra et al., 2010) and, mainly through sexual reproduction, by bringing together high-fitness alleles that reside in different individuals (Whitton et al., 2008). In general, sexual reproduction is predominant in eukaryotes and a nearly universal characteristic of angiosperms. In some groups, however, some sexual plants can reproduce asexually via apomixis, which is the production of clonal seeds in the absence of fertilization (Richards, 1986;Daniel et al., 2001), producing exact genetic replicas of maternal plants (Daniel et al., 2001). Apomixis occurs in fewer than 1% of flowering plant species, with an uneven distribution among lineages (Whitton et al., 2008). Apomixis occurs sporadically among angiosperms (APG, 2003) (Whitton et al., 2008). In some genera (i.e., Taraxacum), apomictic clones can be widely distributed and are temporarily ecologically successful (Van Dijk, 2003;Majesky et al., 2012). However, lack of diversity, the limited possibility of acquiring heritable variability (Richards, 1996), and an increased mutation load leading to the extinction of clones (Van Dijk, 2003), give apomicts an adaptive disadvantage. In contrast, apomicts have lower reproductive costs over than sexuals, a high proportion of loci fixed to heterozygous conditions, and significant advantages over sexuals in colonizing new areas (Majesky et al., 2012). Due to these short-term advantages, natural populations of apomicts are of interest for agricultural development.
There are two major types of apomixis, adventitious embryony and gametophytic apomixis that differ in the way embryos are formed (Whitton et al., 2008;Lo, Stefanovic & Dickinson, 2009). The origin of the former is somatic tissue surrounding the fertilized ovule, and the origin of the latter is an unreduced megagametophyte. Adventitious embryony is widely distributed in nature, and gametophytic apomixis is reported in a few families, e.g., Asteraceae, Poaceae, and Rosaceae. There are two ways for these to spread across space, direct dispersal via apomictic seeds, and indirect transmission via pollen (Whitton et al., 2008). For indirect transmission via pollen, the genes for maternal clonality can be transmitted via male gametes, and this mode of transmission may well be important in the establishment and spread of apomixis (Brock, 2004;Preite et al., 2015). Therefore, the transmission of apomixis genes to sexuals via pollen may be of long-term importance for the spread of apomixis, especially for an agriculturally important tree such as Lindera glauca.
Lindera glauca (Sieb. et Zucc.) Blume (Lauraceae), a deciduous shrub or small tree with both apomixis (asexual reproduction by seeds) (Dupont, 2002) and a sexual reproduction system (Tsui, Xia & Li, 1982;Tsui & Werff, 2008), is distributed extensively in low-altitude montane forests of central and southern mainland China, as well as in Japan, Korea, Vietnam and Taiwan (Wang, 1972;Chang, 1976;Zheng, 1983). As one of the main trees making up the shrubbery and young forest ecosystems in the central and southern areas of mainland China, this species has both economic value and ecological importance. Its fruits are rich in fatty acids and aromatic oils, and they contain terpenoids, flavonoids, and alkaloids, which are used for various applications in traditional medicine. Fruits are also used as raw materials to produce medicines, lubricants, and biochemical products (Zheng, 1983;Wang et al., 1994;Kim et al., 2014;Suh et al., 2015;Qi et al., 2016). Some root extract components, like N-methyllaurotetanine, exhibit significant anti-tumor metastatic activity (Kim et al., 2014;Suh et al., 2015) and some volatile oils from the leaves are used in the industrial production pf spices (Qi et al., 2016). Additionally, the L. glauca species has emerged as a novel potential source of biodiesel in China due to the high quality and quantity of its fruit oil (Lin et al., 2017;Xiong, Dong & Zhang, 2018). There has been increased scientific interest in the species, but relatively little remains known about its reproductive modes and their potential effects on genetic diversity in population dynamics and population differentiation.
L. glauca is native to the mainland China, and diploids (2n = 24) (Yang, 1999) with sexual reproduction and male plants have been known to exist in continental East Asia for several decades (Wang, 1972;Tsui, Xia & Li, 1982;Tsui & Werff, 2008). In a study conducted in Japan, Dupont (2002) found that female L. glauca could asexually reproduce via seeds. Adult population sex ratios of other Lindera species observed in Japan ranged from equal to a strong male bias (including L. obtusiloba, L. umbellata, and L. erythrocarpa) (Dupont, 2002). However, recent empirical studies revealed how L. glauca males are very rare in mainland China, with females reproducing via apomixis. This indicates that natural populations have a mixed reproduction mode that includes apomixis and sexual propagation. Apomixis might play a major role in shaping the genetic structure of the species, by limiting gene flow within populations (Daniel et al., 2001). Interpopulation gene flow in plants is mediated by a combination of pollen and seed dispersal (Ennos, 1994). Some natural populations of apomicts retain residual sexual function as pollen donors and thus have the potential to spread apomixis via male gametes, thereby increasing the genetic diversity observed within apomictic populations (Whitton et al., 2008). In previous studies and records (Wang, 1972;Chang, 1976;Zheng, 1983;Dupont, 2002;Tsui & Werff, 2008), L. glauca is dioecious, and has bisexual or functionally unisexual flowers. However, our survey indicates that there was a very small amount of pollen from the staminode of female flowers (By 2, 3, 5-Triphenyltetrazolium chloride, or TTC method) (Hu, 1993), implying that they have the potential for natural pollination.
Thus, the genetic diversity and structure of natural populations of L. glauca may well be more complex than previously thought. It is essential to study the gene flow and estimate the relative rates of pollen and seed migration among natural populations. Furthermore, population bottleneck effect is thought to be responsible for the very low levels of genetic variation found in a number of species that now have large population sizes (Pannell, 2013). Given that there are very few males of L. glauca in mainland China over the last decade, and many males grown on Taiwan (Zhang, 2007), it is interesting to investigate whether natural populations in mainland China experienced a bottleneck effect or not?
In plants, organelle genomes are often uniparentally inherited (Wills et al., 2005). It is widely believed that plastid genomes are inherited through the maternal parent (Zhang & Sodmergen, 2010). Recently strict maternal inheritance of the plastid was observed in around 82% of the species in two large-scale studies totalling over 500 angiosperms (Corriveau & Coleman, 1988;Zhang, Liu & Sodmergen, 2003;Zhang & Sodmergen, 2010). Strict paternal inheritance is rare, and population size matters as the levels of paternal transmission can be as low as 0.03% in determining the mode of inheritance of the chloroplast genome (Wang et al., 2004). Thus, considering the actual sample size per population, we assume that chloroplast DNA is maternal inheritance for L. glauca.
In the present study, the aims were to (1) investigate the genetic diversity of L. glauca populations in the mainland China, (2) detect genetic variation within and differentiation among natural populations, (3) assess the relative importance of pollen and seeds as agents of gene flow, and (4) determine whether natural populations experienced a decline in size (bottleneck effect). Molecular genetic analyses were performed and individuals in 20 wild populations (and two cultivated populations) of L. glauca were genotyped using 13 nuclear and five chloroplast microsatellite markers developed in the previous work (Xiong et al., 2016;Xiong et al., 2018).

Sample collection
During field expeditions carried out from 2013 to 2017, 300 individuals were sampled from 20 wild populations and two cultivated populations, representing nearly the entire natural distribution of L. glauca in mainland China (Table 1; Fig. 1A). Most L. glauca individuals are able to form clones via vegetative reproduction with stolons (Tsui, Xia & Li, 1982), as found in our field survey. In order to avoid the collection of several ramets from the same genet, a single sample was obtained from each cluster of shrubs in close proximity to a main tree, excluding the surrounding young branches growing on the ground. Each sample (individual) in same population was collected at least 10 m apart. In some smaller populations, fewer than 10 plants of putatively nonclonal origin were available. Overall sample sizes varied from 5 to 30 per wild population (Table 1). There was definite apomixis in the individuals of two cultivated populations (based on our survey for three years), and we added it to global analysis as a reference for wild populations. In the field, fresh leaves were immediately dried in silica gel after collection, and preserved at room temperature until DNA extraction.

DNA extraction and microsatellite genotyping
Genomic DNA was extracted from 100-150 mg of dried leaves per sample using a modified cetyltrimethylammonium bromide (CTAB) method (Doyle & Doyle, 1987). Microsatellite loci of all 300 individuals of L. glauca were screened, including 13 polymorphic nuclear microsatellite markers (EST-based microsatellites; hereafter nSSRs; Table S1) and five polymorphic chloroplast microsatellite markers (hereafter cpSSRs; Table S2). All nSSRs and cpSSRs were labeled with fluorescently labeled nucleotides (forward primer with M13F) and detected by capillary gel electrophoresis. Subsequent\ steps and the PCR assay were conducted according to Xiong et al. (2016). Genotyping was performed using an ABI 3730XL DNA Analyzer (Applied Biosystems, California, USA ) with a GeneScan 500 LIZ Size Standard, and alleles for each locus were manually scored using GeneMarker version 2.2.0 software (SoftGenetics, State College, PA, USA). Notes. a cultivated population. SZY cultivars were introduced from Jiangsu between 1973-1977. HZY cultivars were introduced from Guangdong in before 1985; all Samples accession numbers refer to voucher specimens deposited in the Beijing Forestry University (BJFU); geographic coordinates and elevation were obtained with portable GPS receiver.

Data analyses
Raw data matrices containing information of alleles and haplotypes for 13 nSSR and 5 cpSSR loci were checked for scoring errors. All SSR analyses were conducted with 300 samples. Data editing and formatting were performed using GenAlEx v. 6.502 (Smouse, Whitehead & Peakall, 2015). The related indexes of genetic diversity were calculated as followed: for the nSSR data set, genetic diversity indices, including the number of alleles (N A ), observed heterozygosity (Ho), expected heterozygosity (He), percentage of polymorphic loci (PPB), Wright's inbreeding coefficient (F IS ), and Nei's (Nei, 1978) genetic distances, were estimated using GenAlEx v. 6.502 (Smouse, Whitehead & Peakall, 2015) and POPGENE v. 1.32 (Yeh, Yang & Boyle, 1999). The online package GENEPOP v. 4.1.4 (Rousset, 2008) was used to perform exact Hardy-Weinberg equilibrium (HWE) tests and to test for the presence of private (null) alleles. The differentiation index F ST was computed for pairs of populations using Arlequin v. 3.5.1.3 (Excoffier, Laval & Schneider, 2005). Allelic richness (A R ) was calculated using the software FSTAT v. 1.2 (Goudet, 1995). For the cpSSR data set, number of  Table 1 for details). Each population is represented by a triangle, and pie charts are shown when a population was present in more than one haplotype. The green background shows the provincial-level distribution of the species in China. (B) Haplotype network generated with the TCS program. Each haplotype is represented by a single color, and circle sizes correspond to the relative frequency of a particular haplotype in the total sample.
Full-size DOI: 10.7717/peerj.8304/ fig-1 haplotypes (Nb), genetic diversity (Dv), haplotype richness (Hrs), the number of private alleles (Prv), and the polymorphism information content (PIC) per locus were estimated using HAPLOTYPE v. 1.05 (Eliades & Eliades, 2009). The population genetic structure was analyzed as followed: for the nSSR data set, the genetic structure of the 22 populations (20 wild and 2 cultivated) was analyzed using the Bayesian clustering approach implemented in STRUCTURE v. 2.3.4 (Pritchard, Stephens & Donnelly, 2000), assuming an admixture model. In order to determine the most appropriate number of genetic clusters or groups (K value), K was set from 1 to 20, and the analysis was run with 20 iterations for each K with a burn-in of 1,000,000 generations followed by 50,000 generations for the Markov chain Monte Carlo (MCMC) simulation. The admixture level for each individual (Q) was also inferred. The program STRUCTURE HARVESTER v. 0.6.94 (Earl & Von Holdt, 2012) was used to estimate the number of population clusters based on the K parameter according to Evanno, Regnaut & Goudet (2005). Based on the most appropriate number of clusters suggested by Bayesian clustering, analysis of molecular variance (AMOVA) was performed using Arlequin, with 10,000 iterations for the permutation test. A neighbor-joining (NJ) tree was generated using POPTREE2 (Takezaki, Nei & Tamura, 2010) based on pairwise Nei (1978) genetic distances between populations determined by GenAlEx. For the cpSSR data set, the Arlequin was used to determine pairwise F ST values among all populations. A parsimony network illustrating genetic relationships among haplotypes of L. glauca populations was generated using TCS v.1.1 (Clement, Posada & Crandall, 2000).
Considering that isolation by distance (IBD) can be a key factor keeping populations apart by limiting gene flow (Coyne & Orr, 2004), the IBD of wild L. glauca inter-population in mainland China was tested. In view of the potential importance of pollen in the spread of apomixis (Mogie, 1992;Whitton et al., 2008), the pollen/seed migration ratio (r) was calculated. In order to examine IBD, the Mantel test was performed using GenAlEx, correlating the pairwise genetic distances [F ST /(1 − F ST )] with the pairwise geographic distances (in kilometers). To calculate r, we used the followed formula: (Ennos, 1994;Petit et al., 2005), where mp is the pollen migration rate, ms is the seed migration rate, Population bottlenecks were evaluated using BOTTLENECK v. 1.2.02 (Piry, Luikart & Cornuet, 1990) with the infinite alleles model (IAM) that a single mutation is allocated at a time and the resulting number of alleles is computed, stepwise mutation model (SMM) that is a Bayesian approach and generally more appropriate when testing microsatellite loci, and two-phased model (TPM) that is a modified SMM. According to Piry et al., sign tests, Wilcoxon tests, and mode-shift were applied, excluding standardized differences tests, which are useful when at least 20 polymorphic loci are available.

Genetic variation
A total of 74 alleles and 13 haplotypes were identified at 13 nSSRs and 5 cpSSRs across 300 individuals of L. glauca. For each locus, the number of alleles for 13 nSSR loci ranged from 3 (P-298) to 8 (XBLG-060), with a mean of 5.7 alleles (Table S1). In particular, A R ranged from 1.807 to 2.774, with a mean of 2.329, and PIC ranged from 0.363 to 0.711, with a mean of 0.556. Ho and He varied between 0.210 and 0.563, with a mean of 0.380, and between 0.380 and 0.754, with a mean of 0.602.
The total number of alleles of 20 wild populations across 13 nSSR loci varied from 22 (population SZY) to 44 (population ZJJ), with a mean of 32.9, and allelic richness ranged from 1.231 to 2.011, with a mean of 1.740 (Table 2). Population Ho ranged from 0.108 to 0.708 and population He from 0.106 to 0.477, with means of 0.328 and 0.376. The PPB ranged from 53.85% to 100%, with a mean of 85% (Table 2) (Table S1), and there is no evidence of private alleles within the data set. Among all wild samples (290), 277 individual plants had a unique multi-locus pattern after PCR amplification with 13 nSSRs primers, indicating that these samples were from different individuals. Of the remaining 13 individual plants, 5 pairs exhibited the same multi-locus pattern in pairs, and 3 individuals exhibited the same multi-locus pattern.
Genetic diversity parameters for cpSSR loci are summarized in Table S2. All 5 cpSSR loci exhibited 2-3 alleles per locus across all samples. Dv ranged from 0.115 to 0.249 per locus, and PIC varied from 0.155 to 0.218. Analyzing combinations of all alleles, there were 22 unique haplotypes (hereafter H; Fig. 1). All populations contained several haplotypes, except for populations ATM, YTH, NTB, FHS, TMS, and HZY (Fig. 1A). The network of plastid haplotypes was complex (Fig. 1B). Haplotype H12 exhibited the highest frequency and was detected in 18 of 20 wild populations (Table S3). Of the 22 haplotypes, 7 were private haplotypes excluding cultivated populations, Hrs per population ranged from 0 (populations ATM, NTB, YTH, TMS, and FHS) to 0.964 (population FJS), and the mean genetic distance between individuals (D 2 sh ) varied from 0 to 39.638 (population JGS) ( Table 2).

Genetic clustering and population differentiation
A Bayesian analysis based on 13 nSSRs implemented in STRUCTURE showed the presence of 2 clusters (K = 2), with only slight admixture at the individual level in each population, except for population ATM (Fig. 2A). The K statistic developed by Evanno et al. indicated that the overall differences were not substantial (Fig. 2B). Cluster orange included 13 wild populations (ATM, JGS, LDZ, SJG, NTB, YTH, DBS, HMF, TMS, SQS, LYS, KYS, and WYS), and the remaining 7 wild populations and 2 cultivated populations were assigned to cluster blue. The NJ tree (Fig. 3A) and principal coordinates analysis (Fig. 3B) based on the nSSR dataset supported the results of STRUCTURE analysis, indicating that 22 populations could be grouped into 2 clusters. However, the network diagram of all 22 unique plastid haplotypes revealed by 5 cpSSRs was complex (Fig. 1B), and failed to support the 2 distinct clusters revealed using nuclear data.
The 2 clusters revealed by STRUCTURE analysis were set as groups for AMOVA based on both the nSSR and cpSSR dataset (Table 3). Using the nSSR dataset, the majority of genetic variation was detected within populations (63.82%), indicating a genetic differentiation mostly at the individual level. Nevertheless, a considerable proportion of the total variation (26.42%) was found among populations within groups, and a small amount of variation  (9.76%) occurred among the 2 groups. In contrast to the nSSR results, the AMOVA based on the cpSSR dataset showed that a larger proportion of genetic variation could be attributed to variation within populations (68.84%) and among populations within groups (29.37%), and little variation among groups (1.79%). The overall F ST values calculated by AMOVA were 0.362 (P ≤ 0.0001) for the nSSR dataset and 0.312 (P ≤ 0.0001) for the cpSSR dataset.

Isolation by distance and pollen/seed migration ratios
The estimates of genetic differentiation (F ST value) based on 13 nSSRs ranged from 0.023 (between WJS and GJS) to 0.427 (between LDZ and FHS) (Table S4), excluding 2 cultivated and cpSSR were similar (Table 3), and the r was −1.149, indicating that most gene flow among populations occurs via seed, rather than via pollen.

Population bottleneck effect
Several populations had a significant excess of heterozygosity expected at mutation-drift equilibrium (i.e., He > Heq) (Piry, Luikart & Cornuet, 1990) under the 3 models in the bottleneck analysis, which indicated a deviation from mutation drift equilibrium in wild L. glauca populations (

Genetic variation within populations
In this study, sampling covered a large portion of the natural distribution, and overall genetic diversity across wild L. glauca populations exhibited low levels based on both nSSR (mean A R = 1.74, Ho = 0.33, He = 0.38, F IS = −0.14) and cpSSR (mean Nb = 2.06, Hrs = 0.39) loci. Our estimates of genetic diversity in L. glauca were almost half those of long-lived perennials (Ho = 0.63, He = 0.68), out-crossing species (Ho = 0.63, He = 0.65), and plants with wide distributions (Ho = 0.57, He = 0.62) (Nybom, 2004), which were lower than Laurus nobilis in Lauraceae (A R = 3.22, He = 0.56) (Marzouki et al., 2009). There are several major factors influencing variation that can result, each by itself or in combination, in the low levels of genetic variation observed in wild L. glauca populations. Asexual reproduction through apomictic seeds can decrease genetic variation in a population especially in apomictic populations (Lo, Stefanovic & Dickinson, 2009). Similarly, effective population size could seem like a limitation, due to the population being established from a limited number of individuals and a small sampling quantity. For two cultivated populations, these values were much lower than in wild populations, perhaps because apomixis reduces genetic diversity, or because of the potential confounding effect of small population sample size. Furthermore, the estimated values of genetic diversity therein are also lower than nSSR-based values found in literature (A R = 2.61, He = 0.44, F IS = −0.37) (Zhu et al., 2016). Differences in genetic variation between this study and the previous one are likely to be explained by the number of sampling populations and individuals the study of Zhu et al. (2016) included 6 populations, a total of 96 individual plants, while this study included 20 wild populations, making a total of 290 individual plants.
On the other hand, a marked similarity in the molecular variance revealed by the 2 types of markers (overall F sc = 0.293 and F ST = 0.362 for nSSRs; 0.299 and 0.312 for cpSSRs) was observed, indicating general consistency between chloroplasts and nuclear DNA. More specifically, predominant apomixis in wild L. glauca could explain how F ST observed for nuclear markers is a little higher than that for chloroplast markers, and vice versa. In addition, apomictic reproduction of L. glauca could also affect the results. Therefore the use of clone-corrected data (removing the data of clones from the same parent) is necessary if clones are detected, because L. glauca trees could form clones via vegetative reproduction with stolons.
Nine populations had more heterozygotes than expected. STRUCTURE analysis showed that the 9 populations were grouped into one group and the rest into another, suggesting negative F IS is an important factor affecting group of population difference. Usually, positive F IS values within populations (JGS, LYS, and FJS) indicated inbreeding, and negative F IS of wild populations (KYS, WYS, ZJS, WJS, GJS, NHS, FHS, and ZJJ) suggested outbreeding. However, considering that clonality probably generates significant negative F IS in some wild plant populations with asexual reproduction when considering all individuals (Stoeckel et al., 2006), the observed negative F IS of wild populations (KYS, WYS, ZJS, WJS, GJS, NHS, FHS, and ZJJ), coupled with the result grouped by NJ tree and the cluster pattern of the NJ tree (blue cluster) (Fig. 3A), suggest that apomixis may have been common in these populations. On the other hand, global multilocus tests indicated that sexual reproduction did exist in these populations that were positive F IS values (0 < F IS < 1). There is mixed reproductive system existed in this species, thus explaining the F IS pattern.

Genetic differentiation among populations
Populations can cluster according to habitat type or geographic distance. However, for species with predominantly asexual populations, like Daktulosphaira vitifoliae (Vorwerk & Forneck, 2006), Crataegus douglasii (Lo, Stefanovic & Dickinson, 2009), and Taraxacum officinale (Majesky et al., 2012), no significant correlation exists between genetic distances and geographic distances. In this study, a correlation between genetic distances (as measured by [F ST /(1 − F ST )] values) and geographic distances (in kilometers) was not detected in L. glauca populations, regardless of marker type, suggesting that weak barriers to dispersal between populations and/or the common founders between neighbors, distant populations, and apomictic populations did not completely limit gene flow. Results indicate that sexual dispersion and apomixis co-occurred in the same natural population. However, further research is needed to investigate the extent to which apomixis limits gene flow, and the exact rate of sexual production and apomixis that occurred among and within populations.
According to Ennos' formula (1994), the migration ratio of pollen flow versus seed flow (r) in this study was negative (−1.149), suggesting that gene flow between populations is more likely to involve seeds. Usually, the value of r is positive in almost all angiosperm species (Ennos, 1994). The biased r value could be explained by the reproductive mode in different populations being the apomictic seeds of females, coupled with the result that no staminate flower were found at field observation sites for 5 consecutive years. Considering that a situation where a few individuals can self-pollinate is likely, pollen from even just a few staminodes of female flowers may be associated with the negative r value of pollen flow versus seed flow. However, because the formula is derived for a hermaphrodite species, it needs to be modified to account for the disproportionate maternal contribution from the females to the next generation if the exact r value of L. glauca can be expected to be observed.

Source of evolutionary potential of apomicts
Although 22 haplotypes were observed across 20 wild populations and 2 cultivated populations of L. glauca, H12 accounted for 62.03% of overall haplotypes and was detected in all populations (Table S3), except for FHS, LYS, SZY, and HZY. This haplotype existing in many populations separated by considerable geographical distances (e.g., greater than 1,720 km between KYS and GJS), coupled with the complex network (Fig. 1B), suggested three inferences to account for the observed result. First, these individuals from different populations had a relatively recent shared female founder. Second, there was an apomictic lineage for this species within some populations. Third, the genome of this species has genetic traits that lead to apomixis, which might be induced by some factors (e.g., biological stimulation, environmental influence, climate changing, etc.), and coexists with sexual reproduction in identical individual plants.
The first inference could be explained by a hypothesis that the H12 haplotype may be associated with the migratory patterns of some birds responsible for the dispersion of apomictic seeds over a long distance. However, given that H12 is the most frequent and connected haplotype, probably is ancestral, this inference is not reliable. Besides, despite conducting field observations for 5 consecutive years, few birds were found eating grown fruits of L. glauca, as well as few small mammals (e.g., Paguma larvata). The second inference means having an early maternal ancestor through apomictic reproduction for many individuals. However, all the living species are at the tips of the tree of life, and apomixis is a derived condition (APG, 2003;APG, 2016;Horandl, 2006;Thompson & Ritland, 2006;Lo, Stefanovic & Dickinson, 2009). Consequently, having ancestral asexual angiosperms is almost impossible, because asexuals fail to maintain sex and recombination in populations that are limited in size, therefore are unable to bring together high-fitness alleles that reside in different individuals. The third inference means that apomicts of L. glauca may be of very recent origin and have the ability to apomictically reproduce through mutations or losses of some sexual genes. A similar situation exists in some species, such as Taraxacum officinale that exhibit alternations of asexual and sexual histories of apomicts (Majesky et al., 2012), some hawthorns (Crataegus; Rosaceae) that have a population genetic structure of diploid sexual and polyploid apomicts (Lo, Stefanovic & Dickinson, 2009), and a marbled crayfish (Procambarus virginalis) that reproduces through parthenogenesis (Gerhard et al., 2003;Ewen, 2018). Therefore, when coupled with the results of filed surveys that show apomixis occurred in all sampled populations, this haplotype (H12) is probably a predominant genotype existing in populations that have the ability to apomixis.
Above all, the authors reject the former two inferences that most populations across such long geographical distances have a relatively recent founder or an ancient apomictic ancestor, and believe the last inference that apomixis caused by mutations or losses of related genes makes this maternal haplotype (H12) present in many individuals from different populations more likely.

Conjecture about histories of natural populations
According to the relationship Ne = 4(NmNf )/(Nm + Nf ), where Nm represents the number of males, and Nf represents the number of females (Beerli & Palczewski, 2010), the accurate values of Ne could not be calculated because there were no males in the samples collected. Even so, based on the rarity of male individuals in mainland China, dioecious reproduction reported in the past several decades (Tsui, Xia & Li, 1982;Wang, 1972), some specimens of branches of male individuals stored in the China National Herbarium (PE), and many males grown on Taiwan (Zhang, 2007), it is inferred that some natural populations of L. glauca recently experienced a severe bottleneck and male individuals experienced a decline, most likely resulting from anthropogenic causes. However, in order to test the above hypotheses and obtain more accurate results, further field investigations are necessary, including samples from male and female individuals in a population (especially males found on Taiwan), correcting for clone and apomictic reproduction, across the full range of L. glauca, including Japan, South Korea, and Taiwan.

Implications for conservation
Genetic diversity is recognized as an important population attribute for both conservation and evolutionary purposes (Cena et al., 2006). The purpose of the conservation of endangered and threatened species is to maintain their contributions to overall genetic diversity. People usually focus only on endangered species and provide protection, whereas some species that are reduced in genetic diversity also require protection. This study detected a lower level of genetic diversity in L. glauca than that of some other species in Lauraceae. The conclusion is that wild L. glauca populations have female-skewed sex ratios, which was consistent with our field survey and sampling for 5 consecutive years. The destruction of habitats in order to plant other commercial or medicinal crops and felling by local farmers may explain the low frequency of male L. glauca in mainland China, although the species is common and widely distributed in other regions. The authors propose finding male individuals and promoting sexual reproduction in order to maintain this species' overall genetic diversity.

CONCLUSION
This study has shown low levels of genetic diversity in L. glauca across nearly its the entire natural distribution in mainland China. A complex correlation between populations was revealed by haplotype networks. Genetic structure within and among populations was similar at the nuclear and chloroplast levels. Furthermore, some populations experienced a recent bottleneck, and gene flow between populations is more likely to involve seed. This implies that wild L. glauca in mainland China has highly skewed sex ratios with predominant females.