Southwestern Populations Research Articles

The identity of Tokays (Squamata: Gekkonidae: Gekko) from Hong Kong.

The Tokay gecko (Gekko gecko), found throughout Southeast Asia and India, is a heavily harvested species of high commercial value. Recent studies have supported the elevation of the black morph of Tokay gecko, found only in southern China and northern Vietnam, to the species Gekko reevesii. Previous genetic studies focused on specimens of G. reevesii from southwestern populations. Hong Kong, in southeastern China, has native populations of G. reevesii. To verify the identity of G. reevesii in Hong Kong, we employed three mitochondrial genes (COI, Cytb & ND2) and constructed a matrilineal genealogy using other specimens from Guangxi (southwestern China) and northern Vietnam, as well as G. gecko from a wide range of Southeast Asian countries. Our study confirmed that G. reevesii occurs naturally in Hong Kong, but one exotic population of G. gecko, likely a translocation from international trade, was also revealed. Our study did not reject the species ranking of G. reevesii. Moreover, like previous studies, we recovered a paraphyletic G. gecko, which may reflect a species complex, hybridization or incomplete lineage sorting. More extensive sampling of the two species over a broader range of their asserted distribution together with the use of both mitochondrial and nuclear DNA are required to better investigate their biogeography.

Evolutionary history of Nasutitermes kemneri (Termitidae, Nasutitermitinae), a termite from the South American diagonal of open formations

Little is known about the phylogeography of termites in the Neotropical region. Here, we explored the genetic patterns and phylogeographical processes in the evolutionary history of Nasutitermes kemneri, an endemic termite of the South American diagonal of open formations (DOF) formed by the Chaco, Cerrado, and Caatinga phytogeographic domains. We sampled 60 individuals across the three domains of the DOF, and using the mitochondrial genes 16S, COI, and COII, as well as the nuclear gene ITS, evaluated the genetic diversity and divergence time of the populations, along with their genetic structure. The results show a strong genetic and spatial structure within the samples, evidencing the existence of two well-differentiated genetic groups: the Northeastern and the Southwestern populations, which diverged about 2.5 Mya, during the Pliocene-Pleistocene boundary. The Northeastern population, which encompasses Caatinga and northern portions of Cerrado, has an intricate structure and seems to have suffered repetitive retraction-expansion events due to climactic fluctuations during the Quaternary. The Southwestern population, which ranges from central-south Cerrado to the northeast peripherical portions of the Chaco, displays a star-shaped haplotype structure, indicating that this region may have acted as a refugia during interglacial periods.

Isotopic signature in isolated south-western populations of European brown bear (Ursus arctos)

Stable isotope analysis of animal tissue samples is increasingly used to study the trophic ecology of target species. The isotopic signatures respond to the type of diet, but also to the environmental conditions of their habitat. In the case of omnivorous, seasonal or opportunistic feeding species, the interpretation of isotopic values is more complex, as it is largely determined by food selection, either due to individual choice or because of availability. We analysed C and N isotopes in brown bear (Ursus arctos) hair from four isolated populations of south-western Europe (Cantabrian, Pyrenees, Central Apennines and Alpine) accounting for the geographical and climatic differences among the four areas. We found inter-population differences in isotopic signatures that cannot be attributed to climatic differences alone, indicating that at least some bears from relatively higher altitude populations experiencing higher precipitation (Pyrenees) show a greater consumption of animal foods than those from lower altitudes (Cantabrian and Apennines). The quantification of isotopic niche space using Layman’s metrics identified significant similarities between the Cantabrian and Central Apennine samples that markedly differ from the Pyrenean and Alpine. Our study provides a baseline to allow further comparisons in isotopic niche spaces in a broad ranged omnivorous mammal, whose European distribution requires further conservation attention especially for southern isolated populations.

Calving rates, long-range movements and site fidelity of southern right whales (Eubalaena australis) in south-eastern Australia

