Cultivated flax (Linum usitatissimum L.) is an important plant valuable for industry. Some flax lines can undergo heritable phenotypic and genotypic changes (LIS-1 insertion being the most common) in response to nutrient stress and are called plastic lines. Offspring of plastic lines, which stably inherit the changes, are called genotrophs. MicroRNAs (miRNAs) are involved in a crucial regulatory mechanism of gene expression. They have previously been assumed to take part in nutrient stress response and can, therefore, participate in genotroph formation. In the present study, we performed high-throughput sequencing of small RNAs (sRNAs) extracted from flax plants grown under normal, phosphate deficient and nutrient excess conditions to identify miRNAs and evaluate their expression. Our analysis revealed expression of 96 conserved miRNAs from 21 families in flax. Moreover, 475 novel potential miRNAs were identified for the first time, and their targets were predicted. However, none of the identified miRNAs were transcribed from LIS-1. Expression of seven miRNAs (miR168, miR169, miR395, miR398, miR399, miR408, and lus-miR-N1) with up- or down-regulation under nutrient stress (on the basis of high-throughput sequencing data) was evaluated on extended sampling using qPCR. Reference gene search identified ETIF3H and ETIF3E genes as most suitable for this purpose. Down-regulation of novel potential lus-miR-N1 and up-regulation of conserved miR399 were revealed under the phosphate deficient conditions. In addition, the negative correlation of expression of lus-miR-N1 and its predicted target, ubiquitin-activating enzyme E1 gene, as well as, miR399 and its predicted target, ubiquitin-conjugating enzyme E2 gene, was observed. Thus, in our study, miRNAs expressed in flax plastic lines and genotrophs were identified and their expression and expression of their targets was evaluated using high-throughput sequencing and qPCR for the first time. These data provide new insights into nutrient stress response regulation in plastic flax cultivars.

BackgroundFlax (Linum usitatissimum L.) is a crop plant used for fiber and oil production. Although potentially high-yielding flax varieties have been developed, environmental stresses markedly decrease flax production. Among biotic stresses, Fusarium oxysporum f. sp. lini is recognized as one of the most devastating flax pathogens. It causes wilt disease that is one of the major limiting factors for flax production worldwide. Breeding and cultivation of flax varieties resistant to F. oxysporum is the most effective method for controlling wilt disease. Although the mechanisms of flax response to Fusarium have been actively studied, data on the plant response to infection and resistance gene candidates are currently very limited.ResultsThe transcriptomes of two resistant and two susceptible flax cultivars with respect to Fusarium wilt, as well as two resistant BC2F5 populations, which were grown under control conditions or inoculated with F. oxysporum, were sequenced using the Illumina platform. Genes showing changes in expression under F. oxysporum infection were identified in both resistant and susceptible flax genotypes. We observed the predominant overexpression of numerous genes that are involved in defense response. This was more pronounced in resistant cultivars. In susceptible cultivars, significant downregulation of genes involved in cell wall organization or biogenesis was observed in response to F. oxysporum. In the resistant genotypes, upregulation of genes related to NAD(P)H oxidase activity was detected. Upregulation of a number of genes, including that encoding beta-1,3-glucanase, was significantly greater in the cultivars and BC2F5 populations resistant to Fusarium wilt than in susceptible cultivars in response to F. oxysporum infection.ConclusionsUsing high-throughput sequencing, we identified genes involved in the early defense response of L. usitatissimum against the fungus F. oxysporum. In response to F. oxysporum infection, we detected changes in the expression of pathogenesis-related protein-encoding genes and genes involved in ROS production or related to cell wall biogenesis. Furthermore, we identified genes that were upregulated specifically in flax genotypes resistant to Fusarium wilt. We suggest that the identified genes in resistant cultivars and BC2F5 populations showing induced expression in response to F. oxysporum infection are the most promising resistance gene candidates.Electronic supplementary materialThe online version of this article (10.1186/s12870-017-1192-2) contains supplementary material, which is available to authorized users.

