RNASeq Analysis of the Shoot Apex of Flax (Linum usitatissimum) to Identify Phloem Fiber Specification Genes

Identification of new potential molecular actors related to fiber quality in flax through Omics

One of the biggest challenges for a more widespread utilization of plant fibers is to better understand the different molecular factors underlying the variability in fineness and mechanical properties of both elementary and scutched fibers. Accordingly, we analyzed genome-wide transcription profiling from bast fiber bearing tissues of seven different flax varieties (4 spring, 2 winter fiber varieties and 1 winter linseed) and identified 1041 differentially expressed genes between varieties, of which 97 were related to cell wall metabolism. KEGG analysis highlighted a number of different enriched pathways. Subsequent statistical analysis using Partial Least-Squares Discriminant Analysis showed that 73% of the total variance was explained by the first 3 X-variates corresponding to 56 differentially expressed genes. Calculation of Pearson correlations identified 5 genes showing a strong correlation between expression and morphometric data. Two-dimensional gel proteomic analysis on the two varieties showing the most discriminant and significant differences in morphometrics revealed 1490 protein spots of which 108 showed significant differential abundance. Mass spectrometry analysis successfully identified 46 proteins representing 32 non-redundant proteins. Statistical clusterization based on the expression level of genes corresponding to the 32 proteins showed clear discrimination into three separate clusters, reflecting the variety type (spring-/winter-fiber/oil). Four of the 32 proteins were also highly correlated with morphometric features. Examination of predicted functions for the 9 (5 + 4) identified genes highlighted lipid metabolism and senescence process. Calculation of Pearson correlation coefficients between expression data and retted fiber mechanical measurements (strength and maximum force) identified 3 significantly correlated genes. The genes were predicted to be connected to cell wall dynamics, either directly (Expansin-like protein), or indirectly (NAD(P)-binding Rossmann-fold superfamily protein). Taken together, our results have allowed the identification of molecular actors potentially associated with the determination of both in-planta fiber morphometrics, as well as ex-planta fiber mechanical properties, both of which are key parameters for elementary fiber and scutched fiber quality in flax.

Spatiotemporal Transcriptomic Atlas of Developing Embryos and Vegetative Tissues in Flax

Flax (Linum usitatissimum L.) is an important multipurpose crop widely grown for oil and fiber. Despite recent advances in genomics, detailed gene activities during the important reproductive phase of its development are not well defined. In this study, we employed high-throughput RNA-sequencing methods to generate in-depth transcriptome profiles of flax tissues with emphasis on the reproductive phases of five key stages of embryogenesis (globular embryo, heart embryo, torpedo embryo, cotyledon embryo, and mature embryo), mature seed, and vegetative tissues viz. ovary, anther, and root. These datasets were used to establish the co-expression networks covering 36 gene modules based on the expression patterns for each gene through weighted gene co-expression network analysis (WGCNA). Functional interrogation with Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) of dominantly expressed genetic modules in tissues revealed pathways involved in the development of different tissues. Moreover, the essential genes in embryo development and synthesis of storage reserves were identified based on their dynamic expression patterns. Together, this comprehensive dataset for developing embryos, mature seeds and vegetative tissues provides new insights into molecular mechanisms of seed development with potential for flax crop improvement.

Genome-wide identification of RING finger genes in flax (Linum usitatissimum) and analyses of their evolution

Background

Flax (Linum usitatissimum) is an important crop for its seed oil and stem fiber. Really Interesting New Gene (RING) finger genes play essential roles in growth, development, and biotic and abiotic stress responses in plants. However, little is known about these genes in flax.

Methods

Here, we performed a systematic genome-wide analysis to identify RING finger genes in flax.

Results

We identified 587 RING domains in 574 proteins and classified them into RING-H2 (292), RING-HCa (181), RING-HCb (23), RING-v (53), RING-C2 (31), RING-D (2), RING-S/T (3), and RING-G (2). These proteins were further divided into 45 groups according to domain organization. These genes were located in 15 chromosomes and clustered into three clades according to their phylogenetic relationships. A total of 312 segmental duplicated gene pairs were inferred from 411 RING finger genes, indicating a major contribution of segmental duplications to the RING finger gene family expansion. The non-synonymous/synonymous substitution ratio of the segmentally duplicated gene pairs was less than 1, suggesting that the gene family was under negative selection since duplication. Further, most RING genes in flax were differentially expressed during seed development or in the shoot apex. This study provides useful information for further functional analysis of RING finger genes in flax and to develop gene-derived molecular markers in flax breeding.

