Locus-specific view of flax domestication history

Crop domestication has been inferred genetically from neutral markers and increasingly from specific domestication-associated loci. However, some crops are utilized for multiple purposes that may or may not be reflected in a single domestication-associated locus. One such example is cultivated flax (Linum usitatissimum L.), the earliest oil and fiber crop, for which domestication history remains poorly understood. Oil composition of cultivated flax and pale flax (L. bienne Mill.) indicates that the sad2 locus is a candidate domestication locus associated with increased unsaturated fatty acid production in cultivated flax. A phylogenetic analysis of the sad2 locus in 43 pale and 70 cultivated flax accessions established a complex domestication history for flax that has not been observed previously. The analysis supports an early, independent domestication of a primitive flax lineage, in which the loss of seed dispersal through capsular indehiscence was not established, but increased oil content was likely occurred. A subsequent flax domestication process occurred that probably involved multiple domestications and includes lineages that contain oil, fiber, and winter varieties. In agreement with previous studies, oil rather than fiber varieties occupy basal phylogenetic positions. The data support multiple paths of flax domestication for oil-associated traits before selection of the other domestication-associated traits of seed dispersal loss and fiber production. The sad2 locus is less revealing about the origin of winter tolerance. In this case, a single domestication-associated locus is informative about the history of domesticated forms with the associated trait while partially informative on forms less associated with the trait.     

Genome-wide analysis of genomic imprinting in the endosperm and allelic variation in flax

Genomic imprinting is an epigenetic phenomenon that causes biased expression of maternally and paternally inherited alleles. In flowering plants, genomic imprinting predominantly occurs in the triploid endosperm and plays a vital role in seed development. In this study, we identified 248 candidate imprinted genes including 114 maternally expressed imprinted genes (MEGs) and 134 paternally expressed imprinted genes (PEGs) in flax (Linum usitatissimum L.) endosperm using deep RNA sequencing. These imprinted genes were neither clustered in specific chromosomal regions nor well conserved among flax and other plant species. MEGs tended to be expressed specifically in the endosperm, whereas the expression of PEGs was not tissue-specific. Imprinted single nucleotide polymorphisms differentiated 200 flax cultivars into the oil flax, oil-fiber dual purpose flax and fiber flax subgroups, suggesting that genomic imprinting contributed to intraspecific variation in flax. The nucleotide diversity of imprinted genes in the oil flax subgroup was significantly higher than that in the fiber flax subgroup, indicating that some imprinted genes underwent positive selection during flax domestication from oil flax to fiber flax. Moreover, imprinted genes that underwent positive selection were related to flax functions. Thirteen imprinted genes related to flax seed size and weight were identified using a candidate gene-based association study. Therefore, our study provides information for further exploration of the function and genomic variation of imprinted genes in the flax population.     

Evidence of the domestication history of flax (Linum usitatissimum L.) from genetic diversity of the sad2 locus

A phylogenetic analysis was conducted on 34 alleles of 2.5 kb sized stearoyl-ACP desaturase II (sad2), obtained from 30 accessions of cultivated and pale flax (Linum spp.), to elucidate the history of flax domestication. The analysis supports a single domestication origin for extant cultivated flax. The phylogenetic evidence indicates that flax was first domesticated for oil, rather than fibre. The genetic diversity of the sad2 locus in cultivated flax is low when compared to that of the pale flax assayed. An absolute archaeological date could be applied to the synonymous substitution rate of sad2 in cultivated flax, yielding a high estimate of 1.60–1.71×10⁻⁷ substitutions/site/year. The occurrence of nonsynonymous substitutions at conserved positions of the third exon in alleles from cultivated flax suggests that the locus may have been subjected to an artificial selection pressure. The elevated synonymous substitution rate is also compatible with a population expansion of flax since domestication, followed by a population decline in historic times. These findings provide new insight into flax domestication and are significant for the continuous exploration of the flax germplasm for utilization.     

Genetic polymorphism of flax Linum usitatissimum based on the use of molecular cytogenetic markers

