Evolutionary history of the Persian Jird,Meriones persicus,based on genetics, species distribution modelling and morphometric data

Biogeography of rodents in Iran: species richness, elevational distribution and their environmental correlates

Rodent biogeographic studies are disproportionately scarce in Iran, however, they are an ideal system to understand drivers of biodiversity distributions in the country. The aims of the present research are to determine (i) the pattern of rodent richness across the country, (ii) quantify their elevational distribution patterns, and (iii) explore the underlying mechanisms. To reach these goals, an updated species list was compiled based on the latest taxonomic revisions, published until December 2021. We mapped all 76 rodent species distributions to develop the first map of rodent richness in Iran. We furthermore investigated their elevational distribution patterns in the following four geographic regions based on 100 m intervals; the Zagros Mountains, north of Alborz Mountains, south of Alborz and Kopet-Dagh mountains, and central and east mountains. North-east of Iran, Zagros Mountains, Alborz Mountains, and northwestern of Iran showed the highest richness and were identified as biodiversity hotspots of rodents in the country. This study highlights the importance of past climate change as the key driver of rodent richness in Iran. We showed that rodents’ elevational distribution patterns differ among geographic regions. Areas and elevational zones with the highest species richness should be prioritized for the conservation planning of rodents in Iran.

Pleistocene climate fluctuations as the major driver of genetic diversity and distribution patterns of the Caspian green lizard, Lacerta strigata Eichwald, 1831

Green lizards of the genus Lacerta have served as excellent models for studying the impact of Pleistocene climatic oscillations on genetic structures. The Caspian green lizard, Lacerta strigata, occupies various habitats across the Caucasus and the South Caspian Sea, with the Hyrcanian Forests and north of the Alborz Mountains forming the core of the range. This study aimed to re-examine the phylogenetic relationships of L. strigata with other congeneric members and to assess the genetic structure and historical demography of the species. Furthermore, Species Distribution Models (SDMs) were performed to infer the species' potential habitat suitability and were then projected on climate scenarios reflecting current and past (6 ky and 21 ky before present) conditions. A total of 39 individuals collected from most of the distribution range, together with additional lacertid species sequences from the GenBank database, were examined using mtDNA (Cyt b and 12S ribosomal RNA) and nuclear (C-mos and β-fibrinogen) sequence data. Based on the phylogenetic analyses, L. strigata was found to be a sister taxon to all other members of the genus. The species included two main clades (regional western and eastern) that diverged in a period between the Early and Middle Pleistocene. Based on the BBM and S-Diva analyses, both dispersal and vicariance events explained the phylogeographic structure of the species in the Hyrcanian Forests. The historical demographic analyses using Bayesian skyline plots showed a mild increase in the effective population size from about 120 Kya for the western regional clade. According to phylogeographic structures and SDMs evidence, as in other species within the region, it appears that the south of the Caspian Sea (Hyrcanian Forests), and the Alborz Mountains acted as multiple refugia during cold periods and promoted expansion outwards amid the warm periods. Overall, the results provided evidence that the genetic structure of the species has been influenced by the Pleistocene climatic fluctuations.

Predicted Future Benefits for an Endemic Rodent in the Irano-Turanian Region

Climate change is expected to have an impact on the geographical distribution ranges of species. Endemic species and those with a restricted geographic range may be especially vulnerable. The Persian jird (Meriones persicus) is an endemic rodent inhabiting the mountainous areas of the Irano-Turanian region, where future desertification may form a threat to the species. In this study, the species distribution modelling algorithm MaxEnt was used to assess the impact of future climate change on the geographic distribution range of the Persian jird. Predictions were made under two Representative Concentration Pathways and five different climate models for the years 2050 and 2070. It was found that both bioclimatic variables and land use variables were important in determining potential suitability of the region for the species to occur. In most cases, the future predictions showed an expansion of the geographic range of the Persian jird which indicates that the species is not under immediate threat. There are however uncertainties with regards to its current range. Predictions may therefore be an over or underestimation of the total suitable area. Further research is thus needed to confirm the current geographic range of the Persian jird to be able to improve assessments of the impact of future climate change.

The Zagros Mountains acting as a natural barrier to gene flow in the Middle East: more evidence from the evolutionary history of spiny-tailed lizards (Uromasticinae: Saara)

