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

Local climate and vernalization sensitivity predict the latitudinal patterns of flowering onset in the crop wild relative Linum bienne Mill

BACKGROUND AND AIMS

The timing of flowering onset is often correlated with latitude, indicative of climatic gradients. Flowering onset in temperate species commonly requires exposure to cold temperatures, known as vernalization. Hence, population differentiation of flowering onset with latitude might reflect adaptation to the local climatic conditions experienced by populations.

METHODS

Within its western range, seeds from Linum bienne populations (the wild relative of cultivated Linum usitatissimum) were used to describe the latitudinal differentiation of flowering onset to determine its association with the local climate of the population. A vernalization experiment including different crop cultivars was used to determine how vernalization accelerates flowering onset, in addition to the vernalization sensitivity response among populations and cultivars. Additionally, genetic differentiation of L. bienne populations along the latitudinal range was scrutinized using microsatellite markers.

KEY RESULTS

Flowering onset varied with latitude of origin, with southern populations flowering earlier than their northern counterparts. Vernalization reduced the number of days to flowering onset, but vernalization sensitivity was greater in northern populations compared with southern ones. Conversely, vernalization delayed flowering onset in the crop, exhibiting less variation in sensitivity. In L. bienne, both flowering onset and vernalization sensitivity were better predicted by the local climate of the population than by latitude itself. Microsatellite data unveiled genetic differentiation of populations, forming two groups geographically partitioned along latitude.

CONCLUSIONS

The consistent finding of latitudinal variation across experiments suggests that both flowering onset and vernalization sensitivity in L. bienne populations are under genetic regulation and might depend on climatic cues at the place of origin. The association with climatic gradients along latitude suggests that the climate experienced locally drives population differentiation of the flowering onset and vernalization sensitivity patterns. The genetic population structure suggests that past population history could have influenced the flowering initiation patterns detected, which deserves further work.

Genetic Architecture of Powdery Mildew Resistance Revealed by a Genome-Wide Association Study of a Worldwide Collection of Flax (Linum usitatissimum L.)

Powdery mildew is one of the most important diseases of flax and is particularly prejudicial to its yield and oil or fiber quality. This disease, caused by the obligate biotrophic ascomycete Oïdium lini, is progressing in France. Genetic resistance of varieties is critical for the control of this disease, but very few resistance genes have been identified so far. It is therefore necessary to identify new resistance genes to powdery mildew suitable to the local context of pathogenicity. For this purpose, we studied a worldwide diversity panel composed of 311 flax genotypes both phenotyped for resistance to powdery mildew resistance over 2 years of field trials in France and resequenced. Sequence reads were mapped on the CDC Bethune reference genome revealing 1,693,910 high-quality SNPs, further used for both population structure analysis and genome-wide association studies (GWASs). A number of four major genetic groups were identified, separating oil flax accessions from America or Europe and those from Asia or Middle-East and fiber flax accessions originating from Eastern Europe and those from Western Europe. A number of eight QTLs were detected at the false discovery rate threshold of 5%, located on chromosomes 1, 2, 4, 13, and 14. Taking advantage of the moderate linkage disequilibrium present in the flax panel, and using the available genome annotation, we identified potential candidate genes. Our study shows the existence of new resistance alleles against powdery mildew in our diversity panel, of high interest for flax breeding program.

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

QTL mapping and molecular characterization of the classical D locus controlling seed and flower color in Linum usitatissimum (flax)

The flowers of flax (linseed) are blue-hued, ephemeral and self-pollinating, and the seeds are typically brown. A century-old interest in natural yellow seed variants and a historical model point to recessive alleles in B1, D and G loci being responsible, but the functional aspects had remained unknown. Here, we characterized the "D" locus by quantitative trait loci (QTL) mapping and identified a FLAVONOID 3'5' HYDROXYLASE (F3'5'H) gene therein. It does not belong to the F3'5'H clade, but resembles biochemically characterized F3'Hs (flavonoid 3' hydroxylase) but without F3'H activity. The genome lacks other F3'H or F3'H-like genes. The apparent neo-functionalization from F3'H is associated with a Thr498 → Ser498 substitution in a substrate recognition site (SRS). The yellow seed and white flower phenotypes of the classical d mutation was found to be due to one nucleotide deletion that would truncate the deduced product and remove three of the six potential SRS, negatively impacting delphinidin synthesis. Delphinidin is sporadic in angiosperms, and flax has no known pollination syndrome(s) with functional pollinator group(s) that are attracted to blue flowers, raising questions on the acquisition of F3'5'H. The appearance of d allele is suggestive of the beginning of the loss of F3'5'H in this species.