Gene Expression Patterns for Proteins With Lectin Domains in Flax Stem Tissues Are Related to Deposition of Distinct Cell Wall Types

The genomes of higher plants encode a variety of proteins with lectin domains that are able to specifically recognize certain carbohydrates. Plants are enriched in a variety of potentially complementary glycans, many of which are located in the cell wall. We performed a genome-wide search for flax proteins with lectin domains and compared the expression of the encoding genes in different stem tissues that have distinct cell wall types with different sets of major polysaccharides. Over 400 genes encoding proteins with lectin domains that belong to different families were revealed in the flax genome; three quarters of these genes were expressed in stem tissues. Hierarchical clustering of the data for all expressed lectins grouped the analyzed samples according to their characteristic cell wall type. Most lectins differentially expressed in tissues with primary, secondary, and tertiary cell walls were predicted to localize at the plasma membrane or cell wall. These lectins were from different families and had various architectural types. Three out of four flax genes for proteins with jacalin-like domains were highly upregulated in bast fibers at the stage of tertiary cell wall deposition. The dynamic changes in transcript level of many genes for lectins from various families were detected in stem tissue over the course of gravitropic response induced by plant gravistimulation. The data obtained in this study indicate a large number of lectin-mediated events in plants and provide insight into the proteins that take part in tissue specialization and reaction to abiotic stress.

The Toolbox for Fiber Flax Breeding: A Pipeline From Gene Expression to Fiber Quality

The goal of any plant breeding program is to improve quality of a target crop. Crop quality is a comprehensive feature largely determined by biological background. To improve the quality parameters of crops grown for the production of fiber, a functional approach was used to search for genes suitable for the effective manipulation of technical fiber quality. A key step was to identify genes with tissue and stage-specific pattern of expression in the developing fibers. In the current study, we investigated the relationship between gene expression evaluated in bast fibers of developing flax plants and the quality parameters of technical fibers measured after plant harvesting. Based on previously published transcriptomic data, two sets of genes that are upregulated in fibers during intrusive growth and tertiary cell wall deposition were selected. The expression level of the selected genes and fiber quality parameters were measured in fiber flax, linseed (oil flax) cultivars, and wild species that differ in type of yield and fiber quality parameters. Based on gene expression data, linear regression models for technical stem length, fiber tensile strength, and fiber flexibility were constructed, resulting in the identification of genes that have high potential for manipulating fiber quality. Chromosomal localization and single nucleotide polymorphism distribution in the selected genes were characterized for the efficacy of their use in conventional breeding and genome editing programs. Transcriptome-based selection is a highly targeted functional approach that could be used during the development of new cultivars of various crops.

Targeted Metagenomics of Retting in Flax: The Beginning of the Quest to Harness the Secret Powers of the Microbiota

The mechanical and chemical properties of natural plant fibers are determined by many different factors, both intrinsic and extrinsic to the plant, during growth but also after harvest. A better understanding of how all these factors exert their effect and how they interact is necessary to be able to optimize fiber quality for use in different industries. One important factor is the post-harvest process known as retting, representing the first step in the extraction of bast fibers from the stem of species such as flax and hemp. During this process microorganisms colonize the stem and produce hydrolytic enzymes that target cell wall polymers thereby facilitating the progressive destruction of the stem and fiber bundles. Recent advances in sequencing technology have allowed researchers to implement targeted metagenomics leading to a much better characterization of the microbial communities involved in retting, as well as an improved understanding of microbial dynamics. In this paper we review how our current knowledge of the microbiology of retting has been improved by targeted metagenomics and discuss how related '-omics' approaches might be used to fully characterize the functional capability of the retting microbiome.

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.

Genomic insights into HSFs as candidate genes for high-temperature stress adaptation and gene editing with minimal off-target effects in flax

Flax (Linum usitatissimum) is a cool season crop commercially cultivated for seed oil and stem fibre production. A comprehensive characterization of the heat shock factor (HSF) candidate genes in flax can accelerate genetic improvement and adaptive breeding for high temperature stress tolerance. We report the genome-wide identification of 34 putative HSF genes from the flax genome, which we mapped on 14 of the 15 chromosomes. Through comparative homology analysis, we classified these genes into three broad groups, and sub-groups. The arrangement of HSF-specific protein motifs, DNA-binding domain (DBD) and hydrophobic heptad repeat (HR-A/B), and exon-intron boundaries substantiated the phylogenetic separation of these genes. Orthologous relationships and evolutionary analysis revealed that the co-evolution of the LusHSF genes was due to recent genome duplication events. Digital and RT-qPCR analyses provided significant evidence of the differential expression of the LusHSF genes in various tissues, at various developmental stages, and in response to high-temperature stress. The co-localization of diverse cis-acting elements in the promoters of the LusHSF genes further emphasized their regulatory roles in the abiotic stress response. We further confirmed DNA-binding sites on the LusHSF proteins and designed guide RNA sequences for gene editing with minimal off-target effects. These results will hasten functional investigations of LusHSFs or assist in devising genome engineering strategies to develop high-temperature stress tolerant flax cultivars.

