Maxim S. Belenikin

Identification, Expression Analysis, and Target Prediction of Flax Genotroph MicroRNAs Under Normal and Nutrient Stress Conditions.

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.

Gene expression profiling of flax ( Linum usitatissimum L.) under edaphic stress

BackgroundCultivated flax (Linum usitatissimum L.) is widely used for production of textile, food, chemical and pharmaceutical products. However, various stresses decrease flax production. Search for genes, which are involved in stress response, is necessary for breeding of adaptive cultivars. Imbalanced concentration of nutrient elements in soil decrease flax yields and also results in heritable changes in some flax lines. The appearance of Linum Insertion Sequence 1 (LIS-1) is the most studied modification. However, LIS-1 function is still unclear.ResultsHigh-throughput sequencing of transcriptome of flax plants grown under normal (N), phosphate deficient (P), and nutrient excess (NPK) conditions was carried out using Illumina platform. The assembly of transcriptome was performed, and a total of 34924, 33797, and 33698 unique transcripts for N, P, and NPK sequencing libraries were identified, respectively. We have not revealed any LIS-1 derived mRNA in our sequencing data. The analysis of high-throughput sequencing data allowed us to identify genes with potentially differential expression under imbalanced nutrition. For further investigation with qPCR, 15 genes were chosen and their expression levels were evaluated in the extended sampling of 31 flax plants. Significant expression alterations were revealed for genes encoding WRKY and JAZ protein families under P and NPK conditions. Moreover, the alterations of WRKY family genes differed depending on LIS-1 presence in flax plant genome. Besides, we revealed slight and LIS-1 independent mRNA level changes of KRP2 and ING1 genes, which are adjacent to LIS-1, under nutrition stress.ConclusionsDifferentially expressed genes were identified in flax plants, which were grown under phosphate deficiency and excess nutrition, on the basis of high-throughput sequencing and qPCR data. We showed that WRKY and JAS gene families participate in flax response to imbalanced nutrient content in soil. Besides, we have not identified any mRNA, which could be derived from LIS-1, in our transcriptome sequencing data. Expression of LIS-1 flanking genes, ING1 and KRP2, was suggested not to be nutrient stress-induced. Obtained results provide new insights into edaphic stress response in flax and the role of LIS-1 in these process.

Retrotransposon-based molecular markers for analysis of genetic diversity within the Genus Linum.

SSAP method was used to study the genetic diversity of 22 Linum species from sections Linum, Adenolinum, Dasylinum, Stellerolinum, and 46 flax cultivars. All the studied flax varieties were distinguished using SSAP for retrotransposons FL9 and FL11. Thus, the validity of SSAP method was demonstrated for flax marking, identification of accessions in genebank collections, and control during propagation of flax varieties. Polymorphism of Fl1a, Fl1b, and Cassandra insertions were very low in flax varieties, but these retrotransposons were successfully used for the investigation of Linum species. Species clusterization based on SSAP markers was in concordance with their taxonomic division into sections Dasylinum, Stellerolinum, Adenolinum, and Linum. All species of sect. Adenolinum clustered apart from species of sect. Linum. The data confirmed the accuracy of the separation in these sections. Members of section Linum are not as closely related as members of other sections, so taxonomic revision of this section is desirable. L. usitatissimum accessions genetically distant from modern flax cultivars were revealed in our work. These accessions are of utmost interest for flax breeding and introduction of new useful traits into flax cultivars. The chromosome localization of Cassandra retrotransposon in Linum species was determined.

Excess fertilizer responsive miRNAs revealed in Linum usitatissimum L.

Effective fertilizer application is necessary to increase crop yields and reduce risk of plant overdosing. It is known that expression level of microRNAs (miRNAs) alters in plants under different nutrient concentrations in soil. The aim of our study was to identify and characterize miRNAs with expression alterations under excessive fertilizer in agriculturally important crop - flax (Linum usitatissimum L.). We have sequenced small RNAs in flax grown under normal and excessive fertilizer using Illumina GAIIx. Over 14 million raw reads was obtained for two small RNA libraries. 84 conserved miRNAs from 20 families were identified. Differential expression was revealed for several flax miRNAs under excessive fertilizer according to high-throughput sequencing data. For 6 miRNA families (miR395, miR169, miR408, miR399, miR398 and miR168) expression level alterations were evaluated on the extended sampling using qPCR. Statistically significant up-regulation was revealed for miR395 under excessive fertilizer. It is known that target genes of miR395 are involved in sulfate uptake and assimilation. However, according to our data alterations of the expression level of miR395 could be associated not only with excess sulfur application, but also with redundancy of other macro- and micronutrients. Furthermore expression level was evaluated for miRNAs and their predicted targets. The negative correlation between miR399 expression and expression of its predicted target ubiquitin-conjugating enzyme E2 gene was shown in flax for the first time. So we suggested miR399 involvement in phosphate regulation in L. usitatissimum. Revealed in our study expression alterations contribute to miRNA role in flax response to excessive fertilizer.

miR319, miR390, and miR393 Are Involved in Aluminum Response in Flax (Linum usitatissimum L.)

