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

BACKGROUND: The 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. RESULTS: High-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). CONCLUSIONS: High-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.

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

BACKGROUND: Cultivated 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. RESULTS: High-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. CONCLUSIONS: Differentially 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.

Aerobiological Monitoring and Metabarcoding of Grass Pollen.

Grass pollen is one of the leading causes of pollinosis, affecting 10-30% of the world's population. The allergenicity of pollen from different Poaceae species is not the same and is estimated from moderate to high. Aerobiological monitoring is a standard method that allows one to track and predict the dynamics of allergen concentration in the air. Poaceae is a stenopalynous family, and thus grass pollen can usually be identified only at the family level with optical microscopy. Molecular methods, in particular the DNA barcoding technique, can be used to conduct a more accurate analysis of aerobiological samples containing the DNA of various plant species. This study aimed to test the possibility of using the ITS1 and ITS2 nuclear loci for determining the presence of grass pollen from air samples via metabarcoding and to compare the analysis results with the results of phenological observations. Based on the high-throughput sequencing data, we analyzed the changes in the composition of aerobiological samples taken in the Moscow and Ryazan regions for three years during the period of active flowering of grasses. Ten genera of the Poaceae family were detected in airborne pollen samples. The representation for most of them for ITS1 and ITS2 barcodes was similar. At the same time, in some samples, the presence of specific genera was characterized by only one sequence: either ITS1 or ITS2. Based on the analysis of the abundance of both barcode reads in the samples, the following order could describe the change with time in the dominant species in the air: Poa, Alopecurus, and Arrhenatherum in early mid-June, Lolium, Bromus, Dactylis, and Briza in mid-late June, Phleum, Elymus in late June to early July, and Calamagrostis in early mid-July. In most samples, the number of taxa found via metabarcoding analysis was higher compared to that in the phenological observations. The semi-quantitative analysis of high-throughput sequencing data well reflects the abundance of only major grass species at the flowering stage.

Complete plastome sequencing of Allium paradoxum reveals unusual rearrangements and the loss of the ndh genes as compared to Allium ursinum and other onions.

In this work the complete chloroplast DNAs of Allium paradoxum and Allium ursinum, two edible species of Allium subg. Amerallium (the first lineage), were sequenced, assembled, annotated, and compared with complete Allium plastomes of the second and third evolutionary lines from GenBank database. The A. ursinum plastome contains 90 predicted genes (81 unique) including 5 pseudogenes, while A. paradoxum has 88 predicted genes (79 unique) including 19 pseudogenes. The comparative analysis has revealed that the A. paradoxum plastome differs markedly from those of other species. Due to many deletions, the A. paradoxum plastome is the shortest of known for Allium species, being only 145,819 bp long. The most prominent distinctions are (1) a 4825 bp long local inversion that spans from the ndhE to the rpl32 gene in the small single copy region and (2) pseudogenization, or the loss of all NADH-genes. In contrast, the plastome of A. ursinum - a species from the first evolutionary line (as well as A. paradoxum) - resembles the Allium species of the second and third evolutionary lines, showing no large rearrangements or discrepancies in gene content. It is unclear yet whether only A. paradoxum was affected by some evolutionary events or its close relatives from both sect. Briseis and other sections of Amerallium were altered as well. We speculate the sunlight-intolerant, shade-loving nature of A. paradoxum and the impairment of the ndh genes in its plastome could be interrelated phenomena.

Improved Protocols of ITS1-Based Metabarcoding and Their Application in the Analysis of Plant-Containing Products.

Plants are widely used for food and beverage preparation, most often in the form of complex mixtures of dried and ground parts, such as teas, spices or herbal medicines. Quality control of such products is important due to the potential health risks from the presence of unlabelled components or absence of claimed ones. A promising approach to analyse such products is DNA metabarcoding due to its high resolution and sensitivity. However, this method's application in food analysis requires several methodology optimizations in DNA extraction, amplification and library preparation. In this study, we present such optimizations. The most important methodological outcomes are the following: 1) the DNA extraction method greatly influences amplification success; 2) the main problem for the application of metabarcoding is DNA purity, not integrity or quantity; and 3) the "non-amplifiable" samples can be amplified with polymerases resistant to inhibitors. Using this optimized workflow, we analysed a broad set of plant products (teas, spices and herbal remedies) using two NGS platforms. The analysis revealed the problem of both the presence of extraneous components and the absence of labelled ones. Notably, for teas, no correlation was found between the price and either the absence of labelled components or presence of unlabelled ones; for spices, a negative correlation was found between the price and presence of unlabelled components.

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.

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

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.

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.

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.

Impact of recombination on polymorphism of genes encoding Kunitz-type protease inhibitors in the genus Solanum.

BACKGROUND: The group of Kunitz-type protease inhibitors (KPI) from potato is encoded by a polymorphic family of multiple allelic and non-allelic genes. The previous explanations of the KPI variability were based on the hypothesis of random mutagenesis as a key factor of KPI polymorphism. RESULTS: KPI-A genes from the genomes of Solanum tuberosum cv. Istrinskii and the wild species Solanum palustre were amplified by PCR with subsequent cloning in plasmids. True KPI sequences were derived from comparison of the cloned copies. "Hot spots" of recombination in KPI genes were independently identified by DnaSP 4.0 and TOPALi v2.5 software. The KPI-A sequence from potato cv. Istrinskii was found to be 100% identical to the gene from Solanum nigrum. This fact illustrates a high degree of similarity of KPI genes in the genus Solanum. Pairwise comparison of KPI A and B genes unambiguously showed a non-uniform extent of polymorphism at different nt positions. Moreover, the occurrence of substitutions was not random along the strand. Taken together, these facts contradict the traditional hypothesis of random mutagenesis as a principal source of KPI gene polymorphism. The experimentally found mosaic structure of KPI genes in both plants studied is consistent with the hypothesis suggesting recombination of ancestral genes. The same mechanism was proposed earlier for other resistance-conferring genes in the nightshade family (Solanaceae). Based on the data obtained, we searched for potential motifs of site-specific binding with plant DNA recombinases. During this work, we analyzed the sequencing data reported by the Potato Genome Sequencing Consortium (PGSC), 2011 and found considerable inconsistence of their data concerning the number, location, and orientation of KPI genes of groups A and B. CONCLUSIONS: The key role of recombination rather than random point mutagenesis in KPI polymorphism was demonstrated for the first time.