Quality assessment of flax advanced breeding lines varying in seed coat color and their potential use in the food and industrial applications.

BACKGROUND: With the increasing consumer awareness of the strong relationship between food and health, flax became a promising functional food due to its bioactive nutraceutical composition. Intra-specific crosses of eight contrasting flax genotypes were performed previously, and within segregating F6 progeny families, we investigated a close-up composition of phytochemicals derived from whole seeds. RESULTS: The considerable genetic variation among the flax F6 families suggested that intra-specific hybridization is essential in flax breeding to obtain and broaden genetic variability and largely affirmed the opportunity for selecting promising lines. Also, significant variations in the targeted metabolite contents and antioxidant properties were observed among brown and yellow-seeded families. Notably, brown-seeded families expressed the highest average values of saturated fatty acids, protein, fiber, tocopherol, phenolics, SDG, and SECO lignans. Yellow-seeded families represented the highest average content of unsaturated fatty acids and mucilage. The cultivation year significantly affects flaxseed's composition and functional properties, presumably due to temperature, humidity, and sunshine time differences. Interestingly, the seeds obtained in warmer conditions were more potent and had more chemical constituents. The favorable genetic correlations among all evaluated traits suggest the possibility of joint genetic selection for several nutritional and phytochemical characteristics in flax. The current study highlights the importance and utilization of 19 top families as their seeds and oil play imperative roles in the pharmaceuticals and food industries. The antioxidant capacity of the seeds showed that families 84B, 23B, 35Y, 95Y, 30B, 88B, and 78B serve as a natural source of dietary antioxidants beneficial to human health. To increase the oxidative stability of the flaxseed oil, the quality evaluation identified some families with low levels of linolenic acid. CONCLUSIONS: These findings are essential to improving flaxseed's nutritional quality and therapeutic properties through a bulk breeding program.

Investigation of genetic diversity using molecular and biochemical markers associated with powdery mildew resistance in different flax (Linum usitatissimum L.) genotypes.

Under greenhouse conditions, the resistance of 18 different genotypes of flax to powdery mildew was evaluated. To investigate genetic diversity and identify the molecular and biochemical markers linked to powdery mildew resistance in the tested genotypes, two molecular marker systems-start codon targeted (SCoT) and inter-simple sequence repeat (ISSR)-as well as a biochemical marker (protein profiles, antioxidant enzyme activity, and secondary metabolites) were used. Based on the results, the genotypes were classified into four categories: highly susceptible, susceptible, moderately susceptible, and moderately resistant. The genotypes differed significantly in powdery mildew severity: Polk had a severity of 92.03% and Leona had a severity of 18.10%. Compared to the other genotypes, the moderately resistant genotypes had higher levels of flavonoids, antioxidant enzymes, phenolics, and straw yield; nevertheless, their hydrogen peroxide and malondialdehyde levels were lower. Protein profiles revealed 93.75% polymorphism, although the ISSR marker displayed more polymorphism (78.4%) than the SCoT marker (59.7%). Specific molecular and biochemical markers associated with powdery mildew resistance were identified. The 18 genotypes of flax were divided into two major clusters by the dendrogram based on the combined data of molecular markers. The first main cluster included Leona (genotype number 7), considered moderate resistance to powdery mildew and a separate phenetic line. The second main cluster included the other 17 genotypes, which are grouped together in a sub-cluster. This means that, besides SCoT, ISSR markers can be a useful supplementary technique for molecular flax characterization and for identifying genetic associations between flax genotypes under powdery mildew infection.

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.

Changes in microRNAs expression of flax (Linum usitatissimum L.) planted in a cadmium-contaminated soil following the inoculation with root symbiotic fungi.

Cadmium is one of the most harmful heavy metals that harm agricultural products. Evaluating microRNAs expression is a new and accurate method to study plant response in various environmental conditions. So this study aimed to evaluate the contribution of two symbiotic fungi in improving flax tolerance in a Cd-polluted soil using microRNAs and their target gene expression. A factorial pot experiment in a completely randomized design was conducted with different levels of Cd (0, 20, and 40 mg kg-1) on non-inoculated and inoculated flax with Claroideoglomus etunicatum and Serendipita indica. The results presented that increasing Cd levels caused a constant decline of alkaline phosphatase of soil (from 243 to 210 and 153 µg PNP g-1 h-1), respectively, from control (Cd0) to 20 and 40 mg Cd kg-1. However, the inoculation of flax with fungi significantly enhanced these properties. A negative correlation was observed between the expression level of microRNA 167 and microRNA 398 with their corresponding target genes, auxin response factor 8 and superoxide dismutase zinc/copper 1, respectively. The expression level of both microRNAs and their targets indicated that the inoculation with symbiont fungi could diminish Cd stress and enhance the growth of flax.