We investigated the phylogeographic structure of the genus Saara and studied whether the Zagros Mountain ranges influence the gene flow and dispersal at the landscape scale between the Iranian Saara species, including S. asmussi in the central Iranian plateau and S. loricata in the Mesopotamian lowlands. Phylogenetic analyses clearly show three well-supported species, including S. loricata, S. asmussi and S. hardwickii, that are distinct from Uromastyx species. The S-DIVA and BBM analyses demonstrate that species of Saara originated from an ancestor somewhere in the Iranian Plateau and then dispersed to its current geographical range. Our results indicate that the separation of S. loricata from S. asmussi coincides with the orogenic events of the Zagros Mountains during the Pliocene, and thus it may show a vicariance event. Diversification within populations of S. loricata and S. assmussi are estimated to have occurred during the Pleistocene. The haplotype network indicates one haplogroup for each of the Iranian Saara species. Population genetic analyses shows signals of demographic expansions at the beginning of the Holocene for S. loricata. Our results support the hypothesis that the Zagros Mountains act as a barrier for gene flow and Quaternary climatic oscillations affected intraspecific genetic divergences of S. loricata and S. asmussi.

Reptile species richness associated to ecological and historical variables in Iran

Spatial gradients of species richness can be shaped by the interplay between historical and ecological factors. They might interact in particularly complex ways in heterogeneous mountainous landscapes with strong climatic and geological contrasts. We mapped the distribution of 171 lizard species to investigate species richness patterns for all species (171), diurnal species (101), and nocturnal species (70) separately. We related species richness with the historical (past climate change, mountain uplifting) and ecological variables (climate, topography and vegetation). We found that assemblages in the Western Zagros Mountains, north eastern and north western parts of Central Iranian Plateau have the highest number of lizard species. Among the investigated variables, annual mean temperature explained the largest variance for all species (10%) and nocturnal species (31%). For diurnal species, temperature change velocity shows strongest explained variance in observed richness pattern (26%). Together, our results reveal that areas with annual temperature of 15–20 °C, which receive 400–600 mm precipitation and experienced moderate level of climate change since the Last Glacial Maximum (LGM) have highest number of species. Documented patterns of our study provide a baseline for understanding the potential effect of ongoing climate change on lizard diversity in Iran.

Evolutionary history of water voles revisited: confronting a new phylogenetic model from molecular data with the fossil record

Recent water voles (genus Arvicola) display a prominent morphological diversity with a strong ecotypical background but with unclear taxonomic associations. We provide a novel synthetic view on the evolutionary history and the current taxonomic richness in the genus. Our molecular reconstruction, based on a 1143-bp-long sequence of cytochrome b and a 926-bp interphotoreceptor retinoid binding protein (irbp) confirmed the monophyly of four species (amphibius, sapidus, monticola and italicus) recognized thus far, and retrieved a new deeply divergent lineage from West Iran. Genetic divergence of the Iranian lineage (>9.0%) is inside the range of interspecies distances, exceeding the interspecies divergences between the remaining Arvicola species (range, 4.3–8.7%). The oldest name available for the Iranian phylogroup is Arvicola persicus de Filippi, 1865, with the type locality in Soltaniyeh, Iran. The molecular clock suggests the divergence of A. persicus in the Early Pleistocene (2.545 Ma), and the current radiation of the remaining species between 1.535 Ma (Arvicola sapidus) and 0.671 Ma. While A. sapidus possibly evolved from Arvicola jacobaeus, a fossil ancestor to A. persicus is unknown. The aquatic life-style of Mimomys savini, a direct ancestor to some fossil Arvicola, is retained in recent stem species A. sapidus and A. persicus, while a major shift toward fossorial morphotype characterizes the terminal lineages (amphibius, italicus and monticola). We suggest that habitat-dependent morphological plasticity and positive enamel differentiation in Arvicola amphibius widened its ecological niche that might trigger a range expansion across c. 12 million km2, making it one of the largest among arvicolines.

Population genetic structure and phylogeography of the greater horseshoe bat (Rhinolophus ferrumequinum) along Alborz and Zagros Mts. in Iran

