Next Article in Journal
Colored LED Lights: Use One Color Alone or with Others for Growth in Hedyotis corymbosa In Vitro?
Previous Article in Journal
Photosynthetic Plasticity and Stomata Adjustment in Chromosome Segment Substitution Lines of Rice Cultivar KDML105 under Drought Stress
Previous Article in Special Issue
Development of Microsatellite Markers for Tanacetum cinerariifolium (Trevis.) Sch. Bip., a Plant with a Large and Highly Repetitive Genome
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Liubov V. Povkhova
1,
Elena N. Pushkova
1,
Tatiana A. Rozhmina
1,2,
Alexander A. Zhuchenko
2,3,
Roman I. Frykin
1,4,
Roman O. Novakovskiy
1,
Ekaterina M. Dvorianinova
1,5,
Aleksey A. Gryzunov
6,
Elena V. Borkhert
1,
Elizaveta A. Sigova
1,5,
Gleb N. Vladimirov
7,
Anastasiya V. Snezhkina
1,
Anna V. Kudryavtseva
1,
George S. Krasnov
1,
Alexey A. Dmitriev
1 and
Nataliya V. Melnikova
1,*
1
Engelhardt Institute of Molecular Biology, Russian Academy of Sciences, 119991 Moscow, Russia
2
Federal Research Center for Bast Fiber Crops, 172002 Torzhok, Russia
3
All-Russian Horticultural Institute for Breeding, Agrotechnology and Nursery, 115598 Moscow, Russia
4
Faculty of Biology, Lomonosov Moscow State University, 119234 Moscow, Russia
5
Moscow Institute of Physics and Technology, 141701 Moscow, Russia
6
All-Russian Scientific Research Institute of Refrigeration Industry—Branch of V.M. Gorbatov Federal Research Center for Food Systems of Russian Academy of Sciences, 127422 Moscow, Russia
7
Skolkovo Institute of Science and Technology, 121205 Moscow, Russia
*
Author to whom correspondence should be addressed.
Plants 2023, 12(1), 95; https://doi.org/10.3390/plants12010095
Submission received: 11 October 2022 / Revised: 14 December 2022 / Accepted: 16 December 2022 / Published: 24 December 2022
(This article belongs to the Special Issue Plant Molecular Breeding)

Abstract

:
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.

Graphical Abstract

1. Introduction

Flax (Linum usitatissimum L.) is one of the oldest cultivated plants, its seeds are traditionally used to produce linseed oil [1]. Flax seeds are one of the richest plant sources of healthy polyunsaturated linolenic acid and also contain a number of other valuable compounds (primarily lignans), which prevent the development of cardiovascular, oncological, and a number of other diseases [2,3,4,5,6,7,8,9,10,11]. Flax seeds are used in the production of paints, linoleum, composite materials, pharmaceuticals, food products, dietary supplements, and functional animal feed [3,6,12,13,14,15]. The main fatty acids of linseed oil are the unsaturated linolenic (LIN, ω-3), linoleic (LIO, ω-6), and oleic (OLE, ω-9) acids and saturated stearic and palmitic ones. The fatty acid composition of the oil differs significantly between varieties and determines the direction of the use of flax seeds [2,16,17,18]. Linseed oil with a high content of LIN (more than 50%) is used in the pharmaceutical and polymer industries, however, rapid oxidation limits its use in the food industry, where the seeds with a low (about 5%) and medium (30–40%) content of LIN are preferred [2,3,10,15,19,20].
For the development of flax lines with a low LIN content, mutagenesis was used. Mutant lines with a reduced content of LIN and increased content of LIO were obtained, and it was revealed that the LIN content was determined by the two genes with additive effect [21,22,23]. It is known that desaturases play a key role in the synthesis of oleic, linoleic, and linolenic acids and determine the fatty acid composition of plant oil [24]. Fatty acid desaturases 3 (FAD3) are responsible for the synthesis of LIN by the desaturation of LIO [25]. It was demonstrated that LuFAD3A and LuFAD3B are the major genes responsible for the conversion of LIO to LIN in flax, and mutations in these genes determine a low content of LIN in linseed oil [26].
Several studies aimed at finding the FAD3A and FAD3B mutations resulting in a low LIN content in linseed oil. Nonsense mutations in FAD3A (exon 5) and FAD3B (exon 1) resulted in a low content of LIN in line 593–708 [26]. Nonsense mutation in FAD3A (exon 2) and histidine to tyrosine substitution in FAD3B (exon 2) caused a low LIN content in line SP2047 [16]. Nonsense mutation in FAD3A (exon 5) and glycine to glutamic acid substitution in FAD3B (exon 5) probably prevented the conversion of LIO to LIN in the low-LIN variety TL23 [27]. Two nonsense mutations in FAD3A (exons 1 and 5) and a nonsense mutation and histidine to tyrosine substitution in FAD3B (exons 1 and 2 respectively) were revealed in flax varieties with a reduced LIN content [28]. Nonsense mutation in FAD3A (exon 1) and histidine to tyrosine substitution in FAD3B (exon 2) determined a low LIN content in the study on 40 flax lines [29]. In our previous studies on 279 flax varieties [18,30], three low-LIN mutations were identified: two nonsense mutations in FAD3A (exons 1 and 5) and histidine to tyrosine substitution in FAD3B (exon 2). Varieties with only one of the mutations had a medium content of LIN, while varieties with the mutations in both FAD3A and FAD3B had a low LIN content. Thus, inactivating mutations in the FAD3A and FAD3B genes responsible for a decrease in the content of LIN in linseed oil were identified in a significant number of studies.
Detection of low-LIN mutations in FAD3A and FAD3B can be useful for breeding flax varieties with a desired content of LIN for a particular application, especially for food purposes [20,26,28]. Marker-assisted selection is one of the most effective approaches for selecting valuable genotypes at the early stages of breeding [31]. It significantly increases the accuracy of breeding and reduces the time and material costs for creating varieties of many agricultural plants [32]. However, there is a lack of reliable test systems for the identification of low-LIN mutations in flax genotypes. Our study aimed to solve this issue.

2. Results

2.1. Selected Mutations of FAD3 Genes

The three most common mutations in the FAD3A and FAD3B genes that lead to a decrease in the LIN content and an increase in the LIO content in linseed oil were chosen for the development of test systems for their identification: G-to-A in exon 1 of FAD3A (CP027631.1: 16092348, GCA_000224295.2 ASM22429v2; abbreviated G-to-A/exon 1/FAD3A), resulting in a tryptophan to a stop codon substitution; C-to-T in exon 2 of FAD3B (CP027622.1: 1035655; abbreviated C-to-T/exon 2/FAD3B), resulting in a histidine to tyrosine substitution; and C-to-T in exon 5 of FAD3A (CP027631.1:16090340; abbreviated C-to-T/exon 5/FAD3A), resulting in an arginine to a stop codon substitution. Three approaches were proposed for the detection of these mutations in flax varieties and F3 hybrids: targeted deep sequencing, high-resolution melting (HRM), and cleaved amplified polymorphic sequences (CAPS) markers. The list of flax samples used in the work and approaches used for their analysis is presented in Supplementary Table S1.

2.2. Targeted Deep Sequencing

We used our previously proposed approach [18], which is based on the targeted deep sequencing, for the assessment of the presence/absence of the three mutations in FAD3A and FAD3B in 24 flax lines obtained by crossing flax varieties differing in oil fatty acid composition and in 19 varieties promising for the use in breeding, including those with a slightly reduced content of LIN (45–55%). G-to-A/exon 1/FAD3A was identified in line g.604 (about 25% of plants). C-to-T/exon 2/FAD3B was revealed in lines g.582 (about 65% of plants) and g.592, g.594, g.595, g.598, g.604 (about 25% of plants). C-to-T/exon 5/FAD3A was not found in any sample. Results of the detection of G-to-A/exon 1/FAD3A and C-to-T/exon 2/FAD3B are presented in Supplementary Table S2. Plants with both identified mutations are low-LIN and those with only one mutation are mid-LIN.
Targeted deep sequencing also enabled the screening of 94 F3 families from the Raciol × AGT 427/10 cross and 87 F3 families from the LM 98 × AGT 427/10 cross for the presence of G-to-A/exon 1/FAD3A and C-to-T/exon 2/FAD3B. Results are presented in Supplementary Table S2. F3 families that had homozygous mutations leading to a low or medium LIN content were chosen for further breeding flax varieties for the food industry: 44 F3 families for Raciol × AGT 427/10 with C-to-T/exon 2/FAD3B—medium LIN content; 36 F3 families for LM 98 × AGT 427/10 with G-to-A/exon 1/FAD3A or C-to-T/exon 2/FAD3B—medium LIN content; and 9 F3 families for LM 98 × AGT 427/10 with both G-to-A/exon 1/FAD3A and C-to-T/exon 2/FAD3B—low LIN content.
Moreover, we used targeted deep sequencing to detect G-to-A/exon 1/FAD3A and C-to-T/exon 2/FAD3B for 96 individual plants of line AGT 1535/07 and 96 individual plants of line AGT 987/02, for which high heterogeneity was revealed for the sites of these low-LIN mutations [18]. The results of the identification of the mutations are presented in Supplementary Table S2. The performed analysis enabled distributing the studied plants into groups with particular allelic variants of the FAD3A and FAD3B genes corresponding to a high, medium, or low LIN content and obtaining lines homogeneous by the LIN content for further use in practice and basic research.

