Figures (5)  Tables (2)
    • Figure 1. 

      Distribution of 24 chloroplast haplotypes of Pterocarya hupehensis. The lower right panel shows the minimum spanning network of 24 chlorotypes. Circle sizes are proportional to the number of samples per haplotype. The red dashed line on the map represents the Sino-Himalayan/Sino-Japanese forest boundary. The black dashed line represents the phylogeographic break along the Yangtze River in the eastern Sichuan Basin. The black dotted lines delineate three phylogroups with closely related chlorotypes.

    • Figure 2. 

      BEAST-derived chronogram of 24 Pterocarya hupehensis haplotypes based on six chloroplast DNA (cpDNA) fragments. The divergence time are shown above branches, with blue bars indicating the 95% highest posterior densities (HPDs). Posterior probabilities are labeled at each node.

    • Figure 3. 

      Genetic structure of Pterocarya hupehensis based on 2,889 SNPs dataset. (a) Geographic origins of 18 P. hupehensis populations and their color-coded grouping at close to K = 2. The red dashed line represents the Sino-Himalayan/Sino-Japanese forest boundary. (b) Principal component analysis (PCA), with pink and blue colors representing two clusters. (c) Histogram of the Admixture analysis for P. hupehensis with K = 2.

    • Figure 4. 

      Hybridization among populations of Pterocarya hupehensis. (a) Maximum likelihood (ML) tree inferred using Treemix; gene-flow events are depicted with arrows, and ten migration events were allowed. Migration arrows are colored according to their weight. Horizontal branch lengths are proportional to the amount of genetic drift that occurred on each branch. The scale bar shows 10× the average standard error of the entries in the sample covariance matrix. (b) Heatmaps of residual fit from the ML tree. Residuals white through blue indicate that the corresponding populations are more closely related to each other than on the ML tree, suggesting confounding events between these populations. (c) Demographic history of P. hupehensis inferred with Stairway Plot 2 using unfolded site frequency spectra. The 95% confidence interval for estimated effective population size is shown with gray lines.

    • Figure 5. 

      Results of ecological niche modeling of P. hupehensis in four time periods from past to future. (a) Average projection of the model to present climatic conditions. (b) Average projection of the model to the last glacial maximum (LGM: c. 21 kya BP (before present)). (c) Average projection of the model for the last interglacial (LIG: c. 120–140kya BP). (d) Average projection of the model to the year 2070 (2061–2080) under an intermediate climate warming scenario (RCP 4.5). Colors from blue to red represent the degree of habitat suitability for P. hupehensis survival, from unsuitable to suitable.

    • CodeSiteLongitude
      (°E)
      Latitude
      (°N)
      RAD-SNPscpDNA
      nHOHEπFishπ × 103Haplotype (number
      of individuals)
      HHXian, Shaanxi107.9333.8760.1990.1730.2030.00700H1 (10)
      HXCAnkang, Shaanxi109.4932.7230.2070.1650.2090.0040.420.005H2 (1), H3 (1), H4 (7)
      FLCAnkang, Shaanxi109.431.9550.1920.1430.174−0.03300H5 (9)
      DSPBaoji, Shaanxi107.4933.8440.1830.1440.177−0.01300H6 (8)
      QDZNanyang, Henan111.9733.5260.2100.1750.197−0.0270.390.004H7 (7), H8 (2)
      HKCNanyang, Henan112.0233.5760.1930.1840.2090.0330.50H7 (6), H9 (3)
      TSGNanyang, Henan111.7233.6360.1820.1710.1980.0310.330H7 (5), H10 (1)
      LJLNanyang, Henan111.7033.6340.1910.1640.2070.02700H7 (4)
      LYGTianshui, Gansu106.134.2350.1890.1520.1770.02100H1 (10)
      YPCLongnan, Gansu106.2333.6750.2060.1490.1700.06800H11 (7)
      HJGKangxian, Gansu105.5133.3950.1840.1410.1620.04100H12 (7)
      MYGYangba, Gansu105.7433.0360.1910.1700.1960.01000H13 (6)
      SNJShennongjia, Hubei110.9231.6570.2160.1990.2170.00400H14 (7)
      XJZShennongjia, Hubei110.5831.59110.2090.2090.2200.0330.470.067H15 (7), H16 (3)
      SHJZEnshi, Hubei109.830.16110.2060.1790.194−0.0260.760.067H17 (5), H18 (1), H19 (2),
      H20 (1), H21 (1)
      JSZNanchuan, Chongqing107.1429.0270.2130.2120.2230.029
      DFXBijie, Guizhou105.8827.33120.2240.2120.222−0.0010.360.004H22 (8), H23 (2)
      NYXBijie, Guizhou105.4726.7130.2270.1950.204−0.04700H24 (10)
      AVERAGE0.2010.1740.198−0.0060.190.009
      n, the number of samples used for RAD-seq analysis; HO, observed heterozygosity; HE, expected heterozygosity; π, nucleotide diversity; Fis, inbreeding coefficient; h, haplotype diversity.

      Table 1. 

      Sample codes, sample locations, and genetic diversity indices for 18 populations of Pterocarya hupehensis based on RAD-seq data and cpDNA data.

    • Source of variationcpDNARAD-SNPs
      d.f.Percentage variation (%)Fixation indicesd.f.Percentage variation (%)Fixation indices
      Among lineages175.05FCT = 0.75122.22FCT = 0.22
      Among populations within lineages1522.64FSC = 0.911611.47FSC = 0.15
      Within populations1242.31FST = 0.9822666.31FST = 0.34

      Table 2. 

      Analyses of molecular variance (AMOVA) based on cpDNA data and RAD-seq data for Pterocarya hupehensis from western and eastern lineages.