Background Flax (Linum usitatissimum L.) is grown for fiber and seed in many countries. Flax cultivars differ in the oil composition and, depending on the ratio of fatty acids, are used in pharmaceutical, food, or paint industries. It is known that genes of SAD (stearoyl-ACP desaturase) and FAD (fatty acid desaturase) families play a key role in the synthesis of fatty acids, and some alleles of these genes are associated with a certain composition of flax oil. However, data on genetic polymorphism of these genes are still insufficient. Results On the basis of the collection of the Institute for Flax (Torzhok, Russia), we formed a representative set of 84 cultivars and lines reflecting the diversity of fatty acid composition of flax oil. An approach for the determination of full-length sequences of SAD1, SAD2, FAD2A, FAD2B, FAD3A, and FAD3B genes using the Illumina platform was developed and deep sequencing of the 6 genes in 84 flax samples was performed on MiSeq. The obtained high coverage (about 400x on average) enabled accurate assessment of polymorphisms in SAD1, SAD2, FAD2A, FAD2B, FAD3A, and FAD3B genes and evaluation of cultivar/line heterogeneity. The highest level of genetic diversity was observed for FAD3A and FAD3B genes – 91 and 62 polymorphisms respectively. Correlation analysis revealed associations between particular variants in SAD and FAD genes and predominantly those fatty acids whose conversion they catalyze: SAD – stearic and oleic acids, FAD2 – oleic and linoleic acids, FAD3 – linoleic and linolenic acids. All except one low-linolenic flax cultivars/lines contained both the substitution of tryptophan to stop codon in the FAD3A gene and histidine to tyrosine substitution in the FAD3B gene, while samples with only one of these polymorphisms had medium content of linolenic acid and cultivars/lines without them were high-linolenic. Conclusions Genetic polymorphism of SAD and FAD genes was evaluated in the collection of flax cultivars and lines with diverse oil composition, and associations between particular polymorphisms and the ratio of fatty acids were revealed. The achieved results are the basis for the development of marker-assisted selection and DNA-based certification of flax cultivars.

About 30% of the world's ice-free land area is occupied by acid soils. In soils with pH below 5, aluminum (Al) releases to the soil solution, and becomes highly toxic for plants. Therefore, breeding of varieties that are resistant to Al is needed. Flax (Linum usitatissimum L.) is grown worldwide for fiber and seed production. Al toxicity in acid soils is a serious problem for flax cultivation. However, very little is known about mechanisms of flax resistance to Al and the genetics of this resistance. In the present work, we sequenced 16 transcriptomes of flax cultivars resistant (Hermes and TMP1919) and sensitive (Lira and Orshanskiy) to Al, which were exposed to control conditions and aluminum treatment for 4, 12, and 24 h. In total, 44.9–63.3 million paired-end 100-nucleotide reads were generated for each sequencing library. Based on the obtained high-throughput sequencing data, genes with differential expression under aluminum exposure were revealed in flax. The majority of the top 50 up-regulated genes were involved in transmembrane transport and transporter activity in both the Al-resistant and Al-sensitive cultivars. However, genes encoding proteins with glutathione transferase and UDP-glycosyltransferase activity were in the top 50 up-regulated genes only in the flax cultivars resistant to aluminum. For qPCR analysis in extended sampling, two UDP-glycosyltransferases (UGTs), and three glutathione S-transferases (GSTs) were selected. The general trend of alterations in the expression of the examined genes was the up-regulation under Al stress, especially after 4 h of Al exposure. Moreover, in the flax cultivars resistant to aluminum, the increase in expression was more pronounced than that in the sensitive cultivars. We speculate that the defense against the Al toxicity via GST antioxidant activity is the probable mechanism of the response of flax plants to aluminum stress. We also suggest that UGTs could be involved in cell wall modification and protection from reactive oxygen species (ROS) in response to Al stress in L. usitatissimum. Thus, GSTs and UGTs, probably, play an important role in the response of flax to Al via detoxification of ROS and cell wall modification.