Flax tubulin and CesA superfamilies represent attractive and challenging targets for a variety of genome- and base-editing applications

Flax is both a valuable resource and an interesting model crop. Despite a long history of flax genetic transformation only one transgenic linseed cultivar has been so far registered in Canada. Implementation and use of the genome-editing technologies that allow site-directed modification of endogenous genes without the introduction of foreign genes might improve this situation. Besides its potential for boosting crop yields, genome editing is now one of the best tools for carrying out reverse genetics and it is emerging as an especially versatile tool for studying basic biology. A complex interplay between the flax tubulin family (6 α-, 14 β-, and 2 γ-tubulin genes), the building block of microtubules, and the CesA (15-16 genes), the subunit of the multimeric cellulose-synthesizing complex devoted to the oriented deposition of the cellulose microfibrils is fundamental for the biosynthesis of the cell wall. The role of the different members of each family in providing specificities to the assembled complexes in terms of structure, dynamics, activity, and interaction remains substantially obscure. Genome-editing strategies, recently shown to be successful in flax, can therefore be useful to unravel the issue of functional redundancy and provide evidence for specific interactions between different members of the tubulin and CesA gene families, in relation to different phase and mode of cell wall biosynthesis.

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.

Transcriptome Analysis of Intrusively Growing Flax Fibers Isolated by Laser Microdissection

The intrusive growth, a type of plant cell elongation occurring in the depths of plant tissues, is characterized by the invasion of a growing cell between its neighbours due to a higher rate of elongation. In order to reveal the largely unknown molecular mechanisms of intrusive growth, we isolated primary flax phloem fibers specifically at the stage of intrusive growth by laser microdissection. The comparison of the RNA-Seq data from several flax stem parts enabled the characterization of those processes occurring specifically during the fiber intrusive elongation. The revealed molecular players are summarized as those involved in the supply of assimilates and support of turgor pressure, cell wall enlargement and modification, regulation by transcription factors and hormones, and responses to abiotic stress factors. The data obtained in this study provide a solid basis for developing approaches to manipulate fiber intrusive elongation, which is of importance both for plant biology and the yield of fiber crops.

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

BACKGROUND

The species relationships within the genus Linum have already been studied several times by means of different molecular and phylogenetic approaches. Nevertheless, a number of ambiguities in phylogeny of Linum still remain unresolved. In particular, the species relationships within the sections Stellerolinum and Dasylinum need further clarification. Also, the question of independence of the species of the section Adenolinum still remains unanswered. Moreover, the relationships of L. narbonense and other species of the section Linum require further clarification. Additionally, the origin of tetraploid species of the section Linum (2n = 30) including the cultivated species L. usitatissimum has not been explored. The present study examines the phylogeny of blue-flowered species of Linum by comparisons of 5S rRNA gene sequences as well as ITS1 and ITS2 sequences of 35S rRNA genes.

RESULTS

High-throughput sequencing has been used for analysis of multicopy rRNA gene families. In addition to the molecular phylogenetic analysis, the number and chromosomal localization of 5S and 35S rDNA sites has been determined by FISH. Our findings confirm that L. stelleroides forms a basal branch from the clade of blue-flowered flaxes which is independent of the branch formed by species of the sect. Dasylinum. The current molecular phylogenetic approaches, the cytogenetic analysis as well as different genomic DNA fingerprinting methods applied previously did not discriminate certain species within the sect. Adenolinum. The allotetraploid cultivated species L. usitatissimum and its wild ancestor L. angustifolium (2n = 30) could originate either as the result of hybridization of two diploid species (2n = 16) related to the modern L. gandiflorum and L. decumbens, or hybridization of a diploid species (2n = 16) and a diploid ancestor of modern L. narbonense (2n = 14).

CONCLUSIONS

High-throughput sequencing of multicopy rRNA gene families allowed us to make several adjustments to the phylogeny of blue-flowered flax species and also reveal intra- and interspecific divergence of the rRNA gene sequences.

Genome-wide identification, characterization, and expression profile of aquaporin gene family in flax (Linum usitatissimum)