Within New Zealand and eastern Australia, over 58,000 southern right whales were harvested by commercial whalers between 1790 and 1980, with approximately 19,000 harvested from south-eastern Australia. Local extirpation is believed to have led to a loss of cultural memory of calving areas, contributing to the limited recovery of the south-eastern Australian population. While the number of whales observed using the south-eastern Australian coastline is increasing, there has been no change over three decades in the annual abundance of cow-calf pairs at Logans Beach in Warrnambool, Victoria, the only established calving ground in the region. Knowledge of life history parameters of the south-eastern Australian population is lacking. Here, we examine sightings and photo-identification data from southern Australia to investigate calving intervals, long range movements and fidelity to the Logans Beach nursery ground. Sightings data revealed at least 93 calves were born at Logans Beach between 1980 and 2018 (an average of 2.6 per year) with a mean calving interval of 3.5 ± 1.0 years (± SE, n = 34). Comparison between photo-identification catalogues compiled for south-eastern and south-western Australian populations shows that southern right whales are wide ranging within southern Australian waters. Females can be sighted at locations as far apart as 3,800 km across seasons and there is overlap in the wintering range of the south-eastern and south-western populations, with at least 7% of whales using both regions. We also provide the first report of an Australian southern right whale female with strong site fidelity to a calving area in one region relocating long-term to a calving area in another region. This work highlights several knowledge gaps, such as; the location of feeding and conception grounds for this population as well as the degree of mixing between the two Australian populations outside their wintering areas. In addition, the proportion of female calves born at Logans Beach returning to their natal site to calve remains unclear. Our work provides the first assessment of calving rates, movement and site fidelity within the south-eastern Australian population, critical for understanding constraints to recovery and informing conservation management of southern right whales in Australia. Targeted, long-term monitoring programs across the south-eastern Australian region are needed to provide demographic information on which to base predictions of the impacts of anthropogenic threats such as noise disturbance, entanglement and vessel strike.

Truss network based morphometric and meristic variations among south-western populations of Macrognathus aculeatus in Bangladesh

Aim: Landmark-based truss network system was applied to examine the morphometric and meristic differences among four different populaces of Macrognathus aculeatus namely; Dhakuriabeel, Jashore (DBJ); Arial Kha River, Madaripur (AKRM); Nabaganga River, Jhenaidah (NRJ); Kopotakkha River, Jashore (KRJ) in Bangladesh. Methodology: Completely 80 samples of M. aculeatus were collected (20 individuals from each population) for this study. Six meristic characters, 6 morphometric characters and 15 truss-based morphometric characters were measured for one-way analysis of variance (ANOVA) and Tukey-HSD post hoc test. Discriminant function analyses were also performed by using morphometric and truss-based morphometric characters to detect potential characters responsible for population separation. Results: Completely 80 samples of M. aculeatus were collected (20 individuals from each population) for this study. Six meristic characters, 6 morphometric characters and 15 truss-based morphometric characters were measured for one-way analysis of variance (ANOVA) and Tukey-HSD post hoc test. Discriminant function analyses were also performed by using morphometric and truss-based morphometric characters to detect potential characters responsible for population separation. Interpretation: The baseline information derived from the present study would be useful for appropriate conservation of M. aculeatus populations as well as for the researchers for further investigation.

Genotyping-by-sequencing and ecological niche modeling illuminate phylogeography, admixture, and Pleistocene range dynamics in quaking aspen (Populus tremuloides).

