Transcriptome Analysis of Intrusively Growing Flax Fibers Isolated by Laser Microdissection

Gorshkova, Tatyana; Chernova, Tatyana; Mokshina, Natalia; Gorshkov, Vladimir; Kozlova, Liudmila; Gorshkov, Oleg

doi:10.1038/s41598-018-32869-2

Download PDF

Article
Open access
Published: 01 October 2018

Transcriptome Analysis of Intrusively Growing Flax Fibers Isolated by Laser Microdissection

Tatyana Gorshkova¹,
Tatyana Chernova¹,
Natalia Mokshina¹,
Vladimir Gorshkov¹,
Liudmila Kozlova¹ &
…
Oleg Gorshkov¹

Scientific Reports volume 8, Article number: 14570 (2018) Cite this article

2092 Accesses
33 Citations
1 Altmetric
Metrics details

Subjects

Abstract

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.

A transcriptomic view to wounding response in young Scots pine stems

Article Open access 12 February 2021

Transcriptome analyses suggest that changes in fungal endophyte lifestyle could be involved in grapevine bud necrosis

Article Open access 11 June 2020

Induction of PrMADS10 on the lower side of bent pine tree stems: potential role in modifying plant cell wall properties and wood anatomy

Article Open access 12 December 2019

Introduction

Enlargement of individual cells is a fundamental process of higher plant development. This increase of cell size in a multicellular organism can proceed in three major ways: by symplastic, protrusive, or intrusive growth^1,2,3. During symplastic growth, the involved cells increase their surface simultaneously, retaining the contacts between them through the joint middle lamellae and plasmodesmata. Symplastic growth is the most common type, usually characterizing cells in the elongation zones. Protrusive growth is an outgrowth of a cell part; it takes place only in the epidermal layer and leads to the formation of trichomes, such as root hairs, or hairs of the cotton seed (“cotton fibers”). The third type of cell enlargement, the intrusive growth, is a type characterized by a higher rate of a cell enlargement compared with the adjacent cells. During intrusive elongation, a cell squeezes in between the surrounding cells, splitting their middle lamellae. Several types of plant cells can perform the intrusive growth, with fibers being the classical example^2,4,5,6. In recent years, the intrusive growth has attracted much attention and has been the subject of discussions in several reviews^3,7,8,9, all of which have pointed to the lack of knowledge of the molecular mechanisms underpinning intrusive growth of fibers.