SSR and morphological trait based population structure analysis of 130 diverse flax (Linum usitatissimum L.) accessions

A total of 130 flax accessions of diverse morphotypes and worldwide origin were assessed for genetic diversity and population structure using 11 morphological traits and microsatellite markers (15 gSSRs and 7 EST-SSRs). Analysis performed after classifying these accessions on the basis of plant height, branching pattern, seed size, Indian/foreign origin into six categories called sub-populations viz. fibre type exotic, fibre type indigenous, intermediate type exotic, intermediate type indigenous, linseed type exotic and linseed type indigenous. The study assessed different diversity indices, AMOVA, population structure and included a principal coordinate analysis based on different marker systems. The highest diversity was exhibited by gSSR markers (SI=0.46; He=0.31; P=85.11). AMOVA based on all markers explained significant difference among fibre type, intermediate type and linseed type populations of flax. In terms of variation explained by different markers, EST-SSR markers (12%) better differentiated flax populations compared to morphological (9%) and gSSR (6%) markers at P=0.01. The maximum Nei's unbiased genetic distance (D=0.11) was observed between fibre type and linseed type exotic sub-populations based on EST-SSR markers. The combined structure analysis by using all markers grouped Indian fibre type accessions (63.4%) in a separate cluster along with the Indian intermediate type (48.7%), whereas Indian accessions (82.16%) of linseed type constituted an independent cluster. These findings were supported by the results of the principal coordinate analysis. Morphological markers employed in the study found complementary with microsatellite based markers in deciphering genetic diversity and population structure of the flax germplasm.

Genetic variation of six desaturase genes in flax and their impact on fatty acid composition

Flax (Linum usitatissimum L.) is one of the richest plant sources of omega-3 fatty acids praised for their health benefits. In this study, the extent of the genetic variability of genes encoding stearoyl-ACP desaturase (SAD), and fatty acid desaturase 2 (FAD2) and 3 (FAD3) was determined by sequencing the six paralogous genes from 120 flax accessions representing a broad range of germplasm including some EMS mutant lines. A total of 6 alleles for sad1 and sad2, 21 for fad2a, 5 for fad2b, 15 for fad3a and 18 for fad3b were identified. Deduced amino acid sequences of the alleles predicted 4, 2, 3, 4, 6 and 7 isoforms, respectively. Allele frequencies varied greatly across genes. Fad3a, with 110 SNPs and 19 indels, and fad3b, with 50 SNPs and 5 indels, showed the highest levels of genetic variations. While most of the SNPs and all the indels were silent mutations, both genes carried nonsense SNP mutations resulting in premature stop codons, a feature not observed in sad and fad2 genes. Some alleles and isoforms discovered in induced mutant lines were absent in the natural germplasm. Correlation of these genotypic data with fatty acid composition data of 120 flax accessions phenotyped in six field experiments revealed statistically significant effects of some of the SAD and FAD isoforms on fatty acid composition, oil content and iodine value. The novel allelic variants and isoforms identified for the six desaturases will be a resource for the development of oilseed flax with unique and useful fatty acid profiles.

Genetic variation of six desaturase genes in flax and their impact on fatty acid composition

Flax (Linum usitatissimum L.) is one of the richest plant sources of omega-3 fatty acids praised for their health benefits. In this study, the extent of the genetic variability of genes encoding stearoyl-ACP desaturase (SAD), and fatty acid desaturase 2 (FAD2) and 3 (FAD3) was determined by sequencing the six paralogous genes from 120 flax accessions representing a broad range of germplasm including some EMS mutant lines. A total of 6 alleles for sad1 and sad2, 21 for fad2a, 5 for fad2b, 15 for fad3a and 18 for fad3b were identified. Deduced amino acid sequences of the alleles predicted 4, 2, 3, 4, 6 and 7 isoforms, respectively. Allele frequencies varied greatly across genes. Fad3a, with 110 SNPs and 19 indels, and fad3b, with 50 SNPs and 5 indels, showed the highest levels of genetic variations. While most of the SNPs and all the indels were silent mutations, both genes carried nonsense SNP mutations resulting in premature stop codons, a feature not observed in sad and fad2 genes. Some alleles and isoforms discovered in induced mutant lines were absent in the natural germplasm. Correlation of these genotypic data with fatty acid composition data of 120 flax accessions phenotyped in six field experiments revealed statistically significant effects of some of the SAD and FAD isoforms on fatty acid composition, oil content and iodine value. The novel allelic variants and isoforms identified for the six desaturases will be a resource for the development of oilseed flax with unique and useful fatty acid profiles.