Populus tremuloides is the widest‐ranging tree species in North America and an ecologically important component of mesic forest ecosystems displaced by the Pleistocene glaciations. Using phylogeographic analyses of genome‐wide SNPs (34,796 SNPs, 183 individuals) and ecological niche modeling, we inferred population structure, ploidy levels, admixture, and Pleistocene range dynamics of P. tremuloides, and tested several historical biogeographical hypotheses. We found three genetic lineages located mainly in coastal–Cascades (cluster 1), east‐slope Cascades–Sierra Nevadas–Northern Rockies (cluster 2), and U.S. Rocky Mountains through southern Canadian (cluster 3) regions of the P. tremuloides range, with tree graph relationships of the form ((cluster 1, cluster 2), cluster 3). Populations consisted mainly of diploids (86%) but also small numbers of triploids (12%) and tetraploids (1%), and ploidy did not adversely affect our genetic inferences. The main vector of admixture was from cluster 3 into cluster 2, with the admixture zone trending northwest through the Rocky Mountains along a recognized phenotypic cline (Utah to Idaho). Clusters 1 and 2 provided strong support for the “stable‐edge hypothesis” that unglaciated southwestern populations persisted in situ since the last glaciation. By contrast, despite a lack of clinal genetic variation, cluster 3 exhibited “trailing‐edge” dynamics from niche suitability predictions signifying complete northward postglacial expansion. Results were also consistent with the “inland dispersal hypothesis” predicting postglacial assembly of Pacific Northwestern forest ecosystems, but rejected the hypothesis that Pacific‐coastal populations were colonized during outburst flooding from glacial Lake Missoula. Overall, congruent patterns between our phylogeographic and ecological niche modeling results and fossil pollen data demonstrate complex mixtures of stable‐edge, refugial locations, and postglacial expansion within P. tremuloides. These findings confirm and refine previous genetic studies, while strongly supporting a distinct Pacific‐coastal genetic lineage of quaking aspen.

Phylogeography of Triplophysa stenura (Nemacheilidae): Responded to the Mid-Pleistocene Climate Transition in the Qinghai-Tibetan Plateau.

Pleistocene climatic oscillations have greatly influenced the evolutionary histories and distribution patterns of most extant species. In this study, the phylogeographic patterns and evolutionary history of Triplophysa stenura were investigated. A total of 359 individuals from 19 locations covering the species' distribution range were collected, and two mitochondrial genes (COI and Cyt b) and the recombination activating protein 1 gene (Rag1) were analyzed. Two deeply divergent lineages, corresponding geographically to the northern and southern Tanggula Mountain, were observed, and shallow divergences were found within the southern and central Qinghai-Tibetan plateau (QTP). The estimated divergence time between the northern and southern Tanggula lineages was around 1.07 Mya. Within the southern Tanggula Mountain lineages, the Nu and Subansiri rivers populations were split about 0.74 Mya, and the southwestern and central QTP populations were divided with the southeastern QTP ones approximately 0.67 Mya. The divergence times of the lineages were matched with two major glaciations in QTP (the Xixiabangma Glaciation, 1.17-0.80 Mya and Kunlun Glaciation, 0.72-0.62 Mya). Together with demographic history analyses, our results highlighted that major glacial cycles during the mid-Pleistocene played a major role in sculpting the distribution pattern of T. stenura, and led to the gene homogenization crossing the drainage populations in the southwestern and central QTP.

Isotopic evidence of strong reliance on animal foods and dietary heterogeneity among Early-Middle Neolithic communities of Iberia

Stable carbon and nitrogen isotope research on past populations in the Iberian Neolithic has emphasized the Atlantic and Mediterranean coasts. This study provides the first isotopic insights into the diet and subsistence economy of Early and Middle Neolithic populations from open-air sites in interior north-central Iberia. We present bone collagen carbon (δ13C) and nitrogen (δ15N) isotope ratios for 44 humans and 33 animals recovered from six cemeteries of the Ebro valley and the northern Iberian Plateau. The results obtained are consistent with the C3 terrestrial diets typical of other contemporary south-western European populations, but the spacing between human and herbivore values from Los Cascajos and Paternanbidea sites is higher than expected, and a significant positive correlation is identified between the δ13C and δ15N human values at both. Moreover, the results clearly differ from those of the Late Neolithic/Early Chalcolithic in the same region, which show significantly lower δ13C and δ15N values. These findings contribute to an understanding of the implementation of an agro-pastoral economy in interior Iberia, suggesting a stronger reliance on animal foods among the first Neolithic groups of inner Iberia than in subsequent periods as well as differential access to some resources (possibly suckling herbivores) in the diet, which may point to the existence of early social or economic inequalities that do not seem to be linked to age and sex parameters or to mortuary treatment.