2.3. HRM Analysis

Another approach we chose to identify mutations in the FAD3 genes was the HRM analysis [33]. We designed three primer pairs for the amplification of FAD3A and FAD3B regions containing G-to-A/exon 1/FAD3A, C-to-T/exon 2/FAD3B, and C-to-T/exon 5/FAD3A (sequences are presented in Supplementary Table S3). The primers were specific to the conserved DNA regions chosen on the basis of our previous studies on polymorphisms of the SAD and FAD genes in 279 flax varieties [18,30]. Six oligonucleotide standards were also developed, they corresponded to the variants of FAD3A and FAD3B with and without the analyzed mutations (Supplementary Table S3). PCR and HRM conditions were optimized. Validation of the developed test systems was carried out using oligonucleotide standards (Figure 1, Figure 2 and Figure 3) and individual plants of heterogeneous flax varieties, which were already subjected to the targeted deep sequencing analysis (AGT 1535/07 and AGT 987/02, 92 plants for each) (Figure 4, Figure 5 and Figure 6). Homozygous samples with and without the mutations and heterozygous samples were clearly distinguished from each other by the temperature shift and shape change of the DNA melting curves. The results of the mutation identification in the FAD3A and FAD3B genes obtained using deep sequencing and HRM analysis were in complete agreement. Thus, genotypes with the FAD3A and FAD3B allelic variants corresponding to a high, medium, or low LIN content were chosen from varieties AGT 1535/07 and AGT 987/02. HRM was also used for the selection of proper genotypes for further transcriptomic studies of heterogeneous varieties Lola, AGT 981/05, and AGT 1535/07 (40 individual plants for each variety). It enabled the identification and exclusion of non-typical plants from further research: 9 for Lola, 4 for AGT 981/05, and 15 for AGT 1535/07.

2.4. CAPS Markers

We also developed an approach for the identification of the three FAD3A and FAD3B mutations based on CAPS markers. This analysis includes the amplification of the target gene region with subsequent cleavage with restriction enzymes and evaluation of the DNA fragment lengths by gel electrophoresis [34]. We designed primers for the amplification of FAD3A and FAD3B regions containing the studied mutations (sequences are presented in Supplementary Table S3). As in the HRM analysis, primers were specific to the conserved regions of flax DNA. After PCR, in the absence of G-to-A/exon 1/FAD3A and C-to-T/exon 2/FAD3B (high LIN content), the resulting amplicons contained the HaeIII restriction site (GG∨CC/CC^GG) at the sites of these mutations. Amplicons after the restriction were separated by gel electrophoresis, and two DNA fragments were visualized in each case. However, in the presence of the studied mutations, the HaeIII restriction site in the amplicons was absent, and only one fragment was observed in each case (Figure 7 and Figure 8). In the absence of C-to-T/exon 5/FAD3A, the resulting amplicon contained the BstMBI restriction site (∨GATC/CTAG^), and two DNA fragments were visualized. However, in the presence of this mutation, the BstMBI restriction site was absent in the amplicon, and only one fragment was observed (Figure 9). In the case of heterozygosity of the studied mutations, three DNA fragments were visualized, corresponding to the initial length of the amplicon and two of its parts after restriction (Figure 7, Figure 8 and Figure 9).
CAPS markers were used for the identification of G-to-A/exon 1/FAD3A and C-to-T/exon 2/FAD3B in 185 individual plants of line AGT 1535/07 and 162 individual plants of line AGT 987/02. In these lines, high heterogeneity in the sites carrying low-LIN mutations was previously revealed [18]. From the examined sample set, 96 plants of AGT 1535/07 and 96 plants of AGT 987/02 were already studied for the presence of mutations in the FAD3A and FAD3B genes by the targeted deep sequencing, and the results of the two approaches were concordant. The application of CAPS markers made it possible to select plants of lines AGT 1535/07 and AGT 987/02 with homozygous allelic variants of the FAD3A and FAD3B genes corresponding to a high, medium, or low LIN content. The obtained groups of genotypes provide a perspective for the further development of flax varieties for specialized use in industry, food, and pharmaceuticals.
The developed CAPS markers were also used for the assessment of the presence of G-to-A/exon 1/FAD3A, C-to-T/exon 2/FAD3B, and C-to-T/exon 5/FAD3A in 600 samples of flax varieties Lola, AGT 981/05, and AGT 1535/07, which we used for transcriptomic studies of the mechanisms responsible for the linseed oil composition. It was extremely important as non-typical genotypes could distort the results of the analysis and lead to incorrect conclusions. Typical plants of Lola should have one of the three studied mutations—C-to-T/exon 5/FAD3A—while plants of AGT 981/05 and AGT 1535/07 should have two mutations—G-to-A/exon 1/FAD3A and C-to-T/exon 2/FAD3B. As a result, 42 of 200 Lola samples, 20 of 200 AGT 981/05 samples, and 81 of 200 AGT 1535/07 samples were excluded from further analysis, as they did not have the necessary mutations, or heterozygosity was revealed. Thus, a representative set of flax genotypes was formed for the assessment of gene expression in plants with a reduced LIN content, as required by the planned research.

3. Discussion

The content of unsaturated fatty acids in oil is one of the most important characteristics of flax. In the classical breeding of linseed, estimation of the fatty acid composition of the oil is usually carried out by the chromatographic method [35]. However, this method has several disadvantages: it often requires a significant number of seeds for the analysis; the use of a pool of heterogeneous seeds leads to an incorrect determination of oil fatty acid composition for a variety; the method does not allow differentiation between homozygous and heterozygous genotypes. Therefore, the assessment of fatty acid composition is usually performed only at the later stages of flax breeding, which makes the creation of varieties with the desired content of LIN a complex and lengthy process. The detection of FAD3A and FAD3B mutations that lead to a decrease in the LIN content can solve the indicated problems. It will enable the analysis of individual plants using only a few leaves for DNA isolation, making it possible to perform other studies of this genotype and obtain its seeds. It will also allow distinguishing between the mutations in heterozygous and homozygous states. Often, the use of data on the FAD3 mutations is the only reliable way to reasonably select material for flax breeding or research. For example, the medium content of LIN in linseed oil can be determined by any of the several mutations in the FAD3A or FAD3B genes. By phenotype, it is impossible to resolve which mutation is present in a variety. Choosing mid-LIN varieties for crossing to create low-LIN varieties, it is necessary that one parent carries the mutation in FAD3A and the other in FAD3B. This can be discovered from the sequences of these genes but not from the content of LIN in the oil. In the case of basic research, data on which mutations determine a reduced LIN content are necessary to create a representative set of varieties. This also cannot be carried out based on the analysis of the oil fatty acid composition.
For the identification of FAD3 mutations that result in a low LIN content in linseed oil, only CAPS markers were previously used [26,29,36]. Perhaps, the CAPS analysis was the only method because it is easy-to-use and can be applied in a wide range of laboratories, including those where high-throughput sequencers and real-time PCR systems are unavailable [37,38]. CAPS markers were designed for the detection of nonsense mutations in FAD3A (exon 5) and FAD3B (exon 1) [26], nonsense mutation in FAD3A (exon 5) [39], nonsense mutations in FAD3A (exon 1) and histidine to tyrosine substitution in FAD3B (exon 2) [29]. However, in the study [39], authors failed to identify low-LIN varieties using the developed CAPS markers. Probably, the studied varieties had other low-LIN mutations. In the present work, we developed CAPS markers for three low-LIN mutations of flax, which we identified earlier in the valuable-for-breeding flax varieties: G-to-A/exon 1/FAD3A, C-to-T/exon 2/FAD3B, and C-to-T/exon 5/FAD3A [18,30]. In developing the test system, we took into account the data on polymorphisms in about 300 flax varieties [18,30]. This enabled us to design primers that were specific to the conserved regions of flax DNA. Therefore, they are suitable for a wide range of flax genotypes. In addition, we designed primers in such a way that the amplified fragments were relatively short in length—250–300 bp. The amplification of longer fragments, as suggested in the works [26,29,36], can increase the requirements for DNA quality and PCR conditions. Thus, the developed CAPS markers can be applied in the analysis of diverse flax genotypes and allows the use of even low-quality DNA. Moreover, we developed two throughput approaches for the identification of the three low-LIN mutations in flax. They are based on the HRM analysis and targeted deep sequencing.
HRM analysis is widely used to identify DNA polymorphisms and genotype samples [40,41,42,43]. HRM does not require the use of expensive fluorophore-labeled SNP-specific probes in contrast to TaqMan assays [44,45]. The cost of the HRM analysis is usually slightly lower than that of another widely used fluorescence-based method—kompetitive allele specific PCR (KASP) [46,47,48,49,50]. Therefore, the use of HRM for plant genotyping could be a reasonable approach [51]. In our work, we developed the HRM-based test system for identifying the three mutations in the FAD3A and FAD3B genes that determine the content of LIN. HRM has a lower risk of contamination compared to CAPS since the analysis is carried out on real-time PCR platforms and there is no need to open tubes with amplicons. HRM is also faster than CAPS, it allows the assessment of a large sample set in a short time (only few hours are required).
High-throughput sequencing is a powerful tool for studying genetic diversity, developing molecular markers, and introducing new methods to plant breeding [52,53,54,55,56,57]. This technology is actively used in molecular genetic studies of flax using various approaches [58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89,90,91,92]. Along with microarrays [93,94,95], this method is the most throughput, but the second-generation sequencing (Illumina, BGI, and other platforms [96]) is much more precise than microarrays. Targeted deep sequencing is based on the enrichment of the regions of interest (PCR with target-specific primers in our case) and sequencing of a large number of samples in one run. This approach significantly reduces the cost of sequencing per genotype compared to the whole genome analysis and allows the sequencing of dozens and hundreds of genome regions in hundreds and thousands of samples. It also overcomes the difficulties of Sanger sequencing associated with polyploidy of plant genomes and a high degree of duplication of their genes [18,97,98,99,100,101]. Therefore, the third approach, developed by us, for the identification of the three low-LIN mutations in the FAD3A and FAD3B genes in flax genotypes was based on targeted deep sequencing. Compared to the HRM analysis and CAPS markers, the main advantage of targeted deep sequencing for the identification of the FAD3 mutations was its high accuracy and the ability to assess sample heterogeneity. For this method, a sample can be either an individual plant or a pool of dozens or hundreds of plants of the same variety.
In our work, the costs of HRM-based and CAPS-based approaches were comparable, but the targeted deep sequencing was several times more expensive per sample. The advantage of HRM and CAPS is the capacity to analyze only a few samples at a time. Targeted deep sequencing has an unmatched accuracy, it is designed for large sample sets and allows assessing many polymorphisms at once. In the case of screening studies of several thousand samples in one run, the cost of the analysis per sample performed using targeted deep sequencing will be close to those performed using HRM and CAPS.
The complex application of the three proposed approaches for the identification of the FAD3 mutations makes it possible to increase the efficiency of flax breeding and the accuracy of scientific research on this species. It is advisable to start the analysis of a collection of flax varieties of interest with targeted deep sequencing. For each variety, a DNA pool of several dozens or hundreds of plants should be analyzed to determine the presence of the FAD3 gene mutations and assess the varietal heterogeneity (the mutations may be present only in a subset of plants of variety). On the basis of the obtained data, it is possible to select the most suitable genotypes for further work. These can be promising parental pairs for crossing, heterogeneous varieties whose purity needs to be improved, and genotypes valuable for basic research. Moreover, targeted deep sequencing could be effective to develop DNA certificates of varieties and control the varietal purity when maintaining collections and seed production of flax.
The use of HRM and CAPS markers is advisable for choosing individual plants in the marker-assisted selection of linseed, as well as plants with a desired genotype for basic research. The HRM and CAPS approaches enable the identification of the three most common low-LIN mutations in the FAD3 genes already in F2 hybrids (and hybrids of later generations) and the prediction of the LIN content in seeds of these plants. Therefore, even at the early stages of breeding, it is possible to select genotypes with desirable allelic variants of the FAD3A and FAD3B genes in the homozygous state and continue working only with these genotypes. In addition, the use of HRM and CAPS markers is effective for increasing the purity of varieties for which heterogeneity in the FAD3 gene mutations was identified.
Thus, based on the targeted deep sequencing, HRM analysis, and CAPS markers, we developed reliable approaches for the identification of the three most common mutations in the FAD3A and FAD3B genes that determine a reduced content of LIN in linseed oil: G-to-A/exon 1/FAD3A, C-to-T/exon 2/FAD3B, and C-to-T/exon 5/FAD3A. The choice of one of the methods or their combinations depends on the task being solved, the material being assessed, and the availability of the equipment. Using different combinations of the developed approaches, we successfully analyzed more than a thousand of flax samples and selected plants with the desired stably inherited LIN and LIO content for breeding flax varieties for specialized purposes, as well as studying molecular mechanisms related to the fatty acid synthesis.

