Chromosome-scale pseudomolecules refined by optical, physical and genetic maps in flax

An R2R3-MYB Transcriptional Factor LuMYB314 Associated with the Loss of Petal Pigmentation in Flax (Linum usitatissimum L.)

The loss of anthocyanin pigments is one of the most common evolutionary transitions in petal color, yet the genetic basis for these changes in flax remains largely unknown. In this study, we used crossing studies, a bulk segregant analysis, genome-wide association studies, a phylogenetic analysis, and transgenic testing to identify genes responsible for the transition from blue to white petals in flax. This study found no correspondence between the petal color and seed color, refuting the conclusion that a locus controlling the seed coat color is associated with the petal color, as reported in previous studies. The locus controlling the petal color was mapped using a BSA-seq analysis based on the F2 population. However, no significantly associated genomic regions were detected. Our genome-wide association study identified a highly significant QTL (BP4.1) on chromosome 4 associated with flax petal color in the natural population. The combination of a local Manhattan plot and an LD heat map identified LuMYB314, an R2R3-MYB transcription factor, as a potential gene responsible for the natural variations in petal color in flax. The overexpression of LuMYB314 in both Arabidopsis thaliana and Nicotiana tabacum resulted in anthocyanin deposition, indicating that LuMYB314 is a credible candidate gene for controlling the petal color in flax. Additionally, our study highlights the limitations of the BSA-seq method in low-linkage genomic regions, while also demonstrating the powerful detection capabilities of GWAS based on high-density genomic variation mapping. This study enhances our genetic insight into petal color variations and has potential breeding value for engineering LuMYB314 to develop colored petals, bast fibers, and seeds for multifunctional use in flax.

Genomic prediction for agronomic traits in a diverse Flax (Linum usitatissimum L.) germplasm collection

Breeding programs require exhaustive phenotyping of germplasms, which is time-demanding and expensive. Genomic prediction helps breeders harness the diversity of any collection to bypass phenotyping. Here, we examined the genomic prediction's potential for seed yield and nine agronomic traits using 26,171 single nucleotide polymorphism (SNP) markers in a set of 337 flax (Linum usitatissimum L.) germplasm, phenotyped in five environments. We evaluated 14 prediction models and several factors affecting predictive ability based on cross-validation schemes. Models yielded significant variation among predictive ability values across traits for the whole marker set. The ridge regression (RR) model covering additive gene action yielded better predictive ability for most of the traits, whereas it was higher for low heritable traits by models capturing epistatic gene action. Marker subsets based on linkage disequilibrium decay distance gave significantly higher predictive abilities to the whole marker set, but for randomly selected markers, it reached a plateau above 3000 markers. Markers having significant association with traits improved predictive abilities compared to the whole marker set when marker selection was made on the whole population instead of the training set indicating a clear overfitting. The correction for population structure did not increase predictive abilities compared to the whole collection. However, stratified sampling by picking representative genotypes from each cluster improved predictive abilities. The indirect predictive ability for a trait was proportionate to its correlation with other traits. These results will help breeders to select the best models, optimum marker set, and suitable genotype set to perform an indirect selection for quantitative traits in this diverse flax germplasm collection.

Whole-genome sequencing in medicinal plants: current progress and prospect

Advancements in genomics have dramatically accelerated the research on medicinal plants, and the development of herbgenomics has promoted the "Project of 1K Medicinal Plant Genome" to decipher their genetic code. However, it is difficult to obtain their high-quality whole genomes because of the prevalence of polyploidy and/or high genomic heterozygosity. Whole genomes of 123 medicinal plants were published until September 2022. These published genome sequences were investigated in this review, covering their classification, research teams, ploidy, medicinal functions, and sequencing strategies. More than 1,000 institutes or universities around the world and 50 countries are conducting research on medicinal plant genomes. Diploid species account for a majority of sequenced medicinal plants. The whole genomes of plants in the Poaceae family are the most studied. Almost 40% of the published papers studied species with tonifying, replenishing, and heat-cleaning medicinal effects. Medicinal plants are still in the process of domestication as compared with crops, thereby resulting in unclear genetic backgrounds and the lack of pure lines, thus making their genomes more difficult to complete. In addition, there is still no clear routine framework for a medicinal plant to obtain a high-quality whole genome. Herein, a clear and complete strategy has been originally proposed for creating a high-quality whole genome of medicinal plants. Moreover, whole genome-based biological studies of medicinal plants, including breeding and biosynthesis, were reviewed. We also advocate that a research platform of model medicinal plants should be established to promote the genomics research of medicinal plants.

Fine-mapping of a major locus for Fusarium wilt resistance in flax (Linum usitatissimum L.)

KEY MESSAGE

Fine-mapping of a locus on chromosome 1 of flax identified an S-lectin receptor-like kinase (SRLK) as the most likely candidate for a major Fusarium wilt resistance gene. Fusarium wilt, caused by the soil-borne fungal pathogen Fusarium oxysporum f. sp. lini, is a devastating disease in flax. Genetic resistance can counteract this disease and limit its spread. To map major genes for Fusarium wilt resistance, a recombinant inbred line population of more than 700 individuals derived from a cross between resistant cultivar 'Bison' and susceptible cultivar 'Novelty' was phenotyped in Fusarium wilt nurseries at two sites for two and three years, respectively. The population was genotyped with 4487 single nucleotide polymorphism (SNP) markers. Twenty-four QTLs were identified with IciMapping, 18 quantitative trait nucleotides with 3VmrMLM and 108 linkage disequilibrium blocks with RTM-GWAS. All models identified a major QTL on chromosome 1 that explained 20-48% of the genetic variance for Fusarium wilt resistance. The locus was estimated to span ~ 867 Kb but included a ~ 400 Kb unresolved region. Whole-genome sequencing of 'CDC Bethune', 'Bison' and 'Novelty' produced ~ 450 Kb continuous sequences of the locus. Annotation revealed 110 genes, of which six were considered candidate genes. Fine-mapping with 12 SNPs and 15 Kompetitive allele-specific PCR (KASP) markers narrowed down the interval to ~ 69 Kb, which comprised the candidate genes Lus10025882 and Lus10025891. The latter, a G-type S-lectin receptor-like kinase (SRLK) is the most likely resistance gene because it is the only polymorphic one. In addition, Fusarium wilt resistance genes previously isolated in tomato and Arabidopsis belonged to the SRLK class. The robust KASP markers can be used in marker-assisted breeding to select for this major Fusarium wilt resistance locus.

The Genetic Dissection of Nitrogen Use-Related Traits in Flax (Linum usitatissimum L.) at the Seedling Stage through the Integration of Multi-Locus GWAS, RNA-seq and Genomic Selection

Nitrogen (N), the most important macro-nutrient for plant growth and development, is a key factor that determines crop yield. Yet its excessive applications pollute the environment and are expensive. Hence, studying nitrogen use efficiency (NUE) in crops is fundamental for sustainable agriculture. Here, an association panel consisting of 123 flax accessions was evaluated for 21 NUE-related traits at the seedling stage under optimum N (N+) and N deficiency (N-) treatments to dissect the genetic architecture of NUE-related traits using a multi-omics approach integrating genome-wide association studies (GWAS), transcriptome analysis and genomic selection (GS). Root traits exhibited significant and positive correlations with NUE under N- conditions (r = 0.33 to 0.43, p < 0.05). A total of 359 QTLs were identified, accounting for 0.11% to 23.1% of the phenotypic variation in NUE-related traits. Transcriptomic analysis identified 1034 differentially expressed genes (DEGs) under contrasting N conditions. DEGs involved in N metabolism, root development, amino acid transport and catabolism and others, were found near the QTLs. GS models to predict NUE stress tolerance index (NUE_STI) trait were tested using a random genome-wide SNP dataset and a GWAS-derived QTLs dataset. The latter produced superior prediction accuracy (r = 0.62 to 0.79) compared to the genome-wide SNP marker dataset (r = 0.11) for NUE_STI. Our results provide insights into the QTL architecture of NUE-related traits, identify candidate genes for further studies, and propose genomic breeding tools to achieve superior NUE in flax under low N input.