Soil contamination with Cd affects plant growth.Root symbiotic fungi can improve plant growth in Cd-polluted soils.Examining microRNA expression is a new and accurate method to evaluate plant response to Cd pollution and symbiotic fungi.

From Gene-for-Gene to Resistosomes: Flor's Enduring Legacy.

The gene-for-gene model proposed by H. H. Flor has been one of the fundamental precepts of plant-pathogen interactions that has underpinned decades of research towards our current concepts of plant immunity. The broad validity of this model as an elegant and accurate genetic description of specific recognition events between the products of plant resistance (R) and pathogen avirulence (Avr) genes has been demonstrated many times over in a wide variety of plant disease systems. In recent years detailed molecular and structural analyses have provided a deep understanding of the principles by which plant immune receptors recognize pathogen effectors, including providing molecular descriptions of many of the genetic loci in flax and flax rust characterized by Flor. Recent advances in molecular and structural understanding of immune receptor recognition and activation mechanisms have brought the field to a new level, where rational design of novel receptors through engineering approaches is becoming a realizable goal. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Drivers of diversification in Linum (Linaceae) by means of chromosome evolution: correlations with biogeography, breeding system and habit.

BACKGROUND AND AIMS: Chromosome evolution leads to hybrid dysfunction and recombination patterns and has thus been proposed as a major driver of diversification in all branches of the tree of life, including flowering plants. In this study we used the genus Linum (flax species) to evaluate the effects of chromosomal evolution on diversification rates and on traits that are important for sexual reproduction. Linum is a useful study group because it has considerable reproductive polymorphism (heterostyly) and chromosomal variation (n = 6-36) and a complex pattern of biogeographical distribution. METHODS: We tested several traditional hypotheses of chromosomal evolution. We analysed changes in chromosome number across the phylogenetic tree (ChromEvol model) in combination with diversification rates (ChromoSSE model), biogeographical distribution, heterostyly and habit (ChromePlus model). KEY RESULTS: Chromosome number evolved across the Linum phylogeny from an estimated ancestral chromosome number of n = 9. While there were few apparent incidences of cladogenesis through chromosome evolution, we inferred up to five chromosomal speciation events. Chromosome evolution was not related to heterostyly but did show significant relationships with habit and geographical range. Polyploidy was negatively correlated with perennial habit, as expected from the relative commonness of perennial woodiness and absence of perennial clonality in the genus. The colonization of new areas was linked to genome rearrangements (polyploidy and dysploidy), which could be associated with speciation events during the colonization process. CONCLUSIONS: Chromosome evolution is a key trait in some clades of the Linum phylogeny. Chromosome evolution directly impacts speciation and indirectly influences biogeographical processes and important plant traits.

The pattern of alternative splicing and DNA methylation alteration and their interaction in linseed (Linum usitatissimum L.) response to repeated drought stresses.

BACKGROUND: Drought stress has significantly hampered agricultural productivity worldwide and can also result in modifications to DNA methylation levels. However, the dynamics of DNA methylation and its association with the changes in gene transcription and alternative splicing (AS) under drought stress are unknown in linseed, which is frequently cultivated in arid and semiarid regions. RESULTS: We analysed AS events and DNA methylation patterns in drought-tolerant (Z141) and drought-sensitive (NY-17) linseed under drought stress (DS) and repeated drought stress (RD) treatments. We found that the number of intron-retention (IR) and alternative 3' splice site (Alt3'SS) events were significantly higher in Z141 and NY-17 under drought stress. We found that the linseed response to the DS treatment was mainly regulated by transcription, while the response to the RD treatment was coregulated by transcription and AS. Whole genome-wide DNA methylation analysis revealed that drought stress caused an increase in the overall methylation level of linseed. Although we did not observe any correlation between differentially methylated genes (DMGs) and differentially spliced genes (DSGs) in this study, we found that the DSGs whose gene body region was hypermethylated in Z141 and hypomethylated in NY-17 were enriched in abiotic stress response Gene Ontology (GO) terms. This finding implies that gene body methylation plays an important role in AS regulation in some specific genes. CONCLUSION: Our study is the first comprehensive genome-wide analysis of the relationship between linseed methylation changes and AS under drought and repeated drought stress. Our study revealed different interaction patterns between differentially expressed genes (DEGs) and DSGs under DS and RD treatments and differences between methylation and AS regulation in drought-tolerant and drought-sensitive linseed varieties. The findings will probably be of interest in the future. Our results provide interesting insights into the association between gene expression, AS, and DNA methylation in linseed under drought stress. Differences in these associations may account for the differences in linseed drought tolerance.

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.

