Abstract: Species delimitation is a key foundation for exploring biodiversity. However, the existence of continuous phenotypic variation in widespread species challenges accurate species delimitation based on classical taxonomy. In this study, we investigated the cryptic diversity of a widespread herb (Roscoea tibetica Batalin) in a biodiversity hotspot (the Hengduan Mountains, China) using genotyping by sequencing, examining morphological traits, developing species distribution models, and simulating demographic history. Phylogenomic reconstruction, principal component analysis, and genetic structure inferences indicated that previously reported R. tibetica comprised two monophyletic lineages with a deep divergence. Several morphological diagnostic characteristics were discovered from field and common garden that corresponded to these independent evolutionary lineages. Species distribution models illustrated significant ecological divergence between both lineages. All evidence strongly supported that R. tibetica, as described in previous taxonomy, actually comprises two distinct species. Model test of gene flow and effective population size changes in fastsimcoal2, and a negative Tajima's D-value suggested that recent contact likely occurred between the two lineages. Our results proposed that cryptic diversity in previously reported R. tibetica was possibly associated with phenotypic plasticity in heterogeneous environments and morphological convergence in similar habitats. This study suggests that caution should be exercised when attempting to gain biological insight into species with large-scale morphological variation, and species delimitation should be done in advance.

Key words: cryptic species, genotyping by sequencing, phylogenomics, Roscoea tibetica, species delimitation