Using a set of approaches based on the use of molecular cytogenetic markers (DAPI/C-banding, estimation of the total area of DAPI-positive regions in prophase nuclei, FISH with 26S and 5S rDNA probes) and the microsatellite (SSR-PCR) assay, we studied genomic polymorphism in 15 flax (Linum usitatissimum L.) varieties from different geographic regions belonging to three directions of selection (oil, fiber, and intermediate flax) and in the k-37 × Viking hybrid. All individual chromosomes have been identified in the karyotypes of these varieties on the basis of the patterns of differential DAPI/C-banding and the distribution of 26S and 5S rDNA, and idiograms of the chromosomes have been generated. Unlike the oil flax varieties, the chromosomes in the karyotypes of the fiber flax varieties have, as a rule, pericentromeric and telomeric DAPI-positive bands of smaller size, but contain larger intercalary regions. Two chromosome rearrangements (chromosome 3 inversions) were detected in the variety Luna and in the k-37 × Viking hybrid. In both these forms, no colocalization of 26S rDNA and 5S rDNA on the satellite chromosome was detected. The SSR assay with the use of 20 polymorphic pairs of primers revealed 22 polymorphic loci. Based on the SSR data, we analyzed genetic similarity of the flax forms studied and constructed a genetic similarity dendrogram. The genotypes studied here form three clusters. The oil varieties comprise an independent cluster. The genetically related fiber flax varieties Vita and Luna, as well as the landrace Lipinska XIII belonging to the intermediate type, proved to be closer to the oil varieties than the remaining fiber flax varieties. The results of the molecular chromosome analysis in the fiber and oil flax confirm their very close genetic similarity. In spite of this, the combined use of the chromosome and molecular markers has opened up unique possibilities for describing the genotypes of flax varieties and creating their genetic passports.     

Genetic polymorphism of flax Linum usitatissimum based on use of molecular cytogenetic markers

Using a set of approaches based on the use of molecular cytogenetic markers (DAPI/C-banding, estimation of the total area of DAPI-positive regions in prophase nuclei, FISH with 26S and 5S rDNA probes) and the microsatellite (SSR-PCR) assay, we studied genomic polymorphism in 15 flax (Linum usitatissimum L.) varieties from different geographic regions belonging to three directions of selection (oil, fiber, and intermediate flaxes) and in the k-37 x Viking hybrid. All individual chromosomes have been identified in the karyotypes of these varieties on the basis of the patterns of differential DAPI/C-banding and the distribution of 26S and 5S rDNA, and idiograms of the chromosomes have been generated. Unlike the oil flax varieties, the chromosomes in the karyotypes of the fiber flax varieties have, as a rule, pericentromeric and telomeric DAPI-positive bands of smaller size, but contain larger intercalary regions. Two chromosomal rearrangements (chromosome 3 inversions) were discovered in the variety Luna and in the k-37 x Viking hybrid. In both these forms, no colocalization of 26S rDNA and 5S rDNA on the satellite chromosome was detected. The SSR assay with the use of 20 polymorphic pairs of primers revealed 22 polymorphic loci. Based on the SSR data, we analyzed genetic similarity of the flax forms studied and constructed a genetic similarity dendrogram. The genotypes studied here form three clusters. The oil varieties comprise an independent cluster. The genetically related fiber flax varieties Vita and Luna, as well as the landrace Lipinska XIII belonging to the intermediate type, proved to be closer to the oil varieties than the remaining fiber flax varieties. The results of the molecular chromosomal analysis in the fiber and oil flaxes confirm their very close genetic similarity. In spite of this, the combined use of the chromosomal and molecular markers has opened up unique possibilities for describing the genotypes of flax varieties and creating their genetic passports.     

Developing genomic resources in two Linum species via 454 pyrosequencing and genomic reduction

Recent advances in next-generation DNA sequencing (NGS) have enhanced the development of genomic resources such as contigs or single-nucleotide polymorphisms (SNPs) for evolutionary studies of a nonmodel species with a complex and unsequenced genome. This study presents an application of a NGS technique in combination with genomic reduction and advanced bioinformatics tools to identify contigs and SNPs from multiple samples of two Linum species. A full Roche 454 GS FLX run of 16 diverse Linum samples representing cultivated flax (Linum usitatissimum L.) and its wild progenitor (Linum bienne Mill.) generated approximately 1.6 million sequence reads with a total length of 498 Mbp. Application of the computational pipeline de novo identification of alleles identified 713 contigs and 1067 SNPs. A blast search revealed alignments of all 713 contigs with 491 existing Linum scaffolds and gene annotations associated with 512 contigs. Sanger sequencing confirmed 95% of 79 selected contigs and 94% of 272 SNPs and identified 211 new SNPs and 19 new indels. The scored 454 SNP data were highly imbalanced for assayed samples. These findings not only are useful for evolutionary studies of Linum species but also help to illustrate the utility of NGS technologies in SNP discovery for nonmodel organisms.     

Historical and botanical evidence of distribution, cultivation and utilization of Linum usitatissimum L. (flax) in China

