Microarray analysis of flax (Linum usitatissimum L.) stems identifies transcripts enriched in fibre-bearing phloem tissues

How to explore what is hidden? A review of techniques for vascular tissue expression profile analysis

The evolution of plants to efficiently transport water and assimilates over long distances is a major evolutionary success that facilitated their growth and colonization of land. Vascular tissues, namely xylem and phloem, are characterized by high specialization, cell heterogeneity, and diverse cell components. During differentiation and maturation, these tissues undergo an irreversible sequence of events, leading to complete protoplast degradation in xylem or partial degradation in phloem, enabling their undisturbed conductive function. Due to the unique nature of vascular tissue, and the poorly understood processes involved in xylem and phloem development, studying the molecular basis of tissue differentiation is challenging. In this review, we focus on methods crucial for gene expression research in conductive tissues, emphasizing the importance of initial anatomical analysis and appropriate material selection. We trace the expansion of molecular techniques in vascular gene expression studies and discuss the application of single-cell RNA sequencing, a high-throughput technique that has revolutionized transcriptomic analysis. We explore how single-cell RNA sequencing will enhance our knowledge of gene expression in conductive tissues.

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.

Drought and Oxidative Stress in Flax (Linum usitatissimum L.) Entails Harnessing Non-Canonical Reference Gene for Precise Quantification of qRT-PCR Gene Expression

Flax (Linum usitatissimum L.) is a self-pollinating, annual, diploid crop grown for multi-utility purposes for its quality oil, shining bast fiber, and industrial solvent. Being a cool (Rabi) season crop, it is affected by unprecedented climatic changes such as high temperature, drought, and associated oxidative stress that, globally, impede its growth, production, and productivity. To precisely assess the imperative changes that are inflicted by drought and associated oxidative stress, gene expression profiling of predominant drought-responsive genes (AREB, DREB/CBF, and ARR) was carried out by qRT-PCR. Nevertheless, for normalization/quantification of data obtained from qRT-PCR results, a stable reference gene is mandatory. Here, we evaluated a panel of four reference genes (Actin, EF1a, ETIF5A, and UBQ) and assessed their suitability as stable reference genes for the normalization of gene expression data obtained during drought-induced oxidative stress in flax. Taking together, from the canonical expression of the proposed reference genes in three different genotypes, we report that EF1a as a stand-alone and EF1a and ETIF5A in tandem are suitable reference genes to be used for the real-time visualization of cellular impact of drought and oxidative stress on flax.

Cell Wall Layer Induced in Xylem Fibers of Flax Upon Gravistimulation Is Similar to Constitutively Formed Cell Walls of Bast Fibers

In the fibers of many plant species after the formation of secondary cell walls, cellulose-enriched cell wall layers (often named G-layers or tertiary cell walls) are deposited which are important in many physiological situations. Flax (Linum usitatissimum L.) phloem fibers constitutively develop tertiary cell walls during normal plant growth. During the gravitropic response after plant inclination, the deposition of a cellulose-enriched cell wall layer is induced in xylem fibers on one side of the stem, providing a system similar to that of tension wood in angiosperm trees. Atomic force microscopy (AFM), immunochemistry, and transcriptomic analyses demonstrated that the G-layer induced in flax xylem fibers was similar to the constitutively formed tertiary cell wall of bast (phloem) fibers but different from the secondary cell wall. The tertiary cell walls, independent of tissue of origin and inducibility, were twice as stiff as the secondary cell walls. In the gravitropic response, the tertiary cell wall deposition rate in xylem was higher than that of the secondary cell wall. Rhamnogalacturonan I (RG-I) with galactan side chains was a prominent component in cellulose-rich layers of both phloem and xylem flax fibers. Transcriptomic events underlying G-layer deposition in phloem and xylem fibers had much in common. At the induction of tertiary cell wall deposition, several genes for rhamnosyltransferases of the GT106 family were activated in xylem samples. The same genes were expressed in the isolated phloem fibers depositing the tertiary cell wall. The comparison of transcriptomes in fibers with both inducible and constitutive tertiary cell wall deposition and xylem tissues that formed the secondary cell walls is an effective system that revealed important molecular players involved in the formation of cellulose-enriched cell walls.

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

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

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

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

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.

Plant Soft Rot Development and Regulation from the Viewpoint of Transcriptomic Profiling

Soft rot caused by Pectobacterium species is a devastating plant disease poorly characterized in terms of host plant responses. In this study, changes in the transcriptome of tobacco plants after infection with Pectobacterium atrosepticum (Pba) were analyzed using RNA-Seq. To draw a comprehensive and nontrivially itemized picture of physiological events in Pba-infected plants and to reveal novel potential molecular "players" in plant-Pba interactions, an original functional gene classification was performed. The classifications present in various databases were merged, enriched by "missed" genes, and divided into subcategories. Particular changes in plant cell wall-related processes, perturbations in hormonal and other regulatory systems, and alterations in primary, secondary, and redox metabolism were elucidated in terms of gene expression. Special attention was paid to the prediction of transcription factors (TFs) involved in the disease's development. Herewith, gene expression was analyzed within the predicted TF regulons assembled at the whole-genome level based on the presence of particular cis-regulatory elements (CREs) in gene promoters. Several TFs, whose regulons were enriched by differentially expressed genes, were considered to be potential master regulators of Pba-induced plant responses. Differential regulation of genes belonging to a particular multigene family and encoding cognate proteins was explained by the presence/absence of the particular CRE in gene promoters.

Flax rhamnogalacturonan lyases: phylogeny, differential expression and modeling of protein structure

Rhamnogalacturonan lyases (RGLs; EC 4.2.2.23) degrade the rhamnogalacturonan I (RG-I) backbone of pectins present in the plant cell wall. These enzymes belong to polysaccharide lyase family 4, members of which are mainly from plants and plant pathogens. RGLs are investigated, as a rule, as pathogen 'weapons' for plant cell wall degradation and subsequent infection. Despite the presence of genes annotated as RGLs in plant genomes and the presence of substrates for enzyme activity in plant cells, evidence supporting the involvement of this enzyme in certain processes is limited. The differential expression of some RGL genes in flax (Linum usitatissimum L.) tissues, revealed in our previous work, prompted us to carry out a total revision (phylogenetic analysis, analysis of expression and protein structure modeling) of all the sequences of flax predicted as coding for RGLs. Comparison of the expressions of LusRGL in various tissues of flax stem revealed that LusRGLs belong to distinct phylogenetic clades, which correspond to two co-expression groups. One of these groups comprised LusRGL6-A and LusRGL6-B genes and was specifically upregulated in flax fibers during deposition of the tertiary cell wall, which has complex RG-I as a key noncellulosic component. The results of homology modeling and docking demonstrated that the topology of the LusRGL6-A catalytic site allowed binding to the RG-I ligand. These findings lead us to suggest the presence of RGL activity in planta and the involvement of special isoforms of RGLs in the modification of RG-I of the tertiary cell wall in plant fibers.

Cell wall composition and transcriptomics in stem tissues of stinging nettle (Urtica dioica L.): Spotlight on a neglected fibre crop

