Expression of FAD and SAD Genes in Developing Seeds of Flax Varieties under Different Growth Conditions

Flax seed is one of the richest plant sources of linolenic acid (LIN) and also contains unsaturated linoleic acid (LIO) and oleic acid (OLE). Stearoyl-ACP desaturases (SADs) and fatty acid desaturases (FADs) play key roles in the synthesis of flax fatty acids (FAs). However, there is no holistic view of which genes from the SAD and FAD families and at which developmental stages have the highest expression levels in flax seeds, as well as the influence of genotype and growth conditions on the expression profiles of these genes. We sequenced flax seed transcriptomes at 3, 7, 14, 21, and 28 days after flowering (DAF) for ten flax varieties with different oil FA compositions grown under three temperature/watering conditions. The expression levels of 25 genes of the SAD, FAD2, and FAD3 families were evaluated. FAD3b, FAD3a, FAD2b-2, SAD3-1, SAD2-1, SAD2-2, SAD3-2, FAD2a-1, and FAD2a-2 had the highest expression levels, which changed significantly during seed development. These genes probably play a key role in FA synthesis in flax seeds. High temperature and insufficient watering shifted the maximum expression levels of FAD and SAD genes to earlier developmental stages, while the opposite trend was observed for low temperature and excessive watering. Differences in the FAD and SAD expression profiles under different growth conditions may affect the FA composition of linseed oil. Stop codons in the FAD3a gene, resulting in a reduced LIN content, decreased the level of FAD3a transcript. The obtained results provide new insights into the synthesis of linseed oil.

Editing Metabolism, Sex, and Microbiome: How Can We Help Poplar Resist Pathogens?

Poplar (Populus) is a genus of woody plants of great economic value. Due to the growing economic importance of poplar, there is a need to ensure its stable growth by increasing its resistance to pathogens. Genetic engineering can create organisms with improved traits faster than traditional methods, and with the development of CRISPR/Cas-based genome editing systems, scientists have a new highly effective tool for creating valuable genotypes. In this review, we summarize the latest research data on poplar diseases, the biology of their pathogens and how these plants resist pathogens. In the final section, we propose to plant male or mixed poplar populations; consider the genes of the MLO group, transcription factors of the WRKY and MYB families and defensive proteins BbChit1, LJAMP2, MsrA2 and PtDef as the most promising targets for genetic engineering; and also pay attention to the possibility of microbiome engineering.

Culturable Bacterial Endophytes of Wild White Poplar (Populus alba L.) Roots: A First Insight into Their Plant Growth-Stimulating and Bioaugmentation Potential

The white poplar (Populus alba L.) has good potential for a green economy and phytoremediation. Bioaugmentation using endophytic bacteria can be considered as a safe strategy to increase poplar productivity and its resistance to toxic urban conditions. The aim of our work was to find the most promising strains of bacterial endophytes to enhance the growth of white poplar in unfavorable environmental conditions. To this end, for the first time, we performed whole-genome sequencing of 14 bacterial strains isolated from the tissues of the roots of white poplar in different geographical locations. We then performed a bioinformatics search to identify genes that may be useful for poplar growth and resistance to environmental pollutants and pathogens. Almost all endophytic bacteria obtained from white poplar roots are new strains of known species belonging to the genera Bacillus, Corynebacterium, Kocuria, Micrococcus, Peribacillus, Pseudomonas, and Staphylococcus. The genomes of the strains contain genes involved in the enhanced metabolism of nitrogen, phosphorus, and metals, the synthesis of valuable secondary metabolites, and the detoxification of heavy metals and organic pollutants. All the strains are able to grow on media without nitrogen sources, which indicates their ability to fix atmospheric nitrogen. It is concluded that the strains belonging to the genus Pseudomonas and bacteria of the species Kocuria rosea have the best poplar growth-stimulating and bioaugmentation potential, and the roots of white poplar are a valuable source for isolation of endophytic bacteria for possible application in ecobiotechnology.

Selection of Flax Genotypes for Pan-Genomic Studies by Sequencing Tagmentation-Based Transcriptome Libraries