Genome-wide association studies using multi-models and multi-SNP datasets provide new insights into pasmo resistance in flax

Introduction

Flax (Linum usitatissimum L.) is an economically important crop due to its oil and fiber. However, it is prone to various diseases, including pasmo caused by the fungus Septoria linicola.

Methods

In this study, we conducted field evaluations of 445 flax accessions over a five-year period (2012-2016) to assess their resistance to pasmo A total of 246,035 single nucleotide polymorphisms (SNPs) were used for genetic analysis. Four statistical models, including the single-locus model GEMMA and the multi-locus models FarmCPU, mrMLM, and 3VmrMLM, were assessed to identify quantitative trait nucleotides (QTNs) associated with pasmo resistance.

Results

We identified 372 significant QTNs or 132 tag QTNs associated with pasmo resistance from five pasmo resistance datasets (PAS2012-PAS2016 and the 5-year average, namely PASmean) and three genotypic datasets (the all SNPs/ALL, the gene-based SNPs/GB and the RGA-based SNPs/RGAB). The tag QTNs had R2 values of 0.66-16.98% from the ALL SNP dataset, 0.68-20.54%from the GB SNP dataset, and 0.52-22.42% from the RGAB SNP dataset. Of these tag QTNs, 93 were novel. Additionally, 37 resistance gene analogs (RGAs)co-localizing with 39 tag QTNs were considered as potential candidates for controlling pasmo resistance in flax and 50 QTN-by-environment interactions(QEIs) were identified to account for genes by environmental interactions. Nine RGAs were predicted as candidate genes for ten QEIs.

Discussion

Our results suggest that pasmo resistance in flax is polygenic and potentially influenced by environmental factors. The identified QTNs provide potential targets for improving pasmo resistance in flax breeding programs. This study sheds light on the genetic basis of pasmo resistance and highlights the importance of considering both genetic and environmental factors in breeding programs for flax.

Key FAD2, FAD3, and SAD Genes Involved in the Fatty Acid Synthesis in Flax Identified Based on Genomic and Transcriptomic Data

FAD (fatty acid desaturase) and SAD (stearoyl-ACP desaturase) genes play key roles in the synthesis of fatty acids (FA) and determination of oil composition in flax (Linum usitatissimum L.). We searched for FAD and SAD genes in the most widely used flax genome of the variety CDC Bethune and three available long-read assembled flax genomes-YY5, 3896, and Atlant. We identified fifteen FAD2, six FAD3, and four SAD genes. Of all the identified genes, 24 were present in duplicated pairs. In most cases, two genes from a pair differed by a significant number of gene-specific SNPs (single nucleotide polymorphisms) or even InDels (insertions/deletions), except for FAD2a-1 and FAD2a-2, where only seven SNPs distinguished these genes. Errors were detected in the FAD2a-1, FAD2a-2, FAD3c-1, and FAD3d-2 sequences in the CDC Bethune genome assembly but not in the long-read genome assemblies. Expression analysis of the available transcriptomic data for different flax organs/tissues revealed that FAD2a-1, FAD2a-2, FAD3a, FAD3b, SAD3-1, and SAD3-2 were specifically expressed in embryos/seeds/capsules and could play a crucial role in the synthesis of FA in flax seeds. In contrast, FAD2b-1, FAD2b-2, SAD2-1, and SAD2-2 were highly expressed in all analyzed organs/tissues and could be involved in FA synthesis in whole flax plants. FAD2c-2, FAD2d-1, FAD3c-1, FAD3c-2, FAD3d-1, FAD3d-2, SAD3-1, and SAD3-2 showed differential expression under stress conditions-Fusarium oxysporum infection and drought. The obtained results are essential for research on molecular mechanisms of fatty acid synthesis, FAD and SAD editing, and marker-assisted and genomic selection for breeding flax varieties with a determined fatty acid composition of oil.

Development of genomic resources for Rhodes grass (Chloris gayana), draft genome and annotated variant discovery

Genomic resources for grasses, especially warm-season grasses are limited despite their commercial and environmental importance. Here, we report the first annotated draft whole genome sequence for diploid Rhodes grass (Chloris gayana), a tropical C4 species. Generated using long read nanopore sequencing and assembled using the Flye software package, the assembled genome is 603 Mbp in size and comprises 5,233 fragments that were annotated using the GenSas pipeline. The annotated genome has 46,087 predicted genes corresponding to 92.0% of the expected genomic content present via BUSCO analysis. Gene ontology terms and repetitive elements are identified and discussed. An additional 94 individual plant genotypes originating from three diploid and two tetraploid Rhodes grass cultivars were short-read whole genome resequenced (WGR) to generate a single nucleotide polymorphism (SNP) resource for the species that can be used to elucidate inter- and intra-cultivar relationships across both ploidy levels. A total of 75,777 high quality SNPs were used to generate a phylogenetic tree, highlighting the diversity present within the cultivars which agreed with the known breeding history. Differentiation was observed between diploid and tetraploid cultivars. The WGR data were also used to provide insights into the nature and evolution of the tetraploid status of the species, with results largely agreeing with the published literature that the tetraploids are autotetraploid.

High-Quality Genome Assembly and Genome-Wide Association Study of Male Sterility Provide Resources for Flax Improvement

Flax is an economic crop with a long history. It is grown worldwide and is mainly used for edible oil, industry, and textiles. Here, we reported a high-quality genome assembly for "Neiya No. 9", a popular variety widely grown in China. Combining PacBio long reads, Hi-C sequencing, and a genetic map reported previously, a genome assembly of 473.55 Mb was constructed, which covers ~94.7% of the flax genome. These sequences were anchored onto 15 chromosomes. The N50 lengths of the contig and scaffold were 0.91 Mb and 31.72 Mb, respectively. A total of 32,786 protein-coding genes were annotated, and 95.9% of complete BUSCOs were found. Through morphological and cytological observation, the male sterility of flax was considered dominant nuclear sterility. Through GWAS analysis, the gene LUSG00017705 (cysteine synthase gene) was found to be closest to the most significant SNP, and the expression level of this gene was significantly lower in male sterile plants than in fertile plants. Among the significant SNPs identified in the GWAS analysis, only two were located in the coding region, and these two SNPs caused changes in the protein encoded by LUSG00017565 (cysteine protease gene). It was speculated that these two genes may be related to male sterility in flax. This is the first time the molecular mechanism of male sterility in flax has been reported. The high-quality genome assembly and the male sterility genes revealed, provided a solid foundation for flax breeding.

Identification of new potential molecular actors related to fiber quality in flax through Omics

One of the biggest challenges for a more widespread utilization of plant fibers is to better understand the different molecular factors underlying the variability in fineness and mechanical properties of both elementary and scutched fibers. Accordingly, we analyzed genome-wide transcription profiling from bast fiber bearing tissues of seven different flax varieties (4 spring, 2 winter fiber varieties and 1 winter linseed) and identified 1041 differentially expressed genes between varieties, of which 97 were related to cell wall metabolism. KEGG analysis highlighted a number of different enriched pathways. Subsequent statistical analysis using Partial Least-Squares Discriminant Analysis showed that 73% of the total variance was explained by the first 3 X-variates corresponding to 56 differentially expressed genes. Calculation of Pearson correlations identified 5 genes showing a strong correlation between expression and morphometric data. Two-dimensional gel proteomic analysis on the two varieties showing the most discriminant and significant differences in morphometrics revealed 1490 protein spots of which 108 showed significant differential abundance. Mass spectrometry analysis successfully identified 46 proteins representing 32 non-redundant proteins. Statistical clusterization based on the expression level of genes corresponding to the 32 proteins showed clear discrimination into three separate clusters, reflecting the variety type (spring-/winter-fiber/oil). Four of the 32 proteins were also highly correlated with morphometric features. Examination of predicted functions for the 9 (5 + 4) identified genes highlighted lipid metabolism and senescence process. Calculation of Pearson correlation coefficients between expression data and retted fiber mechanical measurements (strength and maximum force) identified 3 significantly correlated genes. The genes were predicted to be connected to cell wall dynamics, either directly (Expansin-like protein), or indirectly (NAD(P)-binding Rossmann-fold superfamily protein). Taken together, our results have allowed the identification of molecular actors potentially associated with the determination of both in-planta fiber morphometrics, as well as ex-planta fiber mechanical properties, both of which are key parameters for elementary fiber and scutched fiber quality in flax.