Stinging nettle (Urtica dioica L.) produces silky cellulosic fibres, as well as bioactive molecules. To improve the knowledge on nettle and enhance its opportunities of exploitation, a draft transcriptome of the "clone 13" (a fibre clone) is here presented. The transcriptome of whole internodes sampled at the top and middle of the stem is then compared with the core and cortical tissues sampled at the bottom. Young internodes show an enrichment in genes involved in the biosynthesis of phytohormones (auxins and jasmonic acid) and secondary metabolites (flavonoids). The core of internodes collected at the bottom of the stem is enriched in genes partaking in different aspects of secondary cell wall formation (cellulose, hemicellulose, lignin biosynthesis), while the cortical tissues reveal the presence of a C starvation signal probably due to the UDP-glucose demand necessary for the thickening phase of bast fibres. Cell wall analysis indicates a difference in rhamnogalacturonan structure/composition of mature bast fibres, as evidenced by the higher levels of galactose measured, as well as the occurrence of more water-soluble pectins in elongating internodes. The targeted quantification of phenolics shows that the middle internode and the cortical tissues at the bottom have higher contents than top internodes. Ultrastructural analyses reveal the presence of a gelatinous layer in bast fibres with a lamellar structure. The data presented will be an important resource and reference for future molecular studies on a neglected fibre crop.

Flax (Linum usitatissimum L.) Fibers for Composite Reinforcement: Exploring the Link Between Plant Growth, Cell Walls Development, and Fiber Properties

Due to the combination of high mechanical performances and plant-based origin, flax fibers are interesting reinforcement for environmentally friendly composite materials. An increasing amount of research articles and reviews focuses on the processing and properties of flax-based products, without taking into account the original key role of flax fibers, namely, reinforcement elements of the flax stem (Linum usitatissimum L.). The ontogeny of the plant, scattering of fiber properties along the plant, or the plant growth conditions are rarely considered. Conversely, exploring the development of flax fibers and parameters influencing the plant mechanical properties (at the whole plant or fiber scale) could be an interesting way to control and/or optimize fiber performances, and to a greater extent, flax fiber-based products. The first part of the present review synthesized the general knowledge about the growth stages of flax plants and the internal organization of the stem biological tissues. Additionally, key findings regarding the development of its fibers, from elongation to thickening, are reviewed to offer a piece of explanation of the uncommon morphological properties of flax fibers. Then, the slenderness of flax is illustrated by comparison of data given in scientific research on herbaceous plants and woody ones. In the second section, a state of the art of the varietal selection of several main industrial crops is given. This section includes the different selection criteria as well as an overview of their impact on plant characteristics. A particular interest is given to the lodging resistance and the understanding of this undesired phenomenon. The third section reviews the influence of the cultural conditions, including seedling rate and its relation with the wind in a plant canopy, as well as the impact of main tropisms (namely, thigmotropism, seismotropism, and gravitropism) on the stem and fiber characteristics. This section illustrates the mechanisms of plant adaptation, and how the environment can modify the plant biomechanical properties. Finally, this review asks botanists, breeders, and farmers' knowledge toward the selection of potential flax varieties dedicated to composite applications, through optimized fiber performances. All along the paper, both fibers morphology and mechanical properties are discussed, in constant link with their use for composite materials reinforcement.

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

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

Reduced expression of selected FASCICLIN-LIKE ARABINOGALACTAN PROTEIN genes associates with the abortion of kernels in field crops of Zea mays (maize) and of Arabidopsis seeds

Abortion of fertilized ovaries at the tip of the ear can generate significant yield losses in maize crops. To investigate the mechanisms involved in this process, 2 maize hybrids were grown in field crops at 2 sowing densities and under 3 irrigation regimes (well-watered control, drought before pollination, and drought during pollination), in all possible combinations. Samples of ear tips were taken 2-6 days after synchronous hand pollination and used for the analysis of gene expression and sugars. Glucose and fructose levels increased in kernels with high abortion risk. Several FASCICLIN-LIKE ARABINOGALACTAN PROTEIN (FLA) genes showed negative correlation with abortion. The expression of ZmFLA7 responded to drought only at the tip of the ear. The abundance of arabinogalactan protein (AGP) glycan epitopes decreased with drought and pharmacological treatments that reduce AGP activity enhanced the abortion of fertilized ovaries. Drought also reduced the expression of AthFLA9 in the siliques of Arabidopsis thaliana. Gain- and loss-of-function mutants of Arabidopsis showed a negative correlation between AthFLA9 and seed abortion. On the basis of gene expression patterns, pharmacological, and genetic evidence, we propose that stress-induced reductions in the expression of selected FLA genes enhance abortion of fertilized ovaries in maize and Arabidopsis.

Phloem fibres as motors of gravitropic behaviour of flax plants: level of transcriptome

Restoration of stem vertical position after plant inclination is a widely spread version of plant orientation in accordance with gravity vector direction. Gravitropic behaviour of flax plants involves the formation of curvature in stem region that has ceased elongation long in advance of stem inclination. The important participants of such behaviour are phloem fibres with constitutively formed tertiary cell wall (G-layer). We performed the large-scale transcriptome profiling of phloem fibres isolated from pulling and opposite sides of gravitropic curvature and compared with control plant fibres. Significant changes in transcript abundance take place for genes encoding proteins of several ion channels, transcription factors and other regulating elements. The largest number of upregulated genes belonged to the cell wall category; many of those were specifically upregulated in fibres of pulling stem side. The obtained data permit to suggest the mechanism of fibre participation in gravitropic reaction that involves the increase of turgor pressure and the rearrangements of cell wall structure in order to improve contractile properties, and to identify the regulatory elements that operate specifically in the fibres of the pulling stem side making gelatinous phloem fibres an important element of gravitropic response in herbaceous plants.

Evolutionary characterization and transcript profiling of β-tubulin genes in flax (Linum usitatissimum L.) during plant development

BACKGROUND

Microtubules, polymerized from alpha and beta-tubulin monomers, play a fundamental role in plant morphogenesis, determining the cell division plane, the direction of cell expansion and the deposition of cell wall material. During polarized pollen tube elongation, microtubules serve as tracks for vesicular transport and deposition of proteins/lipids at the tip membrane. Such functions are controlled by cortical microtubule arrays. Aim of this study was to first characterize the flax β-tubulin family by sequence and phylogenetic analysis and to investigate differential expression of β-tubulin genes possibly related to fibre elongation and to flower development.

RESULTS

We report the cloning and characterization of the complete flax β-tubulin gene family: exon-intron organization, duplicated gene comparison, phylogenetic analysis and expression pattern during stem and hypocotyl elongation and during flower development. Sequence analysis of the fourteen expressed β-tubulin genes revealed that the recent whole genome duplication of the flax genome was followed by massive retention of duplicated tubulin genes. Expression analysis showed that β-tubulin mRNA profiles gradually changed along with phloem fibre development in both the stem and hypocotyl. In flowers, changes in relative tubulin transcript levels took place at anthesis in anthers, but not in carpels.

CONCLUSIONS

Phylogenetic analysis supports the origin of extant plant β-tubulin genes from four ancestral genes pre-dating angiosperm separation. Expression analysis suggests that particular tubulin subpopulations are more suitable to sustain different microtubule functions such as cell elongation, cell wall thickening or pollen tube growth. Tubulin genes possibly related to different microtubule functions were identified as candidate for more detailed studies.

Identification of fasciclin-like arabinogalactan proteins in textile hemp (Cannabis sativa L.): in silico analyses and gene expression patterns in different tissues

BACKGROUND