Acid soils limit agricultural production worldwide. Major reason of crop losses in acid soils is the toxicity of aluminum (Al). In the present work, we investigated expression alterations of microRNAs in flax (Linum usitatissimum L.) plants under Al stress. Flax seedlings of resistant (TMP1919 and G1071/4_k) and sensitive (Lira and G1071/4_o) to Al cultivars and lines were exposed to AlCl3 solution for 4 and 24 hours. Twelve small RNA libraries were constructed and sequenced using Illumina platform. In total, 97 microRNAs from 18 conserved families were identified. miR319, miR390, and miR393 revealed expression alterations associated with Al treatment of flax plants. Moreover, for miR390 and miR393, the alterations were distinct in sensitive and resistant to Al genotypes. Expression level changes of miR319 and miR390 were confirmed using qPCR analysis. In flax, potential targets of miR319 are TCPs, miR390–TAS3 and GRF5, and miR393–AFB2-coding transcripts. TCPs, TAS3, GRF5, and AFB2 participate in regulation of plant growth and development. The involvement of miR319, miR390, and miR393 in response to Al stress in flax was shown here for the first time. We speculate that these microRNAs play an important role in Al response via regulation of growth processes in flax plants.

Comparative analysis of inverted repeats of polypod fern (Polypodiales) plastomes reveals two hypervariable regions

BackgroundFerns are large and underexplored group of vascular plants (~ 11 thousands species). The genomic data available by now include low coverage nuclear genomes sequences and partial sequences of mitochondrial genomes for six species and several plastid genomes.ResultsWe characterized plastid genomes of three species of Dryopteris, which is one of the largest fern genera, using sequencing of chloroplast DNA enriched samples and performed comparative analysis with available plastomes of Polypodiales, the most species-rich group of ferns. We also sequenced the plastome of Adianthum hispidulum (Pteridaceae). Unexpectedly, we found high variability in the IR region, including duplication of rrn16 in D. blanfordii, complete loss of trnI-GAU in D. filix-mas, its pseudogenization due to the loss of an exon in D. blanfordii. Analysis of previously reported plastomes of Polypodiales demonstrated that Woodwardia unigemmata and Lepisorus clathratus have unusual insertions in the IR region. The sequence of these inserted regions has high similarity to several LSC fragments of ferns outside of Polypodiales and to spacer between tRNA-CGA and tRNA-TTT genes of mitochondrial genome of Asplenium nidus. We suggest that this reflects the ancient DNA transfer from mitochondrial to plastid genome occurred in a common ancestor of ferns. We determined the marked conservation of gene content and relative evolution rate of genes and intergenic spacers in the IRs of Polypodiales. Faster evolution of the four intergenic regions had been demonstrated (trnA- orf42, rrn16-rps12, rps7-psbA and ycf2-trnN).ConclusionsIRs of Polypodiales plastomes are dynamic, driven by such events as gene loss, duplication and putative lateral transfer from mitochondria.

The systematic position of Dryopteris blanfordii subsp. nigrosquamosa (Ching) Fraser-Jenkins within the genus Dryopteris Adans.

Abstract Dryopteris blanfordii (C.Hope) C.Chr. is a member of the Dryopteridaceae, growing in high altitude Picea or Abies forests (2900–3500 m) in China and India. Phylogenetic relationships between D. blanfordii subsp. nigrosquamosa and closely related species of Dryopteris were investigated using a combined analysis of multiple molecular data sets (the protein-coding region of rbcL and matK genes and intergenic spacers psbA-trnH, trnP-petG, rps4-trnS, trnL-trnF and rbcL-accD). An assumption about the position of D. blanfordii subsp. nigrosquamosa within Dryopteris was made by using the Maximum Likelihood and Bayesian Inference approach and chloroplast marker sequences of Dryopteris species from GenBank. The results demonstrated that Asian taxa D. blanfordii subsp. nigrosquamosa and D. laeta as well as two American species D. arguta and D. marginalis belong to the same clade, all four of them being part of Dryopteris section Dryopteris.

Erratum to “miR319, miR390, and miR393 Are Involved in Aluminum Response in Flax (Linum usitatissimum L.)”

[This corrects the article DOI: 10.1155/2017/4975146.].

Evolution of blue-flowered species of genus Linum based on high-throughput sequencing of ribosomal RNA genes

BackgroundThe species relationships within the genus Linum have already been studied several times by means of different molecular and phylogenetic approaches. Nevertheless, a number of ambiguities in phylogeny of Linum still remain unresolved. In particular, the species relationships within the sections Stellerolinum and Dasylinum need further clarification. Also, the question of independence of the species of the section Adenolinum still remains unanswered. Moreover, the relationships of L. narbonense and other species of the section Linum require further clarification. Additionally, the origin of tetraploid species of the section Linum (2n = 30) including the cultivated species L. usitatissimum has not been explored. The present study examines the phylogeny of blue-flowered species of Linum by comparisons of 5S rRNA gene sequences as well as ITS1 and ITS2 sequences of 35S rRNA genes.ResultsHigh-throughput sequencing has been used for analysis of multicopy rRNA gene families. In addition to the molecular phylogenetic analysis, the number and chromosomal localization of 5S and 35S rDNA sites has been determined by FISH.Our findings confirm that L. stelleroides forms a basal branch from the clade of blue-flowered flaxes which is independent of the branch formed by species of the sect. Dasylinum. The current molecular phylogenetic approaches, the cytogenetic analysis as well as different genomic DNA fingerprinting methods applied previously did not discriminate certain species within the sect. Adenolinum. The allotetraploid cultivated species L. usitatissimum and its wild ancestor L. angustifolium (2n = 30) could originate either as the result of hybridization of two diploid species (2n = 16) related to the modern L. gandiflorum and L. decumbens, or hybridization of a diploid species (2n = 16) and a diploid ancestor of modern L. narbonense (2n = 14).ConclusionsHigh-throughput sequencing of multicopy rRNA gene families allowed us to make several adjustments to the phylogeny of blue-flowered flax species and also reveal intra- and interspecific divergence of the rRNA gene sequences.