A LuALS Mutation with High Sulfonylurea Herbicide Resistance in Linum usitatissimum L.

The cultivation of herbicide-resistant crops is an effective tool for weed management in agriculture. Weed control in flax (Linum usitatissimum L.) remains challenging due to the lack of available herbicide-resistant cultivars. In this study, a mutant resistant to acetolactate synthase (ALS)-inhibiting herbicides was obtained by ethyl methanesulphonate (EMS) mutagenesis using an elite cultivar, Longya10. Whole-plant dose-response assays revealed that, compared to Longya10, the mutant was 11.57-fold more resistant to tribenuron-methyl (TBM) and slightly resistant to imazethapyr (resistance index (mutant/Longya10) < 3). In vitro acetolactate synthase assays showed that the relative resistance of the mutant was 12.63 times more than that of Longya10. A biochemical analysis indicated that there was a Pro197Ser (relative to the Arabidopsis thaliana ALS sequence) substitution within the LuALS1, conferring high resistance to sulfonylurea herbicides in the mutant. Additionally, two cleaved amplified polymorphic sequence (CAPS) markers, BsaI-LuALS1 and EcoO109I-LuALS1, were developed based on the mutation site for marker assistant selection in breeding. Moreover, the mutant did not cause losses in natural field conditions. We find a mutant with ALS-inhibiting herbicide resistance chemically induced by EMS mutagenesis, providing a valuable germplasm for breeding herbicide-resistant flax varieties.

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 Identification and Expression Analysis of Auxin Response Factor Gene Family in Linum usitatissimum.

Auxin response factors (ARFs) are critical components of the auxin signaling pathway, and are involved in diverse plant biological processes. However, ARF genes have not been investigated in flax (Linum usitatissimum L.), an important oilseed and fiber crop. In this study, we comprehensively analyzed the ARF gene family and identified 33 LuARF genes unevenly distributed on the 13 chromosomes of Longya-10, an oil-use flax variety. Detailed analysis revealed wide variation among the ARF family members and predicted nuclear localization for all proteins. Nineteen LuARFs contained a complete ARF structure, including DBD, MR, and CTD, whereas the other fourteen lacked the CTD. Phylogenetic analysis grouped the LuARFs into four (I-V) clades. Combined with sequence analysis, the LuARFs from the same clade showed structural conservation, implying functional redundancy. Duplication analysis identified twenty-seven whole-genome-duplicated LuARF genes and four tandem-duplicated LuARF genes. These duplicated gene pairs' Ka/Ks ratios suggested a strong purifying selection pressure on the LuARF genes. Collinearity analysis revealed that about half of the LuARF genes had homologs in other species, indicating a relatively conserved nature of the ARFs. The promoter analysis identified numerous hormone- and stress-related elements, and the qRT-PCR experiment revealed that all LuARF genes were responsive to phytohormone (IAA, GA3, and NAA) and stress (PEG, NaCl, cold, and heat) treatments. Finally, expression profiling of LuARF genes in different tissues by qRT-PCR indicated their specific functions in stem or capsule growth. Thus, our findings suggest the potential functions of LuARFs in flax growth and response to an exogenous stimulus, providing a basis for further functional studies on these genes.

Genome-Wide Analysis of Flax (Linum usitatissimum L.) Growth-Regulating Factor (GRF) Transcription Factors.

