Population Genetic Study on the European Flounder (<i>Platichthys flesus</i>) from the Southern Baltic Sea Using SNPs and Microsatellite Markers

The European flounder (<i>Platichthys flesus</i>), which is closely related to the recently discovered Baltic flounder (<i>Platichthys solemdali</i>), is currently the third most commercially fished species in the Baltic Sea. According to the available data from the Polish Fisheries Monitoring Center and fishermen’s observations, the body condition indices of the species in the Baltic Sea have declined in recent years. The aim of the present study was to obtain information on the current patterns of genetic variability and the population structure of the European flounder and to verify whether the Baltic flounder is present in the southern Baltic Sea. Moreover, we aimed to verify whether the observed decline in the body condition indices of the species in the Baltic Sea might be associated with adaptive alterations in its gene pool due to increased fishing pressure. For this purpose, 190 fish were collected from four locations along the central coastline of Poland, i.e., Mechelinki, Władysławowo, the Vistula Lagoon in 2018, and the Słupsk Bank in 2020. The fish were morphologically analyzed and then genetically screened by the application of nineteen microsatellite DNA and two diagnostic SNP markers. The examined European flounder specimens displayed a high level of genetic diversity (<i>PIC</i> = 0.832–0.903, <i>I</i> = 2.579–2.768). A lack of significant genetic differentiation (<i>Fst</i> = 0.004, <i>p</i> > 0.05) was observed in all the examined fish, indicating that the European flounder in the sampled area constitutes a single genetic cluster. A significant deficiency in heterozygotes (<i>Fis</i> = 0.093, <i>p</i> < 0.05) and overall deviations from Hardy–Weinberg expectations (H-WE) were only detected in fish sampled from the Słupsk Bank. The estimated effective population size (<i>Ne</i>) among the sampled fish groups varied from 712 (Słupsk Bank) to 10,115 (Władysławowo and Mechelinki). However, the recorded values of the Garza–Williamson indicator (<i>M</i> = 0.574–0.600) and the lack of significant (<i>p</i> > 0.05) differences in <i>Heq</i> > <i>He</i> under the <i>SMM</i> model did not support the species’ population size changes in the past. The applied SNP markers did not detect the presence of the Baltic flounder among the fish sampled from the studied area. The analysis of an association between biological traits and patterns of genetic diversity did not detect any signs of directional selection or density-dependent adaptive changes in the gene pool of the examined fish that might be caused by increased fishing pressure.