Genetic dissection of thousand-seed weight in linseed (Linum usitatissimum L.) using multi-locus genome-wide association study

Flaxseed/linseed is an important oilseed crop having applications in the food, nutraceutical, and paint industry. Seed weight is one of the most crucial determinants of seed yield in linseed. Here, quantitative trait nucleotides (QTNs) associated with thousand-seed weight (TSW) have been identified using multi-locus genome-wide association study (ML-GWAS). Field evaluation was carried out in five environments in multi-year-location trials. SNP genotyping information of the AM panel of 131 accessions comprising 68,925 SNPs was employed for ML-GWAS. From the six ML-GWAS methods employed, five methods helped identify a total of 84 unique significant QTNs for TSW. QTNs identified in ≥ 2 methods/environments were designated as stable QTNs. Accordingly, 30 stable QTNs have been identified for TSW accounting up to 38.65% trait variation. Alleles with positive effect on trait were analyzed for 12 strong QTNs with r ² ≥ 10.00%, which showed significant association of specific alleles with higher trait value in three or more environments. A total of 23 candidate genes have been identified for TSW, which included B3 domain-containing transcription factor, SUMO-activating enzyme, protein SCARECROW, shaggy-related protein kinase/BIN2, ANTIAUXIN-RESISTANT 3, RING-type E3 ubiquitin transferase E4, auxin response factors, WRKY transcription factor, and CBS domain-containing protein. In silico expression analysis of candidate genes was performed to validate their possible role in different stages of seed development process. The results from this study provide significant insight and elevate our understanding on genetic architecture of TSW trait in linseed.

Diversity of Linum genetic resources in global genebanks: from agro-morphological characterisation to novel genomic technologies - a review

Linseed or flaxseed is a well-recognized nutritional food with nutraceutical properties owing to high omega-3 fatty acid (α-Linolenic acid), dietary fiber, quality protein, and lignan content. Currently, linseed enjoys the status of a 'superfood' and its integration in the food chain as a functional food is evolving continuously as seed constituents are associated with lowering the risk of chronic ailments, such as heart diseases, cancer, diabetes, and rheumatoid arthritis. This crop also receives much attention in the handloom and textile sectors as the world's coolest fabric linen is made up of its stem fibers which are endowed with unique qualities such as luster, tensile strength, density, bio-degradability, and non-hazardous nature. Worldwide, major linseed growing areas are facing erratic rainfall and temperature patterns affecting flax yield, quality, and response to biotic stresses. Amid such changing climatic regimes and associated future threats, diverse linseed genetic resources would be crucial for developing cultivars with a broad genetic base for sustainable production. Furthermore, linseed is grown across the world in varied agro-climatic conditions; therefore it is vital to develop niche-specific cultivars to cater to diverse needs and keep pace with rising demands globally. Linseed genetic diversity conserved in global genebanks in the form of germplasm collection from natural diversity rich areas is expected to harbor genetic variants and thus form crucial resources for breeding tailored crops to specific culinary and industrial uses. Global genebank collections thus potentially play an important role in supporting sustainable agriculture and food security. Currently, approximately 61,000 germplasm accessions of linseed including 1,127 wild accessions are conserved in genebanks/institutes worldwide. This review analyzes the current status of Linum genetic resources in global genebanks, evaluation for agro-morphological traits, stress tolerance, and nutritional profiling to promote their effective use for sustainable production and nutrition enhancement in our modern diets.

Development and Complex Application of Methods for the Identification of Mutations in the FAD3A and FAD3B Genes Resulting in the Reduced Content of Linolenic Acid in Flax Oil

Flax is grown worldwide for seed and fiber production. Linseed varieties differ in their oil composition and are used in pharmaceutical, food, feed, and industrial production. The field of application primarily depends on the content of linolenic (LIN) and linoleic (LIO) fatty acids. Inactivating mutations in the FAD3A and FAD3B genes lead to a decrease in the LIN content and an increase in the LIO content. For the identification of the three most common low-LIN mutations in flax varieties (G-to-A in exon 1 of FAD3A substituting tryptophan with a stop codon, C-to-T in exon 5 of FAD3A leading to arginine to a stop codon substitution, and C-to-T in exon 2 of FAD3B resulting in histidine to tyrosine substitution), three approaches were proposed: (1) targeted deep sequencing, (2) high resolution melting (HRM) analysis, (3) cleaved amplified polymorphic sequences (CAPS) markers. They were tested on more than a thousand flax samples of various types and showed promising results. The proposed approaches can be used in marker-assisted selection to choose parent pairs for crosses, separate heterogeneous varieties into biotypes, and select genotypes with desired homozygous alleles of the FAD3A and FAD3B genes at the early stages of breeding for the effective development of varieties with a particular LIN and LIO content, as well as in basic studies of the molecular mechanisms of fatty acid synthesis in flax seeds to select genotypes adequate to the tasks.

The flax genome reveals orbitide diversity

Background

Ribosomally-synthesized cyclic peptides are widely found in plants and exhibit useful bioactivities for humans. The identification of cyclic peptide sequences and their precursor proteins is facilitated by the growing number of sequenced genomes. While previous research largely focused on the chemical diversity of these peptides across various species, there is little attention to a broader range of potential peptides that are not chemically identified.

Results

A pioneering study was initiated to explore the genetic diversity of linusorbs, a group of cyclic peptides uniquely occurring in cultivated flax (Linum usitatissimum). Phylogenetic analysis clustered the 5 known linusorb precursor proteins into two clades and one singleton. Preliminary tBLASTn search of the published flax genome using the whole protein sequence as query could only retrieve its homologues within the same clade. This limitation was overcome using a profile-based mining strategy. After genome reannotation, a hidden Markov Model (HMM)-based approach identified 58 repeats homologous to the linusorb-embedded repeats in 8 novel proteins, implying that they share common ancestry with the linusorb-embedded repeats. Subsequently, we developed a customized profile composed of a random linusorb-like domain (LLD) flanked by 5 conserved sites and used it for string search of the proteome, which extracted 281 LLD-containing repeats (LLDRs) in 25 proteins. Comparative analysis of different repeat categories suggested that the 5 conserved flanking sites among the non-homologous repeats have undergone convergent evolution driven by functional selection.

Conclusions

The profile-based mining approach is suitable for analyzing repetitive sequences. The 25 LLDR proteins identified herein represent the potential diversity of cyclic peptides within the flax genome and lay a foundation for further studies on the functions and evolution of these protein tandem repeats.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08735-x.

Isolating Linum usitatissimum L. Nuclear DNA Enabled Assembling High-Quality Genome

High-quality genome sequences help to elucidate the genetic basis of numerous biological processes and track species evolution. For flax (Linum usitatissimum L.)—a multifunctional crop, high-quality assemblies from Oxford Nanopore Technologies (ONT) data were unavailable, largely due to the difficulty of isolating pure high-molecular-weight DNA. This article proposes a scheme for gaining a contiguous L. usitatissimum assembly using Nanopore data. We developed a protocol for flax nuclei isolation with subsequent DNA extraction, which allows obtaining about 5 μg of pure high-molecular-weight DNA from 0.5 g of leaves. Such an amount of material can be collected even from a single plant and yields more than 30 Gb of ONT data in two MinION runs. We performed a comparative analysis of different genome assemblers and polishers on the gained data and obtained the final 447.1-Mb assembly of L. usitatissimum line 3896 genome using the Canu—Racon (two iterations)—Medaka combination. The genome comprised 1695 contigs and had an N50 of 6.2 Mb and a completeness of 93.8% of BUSCOs from eudicots_odb10. Our study highlights the impact of the chosen genome construction strategy on the resulting assembly parameters and its eligibility for future genomic studies.

Adaptive Response and Transcriptomic Analysis of Flax (Linum usitatissimum L.) Seedlings to Salt Stress