4. Materials and Methods

4.1. Plant Material

The following flax samples from the collection of the Institute for Flax (Torzhok, Russia) were used in the work (the list of samples and methods is presented in Supplementary Table S1):
  • A total of 24 lines (g.579, g.580, g.582, g.586, g.587, g.588, g.591, g.592, g.594, g.595, g.597, g.598, g.602, g.603, g.604, g.605, g.610, g.614, g.615, g.616, g.618, g.619, g.620, g.621) obtained from crossing flax varieties differing in the fatty acid composition of the oil and 19 varieties promising for use in breeding, including those with a slightly reduced LIN content (45–55%): Hi de France, L.8709-5-10, N 3860, B-15, N 3895, K-74, Moolzyzan, Cki 10, NP RR, Soletsky kr, S Yanyshina A.A., T Yanyshina A.A., N 3810, N 3841, Stormont Motley, N 3809, Istok, Honkie 49, Focus. For each sample, material was obtained from pools of 50 or more plants;
  • Families of F3 hybrids of flax obtained from crosses Raciol × AGT 427/10 (94 F3 families) and LM 98 × AGT 427/10 (87 F3 families) with a contrast ratio of LIN and LIO in linseed oil and different allelic variants of the FAD3A and FAD3B genes. For each F3 family, material was obtained from a pool of at least 20 plants;
  • Varieties AGT 1535/07 (185 individual plants) and AGT 987/02 (162 individual plants) for which high heterogeneity was revealed for sites of low-LIN mutations [18]. Material was obtained from each individual plant;
  • Varieties Lola, AGT 981/05, and AGT 1535/07 (200 samples for each variety), which we used for transcriptomic studies of the mechanisms of linseed oil synthesis and for which some heterogeneity was revealed for sites of the low-LIN mutations [18]. Material was obtained from each individual plant.

4.2. DNA Isolation

Leaves were used for DNA extraction. They were ground using a MagNA Lyser (Roche, Basel, Switzerland) or TissueLyser II (Qiagen, Chatsworth, CA, USA) homogenizers. Thereafter, the CTAB method (in case of pools of several plants) or DNA-Extran-3 kit (Syntol, Moscow, Russia) (in case of individual plants) were used for DNA extraction. DNA quality was assessed by electrophoresis in 1% agarose gel. DNA concentration was determined with a Qubit fluorometer (Thermo Fisher Scientific, Waltham, MA, USA) or a NanoDrop spectrophotometer (Thermo Fisher Scientific).

4.3. Targeted Deep Sequencing

To prepare DNA libraries for the targeted deep sequencing, an approach based on two sequential PCR was used: the first PCR was to amplify selected gene regions and add universal adapter sequences to them, and the second PCR was to add sequences required for high-throughput sequencing and barcodes for sample identification [18]. For amplification of sites containing the G-to-A mutation in exon 1 of FAD3A (CP027631.1: 16092348, GCA_000224295.2 ASM22429v2; tryptophan to a stop codon substitution), C-to-T mutation in exon 2 of FAD3B (CP027622.1: 1035655; histidine to tyrosine substitution), and C-to-T mutation in exon 5 of FAD3A (CP027631.1: 16090340; arginine to a stop codon substitution), we used primers, as proposed by us earlier (sequences are presented in Supplementary Table S3), and a protocol for DNA library preparation for targeted deep sequencing [18]. One primer pair enabled simultaneous amplification of the FAD3A region carrying the G-to-A mutation and FAD3B region carrying the C-to-T mutation, while another primer pair amplified the FAD3A region with the C-to-T mutation. The quality and concentration of the obtained DNA libraries were evaluated on a 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA) and Qubit fluorometer (Thermo Fisher Scientific). The resulting pools of DNA libraries were sequenced on MiSeq (Illumina, San Diego, CA, USA) using the MiSeq Reagent Kit v3 (600 cycles). For deep sequencing data processing, reads were trimmed and filtered using Trimmomatic [102], mapped to the reference genome of flax variety Bethune (GCA_000224295.2 ASM22429v2) using BWA-MEM [103]. Mapped reads were grouped using Picard, and polymorphisms were searched in the regions of interest using the FreeBayes tool [104].

4.4. HRM Analysis

For the HRM analysis of FAD3A and FAD3B gene regions containing the sites of the three studied mutations determining the content of LIN and LIO (G-to-A in exon 1 of FAD3A (CP027631.1: 16092348, GCA_000224295.2 ASM22429v2), C-to-T in exon 2 of FAD3B (CP027622.1: 1035655), and C-to-T in exon 5 of FAD3A (CP027631.1: 16090340)), primers (Supplementary Table S3) were designed for the conserved regions taking into account the data on polymorphisms obtained by us earlier for 279 flax varieties [18,30]. Reaction mix was prepared in a volume of 20 μL for each sample and standard as follows: 1× Taq Turbo buffer (Evrogen, Moscow, Russia), 250 nM of dNTPs (Evrogen), 350 nM of each primer, 1× EvaGreen dye (Biotium, Fremont, CA, USA), 1× LowROX dye (Evrogen), 2 units of HS Taq DNA polymerase (Evrogen). DNA concentration was about 5 ng/µL for samples and 5 fg/µL for standards. Synthesized oligonucleotides were used as standards (with and without each of the three analyzed mutations, Supplementary Table S3). The following amplification program was used: 95 °C—10 min; 50 cycles of 95 °C—15 s, 60 °C—1 min. Each reaction was performed in triplicate. Then, the HRM analysis was carried out: 65–90 °C, step 0.1 °C. Data processing was performed using the QuantStudio 5 (Thermo Fisher Scientific) software. For each low-LIN mutation, the obtained melting curves for the studied flax samples were compared with the melting curves for two standards (with and without the mutation) and a 1:1 mixture of the two standards. The similarity between the melting curves of the examined sample and the standard carrying the mutation indicated the presence of the mutation in the studied genotype in the homozygous state. The similarity between the melting curves of the examined sample and the standard without the mutation indicated the absence of the mutation in the studied genotype. The similarity between the melting curves of the examined sample and the 1:1 mixture of two standards indicated the presence of the mutation in the studied genotype in the heterozygous state.

