Development and validation of a flax (Linum usitatissimum L.) gene expression oligo microarray
- Stéphane Fenart1,
- Yves-Placide Assoumou Ndong2,
- Jorge Duarte3,
- Nathalie Rivière3,
- Jeroen Wilmer4,
- Olivier van Wuytswinkel2,
- Anca Lucau1,
- Emmanuelle Cariou5,
- Godfrey Neutelings1,
- Laurent Gutierrez6,
- Brigitte Chabbert7,
- Xavier Guillot8,
- Reynald Tavernier9,
- Simon Hawkins1 and
- Brigitte Thomasset10Email author
© Fenart et al; licensee BioMed Central Ltd. 2010
Received: 29 July 2010
Accepted: 21 October 2010
Published: 21 October 2010
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.
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.
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.
Flax (Linum usitatissimum L.) has long been cultivated by man for its cellulose-rich bast fibres and seeds [1, 2]. The long bast fibres are traditionally used in the textile industry for the production of linen or mixed fibre textiles, and, together with shorter xylem fibres are also used in the automobile and construction industries [3, 4]. Flax seeds are widely integrated into animal feeds  and are also important in human health since they are a major source of omega-3 fatty acids  and biologically active lignan . Linseed oil from seeds is equally used in the fabrication of paint, lacquer and varnish, soap, putty and polymers [8, 9].
In order to breed improved flax varieties we need to increase our fundamental knowledge of flax biology (fibre and seed formation, disease resistance, growth etc.) since many important questions remain unanswered. For example, although we know that fibre quality is related to the extremely low lignin levels found in this cell type as compared to wood fibres, we do not know how lignin biosynthesis is regulated in flax fibres . Similarly, although the different biosynthetic pathways involved in fatty acid production in linseed are known, we still know relatively little about the biological mechanisms controlling linolenic acid levels in traditional flax varieties (45-65% omega-3 fatty acid) compared to solin or linola varieties (2%) [11, 12].
The size of the flax genome (686 Mbp) is around four times that of Arabidopsis and recently different research teams have developed reverse genetics and genomic approaches to learn more about fibre and seed formation in this economicly-important species [13–18]. A cDNA custom flax array system has recently been constructed and successfully used to characterize genes expression profiles in fibre-bearing stem tissues [17, 18]. This platform was constructed by spotting 9,600 anonymous cDNA clones obtained from a flax stem-peel cDNA library. In this paper we report the development and validation of a flax-specific high-density oligo-microarray platform using Nimblegen technology. Flax cDNAs generated from different tissues and/or developmental stages from two different flax cultivars (oil-seed and fibre flax respectively) were 454-sequenced and assembled into contigs. These contigs representing genes found in both oil-seed and fibre flax were used to produce a 48K array. The high quality of our array was demonstrated by the high reproducibility of different technical replications, as well as by the platform's capacity to identify differential gene expression profiles in different tissues and flax varieties.
Results and discussion
454 Sequencing and contig assembly
Details of flax samples used for GS FLX sequencing and for microarray validation experiments.
Green capsule phase
Green capsule phase
454 GS FLX sequencing data for the 9 samples
Bases after clip #
After cleaning 881,950 reads were assembled into 59,494 contigs, with 132 singlets giving a total of 59,626 unigenes http://urgi.versailles.inra.fr/index.php/urgi/Species/Flax/Dowload-sequences. The N50 was 569 bp and 17,487 contigs were longer than 500 bp. 17,731 contigs had only 2 reads, >15,000 contigs contained more than 10 reads and >8,500 contigs contained more than 20 reads. The number of reads varied from 1 to more than 3,500 potentially reflecting the relative abundance of the different transcripts. Incorporation of Sanger sequences available in the public domain did not alter our assembly suggesting that the shorter reads generated by 454 sequencing did not represent a problem for correct assembly (data not shown).
Functional Characterization, GO notation
Numbers and percentages of flax unigenes showing blast hits against different databases
Number of flax unigene hits
Analyses of the 59K unigene set allowed for several different probe design possibilities. After examination, we chose a design based upon 48K contigs as being the most representative of the 59K unigene set. Microarrays were designed according to the Roche/Nimblegen protocol based on 3' end cDNA synthesis and both annotated and non-annotated contigs were used. The final design contained 8 non-overlapping oligos (25-30 mers) per contig allowing for both qualitative and quantitative hybridization. This design is available at the GEO web site http://www.ncbi.nlm.nih.gov/projects/geo/ under the platform number GPL104.
Parameters of microarray performance
The Nimblegen array system is based upon the hybridization of a single labeled sample (derived from RNA), followed by one-channel detection. The intensity of the hybridization signal is then used to determine target concentration. In order to check the technical quality of each probe in our flax Nimblegen arrays, we performed 18 independent hybridizations on nine samples representing different tissues and/or different developmental stages (Table 1) with two technical repetitions. In order to avoid any potential bias due to biological variability the hybridizations were performed on the same RNA samples (Table 1) used for 454 sequencing. Technical quality was determined by evaluating the following parameters: hybridization quality, precision and accuracy.
Our results (GEO GSE21868) showed that all probes present on the array were capable of hybridizing successfully (signal > background) when tested on the different samples. The sensitivity of the array was demonstrated by the wide signal dynamic range obtained (log2 values of 4 to 16). The experimental metrics report (NimbleScan v2.5) was used as recommended by Roche/Nimblegen to generate summary statistics that can be used to identify any potential problems during hybridization. Most of these metrics (interquartile density, signal range, uniformity mean, uniformity CV (coefficient of variation), number of empty features on the array, mean empty, the number of random control features present on the array, mean Random) assume probe randomization on the array surface. Thus, deviation from uniformity (i.e. outside the recommended value range) across the array could suggest potential artifacts during hybridization. Metric values (data not shown) for all 18 samples analyzed were within the recommended value range indicating that hybridization quality was satisfactory.
Accuracy describes how close to a true value a measurement lies. It can be estimated in experiments where a number of realistic targets are spiked at known concentrations into relevant RNA populations, or from comparisons with validation experiments, or by correlation of two different methods .
Comparison of expression levels (log2ratio) from qRT-PCR and microarray for selected target genes
Chlorophyll a-b binding protein 3C-like
RuBisCO activase 2
DNA-directed RNA polymerase II
Ubiquitin carboxyl- terminal hydrolase
Fasciclin-like AGP 4
Differential gene expression detection
In order to see whether our flax microarray platform was able to generate biologically-useful information, we analyzed its capacity to detect specific gene expression profiles associated with different tissues, different developmental stages, and different genotypes. Firstly, we used principal component analysis (PCA) to compare expression profiles in the 9 different samples previously hybridized. Secondly, we compared the expression of 4 genes known to be involved in the process of secondary cell wall biosynthesis and wood formation. Thirdly, we directly compared expression profiles between inner and outer stem tissues at two developmental stages. Finally, we compared gene expression profiles in 2 supplementary flax varieties that show differences in fibre quality and disease resistance.
Differential gene expression in flax tissues and developmental stages
Overall, these results are in general agreement with the known physiological processes of the different organs/tissues thereby suggesting that our platform is capable of generating biologically-useful gene-expression data. For example, leaves are well known to be the plant's main photosynthetic organs and leaf structure is closely associated with its photosynthetic function and chloroplast biogenesis [23–31]. The observation that the GO process 'photosynthesis' is entirely composed of genes specifically expressed in the leaves (Figure 5A), as are the GO processes 'cellular component biogenesis', and 'protein modification process' is in agreement with the biological structure and role of this organ. In root tissues the GO processes 'cell growth' and 'response to stress' were the most represented functions in agreement with the biological activities of this organ [32, 33]. As in other plant species, the flax stem is characterized by the differentiation of different specialized tissues [18, 34] and the majority stem-specific gene expression in the development GO process is in keeping with the biological activity of this organ. Finally, cell-cycle and cellular homeostasis GO processes were associated with embryo-specific gene expressions characteristic of embryogenesis .
Similar analyses (Figure 5B) of the relationship between different tissue-specific genes and metabolisms also suggest that our platform is capable of generating biologically-relevant information. Of interest are the relatively high percentages of genes associated with lipid and carbohydrate metabolisms in seeds since direct quantitative measurements  have shown that sucrose and lipids represent the largest metabolite pools in flax seeds. Similarly, secondary metabolism represents the major metabolism in flax stems and is most likely associated with the lignification of xylem secondary cell walls .
Differential gene expression in a targeted process
Differential gene expression between flax inner and outer stem tissues
In order to verify that our oligo-array platform could be used in future studies to provide biologically-relevant information about cell wall formation and fibre development in flax, we compared expression profiles between the heavily lignified inner stem tissues and the outer stem tissues enriched in weakly-lignified bast fibres at two different developmental stages.
The results (Additional files 1 and 2) showed that in vegetative stage samples, 203 genes were over-expressed in inner stem tissues as compared to outer stem tissues (cut-off: log2ratio (SIV/SOV) > 2), and 229 genes were over-expressed in outer stem tissues as compared to inner stem tissues. Examination of the 203 genes over-expressed in the inner stem tissues showed that 56 genes (28%) were potentially associated with cell wall formation and xylem identity. Of these differentially-expressed cell wall genes, 50% (28 genes) were associated with lignification and included both monolignol biosynthetic genes and laccases involved in the oxidative polymerisation of this phenolic polymer. Such an observation is in agreement with the fact that flax inner stem tissues, mainly composed of xylem, are much more heavily lignified than the outer tissues containing cellulose-rich bast fibes . Other cell-wall related genes coded for proteins involved in polysaccharide cell wall polymer synthesis and remodelling. In addition, 3 genes (C1109, C3396, C5556) coding for fasciclin-like arabinogalactan proteins (FLAs) were also over-expressed in inner tissues as compared to outer tissues. Interestingly, FLA genes are highly expressed in flax outer stem tissues, as well as in poplar tension wood where they have been hypothesized to be involved in the formation of cellulose-rich gelatinous fibres (g-fibres) [17, 37, 38]. Our observation (see below) that 3 other FLA genes (C2947, C3576, C5237) are more highly expressed in flax outer stem tissues when compared to inner stem tissues could suggest that these proteins play a role in secondary cell wall formation in both flax inner and outer stem tissues.
Similar analyses (Additional files 1 and 2) of the 229 genes over-expressed in flax outer stem tissues when compared to inner stem tissues revealed that only 16 genes (7%) were associated with cell wall formation. In contrast, 25 genes were potentially associated with lipid and wax metabolism, 8 genes were related to photosynthesis and 8 genes were stress-related. These differences reflect the different physiological status of flax inner and outer stem tissues. For example, the high percentage of 'cell-wall-related' transcripts in inner stem tissues is associated with secondary xylem formation while the relative abundance of photosynthesis-related transcripts reflects the fact that outer stem tissues are green. Nevertheless, it is interesting to note that 6 genes (C1423, C6325, C353, C1410, C8886, C5052) over-expressed in outer stem tissues corresponded to Lipid Transfer Proteins (LTPs). LTPs have been previously associated with diverse aspects of cell wall development and formation and have also been shown to be more highly expressed in flax outer stem tissues [17, 39]. Such an observation could suggest that such proteins might be associated with fibre maturation in flax. As indicated above, 3 FLAs (C2947, C3576, C5237) were also over-expressed in outer stem tissues as compared to inner stem tissues.
Comparison (Additional files 1 and 2) of inner and outer stem tissue expression profiles (cell-wall related genes) at the green capsule stage revealed a very similar differential expression pattern to that observed for the vegetative stage suggesting that little modification/evolution of inner and outer stem transcriptomes occurs between these two developmental stages. This was confirmed by the observation that only 3 cell-wall related genes (C25634, C2733, C11945) were up-regulated in inner stem tissues at the vegetative stage as compared to the green capsule stage, and no cell-wall related genes were up-regulated in inner stem tissues at the green capsule stage as compared to vegetative inner stem tissues. In outer stem tissues, 8 cell wall genes (C31544, C37539, C6312, C59350, C6820, C44241, C5939, C51183) were more highly expressed at the vegetative stage as compared to the green capsule stage. Two of these genes (C44241, C5939) code for xyloglucan endotransglycosylase/hydroalses (XTHs), one gene (C6820) codes for a secondary cell wall associated glycosyltransferase and one gene for a FLA. Since these 4 genes are implicated in cell wall remodelling and assembly events, such an observation could suggest that cell wall formation in outer stem tissues is more active at the vegetative stage as compared to the green capsule stage. No cell wall genes were more highly expressed at the green capsule stage as compared to the vegetative stage.
The vegetative stage (7-8 weeks) corresponds to the fast-growth stage associated with the increase in plant height and fibre lengthening above the snap point, as well as fibre cell wall thickening below the snap point . At the green capsule stage (10-11 weeks) fibres have stopped elongating and further fibre development is restricted to continued fibre cell wall thickening [40, 41]. Our results would suggest that continued fibre thickening (as well as continued secondary xylem formation in flax inner stem tissues) is associated with continued expression of cell wall related genes.
Altogether these analyses confirm that our oligo-array platform represents a powerful tool for investigating cell wall development and fibre formation in flax stems. Our arrays should also prove extremely useful for investigating other interesting biological processes in flax such as oil and seed formation, disease resistance and tolerance to abiotic stress.
Differential gene expression in contrasting flax genotypes
We used our platform to see whether the flax microarrays were able to detect differentially-expressed genes between 2 contrasting flax genotypes - Drakkar and Belinka. Drakkar produces better quality fibres than the variety Belinka, as well as being more resistant to the fungal pathogen Fusarium http://www.lin-itl.com . RNA was isolated at the green capsule stage from outer stem tissues of field-grown flax plants and 3 biological and 2 technical repetitions were used for each variety. Following hybridization, the Log2ratio was calculated as the intensity signal of each gene in the Drakkar genotype vs. Belinka. The Pearson's correlation was > 0.99 for technical repetitions and > 0.97 for biological repetitions for both Belinka and Drakkar. A total of 428 up-regulated genes and 367 down-regulated genes were found. Only 269 genes were annotated or associated with a known biological function. (Additional files 3 and 4).
A total of 14 genes involved in biotic stress response were up-regulated in Drakkar as compared with Belinka, while 9 genes were up-regulated in Belinka as compared with Drakkar (Figure 7B). Interestingly, 4 of the most up-regulated genes (c31122, c59522, c26764, c402) in Drakkar code for cysteine proteases known to play an important role in programmed cell death (PCD) associated with the hypersensitive response (HR) [48, 49]. Two other highly up-regulated genes (c7149 and c48400) encode germin-like proteins potentially associated with ROS production and cell-wall cross-linking in defense . One up-regulated gene (Nbi-D - c37826) has been previously shown to be involved in flax resistance to rust (Melampsora lini) . In Belinka, the most highly up-regulated gene (c505) as compared to Drakkar correspond to a hypersensitive-induced response protein.
Overall, these results would suggest that our flax microarray platform is capable of detecting biologically-relevant differential gene expression between contrasted flax varieties. As such, our microarrays represent a powerful tool for identifying candidate genes potentially associated with quality-related polymorphism and therefore represent a valuable contribution to molecular-based plant improvement programs.
We have developed a powerful and robust high-density oligo-microarray platform for transcriptomics in flax. High correlations were consistently obtained with technical repetitions on a wide range of different samples and results were cross-validated using an independent method (qRT-PCR). The platform is capable of high discrimination and can provide biologically-useful information on specific gene expression profiles of different flax tissues, and developmental stages. Initial studies also enabled the identification of specifically-expressed cell wall- and defence-related genes in 2 different flax varieties showing contrasting fibre quality and resistance towards a fungal pathogen. These results indicate that our microarray platform can make a useful contribution towards understanding the genetic basis of plant quality in flax.
Plant material and tissue collection
Linum usitatissimum (Barbara an oil-seed cultivar and Hermes a fibre cultivar) plants were grown in a growth chamber (light/night cycles 16 h (22°C)/8 h (19°C), 50% humidity and light intensity of 400 μE s-1 m-2). Individual flowers (Barbara) were tagged at anthesis. The developing seeds of Barbara were collected at 10-15, 20-30 and 40-50 Days After Flowering (DAF), immediately flash-frozen in liquid nitrogen and stored at -80°C until used for experiments. 15 cm-long stem samples were recovered from plants (Hermes) harvested at two different developmental stages: (1) vegetative and (2) green capsule. The outer fibre-bearing tissues were peeled off and inner tissues (xylem) cut into short fragments before both tissues were frozen in liquid nitrogen. Leave samples (Hermes) were collected at vegetative stage, immediately frozen in liquid nitrogen and stored at -80°C. For root samples, plants (Hermes) were grown in a mixture of sand and vermiculites in a growth chamber (light/night cycles 16 h (22°C)/8 h (19°C), 50% humidity and light intensity of 150 μE s-1 m-2) for 10 days. Roots were rinsed in cold water (4°C) and the extremity of each root was cut off, frozen in liquid nitrogen, and stored at -80°C.
Total RNA was isolated from the collected tissues. 100 mg of frozen seed samples were ground in liquid nitrogen and polyphenols and polysaccharides were precipitated  prior to RNA extraction by hot phenol, modified after Verwoerd et al. . 100 mg of frozen stem/root/leaf tissues were ground in liquid nitrogen using the Trizol method (Invitrogen, Carlsbad, CA). Total RNAs were purified using the RNeasy Plant kit (Qiagen) according to the manufacturer's instructions. The genomic DNA was eliminated after treatment with Dnase I for 20 min at 37°C using the DNA-free kit (Ambion, Austin, TX, USA). RNA was checked for purity and degradation by capillary electrophoresis using the Bio-Analyzer Experion (Bio-Rad; RNA Standard Sens kit; RNA StdSens chips). RNA concentration was also determined by spectrometry and only RNAs with an OD260:OD280 ratio of > 1.8 and no discernable degradation were used in PCR-based experiments.
454 Sequencing and Bioinformatics
454 GS FLX technology was used to sequence 9 different RNA samples extracted from 2 distinct flax genotypes and different tissues/developmental stages (Table 1). Individual libraries were prepared for each sample using the MID bar-coding system from Roche. Four separate half GS FLX runs were performed (2 total runs) on 4 cDNA library pools: 1) leaf and root; 2) 3 seed developmental stages; 3) Internal/external stem samples, vegetative stage; 4) internal/external stem samples, green capsule stage. Raw data was generated as Sff files. The Sff_extract tool http://bioinf.comav.upv.es/sff_extract/index.html was used to extract reads, qualities and A and B adaptor positions for each sample. A specific search for SMART forward and reverse primers was performed on all reads using cross_match utility.
Assembly and annotation
Reads from each genotype were identified, and a de novo assembly of all reads/samples was performed using the MIRA tool (development version number 2.9.29 × 4). Parameters used were the default ones plus: -job = denovo, est, normal, 454 -GE:not = 8 -GE:kcim = on - LR:mxti = yes -SB:lsd = on -CL:mbc = on. By default, reads smaller than 40 bp are excluded by MIRA. Reads and contigs were then stored in an EST2uni database for data and annotation management.
Microarray design and oligo synthesis
A total of 384,168 oligonucleotides (25-mers long) were designed and used to construct high-density flax microarrays based on the Nimblegen 385K design format (Nimblegen Systems, Inc., Madison, WI, USA). This design enabled an elevated number (8) of distinct oligos to be used for each of the 48,021 contigs selected from the overall total of 59,000 contigs obtained by assembling the GS FLX sequences. Technical specifications and design files of this new platform for high-throughput analysis of gene expression in flax are publicly available on the GEO website under accession number GPL10419.
cDNA synthesis, labeling and hybridization
Double-stranded cDNA (ds-cDNA) was synthesized from 10 μg of total RNA using an Invitrogen SuperScript ds-cDNA synthesis kit in the presence of 250 ng random hexamer primers. ds-cDNA was cleaned and labeled in accordance with the Nimblegen Gene Expression Analysis protocol (Nimblegen Systems, Inc., Madison, WI, USA). Briefly, ds-cDNA was incubated with 4 μg RNase A (Promega) at 37°C for 10 min and cleaned using phenol:chloroform:isoamyl alcohol, followed by ice-cold absolute ethanol precipitation. For Cy3 labeling of cDNA, the Nimblegen One-Color DNA labeling kit was used according to the manufacturer's guideline detailed in the Gene Expression Analysis protocol (Nimblegen Systems, Inc., Madison, WI, USA). One μg ds-cDNA was incubated for 10 min at 98°C with 2 OD of Cy3-9mer primer. Then, 100 pmol of deoxynucleoside triphosphates and 100U of the Klenow fragment (New England Biolabs, Ipswich, MA, USA) were added and the mix incubated at 37°C for 2h30. The reaction was stopped by adding 0.1 volume of 0.5 M EDTA, and the labeled ds-cDNA was purified by isopropanol/ethanol precipitation. Microarrays were hybridized at 38°C during 16 to 18 h with 6 μg of Cy3 labelled ds-cDNA in Nimblegen hybridization buffer/hybridization component A in a hybridization chamber (Hybridization System - Nimblegen Systems, Inc., Madison, WI, USA). Following hybridization, washing was performed using the Nimblegen Wash Buffer kit (Nimblegen Systems, Inc., Madison, WI, USA).
Slides were scanned at 5 μm/pixel resolution using an Axon GenePix 4000 B scanner (Molecular Devices Corporation, Sunnyvale, CA, USA) piloted by GenePix Pro 6.0 software (Axon). Scanned images (TIFF format) were then imported into NimbleScan software (Nimblegen Systems, Inc., Madison, WI, USA) for grid alignment and expression data analyses. Expression data were normalized through quantile normalization  and the Robust Multichip Average (RMA) algorithm  included in the NimbleScan software. Identification of genes displaying a change in expression over repetitions was accomplished with a script utilizing library functions in R with a false discovery rate (FDR) of less than 5%. The SAM  was used to identify differentially expressed genes over different conditions and log2(ratio) ≥ 1 and ≤ -1 were used for filtering gene expression profiles. Analysis was completed with the PCA module of the J-express program . Functional annotation of differentially-expressed genes was based on Gene Ontology http://www.geneontology.org/. All the microarray data have been submitted to the Gene Expression Omnibus (GEO) database http://www.ncbi.nlm.nih.gov/geo with the accession number GSE21868.
Quantitative reverse transcriptase-PCR (qRT-PCR) analysis
5 μg aliquots of total RNA were treated with DNaseI using a TURBO DNA-free Kit (Ambion) according to the manufacturer's instructions, then first-strand cDNA was synthesized using M-MuLV RNase H- reverse transcriptase (Finnzymes) with 2.5 μg of random hexamers and 500 ng of oligo(dT)12 according to the manufacturer's instructions. The reaction was stopped by incubation at 70°C for 10 min, and the reaction mixture treated with RNaseH (BioLabs) according to the manufacturer's instructions before dilution with 600 μL of sterile de-ionized water.
qRT-PCR Experiment Design
Transcript levels were assessed by qRT-PCR, in assays with triplicate reaction mixtures (final volume, 20 μL) containing 5 μL of cDNA, 0.5 μM of both forward and reverse primers, and 1X DyNamo Capillary SYBR Green qRT-PCR mix (Finnzymes). qRT-PCR experiments used a balanced randomized block design, as recently advised . A LightCycler (Roche) was used to acquire the CT values for each sample. The following standard protocol was applied for all amplifications: 10 min at 95°C, followed by 45 cycles of 10 s at 95°C, 15 s at 60°C, and 15 s at 72°C. A melting curve analysis was added to each PCR program and the size of PCR products was assessed by electrophoresis in agarose gels. The primer sequences used for all target genes are presented in Additional file 5.
Relative standard curves describing the PCR efficiencies (E) for each primer pair were generated for each amplicon according to Larionov et al. . Normalization of qRT-PCR was performed using reference genes (R) according to Gutierrez et al. . Twelve genes (Additional file 5) were chosen for their stability across the set of microarray experiments previously shown in this paper. Their expression was assessed by qRT-PCR, and they were ranked according to their stability of expression using geNorm software . The_c3168 (lcl|genolin_c3168 628 nt) and c10916 (lcl|genolin_c10916 813 nt) genes were the most stably expressed ones among the 12 tested and, thus, were used to normalize the qRT-PCR data. The normalized expression patterns obtained using both reference genes were similar, so only the data normalized with c10916 ( highly similar to a Dehydrodolichyl diphosphate synthase 6 of Arabidopsis thaliana) are shown in this article. Gene expression was calculated using CT and E values with the formula, ET(CTtissueA-CTtissueB)/ER((CTtissueA-CTtissueB), where (T) is the target gene and (R) the reference gene, (tissueA) is related to cDNA from the tissue showing the lowest expression and (tissueB) from the tissue showing the highest expression. All qRT-PCR results represent means as calculated from the three technical replicates .
This work was financially supported by ANR (ANR-06-GPAL-017: Genolin project) and by the Pôle of competitivity Industry and Agro-Ressources. SF et YN gratefully acknowledge the financial support of the French national research agency (ANR: Genolin project). SF also acknowledges the financial support of the Nord Pas-de-Calais Region (Plant Teq project 4). BT thanks the EU for the financial contribution towards equipments acquired in the operational regional competitivity program.
- Deyholos MK: Bast fiber of flax (Linum usitatissimum L.): Biological foundations of its ancient and modern uses. Israel Journal of Plant Sciences. 2006, 54 (4): 273-280. 10.1560/IJPS_54_4_273.View ArticleGoogle Scholar
- Vaisey-Genser M, Morris DH: History of the cultivation of flaxseed. Flax, the genus Linum. Edited by: Muir AD, Wescott ND. 2003, CRC Press, Taylor & Francis Group, Florida, 1-21.Google Scholar
- Baley C, Perrot Y, Busnel F, Guezenoc H, Davies P: Transverse tensile behaviour of unidirectional plies reinforced with flax fibres. Materials letters. 2006, 60 (24): 2984-2987. 10.1016/j.matlet.2006.02.028.View ArticleGoogle Scholar
- Kymäläinen HR, Sjöberga AM: Flax and hemp as raw materials for thermal insulations. Building and Environment. 2008, 43: 1261-1269. 10.1016/j.buildenv.2007.03.006.View ArticleGoogle Scholar
- Farmer C, Petit HV, Weiler H, Capuco AV: Effects of dietary supplementation with flax during prepuberty on fatty acid profile, mammogenesis, and bone resorption in gilts. Journal of Animal Science. 2007, 85 (7): 1675-1686. 10.2527/jas.2007-0022.PubMedView ArticleGoogle Scholar
- Dodin S, Cunnane SC, Masse B, Lemay A, Jacques H, Asselin G, Tremblay-Mercier J, Marc I, Lamarche B, Legare F: Flaxseed on cardiovascular disease markers in healthy menopausal women: a randomized, double-blind, placebo-controlled trial. Nutrition. 2008, 24 (1): 23-30.PubMedView ArticleGoogle Scholar
- Felmlee MA, Woo G, Simko E, Krol ES, Muir AD, Alcorn J: Effects of the flaxseed lignans secoisolariciresinol diglucoside and its aglycone on serum and hepatic lipids in hyperlipidaemic rats. The British Journal of Nutrition. 2009, 102 (3): 361-369. 10.1017/S0007114508207488.PubMedView ArticleGoogle Scholar
- Chen J, Soucek MD, Simonsick WJ, Celikay RW: Synthesis and photopolymerization of norbornyl epoxidized linseed oil. Polymer. 2002, 43: 5379-5389. 10.1016/S0032-3861(02)00404-4.View ArticleGoogle Scholar
- Hall C, Tulbek MC, Xu Y: Flaxseed. Adv Food Nutr Res. 2006, 51: 1-97. full_text.PubMedView ArticleGoogle Scholar
- Day A, Neutelings G, Nolin F, Grec S, Habrant A, Cronier D, Maher B, Rolando C, David H, Chabbert B: Caffeoyl coenzyme A O-methyltransferase down-regulation is associated with modifications in lignin and cell-wall architecture in flax secondary xylem. Plant Physiology & Biochemistry. 2009, 47 (1): 9-19.View ArticleGoogle Scholar
- Gutierrez L, Conejero G, Castelain M, Guenin S, Verdeil JL, Thomasset B, Van Wuytswinkel O: Identification of new gene expression regulators specifically expressed during plant seed maturation. Journal of Experimental Botany. 2006, 57 (9): 1919-1932. 10.1093/jxb/erj138.PubMedView ArticleGoogle Scholar
- Rao S, Abdel-Reheem M, Bhella R, McCracken C, Hildebrand D: Characteristics of high alpha-linolenic acid accumulation in seed oils. Lipids. 2008, 43 (8): 749-755. 10.1007/s11745-008-3207-0.PubMedView ArticleGoogle Scholar
- Bennett MD, Leitch IJ: Plant DNA C-values database (release 3.0). 2004, [http://www.rbgkew.org.uk/cval/homepage.html]Google Scholar
- Day A, Addi M, Kim W, David H, Bert F, Mesnage P, Rolando C, Chabbert B, Neutelings G, Hawkins S: ESTs from the fibre-bearing stem tissues of flax (Linum usitatissimum L.): expression analyses of sequences related to cell wall development. Plant Biology. 2005, 7 (1): 23-32. 10.1055/s-2004-830462.PubMedView ArticleGoogle Scholar
- De Pauw MA, Vidmar JJ, Collins J, Bennett RA, Deyholos MK: Microarray analysis of bast fibre producing tissues of Cannabis sativa identifies transcripts associated with conserved and specialised processes of secondary wall development. Functional Plant Biology. 2007, 34 (8): 737-749. 10.1071/FP07014.View ArticleGoogle Scholar
- Hotte NS, Deyholos MK: A flax fibre proteome: identification of proteins enriched in bast fibres. BMC Plant Biology. 2008, 8: 52-10.1186/1471-2229-8-52.PubMed CentralPubMedView ArticleGoogle Scholar
- Roach MJ, Deyholos MK: Microarray analysis of flax (Linum usitatissimum L.) stems identifies transcripts enriched in fibre-bearing phloem tissues. Molecular Genetics & Genomics. 2007, 278 (2): 149-165.View ArticleGoogle Scholar
- Roach MJ, Deyholos MK: Microarray analysis of developing flax hypocotyls identifies novel transcripts correlated with specific stages of phloem fibre differentiation. Annals of Botany. 2008, 102 (3): 317-330. 10.1093/aob/mcn110.PubMed CentralPubMedView ArticleGoogle Scholar
- Moreau Y, Aerts S, De Moor B, De Strooper B, Dabrowski M: Comparison and meta-analysis of microarray data: from the bench to the computer desk. Trends in Genetics. 2003, 19 (10): 570-577. 10.1016/j.tig.2003.08.006.PubMedView ArticleGoogle Scholar
- Nimgaonkar A, Sanoudou D, Butte AJ, Haslett JN, Kunkel LM, Beggs AH, Kohane IS: Reproducibility of gene expression across generations of Affymetrix microarrays. BMC Bioinformatics. 2003, 4: 27-10.1186/1471-2105-4-27.PubMed CentralPubMedView ArticleGoogle Scholar
- Hardiman G: Microarray platforms--comparisons and contrasts. Pharmacogenomics. 2004, 5 (5): 487-502. 10.1517/14622418.104.22.1687.PubMedView ArticleGoogle Scholar
- Wang Y, Barbacioru C, Hyland F, Xiao W, Hunkapiller KL, Blake J, Chan F, Gonzalez C, Zhang L, Samaha RR: Large scale real-time PCR validation on gene expression measurements from two commercial long-oligonucleotide microarrays. BMC genomics. 2006, 7: 59-10.1186/1471-2164-7-59.PubMed CentralPubMedView ArticleGoogle Scholar
- Buick R: When did oxygenic photosynthesis evolve?. Philosophical Transactions of the Royal Society of London B Biological Science. 2008, 363 (1504): 2731-2743. 10.1098/rstb.2008.0041.View ArticleGoogle Scholar
- Campbell NA, Williamson B, Heyden RJ: Biology: Exploring Life. 2006, Boston, Masschussets: Pearson Prentice HallGoogle Scholar
- Douglas SE: Plastid evolution: origins, diversity, trends. Current Opinion in Genetics & Development. 1998, 8 (6): 655-661.View ArticleGoogle Scholar
- Gould SB, Waller RF, McFadden GI: Plastid evolution. Annual Review of Plant Biology. 2008, 59: 491-517. 10.1146/annurev.arplant.59.032607.092915.PubMedView ArticleGoogle Scholar
- Nealson KH, Conrad PG: Life: past, present and future. Philosophical Transactions of the Royal Society of London B Biological Science. 1999, 354 (1392): 1923-1939. 10.1098/rstb.1999.0532.View ArticleGoogle Scholar
- Raven PH, Evert RF, Eichhorn SE: Biology of Plants, 7th Edition. 2005, New York: W.H. Freeman and Company PublishersGoogle Scholar
- Reyes-Prieto A, Weber AP, Bhattacharya D: The origin and establishment of the plastid in algae and plants. Annual Review of Genetics. 2007, 41: 147-168. 10.1146/annurev.genet.41.110306.130134.PubMedView ArticleGoogle Scholar
- Rodriguez-Ezpeleta N, Brinkmann H, Burey SC, Roure B, Burger G, Loffelhardt W, Bohnert HJ, Philippe H, Lang BF: Monophyly of primary photosynthetic eukaryotes: green plants, red algae, and glaucophytes. Current Biology. 2005, 15 (14): 1325-1330. 10.1016/j.cub.2005.06.040.PubMedView ArticleGoogle Scholar
- Smith AL: Oxford dictionary of biochemistry and molecular biology. 1997, Oxford: Oxford University PressGoogle Scholar
- Bengough AG, Bransby MF, Hans J, McKenna SJ, Roberts TJ, Valentine TA: Root responses to soil physical conditions; growth dynamics from field to cell. Journal of Experimental Botany. 2006, 57 (2): 437-447. 10.1093/jxb/erj003.PubMedView ArticleGoogle Scholar
- Okubara PA, Paulitz TC: Root defense responses to fungal pathogens: A molecular perspective. Plant and Soil. 2005, 274 (1-2): 215-226. 10.1007/s11104-004-7328-9.View ArticleGoogle Scholar
- Day A, Ruel K, Neutelings G, Cronier D, David H, Hawkins S, Chabbert B: Lignification in the flax stem: evidence for an unusual lignin in bast fibers. Planta. 2005, 222 (2): 234-245. 10.1007/s00425-005-1537-1.PubMedView ArticleGoogle Scholar
- Dupree P: Plant embryogenesis: cell division forms a pattern. Current Biology. 1996, 6 (6): 683-685. 10.1016/S0960-9822(09)00449-7.PubMedView ArticleGoogle Scholar
- Troufflard S, Roscher A, Thomasset B, Barbotin JN, Rawsthorne S, Portais JC: In vivo 13C NMR determines metabolic fluxes and steady state in linseed embryos. Phytochemistry. 2007, 68 (16-18): 2341-2350. 10.1016/j.phytochem.2007.04.018.PubMedView ArticleGoogle Scholar
- Mellerowicz EJ, Sundberg B: Wood cell walls: biosynthesis, developmental dynamics and their implications for wood properties. Current Opinion in Plant Biology. 2008, 11 (3): 293-300. 10.1016/j.pbi.2008.03.003.PubMedView ArticleGoogle Scholar
- Larfarguette F, Leplé JC, Déjardin A, Laurens F, Costa G, Lesage-Descauses MC, Pilate: Poplar genes encoding fasciclin-like arabinogalactan proteins are highly expressed in tension wood. New Phytol. 2004, 164: 107-121. 10.1111/j.1469-8137.2004.01175.x.View ArticleGoogle Scholar
- Yeats TH, Rose JKC: The biochemistry and biology of extracellular plant lipid-transfer proteins (LTPs). Protein Sci. 2008, 17: 191-198. 10.1110/ps.073300108.PubMed CentralPubMedView ArticleGoogle Scholar
- Gorshkova TA, Sal'nikov VV, Chemikosova SB, Ageeva MV, Pavlencheva NV, van Dam JEG: The snap point: a transition point in Linum usitatissimum bast fiber development, Ind. Crops and Prod. 2003, 18: 213-221. 10.1016/S0926-6690(03)00043-8.View ArticleGoogle Scholar
- His I, Andème-Onzighi C, Morvan C, Driouich A: Microscopic studies on mature flax fibers embedded in LR White: immunogold localization of cell wall matrix polysccharides, 2001. J Histochem and Cytochem. 2001, 49 (12): 1525-1535.View ArticleGoogle Scholar
- Vromans J: Molecular genetic studies in flax (Linum usitatissimum L.). pHD thesis. 2006, Wageningen University, The NetherlandsGoogle Scholar
- Morvan C, Andème-Onzighi C, Girault R, Himmelsbach DS, Driouich A, Akin DE: Building flax fibres: more than one brick in the walls. Plant Physiology and Biochemistry. 2003, 41 (11-12): 935-944. 10.1016/j.plaphy.2003.07.001.View ArticleGoogle Scholar
- Tang GQ, Luscher M, Sturm A: Antisense repression of vacuolar and cell wall invertase in transgenic carrot alters early plant development and sucrose partitioning. Plant Cell. 1999, 11 (2): 177-189. 10.1105/tpc.11.2.177.PubMed CentralPubMedView ArticleGoogle Scholar
- Gorshkova T, Gurjanov O, Mikshina P, Ibragimova N, Mokshina N, Salnikov V, Ageeva M, Amenitskii S, Chernova T, Chemikosova S: Specific type of secondary cell wall formed by plant fibers. Russian Journal of Plant Physiology. 2010, 57 (3): 328-341. 10.1134/S1021443710030040.View ArticleGoogle Scholar
- Hobson N, Roach M, Deyholos M: Gene expression in tension wood and bast fibres. Russian Journal of Plant Physiology. 2010, 57 (3): 321-327. 10.1134/S1021443710030039.View ArticleGoogle Scholar
- Kaku T, Serada S, Baba Ki, Tanaka F, Hayashi T: Proteomic analysis of the G-layer in poplar tension wood. Journal of Wood Science. 2009, 55 (4): 250-257. 10.1007/s10086-009-1032-6.View ArticleGoogle Scholar
- Woltering EJ: Death proteases: alive and kicking. Trends in Plant Science. 2010, 15 (4): 185-188. 10.1016/j.tplants.2010.02.001.PubMedView ArticleGoogle Scholar
- Zhang L, Xu Q, Xing D, Gao C, Xiong H: Real-time detection of caspase-3-like protease activation in vivo using fluorescence resonance energy transfer during plant programmed cell death induced by ultraviolet C overexposure. Plant Physiology. 2009, 150 (4): 1773-1783. 10.1104/pp.108.125625.PubMed CentralPubMedView ArticleGoogle Scholar
- Manosalva PM, Davidson RM, Liu B, Zhu X, Hulbert SH, Leung H, Leach JE: A germin-like protein gene family functions as a complex quantitative trait locus conferring broad-spectrum disease resistance in rice. Plant Physiology. 2009, 149 (1): 286-296. 10.1104/pp.108.128348.PubMed CentralPubMedView ArticleGoogle Scholar
- Islam R, Sheperd KW, Mayo GME: Effect of genotype and temperature on the expression of L genes in flax conferring resistance to rust. Physiology and Molecular Plant Pathology. 1989, 35 (2): 141-150. 10.1016/0885-5765(89)90084-2.View ArticleGoogle Scholar
- Gehrig HH, Winter K, Cushman JC, Borland AM, Taybi T: An improved RNA isolation method for succulent plant species rich in polyphenols and polysaccharides. Plant Molecular Biology Report. 2000, 18: 369-376. 10.1007/BF02825065.View ArticleGoogle Scholar
- Verwoerd TC, Dekker BM, Hoekema A: A small-scale procedure for the rapid isolation of plant RNAs. Nucleic Acids Research. 1989, 17 (6): 2362.-10.1093/nar/17.6.2362.PubMed CentralPubMedView ArticleGoogle Scholar
- Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19 (2): 185-193. 10.1093/bioinformatics/19.2.185.PubMedView ArticleGoogle Scholar
- Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4 (2): 249-264. 10.1093/biostatistics/4.2.249.PubMedView ArticleGoogle Scholar
- Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proceedings of the National Academy of Sciences of the USA. 2001, 98 (9): 5116-5121. 10.1073/pnas.091062498.PubMed CentralPubMedView ArticleGoogle Scholar
- Dysvik B, Jonassen I: J-Express: exploring gene expression data using Java. Bioinformatics. 2001, 17 (4): 369-370. 10.1093/bioinformatics/17.4.369.PubMedView ArticleGoogle Scholar
- Rieu I, Powers SJ: Real-time quantitative RT-PCR: design, calculations, and statistics. Plant Cell. 2009, 21 (4): 1031-1033. 10.1105/tpc.109.066001.PubMed CentralPubMedView ArticleGoogle Scholar
- Larionov A, Krause A, Miller W: A standard curve based method for relative real time PCR data processing. BMC Bioinformatics. 2005, 6: 62-10.1186/1471-2105-6-62.PubMed CentralPubMedView ArticleGoogle Scholar
- Gutierrez L, Mauriat M, Pelloux J, Bellini C, Van Wuytswinkel O: Towards a systematic validation of references in real-time rt-PCR. Plant Cell. 2008, 20 (7): 1734-1735. 10.1105/tpc.108.059774.PubMed CentralPubMedView ArticleGoogle Scholar
- Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biology. 2002, 3 (7): RESEARCH0034-10.1186/gb-2002-3-7-research0034.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.