Soil salinity constrains agricultural development in arid regions. Flax is an economically important crop in many countries, and screening or breeding salinity-resistant flax cultivars is necessary. Based on the previous screening of flaxseed cultivars C71 (salt-sensitive) and C116 (salt-tolerant) as test materials, flax seedlings stressed with different concentrations of NaCl (0, 100, 150, 200, and 250 mmol/L) for 21 days were used to investigate the effects of salt stress on the growth characteristics, osmotic regulators, and antioxidant capacity of these flax seedlings and to reveal the adaptive responses of flax seedlings to salt stress. The results showed that plant height and root length of flax were inhibited, with C116 showing lower growth than C71. The concentrations of osmotic adjustment substances such as soluble sugars, soluble proteins, and proline were higher in the resistant material, C116, than in the sensitive material, C71, under different concentrations of salt stress. Consistently, C116 showed a better rapid scavenging ability for reactive oxygen species (ROS) and maintained higher activities of antioxidant enzymes to balance salt injury stress by inhibiting growth under salt stress. A transcriptome analysis of flax revealed that genes related to defense and senescence were significantly upregulated, and genes related to the growth and development processes were significantly downregulated under salt stress. Our results indicated that one of the important adaptations to tolerance to high salt stress is complex physiological remediation by rapidly promoting transcriptional regulation in flax.

Integrated omics approaches for flax improvement under abiotic and biotic stress: Current status and future prospects

Flax (Linum usitatissimum L.) or linseed is one of the important industrial crops grown all over the world for seed oil and fiber. Besides oil and fiber, flax offers a wide range of nutritional and therapeutic applications as a feed and food source owing to high amount of α-linolenic acid (omega-3 fatty acid), lignans, protein, minerals, and vitamins. Periodic losses caused by unpredictable environmental stresses such as drought, heat, salinity-alkalinity, and diseases pose a threat to meet the rising market demand. Furthermore, these abiotic and biotic stressors have a negative impact on biological diversity and quality of oil/fiber. Therefore, understanding the interaction of genetic and environmental factors in stress tolerance mechanism and identification of underlying genes for economically important traits is critical for flax improvement and sustainability. In recent technological era, numerous omics techniques such as genomics, transcriptomics, metabolomics, proteomics, phenomics, and ionomics have evolved. The advancements in sequencing technologies accelerated development of genomic resources which facilitated finer genetic mapping, quantitative trait loci (QTL) mapping, genome-wide association studies (GWAS), and genomic selection in major cereal and oilseed crops including flax. Extensive studies in the area of genomics and transcriptomics have been conducted post flax genome sequencing. Interestingly, research has been focused more for abiotic stresses tolerance compared to disease resistance in flax through transcriptomics, while the other areas of omics such as metabolomics, proteomics, ionomics, and phenomics are in the initial stages in flax and several key questions remain unanswered. Little has been explored in the integration of omic-scale data to explain complex genetic, physiological and biochemical basis of stress tolerance in flax. In this review, the current status of various omics approaches for elucidation of molecular pathways underlying abiotic and biotic stress tolerance in flax have been presented and the importance of integrated omics technologies in future research and breeding have been emphasized to ensure sustainable yield in challenging environments.

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.

Identification of QTNs Associated With Flowering Time, Maturity, and Plant Height Traits in Linum usitatissimum L. Using Genome-Wide Association Study

Early flowering, maturity, and plant height are important traits for linseed to fit in rice fallows, for rainfed agriculture, and for economically viable cultivation. Here, Multi-Locus Genome-Wide Association Study (ML-GWAS) was undertaken in an association mapping panel of 131 accessions, genotyped using 68,925 SNPs identified by genotyping by sequencing approach. Phenotypic evaluation data of five environments comprising 3 years and two locations were used. GWAS was performed for three flowering time traits including days to 5%, 50%, and 95% flowering, days to maturity, and plant height by employing five ML-GWAS methods: FASTmrEMMA, FASTmrMLM, ISIS EM-BLASSO, mrMLM, and pLARmEB. A total of 335 unique QTNs have been identified for five traits across five environments. 109 QTNs were stable as observed in ≥2 methods and/or environments, explaining up to 36.6% phenotypic variance. For three flowering time traits, days to maturity, and plant height, 53, 30, and 27 stable QTNs, respectively, were identified. Candidate genes having roles in flower, pollen, embryo, seed and fruit development, and xylem/phloem histogenesis have been identified. Gene expression of candidate genes for flowering and plant height were studied using transcriptome of an early maturing variety Sharda (IC0523807). The present study unravels QTNs/candidate genes underlying complex flowering, days to maturity, and plant height traits in linseed.

DNA hypomethylation is the plausible driver of heat stress adaptation in Linum usitatissimum

Heat stress has a significant impact on the climatic adaptation of flax, a cool-season economic crop. Genome-wide DNA methylation patterns are crucial for understanding how flax cultivars respond to heat adversities. It is worth noting that the DNA methylome in flax has yet to be investigated at the nucleotide level. Although heat stress above 40°C caused oxidative damage in flax leaves, 5-azacytidine, a hypomethylating agent, reduced this effect by 15%-24%. Differences in the expression of the LuMET1 (DNA methyltransferase) gene suggested that DNA methylation/demethylation may play a major role in the flax heat stress response. Thus, whole-genome bisulfite sequencing-derived DNA methylation profiles in flax, with or without heat stress and 5-azaC, were developed and analyzed here. In response to heat stress, a high percentage of significant differentially methylated regions (DMRs), particularly hypomethylated DMRs, were identified in the CHH nucleotide sequence context (H = A/T/C). Some of these DMRs overlapped with transposable element insertions. The majority of DMRs were discovered in intergenic regions, but several DMR loci were also found near genes relevant to heat stress response and epigenetic processes. These DMRs, in particular, are linked to CpG islands, implying a possible role in promoter methylation and gene silencing. The DMRs discovered in this study are crucial for understanding and identifying the key players in heat stress response in flax, which will help in developing climate-smart flax varieties.

Insights into the Genetic Architecture and Genomic Prediction of Powdery Mildew Resistance in Flax (Linum usitatissimum L.)

Powdery mildew (PM), caused by the fungus Oidium lini in flax, can cause defoliation and reduce seed yield and quality. To date, one major dominant gene (Pm1) and three quantitative trait loci (QTL) on chromosomes 1, 7 and 9 have been reported for PM resistance. To fully dissect the genetic architecture of PM resistance and identify QTL, a diverse flax core collection of 372 accessions augmented with an additional 75 breeding lines were sequenced, and PM resistance was evaluated in the field for eight years (2010–2017) in Morden, Manitoba, Canada. Genome-wide association studies (GWAS) were performed using two single-locus and seven multi-locus statistical models with 247,160 single nucleotide polymorphisms (SNPs) and the phenotypes of the 447 individuals for each year separately as well as the means over years. A total of 349 quantitative trait nucleotides (QTNs) were identified, of which 44 large-effect QTNs (R² = 10–30%) were highly stable over years. The total number of favourable alleles per accession was significantly correlated with PM resistance (r = 0.74), and genomic selection (GS) models using all identified QTNs generated significantly higher predictive ability (r = 0.93) than those constructed using the 247,160 genome-wide random SNP (r = 0.69), validating the overall reliability of the QTNs and showing the additivity of PM resistance in flax. The QTNs were clustered on the distal ends of all 15 chromosomes, especially on chromosome 5 (0.4–5.6 Mb and 9.4–16.9 Mb) and 13 (4.7–5.2 Mb). To identify candidate genes, a dataset of 3230 SNPs located in resistance gene analogues (RGAs) was used as input for GWAS, from which an additional 39 RGA-specific QTNs were identified. Overall, 269 QTN loci harboured 445 RGAs within the 200 Kb regions spanning the QTNs, including 45 QTNs located within the RGAs. These RGAs supported by significant QTN/SNP allele effects were mostly nucleotide binding site and leucine-rich repeat receptors (NLRs) belonging to either coiled-coil (CC) NLR (CNL) or toll interleukin-1 (TIR) NLR (TNL), receptor-like kinase (RLK), receptor-like protein kinase (RLP), transmembrane-coiled-coil (TM-CC), WRKY, and mildew locus O (MLO) genes. These results constitute an important genomic tool for resistance breeding and gene cloning for PM in flax.

