Species often exhibit different levels of genetic structuring correlated to their environment. However, understanding how environmental heterogeneity influences genetic variation is difficult because the effects of gene flow, drift and selection are confounded. We investigated the genetic variation and its ecological correlates in an endemic and critically endangered stream breeding mountain newt, Neurergus kaiseri, within its entire range in southwestern Iran. We identified two geographic regions based on phylogenetic relationships using Bayesian inference and maximum likelihood of 779 bp mtDNA (D-loop) in 111 individuals from ten of twelve known breeding populations. This analysis revealed a clear divergence between northern populations, located in more humid habitats at higher elevation, and southern populations, from drier habitats at lower elevations regions. From seven haplotypes found in these populations none was shared between the two regions. Analysis of molecular variance (AMOVA) of N. kaiseri indicates that 94.03% of sequence variation is distributed among newt populations and 5.97% within them. Moreover, a high degree of genetic subdivision, mainly attributable to the existence of significant variance among the two regions is shown (θ_CT = 0.94, P = 0.002). The positive and significant correlation between geographic and genetic distances (r = 0.61, P = 0.002) following controlling for environmental distance suggests an important influence of geographic divergence of the sites in shaping the genetic variation and may provide tools for a possible conservation based prioritization policy for the endangered species.

The taxonomic status of the greater mouse-tailed bat (Rhinopoma microphyllum) in Iran is not clear and researchers have reported conflicting results. The initial suggestion of the presence of two subspecies R. m. microphyllum and R. m. harrisoni has been questioned on the basis of small differences between the populations. These differing inferences are based on analysis of morphological characteristics. Here we present a study of the phylogeography of this species using 567 bp of the mitochondrial control region in several localities in Iran and from across its distribution to infer the taxonomic status of this species. Based on the control region sequences, we found high genetic diversity in the Iranian population but variation between populations was not statistically significant. The phylogenetic trees and statistical parsimony network showed all Iranian samples were grouped in the same clade and Levant samples belonged to another clade. These results support the hypothesis that all Iranian specimens belong to one subspecies. We therefore recommend that R. m. microphyllum and R. m. harrisoni can be synonymized as the same subspecies with the name R. m. harrisoni, because molecular results indicate that Iranian samples differ from Levant and Moroccan samples (R. m. microphyllum).

To examine the level of genetic differentiation in the sequences of the mitochondrial D-loop gene of Rhinolophus ferrumequinum, and to evaluate the current taxonomic status of this species, 50 tissue samples of greater horseshoe bats were collected in 2011-2015 from 21 different localities in northwest, northeast, west, central, and south regions of Iran. Twenty-two published D-loop sequences from Europe (Switzerland, United Kingdom, Bulgaria, and Tunisia), and Anatolia (south, west, and east Turkey) were downloaded from GenBank. Molecular genetic analyses revealed remarkable variation among populations of R. ferrumequinum. Two major clades with strong support were identified within the greater horseshoe bat. One of these clades consists of individuals of R. ferrumequinum from Iran and eastern Turkey, and is further subdivided into two subclades. A second clade includes samples from western Turkey and Europe. The two subclades from Iran and Turkey and the second clade from western Turkey and Europe represent three diagnosable categories, which most probably warrant three subspecies for the species. Thus, based on genetic differences, it is clear that two subspecific populations are found in Iran: R. f. irani (southern Iran) and R. f. proximus (northern Iran).

New species of the genus Niphargus is described and named as N. kermanshahi sp. nov. from a karst spring in west of Iran. This species is identified based on the analysis of morphological characters and 28S ribosomal DNA sequences. Taxonomic status and phylogenetic position of this species is discussed in comparison to other Iranian species of Niphargus genus.

Molecular phylogeography and species distribution modelling (SDM) suggest that late Quaternary glacial cycles have portrayed a significant role in structuring current population genetic structure and diversity. Based on phylogenetic relationships using Bayesian inference and maximum likelihood of 535 bp mtDNA (D-loop) and 745 bp mtDNA (Cytb) in 62 individuals of the Mediterranean Horseshoe Bat, Rhinolophus euryale, from 13 different localities in Iran we identified two subspecific populations with differing population genetic structure distributed in southern Zagros Mts. and northern Elburz Mts. Analysis of molecular variance (AMOVA) obtained from D-loop sequences indicates that 21.18% of sequence variation is distributed among populations and 10.84% within them. Moreover, a degree of genetic subdivision, mainly attributable to the existence of significant variance among the two regions is shown (θCT = 0.68, p = .005). The positive and significant correlation between geographic and genetic distances (R² = 0.28, r = 0.529, p = .000) is obtained following controlling for environmental distance. Spatial distribution of haplotypes indicates that marginal population of the species in southern part of the species range have occupied this section as a glacial refugia. However, this genetic variation, in conjunction with results of the SDM shows a massive postglacial range expansion for R. euryale towards higher latitudes in Iran.