Genetic diversity and population structure of trifoliate yam (Dioscorea dumetorum Kunth) in Cameroon revealed by genotyping-by-sequencing (GBS)

BackgroundYams (Dioscorea spp.) are economically important food for millions of people in the humid and sub-humid tropics. Dioscorea dumetorum (Kunth) is the most nutritious among the eight-yam species, commonly grown and consumed in West and Central Africa. Despite these qualities, the storage ability of D. dumetorum is restricted by severe postharvest hardening of the tubers that can be addressed through concerted breeding efforts. The first step of any breeding program is bound to the study of genetic diversity. In this study, we used the Genotyping-By-Sequencing of Single Nucleotide Polymorphism (GBS-SNP) to investigate the genetic diversity and population structure of 44 accessions of D. dumetorum in Cameroon. Ploidy was inferred using flow cytometry and gbs2ploidy.ResultsWe obtained on average 6371 loci having at least information for 75% accessions. Based on 6457 unlinked SNPs, our results demonstrate that D. dumetorum is structured into four populations. We clearly identified, a western/north-western, a western, and south-western populations, suggesting that altitude and farmers-consumers preference are the decisive factors for differential adaptation and separation of these populations. Bayesian and neighbor-joining clustering detected the highest genetic variability in D. dumetorum accessions from the south-western region. This variation is likely due to larger breeding efforts in the region as shown by gene flow between D. dumetorum accessions from the south-western region inferred by maximum likelihood. Ploidy analysis revealed diploid and triploid levels in D. dumetorum accessions with mostly diploid accessions (77%). Male and female accessions were mostly triploid (75%) and diploid (69%), respectively. The 1C genome size values of D. dumetorum accessions were on average 0.333 ± 0.009 pg and 0.519 ± 0.004 pg for diploids and triploids, respectively.ConclusionsGermplasm characterization, population structure and ploidy are an essential basic information in a breeding program as well as for conservation of intraspecific diversity. Thus, results obtained in this study provide valuable information for the improvement and conservation of D. dumetorum. Moreover, GBS appears as an efficient powerful tool to detect intraspecific variation.

Resolving species delimitation within the genus Bunopus Blanford, 1874 (Squamata: Gekkonidae) in Iran using DNA barcoding approach.

Mitochondrial COI sequences were used to investigate species delimitation within the genus Bunopus in Iran. A dataset with a final sequence length of 633 nucleotides including 100 specimens from 31 geographically distant localities across Iran were generated. The result demonstrated that two major clades with strong support can be identified within the genus Bunopus in Iran. Clade A includes Bunopus crassicaudus and two new entities, eastern populations (subclade A2,1) and Shahdad populations (subclade A2,2). The second clade comprises western and southwestern populations (subclade B1,1), Arabian populations (subclade B1,2) and south and southeast populations in Iran, to which Bunopus tuberculatus (subclade B2) is assigned. In addition to Bunopus crassicaudus and B. tuberculatus, three new candidate species in Iran can easily be identified based on the DNA barcoding approach.

Local selection in the presence of high levels of gene flow: Evidence of heterogeneous insecticide selection pressure across Ugandan Culex quinquefasciatus populations.

BackgroundCulex quinquefasciatus collected in Uganda, where no vector control interventions directly targeting this species have been conducted, was used as a model to determine if it is possible to detect heterogeneities in selection pressure driven by insecticide application targeting other insect species.Methodology/Principal findingsPopulation genetic structure was assessed through microsatellite analysis, and the impact of insecticide pressure by genotyping two target-site mutations, Vgsc-1014F of the voltage-gated sodium channel target of pyrethroid and DDT insecticides, and Ace1-119S of the acetylcholinesterase gene, target of carbamate and organophosphate insecticides. No significant differences in genetic diversity were observed among populations by microsatellite markers with HE ranging from 0.597 to 0.612 and low, but significant, genetic differentiation among populations (FST = 0.019, P = 0.001). By contrast, the insecticide-resistance markers display heterogeneous allelic distributions with significant differences detected between Central Ugandan (urban) populations relative to Eastern and Southwestern (rural) populations. In the central region, a frequency of 62% for Vgsc-1014F, and 32% for the Ace1-119S resistant allele were observed. Conversely, in both Eastern and Southwestern regions the Vgsc-1014F alleles were close to fixation, whilst Ace1-119S allele frequency was 12% (although frequencies may be underestimated due to copy number variation at both loci).Conclusions/SignificanceTaken together, the microsatellite and both insecticide resistance target-site markers provide evidence that in the face of intense gene flow among populations, disjunction in resistance frequencies arise due to intense local selection pressures despite an absence of insecticidal control interventions targeting Culex.