Linum usitatissimum L. (flax) has been cultivated for more than 2,000 years in China. Throughout the history of its cultivation, flax was mainly used for oil and medicine. It was not until the last 100 years that fibre flax began to be cultivated and used. This paper discusses the distribution, cultivation and use of flax in China especially in the past, based on literature reviews and ethnobotanical surveys for its traditional uses. The modern distribution of wild flax species in China is also presented in this paper.     

Polymorphism of heterochromatic regions of flax chromosomes

The C-banding technique was used to study flax chromosomes (Linum usitatissimum L., 2n = 30). Heterochromatin was located mainly in pericentromeric regions of chromosomes. In spite of small size (1.5-3.5 microm), all 15 pairs of homologous chromosomes were identified on the basis of the C-banding pattern and morphology. An idiogram of C-banded chromosomes of L usitatissimum L. is presented. Polymorphism of chromosomal heterochromatic regions was studied in karyotypes of three flax samples: L usitatissimum L., accession K-603 (L usitatissimum var. usitatissimum), and accession K-594 (L. usitatissimum var. humile (Mill.)). A common C-banding pattern was observed in all forms studied, although there were some distinctions in the individual band size. The fibre flax (accession K-603) karyotype had the C-banding pattern similar to that of L usitatissimum L., but some intercalary and telomeric C-bands were somewhat larger, and a satellite (NOR) was observed in the short arm of chromosome I. In crown flax, (K-594) chromosomal C-banding pattern exhibited smaller pericentromeric and larger intercalary bands; telomeric bands were present on almost all chromosomes. Thus, the intraspecies polymorphism revealed in the chromosomal C-banding pattern makes possible the use of C-bands as chromosome markers in the studies of genetic and genomic polymorphism of this species.     

Total area of C-heterochromatin in nuclei and chromosomes as a measure of cultivated flax genome variability

Comparative statistical analysis of total C-heterochromatin content (total area of darkly stained C-bands) in prophase nuclei and metaphase chromosomes of eight flax (Linum usitatissimum L.) varieties, viz., three oil flax varieties (Shafir, Oliver, Linotte), three modern (Orshansky-2, Belinka, Baltuchai) and two ancient fiber flax varieties (Zaretsky Kryazh, K-37) has been carried out. Karyotypes of flax varieties selected for different agronomical traits differed from one another in C-heterochromatin content. The total C-band area of oil flax varieties exceeded that of fiber flax varieties. Genomic polymorphism in C-heterochromatin markers is a testimony of different breeding directions. A simple test for determining total C-heterochromatin content in prophase nuclei is recommended as a useful tool for screening specimens with predetermined characteristics during selection of new varieties.     

Osmotic stress alters the balance between organic and inorganic solutes in flax (Linum usitatissimum)

Flax (Linum usitatissimum) is grown for its oil and its fiber. This crop, cultivated in temperate regions, has seen a renewed interest due to the presence of abundant molecules of interest for many applications. Little information is available about the behavior of flax during osmotic stress; yet this is considered a major stress that causes significant yield losses in most crops. To control the presence of this stress better, flax behavior was investigated following the application of osmotic stress and the response was examined by applying increasing concentrations of PEG 8000. This resulted in the reorganization of 32 metabolites and 6 mineral ions in the leaves. The analysis of these two types of solute highlighted the contrasting behavior between a higher metabolite content (particularly fructose, glucose and proline) and a decrease in mineral ions (especially nitrate and potassium) following PEG treatment. However, this reorganization did not lead to a greater accumulation of solutes, with the total amount remaining unchanged in leaves during osmotic stress.     

Genetic diversity of cultivated flax (Linum usitatissimum L.) germplasm assessed by retrotransposon-based markers

Retrotransposon segments were characterized and inter-retrotransposon amplified polymorphism (IRAP) markers developed for cultivated flax (Linum usitatissimum L.) and the Linum genus. Over 75 distinct long terminal repeat retrotransposon segments were cloned, the first set for Linum, and specific primers designed for them. IRAP was then used to evaluate genetic diversity among 708 accessions of cultivated flax comprising 143 landraces, 387 varieties, and 178 breeding lines. These included both traditional and modern, oil (86), fiber (351), and combined-use (271) accessions, originating from 36 countries, and 10 wild Linum species. The set of 10 most polymorphic primers yielded 141 reproducible informative data points per accession, with 52% polymorphism and a 0.34 Shannon diversity index. The maximal genetic diversity was detected among wild Linum species (100% IRAP polymorphism and 0.57 Jaccard similarity), while diversity within cultivated germplasm decreased from landraces (58%, 0.63) to breeding lines (48%, 0.85) and cultivars (50%, 0.81). Application of Bayesian methods for clustering resulted in the robust identification of 20 clusters of accessions, which were unstratified according to origin or user type. This indicates an overlap in genetic diversity despite disruptive selection for fiber versus oil types. Nevertheless, eight clusters contained high proportions (70?C100%) of commercial cultivars, whereas two clusters were rich (60%) in landraces. These findings provide a basis for better flax germplasm management, core collection establishment, and exploration of diversity in breeding, as well as for exploration of the role of retrotransposons in flax genome dynamics.     