The fasciclin-like arabinogalactan proteins (FLAs) belong to the arabinogalactan protein (AGP) superfamily and are known to play different physiological roles in plants. This class of proteins was shown to participate in plant growth, development, defense against abiotic stresses and, notably, cell wall biosynthesis. Although some studies are available on the characterization of FLA genes from different species, both woody and herbaceous, no detailed information is available on the FLA family of textile hemp (Cannabis sativa L.), an economically important fibre crop.

RESULTS

By searching the Cannabis genome and EST databases, 23 CsaFLAs have been here identified which are divided into four phylogenetic groups. A real-time qPCR analysis performed on stem tissues (isolated bast fibres and shivs sampled at three heights), hypocotyls (6-9-12-15-17-20 days-old), whole seedlings, roots, leaves and female/male flowers of the monoecious fibre variety Santhica 27, indicates that the identified FLA genes are differentially expressed. Interestingly, some hemp FLAs are expressed during early phases of fibre growth (elongation), while others are more expressed in the middle and base of the stem and thus potentially involved in secondary cell wall formation (fibre thickening). The bioinformatic analysis of the promoter regions shows that the FLAs upregulated in the younger regions of the stem share a conserved motif related to flowering control and regulation of photoperiod perception. The promoters of the FLA genes expressed at higher levels in the older stem regions, instead, share a motif putatively recognized by MYB3, a transcriptional repressor belonging to the MYB family subgroup S4.

CONCLUSIONS

These results point to the existence of a transcriptional network fine-tuning the expression of FLA genes in the older and younger regions of the stem, as well as in the bast fibres/shivs of textile hemp. In summary, our study paves the way for future analyses on the biological functions of FLAs in an industrially relevant fibre crop.

A Cell Wall Proteome and Targeted Cell Wall Analyses Provide Novel Information on Hemicellulose Metabolism in Flax

Experimentally-generated (nanoLC-MS/MS) proteomic analyses of four different flax organs/tissues (inner-stem, outer-stem, leaves and roots) enriched in proteins from 3 different sub-compartments (soluble-, membrane-, and cell wall-proteins) was combined with publically available data on flax seed and whole-stem proteins to generate a flax protein database containing 2996 nonredundant total proteins. Subsequent multiple analyses (MapMan, CAZy, WallProtDB and expert curation) of this database were then used to identify a flax cell wall proteome consisting of 456 nonredundant proteins localized in the cell wall and/or associated with cell wall biosynthesis, remodeling and other cell wall related processes. Examination of the proteins present in different flax organs/tissues provided a detailed overview of cell wall metabolism and highlighted the importance of hemicellulose and pectin remodeling in stem tissues. Phylogenetic analyses of proteins in the cell wall proteome revealed an important paralogy in the class IIIA xyloglucan endo-transglycosylase/hydrolase (XTH) family associated with xyloglucan endo-hydrolase activity.Immunolocalisation, FT-IR microspectroscopy, and enzymatic fingerprinting indicated that flax fiber primary/S1 cell walls contained xyloglucans with typical substituted side chains as well as glucuronoxylans in much lower quantities. These results suggest a likely central role of xyloglucans and endotransglucosylase/hydrolase activity in flax fiber formation and cell wall remodeling processes.

Spatial regulation of monolignol biosynthesis and laccase genes control developmental and stress-related lignin in flax

BACKGROUND

Bast fibres are characterized by very thick secondary cell walls containing high amounts of cellulose and low lignin contents in contrast to the heavily lignified cell walls typically found in the xylem tissues. To improve the quality of the fiber-based products in the future, a thorough understanding of the main cell wall polymer biosynthetic pathways is required. In this study we have carried out a characterization of the genes involved in lignin biosynthesis in flax along with some of their regulation mechanisms.

RESULTS

We have first identified the members of the phenylpropanoid gene families through a combination of in silico approaches. The more specific lignin genes were further characterized by high throughput transcriptomic approaches in different organs and physiological conditions and their cell/tissue expression was localized in the stems, roots and leaves. Laccases play an important role in the polymerization of monolignols. This multigenic family was determined and a miRNA was identified to play a role in the posttranscriptional regulation by cleaving the transcripts of some specific genes shown to be expressed in lignified tissues. In situ hybridization also showed that the miRNA precursor was expressed in the young xylem cells located near the vascular cambium. The results obtained in this work also allowed us to determine that most of the genes involved in lignin biosynthesis are included in a unique co-expression cluster and that MYB transcription factors are potentially good candidates for regulating these genes.

CONCLUSIONS

Target engineering of cell walls to improve plant product quality requires good knowledge of the genes responsible for the production of the main polymers. For bast fiber plants such as flax, it is important to target the correct genes from the beginning since the difficulty to produce transgenic material does not make possible to test a large number of genes. Our work determined which of these genes could be potentially modified and showed that it was possible to target different regulatory pathways to modify lignification.

Transcriptomic profiling of hemp bast fibres at different developmental stages

Bast fibres are long extraxylary cells which mechanically support the phloem and they are divided into xylan- and gelatinous-type, depending on the composition of their secondary cell walls. The former, typical of jute/kenaf bast fibres, are characterized by the presence of xylan and a high degree of lignification, while the latter, found in tension wood, as well as flax, ramie and hemp bast fibres, have a high abundance of crystalline cellulose. During their differentiation, bast fibres undergo specific developmental stages: the cells initially elongate rapidly by intrusive growth, subsequently they cease elongation and start to thicken. The goal of the present study is to provide a transcriptomic close-up of the key events accompanying bast fibre development in textile hemp (Cannabis sativa L.), a fibre crop of great importance. Bast fibres have been sampled from different stem regions. The developmental stages corresponding to active elongation and cell wall thickening have been studied using RNA-Seq. The results show that the fibres sampled at each stem region are characterized by a specific transcriptomic signature and that the major changes in cell wall-related processes take place at the internode containing the snap point. The data generated also identify several interesting candidates for future functional analysis.

Digital gene expression profiling of flax (Linum usitatissimum L.) stem peel identifies genes enriched in fiber-bearing phloem tissue

To better understand the molecular mechanisms and gene expression characteristics associated with development of bast fiber cell within flax stem phloem, the gene expression profiling of flax stem peels and leaves were screened, using Illumina's Digital Gene Expression (DGE) analysis. Four DGE libraries (2 for stem peel and 2 for leaf), ranging from 6.7 to 9.2 million clean reads were obtained, which produced 7.0 million and 6.8 million mapped reads for flax stem peel and leave, respectively. By differential gene expression analysis, a total of 975 genes, of which 708 (73%) genes have protein-coding annotation, were identified as phloem enriched genes putatively involved in the processes of polysaccharide and cell wall metabolism. Differential expression genes (DEGs) was validated using quantitative RT-PCR, the expression pattern of all nine genes determined by qRT-PCR fitted in well with that obtained by sequencing analysis. Cluster and Gene Ontology (GO) analysis revealed that a large number of genes related to metabolic process, catalytic activity and binding category were expressed predominantly in the stem peels. The Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of the phloem enriched genes suggested approximately 111 biological pathways. The large number of genes and pathways produced from DGE sequencing will expand our understanding of the complex molecular and cellular events in flax bast fiber development and provide a foundation for future studies on fiber development in other bast fiber crops.

Transcriptome portrait of cellulose-enriched flax fibres at advanced stage of specialization

Functional specialization of cells is among the most fundamental processes of higher organism ontogenesis. The major obstacle to studying this phenomenon in plants is the difficulty of isolating certain types of cells at defined stages of in planta development for in-depth analysis. A rare opportunity is given by the developed model system of flax (Linum usitatissimum L.) phloem fibres that can be purified from the surrounding tissues at the stage of the tertiary cell wall deposition. The performed comparison of the whole transcriptome profile in isolated fibres and other portions of the flax stem, together with fibre metabolism characterization, helped to elucidate the general picture of the advanced stage of plant cell specialization and to reveal novel participants potentially involved in fibre metabolism regulation and cell wall formation. Down-regulation of all genes encoding proteins involved in xylan and lignin synthesis and up-regulation of genes for the specific set of transcription factors transcribed during tertiary cell wall formation were revealed. The increased abundance of transcripts for several glycosyltransferases indicated the enzymes that may be involved in synthesis of fibre-specific version of rhamnogalacturonan I.

Development of gravitropic response: unusual behavior of flax phloem G-fibers

The major mechanism of gravitropism that is discussed for herbal plants is based on the nonuniform elongation of cells located on the opposite stem sides, occurring in the growing zone of an organ. However, gravitropic response of flax (Linum usitatissimum L.) is well-pronounced in the lower half of developing stem, which has ceased elongation long in advance of plant inclination. We have analyzed the stem curvature region by various approaches of microscopy and found the undescribed earlier significant modifications in primary phloem fibers that have constitutively developed G-layer. In fibers on the pulling stem side, cell portions were widened with formation of "bottlenecks" between them, leading to the "sausage-like" shape of a cell. Lumen diameter in fiber widening increased, while cell wall thickness decreased. Callose was deposited in proximity to bottlenecks and sometimes totally occluded their lumen. Structure of fiber cell wall changed considerably, with formation of breaks between G- and S-layers. Thick fibrillar structures that were revealed in fiber cell wall by light microscopy got oblique orientation instead of parallel to the fiber axis one in control plants. The described changes occurred at various combinations of gravitational and mechanical stimuli. Thus, phloem fibers with constitutively formed gelatinous cell wall, located in nonelongating parts of herbal plant, are involved in gravitropism and may become an important element in general understanding of the gravity effects on plants. We suggest flax phloem fibers as the model system to study the mechanism of plant position correction, including signal perception and transduction.

Studying Secondary Growth and Bast Fiber Development: The Hemp Hypocotyl Peeks behind the Wall

Cannabis sativa L. is an annual herbaceous crop grown for the production of long extraxylary fibers, the bast fibers, rich in cellulose and used both in the textile and biocomposite sectors. Despite being herbaceous, hemp undergoes secondary growth and this is well exemplified by the hypocotyl. The hypocotyl was already shown to be a suitable model to study secondary growth in other herbaceous species, namely Arabidopsis thaliana and it shows an important practical advantage, i.e., elongation and radial thickening are temporally separated. This study focuses on the mechanisms marking the transition from primary to secondary growth in the hemp hypocotyl by analysing the suite of events accompanying vascular tissue and bast fiber development. Transcriptomics, imaging and quantification of phytohormones were carried out on four representative developmental stages (i.e., 6-9-15-20 days after sowing) to provide a comprehensive overview of the events associated with primary and secondary growth in hemp. This multidisciplinary approach provides cell wall-related snapshots of the growing hemp hypocotyl and identifies marker genes associated with the young (expansins, β-galactosidases, and transcription factors involved in light-related processes) and the older hypocotyl (secondary cell wall biosynthetic genes and transcription factors).

Biotechnology Towards Energy Crops

New crops are gradually establishing along with cultivation systems to reduce reliance on depleting fossil fuel reserves and sustain better adaptation to climate change. These biological assets could be efficiently exploited as bioenergy feedstocks. Bioenergy crops are versatile renewable sources with the potential to alternatively contribute on a daily basis towards the coverage of modern society's energy demands. Biotechnology may facilitate the breeding of elite energy crop genotypes, better suited for bio-processing and subsequent use that will improve efficiency, further reduce costs, and enhance the environmental benefits of biofuels. Innovative molecular techniques may improve a broad range of important features including biomass yield, product quality and resistance to biotic factors like pests or microbial diseases or environmental cues such as drought, salinity, freezing injury or heat shock. The current review intends to assess the capacity of biotechnological applications to develop a beneficial bioenergy pipeline extending from feedstock development to sustainable biofuel production and provide examples of the current state of the art on future energy crops.

Cloning of expansin genes in ramie (Boehmeria nivea L.) based on universal fast walking

Gene cloning is the first step to study the expression profiles and functions of a particular gene; considerable cloning methods have been developed. Expansin, thought to involve in the cell-wall modification events, was not cloned in ramie (Boehmeria nivea L.), which is one of the most important bast fiber crops with little conducted molecular research, especially on its fiber development. Studying the expansin gene family will uncover its possible relationship with ramie fiber development and other growth events. As a result, five expansin genes were cloned with full-length and their sequence information was investigated. Additionally, the phylogenetic analysis was conducted, which suggested that the cloned genes belong to the α-subfamily, and these genes expressed differently during ramie fiber developmental process. In this study, we aimed to apply a strategy for cloning novel full-length genes from genomic DNA of ramie, based on using degenerate primers, touchdown polymerase chain reaction and universal fast walking protocols. By cloning five full-length expansin genes, we believe the polymerase chain reaction-based gene cloning strategy could be applied to general gene studies in ramie and other crops.

Pathways associated with lignin biosynthesis in lignomaniac jute fibres

We generated the bast transcriptomes of a deficient lignified phloem fibre mutant and its wild-type jute (Corchorus capsularis) using Illumina paired-end sequencing. A total of 34,163 wild-type and 29,463 mutant unigenes, with average lengths of 1442 and 1136 bp, respectively, were assembled de novo, ~77-79 % of which were functionally annotated. These annotated unigenes were assigned to COG (~37-40 %) and GO (~22-28 %) classifications and mapped to 189 KEGG pathways (~19-21 %). We discovered 38 and 43 isoforms of 16 and 10 genes of the upstream shikimate-aromatic amino acid and downstream monolignol biosynthetic pathways, respectively, rendered their sequence similarities, confirmed the identities of 22 of these candidate gene families by phylogenetic analyses and reconstructed the pathway leading to lignin biosynthesis in jute fibres. We also identified major genes and bast-related transcription factors involved in secondary cell wall (SCW) formation. The quantitative RT-PCRs revealed that phenylalanine ammonia-lyase 1 (CcPAL1) was co-down-regulated with several genes of the upstream shikimate pathway in mutant bast tissues at an early growth stage, although its expression relapsed to the normal level at the later growth stage. However, cinnamyl alcohol dehydrogenase 7 (CcCAD7) was strongly down-regulated in mutant bast tissues irrespective of growth stages. CcCAD7 disruption at an early growth stage was accompanied by co-up-regulation of SCW-specific genes cellulose synthase A7 (CcCesA7) and fasciclin-like arabinogalactan 6 (CcFLA6), which was predicted to be involved in coordinating the S-layers' deposition in the xylan-type jute fibres. Our results identified CAD as a promising target for developing low-lignin jute fibres using genomics-assisted molecular approaches.

Ectopic lignification in the flax lignified bast fiber1 mutant stem is associated with tissue-specific modifications in gene expression and cell wall composition

Histochemical screening of a flax ethyl methanesulfonate population led to the identification of 93 independent M2 mutant families showing ectopic lignification in the secondary cell wall of stem bast fibers. We named this core collection the Linum usitatissimum (flax) lbf mutants for lignified bast fibers and believe that this population represents a novel biological resource for investigating how bast fiber plants regulate lignin biosynthesis. As a proof of concept, we characterized the lbf1 mutant and showed that the lignin content increased by 350% in outer stem tissues containing bast fibers but was unchanged in inner stem tissues containing xylem. Chemical and NMR analyses indicated that bast fiber ectopic lignin was highly condensed and rich in G-units. Liquid chromatography-mass spectrometry profiling showed large modifications in the oligolignol pool of lbf1 inner- and outer-stem tissues that could be related to ectopic lignification. Immunological and chemical analyses revealed that lbf1 mutants also showed changes to other cell wall polymers. Whole-genome transcriptomics suggested that ectopic lignification of flax bast fibers could be caused by increased transcript accumulation of (1) the cinnamoyl-CoA reductase, cinnamyl alcohol dehydrogenase, and caffeic acid O-methyltransferase monolignol biosynthesis genes, (2) several lignin-associated peroxidase genes, and (3) genes coding for respiratory burst oxidase homolog NADPH-oxidases necessary to increase H2O2 supply.

Transcriptome profiling using pyrosequencing shows genes associated with bast fiber development in ramie (Boehmeria nivea L.)

Background

Ramie (Boehmeria nivea L.), popularly known as “China grass”, is one of the oldest crops in China and the second most important fiber crop in terms of area sown. Ramie fiber, extracted from the plant bast, is important in the textile industry. However, the molecular mechanism of ramie fiber development remains unknown.

Results

A whole sequencing run was performed on the 454 GS FLX + platform using four separately pooled parts of ramie bast. This generated 1,030,057 reads with an average length of 457 bp. Among the 58,369 unigenes (13,386 contigs and 44,983 isotigs) that were generated through de novo assembly, 780 were differentially expressed. As a result, 13 genes that belong to the cellulose synthase gene family (four), the expansin gene family (three) and the xyloglucan endotransglucosylase/hydrolase (XTH) gene family (six) were up-regulated in the top part of the bast, which was in contrast to the other three parts. The identification of these 13 concurrently up-regulated unigenes indicated that the early stage (represented by the top part of the bast) might be important for the molecular regulation of ramie fiber development. Further analysis indicated that four of the 13 unigenes from the expansin (two) and XTH (two) families shared a coincident expression pattern during the whole growth season, which implied they were more relevant to ramie fiber development (fiber quality, etc.) during the different seasons than the other genes.

Conclusions

To the best of our knowledge, this study is the first to characterize ramie fiber development at different developmental stages. The identified transcripts will improve our understanding of the molecular mechanisms involved in ramie fiber development. Moreover, the identified differentially expressed genes will accelerate molecular research on ramie fiber growth and the breeding of ramie with better fiber yields and quality.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-919) contains supplementary material, which is available to authorized users.

Chitinase-like (CTL) and cellulose synthase (CESA) gene expression in gelatinous-type cellulosic walls of flax (Linum usitatissimum L.) bast fibers

Plant chitinases (EC 3.2.1.14) and chitinase-like (CTL) proteins have diverse functions including cell wall biosynthesis and disease resistance. We analyzed the expression of 34 chitinase and chitinase-like genes of flax (collectively referred to as LusCTLs), belonging to glycoside hydrolase family 19 (GH19). Analysis of the transcript expression patterns of LusCTLs in the stem and other tissues identified three transcripts (LusCTL19, LusCTL20, LusCTL21) that were highly enriched in developing bast fibers, which form cellulose-rich gelatinous-type cell walls. The same three genes had low relative expression in tissues with primary cell walls and in xylem, which forms a xylan type of secondary cell wall. Phylogenetic analysis of the LusCTLs identified a flax-specific sub-group that was not represented in any of other genomes queried. To provide further context for the gene expression analysis, we also conducted phylogenetic and expression analysis of the cellulose synthase (CESA) family genes of flax, and found that expression of secondary wall-type LusCESAs (LusCESA4, LusCESA7 and LusCESA8) was correlated with the expression of two LusCTLs (LusCTL1, LusCTL2) that were the most highly enriched in xylem. The expression of LusCTL19, LusCTL20, and LusCTL21 was not correlated with that of any CESA subgroup. These results defined a distinct type of CTLs that may have novel functions specific to the development of the gelatinous (G-type) cellulosic walls.

PT-Flax (phenotyping and TILLinG of flax): development of a flax (Linum usitatissimum L.) mutant population and TILLinG platform for forward and reverse genetics

BACKGROUND

Flax (Linum usitatissimum L.) is an economically important fiber and oil crop that has been grown for thousands of years. The genome has been recently sequenced and transcriptomics are providing information on candidate genes potentially related to agronomically-important traits. In order to accelerate functional characterization of these genes we have generated a flax EMS mutant population that can be used as a TILLinG (Targeting Induced Local Lesions in Genomes) platform for forward and reverse genetics.

RESULTS

A population of 4,894 M2 mutant seed families was generated using 3 different EMS concentrations (0.3%, 0.6% and 0.75%) and used to produce M2 plants for subsequent phenotyping and DNA extraction. 10,839 viable M2 plants (4,033 families) were obtained and 1,552 families (38.5%) showed a visual developmental phenotype (stem size and diameter, plant architecture, flower-related). The majority of these families showed more than one phenotype. Mutant phenotype data are organised in a database and can be accessed and searched at UTILLdb (http://urgv.evry.inra.fr/UTILLdb). Preliminary screens were also performed for atypical fiber and seed phenotypes. Genomic DNA was extracted from 3,515 M2 families and eight-fold pooled for subsequent mutant detection by ENDO1 nuclease mis-match cleavage. In order to validate the collection for reverse genetics, DNA pools were screened for two genes coding enzymes of the lignin biosynthesis pathway: Coumarate-3-Hydroxylase (C3H) and Cinnamyl Alcohol Dehydrogenase (CAD). We identified 79 and 76 mutations in the C3H and CAD genes, respectively. The average mutation rate was calculated as 1/41 Kb giving rise to approximately 9,000 mutations per genome. Thirty-five out of the 52 flax cad mutant families containing missense or codon stop mutations showed the typical orange-brown xylem phenotype observed in CAD down-regulated/mutant plants in other species.

CONCLUSIONS

We have developed a flax mutant population that can be used as an efficient forward and reverse genetics tool. The collection has an extremely high mutation rate that enables the detection of large numbers of independant mutant families by screening a comparatively low number of M2 families. The population will prove to be a valuable resource for both fundamental research and the identification of agronomically-important genes for crop improvement in flax.

Integrated -omics: a powerful approach to understanding the heterogeneous lignification of fibre crops

Lignin and cellulose represent the two main components of plant secondary walls and the most abundant polymers on Earth. Quantitatively one of the principal products of the phenylpropanoid pathway, lignin confers high mechanical strength and hydrophobicity to plant walls, thus enabling erect growth and high-pressure water transport in the vessels. Lignin is characterized by a high natural heterogeneity in its composition and abundance in plant secondary cell walls, even in the different tissues of the same plant. A typical example is the stem of fibre crops, which shows a lignified core enveloped by a cellulosic, lignin-poor cortex. Despite the great value of fibre crops for humanity, however, still little is known on the mechanisms controlling their cell wall biogenesis, and particularly, what regulates their spatially-defined lignification pattern. Given the chemical complexity and the heterogeneous composition of fibre crops' secondary walls, only the use of multidisciplinary approaches can convey an integrated picture and provide exhaustive information covering different levels of biological complexity. The present review highlights the importance of combining high throughput -omics approaches to get a complete understanding of the factors regulating the lignification heterogeneity typical of fibre crops.

LuFLA1PRO and LuBGAL1PRO promote gene expression in the phloem fibres of flax (Linum usitatissimum)

Cell type-specific promoters were identified that drive gene expression in an industrially important product. To identify flax (Linum usitatissimum) gene promoters, we analyzed the genomic regions upstream of a fasciclin-like arabinogalactan protein (LuFLA1) and a beta-galactosidase (LuBGAL1). Both of these genes encode transcripts that have been found to be highly enriched in tissues bearing phloem fibres. Using a beta-glucuronidase (GUS) reporter construct, we found that a 908-bp genomic sequence upstream of LuFLA1 (LuFLA1PRO) directed GUS expression with high specificity to phloem fibres undergoing secondary cell wall development. The DNA sequence upstream of LuBGAL1 (LuBGAL1PRO) likewise produced GUS staining in phloem fibres with developing secondary walls, as well as in tissues of developing flowers and seed bolls. These data provide further evidence of a specific role for LuFLA1 in phloem fibre development, and demonstrate the utility of LuFLA1PRO and LuBGAL1PRO as tools for biotechnology and further investigations of phloem fibre development.

Nutrition metabolism plays an important role in the alternate bearing of the olive tree (Olea europaea L.)

The olive tree (Olea europaea L.) is widely known for its strong tendency for alternate bearing, which severely affects the fruit yield from year to year. Microarray based gene expression analysis using RNA from olive samples (on-off years leaves and ripe-unripe fruits) are particularly useful to understand the molecular mechanisms influencing the periodicity in the olive tree. Thus, we carried out genome wide transcriptome analyses involving different organs and temporal stages of the olive tree using the NimbleGen Array containing 136,628 oligonucleotide probe sets. Cluster analyses of the genes showed that cDNAs originated from different organs could be sorted into separate groups. The nutritional control had a particularly remarkable impact on the alternate bearing of olive, as shown by the differential expression of transcripts under different temporal phases and organs. Additionally, hormonal control and flowering processes also played important roles in this phenomenon. Our analyses provide further insights into the transcript changes between "on year" and "off year" leaves along with the changes from unrpipe to ripe fruits, which shed light on the molecular mechanisms underlying the olive tree alternate bearing. These findings have important implications for the breeding and agriculture of the olive tree and other crops showing periodicity. To our knowledge, this is the first study reporting the development and use of an olive array to document the gene expression profiling associated with the alternate bearing in olive tree.

Intra-platform comparison of 25-mer and 60-mer oligonucleotide Nimblegen DNA microarrays

Background

We performed a Nimblegen intra-platform microarray comparison by assessing two categories of flax target probes (short 25-mers oligonucleotides and long 60-mers oligonucleotides) in identical conditions of target production, design, labelling, hybridization, image analyses, and data filtering. We compared technical parameters of array hybridizations, precision and accuracy as well as specific gene expression profiles.

Results

Comparison of the hybridization quality, precision and accuracy of expression measurements, as well as an interpretation of differential gene expression in flax tissues were performed. Both array types yielded reproducible, accurate and comparable data that are coherent for expression measurements and identification of differentially expressed genes. 60-mers arrays gave higher hybridization efficiencies and therefore were more sensitive allowing the detection of a higher number of unigenes involved in the same biological process and/or belonging to the same multigene family.

Conclusion

The two flax arrays provide a good resolution of expressed functions; however the 60-mers arrays are more sensitive and provide a more in-depth coverage of candidate genes potentially involved in different biological processes.

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

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

Analysis of 2,297 expressed sequence tags (ESTs) from a cDNA library of flax (Linum ustitatissimum L.) bark tissue

Bast fibre crops are the second most important natural fibre crops following cotton. Of these, flax (Linum ustitatissimum L.) is the most widely planted in the world, with its fibre used for high quality linen textile. A cDNA library of flax bark tissues was constructed with the purpose of identifying genes involved in the Bast fibre development. A total of 2,297 unigene sequences were obtained from 3,200 randomly selected clones of the cDNA library. These sequences were grouped into 155 clusters and 2,142 singletons, which have been submitted to the GenBank databases. By putative functional annotation, 23.3% of these sequences were similar to known proteins in GenBank, 44.0% of these sequences were similar to unknown proteins, and 32.7% of these sequences showed no significant similarity to any other protein sequences in existing databases. Classified by the Gene Ontology, 24.8, 23.1 and 14.3% were assigned to molecular function, biological process, and cellular component GO terms, respectively. By further bioinformatics approaches, about 110 ESTs matched cell wall related genes in the MAIZEWALL database, representing 16 functional categories of all 19 categories, of which, the most abundant category was protein synthesis. Based on the PlantTFDB database, 39 of the 64 transcription factor families in the Arabidopsis thaliana genome were identified as being involved in flax cell wall formation. The sequences and bioinformatics analysis data generated in this paper will be useful for gene expression, cloning and genetic engineering studies to characterize bast fibre development and improve the properties of the bast fibres.

The Arabidopsis CORI3 promoter contains two cis-acting regulatory regions required for transcriptional activity in companion cells

Companion cells are metabolically active and functionally specialized cells that behave as terminals for the transport of materials between phloem and the surrounding tissue. Although previous research has clarified the distinct function of companion cells, it is still largely unknown how plants establish and maintain the special identity of these cells. To shed further light on this issue, we carried out expressed sequence tag (EST) analysis. To minimize the difficulty of dissociating and gathering intact companion cells, vascular strings with an abundant content of companion cells were excised from the petioles of Brassica napus. By random sequencing with a string-specific cDNA library derived by suppression subtractive hybridization between the string itself and the petiole from which it had been removed, we identified 377 ESTs and assembled them into 247 EST groups. The most frequent EST was ExBdl-102 (15 of 377 ESTs), which showed the highest sequence similarity to the Arabidopsis CORI3 (CORONATINE INDUCED 3) gene. The CORI3 promoter:GUS showed predominant expression in the vascular tissue of Arabidopsis. Through transient expression assay using Brassica vasculature and gene-gun-mediated transient assay, we found two integrated cis-regulatory regions of the CORI3 promoter. This work has provided not only string-specific EST information and shown that two novel cis-regulatory regions sustain transcriptional activity in companion cells, but also a series of procedures for efficiently examining the transcriptional framework of companion cells by exploiting the histochemical advantage of B. napus as an experimental material.

Development of cellulosic secondary walls in flax fibers requires beta-galactosidase

Bast (phloem) fibers, tension wood fibers, and other cells with gelatinous-type secondary walls are rich in crystalline cellulose. In developing bast fibers of flax (Linum usitatissimum), a galactan-enriched matrix (Gn-layer) is gradually modified into a mature cellulosic gelatinous-layer (G-layer), which ultimately comprises most of the secondary cell wall. Previous studies have correlated this maturation process with expression of a putative β-galactosidase. Here, we demonstrate that β-galactosidase activity is in fact necessary for the dynamic remodeling of polysaccharides that occurs during normal secondary wall development in flax fibers. We found that developing stems of transgenic (LuBGAL-RNAi) flax with reduced β-galactosidase activity had lower concentrations of free Gal and had significant reductions in the thickness of mature cellulosic G-layers compared with controls. Conversely, Gn-layers, labeled intensively by the galactan-specific LM5 antibody, were greatly expanded in LuBGAL-RNAi transgenic plants. Gross morphology and stem anatomy, including the thickness of bast fiber walls, were otherwise unaffected by silencing of β-galactosidase transcripts. These results demonstrate a specific requirement for β-galactosidase in hydrolysis of galactans during formation of cellulosic G-layers. Transgenic lines with reduced β-galactosidase activity also had biochemical and spectroscopic properties consistent with a reduction in cellulose crystallinity. We further demonstrated that the tensile strength of normal flax stems is dependent on β-galactosidase-mediated development of the phloem fiber G-layer. Thus, the mechanical strength that typifies flax stems is dependent on a thick, cellulosic G-layer, which itself depends on β-galactosidase activity within the precursor Gn-layer. These observations demonstrate a novel role for matrix polysaccharides in cellulose deposition; the relevance of these observations to the development of cell walls in other species is also discussed.

Transcript profiling and gene characterization of three fatty acid desaturase genes in high, moderate, and low linolenic acid genotypes of flax (Linum usitatissimum L.) and their role in linolenic acid accumulation

Three genes encoding fatty acid desaturase 3 (fad3a, fad3b, and a novel fad3c) were cloned from four flax genotypes varying in linolenic acid content. Real-time PCR was used to quantify expression levels of the three fad3 genes during seed development. High amounts of both fad3a and fad3b transcripts were observed and reached their peak levels at 20 days after anthesis, except for fad3a from SP2047 where only low level expression was observed throughout seed development. Transcript accumulation of the novel fad3c gene was at similar background levels. The fatty acid composition was analysed for all genotypes and stages of development and compared with the fad3 gene expression patterns. α-Linolenic acid gradually accumulated during seed development, while linoleic acid was transient and decreased in M5791, UGG5-5, and AC McDuff. In contrast, the linolenic acid present in the early stages of development nearly completely disappeared in SP2047, while linoleic acid steadily accumulated. fad3a of the low linolenic acid line SP2047 encoded a truncated protein caused by a premature stop codon resulting from a single point mutation, and the low level of transcript accumulation in this genotype is likely due to nonsense-mediated mRNA decay caused by the premature termination of translation as a result of this early stop codon. Although substantial amounts of transcript accumulation occurred with fad3b of SP2047 genotype, cloning of the gene revealed a mutation in the first histidine box causing an amino acid change. Heterologous expression in yeast of the SP2047 and UGG5-5 fad3b genes showed that the mutation in the histidine box in SP2047 caused the enzyme inactivity. Taken together, these results showed that fad3a and fad3b are responsible for linolenic acid accumulation in flax seeds but did not support a major role for the novel fad3c. These observations were further supported by phenotypic and genotypic assessment of a doubled haploid population. Expression patterns of fad3a and fad3b were highly correlated with linolenic acid accumulation during seed development, with the exception of fad3b in SP2047 whose lack of activity was caused by the histidine box mutation despite its transcript accumulation being similar to that of the fad3b of the other genotypes.

Gene expression analysis of flax seed development

BACKGROUND

Flax, Linum usitatissimum L., is an important crop whose seed oil and stem fiber have multiple industrial applications. Flax seeds are also well-known for their nutritional attributes, viz., omega-3 fatty acids in the oil and lignans and mucilage from the seed coat. In spite of the importance of this crop, there are few molecular resources that can be utilized toward improving seed traits. Here, we describe flax embryo and seed development and generation of comprehensive genomic resources for the flax seed.

RESULTS

We describe a large-scale generation and analysis of expressed sequences in various tissues. Collectively, the 13 libraries we have used provide a broad representation of genes active in developing embryos (globular, heart, torpedo, cotyledon and mature stages) seed coats (globular and torpedo stages) and endosperm (pooled globular to torpedo stages) and genes expressed in flowers, etiolated seedlings, leaves, and stem tissue. A total of 261,272 expressed sequence tags (EST) (GenBank accessions LIBEST_026995 to LIBEST_027011) were generated. These EST libraries included transcription factor genes that are typically expressed at low levels, indicating that the depth is adequate for in silico expression analysis. Assembly of the ESTs resulted in 30,640 unigenes and 82% of these could be identified on the basis of homology to known and hypothetical genes from other plants. When compared with fully sequenced plant genomes, the flax unigenes resembled poplar and castor bean more than grape, sorghum, rice or Arabidopsis. Nearly one-fifth of these (5,152) had no homologs in sequences reported for any organism, suggesting that this category represents genes that are likely unique to flax. Digital analyses revealed gene expression dynamics for the biosynthesis of a number of important seed constituents during seed development.

CONCLUSIONS

We have developed a foundational database of expressed sequences and collection of plasmid clones that comprise even low-expressed genes such as those encoding transcription factors. This has allowed us to delineate the spatio-temporal aspects of gene expression underlying the biosynthesis of a number of important seed constituents in flax. Flax belongs to a taxonomic group of diverse plants and the large sequence database will allow for evolutionary studies as well.

Members of the LATERAL ORGAN BOUNDARIES DOMAIN transcription factor family are involved in the regulation of secondary growth in Populus

Regulation of secondary (woody) growth is of substantial economic and environmental interest but is poorly understood. We identified and subsequently characterized an activation-tagged poplar (Populus tremula × Populus alba) mutant with enhanced woody growth and changes in bark texture caused primarily by increased secondary phloem production. Molecular characterization of the mutation through positioning of the tag and retransformation experiments shows that the phenotype is conditioned by activation of an uncharacterized gene that encodes a novel member of the LATERAL ORGAN BOUNDARIES DOMAIN (LBD) family of transcription factors. Homology analysis showed highest similarity to an uncharacterized LBD1 gene from Arabidopsis thaliana, and we consequently named it Populus tremula × Populus alba (Pta) LBD1. Dominant-negative suppression of Pta LBD1 via translational fusion with the repressor SRDX domain caused decreased diameter growth and suppressed and highly irregular phloem development. In wild-type plants, LBD1 was most highly expressed in the phloem and cambial zone. Two key Class I KNOTTED1-like homeobox genes that promote meristem identity in the cambium were downregulated, while an Altered Phloem Development gene that is known to promote phloem differentiation was upregulated in the mutant. A set of four LBD genes, including the LBD1 gene, was predominantly expressed in wood-forming tissues, suggesting a broader regulatory role of these transcription factors during secondary woody growth in poplar.

Development and validation of a flax (Linum usitatissimum L.) gene expression oligo microarray

BACKGROUND

Flax (Linum usitatissimum L.) has been cultivated for around 9,000 years and is therefore one of the oldest cultivated species. Today, flax is still grown for its oil (oil-flax or linseed cultivars) and its cellulose-rich fibres (fibre-flax cultivars) used for high-value linen garments and composite materials. Despite the wide industrial use of flax-derived products, and our actual understanding of the regulation of both wood fibre production and oil biosynthesis more information must be acquired in both domains. Recent advances in genomics are now providing opportunities to improve our fundamental knowledge of these complex processes. In this paper we report the development and validation of a high-density oligo microarray platform dedicated to gene expression analyses in flax.

RESULTS

Nine different RNA samples obtained from flax inner- and outer-stems, seeds, leaves and roots were used to generate a collection of 1,066,481 ESTs by massive parallel pyrosequencing. Sequences were assembled into 59,626 unigenes and 48,021 sequences were selected for oligo design and high-density microarray (Nimblegen 385K) fabrication with eight, non-overlapping 25-mers oligos per unigene. 18 independent experiments were used to evaluate the hybridization quality, precision, specificity and accuracy and all results confirmed the high technical quality of our microarray platform. Cross-validation of microarray data was carried out using quantitative qRT-PCR. Nine target genes were selected on the basis of microarray results and reflected the whole range of fold change (both up-regulated and down-regulated genes in different samples). A statistically significant positive correlation was obtained comparing expression levels for each target gene across all biological replicates both in qRT-PCR and microarray results. Further experiments illustrated the capacity of our arrays to detect differential gene expression in a variety of flax tissues as well as between two contrasted flax varieties.

CONCLUSION

All results suggest that our high-density flax oligo-microarray platform can be used as a very sensitive tool for analyzing gene expression in a large variety of tissues as well as in different cultivars. Moreover, this highly reliable platform can also be used for the quantification of mRNA transcriptional profiling in different flax tissues.

Selection of reference genes for quantitative gene expression normalization in flax (Linum usitatissimum L.)

BACKGROUND

Quantitative real-time PCR (qRT-PCR) is currently the most accurate method for detecting differential gene expression. Such an approach depends on the identification of uniformly expressed 'housekeeping genes' (HKGs). Extensive transcriptomic data mining and experimental validation in different model plants have shown that the reliability of these endogenous controls can be influenced by the plant species, growth conditions and organs/tissues examined. It is therefore important to identify the best reference genes to use in each biological system before using qRT-PCR to investigate differential gene expression. In this paper we evaluate different candidate HKGs for developmental transcriptomic studies in the economically-important flax fiber- and oil-crop (Linum usitatissimum L).

RESULTS

Specific primers were designed in order to quantify the expression levels of 20 different potential housekeeping genes in flax roots, internal- and external-stem tissues, leaves and flowers at different developmental stages. After calculations of PCR efficiencies, 13 HKGs were retained and their expression stabilities evaluated by the computer algorithms geNorm and NormFinder. According to geNorm, 2 Transcriptional Elongation Factors (TEFs) and 1 Ubiquitin gene are necessary for normalizing gene expression when all studied samples are considered. However, only 2 TEFs are required for normalizing expression in stem tissues. In contrast, NormFinder identified glyceraldehyde-3-phosphate dehydrogenase (GADPH) as the most stably expressed gene when all samples were grouped together, as well as when samples were classed into different sub-groups.qRT-PCR was then used to investigate the relative expression levels of two splice variants of the flax LuMYB1 gene (homologue of AtMYB59). LuMYB1-1 and LuMYB1-2 were highly expressed in the internal stem tissues as compared to outer stem tissues and other samples. This result was confirmed with both geNorm-designated- and NormFinder-designated-reference genes.

CONCLUSIONS

The use of 2 different statistical algorithms results in the identification of different combinations of flax HKGs for expression data normalization. Despite such differences, the use of geNorm-designated- and NormFinder-designated-reference genes enabled us to accurately compare the expression levels of a flax MYB gene in different organs and tissues. Our identification and validation of suitable flax HKGs will facilitate future developmental transcriptomic studies in this economically-important plant.

Suppression of GhAGP4 gene expression repressed the initiation and elongation of cotton fiber

Cotton fibers, important natural raw materials for the textile industry, are trichomes elongated from epidermal cells of cotton ovules. To date, a number of genes have been shown to be critical for fiber development. In this study, the roles of genes encoding fasciclin-like arabinoglactan proteins (FLAs) in cotton fiber were examined by transforming RNA interfering (RNAi) construct. The RNAi according to the sequence of GhAGP4 caused a significant reduction of its mRNA level, and the expression of other three FLAs (GhAGP2, GhAGP3, GhFLA1) were also partially suppressed. The fiber initiation and fiber elongation were inhibited in the transgenic plants. As for the mature fibers of transgenic cotton, the fiber length became significantly shorter and the fiber quality became worse. In addition, the RNAi of GhAGP4 also affected the cytoskeleton network and the cellulose deposition of fiber cells. Through ovule culture, it was found that the expression of cotton FLA genes were upregulated by GA(3), especially for GhAGP2 and GhAGP4. These results indicate that the FLAs are essential for the initiation and elongation of cotton fiber development.

Identification of a cis-regulatory element that acts in companion cell-specific expression of AtMT2B promoter through the use of Brassica vasculature and gene-gun-mediated transient assay

The molecular basis underlying the development, maintenance and function of companion cells in plants is largely unknown. The identification of several genes expressed specifically in companion cells implies the contribution of specific transcriptional elements to the identity of companion cells. However, less is known about the companion cell-specific transcriptional regulation of promoters. We established a novel assay method using gene-gun delivery of partially deleted promoters to string-containing vascular bundles excised from the petiole of Brassica napus for the rapid identification of cis-elements. To test this system, we analyzed the Arabidopsis METALLOTHIONEIN 2B (MT2B) gene, which is expressed in companion cells. The assay revealed a 49-bp region possessing two predicted cis-regulatory elements: a G-box and an evening element-related sequence (EEr), and EEr showed higher activity. We confirmed the reliability of the result with stable transformants harboring a deleted MT2B promoter:GUS transgene. The lack of EEr completely eliminated the MT2B-like expression, but the lack of G-box did not eliminate it. We conclude that EEr is a major cis-regulatory element of the MT2B promoter. Our method will help to explain the transcriptional background of companion cells through the rapid identification of cis-regulatory regions.

Transcript profiling during fiber development identifies pathways in secondary metabolism and cell wall structure that may contribute to cotton fiber quality

A global gene expression profiling study at different stages of fiber development was undertaken on two cotton species cultivated for fiber, Gossypium hirsutum (L.) and G. barbadense (L.). A large proportion of the genome was expressed during both fiber elongation and subsequent secondary cell wall thickening. There was a major shift in abundance of transcripts for gene regulation, cell organization and metabolism between fiber elongation and fiber thickening that was fundamentally similar in both species. Each stage had its own distinctive features represented by specific metabolic and regulatory genes, a number of which have been noted previously. Many of the genes expressed in the fibers were of a similar type and developmental expression to those seen in other fiber-producing plants, indicating a conservation of mechanisms of cell elongation and wall thickening across diverse plant genera. Secondary metabolism and pectin synthesis and modification genes were amongst the most statistically significant differentially expressed categories between the two species during fiber elongation. The gene profiles of the fiber thickening stage, however, were almost identical between the two species, suggesting that their different final fiber quality properties may be established at earlier stages of fiber development. Expression levels of representative phenylpropanoid and pectin modification genes showed high correlations with specific fiber properties in an inter-specific cotton recombinant inbred line (RIL) population, supporting a role in determining fiber quality.