4.5. CAPS Markers

To conduct the CAPS analysis of FAD3A and FAD3B gene regions containing the sites of three studied mutations that determine the content of LIN and LIO, primers (Supplementary Table S3) were designed for the conserved regions of these genes, taking into account the data on polymorphisms obtained by us earlier for 279 flax varieties [18,30]. For amplification in a volume of 15 µL, the following reaction mixture was used: 1× Taq Turbo buffer (Evrogen), 250 nM of dNTPs (Evrogen), 350 nM of each primer, 1 unit of HS Taq DNA polymerase (Eurogen), 20 ng of DNA. The following amplification program was used: 95 °C—3 min; 40 cycles of 95 °C—30 s, 60 °C—30 s, 72 °C—60 s; 72 °C—3 min. Then, restriction was performed for 4 μL of each PCR product. In the case of the G-to-A mutation in exon 1 of FAD3A (CP027631.1: 16092348, GCA_000224295.2 ASM22429v2) and C-to-T mutation in exon 2 of FAD3B (CP027622.1: 1035655), incubation with 10 units of HaeIII (SibEnzyme, Novosibirsk, Russia, restriction site GG∨CC/CC^GG) was performed at 37 °C for 4 h followed by inactivation of the restriction enzyme by heating to 80 °C for 20 min. In the case of the C-to-T mutation in exon 5 of FAD3A (CP027631.1: 16090340), incubation with 10 units of BstMBI (SibEnzyme, restriction site ∨GATC/CTAG^) at 65 °C for 4 h followed by inactivation of the restriction enzyme by heating to 80 °C for 20 min was performed. Restriction products were visualized by 2% agarose gel electrophoresis. In the absence of the G-to-A mutation in exon 1 of FAD3A, the lengths of the fragments were 195 and 107 bp, while in the presence of the mutation, the fragment length was 302 bp. In the absence of the C-to-T mutation in exon 2 of FAD3B, the lengths of the fragments were 157 and 103 bp, while in the presence of the mutation, the fragment length was 260 bp. In the absence of the C-to-T mutation in exon 5 of FAD3A, the lengths of the fragments were 163 and 109 bp, and in the presence of the mutation, the fragment length was 272 bp.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants12010095/s1, Table S1: Flax plant material from the collection of the Institute for Flax that was used for the identification of the G to A mutation in exon 1 of FAD3A, C to T mutation in exon 2 of FAD3B, and C to T mutation in exon 5 of FAD3A using targeted deep sequencing, high resolution melting (HRM) analysis, and cleaved amplified polymorphic sequences (CAPS) markers. Table S2: Results of the identification of the G to A mutation in exon 1 of FAD3A and C to T mutation in exon 2 of FAD3B in flax samples using targeted deep sequencing. Table S3: Sequences of primers and oligonucleotide standards for targeted deep sequencing, HRM analysis, and CAPS markers developed for the identification of the G to A mutation in exon 1 of FAD3A, C to T mutation in exon 2 of FAD3B, and C to T mutation in exon 5 of FAD3A.

Author Contributions

Conceptualization, A.A.D. and N.V.M.; performing experiments, L.V.P., E.N.P., T.A.R., R.I.F., R.O.N., E.A.S., A.A.G., E.V.B. and A.V.S.; data analysis, L.V.P., A.A.Z., R.I.F., E.M.D., G.N.V., A.V.K., G.S.K., A.A.D. and N.V.M.; writing, L.V.P., E.M.D., A.A.D. and N.V.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the Russian Science Foundation, grant 21-16-00111.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw sequencing data have been deposited in the NCBI Sequence Read Archive (SRA) under the BioProject accession number PRJNA625974.

Acknowledgments

