Locus-specific view of flax domestication history

Genetic diversity analysis of a flax (Linum usitatissimum L.) global collection

Background

A sustainable breeding program requires a minimum level of germplasm diversity to provide varied options for the selection of new breeding lines. To maximize genetic gain of the North Dakota State University (NDSU) flax breeding program, we aimed to increase the genetic diversity of its parental stocks by incorporating diverse genotypes. For this purpose, we analyzed the genetic diversity, linkage disequilibrium, and population sub-structure of 350 globally-distributed flax genotypes with 6200 SNP markers.

Results

All the genotypes tested clustered into seven sub-populations (P1 to P7) based on the admixture model and the output of neighbor-joining (NJ) tree analysis and principal coordinate analysis were in line with that of structure analysis. The largest sub-population separation arose from a cluster of NDSU/American genotypes with Turkish and Asian genotypes. All sub-populations showed moderate genetic diversity (average H = 0.22 and I = 0.34). The pairwise F_st comparison revealed a great degree of divergence (F_st > 0.25) between most of the combinations. A whole collection mantel test showed significant positive correlation (r = 0.30 and p < 0.01) between genetic and geographic distances, whereas it was non-significant for all sub-populations except P4 and P5 (r = 0.251, 0.349 respectively and p < 0.05). In the entire collection, the mean linkage disequilibrium was 0.03 and it decayed to its half maximum within < 21 kb distance.

Conclusions

To maximize genetic gain, hybridization between NDSU stock (P5) and Asian individuals (P6) are potentially the best option as genetic differentiation between them is highest (F_st > 0.50). In contrast, low genetic differentiation between P5 and P2 may enhance the accumulation of favorable alleles for oil and fiber upon crossing to develop dual purpose varieties. As each sub-population consists of many genotypes, a Neighbor-Joining tree and kinship matrix assist to identify distantly related genotypes. These results also inform genotyping decisions for future association mapping studies to ensure the identification of a sufficient number of molecular markers to tag all linkage blocks.

Genomic Comparison and Population Diversity Analysis Provide Insights into the Domestication and Improvement of Flax

Flax has been cultivated for its oil and fiber for thousands of years. However, it remains unclear how the modifications of agronomic traits occurred on the genetic level during flax cultivation. In this study, we conducted genome-wide variation analyses on multiple accessions of oil-use, fiber-use, landraces, and pale flax to identify the genomic variations during flax cultivation. Our findings indicate that, during flax domestication, genes relevant to flowering, dehiscence, oil production, and plant architecture were preferentially selected. Furthermore, regardless of origins, the improvement of the modern oil-use flax preceded that of the fiber-use flax, although the dual selection on oil-use and fiber-use characteristics might have occurred in the early flax domestication. We also found that the expansion of MYB46/MYB83 genes may have contributed to the unique secondary cell wall biosynthesis in flax and the directional selections on MYB46/MYB83 may have shaped the morphological profile of the current oil-use and fiber-use flax.

•

Assemblies of genomes, including oil-use flax, fiber-use flax and pale flax

•

Comparative genomic analysis between pale flax and cultivated flax

•

Dual-selection mode on oil-use and fiber-use characteristics might be existing

•

Expansion and selection of MYB46/MYB83 may shape the morphological profile of flax

Flax latitudinal adaptation at LuTFL1 altered architecture and promoted fiber production

After domestication in the Near East around 10,000 years ago several founder crops, flax included, spread to European latitudes. On reaching northerly latitudes the architecture of domesticated flax became more suitable to fiber production over oil, with longer stems, smaller seeds and fewer axillary branches. Latitudinal adaptations in crops typically result in changes in flowering time, often involving the PEBP family of genes that also have the potential to influence plant architecture. Two PEBP family genes in the flax genome, LuTFL1 and LuTFL2, vary in wild and cultivated flax over latitudinal range with cultivated flax receiving LuTFL1 alleles from northerly wild flax populations. Compared to a background of population structure of flaxes over latitude, the LuTFL1 alleles display a level of differentiation that is consistent with selection for an allele III in the north. We demonstrate through heterologous expression in Arabidopsis thaliana that LuTFL1 is a functional homolog of TFL1 in A. thaliana capable of changing both flowering time and plant architecture. We conclude that specialized fiber flax types could have formed as a consequence of a natural adaptation of cultivated flax to higher latitudes.