In this paper, we investigate the genetic structure and phylogeography of Rhinolophus ferrumequinum, using the mitochondrial cytochrome b gene (1017 bp) in Iran and adjacent regions. The total haplotype and nucleotide diversity are 0.63 ± 0.055 and 0.0021 ± 0.00017, respectively which suggest that R. ferrumequinum exhibits low genetic diversity. AMOVA analysis shows that more variation of genetic differentiation is present among populations of phylogenetic groupings than within populations. Our phylogenetic results support the monophyly of R. ferrumequinum and suggest this taxon comprises three allopatric/parapatric phylogroups that are distributed in Europe-western Turkey, eastern Turkey-northern Iran, and southern Iran. The Europe-western Turkey lineage (clade 2) split from the eastern Turkey-Iran lineage (clade 1) during the middle Pleistocene (0.8534 (ca.I)-0.6454 (ca.II) Ma). The divergence time among subclades A and B occurred during the mid-Pleistocene (0.4849 (ca.I)-0.369 (ca.II) Ma). All phylogenetic analyses also indicate that the Iranian and eastern Turkey R. ferrumequinum diverged from Europe and western Turkey R. ferrumequinum, with the mean percentage sequence differences ranging from 0.92%-0.75% between them. We infer that long-term isolation of R. ferrumequinum in spatially distinct refugia in parts of southwestern and northeastern Iran has promoted distinct phylogeographic lineages during the Pleistocene.

Mitochondrial DNA sequence analysis reveals multiple Pleistocene glacial refugia for the Yellow-spotted mountain newt, Neurergus derjugini (Caudata: Salamandridae) in the mid-Zagros range in Iran and Iraq

Phylogeography is often used to investigate the effects of glacial cycles on current genetic structure of various plant and animal species. This approach can also identify the number and location of glacial refugia as well as the recolonization routes from those refugia to the current locations. To identify the location of glacial refugia of the Yellow‐spotted mountain newt, Neurergus derjugini, we employed phylogeography patterns and genetic variability of this species by analyzing partial ND4 sequences (867 bp) of 67 specimens from 15 sampling localities from the whole species range in Iran and Iraq. Phylogenetic trees concordant with haplotype networks showed a clear genetic structure among populations as three groups corresponding to the populations in the north, center, and south. Evolutionary ages of clades north and south ranging from 0.15 to 0.17 Myr, while the oldest clade is the central clade, corresponding to 0.32 Myr. Bayesian skyline plots of population size change through time show a relatively slight increase until about 25 kyr (around the last glacial maximum) and a decline of population size about 2.5 kyr. The presence of geographically structured clades in north, center, and south sections of the species range signifies the disjunct populations that have emerged in three different refugium. This study illustrates the importance of the effect of previous glacial cycles in shaping the genetic structure of mountain species in the Zagros range. These areas are important in terms of long‐term species persistence and therefore valuable areas for conservation of biodiversity.

Hantavirus infection in Iranian patients suspected to viral hemorrhagic fever

BACKGROUND

Hantaviruses are a group of emerging pathogens causing hemorrhagic fever with renal syndrome and Hantavirus cardiopulmonary syndrome in human. This study was conducted to investigate Hantavirus infection among Iranian viral hemorrhagic fever suspected patients.

METHODS

From April 2014 to June 2016, 113 cases from 25 different provinces of Iran were analyzed for Hantavirus infection by IgM/IgG ELISA and pan-Hantavirus RT-PCR tests.

RESULTS

Although, viral genome was detected in none of the subjects, IgM and IgG antibodies were detected in 19 and 4 cases, respectively. Differentiation of the anti-Hantavirus antibodies according to virus species by EUROLINE Anti-Hantavirus Profile Kit revealed three Puumala virus IgM positive, one Hantaan virus IgM positive, one Hantaan virus IgM borderline, and two Puumala virus IgG borderline cases.

CONCLUSIONS

This study demonstrates the circulation of Hantaviruses in Iran and calls for further investigations of these life-threatening viruses in the country.

Historical explanation of genetic variation in the Mediterranean horseshoe bat Rhinolophus euryale (Chiroptera: Rhinolophidae) inferred from mitochondrial cytochrome-b and D-loop genes in Iran

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.