Quantitative trait loci detection of Edwardsiella tarda resistance in Japanese flounder Paralichthys olivaceus using bulked segregant analysis

In recent years, Edwardsiella tarda has become one of the most deadly pathogens of Japanese flounder (Paralichthys olivaceus), causing serious annual losses in commercial production. In contrast to the rapid advances in the aquaculture of P. olivaceus, the study of E. tarda resistance-related markers has lagged behind, hindering the development of a disease-resistant strain. Thus, a marker-trait association analysis was initiated, combining bulked segregant analysis (BSA) and quantitative trait loci (QTL) mapping. Based on 180 microsatellite loci across all chromosomes, 106 individuals from the F1333 (♀: F0768 ×♂: F0915) (Nomenclature rule: F+year+family number) were used to detect simple sequence repeats (SSRs) and QTLs associated with E. tarda resistance. After a genomic scan, three markers (Scaffold 404-21589, Scaffold 404-21594 and Scaffold 270-13812) from the same linkage group (LG)-1 exhibited a significant difference between DNA, pooled/bulked from the resistant and susceptible groups (P <0.001). Therefore, 106 individuals were genotyped using all the SSR markers in LG1 by single marker analysis. Two different analytical models were then employed to detect SSR markers with different levels of significance in LG1, where 17 and 18 SSR markers were identified, respectively. Each model found three resistance-related QTLs by composite interval mapping (CIM). These six QTLs, designated qE1–6, explained 16.0%–89.5% of the phenotypic variance. Two of the QTLs, qE-2 and qE-4, were located at the 66.7 cM region, which was considered a major candidate region for E. tarda resistance. This study will provide valuable data for further investigations of E. tarda resistance genes and facilitate the selective breeding of disease-resistant Japanese flounder in the future.