Mechanisms and control of rapid genomic changes in flax

BACKGROUND AND AIMS: The nuclear DNA of certain varieties of flax (Linum usitatissimum) can vary within a single generation when the plants are grown under specific environmental conditions. This review details the genomic variations that have been identified and associated with this environmental response. CONCLUSIONS: The variation occurs across the whole spectrum of sequence repetition and has been shown to occur in the highly repeated, middle repetitive and low copy number sequences. Although the variation has been shown to be spread throughout the genome it does not occur at random, as similar molecular events have been shown to occur repeatedly. The changes in two labile regions in the nucleus, the ribosomal RNA genes and a site-specific insertion event, have been shown to occur within the period of vegetative growth and over a relatively short period of that growth. The gradual change in total nuclear DNA that has been described would then need to have arisen through an accumulation of changes occurring over the whole, or most of the, period of growth prior to flowering. The polymorphisms that result from these rapidly occurring genomic events have also been observed in many other flax and linseed varieties as well as in the wild progenitors of flax.     

Variability in the composition of tissue-specific galactan from flax fibers

Tissue-specific galactan of sclerenchyma fibers, with cell walls of the gelatinous type, was examined in flax plants (Linum usitatissimum L.) of 23 various genotypes. The content and average degree of polymerization of side chains of galactan were estimated before its deposition into the cell wall. The variability of the analyzed parameters of tissue-specific galactan from flax fibers was high; within the same genotype, the scope of paratypic variability between replicates and years of research was comparable to variability between different genotypes. The average length of side chains in the studied samples ranged from 5 to 41 galactose residues. The average degrees of polymerization of galactan side chains in flax fibers was found to be discrete, which could be explained by block assemblage of the polymer in the Golgi apparatus.     

RAPD analysis of the flax (Linum usitatissimum L.) varieties and hybrids of various productivity

Genetic polymorphism in varieties and hybrids of cultivated flax (Linum usitatissimum L.) has been investigated by RAPD-PCR. Analysis with 15 primers has revealed varietal specificity and hybrid inheritance of RAPD alleles. This allows genetic certification of the original varieties and their hybrids for breeding purposes. Polymorphic amplification products were obtained in RAPD analysis of DNA from two cultivated flax varieties with the use of 10-11 nucleotide primers.     

The genome of flax (Linum usitatissimum) assembled de novo from short shotgun sequence reads

Flax (Linum usitatissimum) is an ancient crop that is widely cultivated as a source of fiber, oil and medicinally relevant compounds. To accelerate crop improvement, we performed whole‐genome shotgun sequencing of the nuclear genome of flax. Seven paired‐end libraries ranging in size from 300 bp to 10 kb were sequenced using an Illumina genome analyzer. A de novo assembly, comprised exclusively of deep‐coverage (approximately 94× raw, approximately 69× filtered) short‐sequence reads (44–100 bp), produced a set of scaffolds with N₅₀ = 694 kb, including contigs with N₅₀ = 20.1 kb. The contig assembly contained 302 Mb of non‐redundant sequence representing an estimated 81% genome coverage. Up to 96% of published flax ESTs aligned to the whole‐genome shotgun scaffolds. However, comparisons with independently sequenced BACs and fosmids showed some mis‐assembly of regions at the genome scale. A total of 43 384 protein‐coding genes were predicted in the whole‐genome shotgun assembly, and up to 93% of published flax ESTs, and 86% of A. thaliana genes aligned to these predicted genes, indicating excellent coverage and accuracy at the gene level. Analysis of the synonymous substitution rates (K_s) observed within duplicate gene pairs was consistent with a recent (5–9 MYA) whole‐genome duplication in flax. Within the predicted proteome, we observed enrichment of many conserved domains (Pfam‐A) that may contribute to the unique properties of this crop, including agglutinin proteins. Together these results show that de novo assembly, based solely on whole‐genome shotgun short‐sequence reads, is an efficient means of obtaining nearly complete genome sequence information for some plant species.     

