Fig. 4 Changes in distribution range from the LGM (minority prediction) to present day (a, d, g, j), predicted elevational range during the LGM and at the present day (b, e, h, k) and historic effective population size (c, f, i, l) of H. tibetana (a, b, c), H. digitata(d, e, f), H. gregoryi (g, h, i), H. platon (j, k, l). The white dots in the distribution range size map show the actual observation of each population. In the predicted elevational range plots (b,e,h,k), the grey line presents the background elevation range of the whole study area, and the dashed line indicates the mean elevation of the suitable area of species at a time scenario (purple: during the LGM; cyan: at present). The maps can be viewed online at the original spatial resolution at https://fred.igb-berlin.de/data/package/812. The y-axis of the PSMC plots for H. gregoryi (i) was re-scaled to filter the noise closer to the present time (~5,000 years before the present).