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

Phenotypic analysis of Longya-10 × pale flax hybrid progeny and identification of candidate genes regulating prostrate/erect growth in flax plants

Flax is a dual-purpose crop that is important for oil and fiber production. The growth habit is one of the crucial targets of selection during flax domestication. Wild hybridization between cultivated flax and wild flax can produce superior germplasms for flax breeding and facilitate the study of the genetic mechanism underlying agronomically important traits. In this study, we used pale flax, Linum grandiflorum, and L. perenne to pollinate Longya-10. Only pale flax interspecific hybrids were obtained, and the trait analysis of the F₁ and F₂ generations showed that the traits analyzed in this study exhibited disparate genetic characteristics. In the F₁ generation, only one trait, i.e., the number of capsules per plant (140) showed significant heterosis, while the characteristics of other traits were closely associated with those of the parents or a decline in hybrid phenotypes. The traits of the F₂ generation were widely separated, and the variation coefficient ranged from 9.96% to 146.15%. The quantitative trait locus underlying growth habit was preliminarily found to be situated on chromosome 2 through Bulked-segregant analysis sequencing. Then linkage mapping analysis was performed to fine-map GH2.1 to a 23.5-kb interval containing 4 genes. Among them, L.us.o.m.scaffold22.109 and L.us.o.m.scaffold22.112 contained nonsynonymous SNPs with Δindex=1. Combined with the qRT-PCR results, the two genes might be possible candidate genes for GH2.1. This study will contribute to the development of important germplasms for flax breeding, which would facilitate the elucidation of the genetic mechanisms regulating the growth habit and development of an ideal architecture for the flax plant.

Transposable Elements in the Revealing of Polymorphism-Based Differences in the Seeds of Flax Varieties Grown in Remediated Chernobyl Area

The nuclear reactor accident in Chernobyl, Ukraine, resulted in effects both locally and farther away. Most of the contaminated areas were the agricultural fields and forests. Experimental fields were established near Chernobyl—radioactively contaminated fields localized 5 km from Chernobyl Nuclear Power Plant as well as the remediated soil that is localized directly in the Chernobyl town. Two flax varieties growing under chronic exposition to ionizing radiation were used for this study—the local Ukrainian variety Kyivskyi and a commercial variety Bethune. The screening of the length polymorphism generated by transposable elements insertions were performed. All known types of common flax transposon, retrotransposons and iPBS approach were used. In the iPBS multiplex analyze, for the Kyivskyi variety, a unique addition was found in the seeds from the radioactive contaminated field and for the Bethune variety, a total of five amplicon additions were obtained and one deletion. For the TRIM Cassandra fingerprints, two amplicon additions were generated in the seeds from radioactive contaminated fields for the Bethune variety. In summary, the obtained data represent the genetic diversity between control and irradiated subgroups of flax seeds from Chernobyl area and the presence of activated transposable elements due to the irradiation stress.

Phylogenetic Analysis and Flower Color Evolution of the Subfamily Linoideae (Linaceae)

The taxonomy of the subfamily Linoideae at the intergeneric and section levels has been questioned throughout the years, and the evolution of floral characters remains poorly understood. In particular, the evolution of flower color is still uncertain, despite its ecological importance and being one of the most variable and striking traits in Angiospermae. We evaluated the phylogenetic relationships of the genera and sections and used the phylogeny to reconstruct the ancestral state of flower color. The results suggest reevaluating the taxonomic status of segregated genera and re-incorporating them into Linum. Four of the five sections currently accepted were recovered as monophyletic (Cathartolinum, Dasylinum, Linum, and Syllinum). We propose accepting the section Stellerolinon and reevaluating Linopsis, whose representatives were recovered in three separate clades. The ancestral flower color for Linoideae was yellow-white. The flower colors purple and yellow-white were recovered at the deepest nodes of the two main clades. Pink, blue, and red colors were the most recent to evolve. These results appear to be related to diversification events, biogeographical history, and ecological aspects of the subfamily. Our reconstruction constitutes the first plausible scenario that explores the evolution of flower color, leading to new testable hypotheses for future research on the flax group.

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.

Sex-specific polymorphism of MET1 and ARR17 genes in Populus × sibirica

The genus Populus is an effective model in tree genetics. This genus includes dioecious species and, recently, whole genome resequencing of P. trichocarpa and P. balsamifera enabled the identification of sex-linked regions and sex-associated single nucleotide polymorphisms (SNPs). These results created new opportunities to study sex determination in poplars. In the present work, we performed deep sequencing of genes encoding METHYLTRANSFERASE1 (MET1) and homolog of ARABIDOPSIS RESPONSE REGULATOR 17 (ARR17), which are localized in a sex-linked region of Populus genome and contain a number of sex-associated SNPs. Amplicon libraries for 38 samples of P. × sibirica (19 males and 19 females) were sequenced on MiSeq Illumina (300 nt paired-end reads) and approximately 4000× coverage was obtained for each sample. In total, from 80 to 179 SNPs were detected in poplar individuals for MET1, and from 16 to 49 SNPs were detected for ARR17. We identified 17 sex-specific SNPs (11 in MET1 and 6 in ARR17) - they were present in all males but absent in all females. For identified sex-specific SNP sites, females were homozygous, while males were heterozygous. Moreover, colocation of sex-specific SNPs confirming the XY sex-determination system of poplars was revealed: in one allelic variant, males had the same nucleotides as females, while in the other, sex-specific SNPs were present. Based on the data obtained, we developed and successfully applied a high-resolution melting-based approach for sex identification in poplars. The developed molecular markers are useful for distinguishing between male and female poplars in scientific research and can also be applied to select male-only genotypes for use in city landscaping and production of paper, pulp, and biofuel.

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.