Landscape genetics in the subterranean rodent Ctenomys “chasiquensis” associated with highly disturbed habitats from the southeastern Pampas region, Argentina

Studies of genetic differentiation in fragmented environments help us to identify those landscape features that most affect gene flow and dispersal patterns. Particularly, the assessment of the relative significance of intrinsic biological and environmental factors affecting the genetic structure of populations becomes crucial. In this work, we assess the current dispersal patterns and population structure of Ctenomys "chasiquensis", a vulnerable and endemic subterranean rodent distributed on a small area in Central Argentina, using 9 polymorphic microsatellite loci. We use landscape genetics approaches to assess the relationship between genetic connectivity among populations and environmental attributes. Our analyses show that populations of C. "chasiquensis" are moderately to highly structured at a regional level. This pattern is most likely the outcome of substantial gene flow on the more homogeneous sand dune habitat of the Northwest of its distributional range, in conjunction with an important degree of isolation of eastern and southwestern populations, where the optimal habitat is surrounded by a highly fragmented landscape. Landscape genetics analysis suggests that habitat quality and longitude were the environmental factors most strongly associated with genetic differentiation/uniqueness of populations. In conclusion, our results indicate an important genetic structure in this species, even at a small spatial scale, suggesting that contemporary habitat fragmentation increases population differentiation.

A range of possibilities: Assessing geographic variation in climate sensitivity of ponderosa pine using tree rings

Climate change has emerged as a major threat to biodiversity, but there have been few assessments of variability in sensitivity to climate across species’ ranges. We examined the sensitivity of ponderosa pine (Pinus ponderosa), a common, commercially important tree species found throughout the western United States, to climate variability using tree rings mostly from the International Tree Ring Data Bank (161 chronologies) and multiple climate products for 1930–1979. We compared interannual tree ring variation to water-year precipitation (wy ppt), May-July vapor pressure deficit (MJJ VPD), April-July precipitation (AMJJ ppt) and average winter (Dec-Feb) minimum temperature (winter tmn). We mapped growth-climate correlations for each variable and applied hierarchical agglomerative cluster analysis to identify five main groups of climate sensitivities, one of which contained mostly climate-insensitive populations. Ponderosa pine growth was positively correlated with wy ppt and AMJJ ppt and negatively correlated with MJJ VPD throughout its range, but sensitivities were greater in the Southern Rockies and Northern Great Plains than in Pacific and Southwestern populations. Sensitivities to winter tmn were weaker than responses to other climate variables, but Pacific and interior populations mostly demonstrated positive and negative responses, respectively. In summary, our analyses revealed the general pattern that interior populations of ponderosa pine were more climate-sensitive than Pacific populations, but there is considerable heterogeneity in sensitivity throughout the species’ range. Sensitivities to climate variables were related to position along corresponding climate gradients for all variables except MJJ VPD, but the unevenness of these relationships suggests a role of local adaptation. Our results illustrate the variation in climate response for a uniquely wide-ranging species and suggest that “leading” and “trailing” range edge concepts can oversimplify patterns of species exposure to climate change.

Genetic variation patterns of American chestnut populations at EST-SSRs