John L, Jones T, Kitchen S, Hulke BS. Assessment of biogeographic variation in traits of Lewis flax (Linum lewisii) for use in restoration and agriculture

Lewis flax (Linum lewisii) is widely distributed across western North America and is currently used in native ecosystem restoration. There is also growing interest in de novo domestication of Lewis flax as a perennial oilseed crop. To better understand this species and facilitate both restoration and domestication, we used common gardens to assess biogeographical variation in a variety of seed and growth traits from 37 flax accessions, consisting of 35 wild populations from the Intermountain West region, the pre-variety germplasm Maple Grove (L. lewisii) and the cultivar 'Appar' (L. perenne) and related this variation to collection site geography and climate. Results from linear mixed models suggest there is extensive phenotypic variation among populations of Lewis flax within the Intermountain West. Using a multivariate approach, we identify a key suite of traits that are related to latitude and climate and may facilitate adaptation, including flowering indeterminacy, seed mass and stem number. These traits should be taken into account when considering the release of new germplasm for restoration efforts. We also find that Lewis flax seed contains desirably high amounts of alpha-linolenic acid and is otherwise mostly indistinguishable in fatty acid composition from oil-type varieties of domesticated flax (L. usitatissimum), making it a strong candidate for domestication. This study provides fundamental knowledge for future research into the ecology and evolution of Lewis flax, which will inform its use in both restoration and agriculture.

Transcriptomic Analysis of Early Flowering Signals in ‘Royal’ Flax

Canada is one of the world’s leading producers and exporters of flax seed, with most production occurring in the Prairie Provinces. However, reduced season length and risk of frost restricts production in the northern grain belt of the Canadian Prairies. To expand the growing region of flax and increase production in Canada, flax breeders need to develop earlier-flowering varieties capable of avoiding the risk of abiotic stress. A thorough understanding of flowering control of flax is essential for the efficient breeding of such lines. We identified 722 putative flax flowering genes that span all major flowering-time pathways. Frequently, we found multiple flax homologues for a single Arabidopsis flowering gene. We used RNA sequencing to quantify the expression of genes in the shoot apical meristem (SAM) at 10, 15, 19, and 29 days after planting (dap) using the ‘Royal’ cultivar. We observed the expression of 80% of putative flax flowering genes and the differential expression of only 30%; these included homologues of major flowering regulators, such as SOC1, FUL, and AP1. We also found enrichment of differentially expressed genes (DEGs) in transcription factor (TF) families involved in flowering. Finally, we identified the candidates’ novel flowering genes amongst the uncharacterized flax genes. Our transcriptomic dataset provides a useful resource for investigating the regulatory control of the transition to flowering in flax and for the breeding of northern-adapted varieties.

Fine-mapping of a putative glutathione S-transferase (GST) gene responsible for yellow seed colour in flax (Linum usitatissimum)

Objective

The brown seed coat colour of flax (Linum ustiatissimum) results from proanthocyanidin synthesis and accumulation. Glutathione S-transferases (GSTs), such as the TT19 protein in Arabidopsis, have been implicated in the transport of anthocyanidins during the synthesis of the brown proanthocyanidins. This study fine mapped the g allele responsible for yellow seed colour in S95407 and identified it as a putative mutated GST.

Results

We developed a Recombinant Inbred Line population with 320 lines descended from a cross between CDC Bethune (brown seed coat) and S95407 (yellow seed) and used molecular markers to fine map the G gene on Chromosome 6 (Chr 6). We used Next Generation Sequencing (NGS) to identify a putative GST was identified in this region and Sanger sequenced the gene from CDC Bethune, S95407 and other yellow seeded genotypes. The putative GST from S95407 had 13 SNPs encoding, including four non-synonymous amino acid changes, compared to the CDC Bethune reference sequence and the other genotypes. The GST encoded by Lus10019895 is a lambda-GST in contrast to the Arabidopsis TT19 which is a phi-GST.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13104-022-05964-x.

Heat-responsive microRNAs and phased small interfering RNAs in reproductive development of flax

Plants will face increased heat stress due to rising global temperatures. Heat stress affects plant reproductive development and decreases productivity; however, the underlying molecular mechanisms of these processes are poorly characterized. Plant small RNAs (sRNAs) have important regulatory roles in plant reproductive development following abiotic stress responses. We generated sRNA transcriptomes of reproductive bud stages at three different time points to identify sRNA-mediated pathways responsive to heat stress in flax (Linum usitatissimum). With added sRNA transcriptomes of vegetative tissues, we comprehensively annotated miRNA and phasiRNA-encoding genes (PHAS) in flax. We identified 173 miRNA genes, of which 42 are newly annotated. Our analysis revealed that 141 miRNA genes were differentially accumulated between tissue types, while 18 miRNA genes were differentially accumulated in reproductive tissues following heat stress, including members of miR482/2118 and miR2275 families, known triggers of reproductive phasiRNAs. Furthermore, we identified 68 21-PHAS flax loci from protein-coding and noncoding regions, four 24-PHAS loci triggered by miR2275, and 658 24-PHAS-like loci with unknown triggers, derived mostly from noncoding regions. The reproductive phasiRNAs are mostly downregulated in response to heat stress. Overall, we found that several previously unreported miRNAs and phasiRNAs are responsive to heat stress in flax reproductive tissues.

Genes Associated with the Flax Plant Type (Oil or Fiber) Identified Based on Genome and Transcriptome Sequencing Data