We thank the Center for Precision Genome Editing and Genetic Technologies for Biomedicine, EIMB RAS for providing the computing power and techniques for the data analysis. This work was performed using the equipment of EIMB RAS “Genome” center (http://www.eimb.ru/ru1/ckp/ccu_genome_ce.php, accessed on 14 December 2022).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Weiss, E.; Zohary, D.; Hopf, M. Domestication of Plants in the Old World—The Origin and Spread of Domesticated Plants in South-West Asia, Europe, and the Mediterranean Basin; Oxford University Press: Oxford, UK, 2012. [Google Scholar]
  2. Muir, A.D.; Westcott, N.D. Flax: The Genus Linum; CRC Press: Boca Raton, FL, USA, 2003. [Google Scholar]
  3. Goyal, A.; Sharma, V.; Upadhyay, N.; Gill, S.; Sihag, M. Flax and flaxseed oil: An ancient medicine & modern functional food. J. Food Sci. Technol. 2014, 51, 1633–1653. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Imran, M.; Ahmad, N.; Anjum, F.M.; Khan, M.K.; Mushtaq, Z.; Nadeem, M.; Hussain, S. Potential protective properties of flax lignan secoisolariciresinol diglucoside. Nutr. J. 2015, 14, 71. [Google Scholar] [CrossRef] [Green Version]
  5. Abdelhamid, A.S.; Martin, N.; Bridges, C.; Brainard, J.S.; Wang, X.; Brown, T.J.; Hanson, S.; Jimoh, O.F.; Ajabnoor, S.M.; Deane, K.H.; et al. Polyunsaturated fatty acids for the primary and secondary prevention of cardiovascular disease. Cochrane Database Syst. Rev. 2018, 7, CD012345. [Google Scholar] [CrossRef] [PubMed]
  6. Campos, J.R.; Severino, P.; Ferreira, C.S.; Zielinska, A.; Santini, A.; Souto, S.B.; Souto, E.B. Linseed Essential Oil-Source of Lipids as Active Ingredients for Pharmaceuticals and Nutraceuticals. Curr. Med. Chem. 2019, 26, 4537–4558. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Kezimana, P.; Dmitriev, A.A.; Kudryavtseva, A.V.; Romanova, E.V.; Melnikova, N.V. Secoisolariciresinol Diglucoside of Flaxseed and Its Metabolites: Biosynthesis and Potential for Nutraceuticals. Front. Genet. 2018, 9, 641. [Google Scholar] [CrossRef] [Green Version]
  8. Locke, A.; Schneiderhan, J.; Zick, S.M. Diets for Health: Goals and Guidelines. Am. Fam. Physician 2018, 97, 721–728. [Google Scholar]
  9. Parikh, M.; Netticadan, T.; Pierce, G.N. Flaxseed: Its bioactive components and their cardiovascular benefits. Am. J. Physiol. Heart Circ. Physiol. 2018, 314, H146–H159. [Google Scholar] [CrossRef]
  10. Cullis, C.A. Genetics and Genomics of Linum; Springer International Publishing: Berlin/Heidelberg, Germany, 2019. [Google Scholar]
  11. Mali, A.V.; Padhye, S.B.; Anant, S.; Hegde, M.V.; Kadam, S.S. Anticancer and antimetastatic potential of enterolactone: Clinical, preclinical and mechanistic perspectives. Eur. J. Pharmacol. 2019, 852, 107–124. [Google Scholar] [CrossRef]
  12. Kouba, M.; Mourot, J. A review of nutritional effects on fat composition of animal products with special emphasis on n-3 polyunsaturated fatty acids. Biochimie 2011, 93, 13–17. [Google Scholar] [CrossRef]
  13. Singh, K.K.; Mridula, D.; Rehal, J.; Barnwal, P. Flaxseed: A Potential Source of Food, Feed and Fiber. Crit. Rev. Food Sci. Nutr. 2011, 51, 210–222. [Google Scholar] [CrossRef]
  14. Corino, C.; Rossi, R.; Cannata, S.; Ratti, S. Effect of dietary linseed on the nutritional value and quality of pork and pork products: Systematic review and meta-analysis. Meat Sci. 2014, 98, 679–688. [Google Scholar] [CrossRef] [PubMed]
  15. Fombuena, V.; Petrucci, R.; Dominici, F.; Jorda-Vilaplana, A.; Montanes, N.; Torre, L. Maleinized Linseed Oil as Epoxy Resin Hardener for Composites with High Bio Content Obtained from Linen Byproducts. Polymers 2019, 11, 301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Banik, M.; Duguid, S.; Cloutier, S. Transcript profiling and gene characterization of three fatty acid desaturase genes in high, moderate, and low linolenic acid genotypes of flax (Linum usitatissimum L.) and their role in linolenic acid accumulation. Genome 2011, 54, 471–483. [Google Scholar] [CrossRef] [PubMed]
  17. Cloutier, S.; Ragupathy, R.; Niu, Z.; Duguid, S. SSR-based linkage map of flax (Linum usitatissimum L.) and mapping of QTLs underlying fatty acid composition traits. Mol. Breed. 2011, 28, 437–451. [Google Scholar] [CrossRef]
  18. Dmitriev, A.A.; Kezimana, P.; Rozhmina, T.A.; Zhuchenko, A.A.; Povkhova, L.V.; Pushkova, E.N.; Novakovskiy, R.O.; Pavelek, M.; Vladimirov, G.N.; Nikolaev, E.N.; et al. Genetic diversity of SAD and FAD genes responsible for the fatty acid composition in flax cultivars and lines. BMC Plant Biol. 2020, 20, 301. [Google Scholar] [CrossRef] [PubMed]
  19. Chen, Y.; Zhou, X.-R.; Zhang, Z.-J.; Dribnenki, P.; Singh, S.; Green, A. Development of high oleic oil crop platform in flax through RNAi-mediated multiple FAD2 gene silencing. Plant Cell Rep. 2015, 34, 643–653. [Google Scholar] [CrossRef]
  20. Brutch, N.; Porokhoviniva, E.; Shelenga, T. Perspectives of the creation of oil flax varieties for the specialized purpose. Agrar. Report. South-East 2016, 1–2, 50–52. [Google Scholar]
  21. Green, A.G.; Marshall, D.R. Isolation of induced mutants in linseed (Linum usitatissimum) having reduced linolenic acid content. Euphytica 1984, 33, 321–328. [Google Scholar] [CrossRef]
  22. Rowland, G. An EMS-induced low-linolenic-acid mutant in McGregor flax (Linum usitatissimum L.). Can. J. Plant Sci. 1991, 71, 393–396. [Google Scholar] [CrossRef] [Green Version]
  23. Green, A.G. Genetic control of polyunsaturated fatty acid biosynthesis in flax (Linum usitatissimum) seed oil. Theor. Appl. Genet. 1986, 72, 654–661. [Google Scholar] [CrossRef]
  24. Shanklin, J.; Cahoon, E.B. Desaturation and Related Modifications of Fatty Acids1. Annu. Rev. Plant Physiol. Plant Mol. Biol. 1998, 49, 611–641. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Fulco, A.J. Metabolic alterations of fatty acids. Annu. Rev. Biochem. 1974, 43, 215–241. [Google Scholar] [CrossRef] [PubMed]
  26. Vrinten, P.; Hu, Z.; Munchinsky, M.A.; Rowland, G.; Qiu, X. Two FAD3 desaturase genes control the level of linolenic acid in flax seed. Plant Physiol. 2005, 139, 79–87. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Rajwade, A.V.; Joshi, R.S.; Kadoo, N.Y.; Gupta, V.S. Sequence characterization and in silico structure prediction of fatty acid desaturases in linseed varieties with differential fatty acid composition. J. Sci. Food Agric. 2016, 96, 4896–4906. [Google Scholar] [CrossRef]
  28. Thambugala, D.; Duguid, S.; Loewen, E.; Rowland, G.; Booker, H.; You, F.M.; Cloutier, S. Genetic variation of six desaturase genes in flax and their impact on fatty acid composition. Theor. Appl. Genet. 2013, 126, 2627–2641. [Google Scholar] [CrossRef] [Green Version]
  29. Porokhovinova, E.A.; Shelenga, T.V.; Matveeva, T.V.; Pavlov, A.V.; Grigorieva, E.A.; Brutch, N.B. Polymorphism of genes controlling low level of linolenic acid in lines from VIR flax genetic collection. Ecol. Genet. 2019, 17, 5–19. [Google Scholar] [CrossRef] [Green Version]
  30. Kezimana, P.; Rozhmina, T.A.; Krasnov, G.S.; Povkhova, L.V.; Novakovskiy, R.O.; Pushkova, E.N.; Zhuchenko, A.A.; Bjelkova, M.; Pavelek, M.; Dmitriev, A.A.; et al. Evaluation of polymorphism of SAD and FAD genes in flax (Linum usitatissimum L.) cultivars and lines using deep sequencing. In Theory and Practice of Adaptive Plant Breeding, Proceedings of the International Scientific and Practical Conference, Krasnodar, Russia, 25 September 2020; Kuban State Agrarian University: Krasnodar, Russia, 2020; pp. 49–52. [Google Scholar]
  31. Sonkar, P.; Sharma, V.; Kumar, H.; Kumar, M.; Negi, H.; Kiran, U.; Abdin, M.; Choubey, A. Marker-Assisted Selection for Value Addition in Crop Plants. In Technologies in Plant Biotechnology and Breeding of Field Crops; Springer: Berlin/Heidelberg, Germany, 2022; pp. 23–39. [Google Scholar]
  32. Hasan, N.; Choudhary, S.; Naaz, N.; Sharma, N.; Laskar, R.A. Recent advancements in molecular marker-assisted selection and applications in plant breeding programmes. J. Genet. Eng. Biotechnol. 2021, 19, 128. [Google Scholar] [CrossRef]
  33. Reed, G.H.; Kent, J.O.; Wittwer, C.T. High-resolution DNA melting analysis for simple and efficient molecular diagnostics. Pharmacogenomics 2007, 8, 597–608. [Google Scholar] [CrossRef] [Green Version]
  34. Konieczny, A.; Ausubel, F.M. A procedure for mapping Arabidopsis mutations using co-dominant ecotype-specific PCR-based markers. Plant J. Cell Mol. Biol. 1993, 4, 403–410. [Google Scholar] [CrossRef]
  35. Bjelkova, M.; Nozkova, J.; Fatrcova-Sramkova, K.; Tejklova, E. Comparison of linseed (Linum usitatissimum L.) genotypes with respect to the content of polyunsaturated fatty acids. Chem. Pap. 2012, 66, 972–976. [Google Scholar] [CrossRef]
  36. Matuszczak, M.; Spasibionek, S.; Gacek, K.; Bartkowiak-Broda, I. Cleaved amplified polymorphic sequences (CAPS) marker for identification of two mutant alleles of the rapeseed BnaA. FAD2 gene. Mol. Biol. Rep. 2020, 47, 7607–7621. [Google Scholar] [CrossRef] [PubMed]
  37. Shavrukov, Y. CAPS markers in plant biology. Russ. J. Genet. Appl. Res. 2016, 6, 279–287. [Google Scholar] [CrossRef]
  38. Shavrukov, Y. Cleaved Amplified Polymorphic Sequences (CAPS) Markers in Plant Biology; Nova Science Publishers, Inc.: Hauppauge, NY, USA, 2014. [Google Scholar]
  39. Walkowiak, M.; Matuszczak, M.; Spasibionek, S.; Liersch, A.; Mikołajczyk, K. Cleaved Amplified Polymorphic Sequences (CAPS) Markers for Characterization of the LuFAD3A Gene from Various Flax (Linum usitatissimum L.) Cultivars. Agronomy 2022, 12, 1432. [Google Scholar] [CrossRef]
  40. Simko, I. High-Resolution DNA Melting Analysis in Plant Research. Trends Plant Sci. 2016, 21, 528–537. [Google Scholar] [CrossRef]
  41. Tamburro, M.; Ripabelli, G. High Resolution Melting as a rapid, reliable, accurate and cost-effective emerging tool for genotyping pathogenic bacteria and enhancing molecular epidemiological surveillance: A comprehensive review of the literature. Ann. Ig. Med. Prev. Comunita 2017, 29, 293–316. [Google Scholar] [CrossRef]
  42. Druml, B.; Cichna-Markl, M. High resolution melting (HRM) analysis of DNA--its role and potential in food analysis. Food Chem. 2014, 158, 245–254. [Google Scholar] [CrossRef]
  43. Montgomery, J.L.; Sanford, L.N.; Wittwer, C.T. High-resolution DNA melting analysis in clinical research and diagnostics. Expert Rev. Mol. Diagn. 2010, 10, 219–240. [Google Scholar] [CrossRef]
  44. Broccanello, C.; Chiodi, C.; Funk, A.; McGrath, J.M.; Panella, L.; Stevanato, P. Comparison of three PCR-based assays for SNP genotyping in plants. Plant Methods 2018, 14, 28. [Google Scholar] [CrossRef] [Green Version]
  45. Ayalew, H.; Tsang, P.W.; Chu, C.; Wang, J.; Liu, S.; Chen, C.; Ma, X.F. Comparison of TaqMan, KASP and rhAmp SNP genotyping platforms in hexaploid wheat. PLoS ONE 2019, 14, e0217222. [Google Scholar] [CrossRef] [Green Version]
  46. Semagn, K.; Babu, R.; Hearne, S.; Olsen, M. Single nucleotide polymorphism genotyping using Kompetitive Allele Specific PCR (KASP): Overview of the technology and its application in crop improvement. Mol. Breed. 2014, 33, 1–14. [Google Scholar] [CrossRef]
  47. Shavrukov, Y. Comparison of SNP and CAPS markers application in genetic research in wheat and barley. BMC Plant Biol. 2016, 16 (Suppl. S1), 11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Ali, A.; Altaf, M.T.; Nadeem, M.A.; Karakoy, T.; Shah, A.N.; Azeem, H.; Baloch, F.S.; Baran, N.; Hussain, T.; Duangpan, S.; et al. Recent advancement in OMICS approaches to enhance abiotic stress tolerance in legumes. Front. Plant Sci. 2022, 13, 952759. [Google Scholar] [CrossRef]
  49. Sahin, Z.N.; Sahin, E.C.; Aydin, Y.; Uncuoglu, A.A. Molecular sexual determinants in Pistacia genus by KASP assay. Mol. Biol. Rep. 2022, 49, 5473–5482. [Google Scholar] [CrossRef]
  50. Makhoul, M.; Obermeier, C. Development of Breeder-Friendly KASP Markers from Genome-Wide Association Studies Results. Methods Mol. Biol. 2022, 2481, 287–310. [Google Scholar] [CrossRef] [PubMed]
  51. Chou, L.; Huang, S.J.; Hsieh, C.; Lu, M.T.; Song, C.W.; Hsu, F.C. A High Resolution Melting Analysis-Based Genotyping Toolkit for the Peach (Prunus persica) Chilling Requirement. Int. J. Mol. Sci. 2020, 21, 1543. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. McCormack, J.E.; Faircloth, B.C. Next-generation phylogenetics takes root. Mol. Ecol. 2013, 22, 19–20. [Google Scholar] [CrossRef] [Green Version]
  53. Unamba, C.I.; Nag, A.; Sharma, R.K. Next Generation Sequencing Technologies: The Doorway to the Unexplored Genomics of Non-Model Plants. Front. Plant Sci. 2015, 6, 1074. [Google Scholar] [CrossRef] [Green Version]
  54. He, J.; Zhao, X.; Laroche, A.; Lu, Z.X.; Liu, H.; Li, Z. Genotyping-by-sequencing (GBS), an ultimate marker-assisted selection (MAS) tool to accelerate plant breeding. Front. Plant Sci. 2014, 5, 484. [Google Scholar] [CrossRef] [Green Version]
  55. Sarwat, M.; Yamdagni, M.M. DNA barcoding, microarrays and next generation sequencing: Recent tools for genetic diversity estimation and authentication of medicinal plants. Crit. Rev. Biotechnol. 2016, 36, 191–203. [Google Scholar] [CrossRef]
  56. Bohra, A.; Jha, U.C.; Kishor, P.B.; Pandey, S.; Singh, N.P. Genomics and molecular breeding in lesser explored pulse crops: Current trends and future opportunities. Biotechnol. Adv. 2014, 32, 1410–1428. [Google Scholar] [CrossRef]
  57. Kilian, B.; Graner, A. NGS technologies for analyzing germplasm diversity in genebanks. Brief. Funct. Genom. 2012, 11, 38–50. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Fu, Y.B.; Peterson, G.W. Developing genomic resources in two Linum species via 454 pyrosequencing and genomic reduction. Mol. Ecol. Resour. 2012, 12, 492–500. [Google Scholar] [CrossRef] [PubMed]
  59. Kumar, S.; You, F.M.; Cloutier, S. Genome wide SNP discovery in flax through next generation sequencing of reduced representation libraries. BMC Genom. 2012, 13, 684. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Wang, Z.; Hobson, N.; Galindo, L.; Zhu, S.; Shi, D.; McDill, J.; Yang, L.; Hawkins, S.; Neutelings, G.; Datla, R.; et al. The genome of flax (Linum usitatissimum) assembled de novo from short shotgun sequence reads. Plant J. Cell Mol. Biol. 2012, 72, 461–473. [Google Scholar] [CrossRef] [Green Version]
  61. Yu, Y.; Huang, W.; Chen, H.; Wu, G.; Yuan, H.; Song, X.; Kang, Q.; Zhao, D.; Jiang, W.; Liu, Y.; et al. Identification of differentially expressed genes in flax (Linum usitatissimum L.) under saline-alkaline stress by digital gene expression. Gene 2014, 549, 113–122. [Google Scholar] [CrossRef]
  62. Galindo-Gonzalez, L.; Pinzon-Latorre, D.; Bergen, E.A.; Jensen, D.C.; Deyholos, M.K. Ion Torrent sequencing as a tool for mutation discovery in the flax (Linum usitatissimum L.) genome. Plant Methods 2015, 11, 19. [Google Scholar] [CrossRef] [Green Version]
  63. Melnikova, N.V.; Dmitriev, A.A.; Belenikin, M.S.; Speranskaya, A.S.; Krinitsina, A.A.; Rachinskaia, O.A.; Lakunina, V.A.; Krasnov, G.S.; Snezhkina, A.V.; Sadritdinova, A.F.; et al. Excess fertilizer responsive miRNAs revealed in Linum usitatissimum L. Biochimie 2015, 109, 36–41. [Google Scholar] [CrossRef]
  64. Dmitriev, A.A.; Krasnov, G.S.; Rozhmina, T.A.; Kishlyan, N.V.; Zyablitsin, A.V.; Sadritdinova, A.F.; Snezhkina, A.V.; Fedorova, M.S.; Yurkevich, O.Y.; Muravenko, O.V.; et al. Glutathione S-transferases and UDP-glycosyltransferases Are Involved in Response to Aluminum Stress in Flax. Front. Plant Sci. 2016, 7, 1920. [Google Scholar] [CrossRef] [Green Version]
  65. Dmitriev, A.A.; Kudryavtseva, A.V.; Krasnov, G.S.; Koroban, N.V.; Speranskaya, A.S.; Krinitsina, A.A.; Belenikin, M.S.; Snezhkina, A.V.; Sadritdinova, A.F.; Kishlyan, N.V.; et al. Gene expression profiling of flax (Linum usitatissimum L.) under edaphic stress. BMC Plant Biol. 2016, 16, 237. [Google Scholar] [CrossRef] [Green Version]
  66. Fu, Y.B.; Dong, Y.; Yang, M.H. Multiplexed shotgun sequencing reveals congruent three-genome phylogenetic signals for four botanical sections of the flax genus Linum. Mol. Phylogenetics Evol. 2016, 101, 122–132. [Google Scholar] [CrossRef]
  67. Galindo-Gonzalez, L.; Deyholos, M.K. RNA-seq Transcriptome Response of Flax (Linum usitatissimum L.) to the Pathogenic Fungus Fusarium oxysporum f. sp. lini. Front. Plant Sci. 2016, 7, 1766. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Melnikova, N.V.; Dmitriev, A.A.; Belenikin, M.S.; Koroban, N.V.; Speranskaya, A.S.; Krinitsina, A.A.; Krasnov, G.S.; Lakunina, V.A.; Snezhkina, A.V.; Sadritdinova, A.F.; et al. Identification, Expression Analysis, and Target Prediction of Flax Genotroph MicroRNAs Under Normal and Nutrient Stress Conditions. Front. Plant Sci. 2016, 7, 399. [Google Scholar] [CrossRef]
  69. Bolsheva, N.L.; Melnikova, N.V.; Kirov, I.V.; Speranskaya, A.S.; Krinitsina, A.A.; Dmitriev, A.A.; Belenikin, M.S.; Krasnov, G.S.; Lakunina, V.A.; Snezhkina, A.V.; et al. Evolution of blue-flowered species of genus Linum based on high-throughput sequencing of ribosomal RNA genes. BMC Evol. Biol. 2017, 17, 253. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Dash, P.K.; Rai, R.; Mahato, A.K.; Gaikwad, K.; Singh, N.K. Transcriptome Landscape at Different Developmental Stages of a Drought Tolerant Cultivar of Flax (Linum usitatissimum). Front. Chem. 2017, 5, 82. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Dmitriev, A.A.; Krasnov, G.S.; Rozhmina, T.A.; Novakovskiy, R.O.; Snezhkina, A.V.; Fedorova, M.S.; Yurkevich, O.Y.; Muravenko, O.V.; Bolsheva, N.L.; Kudryavtseva, A.V.; et al. Differential gene expression in response to Fusarium oxysporum infection in resistant and susceptible genotypes of flax (Linum usitatissimum L.). BMC Plant Biol. 2017, 17, 253. [Google Scholar] [CrossRef]
  72. Dmitriev, A.A.; Kudryavtseva, A.V.; Bolsheva, N.L.; Zyablitsin, A.V.; Rozhmina, T.A.; Kishlyan, N.V.; Krasnov, G.S.; Speranskaya, A.S.; Krinitsina, A.A.; Sadritdinova, A.F.; et al. miR319, miR390, and miR393 Are Involved in Aluminum Response in Flax (Linum usitatissimum L.). BioMed Res. Int. 2017, 2017, 4975146. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Yi, L.; Gao, F.; Siqin, B.; Zhou, Y.; Li, Q.; Zhao, X.; Jia, X.; Zhang, H. Construction of an SNP-based high-density linkage map for flax (Linum usitatissimum L.) using specific length amplified fragment sequencing (SLAF-seq) technology. PLoS ONE 2017, 12, e0189785. [Google Scholar] [CrossRef] [Green Version]
  74. Wu, J.; Zhao, Q.; Zhang, L.; Li, S.; Ma, Y.; Pan, L.; Lin, H.; Wu, G.; Yuan, H.; Yu, Y.; et al. QTL Mapping of Fiber-Related Traits Based on a High-Density Genetic Map in Flax (Linum usitatissimum L.). Front. Plant Sci. 2018, 9, 885. [Google Scholar] [CrossRef] [Green Version]
  75. You, F.M.; Xiao, J.; Li, P.; Yao, Z.; Jia, G.; He, L.; Zhu, T.; Luo, M.C.; Wang, X.; Deyholos, M.K.; et al. Chromosome-scale pseudomolecules refined by optical, physical and genetic maps in flax. Plant J. Cell Mol. Biol. 2018, 95, 371–384. [Google Scholar] [CrossRef]
  76. Zyablitsin, A.V.; Dmitriev, A.A.; Krasnov, G.S.; Bolsheva, N.L.; Rozhmina, T.A.; Muravenko, O.V.; Fedorova, M.S.; Snezhkina, A.V.; Kudryavtseva, A.V.; Melnikova, N.V. CAX3 Gene is Involved in Flax Response to High Soil Acidity and Aluminum Exposure. Mol. Biol. 2018, 52, 514–519. [Google Scholar] [CrossRef]
  77. Bolsheva, N.L.; Melnikova, N.V.; Kirov, I.V.; Dmitriev, A.A.; Krasnov, G.S.; Amosova, A.C.; Samatadze, T.E.; Yurkevich, O.Y.; Zoshchuk, S.A.; Kudryavtseva, A.V.; et al. Characterization of repeated DNA sequences in genomes of blue-flowered flax. BMC Evol. Biol. 2019, 19, 49. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Dmitriev, A.A.; Krasnov, G.S.; Rozhmina, T.A.; Zyablitsin, A.V.; Snezhkina, A.V.; Fedorova, M.S.; Pushkova, E.N.; Kezimana, P.; Novakovskiy, R.O.; Povkhova, L.V.; et al. Flax (Linum usitatissimum L.) response to non-optimal soil acidity and zinc deficiency. BMC Plant Biol. 2019, 19, 54. [Google Scholar] [CrossRef] [PubMed]
  79. Guo, D.; Jiang, H.; Yan, W.; Yang, L.; Ye, J.; Wang, Y.; Yan, Q.; Chen, J.; Gao, Y.; Duan, L.; et al. Resequencing 200 Flax Cultivated Accessions Identifies Candidate Genes Related to Seed Size and Weight and Reveals Signatures of Artificial Selection. Front. Plant Sci. 2019, 10, 1682. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Xie, D.; Dai, Z.; Yang, Z.; Tang, Q.; Deng, C.; Xu, Y.; Wang, J.; Chen, J.; Zhao, D.; Zhang, S.; et al. Combined genome-wide association analysis and transcriptome sequencing to identify candidate genes for flax seed fatty acid metabolism. Plant Sci. Int. J. Exp. Plant Biol. 2019, 286, 98–107. [Google Scholar] [CrossRef] [PubMed]
  81. Dmitriev, A.A.; Novakovskiy, R.O.; Pushkova, E.N.; Rozhmina, T.A.; Zhuchenko, A.A.; Bolsheva, N.L.; Beniaminov, A.D.; Mitkevich, V.A.; Povkhova, L.V.; Dvorianinova, E.M.; et al. Transcriptomes of Different Tissues of Flax (Linum usitatissimum L.) Cultivars With Diverse Characteristics. Front. Genet. 2020, 11, 565146. [Google Scholar] [CrossRef]
  82. Dmitriev, A.A.; Pushkova, E.N.; Novakovskiy, R.O.; Beniaminov, A.D.; Rozhmina, T.A.; Zhuchenko, A.A.; Bolsheva, N.L.; Muravenko, O.V.; Povkhova, L.V.; Dvorianinova, E.M.; et al. Genome Sequencing of Fiber Flax Cultivar Atlant Using Oxford Nanopore and Illumina Platforms. Front. Genet. 2020, 11, 590282. [Google Scholar] [CrossRef]
  83. Lan, S.; Zheng, C.; Hauck, K.; McCausland, M.; Duguid, S.D.; Booker, H.M.; Cloutier, S.; You, F.M. Genomic Prediction Accuracy of Seven Breeding Selection Traits Improved by QTL Identification in Flax. Int. J. Mol. Sci. 2020, 21, 1577. [Google Scholar] [CrossRef] [Green Version]
  84. Dvorianinova, E.M.; Pushkova, E.N.; Novakovskiy, R.O.; Povkhova, L.V.; Bolsheva, N.L.; Kudryavtseva, L.P.; Rozhmina, T.A.; Melnikova, N.V.; Dmitriev, A.A. Nanopore and Illumina Genome Sequencing of Fusarium oxysporum f. sp. lini Strains of Different Virulence. Front. Genet. 2021, 12, 662928. [Google Scholar] [CrossRef]
  85. Povkhova, L.V.; Melnikova, N.V.; Rozhmina, T.A.; Novakovskiy, R.O.; Pushkova, E.N.; Dvorianinova, E.M.; Zhuchenko, A.A.; Kamionskaya, A.M.; Krasnov, G.S.; Dmitriev, A.A. Genes Associated with the Flax Plant Type (Oil or Fiber) Identified Based on Genome and Transcriptome Sequencing Data. Plants 2021, 10, 2616. [Google Scholar] [CrossRef]
  86. Soto-Cerda, B.J.; Aravena, G.; Cloutier, S. Genetic dissection of flowering time in flax (Linum usitatissimum L.) through single- and multi-locus genome-wide association studies. Mol. Genet. Genom. MGG 2021, 296, 877–891. [Google Scholar] [CrossRef]
  87. Bolsheva, N.L.; Melnikova, N.V.; Dvorianinova, E.M.; Mironova, L.N.; Yurkevich, O.Y.; Amosova, A.V.; Krasnov, G.S.; Dmitriev, A.A.; Muravenko, O.V. Clarification of the Position of Linum stelleroides Planch. within the Phylogeny of the Genus Linum, L. Plants 2022, 11, 652. [Google Scholar] [CrossRef] [PubMed]
  88. House, M.A.; Young, L.W.; Robinson, S.J.; Booker, H.M. Transcriptomic Analysis of Early Flowering Signals in ’Royal’ Flax. Plants 2022, 11, 860. [Google Scholar] [CrossRef] [PubMed]
  89. Li, X.; Guo, D.; Xue, M.; Li, G.; Yan, Q.; Jiang, H.; Liu, H.; Chen, J.; Gao, Y.; Duan, L.; et al. Genome-Wide Association Study of Salt Tolerance at the Seed Germination Stage in Flax (Linum usitatissimum L.). Genes 2022, 13, 486. [Google Scholar] [CrossRef]
  90. You, F.M.; Rashid, K.Y.; Zheng, C.; Khan, N.; Li, P.; Xiao, J.; He, L.; Yao, Z.; Cloutier, S. Insights into the Genetic Architecture and Genomic Prediction of Powdery Mildew Resistance in Flax (Linum usitatissimum L.). Int. J. Mol. Sci. 2022, 23, 4960. [Google Scholar] [CrossRef] [PubMed]
  91. Young, L.; Akhov, L.; Kulkarni, M.; You, F.; Booker, H. Fine-mapping of a putative glutathione S-transferase (GST) gene responsible for yellow seed colour in flax (Linum usitatissimum). BMC Res. Notes 2022, 15, 72. [Google Scholar] [CrossRef] [PubMed]
  92. Dvorianinova, E.M.; Bolsheva, N.L.; Pushkova, E.N.; Rozhmina, T.A.; Zhuchenko, A.A.; Novakovskiy, R.O.; Povkhova, L.V.; Sigova, E.A.; Zhernova, D.A.; Borkhert, E.V.; et al. Isolating Linum usitatissimum L. Nuclear DNA Enabled Assembling High-Quality Genome. Int. J. Mol. Sci. 2022, 23, 3244. [Google Scholar] [CrossRef]
  93. Sun, C.; Dong, Z.; Zhao, L.; Ren, Y.; Zhang, N.; Chen, F. The Wheat 660K SNP array demonstrates great potential for marker-assisted selection in polyploid wheat. Plant Biotechnol. J. 2020, 18, 1354–1360. [Google Scholar] [CrossRef]
  94. You, Q.; Yang, X.; Peng, Z.; Xu, L.; Wang, J. Development and Applications of a High Throughput Genotyping Tool for Polyploid Crops: Single Nucleotide Polymorphism (SNP) Array. Front. Plant Sci. 2018, 9, 104. [Google Scholar] [CrossRef] [Green Version]
  95. Rasheed, A.; Hao, Y.; Xia, X.; Khan, A.; Xu, Y.; Varshney, R.K.; He, Z. Crop Breeding Chips and Genotyping Platforms: Progress, Challenges, and Perspectives. Mol. Plant 2017, 10, 1047–1064. [Google Scholar] [CrossRef] [Green Version]
  96. Dmitriev, A.; Pushkova, E.; Melnikova, N. Plant Genome Sequencing: Modern Technologies and Novel Opportunities for Breeding. Mol. Biol. 2022, 56, 495–507. [Google Scholar] [CrossRef]
  97. McCormack, J.E.; Hird, S.M.; Zellmer, A.J.; Carstens, B.C.; Brumfield, R.T. Applications of next-generation sequencing to phylogeography and phylogenetics. Mol. Phylogenet. Evol. 2013, 66, 526–538. [Google Scholar] [CrossRef] [PubMed]
  98. Cruaud, P.; Rasplus, J.Y.; Rodriguez, L.J.; Cruaud, A. High-throughput sequencing of multiple amplicons for barcoding and integrative taxonomy. Sci. Rep. 2017, 7, 41948. [Google Scholar] [CrossRef] [PubMed]
  99. Borkhert, E.; Krasnov, G.; Bolsheva, N.; Kezimana, P.; Yurkevich, O.; Muravenko, O.; Kudryavtseva, A.; Melnikova, N. Genetics polymorphism of poplars from Moscow region based on high-throughput sequencing of ITS. Vavilov J. Genet. Breed. 2018, 22, 531–535. [Google Scholar] [CrossRef] [Green Version]
  100. Melnikova, N.V.; Kudryavtseva, A.V.; Borkhert, E.V.; Pushkova, E.N.; Fedorova, M.S.; Snezhkina, A.V.; Krasnov, G.S.; Dmitriev, A.A. Sex-specific polymorphism of MET1 and ARR17 genes in Populus × sibirica. Biochimie 2019, 162, 26–32. [Google Scholar] [CrossRef]
  101. Novakovskiy, R.O.; Dvorianinova, E.M.; Rozhmina, T.A.; Kudryavtseva, L.P.; Gryzunov, A.A.; Pushkova, E.N.; Povkhova, L.V.; Snezhkina, A.V.; Krasnov, G.S.; Kudryavtseva, A.V.; et al. Data on genetic polymorphism of flax (Linum usitatissimum L.) pathogenic fungi of Fusarium, Colletotrichum, Aureobasidium, Septoria, and Melampsora genera. Data Brief 2020, 31, 105710. [Google Scholar] [CrossRef]
  102. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef]
  103. Li, H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv 2013, arXiv:1303.3997. [Google Scholar]
  104. Garrison, E.; Marth, G. Haplotype-based variant detection from short-read sequencing. arXiv 2012, arXiv:1207.3907. [Google Scholar]