Flax (Linum usitatissimum L.) products are used in the food, pharmaceutical, textile, polymer, medical, and other industries. The creation of a pan-genome will be an important advance in flax research and breeding. The selection of flax genotypes that sufficiently cover the species diversity is a crucial step for the pan-genomic study. For this purpose, we have adapted a method based on Illumina sequencing of transcriptome libraries prepared using the Tn5 transposase (tagmentase). This approach reduces the cost of sample preparation compared to commercial kits and allows the generation of a large number of cDNA libraries in a short time. RNA-seq data were obtained for 192 flax plants (3-6 individual plants from 44 flax accessions of different morphology and geographical origin). Evaluation of the genetic relationship between flax plants based on the sequencing data revealed incorrect species identification for five accessions. Therefore, these accessions were excluded from the sample set for the pan-genomic study. For the remaining samples, typical genotypes were selected to provide the most comprehensive genetic diversity of flax for pan-genome construction. Thus, high-throughput sequencing of tagmentation-based transcriptome libraries showed high efficiency in assessing the genetic relationship of flax samples and allowed us to select genotypes for the flax pan-genomic analysis.

Sex-determining region complements traditionally used in phylogenetic studies nuclear and chloroplast sequences in investigation of Aigeiros Duby and Tacamahaca Spach poplars (genus Populus L., Salicaceae)

Members of the genus Populus L. play an important role in the formation of forests in the northern hemisphere and are used in urban landscaping and timber production. Populus species of closely related sections show extensive hybridization. Therefore, the systematics of the genus is rather complicated, especially for poplars of hybrid origin. We aimed to assess the efficiency of application of the sex-determining region (SDR) in addition to the nuclear and chloroplast genome loci traditionally used in phylogenetic studies of poplars to investigate relationships in sections Aigeiros Duby and Tacamahaca Spach. Targeted deep sequencing of NTS 5S rDNA, ITS, DSH 2, DSH 5, DSH 8, DSH 12, DSH 29, 6, 15, 16, X18, trnG-psbK-psbI, rps2-rpoC2, rpoC2-rpoC1, as well as SDR and ARR17 gene was performed for 379 poplars. The SDR and ARR17 gene together with traditionally used multicopy and single-copy loci of nuclear and chloroplast DNA allowed us to obtain a clustering that is most consistent with poplar systematics based on morphological data and to shed light on several controversial hypotheses about the origin of the studied taxa (for example, the inexpediency of separating P. koreana, P. maximowiczii, and P. suaveolens into different species). We present a scheme of relationships between species and hybrids of sections Aigeiros and Tacamahaca based on molecular genetic, morphological, and geographical data. The geographical proximity of species and, therefore, the possibility of hybridization between them appear to be more important than the affiliation of species to the same section. We speculate that sections Aigeiros and Tacamahaca are distinguished primarily on an ecological principle (plain and mountain poplars) rather than on a genetic basis. Joint analysis of sequencing data for the SDR and chloroplast genome loci allowed us to determine the ancestors of P. × petrovskoe - P. laurifolia (female tree) × P. × canadensis (male tree), and P. × rasumovskoe - P. nigra (female tree) × P. suaveolens (male tree). Thus, the efficiency of using the SDR for the study of poplars of sections Aigeiros and Tacamahaca and the prospects of its use for the investigation of species of the genus Populus were shown.

Isolating Linum usitatissimum L. Nuclear DNA Enabled Assembling High-Quality Genome

High-quality genome sequences help to elucidate the genetic basis of numerous biological processes and track species evolution. For flax (Linum usitatissimum L.)—a multifunctional crop, high-quality assemblies from Oxford Nanopore Technologies (ONT) data were unavailable, largely due to the difficulty of isolating pure high-molecular-weight DNA. This article proposes a scheme for gaining a contiguous L. usitatissimum assembly using Nanopore data. We developed a protocol for flax nuclei isolation with subsequent DNA extraction, which allows obtaining about 5 μg of pure high-molecular-weight DNA from 0.5 g of leaves. Such an amount of material can be collected even from a single plant and yields more than 30 Gb of ONT data in two MinION runs. We performed a comparative analysis of different genome assemblers and polishers on the gained data and obtained the final 447.1-Mb assembly of L. usitatissimum line 3896 genome using the Canu—Racon (two iterations)—Medaka combination. The genome comprised 1695 contigs and had an N50 of 6.2 Mb and a completeness of 93.8% of BUSCOs from eudicots_odb10. Our study highlights the impact of the chosen genome construction strategy on the resulting assembly parameters and its eligibility for future genomic studies.

