Integrated pangenome and population genomics reveal selection on standing genetic variation driving fiber flax-linseed divergence

Abstract

Flax (Linum usitatissimum L.) has been domesticated for dual end uses as linseed and fiber flax, yet the genomic basis of morphotype divergence remains unclear. Here, we constructed a morphotype-resolved pangenome by integrating three newly generated near telomere-to-telomere genome assemblies with 14 previously published ones. Despite substantial variation in assembly size, driven primarily by DNA transposons, gene content was highly conserved, with little evidence for significant morphotype-specific gene presence-absence variation. Population genomic analyses of 407 accessions revealed that fiber flax had reduced nucleotide diversity, extended linkage disequilibrium, and a more compact population structure relative to linseed, consistent with stronger selection and a narrower genetic base. Genome-wide differentiation was heterogeneous and concentrated in discrete regions. Integration of FST, nucleotide diversity ratios, Tajimas D, and genome-wide association signals identified morphotype-enriched genomic blocks distributed across the genome. Many candidate regions are primarily supported by directional shifts in nucleotide diversity rather than extreme differentiation, indicating selection on standing genetic variation. Genome-wide association analyses identified 1,712 unique quantitative trait nucleotides (QTNs), with predominantly small effect sizes and strong enrichment in gene-proximal regions, consistent with a polygenic architecture. Overall, fiber flax traits tend to be controlled by fewer loci with moderate-to-large effects, whereas linseed traits exhibit a more diffuse genetic architecture. Patterns of Tajimas D further support non-classical selection dynamics, with predominantly positive values in linseed and localized negative values in fiber flax, consistent with selection on standing genetic variation. Together, our results suggest that flax morphotype divergence is driven primarily by selection on pre-existing allelic variation within a conserved gene repertoire. This study provides a comprehensive framework linking genome structure, population genomics, and trait architecture, and highlights the importance of standing genetic variation as a key resource for flax breeding and improvement.