Flax is an important cash crop globally with a variety of commercial uses. It has been widely used for fiber, oil, nutrition, feed and in composite materials. Growth regulatory factor (GRF) is a transcription factor family unique to plants, and is involved in regulating many processes of growth and development. Bioinformatics analysis of the GRF family in flax predicted 17 LuGRF genes, which all contained the characteristic QLQ and WRC domains. Equally, 15 of 17 LuGRFs (88%) are predicted to be regulated by lus-miR396 miRNA. Phylogenetic analysis of GRFs from flax and several other well-characterized species defined five clades; LuGRF genes were found in four clades. Most LuGRF gene promoters contained cis-regulatory elements known to be responsive to hormones and stress. The chromosomal locations and collinearity of LuGRF genes were also analyzed. The three-dimensional structure of LuGRF proteins was predicted using homology modeling. The transcript expression data indicated that most LuGRF family members were highly expressed in flax fruit and embryos, whereas LuGRF3, LuGRF12 and LuGRF16 were enriched in response to salt stress. Real-time quantitative fluorescent PCR (qRT-PCR) showed that both LuGRF1 and LuGRF11 were up-regulated under ABA and MeJA stimuli, indicating that these genes were involved in defense. LuGRF1 was demonstrated to be localized to the nucleus as expected for a transcription factor. These results provide a basis for further exploration of the molecular mechanism of LuGRF gene function and obtaining improved flax breeding lines.

Evaluation of Differentially Expressed Genes in Leaves vs. Roots Subjected to Drought Stress in Flax (Linum usitatissimum L.).

Drought stress is a common environmental challenge that plants face, severely constraining plant growth and reducing crop yield and quality. Several studies have highlighted distinct responses between monocotyledonous and dicotyledonous plants. However, the mechanisms underlying flax tolerance to abiotic stress, such as drought, remain unclear. In this study, we investigated the morphological, physiological, and biochemical characteristics and the genome-wide gene expression of oil flax and fiber flax in response to drought stress. The results revealed that drought stress caused significant wilting of flax leaves. Within the first 24 h of stress, various physiological and biochemical characteristics exhibited rapid responses. These included fresh weight, relative water content (RWC), proline, soluble protein, soluble sugar, superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT) in the leaves or roots of flax. Additionally, drought stress led to a significant rise in lignin content in fiber flax. In addition, the transcriptome analysis demonstrated genome-wide variations in gene expression induced by drought stress. Specifically, genes associated with photosynthesis, proline biosynthesis, and phytohormone metabolism exhibited significant differences in expression levels under stress conditions in flax. These findings highlight the rapid response of flax to drought stress within a short-term period. Our experiment also revealed that, although there were variations in the levels of small compound content or gene expression between Longya10 and Fany under drought stress, most stress-resistance responses were similar. Furthermore, the results provide additional evidence supporting the existence of mechanisms underlying the response to drought stress in plants.

Genome-wide analysis of genomic imprinting in the endosperm and allelic variation in flax.

Genomic imprinting is an epigenetic phenomenon that causes biased expression of maternally and paternally inherited alleles. In flowering plants, genomic imprinting predominantly occurs in the triploid endosperm and plays a vital role in seed development. In this study, we identified 248 candidate imprinted genes including 114 maternally expressed imprinted genes (MEGs) and 134 paternally expressed imprinted genes (PEGs) in flax (Linum usitatissimum L.) endosperm using deep RNA sequencing. These imprinted genes were neither clustered in specific chromosomal regions nor well conserved among flax and other plant species. MEGs tended to be expressed specifically in the endosperm, whereas the expression of PEGs was not tissue-specific. Imprinted single nucleotide polymorphisms differentiated 200 flax cultivars into the oil flax, oil-fiber dual purpose flax and fiber flax subgroups, suggesting that genomic imprinting contributed to intraspecific variation in flax. The nucleotide diversity of imprinted genes in the oil flax subgroup was significantly higher than that in the fiber flax subgroup, indicating that some imprinted genes underwent positive selection during flax domestication from oil flax to fiber flax. Moreover, imprinted genes that underwent positive selection were related to flax functions. Thirteen imprinted genes related to flax seed size and weight were identified using a candidate gene-based association study. Therefore, our study provides information for further exploration of the function and genomic variation of imprinted genes in the flax population.

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.

