RADseq-based population genomic analysis and environmental adaptation of rare and endangered recretohalophyte Reaumuria trigyna

Genetic diversity reflects the survival potential, history, and population dynamics of an organism. It underlies the adaptive potential of populations and their response to environmental change. Reaumuria trigyna is an endemic species in the Eastern Alxa and West Ordos desert regions in China. The species has been considered a good candidate to explore the unique survival strategies of plants that inhabit this area. In this study, we performed population genomic analyses based on restriction-site associated DNA sequencing to understand the genetic diversity, population genetic structure, and differentiation of the species. Analyses of 92,719 high-quality single-nucleotide polymorphisms (SNPs) indicated that overall genetic diversity of R. trigyna was low (H-O = 0.249 and H-E = 0.208). No significant genetic differentiation was observed among the investigated populations. However, a subtle population genetic structure was detected. We suggest that this might be explained by adaptive diversification reinforced by the geographical isolation of populations. Overall, 3513 outlier SNPs were located in 243 gene-coding sequences in the R. trigyna transcriptome. Potential sites under diversifying selection occurred in genes (e.g., AP2/EREBP, E3 ubiquitin-protein ligase, FLS, and 4CL) related to phytohormone regulation and synthesis of secondary metabolites which have roles in adaptation of species. Our genetic analyses provide scientific criteria for evaluating the evolutionary capacity of R. trigyna and the discovery of unique adaptions. Our findings extend knowledge of refugia, environmental adaption, and evolution of germplasm resources that survive in the Ordos area.