The objective of this study is to analyze patterns of genetic variation at genic expressed sequence tag – simple sequence repeats (EST-SSRs) and at chloroplast DNA markers in populations of American chestnut (Castanea dentata Borkh.) to assist in conservation and breeding efforts. Allelic diversity at EST-SSRs decreased significantly from southwest to northeast along the Appalachian range, suggesting repeated founder events during postglacial migration. Comparatively high allelic diversity in Ontario, northwest of the Appalachian range, suggested more recent long-distance dispersal. Clinal variation of allele frequencies along the Appalachian axis was also in accordance with postglacial colonization from one refugium southwest of the Appalachian range. We observed clustering of the northwestern population from Ontario with southwestern populations and sharing of a rare chloroplast haplotype among western populations across the whole latitudinal range. This pattern is consistent with a divergence of postglacial migration routes and higher levels of more recent potentially human-mediated gene exchange between populations west of the Appalachian range. Population pairs east and west of the Appalachian axis showed pronounced allele frequency differences over a small geographic range. These patterns of genetic variation should be considered when sampling reproductive material for conservation and breeding.

Age-related fecundity of free-ranging mule deer Odocoileus hemionus Cervidae in south-central, New Mexico, USA

Abstract Mule deer Odocoileus hemionus (Rafinesque 1817) populations in southern New Mexico have declined significantly since the 1980s, similar to trends across the arid Southwestern USA. Because production of fawns is critical to population growth, we evaluated factors influencing fecundity in two mule deer populations in southern New Mexico. Conception, litter size and survival of ≥1 fawn to weaning were all affected by maternal age, with older (age 8 and older) females exhibiting reproductive senescence as compared to prime-aged (age 2–7) females for the latter two traits despite achieving similar condition as did prime-aged females. Litter size and survival to weaning were also positively affected by increasing spring precipitation and survival to weaning was also positively affected by increased maternal condition, and size during late gestation. Unlike most previous work, reproductive senescence was evident in mule deer in our study populations, possibly because deer in both populations were in poor condition and older females produced on average 0.24 (95% CI=0.10–0.42) fawns through weaning compared to 0.76 (95% CI=0.60–0.94) for prime-aged females. The positive effect of precipitation during gestation on litter size and fawn survival also indicated that both income (i.e. nutritional intake) and capital (i.e. body reserves) were important determinants of fecundity in our arid Southwestern populations. The relatively early onset of senescence compared to the lifespan of female mule deer indicates that more intensive management of female age structure may be necessary to enhance population-level productivity.

Populations of aspen (Populus tremuloides Michx.) with different evolutionary histories differ in their climate occupancy.

Quaking aspens (Populus tremuloides Michx.) are found in diverse habitats throughout North America. While the biogeography of aspens' distribution has been documented, the drivers of the phenotypic diversity of aspen are still being explored. In our study, we examined differences in climate between northern and southwestern populations of aspen, finding large‐scale differences between the populations. Our results suggest that northern and southwestern populations live in distinct climates and support the inclusion of genetic and phenotypic data with species distribution modeling for predicting aspens' distribution.

Speciation in Cloudless Sulphurs Gleaned from Complete Genomes.

For 200 years, zoologists have relied on phenotypes to learn about the evolution of animals. A glance at the genotype, even through several gene markers, revolutionized our understanding of animal phylogeny. Recent advances in sequencing techniques allow researchers to study speciation mechanisms and the link between genotype and phenotype using complete genomes. We sequenced and assembled a complete genome of the Cloudless Sulphur (Phoebis sennae) from a single wild-caught specimen. This genome was used as reference to compare genomes of six specimens, three from the eastern populations (Oklahoma and north Texas), referred to as a subspecies Phoebis sennae eubule, and three from the southwestern populations (south Texas) known as a subspecies Phoebis sennae marcellina. While the two subspecies differ only subtly in phenotype and mitochondrial DNA, comparison of their complete genomes revealed consistent and significant differences, which are more prominent than those between tiger swallowtails Pterourus canadensis and Pterourus glaucus. The two sulphur taxa differed in histone methylation regulators, chromatin-associated proteins, circadian clock, and early development proteins. Despite being well separated on the whole-genome level, the two taxa show introgression, with gene flow mainly from P. s. marcellina to P. s. eubule. Functional analysis of introgressed genes reveals enrichment in transmembrane transporters. Many transporters are responsible for nutrient uptake, and their introgression may be of selective advantage for caterpillars to feed on more diverse food resources. Phylogenetically, complete genomes place family Pieridae away from Papilionidae, which is consistent with previous analyses based on several gene markers.

