Transcriptome Analysis of Intrusively Growing Flax Fibers Isolated by Laser Microdissection

Expression of FAD and SAD Genes in Developing Seeds of Flax Varieties under Different Growth Conditions

Flax seed is one of the richest plant sources of linolenic acid (LIN) and also contains unsaturated linoleic acid (LIO) and oleic acid (OLE). Stearoyl-ACP desaturases (SADs) and fatty acid desaturases (FADs) play key roles in the synthesis of flax fatty acids (FAs). However, there is no holistic view of which genes from the SAD and FAD families and at which developmental stages have the highest expression levels in flax seeds, as well as the influence of genotype and growth conditions on the expression profiles of these genes. We sequenced flax seed transcriptomes at 3, 7, 14, 21, and 28 days after flowering (DAF) for ten flax varieties with different oil FA compositions grown under three temperature/watering conditions. The expression levels of 25 genes of the SAD, FAD2, and FAD3 families were evaluated. FAD3b, FAD3a, FAD2b-2, SAD3-1, SAD2-1, SAD2-2, SAD3-2, FAD2a-1, and FAD2a-2 had the highest expression levels, which changed significantly during seed development. These genes probably play a key role in FA synthesis in flax seeds. High temperature and insufficient watering shifted the maximum expression levels of FAD and SAD genes to earlier developmental stages, while the opposite trend was observed for low temperature and excessive watering. Differences in the FAD and SAD expression profiles under different growth conditions may affect the FA composition of linseed oil. Stop codons in the FAD3a gene, resulting in a reduced LIN content, decreased the level of FAD3a transcript. The obtained results provide new insights into the synthesis of linseed oil.

Key FAD2, FAD3, and SAD Genes Involved in the Fatty Acid Synthesis in Flax Identified Based on Genomic and Transcriptomic Data

FAD (fatty acid desaturase) and SAD (stearoyl-ACP desaturase) genes play key roles in the synthesis of fatty acids (FA) and determination of oil composition in flax (Linum usitatissimum L.). We searched for FAD and SAD genes in the most widely used flax genome of the variety CDC Bethune and three available long-read assembled flax genomes-YY5, 3896, and Atlant. We identified fifteen FAD2, six FAD3, and four SAD genes. Of all the identified genes, 24 were present in duplicated pairs. In most cases, two genes from a pair differed by a significant number of gene-specific SNPs (single nucleotide polymorphisms) or even InDels (insertions/deletions), except for FAD2a-1 and FAD2a-2, where only seven SNPs distinguished these genes. Errors were detected in the FAD2a-1, FAD2a-2, FAD3c-1, and FAD3d-2 sequences in the CDC Bethune genome assembly but not in the long-read genome assemblies. Expression analysis of the available transcriptomic data for different flax organs/tissues revealed that FAD2a-1, FAD2a-2, FAD3a, FAD3b, SAD3-1, and SAD3-2 were specifically expressed in embryos/seeds/capsules and could play a crucial role in the synthesis of FA in flax seeds. In contrast, FAD2b-1, FAD2b-2, SAD2-1, and SAD2-2 were highly expressed in all analyzed organs/tissues and could be involved in FA synthesis in whole flax plants. FAD2c-2, FAD2d-1, FAD3c-1, FAD3c-2, FAD3d-1, FAD3d-2, SAD3-1, and SAD3-2 showed differential expression under stress conditions-Fusarium oxysporum infection and drought. The obtained results are essential for research on molecular mechanisms of fatty acid synthesis, FAD and SAD editing, and marker-assisted and genomic selection for breeding flax varieties with a determined fatty acid composition of oil.

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.

Genetic Determinants of Fiber-Associated Traits in Flax Identified by Omics Data Integration

