Links to Exploration step
Le document en format XML
<record><TEI><teiHeader><fileDesc><titleStmt><title xml:lang="en">Refugial isolation and range expansions drive the genetic structure of <italic>Oxyria sinensis</italic> (Polygonaceae) in the Himalaya-Hengduan Mountains</title><author><name sortKey="Meng, Lihua" sort="Meng, Lihua" uniqKey="Meng L" first="Lihua" last="Meng">Lihua Meng</name><affiliation><nlm:aff id="a1"><institution>School of Life Sciences, Engineering Research Center of Sustainable Development and Utilization of Biomass Energy Ministry of Education, Yunnan Normal University</institution>, Kunming 650092, Yunnan,<country>P. R. China</country></nlm:aff></affiliation></author><author><name sortKey="Chen, Gang" sort="Chen, Gang" uniqKey="Chen G" first="Gang" last="Chen">Gang Chen</name><affiliation><nlm:aff id="a2"><institution>Graduate School of Oceanography, University of Rhode Island</institution>, 02881, RI<country>USA</country></nlm:aff></affiliation></author><author><name sortKey="Li, Zhonghu" sort="Li, Zhonghu" uniqKey="Li Z" first="Zhonghu" last="Li">Zhonghu Li</name><affiliation><nlm:aff id="a3"><institution>Key Laboratory of Resource Biology and Biotechnology in Western China (Ministry of Education), College of Life Sciences, Northwest University</institution>, Xi’an 710069, Shanxi,<country>P. R.China</country></nlm:aff></affiliation></author><author><name sortKey="Yang, Yongping" sort="Yang, Yongping" uniqKey="Yang Y" first="Yongping" last="Yang">Yongping Yang</name><affiliation><nlm:aff id="a4"><institution>Kunming Institute of Botany, Chinese Academy of Sciences</institution>, Kunming 650204, Yunnan,<country>P. R. China</country></nlm:aff></affiliation></author><author><name sortKey="Wang, Zhengkun" sort="Wang, Zhengkun" uniqKey="Wang Z" first="Zhengkun" last="Wang">Zhengkun Wang</name><affiliation><nlm:aff id="a1"><institution>School of Life Sciences, Engineering Research Center of Sustainable Development and Utilization of Biomass Energy Ministry of Education, Yunnan Normal University</institution>, Kunming 650092, Yunnan,<country>P. R. China</country></nlm:aff></affiliation></author><author><name sortKey="Wang, Liuyang" sort="Wang, Liuyang" uniqKey="Wang L" first="Liuyang" last="Wang">Liuyang Wang</name><affiliation><nlm:aff id="a5"><institution>Department of Biology, Duke University</institution>, Durham, 27708, NC<country>USA</country></nlm:aff></affiliation></author></titleStmt><publicationStmt><idno type="wicri:source">PMC</idno><idno type="pmid">26013161</idno><idno type="pmc">4445061</idno><idno type="url">http://www.ncbi.nlm.nih.gov/pmc/articles/PMC4445061</idno><idno type="RBID">PMC:4445061</idno><idno type="doi">10.1038/srep10396</idno><date when="2015">2015</date><idno type="wicri:Area/Pmc/Corpus">000045</idno><idno type="wicri:explorRef" wicri:stream="Pmc" wicri:step="Corpus" wicri:corpus="PMC">000045</idno></publicationStmt><sourceDesc><biblStruct><analytic><title xml:lang="en" level="a" type="main">Refugial isolation and range expansions drive the genetic structure of <italic>Oxyria sinensis</italic> (Polygonaceae) in the Himalaya-Hengduan Mountains</title><author><name sortKey="Meng, Lihua" sort="Meng, Lihua" uniqKey="Meng L" first="Lihua" last="Meng">Lihua Meng</name><affiliation><nlm:aff id="a1"><institution>School of Life Sciences, Engineering Research Center of Sustainable Development and Utilization of Biomass Energy Ministry of Education, Yunnan Normal University</institution>, Kunming 650092, Yunnan,<country>P. R. China</country></nlm:aff></affiliation></author><author><name sortKey="Chen, Gang" sort="Chen, Gang" uniqKey="Chen G" first="Gang" last="Chen">Gang Chen</name><affiliation><nlm:aff id="a2"><institution>Graduate School of Oceanography, University of Rhode Island</institution>, 02881, RI<country>USA</country></nlm:aff></affiliation></author><author><name sortKey="Li, Zhonghu" sort="Li, Zhonghu" uniqKey="Li Z" first="Zhonghu" last="Li">Zhonghu Li</name><affiliation><nlm:aff id="a3"><institution>Key Laboratory of Resource Biology and Biotechnology in Western China (Ministry of Education), College of Life Sciences, Northwest University</institution>, Xi’an 710069, Shanxi,<country>P. R.China</country></nlm:aff></affiliation></author><author><name sortKey="Yang, Yongping" sort="Yang, Yongping" uniqKey="Yang Y" first="Yongping" last="Yang">Yongping Yang</name><affiliation><nlm:aff id="a4"><institution>Kunming Institute of Botany, Chinese Academy of Sciences</institution>, Kunming 650204, Yunnan,<country>P. R. China</country></nlm:aff></affiliation></author><author><name sortKey="Wang, Zhengkun" sort="Wang, Zhengkun" uniqKey="Wang Z" first="Zhengkun" last="Wang">Zhengkun Wang</name><affiliation><nlm:aff id="a1"><institution>School of Life Sciences, Engineering Research Center of Sustainable Development and Utilization of Biomass Energy Ministry of Education, Yunnan Normal University</institution>, Kunming 650092, Yunnan,<country>P. R. China</country></nlm:aff></affiliation></author><author><name sortKey="Wang, Liuyang" sort="Wang, Liuyang" uniqKey="Wang L" first="Liuyang" last="Wang">Liuyang Wang</name><affiliation><nlm:aff id="a5"><institution>Department of Biology, Duke University</institution>, Durham, 27708, NC<country>USA</country></nlm:aff></affiliation></author></analytic><series><title level="j">Scientific Reports</title><idno type="eISSN">2045-2322</idno><imprint><date when="2015">2015</date></imprint></series></biblStruct></sourceDesc></fileDesc><profileDesc><textClass></textClass></profileDesc></teiHeader><front><div type="abstract" xml:lang="en"><p>The formation of the Mekong-Salween Divide and climatic oscillations in Pleistocene were the main drivers for the contemporary diversity and genetic structure of plants in the Himalaya-Hengduan Mountains (HHM). To identify the relative roles of the two historical events in shaping population history of plants in HHM, we investigated the phylogeographic pattern of <italic>Oxyria sinensis</italic>, a perennial plant endemic to the HHM. Sixteen chloroplast haplotypes were identified and were clustered into three phylogenetic clades. The age of the major clades was estimated to be in the Pleistocene, falling into several Pleistocene glacial stages and postdating the formation of the Mekong-Salween Divide. Range expansions occurred at least twice in the early and middle Pleistocene, but the spatial genetic distribution rarely changed since the Last Glacial Maximum. Our results suggest that temporary mountain glaciers may act as barriers in promoting the lineage divergence in <italic>O. sinensis</italic> and that subsequential range expansions and secondary contacts might reshape the genetic distribution in geography and blur the boundary of population differentiation created in the earlier glacial stages. This study demonstrates that Pleistocene climatic change and mountain glaciers, rather than the Mekong-Salween Divide, play the primary role in shaping the spatial genetic structure of <italic>O. sinensis</italic>.</p></div></front><back><div1 type="bibliography"><listBibl><biblStruct><analytic><author><name sortKey="Hewitt, G" uniqKey="Hewitt G">G. Hewitt</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hewitt, G M" uniqKey="Hewitt G">G. M. Hewitt</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Avise, J C" uniqKey="Avise J">J. C. Avise</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Myers, N" uniqKey="Myers N">N. Myers</name></author><author><name sortKey="Mittermeier, R A" uniqKey="Mittermeier R">R. A. Mittermeier</name></author><author><name sortKey="Mittermeier, C G" uniqKey="Mittermeier C">C. G. Mittermeier</name></author><author><name sortKey="Da Fonseca, G A" uniqKey="Da Fonseca G">G. A. da Fonseca</name></author><author><name sortKey="Kent, J" uniqKey="Kent J">J. Kent</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Li, X W" uniqKey="Li X">X. W. Li</name></author><author><name sortKey="Li, J" uniqKey="Li J">J. Li</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Mittermeier, R" uniqKey="Mittermeier R">R. Mittermeier</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wilson, E O" uniqKey="Wilson E">E. O. Wilson</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wu, Z" uniqKey="Wu Z">Z. Wu</name></author><author><name sortKey="Boufford, D E" uniqKey="Boufford D">D. E. Boufford</name></author><author><name sortKey="Ohba, H" uniqKey="Ohba H">H. Ohba</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Gao, L M" uniqKey="Gao L">L. M. Gao</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Li, Y" uniqKey="Li Y">Y. Li</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Qiu, Y X" uniqKey="Qiu Y">Y. X. Qiu</name></author><author><name sortKey="Fu, C X" uniqKey="Fu C">C. X. Fu</name></author><author><name sortKey="Comes, H P" uniqKey="Comes H">H. P. Comes</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Liu, Y" uniqKey="Liu Y">Y. Liu</name></author><author><name sortKey="Yang, S" uniqKey="Yang S">S. Yang</name></author><author><name sortKey="Ji, P" uniqKey="Ji P">P. Ji</name></author><author><name sortKey="Gao, L" uniqKey="Gao L">L. Gao</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhang, T C" uniqKey="Zhang T">T. C. Zhang</name></author><author><name sortKey="Comes, H P" uniqKey="Comes H">H. P. Comes</name></author><author><name sortKey="Sun, H" uniqKey="Sun H">H. Sun</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shi, Y F" uniqKey="Shi Y">Y. F. Shi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Clark, M" uniqKey="Clark M">M. Clark</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Liu, J" uniqKey="Liu J">J. Liu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ward, F K" uniqKey="Ward F">F. K. Ward</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhang, T" uniqKey="Zhang T">T. Zhang</name></author><author><name sortKey="Sun, H" uniqKey="Sun H">H. Sun</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Li, X W" uniqKey="Li X">X. W. Li</name></author><author><name sortKey="Li, J" uniqKey="Li J">J. Li</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Li, X W" uniqKey="Li X">X.W. Li</name></author><author><name sortKey="Li, J" uniqKey="Li J">J. Li</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhang, Z Y" uniqKey="Zhang Z">Z.Y. Zhang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Qiu, Y X" uniqKey="Qiu Y">Y. X. Qiu</name></author><author><name sortKey="Guan, B C" uniqKey="Guan B">B. C. Guan</name></author><author><name sortKey="Fu, C X" uniqKey="Fu C">C. X. Fu</name></author><author><name sortKey="Comes, H P" uniqKey="Comes H">H. P. Comes</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Liu, J Q" uniqKey="Liu J">J. Q. Liu</name></author><author><name sortKey="Gao, T G" uniqKey="Gao T">T. G. Gao</name></author><author><name sortKey="Chen, Z D" uniqKey="Chen Z">Z. D. Chen</name></author><author><name sortKey="Lu, A M" uniqKey="Lu A">A. M. Lu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Liu, J Q" uniqKey="Liu J">J. Q. Liu</name></author><author><name sortKey="Wang, Y J" uniqKey="Wang Y">Y. J. Wang</name></author><author><name sortKey="Wang, A L" uniqKey="Wang A">A. L. Wang</name></author><author><name sortKey="Hideaki, O" uniqKey="Hideaki O">O. Hideaki</name></author><author><name sortKey="Abbott, R J" uniqKey="Abbott R">R. J. Abbott</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wang, L" uniqKey="Wang L">L. Wang</name></author><author><name sortKey="Schneider, H" uniqKey="Schneider H">H. Schneider</name></author><author><name sortKey="Zhang, X C" uniqKey="Zhang X">X.-C. Zhang</name></author><author><name sortKey="Xiang, Q P" uniqKey="Xiang Q">Q.-P. Xiang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wang, L" uniqKey="Wang L">L. Wang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sun, Y" uniqKey="Sun Y">Y. Sun</name></author><author><name sortKey="Wang, A" uniqKey="Wang A">A. Wang</name></author><author><name sortKey="Wan, D" uniqKey="Wan D">D. Wan</name></author><author><name sortKey="Wang, Q" uniqKey="Wang Q">Q. Wang</name></author><author><name sortKey="Liu, J" uniqKey="Liu J">J. Liu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Chen, S" uniqKey="Chen S">S. Chen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shi, Y" uniqKey="Shi Y">Y. Shi</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Shi, Y F" uniqKey="Shi Y">Y. F. Shi</name></author><author><name sortKey="Zheng, B X" uniqKey="Zheng B">B. X. Zheng</name></author><author><name sortKey="Li, S J" uniqKey="Li S">S. J. Li</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wang, L Y" uniqKey="Wang L">L. Y. Wang</name></author><author><name sortKey="Ikeda, H" uniqKey="Ikeda H">H. Ikeda</name></author><author><name sortKey="Liu, T L" uniqKey="Liu T">T. L. Liu</name></author><author><name sortKey="Wang, Y J" uniqKey="Wang Y">Y. J. Wang</name></author><author><name sortKey="Liu, J Q" uniqKey="Liu J">J. Q. Liu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Opgenoorth, L" uniqKey="Opgenoorth L">L. Opgenoorth</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jia, D R" uniqKey="Jia D">D. R. Jia</name></author><author><name sortKey="Liu, T L" uniqKey="Liu T">T. L. Liu</name></author><author><name sortKey="Wang, L Y" uniqKey="Wang L">L. Y. Wang</name></author><author><name sortKey="Zhou, D W" uniqKey="Zhou D">D. W. Zhou</name></author><author><name sortKey="Liu, J Q" uniqKey="Liu J">J. Q. Liu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhang, Q" uniqKey="Zhang Q">Q. Zhang</name></author><author><name sortKey="Chiang, T Y" uniqKey="Chiang T">T. Y. Chiang</name></author><author><name sortKey="George, M" uniqKey="George M">M. George</name></author><author><name sortKey="Liu, J Q" uniqKey="Liu J">J. Q. Liu</name></author><author><name sortKey="Abbott, R J" uniqKey="Abbott R">R. J. Abbott</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Meng, L H" uniqKey="Meng L">L. H. Meng</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Yang, F S" uniqKey="Yang F">F. S. Yang</name></author><author><name sortKey="Li, Y F" uniqKey="Li Y">Y. F. Li</name></author><author><name sortKey="Ding, X" uniqKey="Ding X">X. Ding</name></author><author><name sortKey="Wang, X Q" uniqKey="Wang X">X. Q. Wang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Fan, D M" uniqKey="Fan D">D. M. Fan</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Cun, Y Z" uniqKey="Cun Y">Y. Z. Cun</name></author><author><name sortKey="Wang, X Q" uniqKey="Wang X">X. Q. Wang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Avise, J C" uniqKey="Avise J">J. C. Avise</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Liu, J Q" uniqKey="Liu J">J. Q. Liu</name></author><author><name sortKey="Sun, Y S" uniqKey="Sun Y">Y. S. Sun</name></author><author><name sortKey="Ge, X J" uniqKey="Ge X">X. J. Ge</name></author><author><name sortKey="Gao, L M" uniqKey="Gao L">L. M. Gao</name></author><author><name sortKey="Qiu, Y X" uniqKey="Qiu Y">Y. X. Qiu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Kuhle, M" uniqKey="Kuhle M">M. Kuhle</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Ehlers, J" uniqKey="Ehlers J">J. Ehlers</name></author><author><name sortKey="Gibbard, P L" uniqKey="Gibbard P">P. L. Gibbard</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Andersen, B G" uniqKey="Andersen B">B. G. Andersen</name></author><author><name sortKey="Borns, H W" uniqKey="Borns H">H. W. Borns</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhou, S" uniqKey="Zhou S">S. Zhou</name></author><author><name sortKey="Li, J" uniqKey="Li J">J. Li</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zheng, B" uniqKey="Zheng B">B. Zheng</name></author><author><name sortKey="Xu, Q" uniqKey="Xu Q">Q. Xu</name></author><author><name sortKey="Shen, Y" uniqKey="Shen Y">Y. Shen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Owen, L A" uniqKey="Owen L">L. A. Owen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wu, Y" uniqKey="Wu Y">Y. Wu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Li, B Y" uniqKey="Li B">B. Y. Li</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Owen, L A" uniqKey="Owen L">L. A. Owen</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhuo, Z" uniqKey="Zhuo Z">Z. Zhuo</name></author><author><name sortKey="Baoyin, Y" uniqKey="Baoyin Y">Y. Baoyin</name></author><author><name sortKey="Petit Maire, N" uniqKey="Petit Maire N">N. Petit-Maire</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhang, D F" uniqKey="Zhang D">D. F. Zhang</name></author><author><name sortKey="Fengquan, L" uniqKey="Fengquan L">L. Fengquan</name></author><author><name sortKey="Jianmin, B" uniqKey="Jianmin B">B. Jianmin</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Jia, D R" uniqKey="Jia D">D. R. Jia</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Zhao, F" uniqKey="Zhao F">F. Zhao</name></author><author><name sortKey="Yang, Y P" uniqKey="Yang Y">Y. P. Yang</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Thompson, J D" uniqKey="Thompson J">J. D. Thompson</name></author><author><name sortKey="Gibson, T J" uniqKey="Gibson T">T. J. Gibson</name></author><author><name sortKey="Plewniak, F" uniqKey="Plewniak F">F. Plewniak</name></author><author><name sortKey="Jeanmougin, F" uniqKey="Jeanmougin F">F. Jeanmougin</name></author><author><name sortKey="Higgins, D G" uniqKey="Higgins D">D. G. Higgins</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Sanchez, A" uniqKey="Sanchez A">A. Sanchez</name></author><author><name sortKey="Schuster, T M" uniqKey="Schuster T">T. M. Schuster</name></author><author><name sortKey="Kron, K A" uniqKey="Kron K">K. A. Kron</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Drummond, A J" uniqKey="Drummond A">A. J. Drummond</name></author><author><name sortKey="Suchard, M A" uniqKey="Suchard M">M. A. Suchard</name></author><author><name sortKey="Xie, D" uniqKey="Xie D">D. Xie</name></author><author><name sortKey="Rambaut, A" uniqKey="Rambaut A">A. Rambaut</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Muse, S V" uniqKey="Muse S">S. V. Muse</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Richardson, J E" uniqKey="Richardson J">J. E. Richardson</name></author><author><name sortKey="Pennington, R T" uniqKey="Pennington R">R. T. Pennington</name></author><author><name sortKey="Pennington, T D" uniqKey="Pennington T">T. D. Pennington</name></author><author><name sortKey="Hollingsworth, P M" uniqKey="Hollingsworth P">P.M. Hollingsworth</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Wolfe, K H" uniqKey="Wolfe K">K. H. Wolfe</name></author><author><name sortKey="Li, W H" uniqKey="Li W">W. H. Li</name></author><author><name sortKey="Sharp, P M" uniqKey="Sharp P">P. M. Sharp</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Guindon, S" uniqKey="Guindon S">S. Guindon</name></author><author><name sortKey="Gascuel, O" uniqKey="Gascuel O">O. Gascuel</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Darriba, D" uniqKey="Darriba D">D. Darriba</name></author><author><name sortKey="Taboada, G L" uniqKey="Taboada G">G.L. Taboada</name></author><author><name sortKey="Doallo, R" uniqKey="Doallo R">R. Doallo</name></author><author><name sortKey="Posada, D" uniqKey="Posada D">D. Posada</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Rambaut, A" uniqKey="Rambaut A">A. Rambaut</name></author><author><name sortKey="Drummond, A J" uniqKey="Drummond A">A. J. Drummond</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Clement, M" uniqKey="Clement M">M. Clement</name></author><author><name sortKey="Posada, D" uniqKey="Posada D">D. Posada</name></author><author><name sortKey="Crandall, K A" uniqKey="Crandall K">K. A. Crandall</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Templeton, A R" uniqKey="Templeton A">A. R. Templeton</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Pons, O" uniqKey="Pons O">O. Pons</name></author><author><name sortKey="Petit, R J" uniqKey="Petit R">R. J. Petit</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Excoffier, L" uniqKey="Excoffier L">L. Excoffier</name></author><author><name sortKey="Smouse, P E" uniqKey="Smouse P">P. E. Smouse</name></author><author><name sortKey="Quattro, J M" uniqKey="Quattro J">J. M. Quattro</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Excoffier, L" uniqKey="Excoffier L">L. Excoffier</name></author><author><name sortKey="Laval, G" uniqKey="Laval G">G. Laval</name></author><author><name sortKey="Schneider, S" uniqKey="Schneider S">S. Schneider</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Drummond, A J" uniqKey="Drummond A">A. J. Drummond</name></author><author><name sortKey="Rambaut, A" uniqKey="Rambaut A">A. Rambaut</name></author><author><name sortKey="Shapiro, B" uniqKey="Shapiro B">B. Shapiro</name></author><author><name sortKey="Pybus, O G" uniqKey="Pybus O">O. G. Pybus</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Harpending, H C" uniqKey="Harpending H">H. C. Harpending</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Rogers, A R" uniqKey="Rogers A">A. R. Rogers</name></author><author><name sortKey="Harpending, H" uniqKey="Harpending H">H. Harpending</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Tajima, F" uniqKey="Tajima F">F. Tajima</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Fu, Y X" uniqKey="Fu Y">Y. X. Fu</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Manni, F" uniqKey="Manni F">F. Manni</name></author><author><name sortKey="Guerard, E" uniqKey="Guerard E">E. Guerard</name></author><author><name sortKey="Heyer, E" uniqKey="Heyer E">E. Heyer</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Collins, W D" uniqKey="Collins W">W. D. Collins</name></author><author><name sortKey="Lee Taylor, J M" uniqKey="Lee Taylor J">J. M. Lee-Taylor</name></author><author><name sortKey="Edwards, D P" uniqKey="Edwards D">D. P. Edwards</name></author><author><name sortKey="Francis, G L" uniqKey="Francis G">G. L. Francis</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Collins, W D" uniqKey="Collins W">W. D. Collins</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Hasumi, H" uniqKey="Hasumi H">H. Hasumi</name></author><author><name sortKey="Emori, S" uniqKey="Emori S">S. Emori</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Phillips, S J" uniqKey="Phillips S">S. J. Phillips</name></author><author><name sortKey="Dudik, M" uniqKey="Dudik M">M. Dudík</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Araujo, M B" uniqKey="Araujo M">M. B. Araujo</name></author><author><name sortKey="New, M" uniqKey="New M">M. New</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Marmion, M" uniqKey="Marmion M">M. Marmion</name></author><author><name sortKey="Parviainen, M" uniqKey="Parviainen M">M. Parviainen</name></author><author><name sortKey="Luoto, M" uniqKey="Luoto M">M. Luoto</name></author><author><name sortKey="Heikkinen, R K" uniqKey="Heikkinen R">R.K. Heikkinen</name></author><author><name sortKey="Thuiller, W" uniqKey="Thuiller W">W. Thuiller</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Schoener, T W" uniqKey="Schoener T">T. W. Schoener</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Warren, D L" uniqKey="Warren D">D. L. Warren</name></author><author><name sortKey="Glor, R E" uniqKey="Glor R">R. E. Glor</name></author><author><name sortKey="Turelli, M" uniqKey="Turelli M">M. Turelli</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Warren, D L" uniqKey="Warren D">D. L. Warren</name></author><author><name sortKey="Glor, R E" uniqKey="Glor R">R. E. Glor</name></author><author><name sortKey="Turelli, M" uniqKey="Turelli M">M. Turelli</name></author></analytic></biblStruct><biblStruct><analytic><author><name sortKey="Warren, D L" uniqKey="Warren D">D. L. Warren</name></author><author><name sortKey="Seifert, S N" uniqKey="Seifert S">S. N. Seifert</name></author></analytic></biblStruct></listBibl></div1></back></TEI><pmc article-type="research-article"><pmc-dir>properties open_access</pmc-dir>
<front><journal-meta><journal-id journal-id-type="nlm-ta">Sci Rep</journal-id><journal-id journal-id-type="iso-abbrev">Sci Rep</journal-id><journal-title-group><journal-title>Scientific Reports</journal-title></journal-title-group><issn pub-type="epub">2045-2322</issn><publisher><publisher-name>Nature Publishing Group</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="pmid">26013161</article-id><article-id pub-id-type="pmc">4445061</article-id><article-id pub-id-type="pii">srep10396</article-id><article-id pub-id-type="doi">10.1038/srep10396</article-id><article-categories><subj-group subj-group-type="heading"><subject>Article</subject></subj-group></article-categories><title-group><article-title>Refugial isolation and range expansions drive the genetic structure of <italic>Oxyria sinensis</italic> (Polygonaceae) in the Himalaya-Hengduan Mountains</article-title></title-group><contrib-group><contrib contrib-type="author"><name><surname>Meng</surname><given-names>Lihua</given-names></name><xref ref-type="corresp" rid="c1">a</xref><xref ref-type="aff" rid="a1">1</xref></contrib><contrib contrib-type="author"><name><surname>Chen</surname><given-names>Gang</given-names></name><xref ref-type="aff" rid="a2">2</xref></contrib><contrib contrib-type="author"><name><surname>Li</surname><given-names>Zhonghu</given-names></name><xref ref-type="aff" rid="a3">3</xref></contrib><contrib contrib-type="author"><name><surname>Yang</surname><given-names>Yongping</given-names></name><xref ref-type="aff" rid="a4">4</xref></contrib><contrib contrib-type="author"><name><surname>Wang</surname><given-names>Zhengkun</given-names></name><xref ref-type="aff" rid="a1">1</xref></contrib><contrib contrib-type="author"><name><surname>Wang</surname><given-names>Liuyang</given-names></name><xref ref-type="corresp" rid="c2">b</xref><xref ref-type="aff" rid="a5">5</xref></contrib><aff id="a1"><label>1</label><institution>School of Life Sciences, Engineering Research Center of Sustainable Development and Utilization of Biomass Energy Ministry of Education, Yunnan Normal University</institution>, Kunming 650092, Yunnan,<country>P. R. China</country></aff><aff id="a2"><label>2</label><institution>Graduate School of Oceanography, University of Rhode Island</institution>, 02881, RI<country>USA</country></aff><aff id="a3"><label>3</label><institution>Key Laboratory of Resource Biology and Biotechnology in Western China (Ministry of Education), College of Life Sciences, Northwest University</institution>, Xi’an 710069, Shanxi,<country>P. R.China</country></aff><aff id="a4"><label>4</label><institution>Kunming Institute of Botany, Chinese Academy of Sciences</institution>, Kunming 650204, Yunnan,<country>P. R. China</country></aff><aff id="a5"><label>5</label><institution>Department of Biology, Duke University</institution>, Durham, 27708, NC<country>USA</country></aff></contrib-group><author-notes><corresp id="c1"><label>a</label><email>menglihua@mails.gucas.ac.cn</email></corresp><corresp id="c2"><label>b</label><email>lw157@duke.edu</email></corresp></author-notes><pub-date pub-type="epub"><day>27</day><month>05</month><year>2015</year></pub-date><pub-date pub-type="collection"><year>2015</year></pub-date><volume>5</volume><elocation-id>10396</elocation-id><history><date date-type="received"><day>27</day><month>06</month><year>2014</year></date><date date-type="accepted"><day>10</day><month>04</month><year>2015</year></date></history><permissions><copyright-statement>Copyright © 2015, Macmillan Publishers Limited</copyright-statement><copyright-year>2015</copyright-year><copyright-holder>Macmillan Publishers Limited</copyright-holder><license license-type="open-access" xlink:href="http://creativecommons.org/licenses/by/4.0/"><pmc-comment>author-paid</pmc-comment>
<license-p>This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit <ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/4.0/">http://creativecommons.org/licenses/by/4.0/</ext-link></license-p></license></permissions><abstract><p>The formation of the Mekong-Salween Divide and climatic oscillations in Pleistocene were the main drivers for the contemporary diversity and genetic structure of plants in the Himalaya-Hengduan Mountains (HHM). To identify the relative roles of the two historical events in shaping population history of plants in HHM, we investigated the phylogeographic pattern of <italic>Oxyria sinensis</italic>, a perennial plant endemic to the HHM. Sixteen chloroplast haplotypes were identified and were clustered into three phylogenetic clades. The age of the major clades was estimated to be in the Pleistocene, falling into several Pleistocene glacial stages and postdating the formation of the Mekong-Salween Divide. Range expansions occurred at least twice in the early and middle Pleistocene, but the spatial genetic distribution rarely changed since the Last Glacial Maximum. Our results suggest that temporary mountain glaciers may act as barriers in promoting the lineage divergence in <italic>O. sinensis</italic> and that subsequential range expansions and secondary contacts might reshape the genetic distribution in geography and blur the boundary of population differentiation created in the earlier glacial stages. This study demonstrates that Pleistocene climatic change and mountain glaciers, rather than the Mekong-Salween Divide, play the primary role in shaping the spatial genetic structure of <italic>O. sinensis</italic>.</p></abstract></article-meta></front><body><p>Inferring the evolutionary causes that drive genetic diversity of a population over time and space has been one of the central topics in phylogeography<xref ref-type="bibr" rid="b1">1</xref><xref ref-type="bibr" rid="b2">2</xref><xref ref-type="bibr" rid="b3">3</xref>. Over the past two decades, DNA-based phylogeographic studies have greatly advanced our understanding of demographic dynamics and the evolutionary history of plants in response to historical events at a much finer time and spatial scale. These events include geological changes and climate changes, such as the formation of vicariance, the uplift of mountains, and climatic oscillations associated with Pleistocene glaciers. However, how geological and climate changes affect demographic processes, genetic differentiation, and even speciation remains poorly understood. It is of particular interest to study global biodiversity hotspot regions<xref ref-type="bibr" rid="b4">4</xref>, which feature exceptional concentrations of endemic species and habitat changes.</p><p>The Himalaya-Hengduan Mountains (HHM) have been recognized as one of the worldwide global biodiversity “hotspots”<xref ref-type="bibr" rid="b4">4</xref>, containing more than 12,000 plant species, more than 20% of which are endemic<xref ref-type="bibr" rid="b4">4</xref><xref ref-type="bibr" rid="b5">5</xref><xref ref-type="bibr" rid="b6">6</xref><xref ref-type="bibr" rid="b7">7</xref><xref ref-type="bibr" rid="b8">8</xref>. Incredibly high inter-/intra-specific diversification rates of plants have been documented in this region and hypothetically attributed to the geological uplifts of Qinghai-Tibetan Plateau (QTP) and/or recurrent Pleistocene glacial-interglacial climate oscillation<xref ref-type="bibr" rid="b9">9</xref><xref ref-type="bibr" rid="b10">10</xref><xref ref-type="bibr" rid="b11">11</xref>. However, the relative roles of these two historical forces in generating genetic diversity in this region are inconsistent among previous studies, leading to two contrasting hypotheses: a geological hypothesis emphasizing a major role for geological impacts and a climatic hypothesis highlighting the role of Pleistocene climate changes<xref ref-type="bibr" rid="b10">10</xref><xref ref-type="bibr" rid="b12">12</xref><xref ref-type="bibr" rid="b13">13</xref>.</p><p>According to the geological hypothesis, a continuous population might have been broken into spatially isolated subpopulations when new mountains formed during the uplifts of the QTP approximately 40 million years ago (Mya). Those isolated subpopulations eventually evolved into genetically distinct lineages via the accumulation of genetic differences and via adaptation to the local environment. This hypothesis has been strongly supported by the geological, ecological and recent genetic evidence. The HHM serve as the southeast boundary of the QTP and include numerous northwest-southeast high mountains with deeply incised valleys and rivers resulting from the extensive uplift movement of the QTP starting 40 Mya<xref ref-type="bibr" rid="b14">14</xref>. This region has an altitude range of ca. 1,000 to 6,000 m a.s.l., creating abundant ecological heterogeneity and diverse niches for plants<xref ref-type="bibr" rid="b14">14</xref>. Moreover, the geological changes have been episodic rather than steady. For example, the latest uplifts of the QTP between the late Miocene and Pliocene have greatly modified the tectonic morphology of the HHM, thereby creating novel geological disjunctions and rearranging the river drainage system<xref ref-type="bibr" rid="b15">15</xref>. Patterns of phylogeographic disjunctions that are concordant with the geological barriers have been found for a number of different plants<xref ref-type="bibr" rid="b9">9</xref><xref ref-type="bibr" rid="b10">10</xref><xref ref-type="bibr" rid="b16">16</xref>. A classic case was the Mekong-Salween Divide, a geological barrier that arose ca. 4 Mya and has been recognized for nearly a century<xref ref-type="bibr" rid="b17">17</xref>. The Mekong-Salween Divide was believed to have driven the population divergence of two plant species, <italic>Sinopodophyllum hexandrum</italic><xref ref-type="bibr" rid="b10">10</xref> and <italic>Taxus wallichiana</italic><xref ref-type="bibr" rid="b9">9</xref>. In the former, the estimated divergence time between genealogies from the west and east side of Mekong-Salween Divide fell into the late Miocene, agreeing with the formation of the Mekong-Salween Divide. In another case, the present genealogical distribution of the plant species <italic>Terminalia franchetii</italic> was believed to be geographically structured by the paleo-drainage re-arrangements rather than the modern drainage systems, highlighting the role of historical geologic events in shaping contemporary genetic distribution<xref ref-type="bibr" rid="b13">13</xref><xref ref-type="bibr" rid="b18">18</xref>. A second known geological line is the “Tanaka-Kaiyong Line” (TKL), separating floristic subkingdoms into the Sino-Himalayan Forest (west) and Sino-Japanese Forest<xref ref-type="bibr" rid="b19">19</xref><xref ref-type="bibr" rid="b20">20</xref>. Molecular dating in several studies found the phylogeographic disjunction coincided well with the formation of this line (ca. 3 Mya), providing genetic evidence that geological processes might contribute to evolutionary history in plants<xref ref-type="bibr" rid="b9">9</xref><xref ref-type="bibr" rid="b21">21</xref><xref ref-type="bibr" rid="b22">22</xref>. Based on that evidence, the extensive uplifts of the QTP until ca. 4 Mya have been hypothesized to be the important force accelerating speciation<xref ref-type="bibr" rid="b23">23</xref><xref ref-type="bibr" rid="b24">24</xref><xref ref-type="bibr" rid="b25">25</xref> and generating high biodiversity at both the inter- and intra-specific levels in the HHM region<xref ref-type="bibr" rid="b9">9</xref><xref ref-type="bibr" rid="b10">10</xref><xref ref-type="bibr" rid="b26">26</xref><xref ref-type="bibr" rid="b27">27</xref><xref ref-type="bibr" rid="b28">28</xref>.</p><p>In contrast to the geological hypothesis, the climatic hypothesis attributes population genetic diversification and species diversity mainly to Pleistocene glacial cycling. Specifically, it is hypothesized that Pleistocene glaciers may have isolated populations into glacial refugia and promoted population differentiation and intra-specific diversification without apparent geological breaks<xref ref-type="bibr" rid="b1">1</xref>. Fossil and moraine records show that the QTP and the adjacent HHM were subjected to a series of glaciations in the Pleistocene<xref ref-type="bibr" rid="b29">29</xref><xref ref-type="bibr" rid="b30">30</xref>. The ice sheet advanced and receded repeatedly during glacial-interglacial cycles, creating temporal physical barriers and dynamic heterogeneous niches. This process is believed to have facilitated genetic divergence<xref ref-type="bibr" rid="b1">1</xref><xref ref-type="bibr" rid="b3">3</xref>. Under this hypothesis, two different scenarios were proposed in recent molecular phylogeographic studies. The first scenario is high-altitude adaptation, in which some plant species survived the glacial stages at high altitude areas <italic>in situ</italic><xref ref-type="bibr" rid="b26">26</xref><xref ref-type="bibr" rid="b31">31</xref><xref ref-type="bibr" rid="b32">32</xref><xref ref-type="bibr" rid="b33">33</xref> and developed accompanying deep allopatric lineages with the Pleistocene glaciers<xref ref-type="bibr" rid="b11">11</xref><xref ref-type="bibr" rid="b26">26</xref><xref ref-type="bibr" rid="b33">33</xref>, even without large unified ice-sheets existing in the late Pleistocene<xref ref-type="bibr" rid="b14">14</xref><xref ref-type="bibr" rid="b29">29</xref>. The second scenario is migratory colonization, in which some plants in refugia migrated to low-altitude ice-free areas (e.g., the east edge of the QTP and the Hengduan Mountains) and recolonized the inner QTP platform during interglacial or postglacial stages<xref ref-type="bibr" rid="b2">2</xref><xref ref-type="bibr" rid="b34">34</xref><xref ref-type="bibr" rid="b35">35</xref><xref ref-type="bibr" rid="b36">36</xref>. The Hengduan Mountain region were proposed as an important refugium to the QTP and neighboring areas<xref ref-type="bibr" rid="b10">10</xref><xref ref-type="bibr" rid="b36">36</xref>, and several case studies have suggested that mountain glaciers in the Pleistocene rather than geological breaks resulted in allopatric divergence and profoundly affected the intraspecific phylogeographic structure<xref ref-type="bibr" rid="b10">10</xref><xref ref-type="bibr" rid="b37">37</xref>.</p><p>The two contrasting hypotheses to explain the causes of the high biodiversity in the HHM region have each been supported by a number of studies. However, they are not necessarily incompatible with each other. The effects of both historical events (geological and Pleistocene climatic events) may be tracked through the genomic imprinting of extant species in phylogeographic studies. Interestingly, up to now, most phylogeographic studies have focused on plants with shallow lineages, and the evolutionary history of these lineages primarily reflects the effects of Pleistocene climate changes<xref ref-type="bibr" rid="b34">34</xref><xref ref-type="bibr" rid="b35">35</xref><xref ref-type="bibr" rid="b38">38</xref>. Plants originating before the Pleistocene could provide good opportunities to infer the evolutionary forces due to both geologic changes in the Miocene and subsequent Pleistocene climate changes that occurred after their speciation. In this study, we studied the phylogeography of an alpine plant, <italic>Oxyria sinensis</italic> Hemsley (Polygonaceae), which diverged from its sole sister species <italic>O. digyna</italic> in the genus <italic>Oxyria</italic> between 12 and 14 Mya<xref ref-type="bibr" rid="b27">27</xref>. Because the TKL line runs primarily from ca. 33 ° N/102 ° E to 19 ° N/108 ° E in the Yunnan and Sichuan provinces<xref ref-type="bibr" rid="b19">19</xref><xref ref-type="bibr" rid="b20">20</xref>, it has seldom overlapped with the distribution of <italic>O. sinensis</italic>. Thus, in this study, we specifically aimed to test how the Mekong-Salween Divide formation and/or Pleistocene glacial-interglacial fluctuations affected the population history of <italic>O. sinensis</italic>.</p><p>In contrast to the worldwide distribution of <italic>O. digyna</italic>, <italic>O. sinensis</italic> is endemic to southwest China, with a center in the Hengduan Mountains, westward extension to the southeast Himalaya, and eastward extension to the east edge of the QTP (east QTP). There is a clear break in the distribution at the Mekong-Salween Divide, where the Mekong-Salween Divide appears to segregate <italic>O. sinensis</italic> into west and east populations (<xref ref-type="fig" rid="f1">Fig. 1</xref>). Typical habitats of <italic>O. sinensis</italic> include mountain slopes, valleys and riversides, all at altitudes ranging from 1600 to 3800 m a.s.l. Considering the early origination and various ecological environments, <italic>O. sinensis</italic> must have experienced both of the above-mentioned historical events, making it an ideal species to test the evolutionary consequences of these events. In most angiosperms, nuclear DNA is biparentally transmitted and shows strong geographical homogenization because of extensive hybridization and/or introgression<xref ref-type="bibr" rid="b3">3</xref>. However, maternally inherited chloroplast DNA, which is dispersed through seeds without recombination, tends to be geographically structured. Chloroplast DNA has been demonstrated to be a powerful tool for tracing population history, founder effects, and demographic fluctuation<xref ref-type="bibr" rid="b3">3</xref><xref ref-type="bibr" rid="b39">39</xref>. In this study, we surveyed chloroplast (cp) DNA variations in <italic>O. sinensis</italic> using the <italic>matK</italic> fragment from the HHM region. We first estimated the genetic diversity and the divergence times of the different lineages detected. We next examined the population structure among different geographical groups. Then, we attempted to determine the possible factors driving the intraspecific divergence and the population structure. Finally, we used Bayesian inference methods and species distribution modeling to model the distribution changes of <italic>O. sinensis</italic> in response to Pleistocene climate changes. Thus, we had the following aims: 1) reveal the genetic diversity and geographic structure of contemporary <italic>O. sinensis</italic> populations; 2) reconstruct the demographic history of <italic>O. sinensis</italic>; 3) estimate the evolutionary consequences of the two historical events and assess their relative roles in shaping the contemporary genetic pattern of <italic>O. sinensis</italic>; and 4) infer the locations of the refugia in the Last Glacial Maximum (LGM) and regions of over-represented diversity that can help in the conservation and management of biodiversity in the HHM region.</p><sec disp-level="1" sec-type="results"><title>Results</title><sec disp-level="2"><title>Phylogenetic history of <italic>O. sinensis</italic> and molecular dating</title><p>We identified 16 distinct <italic>mat</italic>K haplotypes from 477 individuals in 38 populations (H1-H16; <xref ref-type="supplementary-material" rid="S1">Supplementary Tables S1 & S2</xref>; <xref ref-type="fig" rid="f1">Fig. 1</xref>). Among them, H3 was the most common, occurring in 22 of the 28 populations across the Himalaya, the Hengduan Mountains, and east QTP. Haplotype H8 was the second most common, occurring in ten populations, one in the Hengduan Mountains and nine in the east QTP. Six haplotypes (H5, H12–16) were unique to individual populations, and the remaining haplotypes occurred in more than one population. Fourteen populations (37.8%) had a single haplotype, and other populations (63.2%) harbored two or more haplotypes (<xref ref-type="supplementary-material" rid="S1">Supplementary Table S1</xref>; <xref ref-type="fig" rid="f1">Fig. 1</xref>).</p><p>The distinct <italic>mat</italic>K haplotypes differed at 18 variable sites within the 1173-bp sequence alignment. All variations were polymorphic nucleotide substitutions, and no insertions or deletions were found. The Bayesian phylogenetic analyses identified three well-supported clades (named A, B, and C). Clade A appeared to be a widely distributed clade, consisting of 8 unique haplotypes that were widespread throughout the region. Broadly, Clade B contained 5 haplotypes, including an unresolved haplotype H6 (support value < 50%). After excluding the H6, Clade B was exclusively distributed in the east QTP. The H6 mainly occurred in the Hengduan Mountain region. Clade C contained three haplotypes, all of which occurred only in the Himalaya. The minimum-spanning network of these haplotypes produced a similar grouping pattern (<xref ref-type="fig" rid="f2">Fig. 2</xref>).</p><p>The estimated divergence time between <italic>O. sinensis</italic> and <italic>O. digyna</italic> was approximately 13.79 Mya (95% HPD: 7.54–24.63 Mya) in the Miocene. This estimate is well consistent with an estimate based on the fossil record (12–14 Mya)<xref ref-type="bibr" rid="b27">27</xref><xref ref-type="bibr" rid="b40">40</xref>, suggesting that the substitution rates and models applied in this study were appropriate. Clades A and B diverged from C approximately 1.74 Mya (95% HPD: 0.72–3.46 Mya) in the early Pleistocene, and the recent divergence between clades A and B occurred approximately 0.86 Mya (95% HPD: 0.34–1.73 Mya).</p></sec><sec disp-level="2"><title>Genetic diversity</title><p>Genetic differentiation for the overall populations (<italic>G</italic><sub>ST</sub> = 0.611; <italic>N</italic><sub>ST</sub> = 0.662) was high. The high population differentiation was also supported by AMOVA and BARRIER analyses, as reflected by the distinct haplotype composition among the three geographic groups or the predicted phylogeographic groups. However, <italic>N</italic><sub>ST</sub> was not significantly higher than <italic>G</italic><sub>ST</sub> (<italic>P </italic>> 0.05), indicating a nonsignificant phylogeographic structure. This might be due to the discontinuous distribution of the shared dominant haplotypes (e.g., H3 and H8) across different geographic regions. Furthermore, BARRIER predicted that a major phylogeographic structure existed between the HHM and the east QTP (<xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S2</xref>). AMOVA analysis revealed that genetic differentiation across all populations accounted for 66.82% of the total variation, whereas among three geographic groups (Himalaya, Hengduan Mountains and east QTP), the variance components accounted for a slightly higher proportion of variation (38.95%) compared with the within-population among groups and the within-population variation (<xref ref-type="table" rid="t1">Table 1</xref>). Theoverall genetic diversity was very high (<italic>H</italic><sub>T</sub> = 0.791), contrasting the low degree of within-population diversity (<italic>H</italic><sub>S</sub> = 0.308). The haplotype diversity (<italic>h</italic><sub><italic>e</italic></sub>) within populations varied from 0.000 to 0.833 (<xref ref-type="table" rid="t1">Table 1</xref>; <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S1</xref>). We also computed the average haplotype diversity for each group. The highest diversity (0.6297) was found in the Himalaya, with an intermediate value in the Hengduan Mountains (0.3522), and the lowest value in east QTP (0.1005).</p></sec><sec disp-level="2"><title>Population demographic expansion</title><p>We investigated the population demographic history of <italic>O. sinensis</italic> for the three focal geographic regions: Himalaya, Hengduan Mountains and the east QTP. The Himalaya populations (1–3) showed a multimodal mismatch distribution (<xref ref-type="fig" rid="f3">Fig. 3</xref>), suggesting that no historical expansions occurred in this region. By contrast, unimodal mismatch distributions were found from the Hengduan Mountains and east QTP (<xref ref-type="fig" rid="f3">Fig. 3</xref>), which did not differ significantly from the expected distribution under a sudden population expansion model (P > 0.05, <xref ref-type="table" rid="t2">Table 2</xref>). In addition, theneutrality test for the Hengduan Mountains and east QTP groups showed negative Tajima’s <italic>D</italic> and Fu’s <italic>Fs</italic> values, though these were not statistically significant (except Fu’s <italic>Fs</italic> for east QTP, P < 0.001; <xref ref-type="table" rid="t2">Table 2</xref>), suggesting that the populations in these regions deviated from the neutral equilibrium population model. The expansion time estimated for the Hengduan Mountains and the east QTP based on mismatch analysis was approximately 0.14 and 0.31 Mya, respectively (<xref ref-type="table" rid="t2">Table 2</xref>). A Bayesian skyline plots (BSPs) analysis of total populations indicated a slight decline in population size since 0.8 Mya and an approximately five-fold rapid increase since 0.3 Mya (<xref ref-type="fig" rid="f3">Fig. 3</xref>), which is consistent with the mismatch analysis.</p></sec><sec disp-level="2"><title>Species distribution modeling in the present and the past</title><p>We reconstructed the glaciers under present conditions and the LGM. Because there is still no consensus about the exact locations and range of glaciers at the LGM in the QTP, we used the annual mean temperature layer from both ecological conditions. It is believed that the ice sheet might have reached 2000 meters in thickness and covered the entire plateau regions at an altitude above 3000 meters<xref ref-type="bibr" rid="b41">41</xref>. In general, our reconstructed glaciers are congruent with those of Ehlers and Gibbard<xref ref-type="bibr" rid="b42">42</xref>, who use minimum-glaciations range reconstruction for glacial geology. Most glaciers run along the mountains ridges in the HHM, indicating that mountain glaciers may have played an important role in contributing to the population structure of <italic>O. sinensis</italic>.</p><p>Species distribution models (SDM) were reconstructed for the present and the LGM (<xref ref-type="fig" rid="f4">Fig. 4b–d)</xref>. Model performance showed that all models met the threshold for inclusion in the ensemble, with all three statistical scores > 0.72, indicating good model performance. The realized present distribution of <italic>O. sinensis</italic> was generally well recovered by the consensus models (<xref ref-type="fig" rid="f4">Fig. 4b</xref>). Projection to the environmental niche of the LGM indicates suitable habitats essentially located within the mountain valley of the HHM, a low altitude area of south Hengduan Mountain and east QTP (<xref ref-type="fig" rid="f4">Fig. 4c,d</xref>). A high concordance between the suitable climate areas under the modeled LGM and the present was found (Schoener’s <italic>I</italic> = 0.9585, <italic>D</italic> = 0.7936 for CCSM; <xref ref-type="table" rid="t3">Table 3</xref>), suggesting a stable distribution of <italic>O. sinensis</italic> over the last 21,000 years. The SDM uncovered a dynamic history of geographic shifts in suitable areas for <italic>O. sinensis</italic> (<xref ref-type="fig" rid="f4">Fig. 4b–d</xref>), with mountain valleys and low-altitude regions maintaining suitable habitats, possibly serving as refugia in the glacial periods throughout the glacial-interglacial cycles (<xref ref-type="fig" rid="f5">Fig. 5</xref>).</p></sec></sec><sec disp-level="1" sec-type="discussion"><title>Discussion</title><p>Recently, several studies in plants found that phylogeographic disjunction coincided with the formation of the Mekong-Salween Divide<xref ref-type="bibr" rid="b9">9</xref><xref ref-type="bibr" rid="b10">10</xref><xref ref-type="bibr" rid="b16">16</xref>. Our data do not conform to this scenario but indicate that the main chloroplast genetic divergence of <italic>O. sinensis</italic> coincided with several Pleistocene glacial stages. Much of the variation within the Himalaya and east QTP is likely due to isolation and population expansions induced by glacial–interglacial cycles. In the following sections, we first discuss Pleistocene climate changes rather than the formation of Mekong-Salween Divide as the possible primary evolutionary factors driving intraspecific divergence. We then highlight the impact of Pleistocene glacial-interglacial cycles and range expansions on the maternal evolutionary history of <italic>O. sinensis</italic> and further discuss the implications of our findings for conservation biology.</p><sec disp-level="2"><title>Refugial isolation vs. geological isolation</title><p>In this study, we explored the phylogeography of <italic>O. sinensis.</italic> The estimated divergence time between <italic>O. sinensis</italic> and its closely related species <italic>O. digyna</italic> was approximately 13.8 Mya (<xref ref-type="fig" rid="f2">Fig. 2</xref>). Although the time calibration must be interpreted cautiously, our estimated speciation timing of <italic>O. sinensis</italic> was well consistent with dating based on fossils from a recent study<xref ref-type="bibr" rid="b27">27</xref>. Since <italic>O. sinensis</italic> diverged from its sibling species <italic>O. digyna</italic> approximately 13.8 Mya (<xref ref-type="fig" rid="f2">Fig. 2</xref>), it must have experienced both Pleistocene climate changes and the formation of the Mekong-Salween Divide. Both historical events presumably contributed to population isolation and thus to allopatric genetic differentiation between isolated populations. However, our phylogenetic analyses identified three major clades (<xref ref-type="fig" rid="f2">Fig. 2</xref>), and those clades show spatial clustering into three primary geographic regions. In contrast to early origination, relatively shallow intraspecific lineages were found in our study (<xref ref-type="fig" rid="f2">Fig. 2</xref>). The divergence time for the three major clades was approximately 1.74 to 0.86 Mya, which falls into the early or middle Pleistocene (<xref ref-type="fig" rid="f2">Fig. 2</xref>). The age of the deepest lineages of <italic>O. sinensis</italic> clearly post-dates the formation of a known geological break, the Mekong-Salween Divide, which arose in the late Miocene and early Pliocene<xref ref-type="bibr" rid="b9">9</xref><xref ref-type="bibr" rid="b10">10</xref>. In addition, BARRIER analysis found that a major phylogeographic disjunction occurred between the HHM and east QTP, not following the Mekong-Salween Divide. Taken together, these findings suggest that geological isolation might not be the main physical factor driving lineage divergence in <italic>O. sinensis</italic>.</p><p>By contrast, the estimated divergence times for intraspecific lineages are in concordant with two glacial events: the Eburonian glacial stage (1.2–1.7 Mya) and the Maximum glacial stage (0.6–0.8 Mya) in the early and/or middle Pleistocene<xref ref-type="bibr" rid="b43">43</xref>. On the QTP, the Pleistocene glaciers reached a maximum altitude during the maximum glacial stage<xref ref-type="bibr" rid="b30">30</xref><xref ref-type="bibr" rid="b44">44</xref><xref ref-type="bibr" rid="b45">45</xref>. Ice covered five to seven times more area than at present, and extensive mountain glaciers developed in the eastern Himalaya and Hengduan Mountains<xref ref-type="bibr" rid="b30">30</xref><xref ref-type="bibr" rid="b45">45</xref><xref ref-type="bibr" rid="b46">46</xref><xref ref-type="bibr" rid="b47">47</xref>. Therefore, we suppose that these mountains may harbor individuals of <italic>O. sinensis</italic> that survived <italic>in situ</italic> for a long time and that the temporary mountain glaciers might have acted as physical barriers of gene flow and promoted fragmentation and genetic differentiation. Given a sufficiently long time, population fragmentation and genetic differentiation eventually arose. Although a number of sequential glaciers occurred after the maximum glacial stage, they had minor impacts on the extant spatial genetic distribution of <italic>O. sinensis</italic><xref ref-type="bibr" rid="b14">14</xref>. Bearing this in mind, we propose that the development of mountain glaciers during early/middle Pleistocene glacial stages mostly likely served as the main factor causing lineage divergence of <italic>O. sinensis</italic> due to its relatively low intensity. In agreement with this study, Fan <italic>et. al</italic>.<xref ref-type="bibr" rid="b37">37</xref> also found that climate changes (monsoon) associated with the formation of TKL rather than TKL rising drove the spatial genetic structure of <italic>Sophora davidii</italic>. Deep lineages and similar diversification patterns due to Pleistocene climate changes were also found in other plants in the QTP region<xref ref-type="bibr" rid="b10">10</xref><xref ref-type="bibr" rid="b26">26</xref><xref ref-type="bibr" rid="b31">31</xref><xref ref-type="bibr" rid="b33">33</xref>.</p></sec><sec disp-level="2"><title>Multiple glacial refugia and demographic history</title><p>Glacial refugia were recognized as harboring high intraspecific diversity and major lineages<xref ref-type="bibr" rid="b1">1</xref>. The identification of three phylogenetic lineages suggests that at least three historical refugia existed in the followed recurrent cycles of Pleistocene cooling and warming since the origination of <italic>O. sinensis</italic>. The Himalaya harbored three haplotypes on the basal clade, representing the eldest refugium. A number of populations which harbor high genetic diversity and abundant unique haplotypes were potential refugia during the cooling stages (<xref ref-type="fig" rid="f1">Fig. 1</xref>; <xref ref-type="supplementary-material" rid="S1">Supplementary Fig. S1</xref>); for example, the Himalaya populations, the central Hengduan Mountains (e.g., populations 10, 13 and 28), the southern Hengduan Mountains (e.g., populations 21–23), and the eastern QTP (e.g., population 31, 34). Given the early and middle Pleistocene origination of these unique haplotypes (<xref ref-type="fig" rid="f2">Fig. 2</xref>), it is reasonable that these populations have served as glacial refugia in the cooling stages since their origination and expanded their range rapidly during the following interglacial stages. The SDM also predicted suitable habitats of <italic>O. sinensis</italic> sundered by high mountains resulting in multiple centers of genetic diversity, reflecting the geographic isolates (<xref ref-type="fig" rid="f4">Fig. 4b–d</xref>). Therefore, <italic>O. sinensis</italic> likely survived in those multiple locations <italic>in situ</italic> during the LGM. The SDM predicted that the most suitable habitats under present conditions were highly concordant but slightly larger than during the LGM (<xref ref-type="fig" rid="f4">Fig. 4b,c</xref>) and that the present mountain glaciers were highly similar to those in the LGM but slightly smaller at present (<xref ref-type="fig" rid="f5">Fig. 5</xref>). These results suggest that <italic>O. sinensis</italic> seldom changed its distribution since the LGM. This also indicates that mountain glaciers in this region have not greatly expanded but have slightly shrunk since the LGM. Overall, our findings suggest that the mountain glaciers and range expansions in early and middle Pleistocene had a predominant effect on the spatial genetic structure of <italic>O. sinensis</italic>. The fact that suitable areas have changed relatively little between the LGM and the present indicates that <italic>O. sinensis</italic> has only been slightly affected by ice-sheets and/or mountain glaciers since the LGM<xref ref-type="bibr" rid="b48">48</xref><xref ref-type="bibr" rid="b49">49</xref><xref ref-type="bibr" rid="b50">50</xref>. This finding supports the theory that no large ice sheet/mountain glacier developed in the Hengduan Mountains region during this period.</p><p>The duration of mountain glaciers plays an important role in lineage and population divergence. However, rare gene exchanges most likely occurred among neighboring regions during the retreat of the glaciers. The Himalaya group harbored three haplotypes, all of which located on basal clade C and may represent the earliest diversification center. The presence of haplotype H3 in the Himalaya region could be due to occasional gene flow from the Hengduan Mountains and might have been achieved through dispersion along the Salween River valley. Such long-distance dispersion might also occur between the Yunnan and Sichuan ranges, leading to unclear boundaries between the Hengduan Mountains and east QTP (<xref ref-type="fig" rid="f1">Fig. 1</xref>). Those long-distance migration patterns were also found in <italic>Sinopodophyllum hexandrum</italic><xref ref-type="bibr" rid="b10">10</xref>, an alpine herb distributed sympatrically with <italic>O. sinensis</italic>. An alternative scenario is that those shared haplotypes might be relics and/or close derivatives of ancient haplotypes due to an ancient common ancestry and incomplete lineage sorting. In the network topology, we observed a couple of intermediate states, which may represent those missing ancestral haplotypes and/or unsampled haplotypes.</p><p>The estimates of <italic>G</italic>st and <italic>N</italic>st (<italic>G</italic><sub>ST</sub> = 0.6110; <italic>N</italic><sub>ST</sub> = 0.6620) indicate high population differentiation across the entire distribution of <italic>O. sinensis</italic>, but no significant phylogeographic structure was found based on a comparison of <italic>G</italic><sub>ST</sub> and <italic>N</italic><sub>ST</sub> (P > 0.05). The high geographic genetic differentiation primarily reflects the composition of the distinct haplotypes among the three geographic groups or the phylogeographic disjunction between the HHM and east QTP identified from BARRIER. The nonsignificant phylogeographic signal most likely occurred because of the discontinuous distribution of the common haplotypes (e.g., H3 and H8) across different geographic regions. We detected the signatures of repeated regional expansions based on a combination of multiple statistical analyses, including Fu’s <italic>F</italic>s and Tajima’s <italic>D</italic> statistics, BSPs, and hierarchical mismatch distribution. BSPs indicated that the population size had increased dramatically since 0.5 Mya (<xref ref-type="fig" rid="f3">Fig. 3f</xref>). Further hierarchical mismatch analyses suggested at least two independent population expansions occurred, and the expansion times were dated to before LGM. The first expansion led by haplotype H8 occurred mainly in the east QTP approximately 0.3 Mya and fell into the interglacial period after the Maximum glaciations on the QTP<xref ref-type="bibr" rid="b29">29</xref><xref ref-type="bibr" rid="b46">46</xref><xref ref-type="bibr" rid="b51">51</xref>. The second expansion occurred approximately 0.14 Mya during the Last Interglacial period (LIG, 0.12–0.14 Mya)<xref ref-type="bibr" rid="b14">14</xref><xref ref-type="bibr" rid="b29">29</xref><xref ref-type="bibr" rid="b44">44</xref><xref ref-type="bibr" rid="b45">45</xref>, during which the populations recolonized either westward to Himalaya or eastward to the east QTP. These multiple expansions may have reshuffled the genetic distribution in geography, which, together with genetic drift, led to homogeneity within geographic groups. For example, in east QTP, more than half of the populations were fixed with H8, and more than ten populations (>50%) had H3. The repeated range expansions led to the homogeneity within geographic groups, and subsequent secondary contacts may further blur the boundaries among geographic groups, resulting in inconsistency of two major phylogeographic groups (predicted by BARRIER) and three geographic groups. Taken together, these multiple lines of evidence provide a good understanding of the population history of <italic>O. sinensis</italic>; however, more comprehensive knowledge should be gained through the use of multiple markers, such as nuclear DNA markers. Similar pre-LGM demographic expansions were reported for other plants in the QTP, such as <italic>Pedicularis longiflora</italic><xref ref-type="bibr" rid="b36">36</xref>, <italic>Potentilla glabra</italic><xref ref-type="bibr" rid="b31">31</xref> and <italic>Hippophae tibetana</italic><xref ref-type="bibr" rid="b52">52</xref>.</p></sec><sec disp-level="2"><title>Implications for management of biodiversity conservation</title><p>The understanding of the geographical distribution of cpDNA variations may be used for practical applications, such as preventing the loss of genetic resources, conservation management and breeding strategy. More than 12,000 plant species inhabit in the HHM, and a large proportion is endemic. Efforts to identify localities harboring distinct genetic genealogies and regions rich in biodiversity will provide valuable resources for conservation management in this worldwide biodiversity “hotspot” region. Our genetic survey has at least two important implications. First, the three geographic regions (the Himalaya, Hengduan Mountains and east QTP) could be recognized as distinct genetic units, and merit further conservation recognition from a management perspective. Secondly, we found that the Himalaya harbor the highest genetic diversity, whereas the Hengduan Mountains harbor moderate diversity. This poses an urgent need that is of the highest priority for the conservation management. These locations include the Himalaya (populations 1–3), the central Hengduan Mountains (e.g., 10, 13 and 28), the south Hengduan Mountains (e.g., 21–23), and the east QTP (e.g., population 31, 34). In addition, we should note that only a few field samples were obtained in the Himalaya, implying a rather small natural population and a potential loss of genetic diversity of <italic>O. sinensis</italic>. The future conservation policy for <italic>O. sinensis</italic> should focus on protection of regions of distinct lineages and high diversity. The most effective conservation strategy may combine sampling populations containing different genealogies and field protection of those rich biodiversity regions <italic>in situ</italic>. Overall, our findings revealed the genetic structure of <italic>O. sinensis</italic> and tentatively indicated geographic regions with excess genetic diversity. These findings could facilitate short-term and long-term conservation in this worldwide biodiversity “hotspot” region.</p></sec></sec><sec disp-level="1" sec-type="conclusions"><title>Conclusions</title><p>For plants originated in the Miocene on the QTP, their genetic structure could have been affected by the geological changes due to uplifts of QTP and glacial-interglacial oscillations during the Pleistocene since their speciation. Our phylogeographic analysis of <italic>O. sinensis</italic> suggests that Pleistocene refugial isolation may have driven lineage divergence at spatio-temporal timescales and that later interglacial range expansions blurred the phylogeographic boundaries created in the early Pleistocene. Our results provide one of the few detailed maternal evolutionary histories of endemic plants with deep lineages and give important insight into the origin and maintenance of genetic diversity in a global biodiversity “hotspot” region.</p></sec><sec disp-level="1" sec-type="materials|methods"><title>Materials and Methods</title><sec disp-level="2"><title>Plant material</title><p><italic>Oxyria sinensis</italic> is a dioecious perennial herb endemic to the Himalaya-Hengdua Mountains (HHM) that can reproduce sexually via seeds and asexually via rhizomes<xref ref-type="bibr" rid="b53">53</xref>. This plant is characterized by tiny flowers, densely branched panicles, and winged achenes that favor potential dispersal by wind. Samples of <italic>O. sinensis</italic> were collected from 38 populations across its entire range between 2008 and 2011 (<xref ref-type="fig" rid="f1">Fig. 1</xref> and <xref ref-type="supplementary-material" rid="S1">Supplementary Table S1</xref>), representing three broad geographical groups: the Himalaya, the Hengduan Mountains and the east Qinghai-Tibetan Plateau (QTP). A total of 477 individuals were sampled, with an average of 12 individuals per location (<xref ref-type="supplementary-material" rid="S1">Supplementary Table S1</xref>). By comparing records from the Chinese Virtual Herbarium (CVH; http://www.cvh.ac.cn/), our sampling generally reflects the overall distribution of <italic>O. sinensis</italic> (<xref ref-type="fig" rid="f4">Fig. 4a</xref>; <xref ref-type="supplementary-material" rid="S1">Supplementary Table S3</xref>). All samples were stored at -20 °C in the herbarium of the School of Life Science, Yunnan Normal University, P. R. China.</p></sec><sec disp-level="2"><title>DNA extraction, amplification and sequencing</title><p>Genomic DNA was extracted from 20 mg of dried leaf tissue from each <italic>O. sinensis</italic> individual, using a QIAGEN DNeasy Tissue Kit (Qiagen, BOSITE Biology Co. Ltd., Shanghai, China) following the manufacturer’s protocol. We amplified the chloroplast Muturase K (<italic>mat</italic>K) gene using polymerase chain reactions (PCR), with a primer pair of matKAF (Forward, 5’-CTA TAT CCA CTT ATC TTT CAG GAG-3’) and matK8R (Reverse, 5’-AAA GTT CTA GCA CAA GAA AGT CGA-3’). The PCR reactions were performed in a 25 μL reaction volume, containing 1.5 μL genomic DNA, 50 mM Tris-HCL, 1.5 mM MgCl<sub>2</sub>, 250 μg/mL bovine serum albumin (BSA), 0.5 mM dNTPs, 0.2 μM of each primer, and 0.75 units of <italic>Taq</italic> polymerase. Programmed amplification parameters were as follows: 4 min denaturation at 94 °C, then 35 cycles of 50 sec denaturation at 94 °C, 50 sec annealing at 53 °C, and 1 min 30 sec extension at 72 °C, and a final extension of 7 min at 72 °C. PCR products were purified before sequencing to remove excess primers and deoxynucleotide triphosphates using a TIAN quick Midi Purification Kit (Tiangen Biotechnology Co. Ltd., Beijing, China). Sequencing reactions were performed using ABI Prism Sequencing Ready Reaction Kit with the same primers as PCRs, and analyzed on the ABI 3730 genetic analyzer (Applied Biosystems).</p></sec><sec disp-level="2"><title>Phylogenetic analysis and molecular dating</title><p>DNA sequences were aligned using CLUSTAL X<xref ref-type="bibr" rid="b54">54</xref> with default parameters and checked manually. To root a phylogenetic tree, we additionally sequenced the <italic>mat</italic>K gene from a closely related <italic>Oxyria</italic> species (<italic>O. digyna</italic>)<xref ref-type="bibr" rid="b40">40</xref><xref ref-type="bibr" rid="b55">55</xref> and used this sequence as the outgroup for all the <italic>O. sinensis</italic> sequences. All sequences generated in this study have been deposited in GenBank under accession numbers <ext-link ext-link-type="DDBJ/EMBL/GenBank" xlink:href="KJ159010–KJ159025">KJ159010–KJ159025</ext-link>.</p><p>The genealogical relationships among <italic>mat</italic>K haplotypes and divergence times were jointly estimated through a Bayesian Markov Chain Monte Carlo approach implemented in BEAST v.1.7.2<xref ref-type="bibr" rid="b56">56</xref>. Because the substitution rate of chloroplast DNA (cpDNA) for <italic>Oxyria</italic> is unknown, we used a range of mutation rates from 1.01 × 10<sup>−9</sup> to 2.9 × 10<sup>−9</sup> substitution/site/year (s/s/y) that are believed valid for chloroplast DNA of angiosperms overall<xref ref-type="bibr" rid="b57">57</xref><xref ref-type="bibr" rid="b58">58</xref><xref ref-type="bibr" rid="b59">59</xref>. The divergence time of each bifurcating event was estimated under the strict molecular clock model, assuming no rate variation across the tree branches. We used jModeltest v2.1.4<xref ref-type="bibr" rid="b60">60</xref><xref ref-type="bibr" rid="b61">61</xref> to determine the best-fitting DNA substitution model of sequence evolution using the Akaike Information Criterion (AIC) and comparing –<italic>ln</italic> likelihood scores. The GTR + I model was chosen and applied to all Bayesian inferences. Three independent Markov Chain Monte Carlo (MCMC) runs were performed, with each run starting a random tree for 200,000,000 generations, and sampled for every 1000 generations. Convergence of three independent runs was evaluated using Tracer v.1.5 (<ext-link ext-link-type="uri" xlink:href="http://evolve.zoo.ox.ac.uk/">http://evolve.zoo.ox.ac.uk/</ext-link>), with an effective sample size (ESS) for all parameters larger than 100. The results in log files of multiple runs were combined using LogCombiner v.1.5.4<xref ref-type="bibr" rid="b62">62</xref> and analyzed collectively. The first 20% of the generations were discarded as burn-in, and the rest were retained as valid samples for further analysis. A maximum clade credibility tree and node-specific parameters were computed from the sampled data using TreeAnnotator (<ext-link ext-link-type="uri" xlink:href="http://beast.bio.ed.ac.uk/TreeAnnotator/">http://beast.bio.ed.ac.uk/TreeAnnotator/</ext-link>). The tree was plotted using FigTree v.1.31 (<ext-link ext-link-type="uri" xlink:href="http://tree.bio.ed.ac.uk/software/figtree/">http://tree.bio.ed.ac.uk/software/figtree/</ext-link>) and further graphically modified for publication using graphic design software CorelDRAW X6 (Corel Corp., Ottawa, Canada).</p><p>To explore the genealogical relationships and their phylogeographic pattern, we also constructed the parsimonious networks for all <italic>mat</italic>K haplotypes of <italic>O. sinensis</italic> using TCS 1.21<xref ref-type="bibr" rid="b63">63</xref>, under a 95% criterion of statistical parsimony<xref ref-type="bibr" rid="b64">64</xref>.</p></sec><sec disp-level="2"><title>Phylogeographic and population genetic analyses</title><p>The average haplotype diversity within populations (<italic>H</italic><sub>S</sub>) and among populations (<italic>H</italic><sub>T</sub>), and two population differentiation parameters, <italic>G</italic><sub>ST</sub> and <italic>N</italic><sub>ST</sub>, were computed using the program PERMUT with 1000 permutations<xref ref-type="bibr" rid="b65">65</xref> (<ext-link ext-link-type="uri" xlink:href="http://www.pierroton.inra.fr/genetics/labo/Software/PermutCpSSR">http://www.pierroton.inra.fr/genetics/labo/Software/PermutCpSSR</ext-link>). Geographic structure and genetic differentiation within and between geographical regions were statistically tested through analyses of molecular variance (AMOVA)<xref ref-type="bibr" rid="b66">66</xref> in ARELQUIN v3.1.1<xref ref-type="bibr" rid="b67">67</xref>. The haplotype diversity (<italic>H</italic><sub><italic>E</italic></sub>) and nucleotide diversity (π) for geographic populations were also calculated using ARLEQUIN.</p><p>Demographic dynamics through time were inferred using Bayesian skyline plots (BSPs) in BEAST. The BSPs employ a coalescent approach to estimate the effective population size and its changes over piecewise coalescent intervals<xref ref-type="bibr" rid="b68">68</xref>. The BSPs estimates were obtained through MCMC runs in BEAST, using the same running parameters for phylogenetic and clock time inference described above, including a DNA substitution model, starting tree, running length, and burn-in length.</p><p>Historical population expansions were tested using mismatch distribution analysis implemented in ARLEQUIN. The mismatch analysis compares the observed frequency distribution of pairwise differences among haplotypes with the expected from those under the parametric demographic expansion model. A uni-modal distribution of pairwise haplotypes difference is expected if a population has undergone a recent expansion, whereas a multimodal distribution suggests a population in equilibrium<xref ref-type="bibr" rid="b69">69</xref>. Evidence of recent demographic expansion was also examined using Harpending’s Raggedness index<xref ref-type="bibr" rid="b70">70</xref>, Tajima’s <italic>D</italic><xref ref-type="bibr" rid="b71">71</xref> and Fu’s <italic>F</italic><sub>S</sub> test<xref ref-type="bibr" rid="b72">72</xref>. For a population fitting a sudden expansion model, expansion time (<italic>t</italic>, in years) was estimated from the mutation rate (<italic>u</italic>) and number of generations since expansion (τ), using Rogers and Harpending’s<xref ref-type="bibr" rid="b70">70</xref> approach, τ = 2<italic>ut</italic>. The mutation rate <italic>u</italic> was calculated from the DNA substitution rate (<italic>μ</italic>), sequence length (<italic>k</italic>) and the generation time (<italic>g</italic>) of <italic>O. sinensis</italic> using the formula <italic>u</italic> <italic>=</italic> 2<italic>μkg</italic>. A generation time of 2 years and an average <italic>cp</italic>DNA mutation rate of 2 × 10<sup>−9</sup> s/s/y were assumed<xref ref-type="bibr" rid="b57">57</xref><xref ref-type="bibr" rid="b58">58</xref><xref ref-type="bibr" rid="b59">59</xref>.</p><p>BARRIER v2.2<xref ref-type="bibr" rid="b73">73</xref> was used to predict the potential phylogeographic boundaries exhibiting the largest genetic discontinuities among all population pairs. The program uses a genetic matrix of pairwise Fst values from all of the geographic population pairs to construct a geometric network of populations. Putative genetic boundaries were further identified based on Monmonier’s maximum difference algorithm<xref ref-type="bibr" rid="b73">73</xref>.</p></sec><sec disp-level="2"><title>Species distribution modeling</title><p>Species niche modeligng (SDM) was conducted to examine the influence of Pleistocene climate changes on the geographic distribution of <italic>O. sinensis</italic>. Environmental conditions including 19 bioclimate variables were downloaded from the WorldClim database (<ext-link ext-link-type="uri" xlink:href="http://www.worldclim.org">http://www.worldclim.org</ext-link>). Bioclimate layers were obtained in the 2.5 arc-minute resolution for the present and the LGM (ca. 21,000 years ago). Climate data from two general circulation model (GCM) simulations from the Community Climate System Model (CCSM)<xref ref-type="bibr" rid="b74">74</xref><xref ref-type="bibr" rid="b75">75</xref> and the Palaeoclimate Modeling Intercomparison Project (MIROC)<xref ref-type="bibr" rid="b76">76</xref>, were used in this study.</p><p>Distribution records of <italic>O. sinensis</italic> were compiled from our field study and the Chinese Virtual Herbarium (CVH; <ext-link ext-link-type="uri" xlink:href="http://www.cvh.ac.cn/">http://www.cvh.ac.cn/</ext-link>). We obtained 205 hits by searching “<italic>Oxyria Sinensis</italic>” in the CVH (last visit on 10/10/2014), representing historical collections from 1906 to present. Exact geographic information was obtained from 81 records, i.e., the longitude and latitude coordinates were recorded by collectors. For those without geographical coordinates, we mapped their locations to the nearest county. After removing duplicates, we obtained 64 additional unique records from CVH. After merging the CVH and this study, a total of 102 records were included and subjected to the SDM.</p><p>We applied three algorithms to build the SDM – one regression methods: generalized linear models (GLM) and two machine learning methods: random forests (RF) and maximum entropy (MAXENT). Models were created for each algorithm to predict the effects of topography on species distribution. All three algorithms were implemented in “<italic>Biomod2</italic>” (<ext-link ext-link-type="uri" xlink:href="http://cran.r-project.org/web/packages/biomod2/index.html">http://cran.r-project.org/web/packages/biomod2/index.html</ext-link>) packages in R v. 3.0 (<ext-link ext-link-type="uri" xlink:href="http://www.R-project.org/">http://www.R-project.org/</ext-link>). For the MAXENT method, Maxent v 3.3.3k<xref ref-type="bibr" rid="b77">77</xref> was downloaded from <ext-link ext-link-type="uri" xlink:href="http://www.cs.princeton.edu/~schapire/maxent/">http://www.cs.princeton.edu/~schapire/maxent/</ext-link>. Because all three methods require presence-absence data and <italic>O. sinensis</italic> lacks true absences records, we randomly generated 10,000 background pseudo-absence points. All three SDM modeling approaches were combined to construct the realized distribution model using an ensemble forecasting (EM) techniques to account for prediction uncertainty among different algorithms and outperform single models<xref ref-type="bibr" rid="b78">78</xref><xref ref-type="bibr" rid="b79">79</xref>.</p><p>To evaluate model predictive performance, we randomly split the original dataset into two subsets, using 70% (training dataset) to calibrate each model and 30% (testing dataset) for model evaluation. We repeated the processes 30 times to obtain robust estimates of the SDM. The accuracy of the model performance was evaluated using three indices: the receiver operational characteristic (ROC), the True Skill Statistic (TSS) and Cohen’s Kappa (KAPPA). The ROC curve plots with true positive rate of model simulation as a function of the false positive rate<xref ref-type="bibr" rid="b71">71</xref>. The measure of the area under the ROC curve (AUC) measures the area under the curve, providing a threshold independent accuracy index. The AUC values vary from 0 to 1, with a higher value indicating a more fitting performance. An AUC of greater than 0.5 suggests that the given model has non-random discrimination abilities. The TSS statistic measures the agreement between the expected and observed distribution and is not influenced by prevalence. TSS ranges from –1 (perfect disagreement) to 1 (perfect agreement), and a TSS value less than 0 generally indicates the model’s predictive performance is no better than random. The KAPPA statistic is closely related to TSS, taking into account agreement between the observed and the expected distribution by chance.</p><p>We obtained the projections for the LGM by applying the realized models fitted to the present conditions of <italic>O. sinensis</italic>. To compare the similarity of niches (or niche overlap) between the two ecological layers, we computed Schoener’s <italic>D</italic><xref ref-type="bibr" rid="b80">80</xref> and <italic>I</italic> statistics<xref ref-type="bibr" rid="b81">81</xref> using the ENMTOOLS<xref ref-type="bibr" rid="b82">82</xref><xref ref-type="bibr" rid="b83">83</xref>. The values of the two parameters vary between 0 (no overlaps at all) and 1 (complete overlapping), with a fraction number indicating a certain level of similarity in ecological niches. Finally, we visualized the predicted species distribution using DIVA-GIS v.7.5 (<ext-link ext-link-type="uri" xlink:href="http://www.diva-gis.org/">http://www.diva-gis.org/</ext-link>) and R 3.0.</p></sec></sec><sec disp-level="1"><title>Additional Information</title><p><bold>How to cite this article</bold>: Meng, L. <italic>et al.</italic> Refugial isolation and range expansions drive the genetic structure of <italic>Oxyria sinensis</italic> (Polygonaceae) in the Himalaya-Hengduan Mountains. <italic>Sci. Rep.</italic><bold>5</bold>, 10396; doi: 10.1038/srep10396 (2015).</p></sec><sec sec-type="supplementary-material" id="S1"><title>Supplementary Material</title><supplementary-material id="d33e24" content-type="local-data"><caption><title>Supplementary Information</title></caption><media xlink:href="srep10396-s1.doc"></media></supplementary-material></sec></body><back><ack><p>We thank Dr. Dengmei Fan and Mr. Zhiguo Xue for providing the map and Ms. Yun Wang and Dr. Xinmin Tian for their help in field collections. We thank Dr. Damien Georges for his help on the species distribution modeling in R “Biomod2” packages. We also thank Sarah Jaslow and Benjamin Redelings for English editing. This study is funded by the Natural Science Foundation of China (31160084, 41101058) and Ph.D. Programs Foundation of Ministry of Education of China (20115303120001)</p></ack><ref-list><ref id="b1"><mixed-citation publication-type="journal"><name><surname>Hewitt</surname><given-names>G.</given-names></name><article-title>The genetic legacy of the Quaternary ice ages</article-title>. <source>Nature</source><volume>405</volume>, <fpage>907</fpage>–<lpage>913</lpage> (<year>2000</year>).<pub-id pub-id-type="pmid">10879524</pub-id></mixed-citation></ref><ref id="b2"><mixed-citation publication-type="journal"><name><surname>Hewitt</surname><given-names>G. M.</given-names></name><article-title>Some genetic consequences of ice ages, and their role in divergence and speciation</article-title>. <source>Biol. J. Linn. Soc.</source><volume>58</volume>, <fpage>247</fpage>–<lpage>276</lpage> (<year>1996</year>).</mixed-citation></ref><ref id="b3"><mixed-citation publication-type="journal"><name><surname>Avise</surname><given-names>J. C.</given-names></name><article-title>Phylogeography: retrospect and prospect</article-title>. <source>J. Biogeogr.</source><volume>36</volume>, <fpage>3</fpage>–<lpage>15</lpage> (<year>2009</year>).</mixed-citation></ref><ref id="b4"><mixed-citation publication-type="journal"><name><surname>Myers</surname><given-names>N.</given-names></name>, <name><surname>Mittermeier</surname><given-names>R. A.</given-names></name>, <name><surname>Mittermeier</surname><given-names>C. G.</given-names></name>, <name><surname>da Fonseca</surname><given-names>G. A.</given-names></name> & <name><surname>Kent</surname><given-names>J.</given-names></name><article-title>Biodiversity hotspots for conservation priorities</article-title>. <source>Nature</source><volume>403</volume>, <fpage>853</fpage>–<lpage>858</lpage> (<year>2000</year>).<pub-id pub-id-type="pmid">10706275</pub-id></mixed-citation></ref><ref id="b5"><mixed-citation publication-type="journal"><name><surname>Li</surname><given-names>X. W.</given-names></name> & <name><surname>Li</surname><given-names>J.</given-names></name><article-title>A preliminary floristic study on the seed plants from the region of Hengduan Mountain</article-title>. <source>Acta. Bot. Yunnan</source> . <volume>15</volume>, <fpage>217</fpage>–<lpage>231</lpage> (<year>1993</year>).</mixed-citation></ref><ref id="b6"><mixed-citation publication-type="journal"><name><surname>Mittermeier</surname><given-names>R.</given-names></name><italic>et al.</italic><source>Hotspots revisited: Earth’s biologically richest and most endangered terrestrial ecoregions</source> . (CEMEX & Agrupacion Sierra Madre, Mexico, <year>2004</year>).</mixed-citation></ref><ref id="b7"><mixed-citation publication-type="journal"><name><surname>Wilson</surname><given-names>E. O.</given-names></name><source>The diversity of life</source> . Harvard University Press, Boston, <year>1992</year>).</mixed-citation></ref><ref id="b8"><mixed-citation publication-type="journal"><name><surname>Wu</surname><given-names>Z.</given-names></name><article-title>Delineation and unique features of the Sino-Japanese floristic region</article-title>. In <name><surname>Boufford</surname><given-names>D. E.</given-names></name> & <name><surname>Ohba</surname><given-names>H.</given-names></name> (eds) <source>Sino-Japanese Flora-Its Characteristics and Diversification</source> (University of Tokyo Bulletin, Tokyo, <year>1998</year>).</mixed-citation></ref><ref id="b9"><mixed-citation publication-type="journal"><name><surname>Gao</surname><given-names>L. M.</given-names></name><italic>et al.</italic><article-title>High variation and strong phylogeographic pattern among cpDNA haplotypes in <italic>Taxus wallichiana</italic> (Taxaceae) in China and North Vietnam</article-title>. <source>Mol. Ecol.</source><volume>16</volume>, <fpage>4684</fpage>–<lpage>4698</lpage> (<year>2007</year>).<pub-id pub-id-type="pmid">17908214</pub-id></mixed-citation></ref><ref id="b10"><mixed-citation publication-type="journal"><name><surname>Li</surname><given-names>Y.</given-names></name><italic>et al.</italic><article-title>Glacial survival east and west of the ‘Mekong–Salween Divide’in the Himalaya–Hengduan Mountains region as revealed by AFLPs and cpDNA sequence variation in <italic>Sinopodophyllum hexandrum</italic> (Berberidaceae)</article-title>. <source>Mol. Phylogenet. Evol.</source><volume>59</volume>, <fpage>412</fpage>–<lpage>424</lpage> (<year>2011</year>).<pub-id pub-id-type="pmid">21296173</pub-id></mixed-citation></ref><ref id="b11"><mixed-citation publication-type="journal"><name><surname>Qiu</surname><given-names>Y. X.</given-names></name>, <name><surname>Fu</surname><given-names>C. X.</given-names></name> & <name><surname>Comes</surname><given-names>H. P.</given-names></name><article-title>Plant molecular phylogeography in China and adjacent regions: tracing the genetic imprints of Quaternary climate and environmental change in the world’s most diverse temperate flora</article-title>. <source>Mol. Phylogenet. Evol.</source><volume>59</volume>, <fpage>225</fpage>–<lpage>244</lpage> (<year>2011</year>).<pub-id pub-id-type="pmid">21292014</pub-id></mixed-citation></ref><ref id="b12"><mixed-citation publication-type="journal"><name><surname>Liu</surname><given-names>Y.</given-names></name>, <name><surname>Yang</surname><given-names>S.</given-names></name>, <name><surname>Ji</surname><given-names>P.</given-names></name> & <name><surname>Gao</surname><given-names>L.</given-names></name><article-title>Phylogeography of <italic>Camellia taliensis</italic> (Theaceae) inferred from chloroplast and nuclear DNA: insights into evolutionary history and conservation</article-title>. <source>BMC Evol. Biol.</source><volume>12</volume>, <fpage>92</fpage> (<year>2012</year>).<pub-id pub-id-type="pmid">22716114</pub-id></mixed-citation></ref><ref id="b13"><mixed-citation publication-type="journal"><name><surname>Zhang</surname><given-names>T. C.</given-names></name>, <name><surname>Comes</surname><given-names>H. P.</given-names></name> & <name><surname>Sun</surname><given-names>H.</given-names></name><article-title>Chloroplast phylogeography of <italic>Terminalia franchetii</italic> (Combretaceae) from the eastern Sino-Himalayan region and its correlation with historical river capture events</article-title>. <source>Mol. Phylogenet. Evol.</source><volume>60</volume>, <fpage>1</fpage>–<lpage>12</lpage> (<year>2011</year>).<pub-id pub-id-type="pmid">21545839</pub-id></mixed-citation></ref><ref id="b14"><mixed-citation publication-type="journal"><name><surname>Shi</surname><given-names>Y. F.</given-names></name><source>Uplift and Environmental Changes of Qinghai-Xizang (Tibetan) Plateau in the Late Cenozoic</source> . (Guangdong Science & Technology Press, Guangzhou, <year>1998</year>).</mixed-citation></ref><ref id="b15"><mixed-citation publication-type="journal"><name><surname>Clark</surname><given-names>M.</given-names></name><italic>et al.</italic><article-title>Surface uplift, tectonics, and erosion of eastern Tibet from large-scale drainage patterns</article-title>. <source>Tectonics</source><volume>23</volume>, <fpage>1</fpage>–<lpage>6</lpage> (<year>2004</year>).</mixed-citation></ref><ref id="b16"><mixed-citation publication-type="journal"><name><surname>Liu</surname><given-names>J.</given-names></name><italic>et al.</italic><article-title>Geological and ecological factors drive cryptic speciation of yews in a biodiversity hotspot</article-title>. <source>New Phytol.</source><volume>199</volume>, <fpage>1093</fpage>–<lpage>1108</lpage> (<year>2013</year>).<pub-id pub-id-type="pmid">23718262</pub-id></mixed-citation></ref><ref id="b17"><mixed-citation publication-type="journal"><name><surname>Ward</surname><given-names>F. K.</given-names></name><article-title>The Mekong-Salween Divide as a geographical barrier</article-title>. <source>Geogr. J.</source><volume>58</volume>, <fpage>49</fpage>–<lpage>56</lpage> (<year>1921</year>).</mixed-citation></ref><ref id="b18"><mixed-citation publication-type="journal"><name><surname>Zhang</surname><given-names>T.</given-names></name> & <name><surname>Sun</surname><given-names>H.</given-names></name><article-title>Phylogeographic structure of <italic>Terminalia franchetii</italic> (Combretaceae) in southwest China and its implications for drainage geological history</article-title>. <source>J. Plant Res.</source><volume>124</volume>, <fpage>63</fpage>–<lpage>73</lpage> (<year>2011</year>).<pub-id pub-id-type="pmid">20625791</pub-id></mixed-citation></ref><ref id="b19"><mixed-citation publication-type="journal"><name><surname>Li</surname><given-names>X. W.</given-names></name> & <name><surname>Li</surname><given-names>J.</given-names></name><article-title>On the validity of Tanaka line and its significance viewed from the distribution of eastern Asiatic genera in Yunnan</article-title>. <source>Plant Divers. Resour</source> . <volume>14</volume>, <fpage>1</fpage>–<lpage>12</lpage> (<year>1992</year>).</mixed-citation></ref><ref id="b20"><mixed-citation publication-type="journal"><name><surname>Li</surname><given-names>X.W.</given-names></name> & <name><surname>Li</surname><given-names>J.</given-names></name><article-title>The Tanaka-Kaiyong line–an important floristic line for the study of the flora of East Asia</article-title>. <source>Ann. Mo. Bot. Gard.</source><volume>84</volume>, <fpage>888</fpage>–<lpage>892</lpage> (<year>1997</year>).</mixed-citation></ref><ref id="b21"><mixed-citation publication-type="journal"><name><surname>Zhang</surname><given-names>Z.Y.</given-names></name><italic>et al.</italic><article-title>Alkaloid polymorphism and ITS sequence variation in the <italic>Spiraea japonicacomplex</italic> (Rosaceae) in China: traces of the biological effects of the Himalaya-Tibet Plateau uplift</article-title>. <source><italic>Am. J. Bot</italic>.</source><volume>93</volume>, <fpage>762</fpage>–<lpage>769</lpage> (<year>2006</year>).<pub-id pub-id-type="pmid">21642139</pub-id></mixed-citation></ref><ref id="b22"><mixed-citation publication-type="journal"><name><surname>Qiu</surname><given-names>Y. X.</given-names></name>, <name><surname>Guan</surname><given-names>B. C.</given-names></name>, <name><surname>Fu</surname><given-names>C. X.</given-names></name> & <name><surname>Comes</surname><given-names>H. P.</given-names></name><article-title>Did glacials and/or interglacials promote allopatric incipient speciation in East Asian temperate plants? Phylogeographic and coalescent analyses on refugial isolation and divergence in <italic>Dysosma versipellis</italic></article-title>. <source>Mol. Phylogenet. Evol.</source><volume>51</volume>, <fpage>281</fpage>–<lpage>293</lpage>(<year>2009</year>)<pub-id pub-id-type="pmid">19405195</pub-id></mixed-citation></ref><ref id="b23"><mixed-citation publication-type="journal"><name><surname>Liu</surname><given-names>J. Q.</given-names></name>, <name><surname>Gao</surname><given-names>T. G.</given-names></name>, <name><surname>Chen</surname><given-names>Z. D.</given-names></name> & <name><surname>Lu</surname><given-names>A. M.</given-names></name><article-title>Molecular phylogeny and biogeography of the Qinghai-Tibet Plateau endemic Nannoglottis (Asteraceae)</article-title>. <source>Mol. Phylogenet. Evol.</source><volume>23</volume>, <fpage>307</fpage>–<lpage>325</lpage> (<year>2002</year>).<pub-id pub-id-type="pmid">12099790</pub-id></mixed-citation></ref><ref id="b24"><mixed-citation publication-type="journal"><name><surname>Liu</surname><given-names>J. Q.</given-names></name>, <name><surname>Wang</surname><given-names>Y. J.</given-names></name>, <name><surname>Wang</surname><given-names>A. L.</given-names></name>, <name><surname>Hideaki</surname><given-names>O.</given-names></name> & <name><surname>Abbott</surname><given-names>R. J.</given-names></name><article-title>Radiation and diversification within the <italic>Ligularia–Cremanthodium–Parasenecio</italic> complex (Asteraceae) triggered by uplift of the Qinghai-Tibetan Plateau</article-title>. <source><italic>Mol. Phylogenet. Evol</italic>.</source><volume>38</volume>, <fpage>31</fpage>–<lpage>49</lpage> (<year>2006</year>).<pub-id pub-id-type="pmid">16290033</pub-id></mixed-citation></ref><ref id="b25"><mixed-citation publication-type="journal"><name><surname>Wang</surname><given-names>L.</given-names></name>, <name><surname>Schneider</surname><given-names>H.</given-names></name>, <name><surname>Zhang</surname><given-names>X.-C.</given-names></name> & <name><surname>Xiang</surname><given-names>Q.-P.</given-names></name><article-title>The rise of the Himalaya enforced the diversification of SE Asian ferns by altering the monsoon regimes</article-title>. <source><italic>BMC Plant Biol</italic>.</source><volume>12</volume>, <fpage>210</fpage> (<year>2012</year>).<pub-id pub-id-type="pmid">23140168</pub-id></mixed-citation></ref><ref id="b26"><mixed-citation publication-type="journal"><name><surname>Wang</surname><given-names>L.</given-names></name><italic>et al.</italic><article-title>History and evolution of alpine plants endemic to the Qinghai-Tibetan Plateau: <italic>Aconitum gymnandrum</italic> (Ranunculaceae)</article-title>. <source>Mol. Ecol.</source><volume>18</volume>, <fpage>709</fpage>–<lpage>721</lpage> (<year>2009</year>).<pub-id pub-id-type="pmid">19175501</pub-id></mixed-citation></ref><ref id="b27"><mixed-citation publication-type="journal"><name><surname>Sun</surname><given-names>Y.</given-names></name>, <name><surname>Wang</surname><given-names>A.</given-names></name>, <name><surname>Wan</surname><given-names>D.</given-names></name>, <name><surname>Wang</surname><given-names>Q.</given-names></name> & <name><surname>Liu</surname><given-names>J.</given-names></name><article-title>Rapid radiation of <italic>Rheum</italic> (Polygonaceae) and parallel evolution of morphological traits</article-title>. <source>Mol. Phylogenet. Evol.</source><volume>63</volume>, <fpage>150</fpage>–<lpage>158</lpage> (<year>2012</year>).<pub-id pub-id-type="pmid">22266181</pub-id></mixed-citation></ref><ref id="b28"><mixed-citation publication-type="journal"><name><surname>Chen</surname><given-names>S.</given-names></name><italic>et al.</italic><article-title>Phylogeographic analysis reveals significant spatial genetic structure of <italic>Incarvillea sinensis</italic> as a product of mountain building</article-title>. <source>BMC Plant Biol.</source><volume>12</volume>, <fpage>58</fpage> (<year>2012</year>).<pub-id pub-id-type="pmid">22546007</pub-id></mixed-citation></ref><ref id="b29"><mixed-citation publication-type="journal"><name><surname>Shi</surname><given-names>Y.</given-names></name><article-title>Characteristics of late Quaternary monsoonal glaciation on the Tibetan Plateau and in East Asia</article-title>. <source><italic>Quaternary Int</italic>.</source><volume>97</volume>, <fpage>79</fpage>–<lpage>91</lpage> (<year>2002</year>).</mixed-citation></ref><ref id="b30"><mixed-citation publication-type="journal"><name><surname>Shi</surname><given-names>Y. F.</given-names></name>, <name><surname>Zheng</surname><given-names>B. X.</given-names></name> & <name><surname>Li</surname><given-names>S. J.</given-names></name><article-title>Last glaciation and maximum glaciation in Qinghai-Xizang Plateau: A controversy to M. Kuhle’s Ice Sheet Hypothesis</article-title>. <source>Chinese Geog. Sci</source> . <volume>2</volume>, <fpage>293</fpage>–<lpage>311</lpage> (<year>1990</year>).</mixed-citation></ref><ref id="b31"><mixed-citation publication-type="journal"><name><surname>Wang</surname><given-names>L. Y.</given-names></name>, <name><surname>Ikeda</surname><given-names>H.</given-names></name>, <name><surname>Liu</surname><given-names>T. L.</given-names></name>, <name><surname>Wang</surname><given-names>Y. J.</given-names></name> & <name><surname>Liu</surname><given-names>J. Q.</given-names></name><article-title>Repeated range expansion and glacial endurance of <italic>Potentilla glabra</italic> (Rosaceae) in the Qinghai-Tibetan plateau</article-title>. <source>J. Integr. Plant Biol.</source><volume>51</volume>, <fpage>698</fpage>–<lpage>706</lpage> (<year>2009</year>).<pub-id pub-id-type="pmid">19566648</pub-id></mixed-citation></ref><ref id="b32"><mixed-citation publication-type="journal"><name><surname>Opgenoorth</surname><given-names>L.</given-names></name><italic>et al.</italic><article-title>Tree endurance on the Tibetan Plateau marks the world’s highest known tree line of the Last Glacial Maximum</article-title>. <source>New Phytol.</source><volume>185</volume>, <fpage>332</fpage>–<lpage>342</lpage> (<year>2010</year>).<pub-id pub-id-type="pmid">19761444</pub-id></mixed-citation></ref><ref id="b33"><mixed-citation publication-type="journal"><name><surname>Jia</surname><given-names>D. R.</given-names></name>, <name><surname>Liu</surname><given-names>T. L.</given-names></name>, <name><surname>Wang</surname><given-names>L. Y.</given-names></name>, <name><surname>Zhou</surname><given-names>D. W.</given-names></name> & <name><surname>Liu</surname><given-names>J. Q.</given-names></name><article-title>Evolutionary history of an alpine shrub <italic>Hippophae tibetana</italic> (Elaeagnaceae): allopatric divergence and regional expansion</article-title>. <source>Biol. J. Linn. Soc.</source><volume>102</volume>, <fpage>37</fpage>–<lpage>50</lpage> (<year>2011</year>).</mixed-citation></ref><ref id="b34"><mixed-citation publication-type="journal"><name><surname>Zhang</surname><given-names>Q.</given-names></name>, <name><surname>Chiang</surname><given-names>T. Y.</given-names></name>, <name><surname>George</surname><given-names>M.</given-names></name>, <name><surname>Liu</surname><given-names>J. Q.</given-names></name> & <name><surname>Abbott</surname><given-names>R. J.</given-names></name><article-title>Phylogeography of the Qinghai-Tibetan Plateau endemic <italic>Juniperus przewalskii</italic> (Cupressaceae) inferred from chloroplast DNA sequence variation</article-title>. <source>Mol. Ecol.</source><volume>14</volume>, <fpage>3513</fpage>–<lpage>3524</lpage> (<year>2005</year>).<pub-id pub-id-type="pmid">16156819</pub-id></mixed-citation></ref><ref id="b35"><mixed-citation publication-type="journal"><name><surname>Meng</surname><given-names>L. H.</given-names></name><italic>et al.</italic><article-title>Mitochondrial and chloroplast phylogeography of <italic>Picea crassifolia</italic> Kom. (Pinaceae) in the Qinghai-Tibetan Plateau and adjacent highlands</article-title>. <source>Mol. Ecol.</source><volume>16</volume>, <fpage>4128</fpage>–<lpage>4137</lpage> (<year>2007</year>).<pub-id pub-id-type="pmid">17894761</pub-id></mixed-citation></ref><ref id="b36"><mixed-citation publication-type="journal"><name><surname>Yang</surname><given-names>F. S.</given-names></name>, <name><surname>Li</surname><given-names>Y. F.</given-names></name>, <name><surname>Ding</surname><given-names>X.</given-names></name> & <name><surname>Wang</surname><given-names>X. Q.</given-names></name><article-title>Extensive population expansion of <italic>Pedicularis longiflora</italic> (Orobanchaceae) on the Qinghai-Tibetan Plateau and its correlation with the Quaternary climate change</article-title>. <source>Mol. Ecol.</source><volume>17</volume>, <fpage>5135</fpage>–<lpage>5145</lpage> (<year>2008</year>).<pub-id pub-id-type="pmid">19017264</pub-id></mixed-citation></ref><ref id="b37"><mixed-citation publication-type="journal"><name><surname>Fan</surname><given-names>D. M.</given-names></name><italic>et al.</italic><article-title>Phylogeography of <italic>Sophora davidii</italic> (Leguminosae) across the ‘Tanaka-Kaiyong Line’, an important phytogeographic boundary in Southwest China</article-title>. <source>Mol. Ecol.</source><volume>22</volume>, <fpage>4270</fpage>–<lpage>4288</lpage> (<year>2013</year>).<pub-id pub-id-type="pmid">23927411</pub-id></mixed-citation></ref><ref id="b38"><mixed-citation publication-type="journal"><name><surname>Cun</surname><given-names>Y. Z.</given-names></name> & <name><surname>Wang</surname><given-names>X. Q.</given-names></name><article-title>Plant recolonization in the Himalaya from the southeastern Qinghai-Tibetan Plateau: Geographical isolation contributed to high population differentiation</article-title>. <source>Mol. Phylogenet. Evol.</source><volume>56</volume>, <fpage>972</fpage>–<lpage>982</lpage> (<year>2010</year>).<pub-id pub-id-type="pmid">20472080</pub-id></mixed-citation></ref><ref id="b39"><mixed-citation publication-type="journal"><name><surname>Avise</surname><given-names>J. C.</given-names></name><italic>et al.</italic><article-title>Interspecific phylogeny: the mitochondrial DNA bridge between population genetics and systematics</article-title>. <source><italic>Ann. Rev. Ecol. Syst</italic>.</source><volume>18</volume>, <fpage>489</fpage>–<lpage>522</lpage> (<year>1987</year>).</mixed-citation></ref><ref id="b40"><mixed-citation publication-type="journal"><name><surname>Liu</surname><given-names>J. Q.</given-names></name>, <name><surname>Sun</surname><given-names>Y. S.</given-names></name>, <name><surname>Ge</surname><given-names>X. J.</given-names></name>, <name><surname>Gao</surname><given-names>L. M.</given-names></name> & <name><surname>Qiu</surname><given-names>Y. X.</given-names></name><article-title>Phylogeographic studies of plants in China: Advances in the past and directions in the future</article-title>. <source>J. Syst. Evol.</source><volume>50</volume>, <fpage>267</fpage>–<lpage>275</lpage> (<year>2012</year>).</mixed-citation></ref><ref id="b41"><mixed-citation publication-type="journal"><name><surname>Kuhle</surname><given-names>M.</given-names></name><article-title>The Pleistocene glaciation of Tibet and the onset of ice ages: an autocycle hypothesis</article-title>. <source>GeoJournal</source><volume>17</volume>, <fpage>457</fpage>–<lpage>511</lpage> (<year>1988</year>).</mixed-citation></ref><ref id="b42"><mixed-citation publication-type="journal"><name><surname>Ehlers</surname><given-names>J.</given-names></name> & <name><surname>Gibbard</surname><given-names>P. L.</given-names></name><article-title>Quaternary glaciations–extent and chronology</article-title>. <source>Part III: South America, Asia, Africa, Australasia, Antarctica. Amsterdam</source> . (Elsevier, Amsterdam, <year>2004</year>).</mixed-citation></ref><ref id="b43"><mixed-citation publication-type="journal"><name><surname>Andersen</surname><given-names>B. G.</given-names></name> & <name><surname>Borns</surname><given-names>H. W.</given-names></name><source>The Ice Age world: an introduction to Quaternary history and research with emphasis on North America and northern Europe during the last 2.5 million years</source> . (Scandinavian University Press, Oslo, <year>1994</year>).</mixed-citation></ref><ref id="b44"><mixed-citation publication-type="journal"><name><surname>Zhou</surname><given-names>S.</given-names></name> & <name><surname>Li</surname><given-names>J.</given-names></name><article-title>The sequence of Quaternary glaciation in the Bayan Har Mountains</article-title>. <source><italic>Quaternary Int</italic>.</source><volume>45</volume>, <fpage>135</fpage>–<lpage>142</lpage> (<year>1998</year>).</mixed-citation></ref><ref id="b45"><mixed-citation publication-type="journal"><name><surname>Zheng</surname><given-names>B.</given-names></name>, <name><surname>Xu</surname><given-names>Q.</given-names></name> & <name><surname>Shen</surname><given-names>Y.</given-names></name><article-title>The relationship between climate change and Quaternary glacial cycles on the Qinghai–Tibetan Plateau: review and speculation</article-title>. <source>Quaternary Int.</source><volume>97</volume>, <fpage>93</fpage>–<lpage>101</lpage> (<year>2002</year>).</mixed-citation></ref><ref id="b46"><mixed-citation publication-type="journal"><name><surname>Owen</surname><given-names>L. A.</given-names></name><italic>et al.</italic><article-title>Climatic and topographic controls on the style and timing of late Quaternary glaciation throughout Tibet and the Himalaya defined by 10 Be cosmogenic radionuclide surface exposure dating</article-title>. <source>Quaternary Sci. Rev.</source><volume>24</volume>, <fpage>1391</fpage>–<lpage>1411</lpage> (<year>2005</year>).</mixed-citation></ref><ref id="b47"><mixed-citation publication-type="journal"><name><surname>Wu</surname><given-names>Y.</given-names></name><italic>et al.</italic><article-title>Quaternary geomorphological evolution of the Kunlun Pass area and uplift of the Qinghai-Xizang (Tibet) Plateau</article-title>. <source>Geomorphology</source><volume>36</volume>, <fpage>203</fpage>–<lpage>216</lpage> (<year>2001</year>).</mixed-citation></ref><ref id="b48"><mixed-citation publication-type="journal"><name><surname>Li</surname><given-names>B. Y.</given-names></name><source>Quaternary glacial distribution map of Qinghai-Xizang (Tibet) plateau</source> . (Science Press, Beijing, <year>1991)</year>.</mixed-citation></ref><ref id="b49"><mixed-citation publication-type="journal"><name><surname>Owen</surname><given-names>L. A.</given-names></name><article-title>Latest Pleistocene and Holocene glacier fluctuations in the Himalaya and Tibet</article-title>. <source>Quaternary Sci. Rev.</source><volume>28</volume>, <fpage>2150</fpage>–<lpage>2164</lpage> (<year>2009</year>).</mixed-citation></ref><ref id="b50"><mixed-citation publication-type="journal"><name><surname>Zhuo</surname><given-names>Z.</given-names></name>, <name><surname>Baoyin</surname><given-names>Y.</given-names></name> & <name><surname>Petit-Maire</surname><given-names>N.</given-names></name><article-title>Paleoenvironments in China during the Last Glacial Maximum and the Holocene optimum</article-title>. <source>Episodes</source><volume>21</volume>, <fpage>152</fpage>–<lpage>158</lpage> (<year>1998</year>).</mixed-citation></ref><ref id="b51"><mixed-citation publication-type="journal"><name><surname>Zhang</surname><given-names>D. F.</given-names></name>, <name><surname>Fengquan</surname><given-names>L.</given-names></name> & <name><surname>Jianmin</surname><given-names>B.</given-names></name><article-title>Eco-environmental effects of the Qinghai-Tibet Plateau uplift during the Quaternary in China</article-title>. <source>Environ. Geol.</source><volume>39</volume>, <fpage>1352</fpage>–<lpage>1358</lpage> (<year>2000</year>).</mixed-citation></ref><ref id="b52"><mixed-citation publication-type="journal"><name><surname>Jia</surname><given-names>D. R.</given-names></name><italic>et al.</italic><article-title>Out of the Qinghai-Tibet Plateau: evidence for the origin and dispersal of Eurasian temperate plants from a phylogeographic study of <italic>Hippophae rhamnoides</italic> (Elaeagnaceae)</article-title>. <source>New Phytol.</source><volume>194</volume>, <fpage>1123</fpage>–<lpage>1133</lpage> (<year>2012</year>).<pub-id pub-id-type="pmid">22432741</pub-id></mixed-citation></ref><ref id="b53"><mixed-citation publication-type="journal"><name><surname>Zhao</surname><given-names>F.</given-names></name> & <name><surname>Yang</surname><given-names>Y. P.</given-names></name><article-title>Reproductive allocation in a dioecious perennial <italic>Oxyria sinensis</italic> (Polygonaceae) along altitudinal gradients</article-title>. <source>J. Syst Evol.</source><volume>46</volume>, <fpage>830</fpage>–<lpage>835</lpage> (<year>2008</year>).</mixed-citation></ref><ref id="b54"><mixed-citation publication-type="journal"><name><surname>Thompson</surname><given-names>J. D.</given-names></name>, <name><surname>Gibson</surname><given-names>T. J.</given-names></name>, <name><surname>Plewniak</surname><given-names>F.</given-names></name>, <name><surname>Jeanmougin</surname><given-names>F.</given-names></name> & <name><surname>Higgins</surname><given-names>D. G.</given-names></name><article-title>The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools</article-title>. <source>Nucleic. Acids. Res.</source><volume>25</volume>, <fpage>4876</fpage>–<lpage>4882</lpage> (<year>1997</year>).<pub-id pub-id-type="pmid">9396791</pub-id></mixed-citation></ref><ref id="b55"><mixed-citation publication-type="journal"><name><surname>Sanchez</surname><given-names>A.</given-names></name>, <name><surname>Schuster</surname><given-names>T. M.</given-names></name> & <name><surname>Kron</surname><given-names>K. A.</given-names></name><article-title>A Large-Scale Phylogeny of Polygonaceae Based on Molecular Data</article-title>. <source>Int. J. Plant Sci.</source><volume>170</volume>, <fpage>1044</fpage>–<lpage>1055</lpage> (<year>2009</year>).</mixed-citation></ref><ref id="b56"><mixed-citation publication-type="journal"><name><surname>Drummond</surname><given-names>A. J.</given-names></name>, <name><surname>Suchard</surname><given-names>M. A.</given-names></name>, <name><surname>Xie</surname><given-names>D.</given-names></name> & <name><surname>Rambaut</surname><given-names>A.</given-names></name><article-title>Bayesian phylogenetics with BEAUti and the BEAST 1.7</article-title>. <source>Mol. Biol. Evol.</source><volume>29</volume>, <fpage>1969</fpage>–<lpage>1973</lpage> (<year>2012</year>).<pub-id pub-id-type="pmid">22367748</pub-id></mixed-citation></ref><ref id="b57"><mixed-citation publication-type="journal"><name><surname>Muse</surname><given-names>S. V.</given-names></name><article-title>Examining rates and patterns of nucleotide substitution in plants</article-title>. <source>Plant Mol. Biol.</source><volume>42</volume>, <fpage>25</fpage>–<lpage>43</lpage> (<year>2000</year>).<pub-id pub-id-type="pmid">10688129</pub-id></mixed-citation></ref><ref id="b58"><mixed-citation publication-type="journal"><name><surname>Richardson</surname><given-names>J. E.</given-names></name>, <name><surname>Pennington</surname><given-names>R. T.</given-names></name>, <name><surname>Pennington</surname><given-names>T. D.</given-names></name> & <name><surname>Hollingsworth</surname><given-names>P.M.</given-names></name><article-title>Rapid diversification of a species-rich genus of neotropical rain forest trees</article-title>. <source>Science</source><volume>293</volume>, <fpage>2242</fpage>–<lpage>2245</lpage> (<year>2001</year>).<pub-id pub-id-type="pmid">11567135</pub-id></mixed-citation></ref><ref id="b59"><mixed-citation publication-type="journal"><name><surname>Wolfe</surname><given-names>K. H.</given-names></name>, <name><surname>Li</surname><given-names>W. H.</given-names></name> & <name><surname>Sharp</surname><given-names>P. M.</given-names></name><article-title>Rates of nucleotide substitution vary greatly among plant mitochondrial, chloroplast, and nuclear DNAs</article-title>. <source>Proc. Nat. Acad. Sci., USA</source><volume>84</volume>, <fpage>9054</fpage>–<lpage>9058</lpage> (<year>1987</year>).<pub-id pub-id-type="pmid">3480529</pub-id></mixed-citation></ref><ref id="b60"><mixed-citation publication-type="journal"><name><surname>Guindon</surname><given-names>S.</given-names></name> & <name><surname>Gascuel</surname><given-names>O.</given-names></name><article-title>A simple, fast and accurate method to estimate large phylogenies by maximum-likelihood</article-title>. <source>Syst. Biol.</source><volume>52</volume>, <fpage>696</fpage>–<lpage>704</lpage> (<year>2003</year>).<pub-id pub-id-type="pmid">14530136</pub-id></mixed-citation></ref><ref id="b61"><mixed-citation publication-type="journal"><name><surname>Darriba</surname><given-names>D.</given-names></name><name><surname>Taboada</surname><given-names>G.L.</given-names></name>, <name><surname>Doallo</surname><given-names>R.</given-names></name>& <name><surname>Posada</surname><given-names>D.</given-names></name><article-title>jModelTest 2: more models, new heuristics and parallel computing</article-title>. <source>Nat. Methods</source><volume>9</volume>, <fpage>772</fpage> (<year>2012</year>).<pub-id pub-id-type="pmid">22847109</pub-id></mixed-citation></ref><ref id="b62"><mixed-citation publication-type="journal"><name><surname>Rambaut</surname><given-names>A.</given-names></name> & <name><surname>Drummond</surname><given-names>A. J.</given-names></name><source>LogCombiner v1. 5.4 MCMC Output Combiner</source> . (Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, <year>2010)</year>.</mixed-citation></ref><ref id="b63"><mixed-citation publication-type="journal"><name><surname>Clement</surname><given-names>M.</given-names></name>, <name><surname>Posada</surname><given-names>D.</given-names></name> & <name><surname>Crandall</surname><given-names>K. A.</given-names></name><article-title>TCS: a computer program to estimate gene genealogies</article-title>. <source>Mol. Ecol.</source><volume>9</volume>, <fpage>1657</fpage>–<lpage>1659</lpage> (<year>2000</year>).<pub-id pub-id-type="pmid">11050560</pub-id></mixed-citation></ref><ref id="b64"><mixed-citation publication-type="journal"><name><surname>Templeton</surname><given-names>A. R.</given-names></name><article-title>Human origins and analysis of mitochondrial DNA sequences</article-title>. <source>Science</source><volume>255</volume>, <fpage>737</fpage> (<year>1992</year>).<pub-id pub-id-type="pmid">1590849</pub-id></mixed-citation></ref><ref id="b65"><mixed-citation publication-type="journal"><name><surname>Pons</surname><given-names>O.</given-names></name> & <name><surname>Petit</surname><given-names>R. J.</given-names></name><article-title>Measuring and testing genetic differentiation with ordered versus unordered alleles</article-title>. <source>Genetics</source><volume>144</volume>, <fpage>1237</fpage>–<lpage>1245</lpage> (<year>1996</year>).<pub-id pub-id-type="pmid">8913764</pub-id></mixed-citation></ref><ref id="b66"><mixed-citation publication-type="journal"><name><surname>Excoffier</surname><given-names>L.</given-names></name>, <name><surname>Smouse</surname><given-names>P. E.</given-names></name> & <name><surname>Quattro</surname><given-names>J. M.</given-names></name><article-title>Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data</article-title>. <source>Genetics</source><volume>131</volume>, <fpage>479</fpage>–<lpage>491</lpage> (<year>1992</year>).<pub-id pub-id-type="pmid">1644282</pub-id></mixed-citation></ref><ref id="b67"><mixed-citation publication-type="journal"><name><surname>Excoffier</surname><given-names>L.</given-names></name>, <name><surname>Laval</surname><given-names>G.</given-names></name> & <name><surname>Schneider</surname><given-names>S.</given-names></name><article-title>Arlequin (v. 3.0): an integrated software package for population genetics data analysis</article-title>. <source>Evol. Bioinform. Online</source><volume>1</volume>, <fpage>47</fpage>–<lpage>50</lpage> (<year>2005</year>).<pub-id pub-id-type="pmid">19325852</pub-id></mixed-citation></ref><ref id="b68"><mixed-citation publication-type="journal"><name><surname>Drummond</surname><given-names>A. J.</given-names></name>, <name><surname>Rambaut</surname><given-names>A.</given-names></name>, <name><surname>Shapiro</surname><given-names>B.</given-names></name> & <name><surname>Pybus</surname><given-names>O. G.</given-names></name><article-title>Bayesian coalescent inference of past population dynamics from molecular sequences</article-title>. <source>Mol. Biol. Evol.</source><volume>22</volume>, <fpage>1185</fpage>–<lpage>1192</lpage> (<year>2005</year>).<pub-id pub-id-type="pmid">15703244</pub-id></mixed-citation></ref><ref id="b69"><mixed-citation publication-type="journal"><name><surname>Harpending</surname><given-names>H. C.</given-names></name><article-title>Signature of ancient population growth in a low-resolution mitochondrial DNA mismatch distribution</article-title>. <source>Hum. Biol.</source><volume>66</volume>, <fpage>591</fpage>–<lpage>600</lpage> (<year>1994</year>).<pub-id pub-id-type="pmid">8088750</pub-id></mixed-citation></ref><ref id="b70"><mixed-citation publication-type="journal"><name><surname>Rogers</surname><given-names>A. R.</given-names></name> & <name><surname>Harpending</surname><given-names>H.</given-names></name><article-title>Population growth makes waves in the distribution of pairwise genetic differences</article-title>. <source>Mol. Biol. Evol.</source><volume>9</volume>, <fpage>552</fpage>–<lpage>569</lpage> (<year>1992</year>).<pub-id pub-id-type="pmid">1316531</pub-id></mixed-citation></ref><ref id="b71"><mixed-citation publication-type="journal"><name><surname>Tajima</surname><given-names>F.</given-names></name><article-title>Statistical method for testing the neutral mutation hypothesis by DNA polymorphism</article-title>. <source>Genetics</source><volume>123</volume>, <fpage>585</fpage>–<lpage>595</lpage> (<year>1989</year>).<pub-id pub-id-type="pmid">2513255</pub-id></mixed-citation></ref><ref id="b72"><mixed-citation publication-type="journal"><name><surname>Fu</surname><given-names>Y. X.</given-names></name><article-title>Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection</article-title>. <source>Genetics</source><volume>147</volume>, <fpage>915</fpage>–<lpage>925</lpage> (<year>1997</year>).<pub-id pub-id-type="pmid">9335623</pub-id></mixed-citation></ref><ref id="b73"><mixed-citation publication-type="journal"><name><surname>Manni</surname><given-names>F.</given-names></name>, <name><surname>Guerard</surname><given-names>E.</given-names></name> & <name><surname>Heyer</surname><given-names>E.</given-names></name><article-title>Geographic patterns of (genetic, morphologic, linguistic) variation: how barriers can be detected by using Monmonier’s algorithm</article-title>. <source>Hum. Biol.</source><volume>76</volume>, <fpage>173</fpage>–<lpage>190</lpage> (<year>2004</year>).<pub-id pub-id-type="pmid">15359530</pub-id></mixed-citation></ref><ref id="b74"><mixed-citation publication-type="journal"><name><surname>Collins</surname><given-names>W. D.</given-names></name>, <name><surname>Lee-Taylor</surname><given-names>J. M.</given-names></name>, <name><surname>Edwards</surname><given-names>D. P.</given-names></name> & <name><surname>Francis</surname><given-names>G. L.</given-names></name><article-title>Effects of increased near-infrared absorption by water vapor on the climate system</article-title>. <source><italic>J. Geophys. Res</italic>.</source><volume>111</volume>, <fpage>18</fpage>–<lpage>27</lpage> (<year>2006</year>).</mixed-citation></ref><ref id="b75"><mixed-citation publication-type="journal"><name><surname>Collins</surname><given-names>W. D.</given-names></name><italic>et al.</italic><article-title>The community climate system model version 3 (CCSM3)</article-title>. <source><italic>J. Clim</italic>.</source><volume>19</volume>, <fpage>2122</fpage>–<lpage>2143</lpage> (<year>2006</year>).</mixed-citation></ref><ref id="b76"><mixed-citation publication-type="journal"><name><surname>Hasumi</surname><given-names>H.</given-names></name> & <name><surname>Emori</surname><given-names>S.</given-names></name><source>K-1 Coupled GCM (MIROC) Description</source> . (Center for Climate System Research, University of Tokyo, Tokyo, <year>2004)</year>.</mixed-citation></ref><ref id="b77"><mixed-citation publication-type="journal"><name><surname>Phillips</surname><given-names>S. J.</given-names></name> & <name><surname>Dudík</surname><given-names>M.</given-names></name><article-title>Modeling of species distributions with Maxent: new extensions and a comprehensive evaluation</article-title>. <source>Ecography</source><volume>31</volume>, <fpage>161</fpage>–<lpage>175</lpage> (<year>2008</year>).</mixed-citation></ref><ref id="b78"><mixed-citation publication-type="journal"><name><surname>Araujo</surname><given-names>M. B.</given-names></name> & <name><surname>New</surname><given-names>M.</given-names></name><article-title>Ensemble forecasting of species distributions</article-title>. <source><italic>Trends Ecol. Evol</italic>.</source><volume>22</volume>, <fpage>42</fpage>–<lpage>47</lpage> (<year>2007</year>).<pub-id pub-id-type="pmid">17011070</pub-id></mixed-citation></ref><ref id="b79"><mixed-citation publication-type="journal"><name><surname>Marmion</surname><given-names>M.</given-names></name>, <name><surname>Parviainen</surname><given-names>M.</given-names></name>, <name><surname>Luoto</surname><given-names>M.</given-names></name>, <name><surname>Heikkinen</surname><given-names>R.K.</given-names></name> & <name><surname>Thuiller</surname><given-names>W.</given-names></name><article-title>Evaluation of consensus methods in predictive species distribution modelling</article-title>. <source>Divers. Distrib.</source><volume>15</volume>, <fpage>59</fpage>–<lpage>69</lpage> (<year>2009</year>).</mixed-citation></ref><ref id="b80"><mixed-citation publication-type="journal"><name><surname>Schoener</surname><given-names>T. W.</given-names></name><article-title>The Anolis lizards of Bimini: resource partitioning in a complex fauna</article-title>. <source>Ecology.</source><volume>49</volume>, <fpage>704</fpage>–<lpage>726</lpage> (<year>1968</year>).</mixed-citation></ref><ref id="b81"><mixed-citation publication-type="journal"><name><surname>Warren</surname><given-names>D. L.</given-names></name>, <name><surname>Glor</surname><given-names>R. E.</given-names></name> & <name><surname>Turelli</surname><given-names>M.</given-names></name><article-title>Environmental niche equivalency versus conservatism: quantitative approaches to niche evolution</article-title>. <source>Evolution.</source><volume>62</volume>, <fpage>2868</fpage>–<lpage>2883</lpage> (<year>2008</year>).<pub-id pub-id-type="pmid">18752605</pub-id></mixed-citation></ref><ref id="b82"><mixed-citation publication-type="journal"><name><surname>Warren</surname><given-names>D. L.</given-names></name>, <name><surname>Glor</surname><given-names>R. E.</given-names></name> & <name><surname>Turelli</surname><given-names>M.</given-names></name><article-title>ENMTools: a toolbox for comparative studies of environmental niche models</article-title>. <source>Ecography.</source><volume>33</volume>, <fpage>607</fpage>–<lpage>611</lpage> (<year>2010</year>).</mixed-citation></ref><ref id="b83"><mixed-citation publication-type="journal"><name><surname>Warren</surname><given-names>D. L.</given-names></name> & <name><surname>Seifert</surname><given-names>S. N.</given-names></name><article-title>Ecological niche modeling in Maxent: the importance of model complexity and the performance of model selection criteria</article-title>. <source>Ecol. Appl.</source><volume>21</volume>, <fpage>335</fpage>–<lpage>342</lpage> (<year>2011</year>).<pub-id pub-id-type="pmid">21563566</pub-id></mixed-citation></ref></ref-list><fn-group><fn><p><bold>Author Contributions</bold> L.M. carried out the field work, the molecular genetic studies, statistical analysis, and drafted the manuscript. G.C. designed the study and polished the manuscript. Z.L. participated in the sample collections and the statistical analysis. Y.Y. and Z.W. designed the study. L.W. designed the study and drafted and polished the manuscript. All authors read and reviewed the final manuscript.</p></fn></fn-group></back><floats-group><fig id="f1"><label>Figure 1</label><caption><p>Sampling locations, geographical distribution of 16 matK haplotypes (H1 – H16) and three clades (<bold>A</bold>, <bold>B</bold>, <bold>C</bold>) identified by BEAST in <italic>Oxyria sinensis</italic> (population codes refer to <xref ref-type="supplementary-material" rid="S1">Supplementary Table S1</xref>). The geographic groups for the Himalaya: population 1–3, and the Hengduan Mountains, population 4–28, and east QTP for population 29–38. Map was drawn using the R package “<italic>Rgooglemaps</italic>” and ArcGIS version 9.1 and modified using CorelDRAW X6.</p></caption><graphic xlink:href="srep10396-f1"></graphic></fig><fig id="f2"><label>Figure 2</label><caption><p>The Bayesian phylogenetic tree and minimum-spanning network constructed from all <italic>mat</italic>K haplotypes of <italic>Oxyria sinensis</italic>, with <italic>O. digyna</italic> as outgroup. Values above the tree branches represent the posterior probabilities from the BEAST Bayesian analysis and those at the tree nodes denote the divergence time in million-year units. The 95% HPD for is 13.79 Mya with 95% HPD from 7.54 to 24.63 Mya. The clades A and B diverged from clade C 1.74 Mya with 95% HPD from 0.72 to 3.46 Mya, and the divergence time between clades A and B occurred 0.86 Mya with 95% HPD from 0.34 to 1.73 Mya. In the haplotype network, filled circles represent unique haplotypes and their sizes correspond to their frequencies across all populations. Each crossed bar between two circles depicts one nucleotide difference between the two haplotypes. The map was generated from the BEAST program and modified using the CorelDRAW X6 software.</p></caption><graphic xlink:href="srep10396-f2"></graphic></fig><fig id="f3"><label>Figure 3</label><caption><p>Mismatch distributions (<bold>A–D</bold>) and the Bayesian skyline plot (<bold>E</bold>) of <italic>mat</italic>K sequence data from four different regions. Except for the populations in east Himalaya, all other regional populations showed unimodal mismatch distributions, a sign of population expansion. In the Bayesian skyline plot, the estimated effective population size (<italic>N</italic>e) is depicted on the y-axis is on a <italic>log</italic>2 scale. The median estimate of <italic>N</italic>e is shown by the thick black line, and the 95% highest probability density (HPD) of posterior distribution of <italic>N</italic>e is marked in gray.</p></caption><graphic xlink:href="srep10396-f3"></graphic></fig><fig id="f4"><label>Figure 4</label><caption><p>Distribution records locations, and species potential distribution of <italic>O. sinensis</italic> under current and projected to past climate conditions. (<bold>a</bold>) The distribution records used for SDM were plotted from both this study and the CVH. We obtained the 30 arc-second resolution for altitude data from WorldClim (<ext-link ext-link-type="uri" xlink:href="http://www.worldclim.org/current">http://www.worldclim.org/current</ext-link>). This map provides the information about the locations of mountains and the occurrences of <italic>O. sinensis</italic>. The predicted species occurrences were obtained with an ensemble-forecasting approach: (<bold>b</bold>) current climate, (<bold>c</bold>) the LGM climate for CCSM, and (<bold>d</bold>) the LGM climate for MIROC. The color gradient from white to green and red indicates the suitability levels from low to high. The blue points depict the sampling location in this study, and gray points represent specimen records from the CVH.</p></caption><graphic xlink:href="srep10396-f4"></graphic></fig><fig id="f5"><label>Figure 5</label><caption><p>Glacier locations for current and the LGM. (<bold>a</bold>) The glaciers in present conditions and (<bold>b</bold>) the glaciers in LGM. Because there is no consensus about the exact locations and range of glaciers at LGM on the QTP, we reconstructed the glaciers for both LGM and the present using the Annual Mean Temperature. The sampled <italic>O. sinensis</italic> between latitude 29° and 32° were surrounded by mountain glaciers, indicating that mountain glaciers played an important role in shaping the population structure.</p></caption><graphic xlink:href="srep10396-f5"></graphic></fig><table-wrap position="float" id="t1"><label>Table 1</label><caption><title>Results of the analyses of molecular variance (AMOVA) for <italic>O. sinensis</italic>. Overall, genetic differentiation was examined for a total of 477 individuals from 38 geographic populations, as well as hierarchical geographic structures tests for three geographic groups: the Himalaya, the Hengduan Mountains and east QTP.</title></caption><table frame="hsides" rules="groups" border="1"><colgroup><col align="left"></col><col align="char" char="."></col><col align="char" char="."></col><col align="char" char="."></col><col align="char" char="."></col><col align="char" char=""></col></colgroup><thead valign="bottom"><tr><th align="left" valign="top" charoff="50"><bold>Source of variation</bold></th><th align="char" valign="top" char="." charoff="50"><bold>df</bold></th><th align="char" valign="top" char="." charoff="50"><bold>Sum of Squares</bold></th><th align="char" valign="top" char="." charoff="50"><bold>Variance Components</bold></th><th align="char" valign="top" char="." charoff="50"><bold>Variation (%)</bold></th><th align="char" valign="top" charoff="50"><bold>Fixation index</bold></th></tr></thead><tbody valign="top"><tr><td colspan="6" align="left" valign="top" charoff="50"><italic><bold>All</bold></italic></td></tr><tr><td align="left" valign="top" charoff="50">Among Populations</td><td align="char" valign="top" char="." charoff="50">37</td><td align="char" valign="top" char="." charoff="50">294.526</td><td align="char" valign="top" char="." charoff="50">0.6196 Va</td><td align="char" valign="top" char="." charoff="50">66.82</td><td align="center" valign="top" charoff="50"><italic>F</italic>st = 0.6682***</td></tr><tr><td align="left" valign="top" charoff="50">Within Populations</td><td align="char" valign="top" char="." charoff="50">439</td><td align="char" valign="top" char="." charoff="50">135.047</td><td align="char" valign="top" char="." charoff="50">0.3076 Vb</td><td align="char" valign="top" char="." charoff="50">33.18</td><td align="char" valign="top" charoff="50"> </td></tr><tr><td colspan="6" align="center" valign="top" charoff="50"><italic><bold>Himalaya (1–3) vs. Hengduan Mountains (4–28) vs. east QTP (29–38)</bold></italic></td></tr><tr><td align="left" valign="top" charoff="50">Among groups</td><td align="char" valign="top" char="." charoff="50">2</td><td align="char" valign="top" char="." charoff="50">123.615</td><td align="char" valign="top" char="." charoff="50">0.4355 Va</td><td align="char" valign="top" char="." charoff="50">38.95</td><td align="center" valign="top" charoff="50"><italic>F</italic>sc = 0.5512***</td></tr><tr><td align="left" valign="top" charoff="50">Among populations within groups</td><td align="char" valign="top" char="." charoff="50">35</td><td align="char" valign="top" char="." charoff="50">170.911</td><td align="char" valign="top" char="." charoff="50">0.3778 Vb</td><td align="char" valign="top" char="." charoff="50">33.71</td><td align="center" valign="top" charoff="50">Fst = 0.7256***</td></tr><tr><td align="left" valign="top" charoff="50">Within populations</td><td align="char" valign="top" char="." charoff="50">439</td><td align="char" valign="top" char="." charoff="50">135.047</td><td align="char" valign="top" char="." charoff="50">0.3076 Vc</td><td align="char" valign="top" char="." charoff="50">27.44</td><td align="center" valign="top" charoff="50"><italic>F</italic>ct = 0.3885***</td></tr></tbody></table><table-wrap-foot><fn id="t1-fn1"><p>All tests were based on both molecular distances and haplotype frequencies. Statistical significance: P < 0.001 (***) and P < 0.01 (**).</p></fn></table-wrap-foot></table-wrap><table-wrap position="float" id="t2"><label>Table 2</label><caption><title>Neutrality and mismatch analyses for the <italic>O. sinensis</italic> populations and four designated regional groups,</title></caption><table frame="hsides" rules="groups" border="1"><colgroup><col align="left"></col><col align="char" char="."></col><col align="char" char="."></col><col align="char" char="."></col><col align="char" char="."></col><col align="char" char="."></col><col align="char" char="."></col><col align="char" char="."></col><col align="char" char="."></col><col align="char" char="."></col></colgroup><thead valign="bottom"><tr><th align="center" valign="top" charoff="50"> </th><th align="center" valign="top" char="." charoff="50"><italic><bold>N</bold></italic></th><th align="center" valign="top" char="." charoff="50"><bold>τ</bold></th><th align="center" valign="top" char="." charoff="50"><bold>Θ</bold><sub><bold>0</bold></sub></th><th align="center" valign="top" char="." charoff="50"><bold>Θ</bold><sub><bold>1</bold></sub></th><th align="center" valign="top" char="." charoff="50"><italic><bold>SSD</bold></italic></th><th align="center" valign="top" char="." charoff="50"><italic><bold>RAG (p)</bold></italic></th><th align="center" valign="top" char="." charoff="50"><italic><bold>Tajima’s</bold></italic><bold>D</bold></th><th align="center" valign="top" char="." charoff="50"><italic><bold>Fu’s</bold></italic><bold>Fs</bold></th><th align="center" valign="top" char="." charoff="50"><italic><bold>Expansion Time (year)</bold></italic></th></tr></thead><tbody valign="top"><tr><td align="left" valign="top" charoff="50">Himalaya</td><td align="char" valign="top" char="." charoff="50">26</td><td align="char" valign="top" char="." charoff="50">7.2813</td><td align="char" valign="top" char="." charoff="50">0.0000</td><td align="char" valign="top" char="." charoff="50">3.6719</td><td align="char" valign="top" char="." charoff="50">0.0919</td><td align="char" valign="top" char="." charoff="50">0.1199 (0.2980)</td><td align="char" valign="top" char="." charoff="50">1.5950</td><td align="char" valign="top" char="." charoff="50">2.7058</td><td align="char" valign="top" char="." charoff="50">−</td></tr><tr><td align="left" valign="top" charoff="50">Hengduan Mountains</td><td align="char" valign="top" char="." charoff="50">321</td><td align="char" valign="top" char="." charoff="50">1.3457</td><td align="char" valign="top" char="." charoff="50">2.0000</td><td align="char" valign="top" char="." charoff="50">99999</td><td align="char" valign="top" char="." charoff="50">0.0018</td><td align="char" valign="top" char="." charoff="50">0.0585 (0.1480)</td><td align="char" valign="top" char="." charoff="50">−0.4614</td><td align="char" valign="top" char="." charoff="50">−0.7633</td><td align="char" valign="top" char="." charoff="50">143404</td></tr><tr><td align="left" valign="top" charoff="50">East QTP</td><td align="char" valign="top" char="." charoff="50">130</td><td align="char" valign="top" char="." charoff="50">3.0000</td><td align="char" valign="top" char="." charoff="50">0.0000</td><td align="char" valign="top" char="." charoff="50">0.0808</td><td align="char" valign="top" char="." charoff="50">0.0024</td><td align="char" valign="top" char="." charoff="50">0.7502 (0.7670)</td><td align="char" valign="top" char="." charoff="50">−1.5444*</td><td align="char" valign="top" char="." charoff="50">−5.2660**</td><td align="char" valign="top" char="." charoff="50">319693</td></tr><tr><td align="left" valign="top" charoff="50">Total</td><td align="char" valign="top" char="." charoff="50">477</td><td align="char" valign="top" char="." charoff="50">2.0703</td><td align="char" valign="top" char="." charoff="50">0.0211</td><td align="char" valign="top" char="." charoff="50">11.7047</td><td align="char" valign="top" char="." charoff="50">0.0207*</td><td align="char" valign="top" char="." charoff="50">0.0744 (0.1230)</td><td align="char" valign="top" char="." charoff="50">−0.7842</td><td align="char" valign="top" char="." charoff="50">−2.9243</td><td align="char" valign="top" char="." charoff="50">220621</td></tr></tbody></table><table-wrap-foot><fn id="t2-fn1"><p>Note: <italic>N</italic>, sample size; <italic>τ</italic>, time in generation unit since the population expansion happened; <italic>θ</italic><sub>0</sub> and <italic>θ</italic><sub><italic>1</italic></sub>, nucleotide diversity measures prior to and posterior to expansion scaled with population size; RAG, the Harpending’s Raggedness index; <italic>a</italic>, H3 was excluded for the calculation; SSD, sum of squared deviations.</p></fn></table-wrap-foot></table-wrap><table-wrap position="float" id="t3"><label>Table 3</label><caption><title>Degree of niche overlap based on Schoener’s <italic>I</italic> and <italic>D</italic> statistics. The lower triangles represent the <italic>D</italic> value and the upper represent the <italic>I</italic> value. The statistics indicate high identity among suitable habitats among all pairwise comparisons.</title></caption><table frame="hsides" rules="groups" border="1"><colgroup><col align="left"></col><col align="char" char="."></col><col align="char" char="."></col><col align="char" char="."></col></colgroup><thead valign="bottom"><tr><th align="center" valign="top" charoff="50"><italic><bold>D\I</bold></italic></th><th align="center" valign="top" char="." charoff="50"><bold>Present</bold></th><th align="center" valign="top" char="." charoff="50"><bold>LGM (CCSM)</bold></th><th align="center" valign="top" char="." charoff="50"><bold>LGM (MIROC)</bold></th></tr></thead><tbody valign="top"><tr><td align="left" valign="top" charoff="50">Present</td><td align="char" valign="top" char="." charoff="50">−</td><td align="char" valign="top" char="." charoff="50">0.9585</td><td align="char" valign="top" char="." charoff="50">0.9405</td></tr><tr><td align="left" valign="top" charoff="50">LGM (CCSM)</td><td align="char" valign="top" char="." charoff="50">0.7936</td><td align="char" valign="top" char="." charoff="50">−</td><td align="char" valign="top" char="." charoff="50">0.9627</td></tr><tr><td align="left" valign="top" charoff="50">LGM (MIROC)</td><td align="char" valign="top" char="." charoff="50">0.7447</td><td align="char" valign="top" char="." charoff="50">0.8063</td><td align="char" valign="top" char="." charoff="50">−</td></tr></tbody></table></table-wrap></floats-group></pmc></record>Pour manipuler ce document sous Unix (Dilib)
EXPLOR_STEP=$WICRI_ROOT/Wicri/Linguistique/explor/TamazightV2/Data/Pmc/Corpus
HfdSelect -h $EXPLOR_STEP/biblio.hfd -nk 000045 | SxmlIndent | more
Ou
HfdSelect -h $EXPLOR_AREA/Data/Pmc/Corpus/biblio.hfd -nk 000045 | SxmlIndent | more
Pour mettre un lien sur cette page dans le réseau Wicri
{{Explor lien
|wiki= Wicri/Linguistique
|area= TamazightV2
|flux= Pmc
|étape= Corpus
|type= RBID
|clé=
|texte=
}}
| This area was generated with Dilib version V0.6.33. |  |