The first taxonomic revaluation of the Iranian water frogs of the genus Pelophylax (Anura: Ranidae) using sequences of the mitochondrial genome

The Eurasian water frog species and their geographic ranges have undergone considerable changes in the last four decades, but the Iranian populations have largely remained unknown. All the Iranian populations of water frogs, despite their vast distribution range have attributed to a single species: Rana ridibunda. In order to understand the phylogenetic relationships and taxonomic status of water frogs of Iran, we collected samples from many populations across the country and used the mitochondrial DNA sequence variation. A data set with a final sequence length of 616 nucleotides was generated for Cyt b from 70 individuals of Pelophylax in which there are 422 invariable sites, 174 variable sites of which 123 were parsimony informative. In total, 43 haplotypes were found (Hd: 0.9752). The result demonstrated that, two major clades with strong support can be identified within the Iranian water frogs. One of these clades that include north western and southwestern populations forms a monophyletic group along with P. bedriagae samples from Turkey. The second clade consists of water frog populations of north and northeastern parts of Iran which in turn is subdivided into two subclades. Inclusion of water frog samples from adjacent areas showed that the second clade of our study is, most likely, a distinct taxonomic entity at species rank with its two subclades indicating two diagnosable subspecies for the clade. In conclusion, we suggest that two distinct species, P. bedriagae and Pelophylax sp., with its two subspecies, should be identified as water frogs of Iran. In Addition, another traditionally reported water frog of Iran, P.ridibundus, most likely should be excluded from the Iranian water frog’s checklist.

Geographic distribution of Astroides calycularis (Scleractinia: Dendrophylliidae) as a baseline to assess future human impacts on the Southern Iberian Peninsula

Human activities have increasingly affected biodiversity in the Mediterranean Sea. Data on the distribution and abundance of species allows researchers to assess the possible degradation of wild populations. These data could act as a baseline to assess the magnitude of the effects of human activities on a bioindicator species. The distribution and relative abundance of the south-western populations of the endemic Astroides calycularis in the South Iberian Peninsula were studied to establish a baseline for future studies. The rocky shoreline was studied at a depth range of 0–12 m, including more than 650 km of Spain's Andalusian coastline. The species was present in 135 of the 585 dive points sampled. ANOVA analysis showed differences in depth in the four provinces studied, and there was no interaction between the two factors. As human activities on the Mediterranean coast are reducing the A. calycularis populations, a baseline on marine populations is greatly recommended for monitoring, assessment, and management studies, especially for endangered or bioindicator species. This baseline could be useful as a reference tool to assess the effects of human activities on marine biodiversity, including global change.

Geographic Variation of Morphological Characters in Bunopus tuberculatus Blanford 1874 (Squamata: Gekkonidae) in Iran

A species with populations that occur in variety of habitats might be expected to exhibit variation in morphological traits across its distribution. Based on our examination of 95 specimens of Bunopus tuberculatus (Squamata: Gekkonidae) collected across its range within Iran, we describe the geographic variation in meristic traits and morphology. Values for traits were compared using principal component analysis and discriminant function analysis. None of the populations showed sexual dimorphism in morphometric or meristic characters. One-way analyses of variance revealed that differences exist in most characters among populations. Additionally, canonical variable analyses indicate that the sampled populations of B. tuberculatus are divided into four distinct forms: (1) western and southwestern populations, (2) southern populations, (3) the population of Bampoor, and (4) eastern populations. Our results provide a comparative baseline that will help to resolve the complicated taxonomy of B. tuberculatus across the rest of its distribution in southwest Asia and the Arabian Peninsula.