In this paper, we explore potential genetic factors in control of flax phenotypes associated with fiber by mining a collection of 306 flax accessions from the Federal Research Centre of the Bast Fiber Crops, Torzhok, Russia. In total, 11 traits were assessed in the course of 3 successive years. A genome-wide association study was performed for each phenotype independently using six different single-locus models implemented in the GAPIT3 R package. Moreover, we applied a multivariate linear mixed model implemented in the GEMMA package to account for trait correlations and potential pleiotropic effects of polymorphisms. The analyses revealed a number of genomic variants associated with different fiber traits, implying the complex and polygenic control. All stable variants demonstrate a statistically significant allelic effect across all 3 years of the experiment. We tested the validity of the predicted variants using gene expression data available for the flax fiber studies. The results shed new light on the processes and pathways associated with the complex fiber traits, while the pinpointed candidate genes may be further used for marker-assisted selection.

Laticifer growth pattern is guided by cytoskeleton organization

Laticifers are secretory structures that produce latex, forming a specialized defense system against herbivory. Studies using anatomical approaches to investigate laticifer growth patterns have described their origin; however, their mode of growth, i.e., whether growth is intrusive or diffuse, remains unclear. Studies investigating how cytoskeleton filaments may influence laticifer shape establishment and growth patterns are lacking. In this study, we combined microtubule immunostaining and developmental anatomy to investigate the growth patterns in different types of laticifers. Standard anatomical methods were used to study laticifer development. Microtubules were labelled through immunolocalization of α-tubulin in three types of laticifers from three different plant species: nonanastomosing (Urvillea ulmacea), anastomosing unbranched with partial degradation of terminal cell walls (Ipomoea nil), and anastomosing branched laticifers with early and complete degradation of terminal cell walls (Asclepias curassavica). In both nonanastomosing and anastomosing laticifers, as well as in differentiating meristematic cells, parenchyma cells and idioblasts, microtubules were perpendicularly aligned to the cell growth axis. The analyses of laticifer microtubule orientation revealed an arrangement that corresponds to those cells that grow diffusely within the plant body. Nonanastomosing and anastomosing laticifers, branched or not, have a pattern which indicates diffuse growth. This innovative study on secretory structures represents a major advance in the knowledge of laticifers and their growth mode.

Bouncing back stronger: Diversity, structure, and molecular regulation of gelatinous fiber development

Gelatinous fibers (G-fibers) are specialized contractile cells found in a diversity of vascular plant tissues, where they provide mechanical support and/or facilitate plant mobility. G-fibers are distinct from typical fibers by the presence of an innermost thickened G-layer, comprised mainly of axially oriented cellulose microfibrils. Despite the disparate developmental origins-tension wood fibers from the vascular cambium or primary phloem fibers from the procambium-G-fiber development, composition, and molecular signatures are remarkably similar; however, important distinctions do exist. Here, we synthesize current knowledge of the phylogenetic diversity, compositional makeup, and the molecular profiles that characterize G-fiber development and highlight open questions for future investigation.

Using FIBexDB for In-Depth Analysis of Flax Lectin Gene Expression in Response to Fusarium oxysporum Infection

Plant proteins with lectin domains play an essential role in plant immunity modulation, but among a plurality of lectins recruited by plants, only a few members have been functionally characterized. For the analysis of flax lectin gene expression, we used FIBexDB, which includes an efficient algorithm for flax gene expression analysis combining gene clustering and coexpression network analysis. We analyzed the lectin gene expression in various flax tissues, including root tips infected with Fusarium oxysporum. Two pools of lectin genes were revealed: downregulated and upregulated during the infection. Lectins with suppressed gene expression are associated with protein biosynthesis (Calreticulin family), cell wall biosynthesis (galactose-binding lectin family) and cytoskeleton functioning (Malectin family). Among the upregulated lectin genes were those encoding lectins from the Hevein, Nictaba, and GNA families. The main participants from each group are discussed. A list of lectin genes, the expression of which can determine the resistance of flax, is proposed, for example, the genes encoding amaranthins. We demonstrate that FIBexDB is an efficient tool both for the visualization of data, and for searching for the general patterns of lectin genes that may play an essential role in normal plant development and defense.