Repeatome Analyses and Satellite DNA Chromosome Patterns in Deschampsia sukatschewii, D. cespitosa, and D. antarctica (Poaceae)

Subpolar and polar ecotypes of Deschampsia sukatschewii (Popl.) Roshev, D. cespitosa (L.) P. Beauv, and D. antarctica E. Desv. are well adapted to stressful environmental conditions, which make them useful model plants for genetic research and breeding. For the first time, the comparative repeatome analyses of subpolar and polar D. sukatschewii, D. cespitosa, and D. antarctica was performed using RepeatExplorer/TAREAN pipelines and FISH-based chromosomal mapping of the identified satellite DNA families (satDNAs). In the studied species, mobile genetic elements of class 1 made up the majority of their repetitive DNA; interspecific variations in the total amount of Ty3/Gypsy and Ty1/Copia retroelements, DNA transposons, ribosomal, and satellite DNA were revealed; 12–18 high confident and 7–9 low confident putative satDNAs were identified. According to BLAST, most D. sukatschewii satDNAs demonstrated sequence similarity with satDNAs of D. antarctica and D. cespitosa indicating their common origin. Chromosomal mapping of 45S rDNA, 5S rDNA, and satDNAs of D. sukatschewii allowed us to construct the species karyograms and detect new molecular chromosome markers important for Deschampsia species. Our findings confirmed that genomes of D. sukatschewii and D. cespitosa were more closely related compared to D. antarctica according to repeatome composition and patterns of satDNA chromosomal distribution.

Clarification of the Position of Linum stelleroides Planch. within the Phylogeny of the Genus Linum L