Genetics polymorphism of poplars from Moscow region based on high-throughput sequencing of ITS

Poplars are widely used in landscaping of Moscow due to the ability to effectively purify the air from harmful impurities and to release a large amount of oxygen. The genusPopulusis characterized by a high level of intraspecies polymorphism, as well as the presence of natural interspecies hybrids. The aim of our work was to evaluate the genetic diversity of poplars, which are growing on the territory of Moscow city by high-throughput sequencing of internal transcribed spacers of 45S rRNA genes (ITS sequences). Sequencing of ITS of 40 poplar plants was performed on Illumina platform (MiSeq) and about 3 000 reads were obtained for each sample in average. Bioinformatics analysis was performed using CLC Genomics Workbench tool. The involved set of poplars had a high level of genetic diversity – the number of single nucleotide polymorphisms (SNPs) detected in each genotype relative to the reference ITS1 and ITS2 sequences ofP. trichocarpavarying from 4 to 44. We showed that even trees which had been planted on the same territory and, probably, at the same time had significant genetic differences. It can be speculated that highly polymorphic plant material was used for planting poplars in Moscow. For some sites with SNPs, several variants of nucleotides were found in the same individual and the ratio of SNPs was different. We assume that close to 50/50 ratio is observed in interspecific hybrids due to genetic differences in the ITS sequences between maternal and paternal genotypes. For SNPs with a predominance of one of the variants, the presence of paralogues among numerous genomic copies of ITS sequences is more likely. The results of our work can provide a framework for molecular genetic markers application with the purpose ofPopulusspecies and interspecific hybrids identification, determination the origin of a number of natural hybrids, and monitoring the diversity of genusPopulusin the Moscow city.

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.

Useful parasites: the evolutionary biology and biotechnology applications of transposable elements

Transposable elements usually comprise the most abundant nongenic fraction of eukaryotic genomes. Because of their capacity to selfreplicate and to induce a wide range of mutations, transposable elements have long been considered as 'parasitic' or 'selfish'. Today, we recognize that the findings about genomic changes affected by transposable elements have considerably altered our view of the ways in which genomes evolve and work. Numerous studies have provided evidences that mobile elements have the potential to act as agents of evolution by increasing, rearranging and diversifying the genetic repertoire of their hosts. With large-scale sequencing becoming increasingly available, more and more scientists come across transposable element sequences in their data. I will provide examples that transposable elements, although having signatures of 'selfish' DNA, play a significant biological role in the maintainance of genome integrity and providing novel regulatoty networks. These features, along with the transpositional and mutagenic capacity to produce a raw genetic diversity, make the genome mobile fraction, a key player in species adaptation and microevolution. The last but not least, transposable elements stand as informative DNA markers that may complement other conventional DNA markers. Altogether, transposable elements represent a promising, but still largely unexplored research niche and deserve to be included into the agenda of molecular ecologists, evolutionary geneticists, conservation biologists and plant breeders.

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 AlCl₃ 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.

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.

Multiplexed shotgun sequencing reveals congruent three-genome phylogenetic signals for four botanical sections of the flax genus Linum

A genome-wide detection of phylogenetic signals by next generation sequencing (NGS) has recently emerged as a promising genomic approach for phylogenetic analysis of non-model organisms. Here we explored the use of a multiplexed shotgun sequencing method to assess the phylogenetic relationships of 18 Linum samples representing 16 species within four botanical sections of the flax genus Linum. The whole genome DNAs of 18 Linum samples were fragmented, tagged, and sequenced using an Illumina MiSeq. Acquired sequencing reads per sample were further separated into chloroplast, mitochondrial and nuclear sequence reads. SNP calls upon genome-specific sequence data sets revealed 6143 chloroplast, 2673 mitochondrial, and 19,562 nuclear SNPs. Phylogenetic analyses based on three-genome SNP data sets with and without missing observations showed congruent three-genome phylogenetic signals for four botanical sections of the Linum genus. Specifically, two major lineages showing a separation of Linum-Dasylinum sections and Linastrum-Syllinum sections were confirmed. The Linum section displayed three major branches representing two major evolutionary stages leading to cultivated flax. Cultivated flax and its immediate progenitor were formed as its own branch, genetically more closely related to L. decumbens and L. grandiflorum with chromosome count of eight, and distantly apart from six other species with chromosome count of nine. Five species of the Linastrum and Syllinum sections were genetically more distant from cultivated flax, but they appeared to be more closely related to each other, even with variable chromosome counts. These findings not only provide the first evidence of congruent three-genome phylogenetic pathways within the Linum genus, but also demonstrate the utility of the multiplexed shotgun sequencing in acquisition of three-genome phylogenetic signals of non-model organisms.

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.

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.