Outline
Comptes Rendus

Taxonomy / Taxinomie
Surfing among species, populations and morphotypes: Inferring boundaries between two species of new world silversides (Atherinopsidae)
Comptes Rendus. Biologies, Volume 339 (2016) no. 1, pp. 10-23.

Abstract

Atherinopsidae are widespread freshwater and shallow marine fish with singular economic importance. Morphological, genetical and life cycles differences between marine and estuarine populations were already reported in this family, suggesting ongoing speciation. Also, coexistence and interbreeding between closely related species were documented. The aim of this study was to infer boundaries among: (A) Odontesthes bonariensis and Oargentinensis at species level, and intermediate morphs; (B) the population of Oargentinensis of Mar Chiquita Lagoon and its marine conspecifics. To achieve this, we integrated, meristic, Geometrics Morphometrics and DNA Barcode approaches. Four groups were discriminated and subsequently characterized according to their morphological traits, shape and meristic characters. No shared haplotypes between Obonariensis and Oargentinensis were found. Significative-meristic and body shape differences between the Mar Chiquita and marine individuals of Oargentinensis were found, suggesting they behave as well differentiated populations, or even incipient ecological species. The fact that the Odontesthes morphotypes shared haplotypes with both, Oargentinensis and Obonariensis, but also possess meristic and morphometric distinctive traits open new questions related to the origin of this morphogroup.

Metadata
Received:
Accepted:
Published online:
DOI: 10.1016/j.crvi.2015.11.004
Keywords: Atherinopsidae, Landmarks, Taxonomy, Meristic, Barcoding, Multivariate analysis

Mariano González-Castro 1, 2; Juan José Rosso 1, 2; Ezequiel Mabragaña 1, 2; Juan Martín Díaz de Astarloa 1, 2

1 Grupo de Biotaxonomía Morfológica y molecular de peces, IIMyC-CONICET, Universidad Nacional de Mar del Plata, Mar del Plata, Argentina
2 Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Buenos Aires, Argentina
@article{CRBIOL_2016__339_1_10_0,
     author = {Mariano Gonz\'alez-Castro and Juan Jos\'e Rosso and Ezequiel Mabraga\~na and Juan Mart{\'\i}n D{\'\i}az de Astarloa},
     title = {Surfing among species, populations and morphotypes: {Inferring} boundaries between two species of new world silversides {(Atherinopsidae)}},
     journal = {Comptes Rendus. Biologies},
     pages = {10--23},
     publisher = {Elsevier},
     volume = {339},
     number = {1},
     year = {2016},
     doi = {10.1016/j.crvi.2015.11.004},
     language = {en},
}
TY  - JOUR
AU  - Mariano González-Castro
AU  - Juan José Rosso
AU  - Ezequiel Mabragaña
AU  - Juan Martín Díaz de Astarloa
TI  - Surfing among species, populations and morphotypes: Inferring boundaries between two species of new world silversides (Atherinopsidae)
JO  - Comptes Rendus. Biologies
PY  - 2016
SP  - 10
EP  - 23
VL  - 339
IS  - 1
PB  - Elsevier
DO  - 10.1016/j.crvi.2015.11.004
LA  - en
ID  - CRBIOL_2016__339_1_10_0
ER  - 
%0 Journal Article
%A Mariano González-Castro
%A Juan José Rosso
%A Ezequiel Mabragaña
%A Juan Martín Díaz de Astarloa
%T Surfing among species, populations and morphotypes: Inferring boundaries between two species of new world silversides (Atherinopsidae)
%J Comptes Rendus. Biologies
%D 2016
%P 10-23
%V 339
%N 1
%I Elsevier
%R 10.1016/j.crvi.2015.11.004
%G en
%F CRBIOL_2016__339_1_10_0
Mariano González-Castro; Juan José Rosso; Ezequiel Mabragaña; Juan Martín Díaz de Astarloa. Surfing among species, populations and morphotypes: Inferring boundaries between two species of new world silversides (Atherinopsidae). Comptes Rendus. Biologies, Volume 339 (2016) no. 1, pp. 10-23. doi : 10.1016/j.crvi.2015.11.004. https://comptes-rendus.academie-sciences.fr/biologies/articles/10.1016/j.crvi.2015.11.004/

Version originale du texte intégral

1 Introduction

The order Atheriniformes is a monophyletic group diagnosed by ten characters, and composed by six families and 49 genera of generally small, silvery fish, which belong to the series Atherinomorpha, the most successful fish at the surface layer of the ocean and of many freshwaters habitats [1,2]. The members of this Order commonly share (among others) the following characters: usually two separate dorsal fins, the first, if present, with flexible spines and the second preceded by a supplementary simple spine; anal fin usually preceded by a thorn; lateral line very weak or absent; pectoral fin usually located high in the flanks, abdominal pelvic fins in most species [3]. The New World silversides (Atherinopsidae) are widespread freshwater and marine fish commonly occurring in schools in shallow waters [2]. It is a family of singular economic importance, normally employed for farming and for game fishing, but also commercially exploited by the artisanal and commercial coastal fleets [4]. In Argentina, this family is represented by nine species commonly known as “pejerreyes”, including: Odontesthes argentinensis (Valenciennes, 1835); Obonariensis (Valenciennes 1835); Oincisa (Jenyns 1841); Osmitti (Lahille, 1929); Ohatcheri (Eigenmann 1909); Ohumensis de Buen 1953; Oplatensis (Berg 1895); Onigricans (Richardson 1848) and Operugiae Evermann & Kendall 1906 [5,6].

Although most of Odontesthes species occurring in Argentina inhabit exclusively fresh or marine waters, Oargentinensis and Obonariensis, are commonly found co-existing in brackish waters. These species are valuable resources for regional fisheries, being marketed fresh, but also very appreciated for game fishing in lakes, coastal lagoons and marine coastal areas [7]. Both have been reported in Mar Chiquita coastal lagoon (Buenos Aires, Argentina) and its freshwater tributaries, being in fact Oargentinensis much more abundant than Obonariensis [8–10], which is a conspicuous silverside of the shallows lakes, small rivers and channels of the Pampa plain of Argentina [11]. However, distinguishing both species is not always straightforward for people with not expertise in fish identification.

Odontesthes argentinensis and Obonariensis are genetically [12] and morphologically [13] closely related. Recent molecular evidence suggested that both species could be distinguished using several genetic markers [12,14]. However, interbreeding among both species was also reported [15]. Thus, it is likely to find intermediate forms in environments where both species coexist (i.e. coastal lagoons). Moreover, current identification keys are ambiguous for several characters [5,16,17]. For instance, Dyer [5] indicated that Oargentinensis has 26 to 28 gill rakers (GR) in the lower branch and the first dorsal fin is situated over or posterior to the anus, while Obonariensis has 30–40 GR and first dorsal fin anterior to anus. However, Bemvenuti [16] and Cousseau et al. [17] stated that Oargentinensis has 19–25 GR in the lower branch. Unexpectedly, González-Castro et al. [9] found that most specimens of Oargentinensis collected in Mar Chiquita coastal lagoon shared characteristics of both species (20 to 25 GR and first dorsal fin anterior to anus).

Beheregaray and Sunnucks [18], based on molecular data, suggested that Oargentinensis is ongoing a recent speciation in brackish waters of South America. Moresco and Bemvenuti [19] stated that, in Río Grande do Sul (Brazil), Oargentinensis is represented by two populations: a resident population in the Patos Lagoon estuary and another one in the sea. Moreover, both populations showed evidence of spawning in its respective environment. Bemvenuti [20] found that the estuarine population of Oargentinensis has a spawning period during August–September. Moresco and Bemvenuti [19] characterized the reproductive biology of the marine population of this species and indicated that the spawning period extends from August to December. The ecological behaviour of Oargentinensis at higher latitudes than Brazil seems to be quite similar: there is a growing bulk of evidence about independent events of reproduction of marine and estuarine populations of this species. For example, the spawning of Oargentinensis in the Mar Chiquita coastal lagoon (Argentina) has been documented [9]. The authors found ripe and spent females in the inner zone of the lagoon (where waters are mixo-oligohaline) between June and November, suggesting the spawning of the Oargentinensis inside Mar Chiquita coastal lagoon. Cousseau [6] indicates that Oargentinenesis occur in coastal waters of the Province of Buenos Aires (Argentina) almost all over the year, but in late spring or early summer it “probably” migrates to estuarine waters with reproductive purposes. However, Llompart et al. [21] found that this species has a reproductive period in the coastal area of Bahía San Blas (Argentine, 40°S) between September and November, evidenced by an increase in the gonado-somatic index, with a peak in October corresponding to spawning. Therefore (as demonstrated for Brazil), the possibility of finding both, marine but also estuarine (or coastal lagoon) populations of Oargentinensis should not be ruled out.

All these previous results highlight a more than complex scenario where the accurate discrimination between Oargentinensis and Obonariensis (and its plausible morphs), but also between the co-existing populations of Oargentinensis are desirable. In this respect, we noted that there is a lack of meristic, morphometric and molecular comparative data between Oargentinensis and Obonariensis, but also between the presumptive coastal lagoon and marine populations of Oargentinensis. In this regard, a multidisciplinary approach is highly desirable for achieving valuable and robust results.

The aim of this study was to infer boundaries among: (A) Obonariensis and Oargentinensis at the species level, and the eventual intermediate morphs that could inhabit Mar Chiquita Coastal Lagoon; (B) the putative populations of Oargentinensis from Mar Chiquita Coastal Lagoon and marine environments. To achieve this, we integrated meristic, landmarks-based morphometrics and molecular (DNA Barcode sequences) approaches.

2 Materials and methods

2.1 Study area

Mar Chiquita coastal lagoon is a temperate shallow estuary, separated from the sea by a littoral line of dunes with an inlet joining it to the ocean. It covers approximately 60 km2, with a maximum length of 25 km parallel to the sea. This lagoon, considered a World Reserve of Biosphere by UNESCO [22], is located in the south-west Atlantic (37°32’–37°45'S–57°19’–57°26’W). Salinity fluctuates over a wide range between 0 and 36 PSU, and it is extremely variable and influenced by the freshwater volume present in the lagoon, the tide, and the wind direction/intensity [23]. The fish composition of Mar Chiquita coastal lagoon has been studied during the last fifteen years and several fish species have been reported to make extensive use of the lagoon, in a permanent, seasonal or occasional way [8–10,24–28]. Moreover, González-Castro et al. [9] not only analysed the spatial and temporal patterns in fish assemblage composition and relative abundance of fish species, but also evaluated the relative contribution of some environmental variables over these patterns.

2.2 Fish sampling

A total of 310 adult specimens belonging to genus Odontesthes were collected and analyzed, comprising three main continental, brackish and marine environments of the Province of Buenos Aires (Argentina): (a) fresh water environments, including Los Padres, La Brava, San Lorenzo, Chascomús, Pigué, Alsina and del Sol Lakes; (b) brackish water environments, i.e. Mar Chiquita coastal lagoon; (c) marine environments: Mar del Plata and Miramar coasts (Table 1).

Table 1

Groups, group-code, collection sites and sample size of the specimens used for this study. n, total number of specimens; nCOI, number of barcoded specimens; nM, number of specimens for meristic analyses; nIlD, number of specimens for morphometric analyses based on interlandmark distances; nGM: number of specimens for geometric morphometric analyses.

Groups Locality ST length range (mm) Group-code n n COI n M n IlD n GM
Oargentinensis Mar Chiquita coastal lagoon 133–322 Oarg_Mch 190 19 190 55 40
Oargentinensis Miramar coast; Mar del Plata coast 150–282 Oarg_marine 50 10 50 50 42
O. morphotypes Mar Chiquita coastal lagoon 232–326 Omorph 19 5 19 10 10
Obonariensis Mar Chiquita coastal lagoon. Lake La Brava; Lake Chascomús; Lake Pigué; Lake Alsina; Lake del Sol; Lake San Lorenzo 114–406 Obon 51 17 51 33 29
Total N 310 51 310 121 87

Fish were transported to the laboratory, where they were measured, weighted and sexed macroscopically. A small portion of tissue muscle was excised from representative specimens inhabiting the three sampled environments.

2.2.1 Taxonomic (morphological) identification

Each specimen was, at first, taxonomically identified by means of morphological identification keys of Bemvenutti [13], Cousseau et al. [17] and Dyer [5].

2.2.2 Meristic data analysis

Meristic characters were considered for each specimen as follows: series of lateral scales (Lss), gill rakers count in the lower branch of the first arch (Gr), number of branched and unbranched rays of first dorsal fin (D1), second dorsal fin rays (D2) and anal fin rays (A). Also, the relative position of the first dorsal fin (D1) in relation to the anus was recorded and coded as 1 (D1 anterior to anus), 2 (D1 over to anus), and 3 (D1 posterior to anus).

As the second objective of the present paper is to infer limits between the putative population of Oargentinensis of Mar Chiquita Coastal Lagoon and its marine conspecifics, we further refer to Oargentinensis from oligohaline environment (Mar Chiquita) as Oarg_Mch, and to those from marine localities (Mar del Plata, Miramar) as Oarg_marine.

The uni/multivariate meristic analyses were separately performed for each objective: (A) groups considered for analyses: Oargentinensis, Obonariensis and morphotypes and (B) groups considered for analyses: Oargentinensis from Mar Chiquita coastal Lagoon and marine environments (coast of Mar del Plata and Miramar).

A non-parametric Kruskal–Wallis one-way ANOVA on ranks, followed by the Kruskal–Wallis Multiple-Comparison Z-value test [29], were performed in order to evaluate significant differences for the variables studied among Odontesthes groups included in objective (A). For the case of the objective (B), the two groups were tested for statistical significance by means of the Mann–Whitney test.

On the other hand, the compositional similarity of the meristic variables between Odontesthes groups was determined by a non-parametric multivariate analysis (ANOSIM: analysis of similarities) test [30,31]. The ANOSIM test was used to search for differences in the meristic variables between “groups”. This permutation test analyzes differences among replicates within “groups”, contrasted with differences between “groups”, computing an R statistics under the null hypothesis of “no difference between groups”. The R statistics falls between –1 and 1, so that R is approximately 0 if the null hypothesis is true and R = 1 if all replicates within species are more similar to each other than any replicates from different species. A multidimensional scaling (MDS) analysis (group average sorting of the Bray–Curtis similarity measures based on non-transformed data) was also performed using the PRIMER software [31]. Similarity percentages (SIMPER) were used to identify which meristic variable characterized each “group” and which made the greatest contributions to any dissimilarity between “groups”.

Finally, the accuracy of classification of the specimens among the groups of Odontesthes, using meristic characters, was explored using a Discriminant Analysis (DA). The DA was followed by a Canonical Variate Analysis. This analysis is a special case of a Canonical Correlation Analysis for both, independent (meristic characters) and group (Odontesthes groups) variables that can be graphically displayed. All meristic variables were log transformed prior to this analysis.

2.2.3 Molecular genetics

Cytochrome oxidase subunit I (COI) was employed to test for species boundaries at the molecular level among groups tested in Objective A. Samples of white muscle tissue were excised from 51 individuals belonging to the Odontesthes groups (Table 1). In order to compare these groups with other Odontesthes species reported for the area, specimens of Oplatensis (n = 2), Osmitti (n = 3) and Oincisa (n = 1) were also included. Tissue muscles were preserved in 100% ethanol at –20 °C for genetic analysis. Specimens were labelled, photographed, formalin fixed (with further alcohol long-term preservation) and deposited as vouchers in the fish collection of the Universidad Nacional de Mar del Plata, Argentina.

DNA extraction, polymerase chain reaction (PCR), and sequencing of the COI gene were performed according to standard DNA barcoding protocols [32] and primer cocktails developed for fish [33,34]. Extraction and amplification were performed at the International Barcode of Life Argentinean Reference Barcode Laboratory of CONICET at the Museo Argentino de Ciencias Naturales, Buenos Aires, Argentina. Amplification of the 5′ region of COI, corresponding to base positions 6474 to 7126 of the Danio rerio mitochondrial genome [35], was first attempted using FF2d_t1/FR1d_t1 primer combination and C_FishF1t1/C_FishR1t1 primer cocktails [32]. The primer combinations C_FishF1t1 and C_FishR1t1 both contained two primers (FishF2_t1/VF2_t1 and FishR2_ t1/FR1d_t1, respectively). PCR reactions were performed in 96-well plates. The reaction master mix consisted of 825 μl water, 125 μl 106 buffer, 62.5 μl MgCl2 (25 mM), 6.25 μl dNTP (10 mM), 6.25 μl each primer (0.01 mM) and 6.25 μl Taq DNA polymerase (5 U/μl). This mixture was prepared for each plate, and each well contained 10.5 μl of solution and 2 μl of genomic DNA. The PCR reaction profile was comprised of an initial step of 2 min at 95 °C, and 35 cycles of 30 s at 94 °C, 40 s at 52 °C, and 1 min at 72 °C, with a final extension at 72 °C for 10 min. For specimens that failed to be amplified using the primer combinations above, the primer combinations C_VF1LFt1/C_VR1LRt1 (34) consisting of VF1_t1/VF1d_t1/LepF1_t1/VFli_t1 and VR1_t1/VR1d_t1/LepR1_t1/VRli_t1 primer sets respectively were tried. All primers were appended with M13 tails to facilitate sequencing.

Amplicons were visualized on a 2% agarose E-Gel H 96-well system (Invitrogen). Sequencing reactions applied M13 forward and reverse primers using the BigDyeH Terminator v.3.1 Cycle Sequencing Kit (Applied Biosystems Inc.), and the reaction profile was comprised of an initial step of 2 min at 96 °C and 35 cycles of 30 s at 96 °C, 15 s at 55 °C, and 4 min at 60 °C. Products were directly sequenced using an ABI 3730 capillary sequencer according to the manufacturer's instructions. Sequencing was performed at the Canadian Centre for DNA Barcoding (CCDB) in Ontario, Canada.

2.2.3.1 Molecular data analysis

DNA sequences were aligned using the Muscle [36] tool and further double-checked visually. Each distinct sequence was considered a different haplotype. The K2P + G model was chosen for comparison purposes, as it was determined as the best-fit model under Akaike information criterion for NJ, ML and MP analyses. A neighbour-joining (NJ) analysis was performed to provide a graphic representation of divergences between species [37].

Nevertheless, as distance-based models erase all character-based information [38], the best nucleotide substitution model was also employed to perform a maximum-likelihood (ML) and Maxima–Parsimonia analyses. All these analyses were performed using MEGA version 5.0 [37]. Robustness of trees was tested using bootstrap analysis [39] with 1000 replicates. The Barcode Index Number System (BINs) was used to reinforce species identification. BINs is “an online framework that clusters barcode sequences algorithmically, generating a web page for each cluster. Since clusters show high concordance with species, BINs can be used to verify species identifications as well as document diversity when taxonomic information is lacking.” [40].

Data analysis of COI sequences for the Odontesthes groups was complemented with the nucleotide diagnostic (ND) approach proposed by Wong et al. [41]. Nucleotide diagnostics, either simple (single nucleotide position) or compound (multiple single nucleotide positions; see Sarkar et al. [42]) are molecular characters that are unique for a particular species relative to the pool of species in which the ND was identified. Nucleotide diagnostics have already been shown to be useful in aid for species identification [43,44] and may further enhance DNA Barcode applications by overcoming the ambiguity inherent to the distance-based identification processes. Moreover, it has been recently demonstrated that using the ND approach as a complementary methodology of analysis was useful to reinforce the utility of the DNA barcoding technique to identify species for a large set of Neotropical fish species with low K2P divergence values (< 2%) [45].

All sequence assemblies, electropherogram (trace) files, primer sequences and specimen provenance data were deposited in the “Odontesthes of Argentina” Project (Project-code: OdArg) on BOLD (Barcode of life Data System). This includes digital images of the morphological voucher specimens, sex and ontogenetic stage (juvenile or adult), total and standard body length as well as GPS coordinates for all collection localities. The sequences were deposited in GenBank (accession numbers: pending).

2.2.4 Morphometry

Two types of variables were employed in order to achieve both objectives: (1) interlandmark distances (IlD); (2) coordinates data (landmarks). Accordingly, two different morphometric approaches were performed.

2.2.4.1 Morphometric analysis based on IlD

Twenty-five morphometric variables were taken as interlandmark distances over the left side of all specimens, using a digital calliper (0.05-mm precision). These variables were based on 12 landmarks obtained by truss network [46], defined on the basis of external anatomy and are homologous among the species (Fig. 1). Statistics and mathematics procedures for IlD analysis followed González-Castro et al. [47,48]. The morphometric characters were organized according to the groups defined for each objective: (A) groups considered for analyses, Oargentinensis, Obonariensis and morphotypes, and (B) groups considered for analyses, Oargentinensis from oligohaline environments and Oargentinensis from marine environments.

Fig. 1

Box truss showing the interlandmark distances collected in the silversides analysed, based on 12 anatomical landmarks. (Number of Box truss in Roman). Landmarks: 1, tip of premaxilla; 2, second nostril; 3, caudal extreme of the subopercular; 4, projection of the external border of the opercular towards the mid-dorsal line; 5, anterior insertion of the pelvic fin; 6, anterior insertion of the first dorsal fin; 7, anterior insertion of the anal fin; 8, anterior insertion of the second dorsal fin; 9, posterior insertion of the anal fin; 10, posterior insertion of the second dorsal fin; 11, ventral insertion of the caudal fin; 12, dorsal insertion of the caudal fin.

A normalization technique to scale the data that exhibit an allometric growth was used according to Lleonart et al. [49]. The standard length (SL) was used as the independent variable, whilst the remaining eight morphometric characters were considered as dependent ones. SL0 represents a reference value of size (250 mm in this paper) to which all individuals are reduced (or amplified) [50]. This transformation scales the data that exhibit allometric growth [49]. After transformation, a new matrix was constructed containing the corrected matrices for each species, and Principal Component Analysis (PCA) was performed using MULTIVARIADO® software [51]. Finally, principal components scores (PCs) were submitted to cross-validated discriminant analysis (DA) using SPSS® v.13.0, in order to build a predictive model of group membership based on observed characteristics of each case. This procedure generates a set of discriminant functions based on linear combinations of the predictor variables that provide the best discrimination between the groups. The functions are generated from a sample of cases for which group membership is known; the functions can then be applied to new cases with measurements for the predictor variables but unknown group membership.

2.2.4.2 Geometric morphometric analysis based on coordinate data (landmarks)

The analysis was performed on the Cartesian coordinates of the 12 anatomical landmarks, reconstructed from distance measurements among the landmarks (after its normalization procedure, as previously explained for IlD analysis), based on the proposed box truss scheme (Fig. 1) and using MORPPHEUS® Software [52]. Reconstruction of the form (shape plus size) from truss measures provides Cartesian coordinates for landmarks and allows estimation of, and compensation for, the measurement error [53].

The landmark coordinates for each specimen were scaled, translated and rotated using the generalized Procrustes superimposition (GLS, also called GPA). Scaling, translation and rotation was employed to minimize the Procrustes distance, the sum of squared distances between the corresponding landmarks. In summary, the objects were centred at the origin by subtracting the coordinates of its centroid from the coordinates of each landmark. Then each object was scaled (to unit centroid size) by dividing each coordinate of each landmark by its centroid size [54]. One landmark configuration was used as a reference and all others were rotated to minimize the partial Procrustes distance. The average shape was then calculated and became the new reference to which all objects were rotated again. This step was repeated until rotation ceased to occur.

The thin-plate spline (TPS) procedure was employed to compare shape differences among the groups, using both the uniform and non-uniform shape components and an upward/downward arching effect of the fish's body was observed. This effect was not related to biological factors (size or species) or to the preservation technique (freezing), but was rather due to slight posture differences between fish during interlandmarks distances capture. This distortion associated with the specimen's posture was already addressed in fish by Valentin et al. [55]. These authors proposed a method that effectively removes this artefact from the data, coupling a PCA-based model of the arching with Burnaby's orthogonal projection. This method also has the property of making the correction directly on landmark coordinates. Then, the methodology of Valentin et al. [55] was applied and the new unbiased coordinates were re-subjected to GLS and TPS (using TPSRELW Software (ver 1.46) [56]. The thin-plate spline is a method that projects data from shape space into a tangent space that is Euclidian and generates deformation grids, which depict shape changes over the entire object by interpolating between landmarks [54,57]. This method is not only an effective visualization tool, but its coefficients (partial warp scores) represent the non-uniform shape variation between specimens-consensus and can be used in descriptive and inferential statistical tests as well [58]. A principal components analysis (PCA) of the partial warps matrix was performed (usually named as relative warp analysis, RWA), in order to describe the major trends in shape variation. To examine the potential for differences in shape in classifying unknown specimens, the relative warps scores were submitted to discriminant analysis (SPSS ver. 13.0). This was carried out using cross-validation.

3 Results

3.1 Taxonomic (morphological) identification

The silversides analyzed were morphologically identified as:

  • Odontesthes bonariensis (Obon): all the specimens collected from Pampa Plain Lakes (n = 48) and a some (n = 3) from the inner (oligohaline) zone of Mar Chiquita Coastal Lagoon;
  • Odontesthes argentinensis (Oarg): all the specimens collected in the marine localities (n = 50) and specimens collected from the oligohaline zone of Mar Chiquita Coastal Lagoon (n = 190);
  • • morphotypes: in this morphological identification process, we detected several specimens (n = 19) that failed to be assigned to either Odontesthes bonariensis (Obon) or Odontesthes argentinensis (Oarg). They presented intermediate values for the gill rakers counts in the lower branch of the first gill arch. Moreover, these specimens did not match any diagnosis of the remaining species reported for the area. Consequently, this third group was regarded as morphotypes (Omorph) for further statistical analyses.

3.2 Meristic

The basic statistics on meristic data for the four groups identified in the present work are summarized in Table 2.

Table 2

Basic statistics on meristic data for the four groups of silversides analysed.

Groups GR Lss D1 D2 A D1-anus (%)
Mean (SD) Range Mean (SD) Range Mean (SD) Range Mean (SD) Range Mean (SD) Range Ant Over/post
Oarg_Mch 21.7 (1.1) 17–24 53.2 (2.0) 50–58 4.9 (0.6) 4–6 10.0 (0.6) 9–12 18.3 (0.8) 16–20 64.2 35.8
Oarg_marine 23.8 (1.0) 21–25 55.4 (1.7) 52–59 4.8 (0.6) 4–6 10.0 (0.7) 8–11 19.5 (0.9) 18–21 16.0 84.0
Omorph 27.3 (1.0) 25–29 58.6 (2.6) 53–64 4.8 (0.5) 4–6 11.1 (0.7) 10–12 19.2 (0.5) 19–20 61.1 38.9
Obon 32.9 (1.6) 30–36 58.0 (2.6) 52–66 5.0 (0.6) 4–6 10.6 (0.7) 9–12 17.7 (0.7) 15–19 68.6 31.4

3.2.1 Objective A

The Kruskal–Wallis one-way ANOVA on ranks showed that both species of Odontesthes and the morphotypes significantly differed in all but D1 meristic characters. The largest difference was found in the number of soft rays in the anal fin which were significantly (χ2: 54.58; P < 0.001) different in all pairwise comparisons. Odontesthes bonariensis and morphotypes also displayed a higher number of scales in the lateral line (χ2: 110.95; P < 0.001), a higher number of rays in D2 (χ2: 54.34; P < 0.001) and a higher gill raker count (χ2: 165.77; P < 0.001) than Oargentinensis. The analysis of contingency showed that the D1 fin is more (χ2: 10.77; P < 0.029) frequently inserted anterior to the anus in all three groups.

Analyses of similarities showed that there was significant differences in overall among groups in terms of meristic variables (Global test of the ANOSIM: R = 0.92; P = 0.1%). Pairwise comparisons within ANOSIM showed significant differences between all compared pair groups (Tables 3A and 3B). According to the ANOSIM test, the highest differences occurred between Oarg and Obon (Tables 3A and 3B). These results were also evident in the ordination analysis performed (MDS), in which specimens of the same group were clustered together as well as groups distinguished each other (Fig. 2A). SIMPER identified Gr and Lss as the meristic variables responsible for these differences (Tables 3A and 3B).

Table 3A

The contribution of meristic variables to observed differences among groups determined by SIMPER analyses.

A
ANOSIM SIMPER
Groups pairwise RS SL AA1 AA2 AD D/Sd C% Cu%
Oarg - Obon 0.984 0.1
 GR 22.18 32.94 4.55 5.11 60.09 60.09
 Lss 53.69 58.08 1.97 1.58 25.98 86.07
 A 18.61 17.74 0.46 1.10 6.02 92.08
Oarg - Omorph 0.811 0.1
 GR 22.18 27.58 2.31 2.71 41.36 41.36
 Lss 53.69 58.58 2.18 1.72 38.88 80.24
 D2 10.09 19.16 0.47 1.41 8.42 88.66
 A 18.61 11.11 0.43 1.22 7.63 96.30
Obon- Omorph 0.692 0.1
 GR 32.94 27.58 2.17 2.63 47.15 49.35
 Lss 58.08 58.58 1.11 1.23 25.20 74.55
 A 17.74 19.16 0.61 1.80 13.81 88.36
 D2 10.60 11.11 0.31 1.02 7.05 95.41
Table 3B

The contribution of meristic variables to observed differences among groups determined by SIMPER analyses.

B
SIMPER
Groups pairwise AA1 AA2 AD D/Sd C% Cu%
Lss 55.42 53.23 1.26 1.42 37.55 37.55
GR 23.80 21.75 0.96 1.62 28.81 66.35
A 19.52 18.37 0.58 1.31 17.42 83.77
D2 10.06 10.09 0.28 0.97 8.45 92.22
Fig. 2

a: multidimensional scaling (MDS) performed on meristic data in order to infer boundaries among Obonariensis (O_bon) (white triangles), Oargentinensis (Oarg) (black cruxes), and the intermediate morphs (O_morph) (black squares); b: multidimensional scaling (MDS) performed on meristic data in order to infer boundaries among the putative populations of Oargentinensis from Mar Chiquita coastal lagoon (Oarg_Mch) (white circles) and Oargentinensis from marine coastal areas (Oarg_marine) (black circles).

In the DA, the linear discrimination function that was used to elucidate the a priori classification significantly included all but D1 meristic characters. The DA further showed that individuals from the three groups were correctly classified in 99%. The Canonical Variate Analysis significantly extracted all the information (100%) in the first two variates. The spatial ordination of samples along the first and second variates is depicted in Fig. 3.

Fig. 3

Discriminant analysis performed on meristic data to infer boundaries among Obonariensis, Oargentinensis, and the intermediate morphs. Symbols: Odontesthes bonariensis (O_bon) (white triangles); Odontesthes morphotypes (O_morph) (black squares); Oargentinensis (black cruxes).

3.2.2 Objective B

The Mann–Whitney test showed that lateral line scales (Z-value: 6.25), number of soft rays in the anal fin (Z-value: 7.2) and gill raker counts (Z-value: 8.94) significantly differed between the two putative populations of silversides. Particularly, marine wanderers of Oargentinensis displayed a higher number of lateral line scales, soft anal fin rays and gill rakers than those of Oargentinensis from Mar Chiquita. In addition, marine specimens displayed a more (χ2: 116.99; P < 0.0001) posteriorly inserted first dorsal fin, and only 8 out of 50 specimens presented this fin inserted anterior to anus. Conversely, 122 of 190 specimens from Mar Chiquita displayed the insertion of first dorsal fin anterior to the anus.

The number of rays in both dorsal fins did not significantly contribute to the linear discrimination function of populations of Oargentinensis. DA further discriminated both populations by means of the remaining meristic characters. A percentage of 79.2% of correct classification for group membership was achieved. Overall, group classification reached 90% and 89.5% for marine and Mar Chiquita populations, respectively.

Analyses of similarities showed significant difference among both Oarg (Mar Chiquita vss marine) groups in terms of meristic variables (Global test of the ANOSIM: R = 0.373; P = 0.1%). These results were also observed in the performed ordination analysis (MDS), in which most specimens of the same group were clustered together as well as groups distinguished each other (Fig. 2B). SIMPER identified Gr and Lss as the meristic variables responsible for these differences (Tables 3A and 3B).

3.3 Molecular analysis

As COI is not intended to be used at population level, the results are mainly focused on objective A.

The NJ, ML and MP analyses (based on K2P + G model) generated trees with nearly identical topologies (data not shown). Odontesthes platensis, Osmitti and Oincisa clustered separately and were assigned to three different BINs. Conversely, Oargentinensis and Obonariensis clustered together and were assigned to the same BIN (data not shown). Interestingly, no shared haplotypes were found between Obonariensis and Oargentinensis specimens. Moreover, a compound nucleotide diagnostic allowed us to discriminate between Oargentinensis and Obonariensis by the exclusive occurrence of characters in two informative sites (124 T + 421 A).

Sequence analysis revealed 13 variable sites corresponding to 10 haplotypes, distributed within Obon, Oarg_marine, Oarg_Mch and Omorph groups (Table 4). Indeed, Obon displayed only two different haplotypes, being both shared with Omorph specimens. The remaining eight haplotypes corresponded to Oargentinensis: one of them was shared with a single Omorph specimen (Hp 3), two were unique haplotypes of Oarg_marine (Hp 6 and Hp 9), and three were unique haplotypes for Oarg_Mch (Hp's 7, 8 and 10) (Table 4).

Table 4

Distribution of the haplotypes (Hp) of the four groups analysed. Group-code as in Table 1.

Groups n Hp 1 Hp 2 Hp 3 Hp 4 Hp 5 Hp 6 Hp 7 Hp 8 Hp 9 Hp 10
Obon 17 13 4
Omorph 5 2 2 1
Oarg_marine 10 1 6 1 1 1
Oarg_Mch 19 12 2 2 1 1 1
Total 51 15 6 14 8 3 1 1 1 1 1

3.4 Morphometry

3.4.1 Objective A

3.4.1.1 Morphometric analysis based on IlD

The 25 normalized interlandmark distances, which were analyzed by PCA of the correlation matrix, produced eight eigenvalues greater than one (data not shown). The first four PCs explained more than 63% of the variance in the data. Only correlations (between variables and components) higher than 0.59 were taken as significant (data not shown). PCA based on IlD allowed graphic segregation of the three groups analyzed, with a slightly degree of overlap between them (data not shown). In this respect, Obon was located in the fourth quadrant (data not shown), with higher loadings for the 1–4, 1–3, 2–4 interlandmark distances (which represents the head length) but also 5–7 (the distance between the origins of the ventral and anal fins). Moreover, Obon specimens displayed higher values for the Ild that constitutes the third Box Truss (PC1 vs. PC3, data not shown) and lower values for the 10–11, 8–9, 9–10 variables (interlandmark distances of fourth and third Box- Truss). The specimens of Omorph were located basically in the first quadrant (data not shown) and characterized by higher loadings for the 3–4 (head width) and 6–7/8–9 (variables related to the origin of first and second dorsal fins, with respect to the Anal fin); lower values were obtained for the 4–6 variable, which represent the distance between the end of the head and the origin of the first dorsal fin, thus indicating that this fin originates very near to the head, when compared with the other Odontesthes groups analysed (data not shown). In this respect, Oarg was the most spread group and showed higher loadings for the 4–6 8–9, 9–10 and 10–11 variables and lower ones for 2–4, 3–4, 1–4 (head shape) and 6–7.

The data corresponding to the 25 PCs of the PCA where employed to perform the Discriminant Analysis. DA for the 148 individuals of Odontesthes produced two significant canonical discrimination functions. The first one explained 90.0% of the total variance in the data, (Wilks’lambda = 0.092, P < 0.000). The DA correctly classified 96.6% of the Odontesthes individuals according to the three groups defined a priori, whereas the cross-validated analysis correctly classified 90.2% of the fish according to their body shape (Table 5). Accordingly, group misclassifications were scarce, with a highest rate of 10.0% of Omorph misclassified as Oarg (Table 5). Three groups were defined, accordingly to those defined a priori, and their centroids and individuals were separated both on the first and second discriminant functions (Fig. 4). Remarkably, there was no overlap for the centroids/specimens of Obon, Omorph and Oarg (Fig. 4).

Table 5

Cross-validated Discriminant Analysis for the Objective A, based on the PĆ scores of: (1) interlandmark distances and (2) landmark coordinates. Group-code as in Table 1.

Predicted Group Membership (%)
Morphometric approach Objective Groups Obon Omorph Oarg Total
Percent
1 A Obon 93.94 6.06 0 100
Omorph 10 80 10 100
Oarg 2.86 6.67 90.48 100
Note: 90.5% of cross-validated grouped cases correctly classified
2 A Obon 82.8 13.8 3.4 100
Omorph 20.0 50.5 30.0 100
Oarg 2.4 9.8 87.8 100
Note: 83.5% of cross-validated grouped cases correctly classified
Fig. 4

Discriminant Analysis based on interlandmarks distances, in order to infer boundaries among Obonariensis, Oargentinensis, and the eventual intermediate morphs. Symbols: centroid of each species group (big grey diamonds); Odontesthes bonariensis (O_bon) (white triangles); Odontesthes morphotypes (O_morph) (black squares); Oargentinensis (black cruxes).

3.4.1.2 Geometric morphometric analysis based on coordinates (landmark data)

The first four RW explained 68.13% (31.76, 17.72, 12.30 and 6.36% respectively) of the total variance for the GLS/RWA analysis of the body shape of the three silversides groups analysed in this Objective A. The pattern of morphological variation described by the first two relative warps is shown in Fig. 5.

Fig. 5

Relative Warp analysis (RW1 vs. RW2) based on landmarks coordinates to infer boundaries among Obonariensis, Oargentinensis, and the intermediate morphs. Thin-plate spline transformation grids for the extreme points of RW1/RW2 are shown; they were superimposed on the shapes predicted when the average landmark configuration of all specimens was deformed into that of a hypothetical specimen positioned at the extreme of the RW of interest. Symbols: Odontesthes bonariensis (O_bon) (white triangles); Odontesthes morphotypes (O_morph) (black squares); Oargentinensis (black cruxes).

The shape variation along the first RW was basically expressed by the depression (negative RW1 scores) or expansion (positive RW1 scores) of the body and head along the dorso-ventral axis (i.e. the body height along by the second, third and fourth box trusses), but also the placement and length of the dorsal and anal fins (Fig. 5). Accordingly, the head shape vary from shorter and depressed (RW1−), which is typical of those Oargentinensis specimens located in the third quadrant, to deeper (in height) and longer (RW1+), as can be observed in O_morph and those samples of Oargentinensis located in the fourth quadrant (Fig. 5). The relative position of the dorsal fins (D1 and D2) was in accordance with the morphological traits observed in the previous analyses: the Oargentinensis specimens (those related to RW1) exhibited D1 posterior to anus and, a shorter D2; in the opposite, the individuals of O_morph and Oargentinensis related to (RW1+) showed the D1 located anterior to anus (Fig. 5). Moreover, the third and fourth box trusses of the O_morph (and Oargentinensis related to RW1+) individuals changed dramatically: D2 is forward displaced and the ventral and anal fins (particularly its origin) are downwardly displaced.

The shape of the RW2 basically allowed differentiating O_bon specimens (RW2+) from the remaining groups. The former was characterized by having the shortest caudal peduncle, the longest snout, and a longer head; a forward displacement of the origin of ventral fin, backward displacement of D1 fin and a smaller anal fin base (forward displacement of its insertion) (Fig. 5).

The data corresponding to the 20 RWs of the RWA were employed to perform the Discriminant Analysis. DA for the 121 original grouped cases of silversides classified by groups (Obon, Omorph and Oarg) produced two significant canonical discrimination functions, were the first one explained 87.3% of the total variance in the data (Wilks’ lambda = 0.164, P < 0.000). The DA correctly classified 90.1% of the original grouped cases, whereas the cross-validated analysis correctly classified 83.5% of the fish according to their body shape (Table 5). Three groups were defined, and its centroids and individuals were separated both on the first and second discriminant functions (data not shown). There was no overlap for the centroids of the three groups. Misclassifications (cross-validated analysis) ranged between 9.8–30.0% according the analysed groups (Table 5).

3.4.2 Objective B

3.4.2.1 Morphometric analysis based on IlD

The 25 normalized interlandmark distances, which were analyzed by PCA of the correlation matrix, produced eight eigenvalues greater than one (data not shown). The first four PCs explained more than 64% of the variance in the data. Only correlations (between variables and components) higher than 0.59 were taken as significant (data not shown). PCA based on IlD allowed almost a total graphic segregation of the two groups analyzed (data not shown). In this respect, Oarg_marine was located in the second and third quadrants (data not shown), with higher loadings for the 4–6, 3–6 interlandmark distances (which represents the distance between the posterior part of the head and the first dorsal fin), but also 7–9 (the anal fin base length), and lower values for the 3–4, 1–4, 2–4 and 2–3 variables (interlandmarks distances of first box truss, which represents the head) and 8–9, 8–10 variables (interlandmark distances of fourth box truss). The specimens of Oarg_Mch were located basically in the first and fourth quadrants (data not shown) and characterized by higher loadings of the variables related to the head shape (1–4, 2–3, 2–4, 3–4), and 7–8, 8–9, 8–10 and 9–10 (variables related to the fourth box truss, which represents the body at the level of the second dorsal and anal fins), and lower values for the 4–6 and 3–6 variables, thus indicating that the first dorsal fin is closer to the head, when compared with the marine Oargentinensis specimens.

The data corresponding to the 25 PCs of the PCA were employed to perform the Discriminant Analysis. DA for the 148 individuals of Odontesthes produced one significant canonical discrimination function, which explained 100% of the total variance in the data, (Wilks’ lambda = 0.059, P < 0.000). It is noticeable that both, the original and the cross-validated DA, correctly classified 100% of the Odontesthes individuals of the two groups included in this Objective B, according to their body shape.

3.4.2.2 Geometric morphometric analysis based on coordinates (landmarks data)

The first four RWs explained 75.38% (45.24, 11.65, 10.53 and 7.96% respectively) of the total variance for the GLS/RWA analysis of the body shape of the two Oargentinensis groups (Mar Chiquita and marine) analyzed in this Objective B. The pattern of morphological variation described by the first two relative warps is shown in Fig. 6. A complete segregation of the groups was observed, which explain the dispersion along the first RW axis observed in the RW analysis of Objective A (Fig. 5).

Fig. 6

Relative Warp analysis (RW1 vs. RW2) based on landmark coordinates performed on the populations of Oargentinensis from Mar Chiquita Coastal Lagoon and marine environment. Thin-plate spline transformation grids for the extreme points of RW1 are shown; they were superimposed on the shapes predicted when the average landmark configuration of all specimens was deformed into that of a hypothetical specimen positioned at the extreme of the RW of interest. Symbols: Oargentinensis from Mar Chiquita coastal lagoon (Oarg_Mch) (white circles); Oargentinensis from marine environments (Oarg_marine) (black circles).

The shape variation between the two groups analyzed in this objective B was exclusively related to the first RW. It was basically expressed by the depression (negative RW1 scores) or expansion (positive RW1 scores) of the head (as a consequence of the displacement of the landmarks 4 and 3) and the body (related to the displacement of the dorsal, ventral and anal fins) (Fig. 6). Accordingly, the head shape vary from shorter and depressed (RW1−) typical of Oarg_marine specimens, to deeper (in height) and longer (RW1+), as can be observed in those Oarg_Mch individuals (Fig. 6). The relative position of the dorsal fins (D1 and D2) was in accordance with the morphological traits observed in the previous analyses: the Oarg_marine specimens (RW1) exhibited a backward displacement of the D1 fin and a shorter D2; to the opposite, the Oarg_Mch individuals (RW1+) showed a forward displacement of the D1 fin (Fig. 6). Moreover, the third and fourth box trusses of the Oarg_Mch (RW1+) individuals have dramatically changed: D2 is forward displaced and the ventral and anal fins (its origins) are downwardly displaced; also, the insertion of the anal fin (landmark 9) is forwardly displaced, resulting in a narrower anal fin. At last, Oarg_marine individuals were also characterized by a wider anal and a smaller caudal peduncle (Fig. 6).

The data corresponding to the 20 RWs of the RWA were employed to perform the Discriminant Analysis. DA for the 83 original grouped cases of silversides classified by groups (Oarg_Mch and Oarg_marine) produced one significant canonical discrimination function (Wilks’ lambda = 0.087, P < 0.000). As in the morphometric analysis based on IlD, both the original and the cross-validated DA correctly classified 100% of the individuals of the two groups analyzed in this Objective, according to their body shape.

4 Discussion

The present study represents the first attempt to analyze together the meristic, genetic and landmarks-based morphometric characters of Oargentinensis and Obonariensis. This multidisciplinary approach allowed us to detect and characterize four groups by means of both their body shape and meristic features: Obonariensis, Oargentinensis from marine environments, Oargentinensis from oligohaline environments (Mar Chiquita Lagoon) and Odontesthes morphotypes. Odontesthes bonariensis was characterized by having 30–36 GR, 52–67 Lss, a long head, a deep body at the level of the second and third box trusses, the anal fin located posteriorly in the body and a narrow anal fin base; Oargentinensis from marine environments, have 21–25 Gr, 52–59 Lss, a small head (short and depressed) and dorsal fins (D1–D2) inserted posteriorly in the body; Oargentinensis from Mar Chiquita lagoon, showed 17–24 Gr, 48–58 Lss, the widest head, dorsal fins inserted anteriorly and a robust body at the third and fourth box trusses (Fig. 1). Finally, the group of Odontesthes morphotypes showed intermediate values of Gr/Lss between Obonariensis and Oargentinensis (25–29 and 53–64, respectively), a higher and rather long head too (higher and mid loadings of 3–4/1–4 variables, respectively) and a D1 inserted closer to the head (anterior to the anus), D2 forwardly displaced and anal fin downwardly displaced.

One of the morphological traits employed in the identification keys to distinguish Oargentinensis from Obonariensis is the position of first dorsal fin in relation to the anus [5,13,17]. These authors agreed that Oargentinensis (adult specimens) has the first dorsal fin (D1) situated over or posterior to the anus, while in Obonariensis the D1 is located anterior to the anus. Surprisingly, we found that 64% of O. arg_Mch (n = 190) presented the D1 anterior to the anus, while only 16% of the marine wanderers (O. arg_marine; n = 50) showed this state. Furthermore, 31.4% of the O. bon (n = 51) analyzed had the D1 over (or even posterior) to the anus. Finally, Odontesthes morphotypes (n = 19) presented the latter state in 39% of the specimens, while the remaining 61% displayed D1 anterior to the anus. These results evidenced that using the current taxonomic identification keys could lead to misidentifications. Therefore, we suggest that this character alone should not be included in the identification keys.

Sufficient evidence has been presented to accept the assumption that landmark-based morphometry can discriminate among species, populations and even morphotypes [48,59–62]. With respect to Atheriniformes, ÓReilly and Horn [63] and, more recently, Flucker et al. [64] and Crichigno et al. [65] presented landmarks-based studies, employing the body or head shape to elucidate the phenotypic variation, at the population or specific taxonomical level. In the present work, a multidisciplinary methodological approach was employed, which includes genetic, meristic and morphometric analyses. Moreover, the fact that two complementary morphometric approaches were used contributed to a better understanding of the taxonomic differences related to the body shape of these fish.

Depending on the kind of fish under study, achieving a standard or neutral posture for each individual is not straightforward. The fish body usually is not a rigid structure, and the specimen's shape could be influenced by its posture during landmark capture. This issue was recently stated by Valentin et al. [55], who detected an upward or downward arching effect in the morphometric data set of a multidisciplinary study on redfish (genus Sebastes) in the north-west Atlantic Ocean [66]. These authors proposed an approach, coupling a PCA-based model of the arching with Burnaby's orthogonal projection for removing the artefact from the data. This approach was successfully employed by González-Castro et al. [48] for a morphometric data set of seven species of Mugilidae (Actinopterygii). In the present work, the same kind of arching effect was encountered for the approaches based on landmark coordinate data. Valentine et al.’s [55] methodology was applied, and the arching effect was removed, yielding satisfactory results as evidenced by the RWA and the correct classification rates in the DA (Figs. 5 and 6).

COI sequences are probed to be highly effective to discriminate at the species level [33,67–69]. However, despite so far Oargentinensis and Obonariensis are regarded as two valid species, their COI sequences were assigned to the same BIN. A lack of discrimination among closely related fish species using COI was already reported [68,70,71]. Moreover, Heras and Roldan [12], using several mitochondrial DNA markers, discriminated between Obonariensis and Oargentinensis, but noted that “genetic distances between them for all molecular markers were within the range of Odontesthes intraspecific levels. Moreover, Obonariensis and Oargentinensis comprised a common lineage in all phylogenetical analyses consistent with their shared morphological characters.” In fact, the three specimens of Oargentinensis employed by Heras and Roldan [12] corresponded to the same sampling sites (Mar Chiquita lagoon and Mar del Plata coast) as those of the present work. Although no delimitation between species was obtained using DNA barcoding, no shared haplotypes between Obonariensis and Oargentinensis were found, and the compound nucleotide diagnostic character analysis allowed us to discriminate between Oargentinesis and Obonariensis. The fact that Omorph specimens shared haplotypes with both Oargentinensis and Obonariensis coupled with the meristic and morphometric diagnostic features of this species opens new questions about the origin of this group: could this morphogroup be the result of natural hybridization between Oargentinensis and Obonariensis (as was reported by Tejedor [15])? Should it just be considered a bunch of rare specimens resulted from the intraspecific morphological variability of Oargentinensis, which can be linked to their novel colonization to the estuarine habitat?

Odontesthes argentinensis is considered a widely distributed western Atlantic coastal species, occurring in marine and estuarine environments from the Sao Paulo State, in Brazil, southwards from the province of Chubut, Argentina [5]. Atherinids show a high degree of plasticity, which is inherent within their genetic make-up. This plasticity pre-adapts the species to radiate into habitats showing more stable environmental characteristics, representative of any part of their endemic range [72]. Moreover, the extent and pattern of divergence between estuarine and marine populations of Oargentinensis of southern Brazil indicated that speciation is occurring in the Patos Lagoon estuary [73]. This is an example of speciation associated with significant behavioural and ecological divergence. Indeed, while inshore sheltered waters and estuaries remain the principal habitat of the Atheriniformes, these fish show a striking ability to invade and speciate within vacant freshwater niches [72]. Bloom et al. [74] shown that transitions from marine to freshwater environments result in accelerated speciation and extinction rates, and that these rate differences may help explain the remarkable disparity in species richness between continents and oceans. García et al. [75] pointed out that promiscuous and recent contact between incipient species of silversides blurs species boundaries, yielding complicated taxonomy and species delimitation among silverside genus Odontesthes. Interestingly, in the present work, we found significative-meristic and body shape differences between the Mar Chiquita coastal lagoon and the marine populations of Oargentinensis, as it was already pointed out by Bemvenuti [76] for Oargentinensis of marine and estuarine areas of Rio Grande do Sul State (Brazil). Beheregaray and Sunnucks [18] suggested that ecological shifts due to colonization of estuarine habitats seem to have promoted rapid adaptive divergence and reproductive isolation in estuarine populations of Oargentinensis, which were considered as incipient ecological species. This conclusion is supported by the existence of a set of environmental factors required for successful reproduction of estuarine ecotypes. The estuarine (Oarg_Mch) and marine (Oarg_marine) populations of Oargentinensis studied here are meristically and morphometrically (this study), as well as genetically (Heras and Roldán [12] and this study), distinguishable and appear to behave as well differentiated populations, or even incipient ecological species according to Beheregaray and Sunnucks [18]. Even though COI haplotypes of estuarine and marine populations are shared (and thus no discrimination between both groups were obtained using NDC), a dominant haplotype in each group was observed. Indeed, Heras and Roldán [12], using several markers, found some minor genetic difference between both groups, although the sample used was relatively low (n = 3) to strongly support this. Moreover, González-Castro et al. [9] stated that there is a reproductive isolation of Oargentinensis in Mar Chiquita coastal lagoon. They found ripe and spent females in the inner zone of the lagoon (Zone III in González-Castro et al. [9]) between June and November, confirming reproductive events of Odontesthes argentinensis inside the coastal lagoon. All these results strongly suggest that each population should be treated separately with regard to management and conservation plans, specially taking into account that Mar Chiquita coastal lagoon is a World Reserve of Biosphere. It would be desirable to employ additional molecular markers (i.e., microsatellites or nuclear genes) in order to delimit and characterize this population/incipient species of Oargentinensis that inhabit Mar Chiquita coastal lagoon.

Disclosure of interest

The authors declare that they have no competing interest.

Acknowledgments

The authors would like to thank: Julio Mangiarotti (forest guard of Mar Chiquita Biosphere Reserve), Marcelo Pons (fishermen of Mar Chiquita), Mónica Iza, Florencia Celesia and Gladys Eiras (Mar Chiquita Town Hall) and Town Hall Mar Chiquita authorities (Luis Facca, Jorge Paredi); Matías Delpiani (BIMOPE Laboratory of UNMdP, Argentina) for collaboration in the field sampling; Darío Lijtmaer for assistance in barcode protocols. This work was funded by CIC (CIC 0823/05), UNMdP (15/E525, EXA 577/2) and PIP-CONICET No. 11220130100339CO grants.


References

[1] B.S. Dyer; B. Chernoff Phylogenetic relationships among atheriniform fishes (Teleostei: Atherinomorpha), Zool. J. Linn. Soc. Lond., Volume 117 (1996), pp. 1-69

[2] S. Helfman; B.B. Collette; D.E. Facey; B.W. Bowen The diversity of fishes, biology, evolution and ecology, Wiley-Blackwell, Oxford, 2009

[3] J.S. Nelson Fishes of the world, John Wiley & Sons, New York, 2006

[4] M.B. Cousseau; R.G. Perrotta Peces marinos de Argentina: biología, distribución, pesca, Instituto Nacional de Investigación y Desarrollo Pesquero INIDEP, Mar del Plata, 2013

[5] B.S. Dyer Systematic revision of the South American silversides (Teleostei, Atheriniformes), Biocell, Volume 30 (2006), pp. 69-88

[6] Ictiología. Aspectos Fundamentales. La vida de los peces sudamericanos (M.B. Cousseau, ed.), EUDEM, Mar del Plata, 2010

[7] R.P. Sánchez; G. Navarro; V. Rozycki Estadísticas de la pesca marina en la Argentina: evolución de los desembarques 1898–2010, Ministerio de Agricultura, Ganadería y Pesca, Buenos Aires, 2012

[8] M.B. Cousseau; J.M. Díaz de Astarloa; D.E. Figueroa La Ictiofauna de la Laguna Mar Chiquita (O. Iribarne, ed.), Reserva de Biósfera Mar Chiquita: Características físicas, biológicas y ecológicas, Editorial Martín, Mar del Plata, 2001

[9] M. González-Castro; J.M. Díaz de Astarloa; M.B. Cousseau; D.E. Figueroa; S.M. Delpiani; D.O. Bruno et al. Fish composition in a south-western Atlantic temperate coastal lagoon: spatial-temporal variation and relationships with environmental variables, J. Mar. Biol. Assoc. UK, Volume 89 (2009), pp. 593-604

[10] M.B. Cousseau; M.C. Marchesi; D.E. Figueroa; J.M. Díaz de Astarloa; M. González-Castro Relación íctica entre la laguna costera Mar Chiquita y el mar adyacente, Historia Natural, Volume 1 (2011), pp. 85-100

[11] J.J. Rosso Peces pampeanos: guía y ecología, Literature of Latin America, Buenos Aires, 2006

[12] S. Heras; M.I. Roldán Phylogenetic inference in Odontesthes and Atherina (Teleostei: Atheriniformes) with insights into ecological adaptation, C.R. Biol., Volume 334 (2011), pp. 273-281

[13] M.A. Bemvenuti Diferenciação morfológica das espécies de peixes-rei, Odontesthes Evermann & Kendall (Osteichthyes, Atherinopsidae) no extremo sul do Brasil: morfometria multivariada, Rev. Bras. Zool., Volume 19 (2002), pp. 251-287

[14] R.R. Betancur; R.E. Broughton; E.O. Wiley; K. Carpenter; J.A. López; C. Li et al. The tree of life and a new classification of bony fishes, PloS Curr., Volume 5 (2013) | DOI

[15] D. Tejedor El pejerrey como recurso genético (F. Grosman, ed.), Fundamentos biológicos, económicos y sociales para una correcta gestión del recurso pejerrey, Editorial Astyanax, Azul, 2001

[16] M.A. Bemvenuti Redescriçao do peixe-rei Odontesthes argentinensis (Valenciennes) Pisces: atherinidae, na costa do Rio Grande do Sul, Atlântica, Volume 15 (1993), pp. 17-35

[17] M.B. Cousseau; A.E. Gosztonyi; I. Elıas; M.E. Re Estado actual del conocimiento de los peces de la plataforma continental argentina y adyacencias (R.P. Sanchez; S.I. Bezzi, eds.), El mar argentino y sus recursos pesqueros Tomo 4, Los peces marinos de interés pesquero. Caracterización biológica y evaluación del estado de explotación, INIDEP, Mar del Plata, 2004, pp. 17-38

[18] L.B. Beheregaray; P. Sunnucks Fine-scale genetic structure, estuarine colonization and incipient speciation in the marine silverside fish Odontesthes argentinensis, Mol. Ecol., Volume 10 (2001), pp. 2849-2866

[19] A. Moresco; M.A. Bemvenuti Biologia reprodutiva do peixe-rei Odontesthes argentinensis (Valenciennes) (Atherinopsidae) da região marinha costeira do sul do Brasil, Revista Brasileira de Zoologia, Volume 23 (2006), pp. 1168-1174

[20] M.A. Bemvenuti Abundância, distribuição e reprodução de peixe-rei (Atherinidae) na zona estuarina da Lagoa dos Patos, RS, Brasil, Atlântica, Rio Grande, Volume 9 (1987), pp. 5-32

[21] F.M. Llompart; D.C. Colautti; T. Maiztegui; A.M. Cruz-Jimenez; C.R.M. Baigún Biological traits and growth patterns of pejerrey Odontesthes argentinensis, J. Fish Biol., Volume 82 (2013), pp. 458-474

[22] O. Iribarne Reserva de Biosfera Mar Chiquita. Características Físicas, Biológicas y Ecológicas, Editorial Martín, Mar del Plata, 2001

[23] R. Reta; P. Martos; G.M.E. Perillo; M.C. Piccolo; A. Ferrante Características hidrográficas del estuario de la Laguna Mar Chiquita (O. Iribarne, ed.), Reserva de Biósfera Mar Chiquita: características físicas, biológicas y ecológicas, Editorial Martín, Mar del Plata, 2001, pp. 31-52

[24] J.M. Díaz de Astarloa; D.E. Figueroa; M.B. Cousseau; M. Barragán Occurrence of Trachinotus carolinus (Carangidae) in Laguna Costera Mar Chiquita, with comments on other occasionally recorded fishes in Argentinean waters, B. Mar. Sci. (2000), pp. 399-404

[25] M. González-Castro; J.M. Díaz de Astarloa; M.B. Cousseau First record of a tropical affinity mullet, Mugil curema (Mugilidae), in a temperate southwestern Atlantic coastal lagoon, Cybium, Volume 30 (2006), pp. 90-91

[26] M. González-Castro; G.J. Macchi; M.B. Cousseau Studies on reproduction of the mullet Mugil platanus Günther, 1880 (Actinopterygii, Mugilidae) from the Mar Chiquita coastal lagoon, Argentina: similarities and differences with related species, Ital. J. Zool., Volume 78 (2011), pp. 343-353

[27] M. González-Castro; S.M. Delpiani; D.O. Bruno; J.M. Díaz de Astarloa First occurrence of the Patagonian blennie, Eleginops maclovinus (Cuvier, 1830) and the silverside Odontesthes smitti (Lahille, 1929), in a temperate south-western Atlantic coastal lagoon, J. Appl. Ichthyol., Volume 29 (2013), pp. 940-942

[28] M. González-Castro; J.J. Rosso; N.A. Lajud; D.L. Castellini; J.M. Díaz de Astarloa Occurrence of Cyprinus carpio (Linnaeus, 1758) in a World Biosphere Reserve, the Mar Chiquita Coastal Lagoon (Argentina), J. Appl. Ichthyol., Volume 31 (2015), pp. 749-751

[29] W.W. Daniel Applied Nonparametric Statistics, Houghton Mifflin Company, Boston, 1978

[30] K.R. Clarke; R.H. Green Statistical design and analysis for a “biological effects” study, Mar. Ecol. Prog. Ser., Volume 46 (1988), pp. 213-226

[31] K.R. Clarke; R.M. Warwick Changes in marine communities: an approach to statistical analysis and interpretation, PRIMER-E, Plymouth, 2001

[32] N.V. Ivanova; J.R. de Waard; P.D.N. Hebert An inexpensive, automation friendly protocol for recovering high-quality DNA, Mol. Ecol. Notes, Volume 6 (2006), pp. 998-1002

[33] R.D. Ward; T.S. Zemlak; B.H. Innes; P.R. Last; P.D.N. Hebert DNA barcoding Australia's fish species, Philos. T. Roy. Soc. B, Volume 360 (2005), pp. 1847-1857

[34] N.V. Ivanova; T.S. Zemlak; R.H. Hanner; P.D.N. Hebert Universal primer cocktails for fish DNA barcoding, Mol. Ecol. Notes, Volume 7 (2007), pp. 544-548

[35] R.E. Broughton; J.E. Milam; B.A. Roe The complete sequence of the zebrafish (Danio rerio) mitochondrial genome and evolutionary patterns in vertebrate mitochondrial DNA, Genome Res., Volume 11 (2001), pp. 1958-1967

[36] R.C. Edgar MUSCLE: multiple sequence alignment with high accuracy and high throughput, Nucleic Acids Res., Volume 32 (2004), pp. 1792-1797

[37] K. Tamura; D. Peterson; N. Peterson; G. Stecher; M. Nei; S. Kumar MEGA5: Molecular Evolutionary Genetics Analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods, Mol. Biol. Evol., Volume 28 (2011), pp. 2731-2739

[38] R. DeSalle Species discovery versus species identification in DNA barcoding efforts: response to Rubinoff, Conserv. Biol., Volume 20 (2006), pp. 1545-1547

[39] J. Felsenstein Confidence limits on phylogenies: an approach using the bootstrap, Evolution, Volume 39 (1985), pp. 783-791

[40] S. Ratnasingham; P.D.N. Hebert A DNA-based registry for all animal species: the Barcode Index Number (BIN) System, PLoS ONE, Volume 8 (2013), p. e66213

[41] E.H.K. Wong; M.S. Shivji; R.H. Hanner Identifying sharks with DNA barcodes: assessing the utility of a nucleotide diagnostic approach, Mol. Ecol. Resour., Volume 9 (2009), pp. 243-256

[42] I.N. Sarkar; P.J. Planet; T.E. Bael; S.E. Stanley; M. Siddall; R. DeSallee et al. Characteristic attributes in cancer microarrays, J. Biomed. Infor., Volume 35 (2002), pp. 111-122

[43] R.P. Kelly; I.N. Sarkar; D.J. Eernisse; R. DeSalle DNA barcoding using chitons (genus Mopalia), Mol. Ecol. Notes, Volume 7 (2007), pp. 177-183

[44] J. Rach; R. DeSalle; I.N. Sarkar; B. Schierwater; H. Hadrys Character-based DNA barcoding allows discrimination of genera, species and populations in Odonata, P. Roy. Soc. B Biol. Sci., Volume 275 (2008), pp. 237-247

[45] L.H.G. Pereira; R. Hanner; F. Foresti; C. Oliveira Can DNA barcoding accurately discriminate megadiverse Neotropical freshwater fish fauna?, BMC Genet., Volume 14 (2013), p. 20

[46] F.B. Bookstein; R. Chernoff; J. Elder; Humpries; G.R. Smith; StraussF R.E. Morphometrics in evolutionary biology, Acad. Nat. Sci. Phila. Spec. Publ., Volume 15 (1985), pp. 1-277

[47] M. González-Castro; S. Heras; M.B. Cousseau; M.I. Roldán Assessing species validity of Mugil platanus Günther, 1880 in relation to Mugil cephalus Linnaeus, 1758 (Actinopterygii), Ital. J. Zool., Volume 75 (2008), pp. 319-325

[48] M. González-Castro; A.L. Ibáñez; S. Heras; M.I. Roldán; M.B. Cousseau Assesment of lineal versus landmarks-based morphometry for discriminating species of Mugilidae (Actinopterygii), Zool. Stud., Volume 51 (2012), pp. 1515-1528

[49] J. Lleonart; J. Salat; G.J. Torres Removing allometrics effects of body size in morphological analysis, J. Theor. Biol., Volume 205 (2000), pp. 85-93

[50] A. Lombarte; J. Lleonart Otolith size changes related with body growth, habitat depth and temperature, Environ. Biol. Fish., Volume 37 (1993), pp. 297-306

[51] R. Salomón; S. Melo; N. Winzer MULTIVARIADO software, Departamento de Matemática, Universidad Nacional del Sur, Argentina, 2004

[52] D.E. Slice MORPHEUS software, Department of Ecology and Evolution, State university of New York at Stony Brook, USA, 1994

[53] R.E. Strauss; F.L. Bookstein The truss: body form reconstruccion in morphometrics, Syst. Zool., Volume 131 (1982), pp. 113-135

[54] F.L. Bookstein Morphometric tools for landmark data, Cambridge University Press, New York, 1991

[55] A.E. Valentin; X. Penin; J.-P. Chanut; J.-M. Sévigny; F.J. Rohlf Arching effect on fish body shape in geometric morphometric studies, J. Fish Biol., Volume 73 (2008), pp. 623-638

[56] F.J. Rohlf tpsRelw, version 1.46 Stony Brook, Department of Ecology and Evolution, State University of New York at Stony Brook, NY, 2008 (Available from http://www.morph.bio.sunysb.edu/morph/)

[57] F.J. Rohlf; L.F. Marcus A revolution in morphometrics, Trends Ecol. Evol., Volume 8 (1993), pp. 129-132

[58] F.J. Rohlf Shape statistics: procrustes superimpositions and tangent spaces, J. Classif., Volume 16 (1999), pp. 197-223

[59] S.X. Cadrin Advances in morphometric identification of fishery stocks, Rev. Fish. Biol. Fisher., Volume 10 (2000), pp. 91-112

[60] M.D. Moles; B.W. Robinson; T.A. Johnston; R.A. Cunjak; T.D. Jardine; J.M. Casselman et al. Morphological and trophic differentiation of growth morphotypes of walleye (Sander vitreus) from Lake Winnipeg, Canada, Can. J. Zool., Volume 88 (2010), pp. 950-960

[61] J.M. Díaz de Astarloa; S. Bezzi; M. González-Castro; E. Mabragaña; D. Hernández; S.M. Delpiani et al. Morphological, morphometric, meristic and osteological evidence for two species of hake (Actinopterygii: Gadiformes: Merluccius) in Argentinean waters, J. Fish Biol., Volume 78 (2011), pp. 1336-1358

[62] P. Orlando; M. González-Castro; E. Mabragaña New insights for the discrimination between species of Sympterygia (Rajidae) of the South West Atlantic Ocean: on the use of Geometric Morphometrics and spinulation patterns, J. Appl. Ichthyol., Volume 31 (2015), pp. 381-389

[63] K.M. O’Reilly; M.H. Horn Horn, Phenotypic variation among populations of Atherinops affinis (Atherinopsidae) with insights from a geometric morphometric analysis, J. Fish Biol., Volume 64 (2004), pp. 1117-1135

[64] B.L. Flucker; F. Pezold; R.L. Minton Molecular and morphological divergence in the inland silverside (Menidia beryllina) along a freshwater-estuarine interface, Environ. Biol. Fish., Volume 91 (2011), pp. 311-325

[65] S. Crichigno; C. Conte-Grand; M. Battini; V. Cussac Cephalic morphological variation in freshwater silversides Odontesthes hatcheri and Odontesthes bonariensis in Patagonia: introgression and ecological relationships, J. Fish Biol., Volume 83 (2013), pp. 542-559

[66] S.E. Campana; A.E. Valentin; J.M. Sevigny; D. Power Tracking seasonal migrations of redfish (Sebastes spp.) in and around the Gulf of St Lawrence using otolith elemental fingerprints, Can. J. Fish. Aquat. Sci., Volume 64 (2007), pp. 6-18

[67] R.D. Ward; B.H. Holmes An analysis of nucleotide and amino acid variability in the barcode region of cytochrome c oxidase I (cox1) in fishes, Mol. Ecol. Notes, Volume 7 (2007), pp. 899-907

[68] N. Hubert; R. Hanner; E. Holm; N.E. Mandrak; E. Taylor Identifying Canadian freshwater fishes through DNA barcodes, PloS ONE, Volume 3 (2008), p. e2490

[69] E. Mabragaña; J.M. Díaz de Astarloa; R. Hanner; J. Zhang; M. González-Castro DNA Barcoding identifies Argentine fishes from marine and brackish waters, PLoS ONE, Volume 6 (2011), p. e28655

[70] R.D. Ward; R. Hanner; P.D.N. Hebert The campaign to DNA barcode all fishes, FISH-BOL, J. Fish Biol., Volume 74 (2009), pp. 329-356

[71] J.J. Rosso; E. Mabragaña; M. González-Castro; J.M. Díaz de Astarloa DNA barcoding Neotropical fishes: news from the Pampa Plain, Argentina, Mol. Ecol. Resour., Volume 12 (2012), pp. 999-1011

[72] R.N. Bamber; P.A. Henderson Pre-adaptive plasticity in atherinids and the estuarine seat of teleost evolution, J. Fish Biol., Volume 33 (1988), pp. 17-23

[73] L.B. Beheregaray; J.A. Levy Population Genetics of the Silverside Odontesthes argentinensis (Teleostei, Atherinopsidae): evidence for Speciation in an Estuary of Southern Brazil, Copeia, Volume 2000 (2000), pp. 441-447

[74] D.D. Bloom; J.T. Weir; K.R. Piller; N.R. Lovejoy Do freshwater fishes Diversity faster then marine fishes? A test using state-dependent diversification analyses and molecular phylogenetics of new world silversides (Atherinopsidae), Evolution, Volume 67 (2013), pp. 2040-2057

[75] G. García, N. Ríos, V. Gutiérrez, J. Guerra-Varela, C. Bouza-Fernández, B. Gómez-Pardo, P. Martínez-Portela, Promiscuous Speciation with Gene Flow in Silverside Fish Genus Odontesthes (Atheriniformes, Atherinopsidae) from South Western Atlantic Ocean Basins, PLoS ONE 8 e104659.

[76] M.A. Bemvenuti Diferenciação geográfica do peixe-rei Odontesthes argentinensis (Atherinopsidae), no extremo sul do Brasil, através da morfometria multivariada, Atlántica, Volume 22 (2000), pp. 71-79


Comments - Policy