A multidisciplinary approach to describe population structure of Solea solea in the Mediterranean Sea

Investigating marine species population structure in a multidisciplinary framework can reveal signatures of potential local adaptation and the consequences for management and conservation. In this study we delineate the population structure of common sole (Solea solea) in the Mediterranean Sea using genomic and otolith data, based on single nucleotide polymorphism (SNPs) markers, otolith shape and otolith trace element composition data. We correlated SNPs with environmental and spatial variables to evaluate the impact of the selected features on the actual population structure. Specifically, we used a seascape genetics approach with redundancy (RDA) and genetic environmental association (GEA) analysis to identify loci potentially involved in local adaptation. Finally, putative functional annotation was investigated to detect genes associated with the detected patterns of neutral and adaptive genetic variation. Results from both genetic and otolith data suggested significant divergence among putative populations of common sole, confirming a clear separation between the Western and Eastern Mediterranean Sea, as well as a distinct genetic cluster corresponding to the Adriatic Sea. Evidence of fine-scale population structure in the Western Mediterranean Sea was observed at outlier loci level and further differentiation in the Adriatic. Longitude and salinity variation accounted for most of the wide and fine spatial structure. The GEA detected significant associated outlier loci potentially involved in local adaptation processes under highly structured differentiation. In the RDA both spatial distribution and environmental features could partially explain the genetic structure. Our study not only indicates that separation among Mediterranean sole population is led primarily by neutral processes because of low connectivity due to spatial segregation and limited dispersal, but it also suggests the presence of local adaptation. These results should be taken into account to support and optimize the assessment of stock units, including a review and possible redefinition of fishery management units.