The fiber intrusive elongation is very important for the development of plant architecture as well as for the yield and quality of commercial plant fibers^6,7,10. Fiber bundles develop within plant organs mainly due to the intrusive elongation of individual fibers¹⁰. The key processes that a fiber must accomplish during the intrusive growth include enlarging the cell wall surface (via synthesis and secretion of cell wall polymers and loosening of the existing layers) and increasing the vacuole volume (via water uptake, osmolyte accumulation, synthesis of membrane components). Furthermore, an intrusively growing fiber requires conditions that enable it to squeeze in between neighbouring cells—the dissolution of middle lamellae of the cells on the way of invasion and the abolishment of possible wound effects. Characterization of any of these processes in intrusively growing fibers from a molecular point of view remains quite poor. There is not a single paper with the analysis of fibers isolated at this developmental stage by comprehensive high-throughput techniques, like RNA-Seq. Such experiments are highly demanded to approach understanding of the mechanisms of intrusive elongation and its regulation, but have, at most, been performed on complex samples, which together with the intrusively growing fibers contained other cell types at other developmental stages. The quite limited examples of the studies are the description of transcriptome landscapes in sections of the young xylem from aspen wood, which contain xylem vessels and ray parenchyma together with intrusively growing xylem fibers¹¹, [cluster e1], or flax stem portion (TOP) that includes intensively elongating fibers together with many other tissues (e.g. epidermis, phloem, parenchyma, cambium, xylem¹².

To perform the RNA-Seq analysis specifically for intrusively growing fibers, we have chosen primary phloem fibers of flax due to the well-characterized location within the stem and duration of the process². Primary flax phloem fibers originate from the procambium, close to the apical meristem⁴, and at the earliest stages of development they grow symplastically with the surrounding tissues¹³. The intrusive elongation is initiated at a 1–2-mm distance from the stem apex and continues several centimetres down the stem until the so-called “snap point”, where the transition of fibers to cell wall thickening occurs^13,14. These morphological markers allow the study of stage-specific parameters of fiber development. However, the combination of a significant cell length with only a primary cell wall makes the intrusively growing fibers quite prone to breakage, which hinders their isolation.

To elucidate the molecular events associated with the fiber intrusive elongation, we isolated the intrusively growing fibers through laser microdissection and, for the first time, performed the RNA-Seq analysis for such cells. Comparison of the obtained data with transcriptome landscapes in the tissues from the meristematic stem portions and in fibers at an advanced developmental stage enabled us (1) to characterize the overall picture of the transcriptome in the intrusively growing fibers, and (2) to determine which genes were subjected to up- and down-regulation of transcription at this critical stage of fiber development.

Results and Discussion

General transcriptome characteristics of the analyzed samples

To obtain a full picture of the genes expressed during the intrusive growth of fibers in planta, we performed a comparative RNA-Seq analysis with the mRNAs isolated from the apical part of the stem (APEX) and the intrusively growing phloem fibers (iFIB) (Fig. 1).

A total of 104,827,884 reads with quality scores above Q30 were acquired; of these, approx. 90% (92,277,439) were successfully mapped onto the flax reference genome. Normalized expression values of genes were estimated as the fragments per kilobase of exon per million mapped fragments (FPKM). Of the 43,486 protein-coding genes predicted in the whole-genome assembly of flax, 28,983 were found with FPKM >1 (see Supplementary Table S1). FPKM values for the genes in APEX sample were in good agreement (Pearson correlation coefficient for the log₂ values was r = 0.87) with the data reported by Zhang and Deyholos¹⁵ for a similar sample verifying our obtained results. The analysis in the software package Cufflinks revealed 10,373 genes with FPKM >16 which could be used for further analysis, as recommended¹⁶. Using the fold change below 0.33 (down-regulated) and above 3 (up-regulated) as cut-offs for iFIB versus APEX, and FPKM >16 at least in one sample, a total of 3,731 differentially expressed genes (DEGs) were identified.

To visualize the changes in gene expression upon transition to the intrusive growth, we compared the transcriptome of iFIB with that of APEX in MapMan¹⁷, supplementing it with the results of the search of functional classes that were either over- or underrepresented in DEGs between the APEX and iFIB (Fig. 2A) using overrepresentation analysis by the PageMan software¹⁸. Major metabolic pathways such as photosynthesis, carbohydrate metabolism, cell wall, hormone metabolism, stress, enzyme families, regulation of transcription, protein degradation, signaling, cell organization and transport were overrepresented among expressed iFIB genes. The highest proportion of expressed genes occurred in the pathways designated as lipid and secondary metabolism, regulation of transcription, DNA and protein synthesis, cell cycle.

Twelve genes that had various dynamics of expression between APEX and iFIB were selected for verification by the qRT-PCR. The expression values calculated in RNA-Seq experiments were in good agreement with the data obtained by qRT-PCR (Fig. 2B).

To discriminate the genes that are specifically up- or down-regulated during the intrusive elongation, and those genes that had similar transcript abundance at a later stage of fiber development, we compared our data with the published results of transcriptome analysis for fibers isolated at the stage of the tertiary cell wall formation from 10 cm stem portions taken below the snap point (FIB)¹² (Fig. 1). Many of the up-regulated differentially expressed genes revealed by comparison of APEX and iFIB samples had similar FPKM values in iFIB and FIB, suggesting that while their expression is important for fiber development, they are not stage-specific. A gene was considered as specifically up- or down-regulated in iFIB, if the fold change in transcript abundance exceeded 3 (|log₂FC| > 1.6), in comparison with both APEX and FIB. With such a cut-off, 242 genes were up-regulated, while 1151 were down-regulated stage-specifically in iFIB (see Supplementary Table S2). Next, we present the obtained data considering the most important aspects of cell physiology in the intrusively growing fibers.

Supply of assimilates and support of turgor pressure

In primary phloem flax fibers intrusive growth leads to the increase of cell surface and volume several thousand times. When these fibers initiate in the procambium they are c. 5 µm in width and less than 100 µm in length, while mature fibers are 20–30-µm wide with a modal length of c. 25,000 µm^7,13. The characteristic feature of plant cell growth is that an enlargement of cell volume is mainly achieved by the increasing the vacuole through the accumulation of osmolytes and an inward water flux¹⁹. The latter is largely provided by aquaporins, located in both the plasma membrane (PIP) and the tonoplast (TIP). Different members of multigene families encoding aquaporins were expressed in iFIB and APEX; the genes for PIPs (Lus10014978, Lus10015083) and TIP (Lus10025808) were considerably up-regulated in the intrusively growing fibers (Fig. 3A, see also Lus10032283, Lus10024651, Lus10022611, and Lus10021510 in Supplementary Table S1). This is typical for cells that undergo rapid elongation²⁰, such as seed trichomes²¹, maize root cells²², and cells in the internodes of deepwater rice²³. Activation of aquaporin expression was also detected in samples with mixtures of different cell types collected from the flax stem part containing the elongating fibers (ref.¹², sample TOP). Transcripts for NIPs, SIPs, and XIPs were barely detected in the flax fibers (see Supplementary Table S1) as well as in the flax stem tissues in general²⁴.

The nature of the osmoticum in the elongating fibers has not been elucidated yet. Ions were strong candidates. However, our transcriptomic analysis did not reveal any significantly up-regulated genes for ion transporters. In the agreement, the content of major ions was similar between the fibers and surrounding parenchyma, same as in fibers at an advanced stage of development (Fig. 3B). All salt overly-sensitive (SOS) genes were neither differentially nor highly expressed in iFIB (see Supplementary Table S1). Together with that, the genes for the three types of the carbohydrate transporters located in the plasma membrane—sucrose-proton symporter 2 (Lus10009428, Lus10006066); SWEET, also named MtN3²⁵ (Lus10016742, Lus10022436); and sugar transport protein 1 (Lus10002597)—were up-regulated in the intrusively growing fibers, similarly to the gene Lus10028303 for a putative early-responsive to dehydration likesugar transporter 6 (ERDL6), which is located on the tonoplast²⁶. In addition the specific up-regulation of genes Lus10001087 and Lus10040134 for chloroplast beta-amylase that cleaves starch into low-molecular weight sugars, and for vacuolar invertase (glycosyl hydrolases family 32 protein, Lus10024605) that cuts sucrose into two monosaccharides, rising up osmolarity, occurred (Fig. 3C). Elevated abundance of mRNA for sugar transporters, specific beta-amylases and invertase suggests low-molecular carbohydrates as the most probable osmolytes involved in fiber elongation.

The list of genes for transporters specifically up-regulated in iFIB also included amino acid-proton symporter AAP8 (Lus10012444), member of the transmembrane amino acid transporter family (Lus10005099) and inositol-proton symporters (Lus10042820, Lus10025565) involved in myo-inositol release from vacuoles²⁷ (Fig. 3C). Activation of inositol metabolism in intrusively growing fibers is supported by the increased transcript abundance of genes for myo-inositol-1-phosphate synthase 2 (MIPS2, Lus10000692)—the enzyme that catalyzes the rate-limiting step in inositol synthesis²⁸, and myo-inositol oxygenase 4 (Lus10012923) involved in the biosynthesis of UDP-glucuronic acid. Additionally, Lus10001542 for the protein from Sec14p-like phosphatidylinositol transfer family (Patellin 3) that may be involved in membrane trafficking events and regulation of phosphoinositide signaling²⁹ was specifically up-regulated in iFIB (Fig. 3C).

Plant cell elongation is coupled to acidification of the apoplast¹⁹. In accordance with that, the genes Lus10026224 and Lus10001437 for H(+)-ATPase that serves as a proton pump were up-regulated in intrusively growing fibers (Fig. 3C). However, the abundance of mRNAs for these genes remained similar at further stages of fiber development. Transcription of the genes for vacuolar proton ATPase (Lus10017490, Lus10028796) was activated less substantially, but more specifically for the elongating fibers (see Supplementary Table S1).

Since plasmodesmata in fibers are disrupted during the intrusive elongation¹³, the intercellular transport is possible only through the fiber apoplast. The increased abundance of mRNA for the sucrose/proton importers (Fig. 3C) may indicate the activation of assimilate transport from the apoplast. Assimilates necessary for cell expansion are also supplied by the fibers themselves, which have fully developed chloroplasts¹³. Genes encoding the components of photosynthesis machinery were up-regulated in the intrusively growing fibers when compared with the meristematic region, including genes for the various components of photosystem II (Lus10004895, Lus10001741, Lus10001169, Lus10006593, Lus10006594), of chloroplast ATP synthase complex (Lus10004894, Lus10032827, Lus10009173) and many others (see Supplementary Table S1).

Symplasmic isolation of cells may be caused by callose plug formation at the plasmodesmata^30,31. Significant up-regulation of genes for a glucan synthase-like protein (GSL) similar to callose synthase (Lus10037471, Lus10003917, Lus10042478, Lus10020751, Lus10040892, Lus10003920) was specifically pronounced in the intrusively growing fibers (Fig. 3D). Symplasmic isolation of fibers may be required to maintain a high turgor pressure, as shown for elongating cotton seed trichomes³². Together with the genes for proteins involved in callose biosynthesis, the genes for callose-binding proteins (carbohydrate-binding X8 domain superfamily proteins; Lus10001740, Lus10001167, Lus10039059) were specifically up-regulated at the intrusive elongation stage (Fig. 3D). Such proteins lack a catalytic domain but can bind the chains of β-1,3-linked glucose in either callose or laminarin³³. Possible involvement of the X8 domain superfamily proteins in symplasmic isolation of the intrusively growing fibers is suggested by the increased callose accumulation and reduced diffusion of green fluorescent protein from cell to cell in Arabidopsis lines with an overexpressed gene of a similar protein, called plasmodesmata callose binding (PDCB) proteins³⁴. The partial hydrolysis of callose, if necessary, can be provided by the O-glycosyl hydrolases of family 17; transcription of numerous genes for such proteins was finely adjusted in analyzed samples (Fig. 3D). Genes for the proteins of carbohydrate-binding X8 domain superfamily, O-glycosyl hydrolases family 17, and for the vacuolar invertase similar to those identified in flax fibers were up-regulated in the mixed cell-type samples of aspen wood that contained intrusively growing xylem fibers¹¹, [cluster e1].

Additionally, genes for putative trehalose-phosphatase/synthase 9 (TPS9) (Lus10030635, Lus10030853) were specifically up-regulated in iFIB sample of flax (Fig. 3B) and assigned to gene cluster e1 upon analysis of aspen wood samples¹¹. Trehalose—a non-reducing disaccharide of glucose, and TPS can be involved in cellular morphogenesis³⁵. The revealed similarities (together with other ones listed below) in fiber intrusive growth between different plant species strengthen our conclusions.

Cell wall enlargement and modification

The primary cell wall of flax phloem fibers is similar to the type I primary cell walls typical for dicots³⁶. The major polymers are cellulose, xyloglucan, and pectins. The middle lamella is enriched in pectins, the major constituents of which are homogalacturonan, rhamnogalacturonan I (RG-I), and rhamnogalacturonan II. Among the genes potentially involved in biosynthesis of these non-cellulosic polymers—those encoding corresponding glycosyltransferases and enzymes for interconversion of their monosaccharidic substrates—none were considerably up-regulated in intrusively growing fibers relative to the APEX sample. A moderate up-regulation was detected for several genes of proteins involved in cellulose synthesis machinery, including the CESAs characteristic of the primary cell wall—CESA3 (Lus10007538, Lus10012198) and CESA6 (Lus10041063), cellulase KORRIGAN1 (Lus10007585), and cellulose synthase interactive protein 1 (Lus10019679) (see Supplementary Table S1). At the same time, considerable changes in transcription of genes for the enzymes for the in muro modification of the deposited polysaccharides and for other cell wall proteins were pronounced (Fig. 4).

As the intrusive elongation involves splitting of the middle lamella, the corresponding modification of pectins is anticipated. The families of plant enzymes involved in the post-synthetic modification of pectins include polygalacturonases, pectin methylesterases (PMEs) and their negative regulators—pectin methylesterase inhibitor (PMEI) proteins, pectate lyases, pectin acetylesterases, and rhamnogalacturonan lyases. Different sets of genes from large families encoding enzymes for pectin modification were expressed in iFIB and APEX (Fig. 4A).

A stage-specific up-regulation of expression in intrusively elongating fibers was characteristic for representatives of two gene families—those annotated as pectin lyase-like superfamily proteins (Lus10023679, Lus10013666, Lus10011885, Lus10022817) and pectin acetylesterases (Lus10033301, Lus10034762) (Fig. 4A). All the up-regulated in iFIB genes of pectin lyase-like superfamily proteins belonged to the subfamily I of pectate lyases (EC 4.2.2.2, PF00544)—the enzymes that depolymerized demethylated homogalacturonan by transelimination mechanism³⁷. Of the numerous genes for polygalacturonases in the flax genome, none was substantially up-regulated in the intrusively growing fibers (see Supplementary Table S1). Out of many tens of PME and PMEI genes³⁸, only Lus10026435 had especially abundant mRNA in iFIB as compared to other analyzed samples (Fig. 4A). Herein, several genes for polygalacturonases, PME (both Type I and Type II) and PMEI were down-regulated in iFIB as compared to APEX (Fig. 4A).

Pectin acetylesterases remove acetyl groups that can be present at positions O-2 and/or O-3 of GalA residues in both RG-I and homogalacturonan³⁹. The specific up-regulation of genes for such enzymes in the course of fiber intrusive growth points to the importance of pectin’s acetylation degree for this process. This brings up the need for further studies, since currently the biological function of pectin acetylation and the consequences of acetylesterase action are not well understood³⁹.

Among the 15 annotated genes for the rhamnogalacturonan lyase protein family, none had an FPKM >16 in the APEX and iFIB samples; herewith, several such genes were highly expressed during tertiary cell wall formation¹². Modification of RG-I in the intrusively growing fibers is suggested from the stage-specific up-regulation of genes for α-arabinofuranosidase (Lus10038988) and β-galactosidase (Lus10040557, Lus10006009) (Fig. 4B) that may trim off the neutral side chains of RG-I.

The major cross-linking glycan in the primary cell walls of dicots is xyloglucan. Many studies provided evidences favouring the hypothesis of xyloglucan endotransglucosylase/hydrolases (XET/XTH) mediating cell elongation^40,41,42. However, gene expression for several XTHs was down-regulated in iFIB, as well as of the genes for fasciclin-like arabinogalactan proteins (FLAs), arabinogalactan proteins (AGPs) and extensin-like proteins (Fig. 4B).

The best-known proteins acting as wall-loosening agents are the expansins⁴³. Of the 65 genes encoding expansins and expansin-like proteins present in the flax genome, many had higher expression in APEX than in iFIB. However, some expansin genes (Lus10040801, Lus10019088) were specifically up-regulated in the intrusively growing fibers (Fig. 4B).

In sum, cell wall extension during the intrusive growth is suggested to arise from the action of up-regulated expansins, probable modification of RG-I, and cleavage of homogalacturonan, which occurs by action of pectate lyases rather than of polygalacturonases and may be coupled to deacetylation of the polymer by pectin acetyl esterases.

Regulation by transcription factors and hormones

Key regulatory proteins that play significant roles in many biological processes, by activating and/or repressing transcription, are called transcription factors (TFs). No information about TFs involved specifically in the fiber intrusive growth has yet been obtained. According to PlantTFDB (http://planttfdb.cbi.pku.edu.cn/), the flax genome contains 2481 genes for TFs classified into 57 families (Fig. 5A). Among the up-regulated genes in iFIB relative to APEX, the lists of MYB-related, bZIP and Dof-protein genes were the most extensive, while among the down-regulated genes of TFs, those for the WRKY, AP2/ERF, MYB, GATA, C2H2, bHLH, GRF, NAC, B3, Trihelix, ZF-HD, NF-YC, NF-YB, TCP, YABBY and GeBP families were the most abundant (Fig. 5A). Among the 217 TF genes with a stage-specific character of expression in the intrusively growing fibers, the most represented families were WRKY, C2H2, ERF, MYB, bHLH, GRF, bZIP, Myb-related, NAC and Dof (Fig. 5, see Supplementary Table S2).

The intrusively growing fibers were enriched in transcripts for several genes encoding MYB-related family members. Among those, three genes, Lus10030184, Lus10012602 and Lus10005134, are homologous to At1g01060 (late elongated hypocotyl 1, LHY1), while Lus10031322 and Lus10019902 are homologues of At3g09600 (REVEILLE 8, RVE8) and At2g18330 (RVE7, EPR1). LHY1, RVE8 and RVE7 are known to function as major regulators of the circadian rhythm in plants^44,45,46. Up-regulation of these genes in the intrusively growing fibers may thus indicate the importance of the daily cycle in the process of intrusive elongation. In addition, two more up-regulated genes, Lus10028306 and Lus10041752, homologous to At1g75250 (RAD-like 6, RSM3), were identified (Fig. 5B). RAD-like 6 was also up-regulated in the flax¹² and hemp⁴⁷ stem portions with developing fibers.

The Dof TFs are a group of plant-specific transcription factors that participate in the control of vascular development and functioning⁴⁸. Transcripts of the flax orthologs of AtDof2.3 (Lus10013905, Lus10001896), AtDof3.3 (Lus10019093), and AtDof3.6 (Lus10005956, Lus10016566) were more abundant in iFIB than in APEX (Fig. 5B). AtDof3.6 is supposed to be involved in the control of cell cycle in cambium⁴⁸. AtDof3.6 gene, accumulating in the vascular tissues during the light period, can act as a positive regulator of the phytochrome B signaling pathway and negatively regulate the cryptochrome-mediated seedling photomorphogenesis⁴⁸.

Several families of TFs containing homeodomains, such as HD-ZIP, HB-other, and the closely related TALE, had the up-regulated in iFIB representatives (Fig. 5A,B). Individual members of these families are believed to control the development of the apical meristem and vascular bundles⁴⁹ and to participate in balancing the carbohydrate supply in carbohydrate-consuming tissues⁵⁰. Lus10025232 is a homologue of Arabidopsis HomeoBox 1 (AtHB1) that regulates genes involved in cell wall formation and acts downstream of PIF1 (PHYTOCHROME-INTERACTING FACTOR 1) to promote hypocotyl elongation, especially in response to short-day photoperiods⁵¹.

The transcription factor bZIP34 (Lus10033630), which was found to be specifically up-regulated at the stage of intrusive growth of flax fibers (Fig. 5B), has a close homology to bZIP61 of Populus (Potri.019G091900, 56.3%), which was up-regulated in the mixed cell-type samples of aspen wood that contain intrusively growing xylem fibers¹¹, [cluster e1]. The expression of these transcription factors was down-regulated in Arabidopsis seedlings at boron deficiency, which is coupled to changes in cell wall properties⁵².

One of the largest families of transcriptional regulators in plants is WRKY, and its members modulate many plant processes⁵³. However, these TFs were rather suppressed in the intrusively growing fibers. Only one gene encoding a member of this family (Lus10019898), which is homologous to ATWRKY2, was specifically up-regulated in iFIB (Fig. 5B).

The regulation of gene expression by TFs is closely interconnected with the action of hormones. Major changes in the expression of hormone-related genes of the intrusively growing fibers occurred with the auxin-related ones (Fig. 6A, see Supplementary Table S2). The list of specifically up-regulated genes contains several genes of auxin-induced proteins that belong to Aux/IAA transcriptional regulators (Fig. 6A). Such proteins do not have a DNA-binding domain, but they do interact with the transcription factors, modulating their activity⁵⁴. Other members of the auxin-signaling pathways with an up-regulated expression in the intrusively growing fibers were small auxin up-regulated RNAs (SAURs) (Fig. 6A). SAURs encode small, plant-specific proteins with no catalytic domains that mediate the effect of auxin in various processes⁵⁵. Notably, the genes for the major receptor of auxin, TIR1 (Lus10011262, Lus10018425), and for auxin signaling F-box 3 (AFB3) protein—other auxin receptor that also mediates Aux/IAA proteins proteasomal degradation and auxin-regulated transcription⁵⁶—were up-regulated in iFIB when compared with APEX (Fig. 6A).

Genes for several auxin transporters—PIN4 (Lus10020830) responsible for the efflux of auxin from a cell, AUX1 (Lus10040212) involved in auxin influx, and auxin transporter-like protein LAX3 (Lus10028078, Lus10025628) that also performs auxin influx—were activated in iFIB when compared with APEX (Fig. 6A and Supplementary Table S2). Genes for calossin-like protein required for polar auxin transport⁵⁷—auxin transport protein BIG (Lus10000087, Lus10005090)—were specifically up-regulated at the stage of fiber intrusive elongation (see Supplementary Table S2). The above-mentioned transporters have different modes of regulation, in combination they provide for the polar transport of auxin⁵⁸. Besides, genes for a plant-specific protein recently identified as vacuolar auxin transporter important for the subcellular distribution of auxin (walls are thin 1, WAT1)⁵⁹ (Lus10010647), and for several WAT1-related proteins previously annotated as nodulin MtN21/EamA-like transporter family protein²⁵ (Lus10024424, Lus10003738, Lus10005965, Lus10019970) were up-regulated in iFIB (Fig. 6B).

Genes for two types of proteins that were included by MapMan in the bin of hormonal regulation as auxin induced-regulated-responsive-activated—an aluminum-induced protein with YGL and LRDR motifs, and a dormancy/auxin associated family protein—were highly and specifically expressed in iFIB (Fig. 6B). The same is true for genes (Lus10019299 and Lus10011523) encoding BTB and TAZ domain protein 4 that according to SUBA4 are involved in the response to auxin and gibberellin stimulus. The exact functions of all these proteins have not been elucidated yet.

Other up-regulated genes for hormone-related proteins included four genes (Lus10010167, Lus10010167, Lus10013445, Lus10040999) homologous to the same Arabidopsis gene (At5g08350) that encodes the GRAM domain-containing protein from the abscisic acid-responsive protein family, and genes for plasma membrane localized ABC transporter involved in abscisic acid transport and responses (Lus10006046, Lus10028704) (Fig. 6C).

Reaction to abiotic stress factors

Splitting of the middle lamellae, and squeezing in between neighbouring cells by pushing them aside, would be expected to cause interrelated mechanical, osmotic, and oxidative stresses. As expected, genes of several members of antioxidant and other stress-related systems constituted a large proportion of those genes specifically expressed by the intrusively growing fibers.

A number of significantly activated genes were those for proteins involved in redox regulation, particularly for several types of thioredoxins (Fig. 7A). Thioredoxins are small oxidoreductases that cleave the disulfide bond between two cysteine residues and are involved in the regulation of the cell redox environment⁶⁰. Lus10007630 and Lus10018383 that were highly activated in iFIB as compared to APEX encoded the specific type of thioredoxin, in which the second cysteine of the redox site was replaced by a serine (thioredoxin-like protein CXXS1); such a monocysteinic thioredoxin has a low disulfide reductase but it has efficient disulfide isomerase activity⁶¹. None of the genes for thioredoxin reductases, another component of the thioredoxin system, were activated in iFIB (see Supplementary Table S1). Lus10030070 encoding lactoylglutathione lyase/glyoxalase I family protein—the enzyme that uses the reduced glutathione as cofactor for detoxification of methylglyoxal and is associated with stress responses⁶²—was highly and specifically activated in iFIB (Fig. 7A); the other member of the same gene family was specifically down-regulated (see Supplementary Table S2).

Many of the proteins encoded by differentially expressed genes are located within the chloroplast, according to the SUBA4 database (http://suba.live/suba-app/factsheet.html)⁶³. Transcription of the genes for the chloroplastic atypical CYS HIS rich thioredoxin 4 (Lus10003437 and Lus10026878), which contributes to the fine-tuning of photosynthetic processes⁶⁴, was significantly activated (Fig. 7A). At the same time, genes for the chloroplast-located glutaredoxins (Lus10011915, Lus10013089, and Lus10022844), glutathione peroxidase (Lus10042418), and glutathione S-transferases (Lus10026643, Lus10040347) were specifically down-regulated in iFIB (see Supplementary Table S2), indicating adjustments of redox regulation in the course of full establishment of photosynthetic machinery. mRNAs of genes for the chloroplast-located chaperone, DnaJ, also known as Hsp40 (Lus10025060, Lus10034484, Lus10017263, Lus10024067, Lus10041671), were specifically abundant in iFIB samples (Fig. 7A).

A significant and specific up-regulation was detected for Lus10038265, Lus10011137, Lus10012939, and Lus10034985 (Fig. 7A, see Supplementary Table S2), all of these are homologous to At5g56550, encoding oxidative stress 3 (OXS3) protein, which is known to be involved in plant responses to heavy metals and oxidative stress⁶⁵. Small, (19 kDa) integral membrane protein from AWPM-19-like family was suggested as a suitable molecular marker of osmotic stress and other abiotic stresses⁶⁶. Several genes for members of this plant-specific family (e.g. Lus10001116, Lus10041238, and Lus10011973) were significantly up-regulated in the intrusively growing fibers (Fig. 7A and Supplementary Table S1). The possibility of osmotic stress in the intrusively elongating fibers is also highlighted by the rather specific and considerable up-regulation of genes for an early responsive to dehydration protein 15 (ERD15) (Lus10009847, Lus10040964), much like the genes for members of the rare-cold inducible (RCI) protein family Lus10014028 and Lus10019890 (Fig. 7 and also Supplementary Table S1) which homolog RCI2A was also induced by dehydration and by salt stress⁶⁷. The genes for two types of small cysteine-rich proteins were up-regulated in iFIB: low-molecular-weight cysteine-rich 68 (LCR68) (Lus10004289, Lus10019227) protein that belongs to plant defensins, and cystatin B (Lus10026117)—proteinase inhibitor that can be induced by multiple abiotic stresses⁶⁸ (Fig. 7A).

The detected pronounced changes in the transcription of genes for mitochondrial proteins included the specific up-regulation of Lus10027755 and Lus10035538, which encode the uncoupling mitochondrial protein involved in protecting plant cells against oxidative stress damage. Besides, genes for the members of mitochondrial electron transfer chain—cytochrome c oxidase (Lus10004084 and Lus10009204), and NADH dehydrogenase (Lus10000460, Lus10009720 and Lus10004838)—were highly activated at the stage of fiber intrusive elongation (see Supplementary Tables S1 and S2). Genes for poorly characterized proteins annotated as wound-responsive (Lus10031617, Lus10033728 and Lus10018529), which Arabidopsis homologues, according to SUBA4⁶³, encode proteins localized in mitochondria, were substantially up-regulated in iFIB (Fig. 7A).

Several genes for the so-called metallochaperones, belonging to the heavy metal transport/detoxification superfamily (e.g. Lus10032924, Lus10015583, and Lus10017730), were also up-regulated in the intrusively growing fibers (Fig. 7A). Cu-binding proteins, which genes (Lus10011158, Lus10022800, Lus10043063, and Lus10005229) were highly activated in iFIB (Fig. 7A), had plastocyanin-like domain and belonged to the family of cupredoxins—small proteins containing type-I copper centres and functioning as electron transfer shuttles between proteins⁶⁹. The significant changes were found in the transcription of genes encoding copper transporters (COPT) COPT1 (Lus10017204) and COPT4 (Lus10021108) (see Supplementary Table S1).

Other genes specifically expressed in the intrusively growing fibers

Several genes specifically expressed in the intrusively growing fibers encode proteins of the cytoskeleton. Their function in these fibers may be related to the intensive cell elongation and to the development of multinuclearity: phragmoplast is not formed and karyokinesis is not coupled to cytokinesis¹³. Genes for actin 7 (Lus10005819, Lus10006784) and for actin 11 (Lus10004169), same as for beta-6 tubulin (Lus10008528) (Fig. 7B, see Supplementary Table S2) were specifically up-regulated in the intrusively growing fibers; the abundance of corresponding mRNAs diminished at a more advanced stage of fiber specialization. Similar dynamics of expression was observed for genes of the kinesin-like protein FRA1 (Lus10010047, Lus10018731, and Lus10024805) from the P-loop containing nucleoside triphosphate hydrolases superfamily. FRA1 contains the N-terminal microtubule binding motor domain and is proposed to contribute to transportation of Golgi-derived vesicles along cortical microtubules to secrete cell wall material⁷⁰. The specifically up-regulated in iFIB Lus10022839 (see Supplementary Table S2) encodes CLIP-associated protein—a microtubule-associated protein that is involved in both cell division and cell expansion. It likely promotes microtubule stability⁷¹. Genes Lus10031066 and Lus10035450 for actin-binding protein with high homology to animal villin were also specifically activated in intrusively growing fibers (Fig. 7B). Such proteins are involved in actin filaments bundling in a calcium-dependent manner⁷². Besides, genes of profilin 5 (Lus10006846, Lus10037591) were up-regulated in iFIB as compared to APEX (Fig. 7B). Profilins are considered to be a part of the molecular network that coordinates the actions of microtubules and actin filaments⁷³. Interestingly, the genes for the CTP synthase family protein (Lus10038403 and Lus10001218) were specifically up-regulated in iFIB (Fig. 7B), same as in sections of aspen wood that contained intrusively growing xylem fibers¹¹, [cluster e1]. Such proteins, independently from their role in biosynthesis of pyrimidine nucleotide, are able to polymerize with formation of filaments, called cytoophidia, that are involved in endoreduplication processes in many living organisms⁷⁴.

Genes for several proteins involved in lipid metabolism were also highly and specifically expressed in iFIB (see Supplementary Table S2). Two of them—Lus10024690, and Lus10032316—are predicted to encode the secreted proteins from the GDSL-like Lipase/Acylhydrolase superfamily. Two other considerably up-regulated genes, Lus10010111 and Lus10012615, are annotated as cytochrome B5 isoform B, which is located in the endoplasmic reticulum and takes part in fatty acid metabolism⁷⁵ (Fig. 7B). Besides, Lus10042398 annotated as polyketide cyclase/dehydrase and lipid transport superfamily protein was specifically up-regulated in iFIB sample; its aspen homolog Potri.004G231000 was assigned to gene cluster that describes the mixed cell-type samples of aspen wood containing intrusively growing xylem fibers¹¹, [cluster e1].

A pronounced feature of the transcriptome landscape in the studied intrusively growing fibers was an increased abundance of transcripts for several plant-specific proteins with poorly characterized functions. Together with the mentioned above GRAM domain-containing protein and plant defensin, remorins can be considered as such. They are vascular plant-specific hydrophilic proteins localized in the plasmodesmata areas and in membrane rafts⁷⁶. Remorin was identified as a protein specifically phosphorylated in the presence of oligogalacturonides⁷⁷. Several remorin genes—Lus10027460, Lus10014785, and Lus10030379—were specifically up-regulated in iFIB (Fig. 7B), same as Lus10020498 for major latex protein 28 (MLP28) (Fig. 7B). Additionally, two genes (Lus10027423, Lus10005769) for phloem protein 2 (PP2), described as one of the most abundant and enigmatic proteins in the phloem⁷⁸, were considerably and rather specifically up-regulated in iFIB (Fig. 7B).

The set of genes for cytochromes p450 (CYP) was also highly activated in iFIB as compared to APEX, and kept the similar level of high transcription at later stages of fiber development (Fig. 7B). CYPs comprise the largest enzymatic protein family in plants; various CYPs function most commonly as monooxygenases that use different substrates to provide a wide range of secondary metabolites⁷⁹.

The genes for several ribosomal proteins—the constituents of cytosolic, chloroplastic and mitochondrial ribosomes—were highly specifically up-regulated in intrusively growing fibers (Fig. 7B); many of the genes for ribosomal proteins were down-regulated in iFIB when compared to APEX (see Supplementary Table S2). Changes in the set of ribosomal proteins may influence the translation machinery; however, understating the function of certain ribosomal protein is still rather vague⁸⁰.

Down-regulated genes in the intrusively growing fibers

In total, the expression of 3,145 genes was suppressed in the intrusively growing fibers, when compared with APEX (log₂FC < −1.6 for iFIB/APEX). The most significantly down-regulated genes included those annotated as cyclins, histones, pathogenesis-related thaumatin, plasmodesmata callose-binding protein 3, as well as some other genes (see Supplementary Table S1). Expression of most of them is the characteristic of apical meristems where they are required for normal cell cycle progression and for meristem organization. Correspondingly, the expression of these genes was detected mainly in the APEX sample: the mRNA abundance of several tens of genes for the histones H2A, H2B, H3, and H4, as well as of 14 genes encoding different classes of cyclins, was reduced in the intrusively growing fibers and was not resumed at the later stages of fiber development. A specific down-regulation at the stage of fiber intrusive elongation was detected for 1,171 genes (log₂FC < −1.6 for iFIB/APEX and for iFIB/FIB) (see Supplementary Table S2) that belonged to various functional categories. The most important down-regulated genes are presented in Fig. 8.

Conclusions

The intrusive growth is performed by relatively few cells, which are surrounded by other cells that may grow symplastically, divide, or differentiate, making it very hard to obtain mRNA samples exclusively from the intrusively growing cells. However, by applying a laser microdissection technique, we could progress in filling the knowledge gaps in the molecular mechanisms that rule this peculiar type of plant cell growth. Having in mind that changes in mRNA levels of individual genes are not automatically related to protein synthesis or, especially, enzyme activity, we used the transcriptome profiling as the very indicative approach to realize the patterns of cell metabolism and to identify the key molecular players in the fiber intrusive growth. Based on our obtained data, fiber physiology may be characterized by several major processes that occur simultaneously. Cell elongation goes together with the establishment of active photosynthetic machinery; the latter is coupled to intensive expression of chloroplast-localized chaperonins and thioredoxins. The high proportion of specifically up-regulated genes encodes proteins localized in the cell wall and in the plasma membrane. These include expansins, enzymes for the modification of pectins (pectate-lyases rather than polygalacturonases and glycanases that can cleave off RG-I side chains) and several proteins with poorly characterized functions.

Symplasmic isolation of the fibers within the other stem tissues involves callose and callose-binding proteins, while the intensive water influx to the plant cell is provided by aquaporins (PIPs and TIPs). Our results suggest that low-molecular sugars rather than inorganic ions are the major osmolytes involved in providing turgor essential for fiber growth. The importance of auxin in fiber development is evidenced by the up-regulation of the genes for receptors and transporters of this hormone and for numerous attenuators of its action, such as the SAURs and AUX/IAAs. A specific set of transcription factors is expressed during intrusive fiber elongation. However, to identify the specific roles of each of these transcription factors in inducing and regulating the intrusive elongation per se, will require much further research. The most suitable candidates are those TFs revealed in the current study with a stage-specific character of expression.

The major molecular players important for the fiber intrusive elongation as revealed at the transcriptome level are summarized in Fig. 8. Their further characterization would help not only to advance our theoretical understanding of this important aspect of plant biology, but also as tools for breeding for higher yields and quality of fiber crops.

Materials and Methods

Plant material and collection of stem samples

Flax plants (Linum usitatissimum L., cultivar Mogilevsky from the collection of the All-Russian Flax Research Institute, Torzhok) were grown under open air in boxes with a 50-cm soil layer and received natural daylight and daily watering. Samples for analysis were collected at the fast-growth period of plant development when the stem height was 25–28 cm (30 days after sowing). A set of samples consisted of the manually-cut off apical part of the stem (APEX)—2 mm of the uppermost stem together with leaf primordia—and the intrusively growing phloem fibers (iFIB) that were isolated from 3rd cm from stem apex by laser microdissection from the stem portion above the snap point¹⁴.

Cryosectioning and laser microdissection

Figure 1 shows the sequential steps taken to obtain fibers at the stage of the intrusive elongation for the RNA-Seq analysis. Flax stem pieces that had the intrusively growing fibers were cut off from the stem with a razor blade, then immediately frozen in liquid nitrogen and stored at −80 °C. Cryosectioning and isolation of the fiber bundles by laser microdissection followed the protocol of Abbott et al.⁸¹ but modified for flax stem tissues. For cryosectioning, the stem pieces were transferred to a cryostat chamber (CM3050, Leica Microsystems, Wetzlar, Germany) and placed atop a drop of Tissue-Tek O.C.T. Compound (Sakura Finetek Inc., CA, USA). Longitudinal stem cryosections (60-μm thickness) were taken at an object temperature of −20 °C, transferred to a pool of cold 100% ethanol (for the RNA isolation) or deionized water (for the ion content analysis) on a POL-membrane frame slide (Leica Microsystems), and allowed to dry for an hour inside the cryostat chamber at −20 °C.

Longitudinal stem sections were analyzed by using a laser microdissection microscope (LMD7000, Leica Microsystems). Identification of the intrusively growing fibers was based on their localization in the stem (between the xylem and 5–6 cell layers of cortex parenchyma and epidermis), elongated cell shape, and their occurrence in bundles. The microdissections were made at a magnification of 10×, a laser power of 35, and a laser speed of 4; they were collected into the caps of 0.2-ml PCR tubes containing 20 μl of RNA lysis solution (RNAqueous-Micro RNA Isolation Kit, Ambion (Austin, TX, USA)) or deionized water (for the ion content analysis) and stored at −80 °C. The final sample available for the RNA analysis contained about 500 microdissections obtained from 30 plants; for the ion content analysis, it consisted of 30–50 microdissections from 5–7 plants. Ion content was determined in two biological replicates, each in two technical repeats. RNA-Seq analysis was performed for two independent biological replicates.

RNA extraction and sequencing

Microdissections from the PCR tube with RNA lysis solution were transferred into a microcentrifuge tube with 300 μl of the same lysis buffer for a silica column-based purification. Elution of total RNA was performed with 2 × 10 μl elution buffer preheated to 95 °C. Total RNA from the APEX (20 plants per sample) was isolated using Trizol extraction method combined with RNeasy Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. For all samples, any residual DNA in them was eliminated with a DNA-free kit (Ambion). RNA quantity and quality were analyzed by a Qubit fluorimeter (Invitrogen, Carlsbad, CA, USA) and Agilent 2100 Bioanalyzer (Agilent Technologies). The RNA integrity index (RIN) was more than 6.2, which corresponds to the quality of RNA for further analysis by methods NGS.

For sequencing, total RNA from the APEX (1–3 µg) sample was processed by using TruSeq Sample Prep Kit (Illumina, San Diego, CA, USA) according to the manufacturer’s instruction. To produce normalized cDNA libraries, followed by sequencing using Illumina MiSeq with single-end 75 bp reads. cDNA libraries from total RNA of iFIB samples (up to 1 µg) were prepared with NEBNext Ultra II Directional RNA Library Prep Kit after selective depletion of ribosomal RNA using RiboMinus™ Plant Kit for RNA-Seq according to the manufacturer’s instructions. Sequence was performed on HiSeq2500 with 60 bp single-end reads. Data in the form of raw reads and sample preparation description were deposited in the Sequence Read Archive (SRA) and are available under the Accession Number as the following accessions: Lus_iFIB_1rep: SRR7284881; Lus_iFIB_2rep: SRR7284883; Lus_APEX: SRR7284882; link to BioProject: PRJNA475325 (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA475325).

RNA-Seq data analysis

The obtained sequence reads (i.e., FASTQ files) were filtered using BBDuk utility of BBMap tools (https://sourceforge.net/projects/bbmap). The clean reads for each sample were mapped onto the flax genome sequence scaffolds⁸² by using HISAT2 v2.1.0⁸³ with default parameters using option–dta-cufflinks and then processed by the software Cufflinks⁸⁴. The annotation of flax genes was taken from Phytozome (phytozome.jgi.doe.gov; downloaded file: Lusitatissimum_200_v1.0.readme.txt). Additionally to 43,484 protein-coding genes in the whole-genome assembly of flax⁸² two CESA7 genes⁸⁵ were integrated into the existing annotation and numbered as Lus10043485 and Lus10043486. Normalized expression values of genes were estimated as the fragments per kilobase of exon per million mapped fragments (FPKM), as determined by Cufflinks (v2.2.1) and according to the published protocol⁸⁴. Data from the RNA-Seq analysis were interpreted with the help of MapMan v3.6¹⁷. and PageMan software¹⁸ modified for flax plants. To avoid division by zero, a pseudocount of 1 was added when calculating the fold changes (FC). To distinguish the genes that were differentially expressed in the intrusively growing fibers, as compared with more advanced stage of fiber specialization, the data obtained in this study were compared with the previously published results of the transcriptome analysis for fibers isolated during the stage of cell wall thickening¹² (sample ‘FIB’).

Validation of transcriptome experiments using qRT–PCR

For quantitative reverse transcription PCR (qRT–PCR), the total RNA (about 0.8 µg) was converted to cDNA with M-MuLV Reverse Transcriptase (Fermentas) using oligo (dT)18 primers according to the standard protocol of manufacturer. Twelve genes identified by RNA-Seq were assayed by qRT–PCR (see Supplementary Table S4) using CFX96 Touch Real-Time PCR Detection System (Bio-Rad). The thermal cycling conditions were 95 °C for 5 min, 40 cycles at 95 °C for 15 s, and 60 °C for1 min. A 60–95 °C melting curve was performed to confirm specificity of the products. From each of two biologically independent cDNA samples two independent technical replications were performed and averaged for further calculations. Relative transcript abundance calculations were performed using the ΔΔCt method⁸⁶ with APEX sample as a reference. The genes encoding basic leucine-zipper 44 (Lus10001347), homeobox 1 (Lus10008795), AGAMOUS-like 19 (Lus10006715), RAD-like 6 (Lus10028306), proline-rich protein 2 (Lus10042392), indoleacetic acid-induced protein (Lus10018765), pectin lyase-like (Lus10011885 and Lus10022817), DUF23 (Lus10038387), DUF642 (Lus10008080), cellulose synthases (Lus10007538 and Lus10012198) were selected for verification by the qRT-PCR. The genes of eukaryotic translation initiation factors 1 A, 5 A (LusETIF1, LusETIF5A) and glyceraldehyde-3-phosphate dehydrogenase (LusGAPDH) were used as the housekeeping genes⁸⁷.

Ion content analysis

For this analysis, four samples types were collected by microdissection: the intrusively growing fibers, the fibers at the cell wall thickening stage, the cortex parenchyma cells above and below the snap point. Samples were ground by using a glass pestle in deionized water. The qualitative and quantitative analyses of the ions were performed, respectively, by using an ICP-MS ELAN 9000 (Perkin Elmer, Shelton, CT, USA) and an Optima 8300 ICP-OES Spectrometer (Perkin Elmer). Argon flow rate through the outer burner tube: 15 l/min, the middle burner: 0.5 l/min and the sprayer: 0.55 l/min; the power for the plasma operation: 1400 V. The solution flow rate used during the samples’ analysis: 1.5 ml/min; washing time between the samples: 30 s; the time to sample nebulization: 30 s. The samples were analyzed at the following wavelengths: for Mg²⁺: 280.272 nm, Ca²⁺: 396.846 nm, Na⁺: 589.004 nm, and K⁺:766.493 nm. The analysis of these ions was done with two biological replicates; each analyzed twice; the results are conveyed as means ± standard deviation.

Data Availability

All data generated or analysed during this study are included in this published article (and its Supplementary Information files).

References

Wasteneys, G. O. & Galway, M. E. Remodeling the cytoskeleton for growth and form: an overview with some new views. Annu. Rev. Plant Biol. 54, 691–722 (2003).
Article CAS Google Scholar
Snegireva, A. V. et al. Intrusive growth of sclerenchyma fibres. Russ. J. Plant Physiol. 57, 342–355 (2010).
Article CAS Google Scholar
Guerriero, G., Hausman, J. F. & Cai, G. No stress! Relax! Mechanisms governing growth and shape in plant cells. Int. J. Mol. Sci. 15, 5094–5114 (2014).
Article Google Scholar
Esau, K. Plant anatomy, 2nd edn (John Wiley & Sons, New York, 1965).
Fahn, A. Plant anatomy, 4th edn (Pergamon Press, Oxford, 1990).
Lev-Yadun, S. Intrusive growth—the plant analog of dendrite and axon growth in animals. New Phytol. 150, 508–512 (2001).
Article Google Scholar
Gorshkova, T. et al. Plant fiber formation: state of the art, recent and expected progress, and open questions. Critic. Rev. Plant Sci. 31, 201–228 (2012).
Article CAS Google Scholar
Nezhad, A. S. & Geitmann, A. The cellular mechanics of an invasive lifestyle. J. Exp. Bot. 64, 4709–4728 (2013).
Article Google Scholar
Lev-Yadun, S. Plant development: Cell movement relative to each other is both common and very important. Plant Signaling & Behav. 10, e991566, https://doi.org/10.4161/15592324.2014.991566 (2015).
Article CAS Google Scholar
Snegireva, A., Chernova, T., Ageeva, M., Lev-Yadun, S. & Gorshkova, T. Intrusive growth of primary and secondary phloem fibres in hemp stem determines fibre-bundle formation and structure. AoB Plants. 7, plv061, https://doi.org/10.1093/aobpla/plv061 (2015).
Article CAS PubMed PubMed Central Google Scholar
Sundell, D. et al. 2017. AspWood: high-spatial-resolution transcriptome profiles reveal uncharacterized modularity of wood formation in Populus tremula. Plant Cell 29(1604), 1585 (2017).
Article CAS Google Scholar
Gorshkov, O. et al. Transcriptome portrait of cellulose-enriched flax fibres at advanced stage of specialization. Plant Mol. Biol. 93, 431–449 (2017).
Article CAS Google Scholar
Ageeva, M. V. et al. Intrusive growth of flax phloem fibres is of intercalary type. Planta 222, 565–574 (2005).
Article CAS Google Scholar
Gorshkova, T. A., Sal’nikov, V. V., Chemikosova, S. B., Ageeva, M. V. & Pavlencheva, N. V. The snap point: A transition point in Linum usitatissimum bast fibre development. Ind. Crops and Prods. 18, 213–221 (2003).
Article Google Scholar
Zhang, N. & Deyholos, M. K. RNASeq analysis of the shoot apex of flax (Linum usitatissimum) to identify phloem fiber specification genes. Front. Plant Sci. 7, 950, https://doi.org/10.3389/fpls.2016.00950 (2016).
Article CAS PubMed PubMed Central Google Scholar
SEQC/MAQC-III Consortium. A comprehensive assessment of RNA-seq accuracy, reproducibility and information content by the Sequencing Quality Control Consortium. Nat. Biotechnol. 32, 903–914 (2014).
Thimm, O. et al. MAPMAN: a user-driven tool to display genomics data sets onto diagrams of metabolic pathways and other biological processes. Plant J. 37, 914–939 (2004).
Article CAS Google Scholar
Usadel, B. et al. PageMan: an interactive ontology tool to generate, display, and annotate overview graphs for profiling experiments. BMC Bioinformatics 18, 535, https://doi.org/10.1186/1471-2105-7-535 (2006).
Article CAS Google Scholar
Cosgrove, D. J. Growth of the plant cell wall. Nature Reviews Mol. Cell Biol. 6, 850–861 (2005).
Article CAS Google Scholar
Fricke, W. & Knipfer, T. P Aquaporins and Cell Elongation. In: Plant Aquaporins. (ed. Chaumont, F. & Tyerman, S. D.) 107–131 (Springer, Heidelberg, 2017).
Aslam, U., Khatoon, A., Cheema, H. M. N. & Bashir, A. Identification and characterization of plasma membrane aquaporins isolated from fiber cells of Calotropis procera. J. Zhejiang Univ. Sci. B. 14, 586–595 (2013).
Article CAS Google Scholar
Hukin, D., Doering-Saad, C., Thomas, C. & Pritchard, J. Sensitivity of cell hydraulic conductivity to mercury is coincident with symplasmic isolation and expression of plasmalemma aquaporin genes in growing maize roots. Planta 215, 1047–1056 (2002).
Article CAS Google Scholar
Muto, Y. et al. Vacuolar proton pumps and aquaporins involved in rapid internode elongation of deepwater rice. Biosci., Biotech. Biochem. 75, 114–122 (2011).
Article CAS Google Scholar
Shivaraj, S. M. et al. Genome-wide identification, characterization, and expression profile of aquaporin gene family in flax (Linum usitatissimum). Sci. Rep. 7, 46137, https://doi.org/10.1038/srep46137 (2017).
Article ADS CAS PubMed PubMed Central Google Scholar
Denancé, N., Szurek, B. & Noël, L. D. Emerging functions of nodulin-like proteins in non-nodulating plant species. Plant Cell Physiol. 55, 469–474 (2014).
Article Google Scholar
Büttner, M. The monosaccharide transporter(-like) gene family in Arabidopsis. FEBS Lett. 581, 2318–2324 (2007).
Article Google Scholar
Schneider, S., Beyhl, D., Hedrich, R. & Sauer, N. Functional and physiological characterization of Arabidopsis INOSITOL TRANSPORTER1, a novel tonoplast-localized transporter for myo-inositol. Plant Cell 20, 1073–1087 (2008).
Article CAS Google Scholar
Gillaspy, G. E. The cellular language of myo-inositol signaling. New Phytol. 192, 823–839 (2011).
Article CAS Google Scholar
Huang, J., Ghosh, R. & Bankaitis, V. A. Sec14-like Phosphatidylinositol transfer proteins and the biological landscape of phosphoinositide signaling in plants. Biochim. Biophys. Acta 1861(9PtB), 1352–1364 (2016).
Article CAS Google Scholar
Rinne, P. L. & van der Schoot, C. Symplasmic fields in the tunica of the shoot apical meristem coordinate morphogenetic events. Development 125, 1477–1485 (1998).
CAS PubMed Google Scholar
Chen, X. Y. & Kim, J. Y. Callose synthesis in higher plants. Plant Signaling & Behav. 4, 489–492 (2009).
Article CAS Google Scholar
Ruan, Y. L., Llewellyn, D. J. & Furbank, R. T. The control of single-celled cotton fiber elongation by developmentally reversible gating of plasmodesmata and coordinated expression of sucrose and K+ transporters and expansin. Plant Cell 13, 47–60 (2001).
CAS PubMed PubMed Central Google Scholar
Barral, P. et al. An olive pollen protein with allergenic activity, Ole e 10, defines a novel family of carbohydrate-binding modules and is potentially implicated in pollen germination. Biochem. J. 390, 77–84 (2005).
Article CAS Google Scholar
Simpson, C., Thomas, C., Findlay, K., Bayer, E. & Maule, A. J. An Arabidopsis GPI-anchor plasmodesmal neck protein with callose binding activity and potential to regulate cell-to-cell trafficking. Plant Cell 21, 581–594 (2009).
Article CAS Google Scholar
Chary, S. N., Hicks, G. R., Choi, Y. G., Carter, D. & Raikhel, N. V. Trehalose-6-Phosphate Synthase/Phosphatase regulates cell shape and plant architecture in Arabidopsis. Plant Physiol. 146, 97–107 (2008).
Article CAS Google Scholar
Carpita, N. C. & Gibeaut, D. M. Structural models of primary cell walls in flowering plants: consistency of molecular structure with the physical properties of the walls during growth. Plant J. 3, 1–30 (1993).
Article CAS Google Scholar
Sénéchal, F., Wattier, C., Rustérucci, C. & Pelloux, J. Homogalacturonan-modifying enzymes: structure, expression, and roles in plants. J. Exp. Bot. 65, 5125–5160 (2014).
Article Google Scholar
Pinzón-Latorre, D. & Deyholos, M. K. Pectinmethylesterases (PME) and pectinmethylesterase inhibitors (PMEI) enriched during phloem fiber development in flax (Linum usitatissimum). PLoS One 9, e105386, https://doi.org/10.1371/journal.pone.0105386 (2014).
Article ADS CAS PubMed PubMed Central Google Scholar
Philippe, F., Pelloux, J. & Rayon, C. Plant pectin acetylesterase structure and function: new insights from bioinformatic analysis. BMC Genomics 18, 456, https://doi.org/10.1186/s12864-017-3833-0 (2017).
Article PubMed PubMed Central Google Scholar
Rose, J. K., Braam, J., Fry, S. C. & Nishitani, K. The XTH family of enzymes involved in xyloglucan endotransglucosylation and endohydrolysis: current perspectives and a new unifying nomenclature. Plant and Cell Physiol. 43, 1421–1435 (2002).
Article CAS Google Scholar
Van Sandt, V. S., Suslov, D., Verbelen, J. P. & Vissenberg, K. Xyloglucan endotransglucosylase activity loosens a plant cell wall. Ann. Bot. 100, 1467–1473 (2007).
Article Google Scholar
Maris, A., Suslov, D., Fry, S. C., Verbelen, J. P. & Vissenberg, K. Enzymic characterization of two recombinant xyloglucan endotransglucosylase/hydrolase (XTH) proteins of Arabidopsis and their effect on root growth and cell wall extension. J. Exp. Bot. 60, 3959–3972 (2009).
Article CAS Google Scholar
Cosgrove, D. J. Cell Wall Loosening by Expansins. Plant Physiol. 118, 333–339 (1998).
Article CAS Google Scholar
Kamioka, M. et al. Direct repression of evening genes by CIRCADIAN CLOCK-ASSOCIATED1 in the Arabidopsis circadian clock. Plant Cell 28, 696–711 (2016).
Article CAS Google Scholar
Nguyen, N. H. & Lee, H. MYB-related transcription factors function as regulators of the circadian clock and anthocyanin biosynthesis in Arabidopsis. Plant Signaling & Behav. 11, e1139278, https://doi.org/10.1080/15592324.2016.1139278 (2016).
Article CAS Google Scholar
Graf, A.et al. Parallel analysis of Arabidopsis circadian clock mutants reveals different scales of transcriptome and proteome regulation. Open Biology 7(3) pii: 160333 10.1098/rsob.160333 (2017).
Article Google Scholar
Guerriero, G. et al. Transcriptomic profiling of hemp bast fibres at different developmental stages. Sci. Rep. 7, 4961, https://doi.org/10.1038/s41598-017-05200-8 (2017).
Article ADS CAS PubMed PubMed Central Google Scholar
Le Hir, R. & Bellini, C. The plant-specific Dof transcription factors family: new players involved in vascular system development and functioning in Arabidopsis. Front. in Plant Sci. 4, 164, https://doi.org/10.3389/fpls.2013.00164 (2013).
Article Google Scholar
Ariel, F. D., Manavella, P. A., Dezar, C. A. & Chan, R. L. The true story of the HD-Zip family. Trends Plant Sci. 12, 419–426 (2007).
Article CAS Google Scholar
Satoh, R., Fujita, Y., Nakashima, K., Shinozaki, K. & Yamaguchi-Shinozaki, K. A novel subgroup of bZIP proteins functions as transcriptional activators in hypoosmolarity-responsive expression of the ProDH gene in Arabidopsis. Plant Cell Physiol. 45, 309–317 (2004).
Article CAS Google Scholar
Capella, M., Ribone, P. A., Arce, A. L. & Chan, R. L. Arabidopsis thaliana HomeoBox 1 (AtHB1), a Homedomain-Leucine Zipper I (HD-Zip I) transcription factor, is regulated by PHYTOCHROME-INTERACTING FACTOR 1 to promote hypocotyl elongation. New Phytol. 207, 669–682 (2015).
Article CAS Google Scholar
González-Fontes, A. et al. Transcription factors as potential participants in the signal transduction pathway of boron deficiency. Plant Signal. Behav. 8, e26114, https://doi.org/10.4161/psb.26114 (2013).
Article CAS PubMed PubMed Central Google Scholar
Zhang, Y. & Wang, L. The WRKY transcription factor superfamily: its origin in eukaryotes and expansion in plants. BMC Evol. Biol. 5, 1, https://doi.org/10.1186/1471-2148-5-1 (2005).
Article CAS PubMed PubMed Central Google Scholar
Eyidogan, F., Oz, M. T., Yucel, M. & Oktem, H. A. Signal transduction of phytohormones under abiotic stresses. In: Phytohormones and Abiotic Stress Tolerance in Plants. (ed. Khan, N. A., Nazar, R., Iqbal, N. & Anjum, N. A.) 1–48 (Springer, Heidelberg, 2012).
Google Scholar
Ren, H. & Gray, W. M. SAUR proteins as effectors of hormonal and environmental signals in plant growth. Mol. Plant 8, 1153–1164 (2015).
Article CAS Google Scholar
Dharmasiri, N. et al. Plant development is regulated by a family of auxin receptor F box proteins. Dev. Cell 9, 109–119 (2005).
Article CAS Google Scholar
Gil, P. et al. BIG: a calossin-like protein required for polar auxin transport in Arabidopsis. Genes Dev. 15, 1985–1997 (2001).
Article CAS Google Scholar
Kleine-Vehn, J., Dhonukshe, P., Swarup, R., Bennett, M. & Frimla, J. Subcellular Trafficking of the Arabidopsis Auxin Influx Carrier AUX1 Uses a Novel Pathway Distinct from PIN. Plant Cell 18, 3171–3181 (2006).
Article CAS Google Scholar
Ranocha, P. et al. Arabidopsis WAT1 is a vacuolar auxin transport facilitator required for auxin homoeostasis. Nat. Commun. 4, 2625, https://doi.org/10.1038/ncomms3625 (2013).
Article CAS PubMed PubMed Central Google Scholar
Gelhaye, E., Rouhier, N., Navrot, N. & Jacquot, J. P. The plant thioredoxin system. Cell. Mol. Life Sci. 62, 24–35 (2005).
Article CAS Google Scholar
Serrato, A. J., Guilleminot, J., Meyer, Y. & Vignols, F. AtCXXS: atypical members of the Arabidopsis thaliana thioredoxin h family with a remarkably high disulfide isomerase activity. Physiol. Plant. 133, 611–622 (2008).
Article CAS Google Scholar
Sankaranarayanan, S. et al. Glyoxalase Goes Green: The Expanding Roles of Glyoxalase in Plants. Int. J. Mol. Sci. 18, 898, https://doi.org/10.3390/ijms18040898 (2017).
Article PubMed Central Google Scholar
Hooper, C. M., Castleden, I., Tanz, S. K., Aryamanesh, N. & Millar, A. H. SUBA4: the interactive data analysis centre for Arabidopsis subcellular protein locations. Nucleic Acids Res. 45(D1), D1064–D1074 (2017).
Article CAS Google Scholar
Nikkanen, L., Toivola, J., Guinea Diaz, M. & Rintamäki, E. Chloroplast thioredoxin systems: prospects for improving photosynthesis. Phil. Trans. R. Soc. Lond. B 372 (1730). pii: 20160474;doi:10.1098/rstb.2016.0474 (2017).
Blanvillain, R., Kim, J. H., Wu, S., Lima, A. & Ow, D. W. OXIDATIVE STRESS 3 is a chromatin-associated factor involved in tolerance to heavy metals and oxidative stress. Plant J. 57, 654–665 (2009).
Article CAS Google Scholar
Wanke, D., Berendzen, K.W., Kilian, J. & Harter, K. Insights into the Arabidopsis Abiotic Stress Response from the AtGenExpress Expression Profile Dataset in Plant Stress Biology: From Genomics to Systems Biology (ed. Hirt, H.) (Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim, Germany, 2009).
Medina, J. et al. Arabidopsis mutants deregulated in RCI2Aexpression reveal new signaling pathways in abiotic stress responses. Plant J. 42, 586–597 (2005).
Article CAS Google Scholar
Zhang, X., Liu, S. & Takano, T. Two cysteine proteinase inhibitors from Arabidopsis thaliana, AtCYSa and AtCYSb, increasing the salt, drought, oxidation and cold tolerance. Plant Mol. Biol. 68, 131–143 (2008).
Article CAS Google Scholar
Choi, M. & Davidson, V. L. Cupredoxins—a study of how proteins may evolve to use metals for bioenergetic processes. Metallomics 3, 140–151 (2011).
Article CAS Google Scholar
Zhu, C. et al. The fragile Fiber1 kinesin contributes to cortical microtubule-mediated trafficking of cell wall components. Plant Physiol. 167, 780–792 (2015).
Article CAS Google Scholar
Ambrose, J. C., Shoji, T., Kotzer, A. M., Pighin, J. A. & Wasteneys, G. O. The Arabidopsis CLASP Gene Encodes a Microtubule-Associated Protein Involved in Cell Expansion and Division. Plant Cell 19, 2763–2775 (2007).
Article CAS Google Scholar
Hepler, P. K. The Cytoskeleton and Its Regulation by Calcium and Protons. Plant Physiol. 170, 3–22 (2016).
Article CAS Google Scholar
Nejedla, M. et al. Profilin connects actin assembly with microtubule dynamics. Mol. Biol. Cell 27, 2381–2393 (2016).
Article CAS Google Scholar
Wang, P.-Y. et al. Regulation of CTP Synthase Filament Formation During DNA Endoreplication in Drosophila. Genetics 201, 1511–1523 (2015).
Article CAS Google Scholar
Kumar, R., Tran, L. S., Neelakandan, A. K. & Nguyen, H. T. Higher plant cytochrome b5 polypeptides modulate fatty acid desaturation. PLoS One. 7, e31370. https://doi.org/10.1371/journal.pone.0031370 (2012).
Article ADS Google Scholar
Raffaele, S. et al. Remorin, a Solanaceae protein resident in membrane rafts and plasmodesmata, impairs Potato virus X movement. Plant Cell 21, 1541–1555 (2009).
Article CAS Google Scholar
Farmer, E. E., Moloshok, T. D., Saxton, M. J. & Ryan, C. A. Oligosaccharide signaling in plants. Specificity of oligouronide-enhanced plasma membrane protein phosphorylation. J. Biol. Chem. 266, 3140–3145 (1991).
CAS PubMed Google Scholar
Dinant, S. et al. Diversity of the superfamily of phloem lectins (phloem protein 2) in angiosperms. Plant Physiol. 131, 114–128 (2003).
Article CAS Google Scholar
Xu, J., Wang, X. & Guo, W. The cytochrome P450 superfamily: key players in plant development and defense. J. Integr. Agric. 14, 1673–1686 (2015).
Article CAS Google Scholar
Carroll, A. J. The Arabidopsis Cytosolic Ribosomal Proteome: From form to Function. Front. Plant Sci. 4, 32, https://doi.org/10.3389/fpls.2013.00032 (2013).
Article PubMed PubMed Central Google Scholar
Abbott, E., Hall, D., Hamberger, B. & Bohlmann, J. Laser microdissection of conifer stem tissues: isolation and analysis of high quality RNA, terpene synthase enzyme activity and terpenoid metabolites from resin ducts and cambial zone tissue of white spruce (Picea glauca). BMC Plant Biol. 10, 106, https://doi.org/10.1186/1471-2229-10-106 (2010).
Article CAS PubMed PubMed Central Google Scholar
Wang, Z. et al. The genome of flax (Linum usitatissimum) assembled de novo from short shotgun sequence reads. Plant J. 72, 461–473 (2012).
Article Google Scholar
Kim, D., Langmead, B. & Salzberg, S. L. HISAT: a fast spliced aligner with low memory requirements. Nat. Methods 12, 357–360 (2015).
Article CAS Google Scholar
Trapnell, C. et al. Differential gene and transcript expression analysis of RNA-Seq experiments with TopHat and Cufflinks. Nature Protocols 7, 562–578 (2012).
Article CAS Google Scholar
Mokshina, N., Gorshkova, T. & Deyholos, M. K. Chitinase-Like (CTL) and cellulose synthase (CESA) gene expression in gelatinous-type cellulosic walls of flax (Linum usitatissimumL.) bast fibers. PLoS One 9, e97949, https://doi.org/10.1371/journal.pone.0097949 (2014).
Article ADS CAS PubMed PubMed Central Google Scholar
Livak, K. J. & Schmittgen, T. D. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods 25, 402–408 (2001).
Article CAS Google Scholar
Huis, R., Hawkins, S. & Neutelings, G. Selection of reference genes for quantitative gene expression normalization in flax (Linum usitatissimum L. ). BMC Plant Biol. 10, 71, https://doi.org/10.1186/1471-2229-10-71 (2010).
Article CAS PubMed PubMed Central Google Scholar

Download references

Acknowledgements

This study was partially supported by the Russian Science Foundation [16–14–10256]. We thank the Research Resource Center “Molecular and Cell Technologies”, of Saint Petersburg State University (Saint Petersburg, Russia), for the access to cryosection and laser microdissection facilities, as well as the Central Scientific Research Institute of Geology of Industrial Minerals (Kazan, Russia) for helping with the ion analysis. We also would like to thank joint Kazan Federal University - RIKEN laboratory (Kazan, Russia) and personally Natalia Gogoleva for the assistance in RNA-Seq analysis.

Author information

Authors and Affiliations

Kazan Institute of Biochemistry and Biophysics, Federal Research Center “Kazan Scientific Center of RAS” 420111, Lobachevsky Str., 2/31, Kazan, Russian Federation
Tatyana Gorshkova, Tatyana Chernova, Natalia Mokshina, Vladimir Gorshkov, Liudmila Kozlova & Oleg Gorshkov

Authors

Tatyana Gorshkova
View author publications
You can also search for this author in PubMed Google Scholar
Tatyana Chernova
View author publications
You can also search for this author in PubMed Google Scholar
Natalia Mokshina
View author publications
You can also search for this author in PubMed Google Scholar
Vladimir Gorshkov
View author publications
You can also search for this author in PubMed Google Scholar
Liudmila Kozlova
View author publications
You can also search for this author in PubMed Google Scholar
Oleg Gorshkov
View author publications
You can also search for this author in PubMed Google Scholar

Contributions

T.G. conceived the study topic, organized, wrote, and reviewed the manuscript. T.Ch. performed cryosectioning and laser microdissection, took part in ion content analysis, prepared Figures 1, 8. N.M. prepared cDNA libraries for transcriptome sequencing, prepared Figures 3–7. V.G. conducted transcriptome sequencing using Illumina platform. L.K. analyzed RNA-Seq data and literature data. O.G. took part in cryosectioning and laser microdissection, isolated RNA, performed all bioinformatics work, prepared Figure 2. All authors reviewed the manuscript.

Corresponding author

Correspondence to Tatyana Gorshkova.

Ethics declarations

Competing Interests

The authors declare no competing interests.

Additional information

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Supplementary Dataset 1

Supplementary Dataset 2

Supplementary Dataset 3

Supplementary Dataset 4

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Cite this article

Gorshkova, T., Chernova, T., Mokshina, N. et al. Transcriptome Analysis of Intrusively Growing Flax Fibers Isolated by Laser Microdissection. Sci Rep 8, 14570 (2018). https://doi.org/10.1038/s41598-018-32869-2

Download citation

Received: 23 March 2018
Accepted: 18 September 2018
Published: 01 October 2018
DOI: https://doi.org/10.1038/s41598-018-32869-2

Keywords

This article is cited by

Genome-wide Identification and Characterization of Expansin Genes in Jute
- Md Sabbir Hossain
- Borhan Ahmed
- Md. Shahidul Islam
Tropical Plant Biology (2022)
Genome-wide identification and expression analysis of the WRKY transcription factor family in flax (Linum usitatissimum L.)
- Hongmei Yuan
- Wendong Guo
- Guangwen Wu
BMC Genomics (2021)
Genome-wide identification of fasciclin-like arabinogalactan proteins in jute and their expression pattern during fiber formation
- Md. Sabbir Hossain
- Borhan Ahmed
- Md. Shahidul Islam
Molecular Biology Reports (2020)
Expression of cellulose synthase-like genes in two phenotypically distinct flax (Linum usitatissimum L.) subspecies
- Dmitry Galinousky
- Tsimafei Padvitski
- Alexandr Kilchevsky
Genetic Resources and Crop Evolution (2020)

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.