Figure 1. HRM analysis for the G-to-A mutation in exon 1 of FAD3A gene for oligonucleotide standards. G—standard without the mutation; A—standard with the mutation; G/A—standards with and without the mutation in 1:1 ratio.
Figure 1. HRM analysis for the G-to-A mutation in exon 1 of FAD3A gene for oligonucleotide standards. G—standard without the mutation; A—standard with the mutation; G/A—standards with and without the mutation in 1:1 ratio.
Plants 12 00095 g001
Figure 2. HRM analysis for the C-to-T mutation in exon 2 of FAD3B gene for oligonucleotide standards. C—standard without the mutation; T—standard with the mutation; C/T—standards with and without the mutation in 1:1 ratio.
Figure 2. HRM analysis for the C-to-T mutation in exon 2 of FAD3B gene for oligonucleotide standards. C—standard without the mutation; T—standard with the mutation; C/T—standards with and without the mutation in 1:1 ratio.
Plants 12 00095 g002
Figure 3. HRM analysis for the C-to-T mutation in exon 5 of FAD3A gene for oligonucleotide standards. C—standard without the mutation; T—standard with the mutation; C/T—standards with and without the mutation in 1:1 ratio.
Figure 3. HRM analysis for the C-to-T mutation in exon 5 of FAD3A gene for oligonucleotide standards. C—standard without the mutation; T—standard with the mutation; C/T—standards with and without the mutation in 1:1 ratio.
Plants 12 00095 g003
Figure 4. HRM analysis for the G-to-A mutation in exon 1 of FAD3A gene for individual plants of line AGT 1535/07. G—samples without the mutation, G/A—samples heterozygous for the mutation; A—samples with the mutation.
Figure 4. HRM analysis for the G-to-A mutation in exon 1 of FAD3A gene for individual plants of line AGT 1535/07. G—samples without the mutation, G/A—samples heterozygous for the mutation; A—samples with the mutation.
Plants 12 00095 g004
Figure 5. HRM analysis for the C-to-T mutation in exon 2 of FAD3B gene for individual plants of line AGT 981/05. C—samples without the mutation, C/T—samples heterozygous for the mutation; T—samples with the mutation.
Figure 5. HRM analysis for the C-to-T mutation in exon 2 of FAD3B gene for individual plants of line AGT 981/05. C—samples without the mutation, C/T—samples heterozygous for the mutation; T—samples with the mutation.
Plants 12 00095 g005
Figure 6. HRM analysis for the C-to-T mutation in exon 5 of FAD3A gene for individual plants of variety Lola. C—samples without the mutation, C/T—samples heterozygous for the mutation; T—samples with the mutation.
Figure 6. HRM analysis for the C-to-T mutation in exon 5 of FAD3A gene for individual plants of variety Lola. C—samples without the mutation, C/T—samples heterozygous for the mutation; T—samples with the mutation.
Plants 12 00095 g006
Figure 7. Results of CAPS analysis for the G-to-A mutation in exon 1 of FAD3A for individual plants of line AGT 1535/07. If the mutation is absent, two fragments with the lengths of 195 and 107 bp were visualized (lines 5, 16, 19, 20), while if the mutation is present, one fragment with a length of 302 bp was visualized (lines 1–4, 6–8, 10–15, 17, 18, 21, 22). In the case of heterozygosity for the studied mutation, three DNA fragments were visualized, corresponding to the initial length of the amplicon and two of its parts after restriction (line 9). M—50 bp ladder.
Figure 7. Results of CAPS analysis for the G-to-A mutation in exon 1 of FAD3A for individual plants of line AGT 1535/07. If the mutation is absent, two fragments with the lengths of 195 and 107 bp were visualized (lines 5, 16, 19, 20), while if the mutation is present, one fragment with a length of 302 bp was visualized (lines 1–4, 6–8, 10–15, 17, 18, 21, 22). In the case of heterozygosity for the studied mutation, three DNA fragments were visualized, corresponding to the initial length of the amplicon and two of its parts after restriction (line 9). M—50 bp ladder.
Plants 12 00095 g007
Figure 8. Results of CAPS analysis for the C-to-T mutation in exon 2 of FAD3B for individual plants of line AGT 981/05. If the mutation is absent, two fragments with the lengths of 157 and 103 bp were visualized (lines 5, 9, 13), while if the mutation is present, one fragment with a length of 260 bp was visualized (lines 1–4, 6, 8, 10–12, 14–21). In the case of heterozygosity for the studied mutation, three DNA fragments were visualized, corresponding to the initial length of the amplicon and two of its parts after restriction (lines 7, 22). M—50 bp ladder.
Figure 8. Results of CAPS analysis for the C-to-T mutation in exon 2 of FAD3B for individual plants of line AGT 981/05. If the mutation is absent, two fragments with the lengths of 157 and 103 bp were visualized (lines 5, 9, 13), while if the mutation is present, one fragment with a length of 260 bp was visualized (lines 1–4, 6, 8, 10–12, 14–21). In the case of heterozygosity for the studied mutation, three DNA fragments were visualized, corresponding to the initial length of the amplicon and two of its parts after restriction (lines 7, 22). M—50 bp ladder.
Plants 12 00095 g008
Figure 9. Results of CAPS analysis for the C-to-T mutation in exon 5 of FAD3A for individual plants of variety Lola. If the mutation is absent, two fragments with the lengths of 163 and 109 bp were visualized (lines 2, 7, 8, 21), while if the mutation is present, one fragment with a length of 272 bp was visualized (1, 4–6, 9, 12–20, 22). In the case of heterozygosity for the studied mutation, three DNA fragments were visualized, corresponding to the initial length of the amplicon and two of its parts after restriction (3, 10, 11). M—50 bp ladder.
Figure 9. Results of CAPS analysis for the C-to-T mutation in exon 5 of FAD3A for individual plants of variety Lola. If the mutation is absent, two fragments with the lengths of 163 and 109 bp were visualized (lines 2, 7, 8, 21), while if the mutation is present, one fragment with a length of 272 bp was visualized (1, 4–6, 9, 12–20, 22). In the case of heterozygosity for the studied mutation, three DNA fragments were visualized, corresponding to the initial length of the amplicon and two of its parts after restriction (3, 10, 11). M—50 bp ladder.
Plants 12 00095 g009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Povkhova, L.V.; Pushkova, E.N.; Rozhmina, T.A.; Zhuchenko, A.A.; Frykin, R.I.; Novakovskiy, R.O.; Dvorianinova, E.M.; Gryzunov, A.A.; Borkhert, E.V.; Sigova, E.A.; et al. 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. Plants 2023, 12, 95. https://doi.org/10.3390/plants12010095