Comparative study of genomes of two species of flax by C-banding of chromosomes

C-banding patterns of the karyotypes of two closely related wild flax species, Linum austriacum L. (2n = 18) and Linum grandiflorum Desf. (2n = 16), were studied. The karyotypes of both species were similar in the chromosome morphology and size. In each species, metacentric and acrocentric chromosomes (1.7-4.3 microns) and one satellite chromosome were observed. In the karyotypes of the species studied, all homologous chromosome pairs were identified, and quantitative ideograms were constructed. Eight chromosome pairs in the two species had similar C-banding patterns. A low level of intraspecific polymorphism in the intercalary and telomeric C-bands was shown in both species. The results indicate that the genomes of two flax species originated from one ancestral genome with the main chromosome number of 8 or 9. Apparently, the doubling of chromosome number or loss of one chromosome with subsequent redistribution of the chromosome material in the ancestral form resulted in the divergence into two species, L. austriacum L. and L. grandiflorum Desf. A considerable similarity of chromosomes in these species provides evidence for their close phylogenetic relatedness, which makes it possible to place them in one section within the Linum genus.     

The effect of soil drought on the phloem fiber development in long-fiber flax

The effects of soil drought on various stages of phloem fiber development during the period of flax (Linum usitatissimum L.) rapid growth were assessed. The formation of the secondary cell wall was shown to be most retarded. The content of a tissue-specific galactan was reduced especially sharply, and its side chains were changed. Under conditions of pronounced stress-induced plant growth retardation, fiber intrusive growth was suppressed relatively softly: their number on the stem transverse sections was reduced only by 16%. However, this determined irreversible diversity in the fiber length in various stem regions. Such insignificant suppression of intrusive growth under osmotic stress (simultaneously with substantial retardation of plant growth and metabolism inhibition) indicates the functioning of special mechanisms of its regulation.     

Structural organization and classification of cytochrome P450 genes in flax (Linum usitatissimum L.)

Flax CYPome analysis resulted in the identification of 334 putative cytochrome P450 (CYP450) genes in the cultivated flax genome. Classification of flax CYP450 genes based on the sequence similarity with Arabidopsis orthologs and CYP450 nomenclature, revealed 10 clans representing 44 families and 98 subfamilies. CYP80, CYP83, CYP92, CYP702, CYP705, CYP708, CYP728, CYP729, CYP733 and CYP736 families are absent in the flax genome. The subfamily members exhibited conserved sequences, length of exons and phasing of introns. Similarity search of the genomic resources of wild flax species Linum bienne with CYP450 coding sequences of the cultivated flax, revealed the presence of 127 CYP450 gene orthologs, indicating amplification of novel CYP450 genes in the cultivated flax. Seven families CYP73, 74, 75, 76, 77, 84 and 709, coding for enzymes associated with phenylpropanoid/fatty acid metabolism, showed extensive gene amplification in the flax. About 59% of the flax CYP450 genes were present in the EST libraries.     

Environmental factors affecting the cyanogenic potential of flax seedlings

The levels of cyanogenic glucosides (linamarin and lotaustralin) and the activity of linamarase were studied in 5-day old seedlings of oil flax (Linum usitatissimum L., cv. LCSD 200) under different environmental conditions. White light enhanced the cyanoglucosides content, and this effect depended on its intensity and the time of exposure. The level of cyanoglucosides rose with temperature, and it reached the highest level at the highest temperature (30 °C). Linamarase (EC. 3.2.1.21) activity was the highest at 20°C, especially in light-grown seedlings. Lower enzyme activity at the extreme temperature (15 and 30 °C) was observed. Water stress (low water potential, ω=−0.34 MPa) reduced by more than twice the cyanoglucoside level and linamarase activity. The possible protective, or/and regulatory roles of cyanogenic glucosides was discussed.     

Botanical analysis of a bundle of flax (Linum usitatissimum L.) from an early medieval site in northern Poland; a contribution to the history of flax cultivation and its field weeds

A bundle of flax (Linum usitatissimum L.) radiocarbon dated to 1210±70 uncal B.P. (830±90 cal A.D.) was analysed for its macrofossil content. Apart from stems, capsules and seeds of flax., a large number of diaspores (fruits and seeds) from other plants was identified. Field weeds were the most numerous taxa present, among them three flax field weeds,Spergula maxima, Camelina alyssum andCuscuta cpilinum. Development of the specific flax weed community is discussed. Indicator values are used to characterize the edaphic conditions of this early medieval flax field. The field weeds spectrum also suggests that this flax was sown as a summer crop after an earlier crop of millet.