Do woody vines use gelatinous fibers to climb?

Many plant movements are facilitated by contractile cells called gelatinous fibers (G-fibers), but how G-fibers function in the climbing movements of woody vines remains underexplored. In this Insight, we compare the presence and distribution of G-fibers in the stems of stem-twiners, which wrap around supports, with non-stem-twiners, which attach to supports via tendrils or adventitious roots. An examination of 164 species spanning the vascular plant phylogeny reveals that G-fibers are common in stem-twiners but scarce in non-stem-twiners, suggesting that G-fibers are preferentially formed in the organ responsible for movement. When present, G-fibers are in the xylem, phloem, pericycle, and/or cortex. We discuss the hypothesis that G-fibers are foundational to plant movement and highlight research opportunities concerning G-fiber development and function.

Genes Associated with the Flax Plant Type (Oil or Fiber) Identified Based on Genome and Transcriptome Sequencing Data

As a result of the breeding process, there are two main types of flax (Linum usitatissimum L.) plants. Linseed is used for obtaining seeds, while fiber flax is used for fiber production. We aimed to identify the genes associated with the flax plant type, which could be important for the formation of agronomically valuable traits. A search for polymorphisms was performed in genes involved in the biosynthesis of cell wall components, lignans, fatty acids, and ion transport based on genome sequencing data for 191 flax varieties. For 143 of the 424 studied genes (4CL, C3'H, C4H, CAD, CCR, CCoAOMT, COMT, F5H, HCT, PAL, CTL, BGAL, ABC, HMA, DIR, PLR, UGT, TUB, CESA, RGL, FAD, SAD, and ACT families), one or more polymorphisms had a strong correlation with the flax type. Based on the transcriptome sequencing data, we evaluated the expression levels for each flax type-associated gene in a wide range of tissues and suggested genes that are important for the formation of linseed or fiber flax traits. Such genes were probably subjected to the selection press and can determine not only the traits of seeds and stems but also the characteristics of the root system or resistance to stresses at a particular stage of development, which indirectly affects the ability of flax plants to produce seeds or fiber.

Elongation of wood fibers combines features of diffuse and tip growth

Xylem fibers are highly elongated cells that are key constituents of wood, play major physiological roles in plants, comprise an important terrestrial carbon reservoir, and thus have enormous ecological and economic importance. As they develop, from fusiform initials, their bodies remain the same length while their tips elongate and intrude into intercellular spaces. To elucidate mechanisms of tip elongation, we studied the cell wall along the length of isolated, elongating aspen xylem fibers and used computer simulations to predict the forces driving the intercellular space formation required for their growth. We found pectin matrix epitopes (JIM5, LM7) concentrated at the tips where cellulose microfibrils have transverse orientation, and xyloglucan epitopes (CCRC-M89, CCRC-M58) in fiber bodies where microfibrils are disordered. These features are accompanied by changes in cell wall thickness, indicating that while the cell wall elongates strictly at the tips, it is deposited all over fibers. Computer modeling revealed that the intercellular space formation needed for intrusive growth may only require targeted release of cell adhesion, which allows turgor pressure in neighboring fiber cells to 'round' the cells creating spaces. These characteristics show that xylem fibers' elongation involves a distinct mechanism that combines features of both diffuse and tip growth.

OMICs, Epigenetics, and Genome Editing Techniques for Food and Nutritional Security

The incredible success of crop breeding and agricultural innovation in the last century greatly contributed to the Green Revolution, which significantly increased yields and ensures food security, despite the population explosion. However, new challenges such as rapid climate change, deteriorating soil, and the accumulation of pollutants require much faster responses and more effective solutions that cannot be achieved through traditional breeding. Further prospects for increasing the efficiency of agriculture are undoubtedly associated with the inclusion in the breeding strategy of new knowledge obtained using high-throughput technologies and new tools in the future to ensure the design of new plant genomes and predict the desired phenotype. This article provides an overview of the current state of research in these areas, as well as the study of soil and plant microbiomes, and the prospective use of their potential in a new field of microbiome engineering. In terms of genomic and phenomic predictions, we also propose an integrated approach that combines high-density genotyping and high-throughput phenotyping techniques, which can improve the prediction accuracy of quantitative traits in crop species.

Genome-wide identification and expression analysis of the WRKY transcription factor family in flax (Linum usitatissimum L.)

BACKGROUND

Members of the WRKY protein family, one of the largest transcription factor families in plants, are involved in plant growth and development, signal transduction, senescence, and stress resistance. However, little information is available about WRKY transcription factors in flax (Linum usitatissimum L.).

RESULTS

In this study, comprehensive genome-wide characterization of the flax WRKY gene family was conducted that led to prediction of 102 LuWRKY genes. Based on bioinformatics-based predictions of structural and phylogenetic features of encoded LuWRKY proteins, 95 LuWRKYs were classified into three main groups (Group I, II, and III); Group II LuWRKYs were further assigned to five subgroups (IIa-e), while seven unique LuWRKYs (LuWRKYs 96-102) could not be assigned to any group. Most LuWRKY proteins within a given subgroup shared similar motif compositions, while a high degree of motif composition variability was apparent between subgroups. Using RNA-seq data, expression patterns of the 102 predicted LuWRKY genes were also investigated. Expression profiling data demonstrated that most genes associated with cellulose, hemicellulose, or lignin content were predominantly expressed in stems, roots, and less in leaves. However, most genes associated with stress responses were predominantly expressed in leaves and exhibited distinctly higher expression levels in developmental stages 1 and 8 than during other stages.

CONCLUSIONS

Ultimately, the present study provides a comprehensive analysis of predicted flax WRKY family genes to guide future investigations to reveal functions of LuWRKY proteins during plant growth, development, and stress responses.

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.

Genome-wide identification of fasciclin-like arabinogalactan proteins in jute and their expression pattern during fiber formation

Fasciclin-like arabinogalactan proteins (FLAs), a class of arabinogalactan proteins (AGPs) are involved in plant growth and development via cell communication and adhesion. FLAs were also associated with fiber and wood formation in plants but no information is available about the roles of FLA proteins during fibre development of jute. Here, we performed molecular characterization, evolutionary relationship and expression profiling of FLAs proteins in jute (Corchorus olitorius). In total, nineteen CoFLA genes have been identified in jute genome, which were divided into four classes like FLAs of other species based on protein structure and similarity. All CoFLAs have N-terminal signal peptide and one or two FAS domain while two FLAs lack well defined AGP region and eight FLAs were devoid of C-terminal glycosylphosphatidylinositol (GPI) anchor. Expression analysis of different regions of jute stem suggested their involvement in different fiber development stages. Four genes CoFLA 11, 12, 20, and 23 were highly or predominately expressed in fiber containing bark tissues while the expression levels of six CoFLA genes 02, 03, 04, 06, 14 and 19 were comparatively higher in stick. Higher transcripts levels of CoFLA 12 and 20 in the middle bark tissues suggest their involvement in fiber elongation. In contrast, the CoFLA 11 and 23 were more expressed in bottom bark tissues suggesting their potential involvement in secondary cell wall synthesis. Our study can serve as solid foundation for further functional exploration of FLAs and in future breeding program of jute aiming fiber improvement.

Transcriptome Analysis Reveals a Gene Expression Pattern Associated with Fuzz Fiber Initiation Induced by High Temperature in Gossypium barbadense

Gossypium barbadense is an important source of natural textile fibers, as is Gossypium hirsutum. Cotton fiber development is often affected by various environmental factors, such as abnormal temperature. However, little is known about the underlying mechanisms of temperature regulating the fuzz fiber initiation. In this study, we reveal that high temperatures (HT) accelerate fiber development, improve fiber quality, and induced fuzz initiation of a thermo-sensitive G. barbadense variety L7009. It was proved that fuzz initiation was inhibited by low temperature (LT), and 4 dpa was the stage most susceptible to temperature stress during the fuzz initiation period. A total of 43,826 differentially expressed genes (DEGs) were identified through comparative transcriptome analysis. Of these, 9667 were involved in fiber development and temperature response with 901 transcription factor genes and 189 genes related to plant hormone signal transduction. Further analysis of gene expression patterns revealed that 240 genes were potentially involved in fuzz initiation induced by high temperature. Functional annotation revealed that the candidate genes related to fuzz initiation were significantly involved in the asparagine biosynthetic process, cell wall biosynthesis, and stress response. The expression trends of sixteen genes randomly selected from the RNA-seq data were almost consistent with the results of qRT-PCR. Our study revealed several potential candidate genes and pathways related to fuzz initiation induced by high temperature. This provides a new view of temperature-induced tissue and organ development in Gossypium barbadense.

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.

The cinnamyl alcohol dehydrogenase gene family is involved in the response to Fusarium oxysporum in resistant and susceptible flax genotypes

Flax (Linum usitatissimum L.) is used for the production of textile, oils, pharmaceuticals, and composite materials. Fusarium wilt, caused by the fungus Fusarium oxysporum f. sp. lini, is a very harmful disease that reduces flax production. Flax cultivars that are resistant to Fusarium wilt have been developed, and the genes that are involved in the host response to F. oxysporum have been identified. However, the mechanisms underlying resistance to this pathogen remain unclear. In the present study, we used transcriptome sequencing data obtained from susceptible and resistant flax genotypes grown under control conditions or F. oxysporum infection. Approximately 250 million reads, generated with an Illumina NextSeq instrument, were analyzed. After filtering to exclude the F. oxysporum transcriptome, the remaining reads were mapped to the L. usitatissimum genome and quantified. Then, the expression levels of cinnamyl alcohol dehydrogenase (CAD) family genes, which are known to be involved in the response to F. oxysporum, were evaluated in resistant and susceptible flax genotypes. Expression alterations in response to the pathogen were detected for all 13 examined CAD genes. The most significant differences in expression between control and infected plants were observed for CAD1B, CAD4A, CAD5A, and CAD5B, with strong upregulation of CAD1B, CAD5A, and CAD5B and strong downregulation of CAD4A. When plants were grown under the same conditions, the expression levels were similar in all studied flax genotypes for most CAD genes, and statistically significant differences in expression between resistant and susceptible genotypes were only observed for CAD1A. Our study indicates the strong involvement of CAD genes in flax response to F. oxysporum but brings no evidence of their role as resistance gene candidates. These findings contribute to the understanding of the mechanisms underlying the response of flax to F. oxysporum infection and the role of CAD genes in stress resistance.

Intrusive Growth of Phloem Fibers in Flax Stem: Integrated Analysis of miRNA and mRNA Expression Profiles

Phloem fibers are important elements of plant architecture and the target product of many fiber crops. A key stage in fiber development is intrusive elongation, the mechanisms of which are largely unknown. Integrated analysis of miRNA and mRNA expression profiles in intrusivelygrowing fibers obtained by laser microdissection from flax (Linum usitatissimum L.) stem revealed all 124 known flax miRNA from 23 gene families and the potential targets of differentially expressed miRNAs. A comparison of the expression between phloem fibers at different developmental stages, and parenchyma and xylem tissues demonstrated that members of miR159, miR166, miR167, miR319, miR396 families were down-regulated in intrusively growing fibers. Some putative target genes of these miRNA families, such as those putatively encoding growth-regulating factors, an argonaute family protein, and a homeobox-leucine zipper family protein were up-regulated in elongating fibers. miR160, miR169, miR390, and miR394 showed increased expression. Changes in the expression levels of miRNAs and their target genes did not match expectations for the majority of predicted target genes. Taken together, poorly understood intrusive fiber elongation, the key process of phloem fiber development, was characterized from a miRNA-target point of view, giving new insights into its regulation.