AUTHOR=Liu Detuan , Zhang Lu , Wang Jihua , Ma Yongpeng 

TITLE=Conservation Genomics of a Threatened Rhododendron: Contrasting Patterns of Population Structure Revealed From Neutral and Selected SNPs

JOURNAL=Frontiers in Genetics

VOLUME=11

YEAR=2020

URL=https://www.frontiersin.org/journals/genetics/articles/10.3389/fgene.2020.00757

DOI=10.3389/fgene.2020.00757

ISSN=1664-8021

ABSTRACT=<p>Though it is well-acknowledged that next generation sequencing (NGS) technologies can provide further insights into plant conservation management than traditional molecular markers, studies employing NGS to address conservation genomics and subsequent conservation strategies for threatened plants are still rare. <italic>Rhododendron</italic> is the largest genus of woody plants in China, and many species are threatened, however, to date there has been no conservation genetic research using NGS in this genus. In the present study, we investigated the conservation genetics of <italic>R. cyanocarpum</italic>, a threatened species endemic to the Cangshan Mountains in Yunnan, China, using a double digest restriction-site-associated DNA-sequencing (ddRAD-seq) approach. Due to the availability of sufficient SNPs, we were able to distinguish between neutral and putatively selected SNPs and were able to further investigate the genetic diversity, population structure, and differentiation in <italic>R. cyanocarpum</italic>, as well as make an estimation of its demographic history. A total of 6,584 SNPs were obtained, of which 5,729 were neutral (detected using Tajima’s D). In terms of the 5,729 neutral SNPs, <italic>R. cyanocarpum</italic> had a higher genetic diversity (π = 0.0702 ± 0.0017, <italic>H<sub><italic>e</italic></sub></italic> = 0.0675 ± 0.0016) than other plant species assessed using Rad-seq methods, while population differentiation (<italic>F</italic><sub><italic>st</italic></sub> from 0.0314 to 0.0452) was weak. Interestingly, contrasting patterns of population structure were revealed from all neutral and selected SNPs, with distinct genetic clusters forming for all SNPs and neutral SNPs, but no distinct subgroups for selected ones. Moreover, we were able to detect changes in effective population size (<italic>N</italic><sub><italic>e</italic></sub>) of <italic>R. cyanocarpum</italic> from 150,000 years ago, including a bottleneck event ca. 60,000 years ago, followed by recovery of <italic>N</italic><sub><italic>e</italic></sub> over a short period, and a subsequent gradual decline in <italic>N</italic><sub><italic>e</italic></sub> to date. Implications for conserving <italic>R. cyanocarpum</italic> based on these main results are then discussed.</p>