The phylogeny of members of the family Linaceae DC. ex Perleb has not been adequately studied. In particular, data on the phylogenetic relationship between Linum stelleroides Planch. and other representatives of the blue-flowered flax are very controversial. In the present work, to clarify this issue, we obtained DNA sequences of three nuclear loci (IGS and ITS1 + 5.8S rDNA + ITS2 of the 35S rRNA gene and the 5S rRNA gene) and eight chloroplast loci (rbcL, the trnL-trnF intergenic spacer, matK, the 3' trnK intron, ndhF, trnG, the psbA-trnH intergenic spacer, and rpl16) of 10 Linum L. species (L. stelleroides, L. hirsutum, L. perenne, L. leonii, L. lewisii, L. narbonense, L. decumbens, L. grandiflorum, L. bienne (syn. L. angustifolium), and L. usitatissimum) using high-throughput sequencing data. The phylogenetic analysis showed that L. stelleroides forms a basal branch in the blue-flowered flax clade. Previously found inconsistencies in the position of L. stelleroides and some other species in the Linaceae phylogenetic tree resulted from the erroneous species identification of some of the studied plant samples.

Genome Assembly and Sex-Determining Region of Male and Female Populus × sibirica

The genus Populus is presented by dioecious species, and it became a promising object to study the genetics of sex in plants. In this work, genomes of male and female Populus × sibirica individuals were sequenced for the first time. To achieve high-quality genome assemblies, we used Oxford Nanopore Technologies and Illumina platforms. A protocol for the isolation of long and pure DNA from young poplar leaves was developed, which enabled us to obtain 31 Gb (N50 = 21 kb) for the male poplar and 23 Gb (N50 = 24 kb) for the female one using the MinION sequencer. Genome assembly was performed with different tools, and Canu provided the most complete and accurate assemblies with a length of 818 Mb (N50 = 1.5 Mb) for the male poplar and 816 Mb (N50 = 0.5 Mb) for the female one. After polishing with Racon and Medaka (Nanopore reads) and then with POLCA (Illumina reads), assembly completeness was 98.45% (87.48% duplicated) for the male and 98.20% (76.77% duplicated) for the female according to BUSCO (benchmarking universal single-copy orthologs). A high proportion of duplicated BUSCO and the increased genome size (about 300 Mb above the expected) pointed at the separation of haplotypes in a large part of male and female genomes of P. × sibirica. Due to this, we were able to identify two haplotypes of the sex-determining region (SDR) in both assemblies; and one of these four SDR haplotypes, in the male genome, contained partial repeats of the ARR17 gene (Y haplotype), while the rest three did not (X haplotypes). The analysis of the male P. × sibirica SDR suggested that the Y haplotype originated from P. nigra, while the X haplotype is close to P. trichocarpa and P. balsamifera species. Moreover, we revealed a Populus-specific repeat that could be involved in translocation of the ARR17 gene or its part to the SDR of P. × sibirica and other Populus species. The obtained results expand our knowledge on SDR features in the genus Populus and poplar phylogeny.

Aluminum Responsive Genes in Flax (Linum usitatissimum L.)

Flax (Linum usitatissimum L.) is a multipurpose crop which is used for the production of textile, oils, composite materials, pharmaceuticals, etc. Soil acidity results in a loss of seed and fiber production of flax, and aluminum toxicity is a major factor that depresses plant growth and development in acid conditions. In the present work, we evaluated gene expression alterations in four flax genotypes with diverse tolerance to aluminum exposure. Using RNA-Seq approach, we revealed genes that are differentially expressed under aluminum stress in resistant (Hermes, TMP1919) and sensitive (Lira, Orshanskiy) cultivars and selectively confirmed the identified alterations using qPCR. To search for differences in response to aluminum between resistant and sensitive genotypes, we developed the scoring that allowed us to suggest the involvement of MADS-box and NAC transcription factors regulating plant growth and development and enzymes participating in cell wall modifications in aluminum tolerance in flax. Using Gene Ontology (GO) enrichment analysis, we revealed that glutathione metabolism, oxidoreductase, and transmembrane transporter activities are the most affected by the studied stress in flax. Thus, we identified genes that are involved in aluminum response in resistant and sensitive genotypes and suggested genes that contribute to flax tolerance to the aluminum stress.

Characterization of repeated DNA sequences in genomes of blue-flowered flax

BACKGROUND

Members of different sections of the genus Linum are characterized by wide variability in size, morphology and number of chromosomes in karyotypes. Since such variability is determined mainly by the amount and composition of repeated sequences, we conducted a comparative study of the repeatomes of species from four sections forming a clade of blue-flowered flax. Based on the results of high-throughput genome sequencing performed in this study as well as available WGS data, bioinformatic analyses of repeated sequences from 12 flax samples were carried out using a graph-based clustering method.

RESULTS

It was found that the genomes of closely related species, which have a similar karyotype structure, are also similar in the repeatome composition. In contrast, the repeatomes of karyologically distinct species differed significantly, and no similar tandem-organized repeats have been identified in their genomes. At the same time, many common mobile element families have been identified in genomes of all species, among them, Athila Ty3/gypsy LTR retrotransposon was the most abundant. The 30-chromosome members of the sect. Linum (including the cultivated species L. usitatissimum) differed significantly from other studied species by a great number of satellite DNA families as well as their relative content in genomes.

CONCLUSIONS

The evolution of studied flax species was accompanied by waves of amplification of satellite DNAs and LTR retrotransposons. The observed inverse correlation between the total contents of dispersed repeats and satellite DNAs allowed to suggest a relationship between both classes of repeating sequences. Significant interspecific differences in satellite DNA sets indicated a high rate of evolution of this genomic fraction. The phylogenetic relationships between the investigated flax species, obtained by comparison of the repeatomes, agreed with the results of previous molecular phylogenetic studies.

Flax (Linum usitatissimum L.) response to non-optimal soil acidity and zinc deficiency

BACKGROUND

Flax (Linum usitatissimum L.) is grown for fiber and seed production. Unfavorable environments, such as nutrient deficiency and non-optimal soil acidity, decrease the quantity and quality of yield. Cultivation of tolerant to stress varieties can significantly reduce the crop losses. Understanding the mechanisms of flax response to the stresses and identification of resistance gene candidates will help in breeding of improved cultivars. In the present work, the response of flax plants to increased pH level and zinc (Zn) deficiency was studied.

RESULTS

We performed high-throughput transcriptome sequencing of two flax cultivars with diverse tolerance to increased pH level and Zn deficiency: Norlin (tolerant) and Mogilevsky (sensitive). Sixteen cDNA libraries were created from flax plants grown under control conditions, increased pH level, Zn deficiency, and both stresses simultaneously, and about 35 million reads were obtained for each experiment type. Unfavorable pH resulted in significantly stronger gene expression alterations compared to Zn deficiency. Ion homeostasis, oxidoreductase activity, cell wall, and response to stress Gene Ontology terms were the most affected by unfavorable pH and Zn deficiency both in tolerant and sensitive flax cultivars. Upregulation of genes encoding metal transporters was identified under increased pH level, Zn deficiency, and both stresses simultaneously. Under Zn deficiency, only in tolerant cultivar Norlin, we revealed the induction of several photosynthesis-related genes and, in this way, this tolerant genotype could overcome unfavorable effects of reduced Zn content.

CONCLUSIONS

We identified genes with expression alterations in flax under non-optimal soil acidity and Zn deficiency based on high-throughput sequencing data. These genes are involved in diverse processes, including ion transport, cell wall biogenesis, and photosynthesis, and could play an important role in flax response to the studied stresses. Moreover, genes with distinct expression changes between examined tolerant and sensitive genotypes could determine the mechanisms of flax tolerance to non-optimal soil acidity and Zn deficiency.

[CAX3 Gene is Involved in Flax Response to High Soil Acidity and Aluminum Exposure]

Understanding the molecular mechanisms of plant response to unfavorable conditions is necessary for the effective selection of tolerant genotypes. Earlier, using high-throughput transcriptome sequencing of flax plants after exposure to aluminum ions (Al^(3+)) and high soil acidity, we detected stress-induced alteration in the expression of several genes, including CAX3, which encodes Ca^(2+)/H^(+)-exchanger involved in calcium ion transport. Here we describe CAX3 mRNA levels in flax cultivars either tolerant (Hermes and TMP1919) or sensitive (Lira and Orshanskiy) to Al^(3+). Stress-induced increased expression of CAX3 was detected only in aluminum-tolerant flax cultivars. The product of CAX3 gene may participate in flax response to high soil acidity and high Al^(3+) concentration through Ca^(2+)-mediated intracellular regulation.

Differential gene expression in response to Fusarium oxysporum infection in resistant and susceptible genotypes of flax (Linum usitatissimum L.)

BACKGROUND

Flax (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.

RESULTS

The transcriptomes of two resistant and two susceptible flax cultivars with respect to Fusarium wilt, as well as two resistant BC₂F₅ 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 BC₂F₅ populations resistant to Fusarium wilt than in susceptible cultivars in response to F. oxysporum infection.

CONCLUSIONS

Using 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 BC₂F₅ populations showing induced expression in response to F. oxysporum infection are the most promising resistance gene candidates.

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.

Integration of Physical, Genetic, and Cytogenetic Mapping Data for Cellulose Synthase (CesA) Genes in Flax (Linum usitatissimum L.)

Flax, Linum usitatissimum L., is a valuable multi-purpose plant, and currently, its genome is being extensively investigated. Nevertheless, mapping of genes in flax genome is still remaining a challenging task. The cellulose synthase (CesA) multigene family involving in the process of cellulose synthesis is especially important for metabolism of this fiber crop. For the first time, fluorescent in situ hybridization (FISH)-based chromosomal localization of the CesA conserved fragment (KF011584.1), 5S, and 26S rRNA genes was performed in landrace, oilseed, and fiber varieties of L. usitatissimum. Intraspecific polymorphism in chromosomal distribution of KF011584.1 and 5S DNA loci was revealed, and the generalized chromosome ideogram was constructed. Using BLAST analysis, available data on physical/genetic mapping and also whole-genome sequencing of flax, localization of KF011584.1, 45S, and 5S rRNA sequences on genomic scaffolds, and their anchoring to the genetic map were conducted. The alignment of the results of FISH and BLAST analyses indicated that KF011584.1 fragment revealed on chromosome 3 could be anchored to linkage group (LG) 11. The common LG for 45S and 5S rDNA was not found probably due to the polymorphic localization of 5S rDNA on chromosome 1. Our findings indicate the complexity of integration of physical, genetic, and cytogenetic mapping data for multicopy gene families in plants. Nevertheless, the obtained results can be useful for future progress in constructing of integrated physical/genetic/cytological maps in L. usitatissimum which are essential for flax breeding.

Comparative molecular cytogenetic characterization of seven Deschampsia (Poaceae) species

The genus Deschampsia P. Beauv (Poaceae) involves a group of widespread polymorphic species. Some of them are highly tolerant to stressful and variable environmental conditions, and D. antarctica is one of the only two vascular plants growing in Antarctic. This species is a source of useful for selection traits and a valuable model for studying an environmental stress tolerance in plants. Genome diversity and comparative chromosomal phylogeny within the genus have not been studied yet as karyotypes of most Deschampsia species are poorly investigated. We firstly conducted a comparative molecular cytogenetic analysis of D. antarctica (Antarctic Peninsula) and related species from various localities (D. cespitosa, D. danthonioides, D. elongata, D. flexuosa (= Avenella flexuosa), D. parvula and D. sukatschewii by fluorescence in situ hybridization with 45S and 5S rDNA, DAPI-banding and sequential rapid in situ hybridization with genomic DNA of D. antarctica, D. cespitosa, and D. flexuosa. Based on patterns of distribution of the examined markers, chromosomes of the studied species were identified. Within these species, common features as well as species peculiarities in their karyotypic structure and chromosomal distribution of molecular cytogenetic markers were characterized. Different chromosomal rearrangements were detected in D. antarctica, D. flexuosa, D. elongata and D. sukatschewii. In karyotypes of D. antarctica, D. cespitosa, D. elongata and D. sukatschewii, 0-3 B chromosomes possessed distinct DAPI-bands were observed. Our findings suggest that the genome evolution of the genus Deschampsia involved polyploidy and also different chromosomal rearrangements. The obtained results will help clarify the relationships within the genus Deschampsia, and can be a basis for the further genetic and biotechnological studies as well as for selection of plants tolerant to extreme habitats.

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

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

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.

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.

[Genome: Origins and evolution of the term]

The appearance of a new scientific term is a significant event in the human cognitive process and the result of the realization of the separateness of an object or a phenomenon. Our article concentrates on the origins of basic genetic terms, such as genetics, gene, genotype, genome, gene pool, and genomics. We propose using the term karyogenomics for the special direction of genomics related to the study of the organization and evolution of eukaryotic genomes by means of modern chromosome analysis, as well as by full genome sequencing.

A comparative cytogenetic study of Drosophila parasitoids (Hymenoptera, Figitidae) using DNA-binding fluorochromes and FISH with 45S rDNA probe

Karyotypes of Leptopilina boulardi (Barbotin, Carton et Keiner-Pillault, 1979) (n = 9), L. heterotoma (Thomson, 1862) (n = 10), L. victoriae Nordlander, 1980 (n = 10) and Ganaspis xanthopoda (Ashmead, 1896) (n = 9) (Hymenoptera, Figitidae) were studied using DNA-binding ligands with different base specificity [propidium iodide (PI), chromomycin A3 (CMA3) and 4',6-diamidino-2-phenylindole (DAPI)], and fluorescence in situ hybridization (FISH) with a 45S rDNA probe. Fluorochrome staining was similar between the different fluorochromes, except for a single CMA3- and PI-positive and DAPI-negative band per haploid karyotype of each species. FISH with 45S rDNA probe detected a single rDNA site in place of the bright CMA3-positive band, thus identifying the nucleolus organizing region (NOR). Chromosomal locations of NORs were similar for both L. heterotoma and L. victoriae, but strongly differed in L. boulardi as well as in G. xanthopoda. Phylogenetic aspects of NOR localization in all studied species are briefly discussed.

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.

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.

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.

Comparative cytogenetic study on two species of the genus Entedon Dalman, 1820 (Hymenoptera, Eulophidae) using DNA-binding fluorochromes and molecular and immunofluorescent markers

Karyotypes of Entedon cionobius Thomson, 1878 and Entedon cioni Thomson, 1878 (Hymenoptera: Eulophidae) were studied using DNA-binding ligands with different base specificity (propidium iodide, chromomycin A3, methyl green and DAPI; all these ligands, except for the last one, were used for the first time in parasitic wasps), C-banding, fluorescence in situ hybridization (FISH) with a 45S rDNA probe and 5-methylcytosine immunodetection. Female karyotypes of both species contain five pairs of relatively large metacentric chromosomes and a pair of smaller acrocentric chromosomes (2n = 12). As in many other Hymenoptera, males of both Entedon Dalman, 1820 species have haploid chromosome sets (n = 6). Fluorochrome staining revealed chromosome-specific banding patterns that were similar between the different fluorochromes, except for the CMA3- and PI-positive and DAPI-negative band in the pericentromeric regions of the long arms of both acrocentric chromosomes. The obtained banding patterns were virtually identical in both species and allowed for the identification of each individual chromosome. C-banding revealed a pattern similar to DAPI staining, although centromeric and telomeric regions were stained more intensively using the former technique. FISH detected a single rDNA site in the same position on the acrocentric chromosomes as the bright CMA3-positive band. Immunodetection of 5-methylcytosine that was performed for the first time in the order Hymenoptera revealed 5-methylcytosine-rich sites in the telomeric, centromeric and certain interstitial regions of most of the chromosomes.

Comparison of genomes of eight species of sections Linum and Adenolinum from the genus Linum based on chromosome banding, molecular markers and RAPD analysis

Karyotypes of species sects. Linum and Adenolinum have been studied using C/DAPI-banding, Ag-NOR staining, FISH with 5S and 26S rDNA and RAPD analysis. C/DAPI-banding patterns enabled identification of all homologous chromosome pairs in the studied karyotypes. The revealed high similarity between species L. grandiflorum (2n = 16) and L. decumbens by chromosome and molecular markers proved their close genome relationship and identified the chromosome number in L. decumbens as 2n = 16. The similarity found for C/DAPI-banding patterns between species with the same chromosome numbers corresponds with the results obtained by RAPD-analysis, showing clusterization of 16-, 18- and 30-chromosome species into three separate groups. 5S rDNA and 26S rDNA were co-localized in NOR-chromosome 1 in the genomes of all species investigated. In 30-chromosome species, there were three separate 5S rDNA sites in chromosomes 3, 8 and 13. In 16-chromosome species, a separate 5S rDNA site was also located in chromosome 3, whereas in 18-chromosome species it was found in the long arm of NOR-chromosome 1. Thus, the difference in localization of rDNA sites in species with 2n = 16, 2n = 30 and 2n = 18 confirms taxonomists opinion, who attributed these species to different sects. Linum and Adenolinum, respectively. The obtained results suggest that species with 2n = 16, 2n = 18 and 2n = 30 originated from a 16-chromosome ancestor.

Reverse transcriptase inhibitors suppress telomerase function and induce senescence-like processes in cultured mouse fibroblasts

Spontaneous transformation of mouse embryonic fibroblasts in the presence of the reverse transcriptase inhibitors azidothymidine and carbovir led to the formation of telomerase-free clones. After prolonged cultivation of fibroblasts in the presence of carbovir, resistant cells with a very high level of telomerase activity were obtained. Azidothymidine and carbovir, but not dideoxycytidine, induced senescence-like processes in cultures of immortal mouse fibroblasts. After long-term incubation, cell proliferation gradually decreased, their morphology becoming similar to that of the senescent ones. The process was reversible: after inhibitor removal, the cells, including the giant ones, entered mitoses. All these data suggest that reverse transcriptase inhibitors block telomerase function in mouse cells.

Formation of a synthetic karyotype of tetraploid triticale

Thirteen forms of 4x triticale were analyzed using a C-banding technique. All forms contained 14 R-genome chromosomes and 69 different combinations of A-, B-, and D-genome chromosomes from wheat were identified. The following general features of the formation of the synthetic karyotype of 4x triticale were found: there was a different rate at which homoelogous groups were stabilized; there was a selection against "hétérozygotes" (AB); and there was an effect exerted by the chromosomal constitution of some homoeologous groups on the rate and direction of stabilization in other groups. The two most probable variants of the wheat genome in 4x triticale are 1B, 2B, 3A, 4A, 5A, 6A, 7B and the reverse. The selection occurred at the level of homoeologues, and the selective advantages are determined by genotype-environment interactions. The environmental factors included both natural and artificial selection pressures. The data presented here can be used for the creation of a hypothesis concerning the practical approaches in chromosomal engineering of cereals.Key words: 4x triticale, general features of karyotype formation, C-banding.

Chromosomal alterations in the karyotype of triticale in comparison with the parental forms : 2. Heterochromatin of the wheat chromosomes

Wheat chromosomes of the primary winter hexaploid and octoploid triticales and of the parental durum and common wheat varieties were studied using morphometric analysis. The size of some heterochromatic segments was shown to change in triticale. Telomeric and intercalary C-bands both increased and decreased in size whereas centromeric bands only increased. The size variability of C-bands in triticale B-genome chromosomes decreased in most of the cases and increased only for several specific C-bands. The C-bands of homologous B-genome chromosomes changed in the same direction in both triticale forms. The changes in size of the C-bands found in R-genome chromosomes detected earlier in these triticale forms (Badaeva et al. 1986) were shown to coincide in their pattern with the size changes of C-bands in homeological B-genome chromosomes. Our data are indicative of regular, directed chromosomal changes in the triticale karyotype.

Chromosome alterations in the karyotype of triticale in comparison with the parental forms : 1. Heterochromatic regions of R genome chromosomes

R genome chromosomes were studied in two forms of primary triticales (hexaploid 'TPG-1/1-78' and octoploid 'AD 825') and in their parent rye forms (Secale cereale L. cv. 'Kharkovskaya 60' and 'VSKhI', respectively) using the methods of C-banding and morphometric analysis. The size of some heterochromatic segments was shown to alter in the karyotype of triticale. An increase in size was detected approximately in half of all telomeric C-bands; the size of the other C-bands either decreased or did not change. The frequencies of these alterations were 1∶1. The variability in the size of telomeric C-bands in rye chromosomes diminished in both triticales studied. The two triticale forms inherited variants of R genome chromosome polymorphism predominantly with the medium size range of telomeric C-bands. The centromeric C-bands in both triticale forms either enlarged or did not alter. Possible mechanisms responsible for the observed pattern of alterations are discussed.

Cytogenetic analysis of forms produced by crossing hexaploid triticale with common wheat

Hexaploid triticales were crossed with common wheats, and the resultant froms were selected for either triticale (AD 213/5-80) or common wheat (lines 381/80, 391/80, 393/80). The cytogenetic analysis showed that all forms differ in their chromosome composition. Triticale AD 213/5-80 and wheat line 381/80 were stable forms with 2n = 6x = 42. Lines 391/80 and 393/80 were cytologically unstable. In triticale AD 213/5-80, a 2R (2D) chromosome substitution was found. Each of the three wheat lines had a chromosome formed by the translocation of the short arm of IR into the long arm of the IB chromosome. In line 381/80, this chromosome seems to be inherited from the 'Kavkaz' wheat variety. In lines 391/80 and 393/80, this chromosome apparently formed de novo since the parent forms did not have it. The karyotype of line 381/80 was found to contain rye chromosomes 4R/7R, 5R and 7R/4R. About 15% of the cells in line 391/80 contained an isochromosome for the 5R short arm and also a chromosome which arose from the translocation of the long arms of the 5D and 5R chromosomes. About one-third of the cells in the common wheat line 393/80 contained the 5R chromosome. This chromosome was normal or rearranged. Practical applications of the C-banding technique in the breeding of triticale is discussed.