Aphanius arakensis, a new species of tooth-carp (Actinopterygii, Cyprinodontidae) from the endorheic Namak Lake basin in Iran

Abstract A new species of tooth-carp, Aphanius arakensis sp. n., is described from the Namak Lake basin in Iran. The new species is distinguished by the congeners distributed in Iran by the following combination of characters: 10–12 anal fin rays, 28–32 lateral line scales, 10–13 caudal peduncle scales, 8–10 gill rakers, 12–19, commonly 15–16, clearly defined flank bars in males, a more prominent pigmentation along the flank added by relatively big blotches in the middle and posterior flank segments in females, a short but high antirostrum of the otolith that has a wide excisura, and a ventral rim with some small, drop-like processes, and 19 molecular apomorphies (17 transitions, two transversions) in the cytochrome b gene. It was suggested based on the phylogenetic analysis that the new species is sister to Aphanius sophiae from the Kor River and that Aphanius farsicus from the Maharlu Lake basin is sister to Aphanius arakensis plus Aphanius sophiae. A noticeable feature of the Aphanius diversity in Iran is the conservatism of the external morphology as well as morphometric and meristic characters, while distinctive differences are present in genetic characters, otolith morphology, and male color pattern. Transformation of the latter was probably driven by sexual selection.


Introduction
Aphanius is the only representative of the Cyprinodontidae (Teleostei, Cyprinodontiformes) in Eurasia. The genus occurs in coastal (brackish) and landlocked (freshwater to saline) water bodies in the Mediterranean and Persian Gulf basins from Iberian Peninsula as far eastwards as Iran and Pakistan (Wildekamp 1993). Aphanius species diversity is highest in the endorheic basins of the mountainous regions of central Anatolia and the Iranian plateau (Coad 2000;Hrbek and Meyer 2003, Hrbek et al. 2006, Esmaeili et al. 2012. Though central Anatolia is believed to represent the center of Aphanius speciation (Wildekamp et al. 1999), a high number of Aphanius species also occurs in Iran. Apart from the widely distributed A. dispar (Rüppell, 1829), seven endemic Aphanius species have been described from Iran to date, namely A. ginaonis (Holly, 1929) from the Genow hot spring near the Persian Gulf; A. isfahanensis Hrbek, Keivany & Coad, 2006 from the endorheic Esfahan basin; A. farsicus Teimori, Esmaeili and Reichenbacher, 2011 from the endorheic Maharlu Lake basin [A. farsicus is a replacement name for the previous A. persicus (Jenkins, 1910) because this name has been recognized as a homonym of the fossil A. persicus (Priem, 1908) (Gaudant 2011, Teimori et al. 2011; A. sophiae (Heckel, 1849) from the endorheic Kor River Basin; A. vladykovi Coad, 1988 from the upper reaches of the Karoun basin; A. mesopotamicus Coad, 2009 from the Tigris-Euphrates drainage; and the recently re-established A. pluristriatus (Jenkins, 1910) from the Mond River drainage. In addition to the species listed above, Lebias punctatus and Lebias crystallodon were originally described from the Nemek Deria near Shiraz by Heckel (1846Heckel ( -1849. Berg (1949) and Coad (1996) considered L. punctatus to be a synonym of A. sophiae but at that time most of now valid species distributed in Iran were thought to be synonyms of the widely distributed A. sophiae. Coad (1996) strongly suggested that the type locality of L. punctatus is not the Lake Maharlu but some other lake nearby as a name Nemek Deria is a very common name in Farsi for a salt lake. However, later, the Kotschy's itinerary in southern Iran in 1841 and 1842 was studied in detail based on botanical labels and it was clearly shown that collections by Kotschy studied by Heckel indeed came from a lake now called Maharlu (Edmondson and Lack 2006). This aspect is not in the focus of this very paper; we tentatively consider L. punctatus to be a synonym of A. sophiae until a proper examination of the extant syntypes of Lebias punctatus is done.
A number of isolated Aphanius populations that might deserve species status have been reported from endorheic drainages in Iran, but have not yet been investigated in detail (Coad and Abdoli 2000;Hrbek et al. 2006;Esmaeili et al. 2010).
They were commonly identified as A. sophiae (Heckel, 1849) (Coad and Abdoli 2000;Kamal et al. 2009); however, it was shown that the true A. sophiae is restricted to the endorheic Kor River basin near Shiraz (Fars Province) (Coad 2009;Esmaeili et al. 2012). This study describes a newly discovered Aphanius population from the Namak Lake basin in northern central Iran (Fig. 1). The specimens were collected in 2007 because they appeared to be different from other Iranian Aphanius species by a specific coloration. Here it is shown that the population from the Namak Lake basin in fact represents a new species, Aphanius arakensis. Our study is based on a total-evidence approach including morphometric and meristic characters, otolith morphology, and molecular data.

Material and methods
Institutional acronyms: ZM-CBSU, Zoological Museum of Shiraz University, Collection of Biology Department; ZSM, Zoological State Collection, Munich.

Morphological analysis
Based on the morphometric schemes introduced in Holcik et al. (1989) and Doadrio et al. (2002), 18 morphometric parameters were measured using a Vernier calliper and recorded to the nearest 0.5 mm. The standard length was measured from the most anterior part of the snout to the base of the caudal fin rays. In total, 21 relative variables were calculated from the measurements (Table 1). Scales removed from the left side of each fish, from the 3rd or 4th row below the dorsal fin, were mounted between microscope slides, and length and width of scales were measured to the nearest 0.1 mm by using a scale reader (Xerox 320). For each individual, scale length and scale width measurements were averaged to obtain a single length value and a single width value per individual and relative width and length of scales were calculated following Esmaeili (2001).
The meristic characters were counted under a stereomicroscope and consist of the numbers of (i) dorsal (ii) pectoral (iii) pelvic and (vi) anal fin rays, (v) lateral line series scales, (vi) caudal peduncle scales (the numbers of scales along the caudal peduncle, i.e. from the base of the last anal fin ray to the base of the caudal fin rays in a direct line), (vii) gill rakers and (viii) flank bars of males. Two posteriormost rays in dorsal and anal fins were calculated as one ray.
For examination of otolith morphology fish skulls were opened ventrally in order to remove the right and left otoliths. Otoliths were cleaned from tissue remains in 1% potassium hydroxide solution for 3-6 h, washed several times and finally rinsed in distilled water for 12 h. Otolith morphology was analyzed under a stereo microscope. In addition, five or six otoliths from each population were examined by a scanning electron microscope (SEM) with a LEO 1430 VP at ZSM.
Univariate analysis of variance (ANOVA, with Duncan's post hoc test, p < 0.05) was used to test the significance of phenotypic differences among species and also between sexes. Canonical discriminant analysis (CDA) was used for multivariate analyses in order to document the classification success of the groups. The statistical analyses were carried out using PASW 19.00 (SPSS Inc 2011) and PAST (Hammer et al. 2001: PAlaeontological STatistics, version 1.81).

Laboratory protocols and molecular analyses
Total genomic DNA was extracted according to phenol/chloroform procedures (Sambrook et al. 1989). A 900 base pairs (bp) fragment of the cytochrome b gene was successfully amplified via PCR using the primers (forward: Glu-F, 5' -AACCAC-  CGTTGTATTCAACTACAA-3'; reverse: ThrR, 5'-CCTCCGATCTTCGGATTA-CAAGACCG-3' (Machordom and Doadrio 2001). Amplification was performed in a thermal cycler programmed as follows: initial 94°C for 3 min, 35 cycles at 94°C for 50 s, 56°C for 45s, 72°C for 1 min, followed by a final extension at 72°C for 5 min. Sequencing was performed by Macrogen company, South Korea. Cytochrome b nucleotide sequences were edited with BioEdit and aligned through Geneious pro v5.4 (Drummond et al. 2011). Additional Aphanius sequences were obtained from the NCBI GenBank (http://www.ncbi.nlm.nih.gov) and included in the analyses (see above). The achieved cytb sequences for the here studied Aphanius populations were deposited in GenBank under numbers JX154880-JX154898. Maximum likelihood-based phylogenetic relationships were estimated by using the program SeaView version 4 (Gouy et al. 2010). The best-fit model of nucleotide substitution was obtained using the program JmodelTest 0.1.1 (Posada 2008). Accordingly, the GTR + I + G model (= General Time Reversible model + proportion of Invariable sites + Gamma-shaped distribution of rates across sites) was chosen.
Maximum parsimony based phylogenetic relationships were estimated using the program SeaView version 4 (Gouy et al. 2010) with 100 heuristic searches using random additions of sequences and implementing the Close-Neighbor-Interchange (CNI) on random tree algorithm. To test this phylogeny, bootstrap method using 2000 replication was used. To document the degree of homoplasy and degree to which potential synapomorphy is exhibited on the tree, the Consistency Index (CI) and the Retention Index (RI) were calculated by using the parsimony model within the Mesquite system for phylogenetic computing (Maddison and Maddison 2011).
The Neighbor Joining (NJ) distance-based phylogenetic relationships were estimated by using the computer program Geneious pro v5.4 (Drummond et al. 2011). The HKY85 model (Hasegawa et al. 1985) of molecular evolution was used with gamma distributed among site rate variation. There were a total of 771 positions in the final dataset.

Results
Aphanius arakensis sp. n. urn:lsid:zoobank.org:act:D9995F4C-AF0A-4791-9D80-D759EFEDA569 http://species-id.net/wiki/Aphanius_arakensis Figure 2A Diagnosis. The new species is distinguished by the congeners distributed in Iran by the following combination of characters: 10-12 anal fin rays, 28-32 lateral line scales, 10-13 caudal peduncle scales, 8-10 gill rakers, 12-19, commonly 11-13, clearly defined flank bars in males, a more prominent pigmentation along the flank added by relatively big blotches in the middle and posterior flank segments in females, a short but high antirostrum of the otolith that has a wide excisura, and a ventral rim with some small, drop-like processes and 19 molecular apomorphies (17 transitions, two transversions) in the cytochrome b gene.
Description of the holotype. The males of the new species reach approximately 32 mm SL and have 12-19 flank bars, the females are usually larger than the males and reach approximately 34 mm SL.
The morphometric characters are summarized in Table 1. Compared to the other examined Aphanius species, A. arakensis sp. n. shows higher mean values of the minimum body depth, width and length of scales, distances between the pectoral and pelvic fins and the interorbital distance, but significantly lower mean values for the eye diameter and the caudal peduncle length (differences are statistically significant, p < 0.05).
The meristic characters are summarized in Table 2. The dorsal fin is characterized by a somewhat curved superior border, and has 11-14 rays; the anal fin shows a round superior border and includes 10-12 rays; the pectoral fin is rounded and consists of 14-18 rays; the pelvic fin is relatively short, positioned just anteriorly to the anal fin and comprises 6-8 rays. The caudal fin is rounded; the caudal peduncle possesses 10-13 scales. The number of lateral line series scales is 27-32. However, the ANOVA analysis reveals that only the numbers of lateral line series scales and caudal  peduncle scales (in males and females), as well as the numbers of flank bars (in males), significantly differ from the values obtained for the other examined species. Moreover, there is a significant correlation between SL and numbers of flank bars (Pearson Correlation r = 0.455, p < 0.05*). The otolith is rounded-trapezoid and characterized by a very wide excisura, a medium-sized and pointed rostrum, and a quite short antirostrum. The ventral and dorsal rims are slightly curved; the ventral rim may bear small irregular processes; the dorsal rim may show a fine crenulation; the posterior rim is steep (Fig. 3W-Aa).
The flank bars in males (Fig. 2a) are narrow and the interspaces are broader than the bars. The first bar is located above the operculum, while the posteriormost bar is located at the base of the caudal fin; the interspaces are wider at the caudal peduncle than in the anterior body part. Dorsally, the head is gray and the body is dark due to a strong melanophore pigmentation. The ventral body portion does not usually show any dark pigmentation. The dorsal, anal and caudal fins have white margins; the first rays of the dorsal fin are dark. The pectoral fins are somewhat yellowish. The pelvic fin is yellowish. Most specimens are characterized by dark blotches at the base of the dorsal and anal fins. Females (Fig. 2b) are characterized by a grayish pigmentation of the back. The lateral flanks of the body are covered by dark pigmentations; series of blotches are present from the middle of the body to the caudal peduncle. The ventral part of the head and belly are light. The chin and sides of the head are speckled with melanophores. Below the eye there is a line of relatively dark melanophores. All fins are white.
Distribution and habitat. The species has been collected from a small natural shallow pond (Fig. 4) in the Namak Lake basin, 5 km south east of the city of Arak (Fig. 1). This pond, which is about 6 x 4 m in size, is fed by the drainage of a nearby natural spring. During sampling, the water body was almost stagnant and water temperature was 23°C. There was no vegetation in the pond, but the surrounding area was covered with Juncus sp. and Typha sp. The bottom of the pond was generally muddy with small gravels. The habitat was in a bad condition due to anthropogenic pollution. Around collection time, the new Aphanius species was the only fish observed living in the pond. In addition, the new species can be found in several springs located in close proximity to the type locality (Fig. 5).
Etymology. The species name refers to the city of Arak, which is located in close proximity to the type locality. Arak is the capital of the Markazi province in northcentral Iran. A proposed common name is Arak tooth-carp. Farsi name is Kapour-edandandar-e-Arak. . Natural shallow pond and type locality of Aphanius arakensis sp. n., in the Namak Lake Basin, 5 km SE of Arak city, Iran (see Fig. 1).

Phylogenetic relationships
The parameters for the maximum likelihood are ln(L) = -85.11.91237, gamma shape parameter of 1.000, proportion of invariant sites of 0.097 and parsimony = 1556. The maximum parsimony phylogeny has a CI of 0.462 and RI of 0.747. The initial tree for the maximum likelihood analysis was obtained by the BIONJ algorithm. The trees of the maximum likelihood and maximum parsimony phylogenies (Fig. 6) are not significantly different in topology (Templeton test, P > 0.05). They support the hypothesis that Aphanius arakensis diverged from the clade leading to the present-day A. sophiae and is sister to this species. Moreover, A. farsicus is sister to A. arakensis + A. sophiae; sister to these taxa is A. isfahanensis, and sister to all previously mentioned species is A. vladykovi. The same topology (Templeton test, P > 0.05) is observed for the tree of the Neighbor Joining (NJ) distance-based analysis. Table 4 shows the estimation of evolutionary divergence between the sequences of the new species and its relatives.

Probable reasons for morphological similarities between endemic Aphanius species
Several endemic Aphanius species are known that are soundly circumscribed by genetic differentiation and specific otolith morphology (see below), whereas they differ only weakly (or only in multivariate space) with regard to morphometry and meristics. Examples are A. isfahanensis from central Iran, A. sophiae and A. farsicus from southern Iran (Hrbek et al. 2006, this study); another example from the Mediterranean area is A. baeticus from Spain (Doadrio et al. 2002). Aphanius arakensis sp. n., from the Namak Lake basin represents yet another example for a species that is difficult to distinguish from its relatives based on external characters (with the exception of the features mentioned above). It is likely that the overall morphological similarity between these taxa are a result of the similar habitats, in which the various endemic Aphanius species are thriving. Thus, common environmental variables may have acted as a stabilizing selection on morphological characters (see also Hrbek et al. 2006). This offers an explanation as to why speciation events in Aphanius have affected genetic characters, rather than morphology, and why rapid genetic diversification can occur with little morphological change in this taxon (see also Adams et al. 2009).

Probable reasons for otolith differences between endemic Aphanius species
Otolith morphology is known to support the distinctive taxonomic state of several Aphanius species (Reichenbacher et al. 2007(Reichenbacher et al. , 2009a  and A. sophiae (Fig. 3) to show that these species are clearly different with regard to otolith morphology. Also A. arakensis shows clear divergence of its otolith morphology in comparison to the other inland Aphanius species, in particular with regard to the weakly pronounced antirostrum (Fig. 3W-Aa). Notably, the otoliths of A. vladykovi are most distinctive in comparison to those of the other studied species as they are characterized by a long ventral part, angular overall shape and long rostrum (Fig. 3R-V). This uniqueness of the A. vladykovi otoliths corresponds well to our and previous phylogenetic analyses, which have established A. vladykovi as being sister to all other Iranian inland species that diverged approximately 10 Ma ago (Hrbek et al. 2003). As a result, Aphanius likely has a higher rate of divergence in otolith morphology than in overall morphology. This difference in divergence rate may be related to the function of the otoliths as parts of the inner ear. In general, otoliths provide a mechanism for measuring motion and position of the head relative to gravity (Manley et al. 2004). However, it is quite important for a fish to know from where a sound is coming, so as to be able to distinguish between different sounds and pick out the biologically most relevant sounds (Popper et al. 2005). In addition, differences in otolith morphology are related to the balance and orientation of a fish (Popper et al. 2005). This means that differences in otolith morphology can reflect changes in intraspecific communication and behavior in fishes, that may have acted as evolutionary pressures.

Role of coloration pattern (flank bar numbers) in Aphanius diversification
Coloration and flank bar numbers are significant characters for the identification of Aphanius species, in particular for the identification of male individuals. Among the allopatric Iranian Aphanius species, males of A. arakensis have the largest number of flank bars, and flank bars are non-overlapping, whereas the number of flank bars is lowest in A. sophiae. Also the flank bars of the central Anatolian Aphanius species vary in thickness and number between species (Hrbek et al. 2002). However, the mechanisms underlying male flank bar variation have not been studied. We hypothesize that flank bar patterns play an important role in sexual selection, and thus represent important factors in the evolutionary history and speciation of Aphanius.
Sexual selection has long been believed to promote species divergence among groups of animals (see Kraaijeveld et al. 2010 for a review). Sexual selection may facilitate speciation because it can cause rapid evolutionary diversification of male mating signals and female preferences (Boughman 2001). Divergence in these traits may then contribute to reproductive isolation.
Several studies indicate that fishes can adapt to variation in underwater light environments by changing their colour, most likely as a result of a more effective intraspecific communication (Boughman 2001(Boughman , 2002Fuller 2002;Seehausen et al. 2008). Adding support to this interpretation is provided by studies on cichlids from the Victoria Lake (Seehausen et al. 2008) and African elephant fishes (Leal and Losos 2010). These studies indicate that variation in male nuptial coloration due to specific light conditions in different environments can result in ecological, phenotypic, genetic and behavioral differentiation. Additionally, color contrast with the visual background was found to be more important for effective intraspecific communication than color brightness (Fuller 2002). Thus, our conclusion is that the specific male flank bar patterns in different Aphanius species may have evolved as a response to different light regimes prevalent in respective habitats for increasing contrast and optimizing intraspecific communication. It can therefore be suggested that sensory-driven speciation might have played a prominent role in Aphanius speciation.

Conclusion
The noticeable features of the present-day diversity of the endemic Aphanius species in Iran include high genetic divergence and clear differences in otolith morphology, but only weak differences in general external morphology, morphometry and meristics. These patterns are probably caused by different rates of evolution in the mentioned characters that may be linked to the similarity of the individual environments, intraspecies communication, and vicariance events. It is likely that additional Aphanius species are present in remote areas of Iran, especially in the Zagros and Alburz Mountains.