Evolution of blue-flowered species of genus Linum based on high-throughput sequencing of ribosomal RNA genes
- Nadezhda L. Bolsheva†1,
- Nataliya V. Melnikova†1,
- Ilya V. Kirov2,
- Anna S. Speranskaya3,
- Anastasia A. Krinitsina3,
- Alexey A. Dmitriev1,
- Maxim S. Belenikin3,
- George S. Krasnov1,
- Valentina A. Lakunina1,
- Anastasiya V. Snezhkina1,
- Tatiana A. Rozhmina1, 4,
- Tatiana E. Samatadze1,
- Olga Yu. Yurkevich1,
- Svyatoslav A. Zoshchuk1,
- Аlexandra V. Amosova1Email author,
- Anna V. Kudryavtseva1 and
- Olga V. Muravenko1
© The Author(s). 2017
Published: 28 December 2017
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.
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).
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.
The genus Linum (Linacea), which involves about 200 species, is subdivided into a group of sections based on morphological traits [1–3]. The investigation of the phylogeny of the genus Linum has been performed several times by means of different molecular and phylogenetic approaches. Particularly, due to the similar results obtained by AFLP, RAPD and SSAP methods, the members of the genus could be subdivided into the groups of species with the similar or even identical karyotype structure [4–7]. In some instances, the groups species revealed according to the genome similarities coincide with the separate sections of the genus with the exception of the sect. Linum and Linastrum involving several groups. Unfortunately, DNA fingerprinting is not an adequate method for studying the phylogenetic relationships between the separate groups of species. Mostly, these relationships were revealed by molecular phylogenetic investigations performed with the use of nuclear ibosomal ITS1 and ITS2 and plastid sequences of flax . The results of these studies showed that the genus Linum was not monophyletic but combined the members of two sister clades: yellow-flowered and blue-flowered flaxes. Later, the analysis of the transcriptomes of 11 flax species representing different sections was performed with the use of high-throughput sequencing . Based on the analysis of 413 nuclear genes, the phylogeny of these 11 species was constructed, and it was almost identical to Linum phylogeny created with the use of nuclear ribosomal ITS and chloroplast genes . The similar results were obtained by multiplexed shotgun sequencing of a great number of nuclear and chloroplast genes of 16 flax species . Thus, the phylogenetic relationships within the genus Linum have been adequately investigated, however, a number of details are still unresolved. Particularly, the species relationships within the monotypic sect. Stellerolinum and the sect. Dasylinum need further clarification. According to the phylogeny of rbcL, L. stelleroides (sect. Stellerolinum) is a sister to sect. Dasylinum. Besides, the analysis of ITS and other chloroplast data showed that this species is a sister to Dasylinum plus the remaining members of sect. Linum including the L. perenne group . The phylogeny of the sect. Adenolinum (syn. L. perenne group) also needs some refinements. Based on the morphological traits, about 20 species can be separated within the sect. Adenolinum. However, this identification of the species are not unambiguous as the morphological traits are mainly the quantitative characters which depend on the environment conditions .
It was found that all members of the section have similar karyotypes. Moreover, genomic differences among the certain species of the section were not revealed by AFLP, RAPD or SSAP methods [4, 5, 7]. As the analysis of the sequences of nuclear and cytoplasmic genes was performed only for a limited number of species [8–10], it is not possible to estimate accurately the intra- or interspecific divergence of these sequences among the members of the sect. Adenolinum, and the question of species independence within this section still remains unanswered.
In the present study, we investigated the relationships within the clade of blue-flowered flaxes based on the analysis of ITS1-5.8S-ITS2 sequences of 35S rRNA genes and also 5S rRNA genes. Moreover, in angiosperm genomes, 35S and 5S rRNA genes are usually localized separately as multiple tandem repeats. The coding sequences of rRNA genes are highly conserved unlike the sequences of non-coding spacers which are polymorphic and therefore, can be widely used in phylogenetic studies [12, 13]. The distinct advantage of using of 5S rRNA genes as molecular markers in phylogenetic studies is the fact that these genes are inherited from the both parent species in hybrid plants. In contrast, 35S rRNA genes of one of the parents can be eliminated in the hybrid genome [14, 15]. For analysis of the selected gene sequences, we used the high-throughput target sequencing. This is a more adequate technique for estimation of intraspecific variability and, consequently, a degree of interspecific differences compared to the traditional approach which is based on the sequencing of few cloned sequences. For clarification of the phylogenetic relationships within the sect. Adenolinum, we increased significantly the number of the studied species and specimens. Additionally, for all the examined plant samples, the karyotypic studies and chromosome mapping of 5S and 35S rRNA genes were carried out.
In the present study, 43 flax accessions belonging to 23 species and subspecies from the sections Linum, Adenolinum, Dasylinum, Stellerolinum and Syllinum were used (Additional file 1). Most of these accessions were obtained from the genebank of Leibniz Institute of Plant Genetics and Crop Plant Research (IPK) (Gatersleben, Germany) and the seed collection of All-Russian Flax Institute (VNIIL) (Torzhok, Russian Federation). The accessions of L. amurense and L. stelleroides were kindly provided by Dr. L.N. Mironova, Botanic Garden Institute of the Far-Eastern Branch of the Russian Academy of Sciences (BGI of RAS) (Vladivostok, Russian Federation). Several accessions were collected from natural populations by Dr. A.A. Svetlova, Komarov Botanical Institute RAS (St. Petersburg, Russian Federation), by Dr. N.L. Bolsheva, Engelhardt Institute of Molecular Biology RAS (Moscow, Russian Federation), and by Dr. M. Pavelka, Euroseeds (Novy Jicin, Chech Republic).
Chromosome preparation and FISH
Preparation and sequencing of DNA libraries
Total DNA was extracted from flax leaves as described earlier .
First-stage PCR specific primers
Second-stage PCR primers 
The PCR products were separated in 1.5% agarose gel for 5 h at 45 V using ТВЕ buffer and then stained with ethidium bromide. The Bio-Rad Gel Doc system was used for visual detection of the presence of the PCR products. Bands of the expected sizes (~350 bp for 5S rDNA and ~950 for ITS1-5.8S rDNA-ITS2) were excised from agarose gel and purified using Agencourt AMPure beads (Beckman Coulter, USA). Purified amplicon libraries concentration was evaluated using Quant-iT PicoGreen dsDNA Assay Kit (Life technologies, USA) with the use of QuantiFlour fluorometer (Promega, USA). Then, equimolar pooling of barcoded amplicons for 5S rDNA and ITS1-5.8S rDNA-ITS2 libraries was performed and library quality was evaluated using Agilent 2100 Bioanalyzer (Agilent Technologies, USA). After library dilution, emulsion PCR for Lib-A library and 454 sequencing by GS Junior were performed according to the protocols provided by the manufacturer (Roche, Switzerland).
Phylogenetic analysis and statistical evaluation
The obtained reads of each accession were aligned and then clustered with the use of Allele Builder (http://sourceforge.net/p/allelebuilder). For the clusters representing more than 5% of the obtained reads, consensus sequences were constructed using the program CLC Genomics Workbench v6.5 (http://www.clcbio.com). The alignment of the obtained consensus sequences was performed using the program Clustal W . For ITS1-5.8S rDNA-ITS2 sequences, the parameters for Pairwise Sequence Alignment were: Gap Opening Penalty - 15, Gap Extension Penalty - 6.66; the parameters for Multiple Sequence Alignment were: Gap Opening Penalty - 15; Gap Extension Penalty - 6.66; DNA Weight Matrix - JUB; Transition Weight - 0.5; Delay Divergent Cutoff - 50%. For 5S rDNA, the parameters for Pairwise Sequence Alignment were: Gap Opening Penalty - 15, Gap Extension Penalty - 6.66; parameters for Multiple Sequence Alignment were: Gap Opening Penalty - 12; Gap Extension Penalty - 6.66; DNA Weight Matrix - JUB; Transition Weight - 0.4; Delay Divergent Cutoff −70%.
Phylogenetic analysis was performed using Maximum Likelihood (ML) by PhyML software . The best fit model of nucleotide substitution was determined by JModelTest 220.127.116.11 software . The selected model was YIM3ef + G and TPM3 + I + G for ITS and 5S regions, respectively, established by corrected Akaike Information Criterion (AICc) and Bayesian Information Criterion (BIC). The tree was visualized by iTOL tool . The node support was assessed by the Bayesian-like modification of the approximate likelihood ratio test (aBayes) . Statistical calculations were conducted in MEGA version 6 .
Karyotype structure and chromosomal mapping of rRNA genes
In the present study, we describe mostly the results of cytogenetic studies of the specimens which have not been published yet as the karyotype structure and localization of ribosomal genes of many flax species had been previously published [5, 6, 16].
The members of the section Dasylinum had the largest chromosomes found in flax (3.5-6 μm). Among three studied specimens, L. hirsutum ssp. hirsutum was diploid with 2n = 16, whereas L. hirsutum ssp. anatolicum and L. hirsutum ssp. pseudoanatolicum were tetraploid (2n = 32) containing two similar sets of chromosomes. In genomes of these plants, two loci of 35S rDNA were detected in the distal regions of the short arms of chromosomes 1 and 5, and one locus of 5S rDNA was found in the median region of the short arm of chromosome 4 (Fig. 1b, c).
The karyotype of L. narbonense possessed two similar chromosome sets and was probably tetraploid with 2n = 4× = 28. The sizes of chromosomes of this species were within 3-5 μm range. The clusters of 35S rDNA were localized in the proximal regions of short arms of two homeologous pairs of chromosomes. Large 5S rDNA clusters were found in the distal regions of two chromosome pairs and also in the pericentromeric heterochromatic regions of two pairs of metacentric chromosomes (Fig. 1d). In some metaphase cells, minor sites of 5S rDNA were detected in the pericentromeric heterochromatin of 2-4 chromosome pairs.
Closely related species L. grandiflorum and L. decumbens (2n = 16) had the similar structure of karyotypes, but their chromosomes were smaller (2-3.5 μm) than chromosomes of L. narbonense. The site of 35S rDNA was localized in the pericentromeric region of the long arm of chromosome 1, and the site of 5S rDNA was revealed in the distal region of the long arm of chromosome 3 (Fig. 2c, d).
The cultivated species L. usitatissimum as well as its wild ancestor L. angustifolium were mesotetraploids with 2n = 30 and had significantly smaller sizes of chromosomes (1.2-3 μm) compared to the species mentioned above. The only site of 35S rDNA was observed in the proximal region of the short arm of metacentric chromosome 1. Sites of 5S rDNA were found in the proximal region of the long arm of chromosome 3 and also in the distal regions of long arms of chromosomes 8 and 13 (Fig. 2e, f).
Structure and intragenomic variability of rDNA repeats
The PCR amplified region of 35S rRNA gene involved 33 bp of 5’end of 18S rRNA gene, ITS1 sequence, the sequence coding 5.8S rRNA gene, and ITS2 sequence (Additional file 2 ). It was found that intraspecific variability of the sequences of intergenic spacers ITS1 and ITS2 was significantly higher than the variability of coding regions of 5.8S rRNA gene, and it was caused mostly by single nucleotide substitutions (SNPs) and indels. Mutations of more extended regions of the sequences occurred rarely with less than 1% frequency. The analysis of the SNPs showed that the number of transitions exceeded the number of transversions. The clustering of ITS1-5.8S rDNA-ITS2 sequences revealed one major cluster including more than 5% of reads for each studied sample.
The PCR amplified sequences of 5S rRNA genes contained 98 bp of 5′ end of coding region, intergenic spacer sequence and 22 bp of 3′ end of coding region. In all studied specimens, the length of coding regions of 5S rRNA genes was 120 bp, and the length of intergenic spacers of 5S rRNA genes ranged from 136 to 239 (Additional file 3). Intraspecific variability of the spacer sequences was significantly higher than the variability of the coding regions. The variability was caused mostly by SNPs, indels and also by deletions or duplications of more extended regions of spacer sequences. The analysis of the SNPs showed that the number of transitions far exceeded the number of transversions (Additional file 4). For 5S rRNA sequences, the number of clusters per sample varied from 1 to 7. The analysis of the nature of mutations that defined clusters showed that the SNPs, in particular C➝T and G➝A transitions, were the most frequent polymorphisms in the majority of the specimens with the exception of the members of the section Dasylinum in which the variability of spacer sequences of 5S rRNA genes were caused mostly by deletions of different spacer regions and also by replacement of CGGAATGGGAA sequence in position 309-319 to TAAAATAATAT sequence (a subrepeat of the spacer) in all studied members of the section Dasylinum.
Interspecies variability of rDNA repeats
Comparison of ITS1-5.8S rDNA-ITS2 sequences showed that 5.8S rDNA sequences were rather conservative in all studied flax specimens (Average Evolutionary Divergence 2,92 + 1,18) whereas strong sequence divergence was found in ITS1 and ITS2 sequences (Average Evolutionary Divergence for ITS1and ITS2: 35,19 + 2,49 and 34,25 + 3,14, respectively). The ITS1 and ITS2 sequences were more similar to each other inside the groups of species having similar karyotype structure.
In all studied species, the coding sequences of 5S rDNA genes were similar. They were also similar to L. usitatissimum sequences described earlier [26–28]. The coding sequence of flax 5S rRNA genes, in common with the sequences of other plants studied earlier, contained a highly conserved regulatory region known as the Intragenic Control Region that included A and C boxes [29, 30]. In general, the coding sequence of flax 5S rDNA genes common to blue-flowered flax species differed by only four nucleotide replacements from the consensus sequence of angiosperms . In contrary, nontranscribed spacers of 5S rDNA genes were highly variable. Average Evolutionary Divergence for the spacers were 41,52± 2,13 compared to 0,78 ± 0,15 for the coding region. Similar to ITS regions, spacers of 5S rDNA genes were more similar to each other inside the groups of species having similar karyotype structure. It has been established that, in many plants, upstream and downstream spacer regions of coding sequence of 5S rRNA gene were conservative as the regulatory elements important to transcriptional mechanisms were localized there [31, 32].
In the studied flax species, these spacer regions were also conservative. So, for example, for 54 bp of the alignment upstream of the coding region, Average Evolutionary Divergence was 6.75 + 1.12; also, in four position, their nucleotide composition was similar in all group of reads of all the samples (Additional file 3). Besides, all examined samples possessed a poly-C sequence followed by a poly-T structure immediately downstream of the 5S gene coding region.
The study of the phylogenetic relationships within the clade of blue-flowered flaxes is especially important as it includes widely cultivated species L. usitatissimum. Currently, the genome of flax has been sequenced, and continued to be intensively investigated. [33–47]. However, the origin of the genome of cultivated flax is still remained uncleared.
Evolution of karyotypes
The karyotype structure is one of the significant taxonomic characters. It was found that within the section, the flax species have the similar karyotypes (excepting the members of the sect. Linum) but the members of different sections vary considerably in chromosome number and morphology. The differences in chromosome number can be caused by ascending/descending aneuploidy or chromosome fusions. These phenomena are involved into the process of evolution of karyotypes of blue-flowered flaxes. It was believed that genome duplication of the ancient ancestor of the blue-flowered clade took place about 30 million years ago . Perhaps, the polyploid origin of blue-flowered flaxes enabled them to loss chromosomes and to maintain the vitality during genome evolution. The intercalary bands of telomere repeats revealed in karyotypes of several species [17, 48] indicate that chromosome fusions, could also occur during the evolution of flax. Since the event of duplication of the ancestral genome, the significant diploidization of genomes has occurred, and most modern blue-flowered flaxes look like diploid. At the same time, the karyotypes of L. hirsutum ssp. anatolicum, L. hirsutum ssp. pseudoanatolicum (2n = 32), L. perenne ssp. anglicum, L. perenne ssp. extraailare, L. perenne k5500 (2n = 36) and L. narbonense (2n = 28) possess two similar diploid chromosome sets indicating that relatively recently, another duplication of their genomes has taken place. This is not surprising as polyploidization occurring by duplication of chromosomal set is widespread in angiosperms. In plant systematics, the forms having similar morphological traits but differed in ploidy levels are considered to be the cytotypes of one biological species. Changes in ploidy status lead to the whole or partial reproductive isolation, and therefore, in some cases, the appearance of cytotypes is probably the starting point of speciation .
The investigation of the cultivated flax (L. usitatissimum, 2n = 30) showed that this species is also tetraploid. The degree of divergence of the orthologous genes indicated that the origin time of this tetraploid was 5-9 million years [35, 50]. However, unlike the tetraploid species mentioned above, its chromosome set could not be subdivided in two parental genomes though similar features between several pairs of chromosomes were detected . These findings indicate that either multiple chromosomal rearrangements have occurred since the event of duplication of the genome of cultivated flax or the tetraploid flax ancestor could appear as the result of distant hybridization between two diploid species having different karyotypes. Perhaps, L. marginale having a considerable number of chromosomes (2n = 84) is a high-ploidy-level species, but the detailed analysis of its karyotype is impossible since its chromosomes are very small.
FISH analysis showed that most studied species possessed single 5S and 35S rDNA loci in their haploid chromosome sets with the exception of members of the sect. Dasylinum having two 35S rDNA loci and L. narbonense with two 5S rDNA loci in their haploid chromosome sets. Probably, these additional loci appeared de novo during evolution but they also could result from ancient polyploidization. The tetraploid species L. angustifolium and L. usitatissimum and also polyploid L. marginale have one pair of satellite chromosomes in their karyotypes. One NOR is often observed in polyploid genomes of allopolyploids as a consequence of nucleolar dominance, and therefore, L. angustifolium, L. usitatissimum and L. marginale have probably hybrid origins. Besides, karyotypes of L. angustifolium and L. usitatissimum possess three pairs of 5S rDNA loci. It is likely that one of the parental species had two 5S rDNA loci or the additional locus appeared later as a result of dispersion of 5S rDNA throughout the genome. Interestingly, we detected only three pairs of 5S rDNA loci in multipolyploid L. marginale (2n = 84). The other 5S rDNA repeats were probably lost during formation of the polyploid genome, and they could not be revealed by FISH. We previously observed the same phenomenon in some yellow-flowered flaxes of sect. Sillinum .
Variability and evolution of rRNA genes
The reported here ITS1-5.8S rDNA-ITS2 sequences are in agreement with the sequences obtained early by the cloning method . In the studied flaxes similar to many other plants and animals , low levels of intraspecific variability of ITS1-5.8S rDNA-ITS2 sequences revealed, and the character of their variability agreed with the model of concerted evolution. It is believed that high levels of homogeneity of these repeats provides the presence of a quality control system that prevents the accumulation of rDNA mutations [52–54]. The results of high-throughput sequencing for the species of the genus Nicotiana revealed correlation between the number of 35S rDNA loci and variability of ITS1-5.8S rDNA-ITS2 sequences . However, we did not find a significant increase in variability of these sequences within diploid and tetraploid members of the sect. Dasylinum having two NOR loci in a haploid set as well as in tetraploids from the sect. Adenolinum and L. narbonense from the sect. Linum bearing active NOR in each of their four subgenomes compared to the other species having one NOR pair. It is quite possible that this discrepancy between the members of Nicotiana and Linum was due to the fact that multiple NOR loci in Linum could formed by the duplication of a single loci of the ancestor species and therefore, they were initially similar. In the present study, the high-throughput sequencing was used for the first time for analysis of 5S rRNA genes. In flax, the intraspecific variability of the spacers of 5S rRNA genes was considerably higher than the variability of ITS1 and ITS2 sequences of 35S rRNA genes. The highest variability was found in L. angustifolium, L. usitatissimum and L. stelleroides. The heterogeneity of 5S rRNA genes was observed earlier in a lot of other species . Unlike 35S rRNA genes (their variability is usually in good agreement with the model of concerted evolution), the type of variation of 5S rRNA genes often agrees with the model of birth-and-death evolution or with the mixed process of concerted evolution and birth-and-death evolution [51, 57, 58].
It was supposed that the newly generated duplicate genes or gene families may evolve to interact with other existing gene families and promote the adaptation of organisms to new environments .
Despite the question if the heterogeneity of 5S rRNAs is somehow associated with their functions remains unresolved in most cases, there is a great deal of evidence that different variants of 5S rRNA genes could be differentially expressed during development [59–61]. In Arabidopsis thaliana, different variants of 5S rRNA were found in the cells of seeds and roots .
Phylogenetic reconstruction of blue-flowered flaxes
It should be highlighted that the results of phylogenetic analyses based on the 5S rDNA and ITS sequences were analogous with each other. They were in good agreement with the phylogenetic studies reported earlier [8–10, 63] and also with the results of the comparative katyotypic analyses performed in the present work. Particularly, according to our results, the topology of the tree constructed with the use of ITS sequences almost coincided with the topology described earlier . The results of phylogenetic analyses reported here are agreed with the type of variation of some morphological traits revealed among the members of different sections. Particularly, in L. stelleroides, margins of sepals are membranous with large black stipitate glands. The members of the sect. Dasylinum have glandular-ciliate sepals and stipitate glands are small and white. Most species of the sect. Adenolinum and Linum have eglandular sepals with the exception of L. nervosum Waldst. et Kit (not examined in this study) which margins of sepals are slightly serrulate-ciliate with small short stalked glands. Interestingly, for a long time, L. nervosum was considered to be a primitive member of the sect. Linum  but molecular phylogenetic studies  showed that this species represented a basal branch of the clade formed by the members of the sect. Adenolinum. Sepals of the other members of the sect. Adenolinum have entire, scarious margins. The species of the sect. Linum possess sepals with ciliate margins. Molecular phylogenetic studies also indicated that the L. narbonense was the most ancient species among the members of the sect. Linum. This species has sepals with minutely serrulate-ciliate, scarious margins. Сonsidering the fact that secretome of sepal glands attracts pollinating insects, there is probably any relationship between the evolution of sepals in flax and a species composition of the pollinating insects though this problem is not adequately investigated.
Besides, there is a tendency for increase of sizes of seeds and seed capsuls among the sect. Stellerolinum - Dasylinum - Adenolinum - Linum with the exception of L. marginale which has small seeds. It is important to note that the structure of pollen grains in the most ancient species L. stelleroides differs greatly compared to the other species. L. stelleroides has pantoporate pollen grains with 12 pores whereas the members of the other three sections have tricolpate pollen with the exception of few popyploid species having polycolpate pollen . The relationship between the revealed phylogeny and geographical distribution of the species is not quite clear. The centre of origin of genus Linum was believed to locate in the Mediterranean region and western Asia where the most species diversity of wild flaxes of the Old World is observed . Inside this region, the areas of distribution of the separate sections are mostly overlapped. However, the areas of distribution of the evolutionary advanced species, such as members of the sect. Adenolinum and also L. angustifolium and L. usitatissimum from the sect. Linum, are wider and spread farther to the North. The progenitor species of Australian L. marginale and its close relative from New Zealand L. monogynum Rchb. (sect. Linum) were probably arrived from the Northern Hemisphere about 5-1 mya similar to many other members of other families of angiosperms . Interestingly, among the species originated in Northern Hemisphere and migrated to Australia and New Zealand regions, a great number of polyploids including the studied L. marginale and also L. monogynum (2n = 84) were revealed. Unlike the most species of the Old World, the oldest member of blue-flowered flaxes, L. stelleroides, grows in East Asia. Probably, the progenitor of L. stelleroides was distributed across the entire Asia region, but it could survive to nowadays only in the east where where it ‘had waited out’ glaciations in China refugia. Alternatively, the Mediterranean region and western Asia were not the primary but secondary centre of origin of blue-flowered flaxes, and the primary centre was located much farther east than it was thought.
Unlike the previous studies [8–10, 63], we used more species and specimens of the sect. Adenolinum and this allowed us to examine the relationships within this section in greater detail. We detected that the specimens belonging to one and the same species were not always clustered with each other but interspersed with the specimens of other species. Besides, in the present study, the full ITS sequence for L. narbonense (sect. Linum) was obtained, and this allowed us to clarify phylogenetic relationships of this species. It was found that this species was clustered with the other members of the sect. Linum despite the considerable karyotypic differences.
The sequences of 5SrRNA genes were successfully used in phylogenetic analyses of various groups of plants and animals. Along with this, a number of variants of 5S rRNA genes differed in DNA sequences were revealed in genomes of many organisms . In particular, we revealed this phenomenon in L. usitatissimum. High variability of the sequences of 5S rRNA genes in L. usitatissimum was also described early .
Besides, we revealed that in most studied species intraspecific variability of the sequences of 5S rRNA genes was usually less than interspecific variability with the exception of cultivated species L. usitatissimum and its wild ancestor L. angustifolium from the sect. Linum as well as the members of sect. Adenolinum.
According to the obtained ITS phylogeny and the phylogeny based on the mitochondrial rcbL , L. stelleroides was clustered with the members of the sect. Dasyllinum. At the same time, in accordance with the obtained phylogeny of 5S rDNA as well as the phylogeny based on chloroplast (ndhF + trnK intron + trnL-F) and also ITS phylogenies , L. stelleroides formed a basal clade which was independent from the clade of sect. Dasyllinum. These data and also the detected karyological and morphological differences  indicate that L. stelleroides is not a close relative to the species of the sect. Dasyllinum.
As already mentioned, using 5S rDNA or ITS sequences as well as different DNA fingerprinting methods [4, 7, 16] did not allow us to reveal the specific differences within the sect. Adenolinum. Besides, these species were also difficult to distinguish morphologically. The classical taxonomical treatment of this group is based on the quantitative characters which are unstable and depend on the environmental conditions [11, 66]. As up to the present moment the stable and also taxonomically significant characters, which are adequate for speciation within this group, have not been revealed. Therefore, we consider the members of the sect. Adenolinum as a single highly polymorphic species widespread in the North Hemisphere. Interestingly, the same treatment of L. perenne group (syn. sect. Adenolinum) was proposed by Planchon as early as in 1848 (cited according to ). The taxonomic classification of the species L. narbonense from the sect. Linum was controversial as this species differs considerably in chromosome number and sizes, chromosomal distribution of rDNA sites and also due to a low level of similarity with EST-SSRs sequences from the other members of the section, particularly, from L. usitatissimum . However, according to our results, phylogenies of 5S rDNA and ITS as well as the phylogeny based on the mitochondrial rbcL genes  this species is clustered significantly with the other members of the sect. Linum.
Phylogeny of L. usitatissimum and its possible origin
The cultivated flax (L. usitatissimum, 2n = 30) was shown to be an allotetraploid species [35, 50]. Based on both 5S rDNA and ITS phylogenies, we revealed close relationships between L. usitatissimum and L. angustifolium which is considered to be a wild ancestor of the cultivated flax. Both species have identical karyotypes (2n = 30) and similar distribution of 5S and 35S rDNA sites. High degree of relationship between these species is also confirmed by the results of whole genomes and transcriptomes sequencing [9, 10, 37]. Moreover, some taxonomists do not consider them as separate species but subspecies: L. usitatissimum subsp. usitatissimum L. and L. usitatissimum L. subsp. angustifolium (Huds.) Thell. According to morphological characters and molecular phylogenetic investigations, the multipolyploid L. marginale (2n = 84) is relative to tetraploid L. usitatissimum and L. angustifolium. Despite the differences in ploidy levels, these species have similar chromosomal localization of 5S and 45S rRNA genes, and this allowed us to assume that their polyploid genomes could comprise common subgenomes. Unfortunately, it is impossible to verify this assumption as the genome of L. marginale has not been sequenced yet. Polyploid species L. usitatissimum, L.angustifolium and L. marginale are related to the karyologically similar diploids L. grandiflorum (2n = 16) and L. decumbens (2n = 16) and also autotetraploid L. narbonense (2n = 28). All these species are similar in a number of morphological traits (alternate, glabrous leaves, without stipular glands; eglandular, ciliate sepals; blue, purple, pink or white, free petals; linear stigmas) and they are included in the sect. Linum. Close relationship between genomes of L. usitatissimum and L.angustifolium and genomes of L. grandiflorum and L. decumbens is also confirmed by the results of whole genomes and transcriptomes sequencing [35, 50]. These data together with the revealed common structure of satellite chromosomes allowed us to suggest that the ancestor of the modern allotetraploids L. angustifolium and L. usitatissimum could appear as a result of hybridization between two ancient members of the lineage with the basic chromosome number x = 8 (a relative of modern L. grandiflorum and L. decumbens) followed by diploidization of the hybrid genome. In the genome of this newly formed allotetraploid, the reduction of chromosomes from 2n = 32 to 2n = 30 could occur during evolution. Allotetraploid with 2n = 30 could also appear after hybridization between the ancestor of modern L. narbonense with the basic chromosome number x = 7 and the ancestor of modern L. grandiflorum and L. decumbens with the basic chromosome number x = 8; and this hybridization led to the formation of a new chromosome number n = 7 + 8 = 15. An argument in favor of the second hypothesis could be the fact that the number of 5S rDNA sites in karyotypes of L. angustifolium and L. usitatissimum (three pairs) was equal to the total number observed in the karyotypes of L. grandiflorum or L. decumbens (one pair) and in the diploid karyotype of L. narbonense (two pairs). However, the both hypotheses obviously need further confirmations.
Combined molecular and cytogenetic studies and also 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 the intra- and interspecific diversity of the rRNA gene sequences. In particular, the basal position of L. stelleroides within the clade of blue-flowered species and also weak divergence among members of the sect. Adenolinum were confirmed. The specification of the phylogenetic relations among the members of the sect. Linum can contribute to a better understanding of the processes underlying the formation and evolution of the cultivated flax – L. usitatissimum.
This work was performed using the equipment of “Genome” center of Engelhardt Institute of Molecular Biology (http://www.eimb.ru/rus/ckp/ccu_genome_c.php).
This work and the publications costs were financially supported by the Russian Foundation for Basic Research, grants 16-04-01239 and 15-04-06198, and the RAS Presidium Program “Biodiversity of natural systems”, subprogram “Gene pools of nature and their conservation”.
Availability of data and materials
The datasets generated during the current study are available in the Additional files.
About this supplement
This article has been published as part of BMC Evolutionary Biology Volume 17 Supplement 2, 2017: Selected articles from Belyaev Conference 2017: evolutionary biology. The full contents of the supplement are available online at https://bmcevolbiol.biomedcentral.com/articles/supplements/volume-17-supplement-2.
NB, NM, AD, TR, and OM conceived and designed the work; AVK, ASS, AAK, AVS, TR, OY, NB, VL, TS, SZ, and NM performed the experiments; AD, GK, MB, IK, NB, and NM analyzed the data; NB, NM, and АA wrote the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Ockendon DJ, Walters SM. Linaceae, vol. 2. Cambridge: Cambridge University Press; 1968.Google Scholar
- Yuzepchuk SA. Genus Linum. In: Shishkin BK, Ebrov EG, editors. Flora SSSR. Izdatel’stvo Akademii Nauk: Leningrad; 1949.Google Scholar
- Egorova TV. Genus Linum, Linaceae, vol. 9. St Petersburg Mir i Zemia: St Petersburg Publishing House; 1996.Google Scholar
- Vromans J. Molecular genetic studies in flax (Linum usitatissimum L.). Wageningen: Wageningen University; 2006.Google Scholar
- Muravenko OV, Yurkevich OY, Bolsheva NL, Samatadze TE, Nosova IV, Zelenina DA, et al. Comparison of genomes of eight species of sections Linum and Adenolinum from the genus Linum based on chromosome banding, molecular markers and RAPD analysis. Genetica. 2009;135(2):245–55.PubMedView ArticleGoogle Scholar
- Muravenko OV, Bol'sheva NL, Iurkevich O, Nosova IV, Rachinskaia OA, Samatadze TE, et al. Karyogenomics of species of the genus Linum L. Genetika. 2010;46(10):1339–42.PubMedGoogle Scholar
- Melnikova NV, Kudryavtseva AV, Zelenin AV, Lakunina VA, Yurkevich OY, Speranskaya AS, et al. Retrotransposon-based molecular markers for analysis of genetic diversity within the genus Linum. Biomed Res Int. 2014;2014:231589.PubMedPubMed CentralView ArticleGoogle Scholar
- McDill J, Repplinger M, Simpson BB, Kadereit JW. The phylogeny of Linum and Linaceae subfamily linoideae, with implications for their systematics, biogeography, and evolution of heterostyly. Syst Bot. 2009;34(2):386–405.View ArticleGoogle Scholar
- Sveinsson S, McDill J, Wong GK, Li J, Li X, Deyholos MK, et al. Phylogenetic pinpointing of a paleopolyploidy event within the flax genus (Linum) using transcriptomics. Ann Bot. 2014;113(5):753–61.PubMedView ArticleGoogle Scholar
- YB F, Dong Y, Yang MH. Multiplexed shotgun sequencing reveals congruent three-genome phylogenetic signals for four botanical sections of the flax genus Linum. Mol Phylogenet Evol. 2016;101:122–32.View ArticleGoogle Scholar
- Ockendon DJ. Taxonomy of the Linum Perenne group in Europe. Watsonia. 1971;8:205–35.Google Scholar
- Alvarez I, Wendel JF. Ribosomal ITS sequences and plant phylogenetic inference. Mol Phylogenet Evol. 2003;29(3):417–34.PubMedView ArticleGoogle Scholar
- Wicke S, Costa A, Munoz J, Quandt D. Restless 5S: the re-arrangement(s) and evolution of the nuclear ribosomal DNA in land plants. Mol Phylogenet Evol. 2011;61(2):321–32.PubMedView ArticleGoogle Scholar
- Navashin M. Chromosomal alterations caused by hybridization and their bearing upon certain general genetic problems. Cytologia. 1934;5:169–203.View ArticleGoogle Scholar
- Pikaard CS. Nucleolar dominance: uniparental gene silencing on a multi-megabase scale in genetic hybrids. Plant Mol Biol. 2000;43(2-3):163–77.PubMedView ArticleGoogle Scholar
- Yurkevich OY, Naumenko-Svetlova AA, Bolsheva NL, Samatadze TE, Rachinskaya OA, Kudryavtseva AV, et al. Investigation of genome polymorphism and seed coat anatomy of species of section Adenolinum from the genus Linum. Genet Resour Crop Ev. 2013;60(2):661–76.View ArticleGoogle Scholar
- Bolsheva NL, Zelenin AV, Nosova IV, Amosova AV, Samatadze TE, Yurkevich OY, et al. The diversity of karyotypes and genomes within section Syllinum of the genus Linum (Linaceae) revealed by molecular cytogenetic markers and RAPD analysis. PLoS One. 2015;10(4):e0122015.PubMedPubMed CentralView ArticleGoogle Scholar
- Amosova AV, Bolsheva NL, Samatadze TE, Twardovska MO, Zoshchuk SA, Andreev IO, et al. Molecular cytogenetic analysis of Deschampsia Antarctica Desv. (Poaceae), Maritime Antarctic. PLoS One. 2015;10(9):e0138878.PubMedPubMed CentralView ArticleGoogle Scholar
- Gholami M, Bekele WA, Schondelmaier J, Snowdon RJ. A tailed PCR procedure for cost-effective, two-order multiplex sequencing of candidate genes in polyploid plants. Plant Biotechnol J. 2012;10(6):635–45.PubMedView ArticleGoogle Scholar
- Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, et al. Clustal W and Clustal X version 2.0. Bioinformatics. 2007;23(21):2947–8.PubMedView ArticleGoogle Scholar
- Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O. New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol. 2010;59(3):307–21.PubMedView ArticleGoogle Scholar
- Darriba D, Taboada GL, Doallo R, Posada D. jModelTest 2: more models, new heuristics and parallel computing. Nat Methods. 2012;9(8):772.PubMedPubMed CentralView ArticleGoogle Scholar
- Letunic I, Bork P. Interactive tree of life (iTOL) v3: an online tool for the display and annotation of phylogenetic and other trees. Nucleic Acids Res. 2016;44(W1):W242–5.PubMedPubMed CentralView ArticleGoogle Scholar
- Anisimova M, Gil M, Dufayard JF, Dessimoz C, Gascuel O. Survey of branch support methods demonstrates accuracy, power, and robustness of fast likelihood-based approximation schemes. Syst Biol. 2011;60(5):685–99.PubMedPubMed CentralView ArticleGoogle Scholar
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Goldsbrough PB, Ellis TH, Lomonossoff GP. Sequence variation and methylation of the flax 5S RNA genes. Nucleic Acids Res. 1982;10(15):4501–14.PubMedPubMed CentralView ArticleGoogle Scholar
- Schneeberger RG, Creissen GP, Cullis CA. Chromosomal and molecular analysis of 5S RNA gene organization in the flax, Linum Usitatissimum. Gene. 1989;83(1):75–84.PubMedView ArticleGoogle Scholar
- Schneeberger RG, Cullis CA. Intraspecific 5S rRNA gene variation in flax, Linum Usitatissimum (Linaceae). Plant Syst Evol. 1992;183(3/4):265–80.View ArticleGoogle Scholar
- Kellogg EA, Appels R. Intraspecific and interspecific variation in 5S RNA genes are decoupled in diploid wheat relatives. Genetics. 1995;140(1):325–43.PubMedPubMed CentralGoogle Scholar
- Cloix C, Tutois S, Yukawa Y, Mathieu O, Cuvillier C, Espagnol MC, et al. Analysis of the 5S RNA pool in Arabidopsis Thaliana: RNAs are heterogeneous and only two of the genomic 5S loci produce mature 5S RNA. Genome Res. 2002;12(1):132–44.PubMedPubMed CentralView ArticleGoogle Scholar
- Brown G, Carlson J. Molecular cytogenetics of the genes encoding 18s-5.8s-26s rRNA and 5s rRNA in two species of spruce (Picea). Theor Appl Genet. 1997;95(1):1–9.View ArticleGoogle Scholar
- Trontin JF, Grandemange C, Favre JM. Two highly divergent 5S rDNA unit size classes occur in composite tandem array in European larch (Larix Decidua Mill.) and Japanese larch (Larix Kaempferi (Lamb.) Carr.). Genome. 1999;42(5):837–48.PubMedView ArticleGoogle Scholar
- Muir AD, Westcott ND. Flax: the genus Linum. London: Taylor & Francis; 2003.Google Scholar
- Johnson C, Moss T, Cullis C. Environmentally induced heritable changes in flax. J Vis Exp. 2011;47:3223.Google Scholar
- Wang Z, Hobson N, Galindo L, Zhu S, Shi D, McDill J, et al. The genome of flax (Linum Usitatissimum) assembled de novo from short shotgun sequence reads. Plant J. 2012;72(3):461–73.PubMedView ArticleGoogle Scholar
- Dmitriev AA, Kudryavtseva AV, Bolsheva NL, Zyablitsin AV, Rozhmina TA, Kishlyan NV, et al. miR319, miR390, and miR393 are involved in aluminum response in flax (Linum Usitatissimum L.). Biomed Res Int. 2017;2017:4975146.PubMedPubMed CentralGoogle Scholar
- YB F, Peterson GW. Developing genomic resources in two Linum species via 454 pyrosequencing and genomic reduction. Mol Ecol Resour. 2012;12(3):492–500.View ArticleGoogle Scholar
- Melnikova NV, Dmitriev AA, Belenikin MS, Koroban NV, Speranskaya AS, Krinitsina AA, et al. Identification, expression analysis, and target prediction of flax Genotroph MicroRNAs under normal and nutrient stress conditions. Front Plant Sci. 2016;7:399.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhang N, Deyholos MK. RNASeq analysis of the shoot apex of flax (Linum Usitatissimum) to identify phloem fiber specification genes. Front Plant Sci. 2016;7:950.PubMedPubMed CentralGoogle Scholar
- Kumar S, You FM, Cloutier S. Genome wide SNP discovery in flax through next generation sequencing of reduced representation libraries. BMC Genomics. 2012;13:684.PubMedPubMed CentralView ArticleGoogle Scholar
- Dmitriev AA, Kudryavtseva AV, Krasnov GS, Koroban NV, Speranskaya AS, Krinitsina AA, et al. Gene expression profiling of flax (Linum Usitatissimum L.) under edaphic stress. BMC Plant Biol. 2016;16(Suppl 3):237.PubMedView ArticleGoogle Scholar
- Dmitriev AA, Krasnov GS, Rozhmina TA, Kishlyan NV, Zyablitsin AV, Sadritdinova AF, et al. Glutathione S-transferases and UDP-glycosyltransferases are involved in response to aluminum stress in flax. Front Plant Sci. 2016;7:1920.PubMedPubMed CentralView ArticleGoogle Scholar
- Dash PK, Cao Y, Jailani AK, Gupta P, Venglat P, Xiang D, et al. Genome-wide analysis of drought induced gene expression changes in flax (Linum Usitatissimum). GM Crops Food. 2014;5(2):106–19.PubMedPubMed CentralView ArticleGoogle Scholar
- Yu Y, Wu G, Yuan H, Cheng L, Zhao D, Huang W, et al. Identification and characterization of miRNAs and targets in flax (Linum Usitatissimum) under saline, alkaline, and saline-alkaline stresses. BMC Plant Biol. 2016;16(1):124.PubMedPubMed CentralView ArticleGoogle Scholar
- Wojtasik W, Kulma A, Dyminska L, Hanuza J, Czemplik M, Szopa J. Evaluation of the significance of cell wall polymers in flax infected with a pathogenic strain of Fusarium oxysporum. BMC Plant Biol. 2016;16:75.PubMedPubMed CentralView ArticleGoogle Scholar
- Wojtasik W, Kulma A, Namysl K, Preisner M, Szopa J. Polyamine metabolism in flax in response to treatment with pathogenic and non-pathogenic Fusarium strains. Front Plant Sci. 2015;6:291.PubMedPubMed CentralView ArticleGoogle Scholar
- Galindo-Gonzalez L, Pinzon-Latorre D, Bergen EA, Jensen DC, Deyholos MK. Ion torrent sequencing as a tool for mutation discovery in the flax (Linum Usitatissimum L.) genome. Plant Methods. 2015;11:19.PubMedPubMed CentralView ArticleGoogle Scholar
- Bolsheva NL, Semenova OY, Muravenko OV, Nosova IV, Popov KV, Zelenin AV. Localization of telomere sequences in chromosomes of two flax species. Biol Membr. 2005;22(3):227–31.Google Scholar
- Soltis DE, Soltis PS, Schemske DW, Hancock JF, Thompson JN, Husband BC, et al. Autopolyploidy in angiosperms: have we grossly underestimated the number of species? Taxon. 2007;56(1):13–30.Google Scholar
- Cullis CA. DNA sequence organization in the flax genome. Biochim Biophys Acta. 1981;652(1):1–15.PubMedView ArticleGoogle Scholar
- Nei M, Rooney AP. Concerted and birth-and-death evolution of multigene families. Annu Rev Genet. 2005;39:121–52.PubMedPubMed CentralView ArticleGoogle Scholar
- Ganley AR, Kobayashi T. Highly efficient concerted evolution in the ribosomal DNA repeats: total rDNA repeat variation revealed by whole-genome shotgun sequence data. Genome Res. 2007;17(2):184–91.PubMedPubMed CentralView ArticleGoogle Scholar
- Kobayashi T. Regulation of ribosomal RNA gene copy number and its role in modulating genome integrity and evolutionary adaptability in yeast. Cell Mol Life Sci. 2011;68(8):1395–403.PubMedPubMed CentralView ArticleGoogle Scholar
- Tchurikov NA, Fedoseeva DM, Sosin DV, Snezhkina AV, Melnikova NV, Kudryavtseva AV, et al. Hot spots of DNA double-strand breaks and genomic contacts of human rDNA units are involved in epigenetic regulation. J Mol Cell Biol. 2015;7(4):366–82.PubMedView ArticleGoogle Scholar
- Matyasek R, Renny-Byfield S, Fulnecek J, Macas J, Grandbastien MA, Nichols R, et al. Next generation sequencing analysis reveals a relationship between rDNA unit diversity and locus number in Nicotiana diploids. BMC Genomics. 2012;13:722.PubMedPubMed CentralView ArticleGoogle Scholar
- Szymanski M, Zielezinski A, Barciszewski J, Erdmann VA, Karlowski WM. 5SRNAdb: an information resource for 5S ribosomal RNAs. Nucleic Acids Res. 2016;44(D1):D180–3.PubMedView ArticleGoogle Scholar
- Freire R, Arias A, Insua AM, Mendez J, Eirin-Lopez JM. Evolutionary dynamics of the 5S rDNA gene family in the mussel Mytilus: mixed effects of birth-and-death and concerted evolution. J Mol Evol. 2010;70(5):413–26.PubMedView ArticleGoogle Scholar
- Eirin-Lopez JM, Rebordinos L, Rooney AP, Rozas J. The birth-and-death evolution of multigene families revisited. Genome Dyn. 2012;7:170–96.PubMedView ArticleGoogle Scholar
- Komiya H, Hasegawa M, Takemura S. Differentiation of oocyte- and somatic-type 5S rRNAs in animals. J Biochem. 1986;100(2):369–74.PubMedView ArticleGoogle Scholar
- Dimarco E, Cascone E, Bellavia D, Caradonna F. Functional variants of 5S rRNA in the ribosomes of common sea urchin Paracentrotus Lividus. Gene. 2012;508(1):21–5.PubMedView ArticleGoogle Scholar
- Hariharan N, Reddy PS, Padayatty JD. 5S-rRNA genes in rice embryos. Plant Mol Biol. 1987;9(5):443–51.PubMedView ArticleGoogle Scholar
- Mathieu O, Jasencakova Z, Vaillant I, Gendrel AV, Colot V, Schubert I, et al. Changes in 5S rDNA chromatin organization and transcription during heterochromatin establishment in Arabidopsis. Plant Cell. 2003;15(12):2929–39.PubMedPubMed CentralView ArticleGoogle Scholar
- McDill J, Simpson BB. Molecular phylogenetics of Linaceae with complete generic sampling and data from two plastid genes. Bot J Linn Soc. 2011;165:64–83.View ArticleGoogle Scholar
- Rogers CM, Xavier KS. Pollen morphology as an aid in determining relationships among some widely separated old world species of Linum. Grana. 1971;11(1):55–7.View ArticleGoogle Scholar
- Winkworth RC, Wagstaff SJ, Glenny D, Lockhart PJ. Evolution of the New Zealand mountain flora: origins, diversification and dispersal. Org Divers Evol. 2005;5(3):237–47.View ArticleGoogle Scholar
- Ockendon DJ. Biosystematic studies in the Linum Perenne group. New Phytol. 1968;67(4):787–813.View ArticleGoogle Scholar
- Soto-Cerda BJ, Saavedra HU, Navarro CN, Ortega PM. Characterization of novel genic SSR markers in Linum Usitatissimum (L.) and their transferability across eleven Linum species. Electron J Biotechnol. 2011;14(2):4.View ArticleGoogle Scholar