Nadezhda L. Bolsheva

The unique genome of two-chromosome grasses Zingeria and Colpodium, its origin, and evolution

Chromosome C-banding and two-color fluorescent in situ hybridization (FISH) were used to compare the chromosomes, to identify the chromosomal localization of the 45S and 5S rRNA genes, and to analyze the sequences of internal transcribed spacers 1 and 2 (ITS1 and ITS2) of the 45S rRNA genes in the genomes of grasses Zingeria biebersteiniana (2n = 4), Z. pisidica, Z. trichopoda (2n = 8), Colpodium versicolor (2n = 4), and Catabrosella variegata (syn. Colpodium variegatum) (2 n = 10). Differences in C-banding pattern were observed for two Z. biebersteiniana accessions from different localities. Similar C-banding patterns of chromosomes 1 and 2 were demonstrated for the Z. pisidica and Z. biebersteininana karyotypes. Chromosome C banding and localization of the 45S and 5S rRNA genes on the chromosomes of the two Zingeria species confirmed the assumption that Z. pisidica is an allotetraploid with one of the subgenomes similar to the Z. biebersteiniana genome. ITS comparisons showed that the unique two-chromosome grasses (x = 2)—Z. biebersteiniana (2n = 4), Z. trichopoda (2n = 8), Z. pisidica (2n = 8), and C. versicolor (2n = 4), which were earlier assigned to different tribes of subtribes of the family Poaceae—represent two closely related genera, the genetic distance (p-distance) between their ITSs being only 1.2–4.4%. The Zingeria species and C. versicolor formed a common clade with Catabrosella araratica (2n = 42, x = 7) on a molecular phylogenetic tree. Thus, the karyotypes of Zingeria and Colpodium, which have the lowest known basic chromosome number (x = 2), proved to be monophyletic, rather than originating from different phylogenetic lineages.

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.

Karyogenomics of species of the genus Linum L.

Using the molecular cytogenetic and RAPD methods of analysis, we studied genomes of 22 cultivated flax varieties and 24 wild species from six sections of the genus Linum L. The chromosome numbers were exactly determined in the karyotypes of all studied species, and all individual chromosomes were identified by the C/DAPI-banding pattern and localization of 26S rDNA and 5S rDNA. B chromosomes were discovered and studied for the first time in species of the section Syllinum Griseb. According to the data obtained, the species studied were divided into eight groups on the basis of similarity of their karyotypes, which corresponded in general to their clustering based on the RAPD results. The systematic positions and phylogenetic relationships of the flax species were verified.

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.

Glutathione S-transferases and UDP-glycosyltransferases Are Involved in Response to Aluminum Stress in Flax

About 30% of the worlds 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.

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.

The diversity of karyotypes and genomes within section Syllinum of the Genus Linum (Linaceae) revealed by molecular cytogenetic markers and RAPD analysis.

The wide variation in chromosome number found in species of the genus Linum (2n = 16, 18, 20, 26, 28, 30, 32, 36, 42, 72, 84) indicates that chromosomal mutations have played an important role in the speciation of this taxon. To contribute to a better understanding of the genetic diversity and species relationships in this genus, comparative studies of karyotypes and genomes of species within section Syllinum Griseb. (2n = 26, 28) were carried out. Elongated with 9-aminoacridine chromosomes of 10 species of section Syllinum were investigated by C- and DAPI/С-banding, CMA and Ag-NOR-staining, FISH with probes of rDNA and of telomere repeats. RAPD analysis was also performed. All the chromosome pairs in karyotypes of the studied species were identified. Chromosome DAPI/C-banding patterns of 28-chromosomal species were highly similar. Two of the species differed from the others in chromosomal location of rDNA sites. B chromosomes were revealed in all the 28-chromosomal species. Chromosomes of Linum nodiflorum L. (2n = 26) and the 28-chromosomal species were similar in DAPI/C-banding pattern and localization of several rDNA sites, but they differed in chromosomal size and number. The karyotype of L. nodiflorum was characterized by an intercalary site of telomere repeat, one additional 26S rDNA site and also by the absence of B chromosomes. Structural similarities between different chromosome pairs in karyotypes of the studied species were found indicating their tetraploid origin. RAPD analysis did not distinguish the species except L. nodiflorum. The species of section Syllinum probably originated from a common tetraploid ancestor. The 28-chromosomal species were closely related, but L. nodiflorum diverged significantly from the rest of the species probably due to chromosomal rearrangements occurring during evolution.

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.

The role of microRNA in abiotic stress response in plants

Regulation of gene expression via microRNA is the key mechanism of response to biotic and abiotic stresses in plants. There are a lot of experimental data on the biological function of microRNAs in response to different stresses in various plant species. This review contains up-to-date information on molecular mechanisms of microRNA action in plants in response to abiotic stresses, including drought, salinity, mineral nutrient deficiency or imbalance.

Molecular Cytogenetic Analysis of Deschampsia antarctica Desv. (Poaceae), Maritime Antarctic

Deschampsia antarctica Desv. (Poaceae) (2n = 26) is one of the two vascular plants adapted to the harshest environment of the Antarctic. Although the species is a valuable model for study of environmental stress tolerance in plants, its karyotype is still poorly investigated. We firstly conducted a comprehensive molecular cytogenetic analysis of D. antarctica collected on four islands of the Maritime Antarctic. D. antarctica karyotypes were studied by Giemsa C- and DAPI/C-banding, Ag-NOR staining, multicolour fluorescence in situ hybridization with repeated DNA probes (pTa71, pTa794, telomere repeats, pSc119.2, pAs1) and the GAA simple sequence repeat probe. We also performed sequential rapid in situ hybridization with genomic DNA of D. caespitosa. Two chromosome pairs bearing transcriptionally active 45S rDNA loci and five pairs with 5S rDNA sites were detected. A weak intercalary site of telomere repeats was revealed on the largest chromosome in addition to telomere hybridization signals at terminal positions. This fact confirms indirectly the hypothesis that chromosome fusion might have been the cause of the unusual for cereals chromosome number in this species. Based on patterns of distribution of the examined molecular cytogenetic markers, all chromosomes in karyotypes were identified, and chromosome idiograms of D. antarctica were constructed. B chromosomes were found in most karyotypes of plants from Darboux Island. A mixoploid plant with mainly triploid cells bearing a Robertsonian rearrangement was detected among typical diploid specimens from Great Jalour Island. The karyotype variability found in D. antarctica is probably an expression of genome instability induced by environmental stress factors. The differences in C-banding patterns and in chromosome distribution of rDNA loci as well as homologous highly repeated DNA sequences detected between genomes of D. antarctica and its related species D. caespitosa indicate that genome reorganization involving coding and noncoding repeated DNA sequences had occurred during the divergence of these species.