Membrane intrinsic proteins (MIPs) form transmembrane channels and facilitate transport of myriad substrates across the cell membrane in many organisms. Majority of plant MIPs have water transporting ability and are commonly referred as aquaporins (AQPs). In the present study, we identified aquaporin coding genes in flax by genome-wide analysis, their structure, function and expression pattern by pan-genome exploration. Cross-genera phylogenetic analysis with known aquaporins from rice, arabidopsis, and poplar showed five subgroups of flax aquaporins representing 16 plasma membrane intrinsic proteins (PIPs), 17 tonoplast intrinsic proteins (TIPs), 13 NOD26-like intrinsic proteins (NIPs), 2 small basic intrinsic proteins (SIPs), and 3 uncharacterized intrinsic proteins (XIPs). Amongst aquaporins, PIPs contained hydrophilic aromatic arginine (ar/R) selective filter but TIP, NIP, SIP and XIP subfamilies mostly contained hydrophobic ar/R selective filter. Analysis of RNA-seq and microarray data revealed high expression of PIPs in multiple tissues, low expression of NIPs, and seed specific expression of TIP3 in flax. Exploration of aquaporin homologs in three closely related Linum species bienne, grandiflorum and leonii revealed presence of 49, 39 and 19 AQPs, respectively. The genome-wide identification of aquaporins, first in flax, provides insight to elucidate their physiological and developmental roles in flax.

miR319, miR390, and miR393 Are Involved in Aluminum Response in Flax (Linum usitatissimum L.)

Acid soils limit agricultural production worldwide. Major reason of crop losses in acid soils is the toxicity of aluminum (Al). In the present work, we investigated expression alterations of microRNAs in flax (Linum usitatissimum L.) plants under Al stress. Flax seedlings of resistant (TMP1919 and G1071/4_k) and sensitive (Lira and G1071/4_o) to Al cultivars and lines were exposed to AlCl₃ solution for 4 and 24 hours. Twelve small RNA libraries were constructed and sequenced using Illumina platform. In total, 97 microRNAs from 18 conserved families were identified. miR319, miR390, and miR393 revealed expression alterations associated with Al treatment of flax plants. Moreover, for miR390 and miR393, the alterations were distinct in sensitive and resistant to Al genotypes. Expression level changes of miR319 and miR390 were confirmed using qPCR analysis. In flax, potential targets of miR319 are TCPs, miR390–TAS3 and GRF5, and miR393–AFB2-coding transcripts. TCPs, TAS3, GRF5, and AFB2 participate in regulation of plant growth and development. The involvement of miR319, miR390, and miR393 in response to Al stress in flax was shown here for the first time. We speculate that these microRNAs play an important role in Al response via regulation of growth processes in flax plants.

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

About 30% of the world's ice-free land area is occupied by acid soils. In soils with pH below 5, aluminum (Al) releases to the soil solution, and becomes highly toxic for plants. Therefore, breeding of varieties that are resistant to Al is needed. Flax (Linum usitatissimum L.) is grown worldwide for fiber and seed production. Al toxicity in acid soils is a serious problem for flax cultivation. However, very little is known about mechanisms of flax resistance to Al and the genetics of this resistance. In the present work, we sequenced 16 transcriptomes of flax cultivars resistant (Hermes and TMP1919) and sensitive (Lira and Orshanskiy) to Al, which were exposed to control conditions and aluminum treatment for 4, 12, and 24 h. In total, 44.9-63.3 million paired-end 100-nucleotide reads were generated for each sequencing library. Based on the obtained high-throughput sequencing data, genes with differential expression under aluminum exposure were revealed in flax. The majority of the top 50 up-regulated genes were involved in transmembrane transport and transporter activity in both the Al-resistant and Al-sensitive cultivars. However, genes encoding proteins with glutathione transferase and UDP-glycosyltransferase activity were in the top 50 up-regulated genes only in the flax cultivars resistant to aluminum. For qPCR analysis in extended sampling, two UDP-glycosyltransferases (UGTs), and three glutathione S-transferases (GSTs) were selected. The general trend of alterations in the expression of the examined genes was the up-regulation under Al stress, especially after 4 h of Al exposure. Moreover, in the flax cultivars resistant to aluminum, the increase in expression was more pronounced than that in the sensitive cultivars. We speculate that the defense against the Al toxicity via GST antioxidant activity is the probable mechanism of the response of flax plants to aluminum stress. We also suggest that UGTs could be involved in cell wall modification and protection from reactive oxygen species (ROS) in response to Al stress in L. usitatissimum. Thus, GSTs and UGTs, probably, play an important role in the response of flax to Al via detoxification of ROS and cell wall modification.