Gaz-angubin, the manna used in the production of famous sweet GAZ of Isfahan (or Khunsar) for a long time, has had special economic and medicinal value in Iran. In creating this manna ( Cyamophila astragalicola ; Gegechkori 1977), recent years, due to the significant decline in the population of the insect production and harvest of this valuable manna have decreased. Climate change, as a threatening factor of biodiversity, has greatly affected species distribution. Further, Genetic diversity in species is required for evolutionary potential and adaptation to and dispersal ability by studying the gene flow between populations. The environmental changes that we can get an insight into the evolutionary history current study aimed to model the current and future potential distribution of C. astragalicola under the climate change scenario (CCSM4), and to evaluate genetic variation of the species, using the mitochondrial CO1 sequence. Phylogenetic relationships were reconstructed using IQ-TREE software, using 27 sequences of C .astragalicola and sequences of other psyllids from superfamily Psylloidea. Distribution modelling was performed, using MaxEnt model, to select the most important variable based on the Jackkife test. Graphical-statistical environment of R software and Biomod2 package were, then, used to create an ensemble model. Among the applicable models, GBM, GLM, GAM, and MaxEnt were changes in suitability under the climate change scenario. The ensemble model used in the ensemble model to predict current habitat suitability and possible revealed a higher performance (AUC=0.97) than single models. The important (NDSI) respectively. Overall, climate change modeling indicates the sensitivity variables in habitat selection include Annual Mean Temperature (BIO1), Precipitation of Driest Quarter (BIO17), and Normalized Difference Snow Index of C .astragalicola to climate change, as such that, in all Representative Concentration Pathways (e.g. RCP 2.6 and RCP 8.5) for the year 2050 and 2070, loss in suitable habitats was evident. In addition, results revealed a negative effect of C .astragalico, as such that, by increasing BIO1 in each scenario, of climate change, especially an increase in temperature on habitat suitability the predicted amount of suitable habitats was reduced. Results of genetic variation in mitochondrial CO1 sequence confirmed the identification of C .astragalicola by the morphological method and all samples collected from populations in Isfahan and Chaharmahal va Bakhtiari Provinces belonged to one species. Three haplotypes were observed in this study that all, but two samples (Owned sample by Khunsar and Koohrang), belonged to the same haplotype. Key words: Climate change, Genetic diversity, Gaz-angubin, Habitat suitability, COI