Effect of chronic radiation on the flax (Linum usitatissimum L.) genome grown for six consecutive generations in the radioactive Chernobyl area.

The growth of plants under chronic radiation stress in the Chernobyl area may cause changes in the genome of plants. To assess the extent of genetic and epigenetic changes in nuclear DNA, seeds of the annual crop flax (Linum usitatissimum L.) of the Kyivskyi variety, sown 21 years after the accident and grown for six generations in radioactive (RAD) and remediated (REM) fields were analysed. Flaxseed used for sowing first generation, which served as a reference (REF), was also analysed. The AFLP (Amplified Fragment Length Polymorphism) revealed a higher number of specific EcoRI-MseI loci (3.4-fold) in pooled flaxseed samples harvested from the RAD field compared with the REM field, indicating a link between the mutation process in the flax genome and the ongoing adaptation process. MSAP (Methylation-Sensitive Amplified Polymorphism) detecting EcoRI-MspI and EcoRI-HpaII loci in flax nuclear DNA genome showed no significant differences in methylation level, reaching about 33% in each of the groups studied. On the other hand, significant changes in the DNA methylation pattern of flaxseed samples harvested from the RAD field compared with controls were detected. Pairwise FST comparison revealed within both, EcoRI-MspI and transformed methylation-Sensitive data sets more than a 3-fold increase of genetic divergence in the RAD field compared with both controls. These results indicate that the nuclear genome of flax exposed to chronic radiation for six generations has more mutations and uses DNA methylation as one of the adaptation mechanisms for sustainability under adverse conditions.

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.

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.

Genetic Determinants of Fiber-Associated Traits in Flax Identified by Omics Data Integration.

In this paper, we explore potential genetic factors in control of flax phenotypes associated with fiber by mining a collection of 306 flax accessions from the Federal Research Centre of the Bast Fiber Crops, Torzhok, Russia. In total, 11 traits were assessed in the course of 3 successive years. A genome-wide association study was performed for each phenotype independently using six different single-locus models implemented in the GAPIT3 R package. Moreover, we applied a multivariate linear mixed model implemented in the GEMMA package to account for trait correlations and potential pleiotropic effects of polymorphisms. The analyses revealed a number of genomic variants associated with different fiber traits, implying the complex and polygenic control. All stable variants demonstrate a statistically significant allelic effect across all 3 years of the experiment. We tested the validity of the predicted variants using gene expression data available for the flax fiber studies. The results shed new light on the processes and pathways associated with the complex fiber traits, while the pinpointed candidate genes may be further used for marker-assisted selection.

Investigation of the Properties of Linen Fibers and Dressings.

In antiquity, flax was used as a dressing for healing wounds. Currently, work is underway on the genetic modification of flax fibers to improve their properties. Genetic modifications have resulted in an increased content of antioxidants and more favorable mechanical properties. The works published so far have presented independent tests of fibers and dressings after appropriate technological treatments in cell cultures. This study aimed to compare the properties of the fibers and the dressing produced in cell cultures-hamster fibroblasts-V79. The research material was traditional NIKE fibers; genetically modified M, B, and MB fibers; and linen dressings obtained from these fibers. The extract from 48-h incubation of 40 mg of fiber in the culture medium, which was desolved into 10, 20, and 30 mg, was administered to the cell culture. On the other hand, a linen dressing was placed on cells with an area of 0.5 cm2, 1 cm2, 1.5 cm2, and 2 cm2. Cells with fiber or dressing were incubated for 48 h, and then, biological tests were performed, including cell viability (in propidium iodide staining), cell proliferation (in the SRB assay), evaluation of the intracellular free radical level (in the DCF-DA assay), genotoxicity (in the comet assay), assessment of the apoptotic and necrotic cells (in staining anexin-V and iodide propidium), the course of the cell cycle, and the scratch test. The correlation between apoptosis and genotoxicity and the levels of free radicals and genotoxicity were determined for the tested linen fibers and fabrics. The tests presented that the fibers are characterized by the ability to eliminate damaged cells in the elimination phase. However, the obtained fabrics gain different properties during the technological processing of the fibers into linen dressings. Linen fabrics have better regenerative properties for cells than fibers. The linseed dressing made of MB fiber has the most favorable regenerative properties.