Genome-Wide Association Studies for Pasmo Resistance in Flax (Linum usitatissimum L.)

Pasmo is one of the most widespread diseases threatening flax production. To identify genetic regions associated with pasmo resistance (PR), a genome-wide association study was performed on 370 accessions from the flax core collection. Evaluation of pasmo severity was performed in the field from 2012 to 2016 in Morden, MB, Canada. Genotyping-by-sequencing has identified 258,873 single nucleotide polymorphisms (SNPs) distributed on all 15 flax chromosomes. Marker-trait associations were identified using ten different statistical models. A total of 692 unique quantitative trait nucleotides (QTNs) associated with 500 putative quantitative trait loci (QTL) were detected from six phenotypic PR datasets (five individual years and average across years). Different QTNs were identified with various statistical models and from individual PR datasets, indicative of the complementation between analytical methods and/or genotype × environment interactions of the QTL effects. The single-locus models tended to identify large-effect QTNs while the multi-loci models were able to detect QTNs with smaller effects. Among the putative QTL, 67 had large effects (3-23%), were stable across all datasets and explained 32-64% of the total variation for PR in the various datasets. Forty-five of these QTL spanned 85 resistance gene analogs including a large toll interleukin receptor, nucleotide-binding site, leucine-rich repeat (TNL) type gene cluster on chromosome 8. The number of QTL with positive-effect or favorite alleles (NPQTL) in accessions was significantly correlated with PR (R 2 = 0.55), suggesting that these QTL effects are mainly additive. NPQTL was also significantly associated with morphotype (R 2 = 0.52) and major QTL with positive effect alleles were present in the fiber type accessions. The 67 large effect QTL are suited for marker-assisted selection and the 500 QTL for effective genomic prediction in PR molecular breeding.

Sequence characterization and in silico structure prediction of fatty acid desaturases in linseed varieties with differential fatty acid composition

BACKGROUND

Linseed is the richest agricultural source of α-linolenic acid (ALA), an ω-3 fatty acid (FA) that offers several nutritional benefits. In the present study, sequence characterization of six desaturase genes (SAD1, SAD2, FAD2, FAD2-2, FAD3A and FAD3B) and 3D structure prediction of their proteins from ten Indian linseed varieties differing in ALA content were performed to determine whether the nucleotide and amino acid (AA) sequence variants have any functional implications in differential accumulation of ALA or other FAs in linseed.

RESULTS

The SAD and FAD2 genes exhibited few sequence variations among the ten varieties, forming only one or two protein isoforms. In contrast, the FAD3A and FAD3B genes showed more sequence variations and three or four protein isoforms. Interestingly, the two high-ALA varieties NL260 and Padmini had the same FAD3B nucleotide and protein isoforms, which differed from all other varieties. Surprisingly, no AA changes altered the 3D structures of the desaturase proteins.

CONCLUSION

Several nucleotide and AA sequence variations in desaturase genes were observed; however, they did not alter the 3D structure of any desaturase protein and were not correlated with FA levels among the ten linseed varieties, which had different ALA contents. This suggests a complex regulatory process of biosynthesis of FAs in linseed. © 2016 Society of Chemical Industry.

Population-based resequencing revealed an ancestral winter group of cultivated flax: implication for flax domestication processes

Cultivated flax (Linum usitatissimum L.) is the earliest oil and fiber crop and its early domestication history may involve multiple events of domestication for oil, fiber, capsular indehiscence, and winter hardiness. Genetic studies have demonstrated that winter cultivated flax is closely related to oil and fiber cultivated flax and shows little relatedness to its progenitor, pale flax (L. bienne Mill.), but winter hardiness is one major characteristic of pale flax. Here, we assessed the genetic relationships of 48 Linum samples representing pale flax and four trait-specific groups of cultivated flax (dehiscent, fiber, oil, and winter) through population-based resequencing at 24 genomic regions, and revealed a winter group of cultivated flax that displayed close relatedness to the pale flax samples. Overall, the cultivated flax showed a 27% reduction of nucleotide diversity when compared with the pale flax. Recombination frequently occurred at these sampled genomic regions, but the signal of selection and bottleneck was relatively weak. These findings provide some insight into the impact and processes of flax domestication and are significant for expanding our knowledge about early flax domestication, particularly for winter hardiness.