Evolution of blue-flowered species of genus Linum based on high-throughput sequencing of ribosomal RNA genes

Drivers of diversification in Linum (Linaceae) by means of chromosome evolution: correlations with biogeography, breeding system and habit

BACKGROUND AND AIMS

Chromosome evolution leads to hybrid dysfunction and recombination patterns and has thus been proposed as a major driver of diversification in all branches of the tree of life, including flowering plants. In this study we used the genus Linum (flax species) to evaluate the effects of chromosomal evolution on diversification rates and on traits that are important for sexual reproduction. Linum is a useful study group because it has considerable reproductive polymorphism (heterostyly) and chromosomal variation (n = 6-36) and a complex pattern of biogeographical distribution.

METHODS

We tested several traditional hypotheses of chromosomal evolution. We analysed changes in chromosome number across the phylogenetic tree (ChromEvol model) in combination with diversification rates (ChromoSSE model), biogeographical distribution, heterostyly and habit (ChromePlus model).

KEY RESULTS

Chromosome number evolved across the Linum phylogeny from an estimated ancestral chromosome number of n = 9. While there were few apparent incidences of cladogenesis through chromosome evolution, we inferred up to five chromosomal speciation events. Chromosome evolution was not related to heterostyly but did show significant relationships with habit and geographical range. Polyploidy was negatively correlated with perennial habit, as expected from the relative commonness of perennial woodiness and absence of perennial clonality in the genus. The colonization of new areas was linked to genome rearrangements (polyploidy and dysploidy), which could be associated with speciation events during the colonization process.

CONCLUSIONS

Chromosome evolution is a key trait in some clades of the Linum phylogeny. Chromosome evolution directly impacts speciation and indirectly influences biogeographical processes and important plant traits.

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

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

Discrimination of Green Coffee (Coffea arabica and Coffea canephora) of Different Geographical Origin Based on Antioxidant Activity, High-Throughput Metabolomics, and DNA RFLP Fingerprinting

The genus Coffea is known for the two species C. arabica (CA) and C. canephora (CC), which are used to prepare the beverage coffee. Proper identification of green beans of coffee varieties is based on phenotypic and phytochemical/molecular characteristics. In this work, a combination of chemical (UV/Vis, HPLC-DAD-MS/MS, GC-MS, and GC-FID) and molecular (PCR-RFLP) fingerprinting was used to discriminate commercial green coffee accessions from different geographical origin. The highest content of polyphenols and flavonoids was always found in CC accessions, whereas CA showed lower values. ABTS and FRAP assays showed a significant correlation between phenolic content and antioxidant activity in most CC accessions. We identified 32 different compounds, including 28 flavonoids and four N-containing compounds. The highest contents of caffeine and melatonin were detected in CC accessions, whereas the highest levels of quercetin and kaempferol derivatives were found in CA accessions. Fatty acids of CC accessions were characterized by low levels of linoleic and cis octadecenoic acid and high amounts of elaidic acid and myristic acid. Discrimination of species according to their geographical origin was achieved using high-throughput data analysis, combining all measured parameters. Lastly, PCR-RFLP analysis was instrumental for the identification of recognition markers for the majority of accessions. Using the restriction enzyme AluI on the trnL-trnF region, we clearly discriminated C. canephora from C. arabica, whereas the cleavage performed by the restriction enzymes MseI and XholI on the 5S-rRNA-NTS region produced specific discrimination patterns useful for the correct identification of the different coffee accessions. This work extends our previous studies and provides new information on the complete flavonoid profile, combining high-throughput data with DNA fingerprinting to assess the geographical discrimination of green coffee.

Development and Complex Application of Methods for the Identification of Mutations in the FAD3A and FAD3B Genes Resulting in the Reduced Content of Linolenic Acid in Flax Oil

Flax is grown worldwide for seed and fiber production. Linseed varieties differ in their oil composition and are used in pharmaceutical, food, feed, and industrial production. The field of application primarily depends on the content of linolenic (LIN) and linoleic (LIO) fatty acids. Inactivating mutations in the FAD3A and FAD3B genes lead to a decrease in the LIN content and an increase in the LIO content. For the identification of the three most common low-LIN mutations in flax varieties (G-to-A in exon 1 of FAD3A substituting tryptophan with a stop codon, C-to-T in exon 5 of FAD3A leading to arginine to a stop codon substitution, and C-to-T in exon 2 of FAD3B resulting in histidine to tyrosine substitution), three approaches were proposed: (1) targeted deep sequencing, (2) high resolution melting (HRM) analysis, (3) cleaved amplified polymorphic sequences (CAPS) markers. They were tested on more than a thousand flax samples of various types and showed promising results. The proposed approaches can be used in marker-assisted selection to choose parent pairs for crosses, separate heterogeneous varieties into biotypes, and select genotypes with desired homozygous alleles of the FAD3A and FAD3B genes at the early stages of breeding for the effective development of varieties with a particular LIN and LIO content, as well as in basic studies of the molecular mechanisms of fatty acid synthesis in flax seeds to select genotypes adequate to the tasks.

A legume kinesin controls vacuole morphogenesis for rhizobia endosymbiosis

Symbioses between legumes and rhizobia require establishment of the plant-derived symbiosome membrane, which surrounds the rhizobia and accommodates the symbionts by providing an interface for nutrient and signal exchange. The host cytoskeleton and endomembrane trafficking systems play central roles in the formation of a functional symbiotic interface for rhizobia endosymbiosis; however, the underlying mechanisms remain largely unknown. Here we demonstrate that the nodulation-specific kinesin-like calmodulin-binding protein (nKCBP), a plant-specific microtubule-based kinesin motor, controls central vacuole morphogenesis in symbiotic cells in Medicago truncatula. Phylogenetic analysis further indicated that nKCBP duplication occurs solely in legumes of the clade that form symbiosomes. Knockout of nKCBP results in central vacuole deficiency, defective symbiosomes and abolished nitrogen fixation. nKCBP decorates linear particles along microtubules, and crosslinks microtubules with the actin cytoskeleton, to control central vacuole formation by modulating vacuolar vesicle fusion in symbiotic cells. Together, our findings reveal that rhizobia co-opted nKCBP to achieve symbiotic interface formation by regulating cytoskeletal assembly and central vacuole morphogenesis during nodule development.

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

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

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.

High-throughput sequencing techniques to flax genetics and breeding

Flax (Linum usitatissimum L.) is an important oil and fiber crop. Using modern methods for flax breeding allows accelerating the introduction of some desired genes into the genotypes of future varieties. Today, an important condition for their creation is the development of research, that is based on next-generation sequencing (NGS). This review summarizes the results obtained using NGS in flax research. To date, a linkage map with a high marker density has been obtained for L. usitatissimum, which is already being used for a more efficient search for quantitative traits loci. Comparative studies of transcriptomes and miRNomes of flax under stress and in control conditions elucidated molecular-genetic mechanisms of abiotic and biotic stress responses. The very accurate model for genomic selection of flax resistant to pasmo was constructed. Based on NGS-sequencing also some details of the genus Linum evolution were clarified. The knowledge systematized in the review can be useful for researchers working in flax breeding and whereas fundamental interest for understanding the phylogenetic relationships within the genus Linum, the ontogenesis, and the mechanisms of the response of flax plants to various stress factors.

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.