As a result of the breeding process, there are two main types of flax (Linum usitatissimum L.) plants. Linseed is used for obtaining seeds, while fiber flax is used for fiber production. We aimed to identify the genes associated with the flax plant type, which could be important for the formation of agronomically valuable traits. A search for polymorphisms was performed in genes involved in the biosynthesis of cell wall components, lignans, fatty acids, and ion transport based on genome sequencing data for 191 flax varieties. For 143 of the 424 studied genes (4CL, C3'H, C4H, CAD, CCR, CCoAOMT, COMT, F5H, HCT, PAL, CTL, BGAL, ABC, HMA, DIR, PLR, UGT, TUB, CESA, RGL, FAD, SAD, and ACT families), one or more polymorphisms had a strong correlation with the flax type. Based on the transcriptome sequencing data, we evaluated the expression levels for each flax type-associated gene in a wide range of tissues and suggested genes that are important for the formation of linseed or fiber flax traits. Such genes were probably subjected to the selection press and can determine not only the traits of seeds and stems but also the characteristics of the root system or resistance to stresses at a particular stage of development, which indirectly affects the ability of flax plants to produce seeds or fiber.

Genome-wide identification of RING finger genes in flax (Linum usitatissimum) and analyses of their evolution

Background

Flax (Linum usitatissimum) is an important crop for its seed oil and stem fiber. Really Interesting New Gene (RING) finger genes play essential roles in growth, development, and biotic and abiotic stress responses in plants. However, little is known about these genes in flax.

Methods

Here, we performed a systematic genome-wide analysis to identify RING finger genes in flax.

Results

We identified 587 RING domains in 574 proteins and classified them into RING-H2 (292), RING-HCa (181), RING-HCb (23), RING-v (53), RING-C2 (31), RING-D (2), RING-S/T (3), and RING-G (2). These proteins were further divided into 45 groups according to domain organization. These genes were located in 15 chromosomes and clustered into three clades according to their phylogenetic relationships. A total of 312 segmental duplicated gene pairs were inferred from 411 RING finger genes, indicating a major contribution of segmental duplications to the RING finger gene family expansion. The non-synonymous/synonymous substitution ratio of the segmentally duplicated gene pairs was less than 1, suggesting that the gene family was under negative selection since duplication. Further, most RING genes in flax were differentially expressed during seed development or in the shoot apex. This study provides useful information for further functional analysis of RING finger genes in flax and to develop gene-derived molecular markers in flax breeding.

Chromosome identification in oil palm (Elaeis guineensis) using in situ hybridization with massive pools of single copy oligonucleotides and transferability across Arecaceae species

Chromosome identification is essential for linking sequence and chromosomal maps, verifying sequence assemblies, showing structural variations and tracking inheritance or recombination of chromosomes and chromosomal segments during evolution and breeding programs. Unfortunately, identification of individual chromosomes and chromosome arms has been a major challenge for some economically important crop species with a near-continuous chromosome size range and similar morphology. Here, we developed oligonucleotide-based chromosome-specific probes that enabled us to establish a reference chromosome identification system for oil palm (Elaeis guineensis Jacq., 2n = 32). Massive oligonucleotide sequence pools were anchored to individual chromosome arms using dual and triple fluorescent in situ hybridization (EgOligoFISH). Three fluorescently tagged probe libraries were developed to contain, in total 52,506 gene-rich single-copy 47-mer oligonucleotides spanning each 0.2-0.5 Mb across strategically placed chromosome regions. They generated 19 distinct FISH signals and together with rDNA probes enabled identification of all 32 E. guineensis chromosome arms. The probes were able to identify individual homoeologous chromosome regions in the related Arecaceae palm species: American oil palm (Elaeis oleifera), date palm (Phoenix dactylifera) and coconut (Cocos nucifera) showing the comparative organization and concerted evolution of genomes in the Arecaceae. The oligonucleotide probes developed here provide a valuable approach to chromosome arm identification and allow tracking chromosome transfer in hybridization and breeding programs in oil palm, as well as comparative studies within Arecaceae.

Efficient construction of a linkage map and haplotypes for Mentha suaveolens using sequence capture

The sustainability of many crops is hindered by the lack of genomic resources and a poor understanding of natural genetic diversity. Particularly, application of modern breeding requires high-density linkage maps that are integrated into a highly contiguous reference genome. Here, we present a rapid method for deriving haplotypes and developing linkage maps, and its application to Mentha suaveolens, one of the diploid progenitors of cultivated mints. Using sequence-capture via DNA hybridization to target single nucleotide polymorphisms (SNPs), we successfully genotyped ∼5000 SNPs within the genome of >400 individuals derived from a self cross. After stringent quality control, and identification of nonredundant SNPs, 1919 informative SNPs were retained for linkage map construction. The resulting linkage map defined a total genetic space of 942.17 cM divided among 12 linkage groups, ranging from 56.32 to 122.61 cM in length. The linkage map is in good agreement with pseudomolecules from our preliminary genome assembly, proving this resource effective for the correction and validation of the reference genome. We discuss the advantages of this method for the rapid creation of linkage maps.

Loci harboring genes with important role in drought and related abiotic stress responses in flax revealed by multiple GWAS models

QTNs associated with drought tolerance traits and indices were identified in a flax mini-core collection through multiple GWAS models and phenotyping at multiple locations under irrigated and non-irrigated field conditions. Drought is a critical phenomenon challenging today's agricultural sector. Crop varieties adapted to moisture deficit are becoming vital. Flax can be greatly affected by limiting moisture conditions, especially during the early development and reproductive stages. Here, a mini-core collection comprising genotypes from more than 20 major growing countries was evaluated for 11 drought-related traits in irrigated and non-irrigated fields for 3 years. Heritability of the traits ranged from 44.7 to 86%. Six of the 11 traits showed significant phenotypic difference between irrigated and non-irrigated conditions. A genome-wide association study (GWAS) was performed for these six traits and their corresponding stress indices based on 106 genotypes and 12,316 single nucleotide polymorphisms (SNPs) using six multi-locus and one single-locus models. The SNPs were then assigned to 8050 linkage disequilibrium (LD) blocks to which a restricted two-stage multi-locus multi-allele GWAS was applied. A total of 144 quantitative trait nucleotides (QTNs) and 13 LD blocks were associated with at least one trait or stress index. Of these, 16 explained more than 15% of the genetic variance. Most large-effect QTN loci harbored gene(s) previously predicted to play role(s) in the associated traits. Genes mediating responses to abiotic stresses resided at loci associated with stress indices. Flax genes Lus10009480 and Lus10030150 that are predicted to encode WAX INDUCER1 and STRESS-ASSOCIATED PROTEIN (SAP), respectively, are among the important candidates detected. Accessions with multiple favorable alleles outperformed others for grain yield, thousand seed weight and fiber/biomass in non-irrigated conditions, suggesting their potential usefulness in breeding and genomic selection.

The Toolbox for Fiber Flax Breeding: A Pipeline From Gene Expression to Fiber Quality

The goal of any plant breeding program is to improve quality of a target crop. Crop quality is a comprehensive feature largely determined by biological background. To improve the quality parameters of crops grown for the production of fiber, a functional approach was used to search for genes suitable for the effective manipulation of technical fiber quality. A key step was to identify genes with tissue and stage-specific pattern of expression in the developing fibers. In the current study, we investigated the relationship between gene expression evaluated in bast fibers of developing flax plants and the quality parameters of technical fibers measured after plant harvesting. Based on previously published transcriptomic data, two sets of genes that are upregulated in fibers during intrusive growth and tertiary cell wall deposition were selected. The expression level of the selected genes and fiber quality parameters were measured in fiber flax, linseed (oil flax) cultivars, and wild species that differ in type of yield and fiber quality parameters. Based on gene expression data, linear regression models for technical stem length, fiber tensile strength, and fiber flexibility were constructed, resulting in the identification of genes that have high potential for manipulating fiber quality. Chromosomal localization and single nucleotide polymorphism distribution in the selected genes were characterized for the efficacy of their use in conventional breeding and genome editing programs. Transcriptome-based selection is a highly targeted functional approach that could be used during the development of new cultivars of various crops.

Genome-wide identification of ATP binding cassette (ABC) transporter and heavy metal associated (HMA) gene families in flax (Linum usitatissimum L.)

Background

The recent release of the reference genome sequence assembly of flax, a self-pollinated crop with 15 chromosome pairs, into chromosome-scale pseudomolecules enables the characterization of gene families. The ABC transporter and HMA gene families are important in the control of cadmium (Cd) accumulation in crops. To date, the genome-wide analysis of these two gene families has been successfully conducted in some plant species, but no systematic evolutionary analysis is available for the flax genome.

Results

Here we describe the ABC transporter and HMA gene families in flax to provide a comprehensive overview of its evolution and some support towards the functional annotation of its members. The 198 ABC transporter and 12 HMA genes identified in the flax genome were classified into eight ABC transporter and four HMA subfamilies based on their phylogenetic analysis and domains’ composition. Nine of these genes, i.e., LuABCC9, LuABCC10, LuABCG58, LuABCG59, LuABCG71, LuABCG72, LuABCG73, LuHMA3, and LuHMA4, were orthologous with the Cd associated genes in Arabidopsis, rice and maize. Ten motifs were identified from all ABC transporter and HMA genes. Also, several motifs were conserved among genes of similar length, but each subfamily each had their own motif structures. Both the ABC transporter and HMA gene families were highly conserved among subfamilies of flax and with those of Arabidopsis. While four types of gene duplication were observed at different frequencies, whole-genome or segmental duplications were the most frequent with 162 genes, followed by 29 dispersed, 14 tandem and 4 proximal duplications, suggesting that segmental duplications contributed the most to the expansion of both gene families in flax. The rates of non-synonymous to synonymous (Ka/Ks) mutations of paired duplicated genes were for the most part lower than one, indicative of a predominant purifying selection. Only five pairs of genes clearly exhibited positive selection with a Ka/Ks ratio greater than one. Gene ontology analyses suggested that most flax ABC transporter and HMA genes had a role in ATP binding, transport, catalytic activity, ATPase activity, and metal ion binding. The RNA-Seq analysis of eight different organs demonstrated diversified expression profiling patterns of the genes and revealed their functional or sub-functional conservation and neo-functionalization.

Conclusion

Characterization of the ABC transporter and HMA gene families will help in the functional analysis of candidate genes in flax and other crop species.

Genetic diversity of SAD and FAD genes responsible for the fatty acid composition in flax cultivars and lines

BACKGROUND

Flax (Linum usitatissimum L.) is grown for fiber and seed in many countries. Flax cultivars differ in the oil composition and, depending on the ratio of fatty acids, are used in pharmaceutical, food, or paint industries. It is known that genes of SAD (stearoyl-ACP desaturase) and FAD (fatty acid desaturase) families play a key role in the synthesis of fatty acids, and some alleles of these genes are associated with a certain composition of flax oil. However, data on genetic polymorphism of these genes are still insufficient.

RESULTS

On the basis of the collection of the Institute for Flax (Torzhok, Russia), we formed a representative set of 84 cultivars and lines reflecting the diversity of fatty acid composition of flax oil. An approach for the determination of full-length sequences of SAD1, SAD2, FAD2A, FAD2B, FAD3A, and FAD3B genes using the Illumina platform was developed and deep sequencing of the 6 genes in 84 flax samples was performed on MiSeq. The obtained high coverage (about 400x on average) enabled accurate assessment of polymorphisms in SAD1, SAD2, FAD2A, FAD2B, FAD3A, and FAD3B genes and evaluation of cultivar/line heterogeneity. The highest level of genetic diversity was observed for FAD3A and FAD3B genes - 91 and 62 polymorphisms respectively. Correlation analysis revealed associations between particular variants in SAD and FAD genes and predominantly those fatty acids whose conversion they catalyze: SAD - stearic and oleic acids, FAD2 - oleic and linoleic acids, FAD3 - linoleic and linolenic acids. All except one low-linolenic flax cultivars/lines contained both the substitution of tryptophan to stop codon in the FAD3A gene and histidine to tyrosine substitution in the FAD3B gene, while samples with only one of these polymorphisms had medium content of linolenic acid and cultivars/lines without them were high-linolenic.

CONCLUSIONS

Genetic polymorphism of SAD and FAD genes was evaluated in the collection of flax cultivars and lines with diverse oil composition, and associations between particular polymorphisms and the ratio of fatty acids were revealed. The achieved results are the basis for the development of marker-assisted selection and DNA-based certification of flax cultivars.

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.

High-throughput sequencing techniques to flax genetics and breeding

Flax (Linum usitatissimum L.) is an important oil and fiber crop. Using modern methods for flax breeding allows accelerating the introduction of some desired genes into the genotypes of future varieties. Today, an important condition for their creation is the development of research, that is based on next-generation sequencing (NGS). This review summarizes the results obtained using NGS in flax research. To date, a linkage map with a high marker density has been obtained for L. usitatissimum, which is already being used for a more efficient search for quantitative traits loci. Comparative studies of transcriptomes and miRNomes of flax under stress and in control conditions elucidated molecular-genetic mechanisms of abiotic and biotic stress responses. The very accurate model for genomic selection of flax resistant to pasmo was constructed. Based on NGS-sequencing also some details of the genus Linum evolution were clarified. The knowledge systematized in the review can be useful for researchers working in flax breeding and whereas fundamental interest for understanding the phylogenetic relationships within the genus Linum, the ontogenesis, and the mechanisms of the response of flax plants to various stress factors.

Orbitides isolated from flaxseed induce apoptosis against SGC-7901 adenocarcinoma cells

Unique plant-derived cyclic peptides family exhibiting various key biological activities has great possibility for anticancer therapy. In this study, we investigated the effects of orbitides isolated from flax (Linum usitatissimum L.) on the growth of SGC-7901 cancer cells and the potential mechanism. Results showed that flaxseed orbitides killed off cancer cells by inducing apoptosis in a dose-dependent manner, which was confirmed by the appearance of nuclear shrinkage and DNA fragmentation, and the inhibitory effect was stronger than that of pure orbitide [1-9-NαC]-linusorb B2 or [1-9-NαC]-linusorb B3. Besides, the mitochondrial apoptosis pathway-related protein cytochrome C (Cyt C) was released from mitochondria to cytosol, associated with the activation of caspases 9 and 3, and the cleavage of PARP. Taken together, these results indicated that flaxseed orbitides induced apoptosis via the mitochondrial pathway, releasing Cyt C, increasing Bax/Bcl-2 ratio and elevating the expression of cleaved caspase 9 and 3 in SGC-7901 cells.

Mapping Quantitative Trait Loci onto Chromosome-Scale Pseudomolecules in Flax

Quantitative trait loci (QTL) are genomic regions associated with phenotype variation of quantitative traits. To date, a total of 313 QTL for 31 quantitative traits have been reported in 14 studies on flax. Of these, 200 QTL from 12 studies were identified based on genetic maps, the scaffold sequences, or the pre-released chromosome-scale pseudomolecules. Molecular markers for QTL identification differed across studies but the most used ones were simple sequence repeats (SSRs) or single nucleotide polymorphisms (SNPs). To uniquely map the SSR and SNP markers from different references onto the recently released chromosome-scale pseudomolecules, methods with several scripts and database files were developed to locate PCR- and SNP-based markers onto the same reference, co-locate QTL, and scan genome-wide candidate genes. Using these methods, 195 out of 200 QTL were successfully sorted onto the 15 flax chromosomes and grouped into 133 co-located QTL clusters; the candidate genes that co-located with these QTL clusters were also predicted. The methods and tools presented in this article facilitate marker re-mapping to a new reference, genome-wide QTL analysis, candidate gene scanning, and breeding applications in flax and other crops.

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

Genomic Prediction Accuracy of Seven Breeding Selection Traits Improved by QTL Identification in Flax

Molecular markers are one of the major factors affecting genomic prediction accuracy and the cost of genomic selection (GS). Previous studies have indicated that the use of quantitative trait loci (QTL) as markers in GS significantly increases prediction accuracy compared with genome-wide random single nucleotide polymorphism (SNP) markers. To optimize the selection of QTL markers in GS, a set of 260 lines from bi-parental populations with 17,277 genome-wide SNPs were used to evaluate the prediction accuracy for seed yield (YLD), days to maturity (DTM), iodine value (IOD), protein (PRO), oil (OIL), linoleic acid (LIO), and linolenic acid (LIN) contents. These seven traits were phenotyped over four years at two locations. Identification of quantitative trait nucleotides (QTNs) for the seven traits was performed using three types of statistical models for genome-wide association study: two SNP-based single-locus (SS), seven SNP-based multi-locus (SM), and one haplotype-block-based multi-locus (BM) models. The identified QTNs were then grouped into QTL based on haplotype blocks. For all seven traits, 133, 355, and 1208 unique QTL were identified by SS, SM, and BM, respectively. A total of 1420 unique QTL were obtained by SS+SM+BM, ranging from 254 (OIL, LIO) to 361 (YLD) for individual traits, whereas a total of 427 unique QTL were achieved by SS+SM, ranging from 56 (YLD) to 128 (LIO). SS models alone did not identify sufficient QTL for GS. The highest prediction accuracies were obtained using single-trait QTL identified by SS+SM+BM for OIL (0.929 ± 0.016), PRO (0.893 ± 0.023), YLD (0.892 ± 0.030), and DTM (0.730 ± 0.062), and by SS+SM for LIN (0.837 ± 0.053), LIO (0.835 ± 0.049), and IOD (0.835 ± 0.041). In terms of the number of QTL markers and prediction accuracy, SS+SM outperformed other models or combinations thereof. The use of all SNPs or QTL of all seven traits significantly reduced the prediction accuracy of traits. The results further validated that QTL outperformed high-density genome-wide random markers, and demonstrated that the combined use of single and multi-locus models can effectively identify a comprehensive set of QTL that improve prediction accuracy, but further studies on detection and removal of redundant or false-positive QTL to maximize prediction accuracy and minimize the number of QTL markers in GS are warranted.

The cinnamyl alcohol dehydrogenase gene family is involved in the response to Fusarium oxysporum in resistant and susceptible flax genotypes

Flax (Linum usitatissimum L.) is used for the production of textile, oils, pharmaceuticals, and composite materials. Fusarium wilt, caused by the fungus Fusarium oxysporum f. sp. lini, is a very harmful disease that reduces flax production. Flax cultivars that are resistant to Fusarium wilt have been developed, and the genes that are involved in the host response to F. oxysporum have been identified. However, the mechanisms underlying resistance to this pathogen remain unclear. In the present study, we used transcriptome sequencing data obtained from susceptible and resistant flax genotypes grown under control conditions or F. oxysporum infection. Approximately 250 million reads, generated with an Illumina NextSeq instrument, were analyzed. After filtering to exclude the F. oxysporum transcriptome, the remaining reads were mapped to the L. usitatissimum genome and quantified. Then, the expression levels of cinnamyl alcohol dehydrogenase (CAD) family genes, which are known to be involved in the response to F. oxysporum, were evaluated in resistant and susceptible flax genotypes. Expression alterations in response to the pathogen were detected for all 13 examined CAD genes. The most significant differences in expression between control and infected plants were observed for CAD1B, CAD4A, CAD5A, and CAD5B, with strong upregulation of CAD1B, CAD5A, and CAD5B and strong downregulation of CAD4A. When plants were grown under the same conditions, the expression levels were similar in all studied flax genotypes for most CAD genes, and statistically significant differences in expression between resistant and susceptible genotypes were only observed for CAD1A. Our study indicates the strong involvement of CAD genes in flax response to F. oxysporum but brings no evidence of their role as resistance gene candidates. These findings contribute to the understanding of the mechanisms underlying the response of flax to F. oxysporum infection and the role of CAD genes in stress resistance.

The Complex Genetic Architecture of Early Root and Shoot Traits in Flax Revealed by Genome-Wide Association Analyses

Roots are fundamental organs for water and nutrient uptake as well as for signal transduction in response to biotic and abiotic stresses. Flax has a shallow tap root system that relies mostly on top soil nutrient and moisture resources. The crop can easily be outcompeted by weeds or other crops in intercropping systems, especially in moisture deficit conditions. However, there is a wide range of variation among genotypes in terms of performance under scarce resources such as moisture limitation. Here we phenotyped 15 root, two shoot traits and shoot to root dry weight ratio on 115 flax accessions grown in a hydroponic pouch system and performed a genome-wide association study (GWAS) based on seven different models to identify quantitative trait loci underlying these traits. Significant variation among genotypes was observed for the two shoot and 12 of the 14 root traits. Shoot dry weight was correlated with root network volume, length, surface area, and root dry weight (r > 0.5, P < 0.001) but not significantly correlated with root depth (r = 0.033, P > 0.05). The seven GWAS models detected a total of 228 quantitative trait nucleotides (QTNs) for 16 traits. Most loci, defined by an interval of 100 kb up and downstream of the QTNs, harbored genes known to play role(s) in root and shoot development, suggesting them as candidates. Examples of candidate genes linked to root network QTNs included genes encoding GRAS transcription factors, mitogen-activated protein kinases, and auxin related lateral organ boundary proteins while QTN loci for shoot dry weight harbored genes involved in photomorphogenesis and plant immunity. These results provide insights into the genetic bases of early shoot and root development traits in flax that could be capitalized upon to improve its root architecture, particularly in view of better withstanding water limiting conditions during the cropping season.

Genome-wide identification and expression analysis of R2R3, 3R- and 4R-MYB transcription factors during lignin biosynthesis in flax (Linum usitatissimum)

MYB transcription factors (TFs) have vital roles in regulating lignin or fiber development. Flax (Linum usitatissimum) is known as one of the plants with high fiber production capacity. However, no studies have been conducted to identify and characterize MYB TFs in the flax genome. Results showed that flax genome harbours 167 R2R3, seven 3R, and one 4R-type MYB TFs. 22 MYB genes (%13) were estimated to be tandem duplicated dated around 13.3-86.98 Mya. 130 flax MYB members have apparent orthologous with Arabidopsis, in which 17 R2R3 MYBs are associated with lignin biosynthesis. MYB062, MYB072, MYB096, MYB141, and MYB146 genes were up-regulated in tissues having higher lignin production capacity. In opposite, MYB012 and MYB113 genes were down-regulated which points out the involvement of those genes in the lignin biosynthesis mechanism. This comprehensive study can provide a basis for understanding the role of MYBs in fiber or lignin production in flax.

The genetic structure of flax illustrates environmental and anthropogenic selections that gave rise to its eco-geographical adaptation

Flax, one of the eight founder crops of agriculture, has been cultivated for its oil and/or fiber for millennia. Understanding genetic divergence and geographic origins of germplasm in line with their cultivation history and ecological adaptation are essential for conservation and breeding. Here we performed a genome-wide assessment based on more than 51,000 single nucleotide polymorphic sites defining 383 flax accessions from a core collection representing 37 flax growing countries. Population structure analysis resulted in a total of 12 populations that were pooled into four major groups: Temperate, South Asian, Abyssinian and Mediterranean. The vast majority (n = 335) belonged to the Temperate group that comprised eight populations including one dominated by fiber flax. Genetic variation between fiber and oil morphotypes was less pronounced than variation within morphotypes. The genetic variation among groups and populations was attributed in part to eco-geographic and anthropogenic factors. Genetic signatures indicated loci under strong selection by environmental factors such as day length. A high concentration of private haplotypes were observed in the South Asian, Mediterranean and Abyssinian populations despite their low genotype representation, hinting at the long history of the crop in these regions. The addition of genotypes from these three regions would enrich the core collection by capturing a wider genetic breadth for breeding and conservation.

Genomic insights into HSFs as candidate genes for high-temperature stress adaptation and gene editing with minimal off-target effects in flax

Flax (Linum usitatissimum) is a cool season crop commercially cultivated for seed oil and stem fibre production. A comprehensive characterization of the heat shock factor (HSF) candidate genes in flax can accelerate genetic improvement and adaptive breeding for high temperature stress tolerance. We report the genome-wide identification of 34 putative HSF genes from the flax genome, which we mapped on 14 of the 15 chromosomes. Through comparative homology analysis, we classified these genes into three broad groups, and sub-groups. The arrangement of HSF-specific protein motifs, DNA-binding domain (DBD) and hydrophobic heptad repeat (HR-A/B), and exon-intron boundaries substantiated the phylogenetic separation of these genes. Orthologous relationships and evolutionary analysis revealed that the co-evolution of the LusHSF genes was due to recent genome duplication events. Digital and RT-qPCR analyses provided significant evidence of the differential expression of the LusHSF genes in various tissues, at various developmental stages, and in response to high-temperature stress. The co-localization of diverse cis-acting elements in the promoters of the LusHSF genes further emphasized their regulatory roles in the abiotic stress response. We further confirmed DNA-binding sites on the LusHSF proteins and designed guide RNA sequences for gene editing with minimal off-target effects. These results will hasten functional investigations of LusHSFs or assist in devising genome engineering strategies to develop high-temperature stress tolerant flax cultivars.

Characterization of repeated DNA sequences in genomes of blue-flowered flax

BACKGROUND

Members of different sections of the genus Linum are characterized by wide variability in size, morphology and number of chromosomes in karyotypes. Since such variability is determined mainly by the amount and composition of repeated sequences, we conducted a comparative study of the repeatomes of species from four sections forming a clade of blue-flowered flax. Based on the results of high-throughput genome sequencing performed in this study as well as available WGS data, bioinformatic analyses of repeated sequences from 12 flax samples were carried out using a graph-based clustering method.

RESULTS

It was found that the genomes of closely related species, which have a similar karyotype structure, are also similar in the repeatome composition. In contrast, the repeatomes of karyologically distinct species differed significantly, and no similar tandem-organized repeats have been identified in their genomes. At the same time, many common mobile element families have been identified in genomes of all species, among them, Athila Ty3/gypsy LTR retrotransposon was the most abundant. The 30-chromosome members of the sect. Linum (including the cultivated species L. usitatissimum) differed significantly from other studied species by a great number of satellite DNA families as well as their relative content in genomes.

CONCLUSIONS

The evolution of studied flax species was accompanied by waves of amplification of satellite DNAs and LTR retrotransposons. The observed inverse correlation between the total contents of dispersed repeats and satellite DNAs allowed to suggest a relationship between both classes of repeating sequences. Significant interspecific differences in satellite DNA sets indicated a high rate of evolution of this genomic fraction. The phylogenetic relationships between the investigated flax species, obtained by comparison of the repeatomes, agreed with the results of previous molecular phylogenetic studies.