AMA Style

Povkhova LV, Pushkova EN, Rozhmina TA, Zhuchenko AA, Frykin RI, Novakovskiy RO, Dvorianinova EM, Gryzunov AA, Borkhert EV, Sigova EA, et al. 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. Plants. 2023; 12(1):95. https://doi.org/10.3390/plants12010095

Chicago/Turabian Style

Povkhova, Liubov V., Elena N. Pushkova, Tatiana A. Rozhmina, Alexander A. Zhuchenko, Roman I. Frykin, Roman O. Novakovskiy, Ekaterina M. Dvorianinova, Aleksey A. Gryzunov, Elena V. Borkhert, Elizaveta A. Sigova, and et al. 2023. "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" Plants 12, no. 1: 95. https://doi.org/10.3390/plants12010095

APA Style

Povkhova, L. V., Pushkova, E. N., Rozhmina, T. A., Zhuchenko, A. A., Frykin, R. I., Novakovskiy, R. O., Dvorianinova, E. M., Gryzunov, A. A., Borkhert, E. V., Sigova, E. A., Vladimirov, G. N., Snezhkina, A. V., Kudryavtseva, A. V., Krasnov, G. S., Dmitriev, A. A., & Melnikova, N. V. (2023). 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. Plants, 12(1), 95. https://doi.org/10.3390/plants12010095

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop