Iterative taxonomic study of Pareiorhaphis hystrix (Siluriformes, Loricariidae) suggests a single, yet phenotypically variable, species in south Brazil

Pareiorhaphis hystrix is a widely distributed species, occurring in the upper and middle Uruguay River and in the Taquari River basin, Patos Lagoon system, southern Brazil. Morphological variation has been detected throughout the distribution of P. hystrix, and this work seeks to test the conspecific nature of populations in several occurrence areas. Specimens from six areas in the Uruguay River basin and three in the Taquari River basin were compared. Variance analysis (ANOVA) was performed for the meristic data, and Principal Component Analysis (PCA) and Linear Discriminant Analysis (LDA) were conducted for morphometric data. Molecular analyses used coI, cytb, 12S and 16S mitochondrial genes, examining nucleotide diversity, haplotype diversity, genetic distance, and delimitation of possible multiple species through the Generalized Mixed Yule Coalescent (GMYC) method. Phylogenetic relationships of studied populations were also investigated through Bayesian inference. While PCA indicated a tendency of overlap between areas, ANOVA and LDA detected a subtle differentiation between populations from the two hydrographic basins. Yet, both latter analyses recovered the population from Pelotas River, a tributary to Uruguay River, as more similar to populations from Taquari River, which is congruent to morphological observations of anterior abdominal plates. The molecular data indicated a nucleotide diversity lower than the haplotypic diversity, suggestive of recent expansion. The concatenated haplotype network points to slight differentiation between areas, with each locality presenting unique and non-shared haplotypes, although with few mutational steps in general. The species delimitation by coalescence analysis suggested the presence of a variable number of OTUs depending on the inclusion or exclusion of an outgroup. In general, the morphological data suggest a subtle variation by river basin, while the genetic data indicates a weak population structuration by hydrographic areas, especially the Chapecó and Passo Fundo rivers. However, there is still not enough differentiation between the specimens to suggest multiple species. The iterative analyses indicate that Pareiorhaphis hystrix is composed of a single, although variable, species.