Selection of reference genes for quantitative gene expression normalization in flax (Linum usitatissimum L.)
© Huis et al; licensee BioMed Central Ltd. 2010
Received: 14 July 2009
Accepted: 19 April 2010
Published: 19 April 2010
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).
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.
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.
Flax is one of the earliest cultivated plants known to man and the remains of seeds found in different archeological sites suggest that it was also used before domestication, most probably for textiles . Cultivated flax (Linum usitatissimum L.) is an annual, diploid, self-pollinated crop used as a source of both high-quality, cellulose-rich bast fibers and oil. The main areas of fiber production are currently found in France (77,000 ha) with an annual turnover of over 200 million euros followed by Russia and China (80,000 ha). In contrast, oil varieties (linseed) are mainly grown in India (930,000 ha), Canada (811,000 ha) and China (570,000 ha) .
Although flax fibers have now been largely replaced in textiles by cotton or synthetic fibers, they still are used in the fabrication of high-quality linen material and are also incorporated in the polymeric matrix during the fabrication of biocomposites in order to improve their mechanical properties .
In the flax stem, fibers are derived from the procambial cells of the protophloem  and are referred to as bast fibers. Individual bast fibers are extremely-long (≤ 77 mm) single cells  characterized by the presence of a very thick secondary cell-wall containing high amounts of cellulose and very low amounts of lignins as compared to the shorter fibers found in the xylem . Bast fibers provide mechanical support to the plant and are organized in bundles that occupy the great majority of the external tissues between the epidermis and the vascular cambium. Flax fiber extraction is initiated during the retting process when stems are left on the ground in order to promote microbiological-mediated separation of the fiber bundles from the surrounding tissues .
Secondary cell-wall formation is most often studied in woody species such as poplar, pine and eucalyptus, or else in model species such as Arabidopsis and tobacco. More recently, growing interest in the unusual structure of the secondary cell wall of bast fibers (high crystalline cellulose, low lignin) is stimulating transcriptomic and genomic studies in flax and other bast fiber species [8–10]. A better understanding of the molecular regulation controlling fiber cell-wall metabolism should increase our knowledge of the very complex coordination necessary for the production and assembly of different polymers within the cell wall.
Flax is a suitable model for genetic and functional genomics as it is autogamous, easy to grow, has a relatively short vegetative stage, and can be genetically transformed [11–13]. The first studies on flax DNA structure were undertaken almost thirty years ago [14, 15] and the first cDNA sequence was deposited in Genbank in 1993 . Previous analyses of the flax transcriptome have mainly involved studies of single gene expression by Northern blot [17–23] or multiple gene expression by classical RT-PCR . Only relatively few studies [10, 25] have reported the use of qRT-PCR in flax. Quantitative real-time RT-PCR (qRT-PCR) is currently the most accurate method for detecting low abundant mRNAs and is capable of detecting slight variations in gene expression in different tissues of the same plant . When compared to traditional methods used to evaluate transcript accumulation, the main advantages of qRT-PCR are its very high sensitivity and specificity. However, in order to obtain reliable and reproducible results, it is necessary to include suitable internal controls. Indeed, the high sensitivity of qRT-PCR can lead to misinterpretation of expression data , especially when biological samples show similar gene expression levels. Slight variations in the RNA integrity or quantity, in the reverse transcription yield, or more generally in the transcriptional activity of the studied organ can all have a major impact on quantification. It is therefore absolutely necessary to develop a normalisation strategy before undertaking gene expression studies. The identification of suitable control genes is thus an important challenge in transcriptome analyses and the major aim is to identify genes ubiquitously expressed in every tissue independently of the experimental context .
In the past a number of different genes were commonly used for normalising expression, especially in classical RT-PCR approaches. Most of these genes proved to be suitable for such studies at the time, mainly because of the low sensibility of RT-PCR. More recently, the same reference genes have also been used for qRT-PCR and it has been assumed that their expression is perfectly stable even though this has not necessarily been experimentally verified. Examples of such commonly-used genes include GAPDH (glyceraldehyde-3-phosphate dehydrogenase) , cyclophilin and 18S rRNA , actin , EF1alpha , ubiquitin  and beta tubulin .
Surprisingly, there are relatively few reports concerning the study of stably-expressed genes, referred to as housekeeping genes (HKGs), in different plant species, even though qRT-PCR is routinely used in many research projects. Suitable HKGs have been identified in model species with sequenced genomes such as rice , Arabidopsis , grapevine  and poplar , as well as in economically-important crops such as wheat , barley , Coffee , tomato , potato  and soybean . In most cases, gene expression stability was determined by using statistical approaches such as geNorm  and NormFinder . The geNorm algorithm allows the identification of the most suitable reference gene(s) and the optimal number of genes that should be used. It relies on a pairwise comparison and evaluates the variation of relative quantity ratios for each gene pair in a set of expression data. The NormFinder algorithm depends on a statistical and mathematical model that not only estimates the overall expression variation of a candidate gene, but also considers the variation between the chosen subgroups.
In the present study, we have evaluated the expression profiles of 13 putative HKGs during the development of flax plants with a special focus on the stem which contains two major cell-wall tissue types: 1) inner tissues (lignin-rich secondary cell walls of xylem) and 2) outer tissues (lignin-poor, cellulose-rich secondary cell walls of bast fibers). The two algorithms geNorm and NormFinder were used in order to determine the best reference genes needed for normalisation. The selected HKGs were then used to normalise gene expression in flax for an investigation of the expression of LuMYB1 - a flax MYB transcription factor highly similar to AtMYB59 .
MYB transcription factors are of particular interest in studies on secondary cell wall development since they have been shown to regulate the expression of genes coding enzymes of the phenylpropanoid biosynthetic pathway responsible for the production of lignin monomers [48, 49]. We decided to investigate the expression profiles of LuMYB1 since initial transcriptome studies in our laboratory (data not shown) had indicated that the gene was most highly expressed in inner stem tissues and therefore potentially associated with secondary cell wall formation in flax xylem tissues. In addition, AtMYB59 has been shown to undergo alternative splicing leading to different transcripts in rice and Arabidopsis, and therefore represents an interesting model for gene regulation studies.
Selection of candidate reference genes and amplification specificity
Flax candidate HKG description and comparison with Arabidopsis orthologs.
GenBank accession number
Arabidopsis ortholog locus
Arabidopsis BlastX E-value
Elongation Factor 1-α
Elongation Factor 2
Eukaryotic translation initiation Factor 1
Eukaryotic translation initiation Factor 3 E
Eukaryotic translation initiation Factor 3 H
Eukaryotic translation initiation Factor 4 E
Eukaryotic translation initiation Factor 5 A
Glyceraldehyde 3-phosphate dehydrogenase
Ubiquitin extension protein
Description of candidate reference genes, primers and amplicons for flax housekeeping gene selection.
Amplicon length (bp)
Elongation Factor 1-α
Elongation Factor 2
Eukaryotic translation initiation Factor 1
Eukaryotic translation initiation Factor 3 E
Eukaryotic translation initiation Factor 3 H
Eukaryotic translation initiation Factor 4 E
Eukaryotic translation initiation Factor 5 A
Glyceraldehyde 3-phosphate dehydrogenase
Ubiquitin extension protein
Expression profile of the reference genes
With regards to gene expression variation, TUA showed the highest value (6.98 cycles) and ETIF4E showed the lowest (3.67 cycles) when all tissues are considered. As might be expected, expression variation was lower when only stem samples are considered. Once again TUA showed the highest value (6.1 cycles) whereas the lowest value was observed for EF1A (3.2 cycles). These results indicated that none of the selected genes showed a near constant expression level and it was therefore necessary to evaluate the best candidates for gene expression normalization.
Expression stability analysis
Since the 13 candidate HKGs show wide variations in expression levels in different flax tissues, it is necessary to use statistical approaches to rank the expression levels and determine the number of housekeeping genes necessary for accurate gene-expression profiling in the selected tissues. We decided to use the two most widely-used algorithms, geNorm and NormFinder.
We then used the NormFinder algorithm  to identify the most stable gene among our candidates and test whether the use of only one HKG would be suitable for gene expression studies. The NormFinder algorithm requires the analysis of a minimum of 3 genes and a minimum of 2 samples per group. It can analyze the expression data obtained from quantitative methods and choose the best normalization gene from a set of candidates. The software calculates and ranks the stability values for candidate genes within the sample set investigated according to their expression profile. The lowest stability value indicates the most stably expressed gene.
Ranking of candidate genes according to NormFinder.
Best combination of two genes
Interestingly, both the rice and flax (but not Arabidopsis) type 2 transcripts contain a stop codon close to (and in frame with) the corresponding ATG start codon in Arabidopsis. As a result in both rice and flax, translation is initiated at the next ATG that corresponds to the functional ATG in type 1 transcripts. In consequence, the flax type 1 and 2 ORFs, as in rice, are identical. We have not yet detected a flax type 3 transcript (LuMYB1-3) in our ESTs even though sequence data predicts the existence of such a transcript. Nevertheless, we designed qRT-PCR primers to specifically amplify the two longest flax transcripts. LuMYB1-1 was detected with the M1R primer located in the intron, and LuMYB1-2 with the M2F primer overlapping the 3' splicing site (Figure 5B).
Over the last years several very powerful techniques have been developed to detect differences in gene expression levels between different cell types, tissues and organs. Among these, qRT-PCR is now commonly used in many laboratories to undertake accurate expression profiling of different candidate genes that have been previously identified, either because they are known to be involved in a specific biological process, or else because they have shown an interesting expression profile in global expression analyses (microarrays).
However, reliable relative quantification can only be performed if an accurate normalization is performed following the choice and verification of suitable (i.e. stable) HKGs. Although microarray data have been analyzed to determine the most stably-expressed genes in Arabidopsis  and rice , it is not always possible to directly transpose suitable HKGs identified in one species to another species. For example UBQ10 is very stable in Arabidopsis  but not in rice , nor in soybean  nor in Brachypodium . Therefore, since there are no universally-suitable reference genes, it is necessary to verify the expression levels of the candidate reference genes under the same experimental conditions as those used for the gene of interest.
In this study, we have measured the expression levels and stability of candidate genes in different tissues at three developmental stages in flax plants. Our major research interest in this species concerns the cell wall formation and development of the phloem (bast) fibers located in the external tissues of the stem. The polymer composition of the secondary cell walls in these long fiber cells (as well as that of the shorter xylem cells located in the inner tissues) has been previously analyzed  and we have shown that the abundance and structure of the phenolic polymer lignin varies according to the age and physiological state of the plant. In order to further investigate cell wall formation (including lignification) in flax, we are currently undertaking gene expression profiling by qRT-PCR on different tissues. However, due to important differences between the structure and the metabolic state of cells in flax inner and outer tissues, it is necessary to identify the most relevant combination of reference genes. Normalization of gene expression in flax has previously been done with either a homolog to an Arabidopsis NADH Ubiquinone oxidoreductase determined by cDNA-AFLP , or with an Elongation Factor EF1A homolog .
We originally selected 20 candidate HKGs from the literature and after determination of amplification efficiencies for the selected primers we retained 13 candidates. We then evaluated their normalization potential using 2 commonly-used algorithms (geNorm  and NormFinder ) in a systematic study of their expression stability in different flax tissues including stems, flowers, roots and leaves. The NormFinder algorithm  generally ranked GAPDH as the most stably expressed gene in all the flax samples irrespective of whether the samples were assembled into 1 main group or divided into several sub-groups. In only 1 case did NormFinder rank GAPDH as the second most stable gene. These results are consistent with those observed recently in Brachypodium distachyon  during cold/heat stress and in Coffea arabica  in different organs and tissues. In contrast to NormFinder, the geNorm algorithm indicated that we should use 3 HKGs (EFIA, ETIF5A, Ubiquitin) for studies involving a mix of root, leaf, stem and flower tissues. However, the same program indicated that only 2 genes (ETIF1, ETIF4F) should be used when just stem tissues were considered. While EF1A is often described as a stable gene and has been used as a reference in many species [37, 43, 52, 54], there are only relatively few reports of Eucaryotic Translation Initiation Factor genes being used as reference genes. For example in soybean, ETIF1A and ETIF1B genes were ranked as the most stable genes tested when 21 samples were pooled, or when photoperiodic treatments were modified . In contrast, expression patterns of several translation initiation factors were shown to be unstable in wheat , and have been excluded as good candidate reference genes in poplar  and Darnel ryegrass .
We also evaluated the stability of tubulin and actin genes that are very often used to normalize expression data. Although actin has been shown to be a suitable normalization gene in developmental studies , it also appears to be unstable in many biological processes . In the same way, tubulin has been shown to be stable during development in orobranche , but is apparently unstable (geNorm analyses) during development  and abiotic stress . Our results would also suggest that these 2 genes are among the most instable for expression profiling during development in flax.
Altogether, our results showed that the geNorm and NormFinder algorithms came to different conclusions concerning the best candidate HKG(s) to use for expression normalization in flax. In order to see whether the use of these different HKGs modified the determined expression profiles (and hence possible interpretations of biological role) of a gene of interest, we analyzed the expression of LuMYB1 in different tissues and organs of Flax. LuMYB1 is a homolog of AtMYB59, a gene that has been shown to undergo alternative splicing in Arabidopsis and rice . We decided to investigate the expression profile of LuMYB1 since our team is interested in gaining a better understanding of the regulation of cell wall formation in flax and different MYB transcription factors have been previously shown to play an important role in controlling this process . In addition, preliminary transcriptome data in our lab had indicated that LuMYB1 was highly expressed in stem inner tissues. Such an observation could suggest that this gene is involved in the transcriptional regulation of cell-wall genes during the formation of xylem.
Our results showed that the 2 different flax spliceforms (LuMYB1-1 and LuMYB1-2) are strongly expressed in the inner stem tissues, but only expressed at low levels in external stem tissues, leaves and roots. In addition, only LuMYB1-2 is expressed in flowers. Such a result is interesting since previous work  has shown that the corresponding Arabidopsis gene is mainly expressed in roots, leaves and seedlings, but only poorly expressed in stem tissues. More recently , AtMYB59 has been shown to be specifically expressed during the S phase of the cell cycle in Arabidopsis cell suspensions, suggesting that this gene plays a key role in cell cycle regulation. Our observation that the corresponding flax gene (both spliceforms) are highly expressed in inner stem tissues could also suggest a role for LuMYB1 in cell cycle events since these flax tissues are the site of intense mitotic activity associated with the division of vascular cambium initials to form new xylem cells. However, further work is obviously necessary to confirm this hypothesis.
In addition, our results also revealed that both the choice of HKGs and the decision or not to analyze samples together or in sub-groups modified the determined expression profiles for the gene of interest LuMYB1. Although such choices did not affect the identity of the 4 (stem) tissues showing the highest expression levels, they did affect the ranking of these different tissues. Such an observation suggests that care should be taken when interpreting the biological significance of small differences in gene expression, and that in this case, the researchers should probably consider using different normalizations to verify their data. In addition, it is clear that functional approaches are also necessary to fully understand the role of a given gene in a particular biological process. Nevertheless, our results showing important differences in LuMYB1 expression levels (regardless of the normalization type) would suggest that this gene is potentially associated with the formation of xylem tissue in flax stems.
A considerable quantity of transcriptional data is currently available for major model plant species and in silico analyses can therefore be used to identify suitable HKGs for gene expression normalization in these species. However, for most other plant species, the suitability of potential HKGs identified in the literature or in heterologous databases must be verified by qRT-PCR. In this study, we have identified several suitable reference genes for studying developmental gene expression in flax, with a special focus on stem tissues. The use of HKGs identified by both geNorm and NormFinder algorithms allowed us to determine the expression profiles of a gene of interest (LuMYB1) in a large range of different flax tissues. Our results also showed that certain classically-used reference genes such as actin and tubulin are not necessarily the most suitable for studies of quantitative gene expression in flax. In conclusion, our study has identified suitable HKGs for future developmental transcriptome studies in flax. These genes also represent potentially interesting targets for normalizing gene expression in flax under stress conditions, although it would obviously be necessary to verify their stabilities in this case.
Flax plants (Linum usitatissinum L. cv Barbara were grown in a greenhouse under 16 h/20°C day and 8 h/18°C night conditions. The plants were harvested at three developmental stages: (1) vegetative (56 days after sowing), (2) flowering (131 days after sowing, 50% open flowers) and (3) green capsule stage (159 days after sowing). Following harvest, organs were dissected and the stems were divided into three equal parts: - apical, medium and basal (the latter was discarded because the secondary cell walls of bast fibers have obtained their maximum thickness in this part of the stem and residual metabolism is low and unrelated to fiber maturation events).
For the flowering and green capsule stages, the outer fiber-bearing tissues were separated from the inner woody tissues. The tissues were immediately frozen in liquid nitrogen and stored at -80°C.
RNA isolation, quality control and cDNA synthesis
Frozen tissues were ground in liquid nitrogen using a mortar and a pestle. Total RNA was extracted using the RNeasy Plant Mini Kit (Qiagen) according to the manufacturer's instructions. RNA purity was assessed on a biophotometer (Eppendorf) by determining the OD260/OD280 and OD260/OD230 ratios, which were between 1.8 and 2. Potentially-contaminating DNA was eliminated by treatment with DNAse I using the DNA-free kit (Ambion). RNA concentration and quality were determined by capillary electrophoresis on an Experion labchip electrophoresis system (Bio-Rad) One microgram of total RNA was reverse-transcribed using the Iscript cDNA synthesis kit (Bio-Rad) according to the manufacturer's instructions. The cDNAs were diluted 1:256 with nuclease free water.
Design and validation of reference gene primers
Specific primer pairs were first designed for 20 commonly used housekeeping genes representing distinct functional classes and gene families. These include ubiquitins, actin, tubulin, elongation factors, cytosolic cylophilins and translational initiation factors (Table 2). The sequence of these genes were identified by BLAST searches (Table 1) in a flax EST database obtained from a cDNA library derived from the outer fiber-bearing tissues of flax . For primer design, Primer3 software  was used (160 bp maximum length, optimal Tm at 60°C, GC % between 20% and 80%). The absence of secondary structure in the primer annealing fragments was verified using mfold . Five-point standard curves of a 4-fold dilution series (1:16 to 1:4096) were used to calculate the PCR efficiency (E) of each primer pair. The PCR efficiency is given by the equation E = (10(-1/m) -1) × 100 where m is the slope of the linear regression model fitted over log transformed data of the input cDNA concentration versus Ct values according to the linear equation y = m log(x) + b. After evaluation of the 20 primer pairs, 13 were retained on the basis of their amplification efficiency value (93-105%).
qRT-PCR conditions and analyses
The qRT-PCRs were carried out in 96-wells plates with a MyIQ real time PCR detection system (Bio-Rad) using Quantitect SYBR Green PCR Kit (Qiagen) in a reaction volume of 20 μL (5 μL diluted cDNAs, 10 μL of 2× SYBR Green mix and primer pairs at 0.4 μM). Aliquots from the same cDNA solutions were used with all primer sets in each experiment. All PCR reactions were performed under the following conditions: 95°C for 15 min, 40 cycles of 10 s at 95 °C and 30 s at 60°C. For each primer pair, a melting curve was generated in order to confirm the specificity of the amplification and the PCR products were checked on a 4% agarose gel. Data were analysed using Bio-Rad iQ5 software. The efficiency (E) value of each reaction was between 1.93 and 2.04 with R2 values higher than 0.998. Each experiment was repeated three times on two biological replicates, each one represented by three technical repetitions. PCR reactions on samples lacking the cDNA template or the reverse transcriptase during the cDNA synthesis were also performed as negative controls for each primer pair.
Statistical analyses of gene expression stability
The stability of the candidate reference genes was evaluated by two statistical approaches. In both approaches, expression levels were expressed relative to the sample with the highest expression. Ct values were converted into relative quantities and imported into the geNorm v3,5 software http://medgen.ugent.be/~jvdesomp/genorm/ and into the NormFinder software http://www.mdl.dk/publicationsnormfinder.htm.
The geNorm algorithm first calculates an expression stability value (M) for each gene and then the pairwise variation (V) of this gene with the others. All the tested genes are ranked according to their stability in the tested tissues and the number of HKGs necessary for an optimal normalization is indicated. NormFinder also ranks the stability of the tested genes, but independently of each other.
Expression of LuMYB1
The authors gratefully acknowledge S. Fénart for his helpful advice on statistics. RH is supported by the French Ministère de la Recherche. This work was financed by the French Nord-Pas de Calais regional project ARCir Plant Teq 4: 'Identification des determinants moléculaires de la qualité des fibres de lin'.
- Zohary D, Hopf M: Domestication of plants in the old world Third edition. Oxford University Press, 2004,Google Scholar
- ITL: [http://www.lin-itl.com/]
- Bodros E, Pillin I, Montrley N, Baley C: Could biopolymers reinforced by randomly scattered flax fibre be used in structural applications?. Composites Science and Technology. 2007, 67: 462-470. 10.1016/j.compscitech.2006.08.024.View ArticleGoogle Scholar
- Gorshkova TA, Wyatt SE, Salnikov VV, Gibeaut DM, Ibragimov MR, Lozovaya VV, Carpita NC: Cell-Wall Polysaccharides of Developing Flax Plants. Plant Physiol. 1996, 110: 721-729.PubMedPubMed CentralGoogle Scholar
- Reddy N, Yang Y: Biofibers from agricultural byproducts for industrial applications. Trends in Biotechnology. 2005, 23: 22-27. 10.1016/j.tibtech.2004.11.002.PubMedView 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: 234-245. 10.1007/s00425-005-1537-1.PubMedView ArticleGoogle Scholar
- Akin DE, Morrison WH, Rigsby LL, Dodd RB: Plant factors influencing enzyme retting of fiber and seed flax. J Agric Food Chem. 2001, 49: 5778-5784. 10.1021/jf010804d.PubMedView ArticleGoogle Scholar
- Wrobel-Kwiatkowska M, Starzycki M, Zebrowski J, Oszmianski J, Szopa J: Lignin deficiency in transgenic flax resulted in plants with improved mechanical properties. J Biotechnol. 2007, 128: 919-934. 10.1016/j.jbiotec.2006.12.030.PubMedView ArticleGoogle Scholar
- Ebskamp MJ: Engineering flax and hemp for an alternative to cotton. Trends Biotechnol. 2002, 20: 229-230. 10.1016/S0167-7799(02)01953-4.PubMedView ArticleGoogle Scholar
- Roach MJ, Deyholos MK: Microarray analysis of flax (Linum usitatissimum L.) stems identifies transcripts enriched in fibre-bearing phloem tissues. Mol Genet Genomics. 2007, 278: 149-165. 10.1007/s00438-007-0241-1.PubMedView ArticleGoogle Scholar
- Lamblin F, Saladin G, Dehorter B, Cronier D, Grenier E, Lacoux J, Bruyant P, Laine E, Chabbert B, Girault F, et al: Overexpression of a heterologous sam gene encoding S-adenosylmethionine synthetase in flax (Linum usitatissimum) cells: Consequences on methylation of lignin precursors and pectins. Physiol Plant. 2001, 112: 223-232. 10.1034/j.1399-3054.2001.1120211.x.PubMedView ArticleGoogle Scholar
- Vrinten P, Hu Z, Munchinsky MA, Rowland G, Qiu X: Two FAD3 desaturase genes control the level of linolenic acid in flax seed. Plant Physiol. 2005, 139: 79-87. 10.1104/pp.105.064451.PubMedPubMed CentralView ArticleGoogle Scholar
- Caillot S, Rosiau E, Laplace C, Thomasset B: Influence of light intensity and selection scheme on regeneration time of transgenic flax plants. Plant Cell Rep. 2009, 28: 359-371. 10.1007/s00299-008-0638-2.PubMedView ArticleGoogle Scholar
- Cullis CA: DNA sequence organization in the flax genome. Biochim Biophys Acta. 1981, 652: 1-15.PubMedView ArticleGoogle Scholar
- Goldsbrough PB, Cullis CA: Characterisation of the genes for ribosomal RNA in flax. Nucleic Acids Res. 1981, 9: 1301-1309. 10.1093/nar/9.6.1301.PubMedPubMed CentralView ArticleGoogle Scholar
- Song W, Funk C, Brash A: Molecular cloning of an allene oxide synthase: a cytochrome P450 specialized for the metabolism of fatty acid hydroperoxides. Proc Natl Acad Sci USA. 1993, 90: 8519-8523. 10.1073/pnas.90.18.8519.PubMedPubMed CentralView ArticleGoogle Scholar
- Dickinson MJ, Zhang R, Pryor A: Nucleotide sequence relationships of double-stranded RNAs in flax rust, Melampsora lini. Curr Genet. 1993, 24: 428-432. 10.1007/BF00351852.PubMedView ArticleGoogle Scholar
- Omann F, Beaulieu N, Tyson H: cDNA sequence and tissue-specific expression of an anionic flax peroxidase. Genome. 1994, 37: 137-147. 10.1139/g94-018.PubMedView ArticleGoogle Scholar
- Roberts JK, Pryor A: Isolation of a flax (Linum usitatissimum) gene induced during susceptible infection by flax rust (Melampsora lini). Plant J. 1995, 8: 1-8. 10.1046/j.1365-313X.1995.08010001.x.PubMedView ArticleGoogle Scholar
- Day A, Dehorter B, Neutelings G, Czeszak X, Chabbert B, Belingheri L, David H: Caffeoyl-coenzyme A 3-O-methyltransferase enzyme activity, protein and transcript accumulation in flax (Linum usitatissimum) stem during development. Physiol Plant. 2001, 113: 275-284. 10.1034/j.1399-3054.2001.1130216.x.PubMedView ArticleGoogle Scholar
- Wrobel M, Zebrowski J, Szopa J: Polyhydroxybutyrate synthesis in transgenic flax. J Biotechnol. 2004, 107: 41-54. 10.1016/j.jbiotec.2003.10.005.PubMedView ArticleGoogle Scholar
- Harms K, Atzorn R, Brash A, Kuhn H, Wasternack C, Willmitzer L, Pena-Cortes H: Expression of a Flax Allene Oxide Synthase cDNA Leads to Increased Endogenous Jasmonic Acid (JA) Levels in Transgenic Potato Plants but Not to a Corresponding Activation of JA-Responding Genes. Plant Cell. 1995, 7: 1645-1654. 10.1105/tpc.7.10.1645.PubMedPubMed CentralView 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: 749-755. 10.1007/s11745-008-3207-0.PubMedView ArticleGoogle 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 Biol (Stuttg). 2005, 7: 23-32. 10.1055/s-2004-830462.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. J Exp Bot. 2006, 57: 1919-1932. 10.1093/jxb/erj138.PubMedView ArticleGoogle Scholar
- Gachon C, Mingam A, Charrier B: Real-time PCR: what relevance to plant studies?. J Exp Bot. 2004, 55: 1445-1454. 10.1093/jxb/erh181.PubMedView ArticleGoogle Scholar
- Guenin S, Mauriat M, Pelloux J, Van Wuytswinkel O, Bellini C, Gutierrez L: Normalization of qRT-PCR data: the necessity of adopting a systematic, experimental conditions-specific, validation of references. J Exp Bot. 2009, 60: 487-493. 10.1093/jxb/ern305.PubMedView ArticleGoogle Scholar
- Sturzenbaum SR, Kille P: Control genes in quantitative molecular biological techniques: the variability of invariance. Comp Biochem Physiol B Biochem Mol Biol. 2001, 130: 281-289. 10.1016/S1096-4959(01)00440-7.PubMedView ArticleGoogle Scholar
- Magneschi L, Kudahettige RL, Alpi A, Perata P: Expansin gene expression and anoxic coleoptile elongation in rice cultivars. J Plant Physiol. 2009, 166 (14): 1576-80. 10.1016/j.jplph.2009.03.008.PubMedView ArticleGoogle Scholar
- Brentner LB, Mukherji ST, Merchie KM, Yoon JM, Schnoor JL, Van Aken B: Expression of glutathione S-transferases in poplar trees (Populus trichocarpa) exposed to 2,4,6-trinitrotoluene (TNT). Chemosphere. 2008, 73: 657-662. 10.1016/j.chemosphere.2008.07.059.PubMedView ArticleGoogle Scholar
- Bracha-Drori K, Shichrur K, Lubetzky TC, Yalovsky S: Functional analysis of Arabidopsis postprenylation CaaX processing enzymes and their function in subcellular protein targeting. Plant Physiol. 2008, 148: 119-131. 10.1104/pp.108.120477.PubMedPubMed CentralView ArticleGoogle Scholar
- Bomal C, Bedon F, Caron S, Mansfield SD, Levasseur C, Cooke JE, Blais S, Tremblay L, Morency MJ, Pavy N, et al: Involvement of Pinus taeda MYB1 and MYB8 in phenylpropanoid metabolism and secondary cell wall biogenesis: a comparative in planta analysis. J Exp Bot. 2008, 59: 3925-3939. 10.1093/jxb/ern234.PubMedPubMed CentralView ArticleGoogle Scholar
- Aquea F, Gutierrez F, Medina C, Arce-Johnson P: A novel Otubain-like cysteine protease gene is preferentially expressed during somatic embryogenesis in Pinus radiata. Mol Biol Rep. 2008, 35: 567-573. 10.1007/s11033-007-9124-0.PubMedView ArticleGoogle Scholar
- Diretto G, Welsch R, Tavazza R, Mourgues F, Pizzichini D, Beyer P, Giuliano G: Silencing of beta-carotene hydroxylase increases total carotenoid and beta-carotene levels in potato tubers. BMC Plant Biol. 2007, 7: 11-10.1186/1471-2229-7-11.PubMedPubMed CentralView ArticleGoogle Scholar
- Kim BR, Nam HY, Kim SU, Kim SI, Chang YJ: Normalization of reverse transcription quantitative-PCR with housekeeping genes in rice. Biotechnol Lett. 2003, 25: 1869-1872. 10.1023/A:1026298032009.PubMedView ArticleGoogle Scholar
- Czechowski T, Stitt M, Altmann T, Udvardi MK, Scheible WR: Genome-wide identification and testing of superior reference genes for transcript normalization in Arabidopsis. Plant Physiol. 2005, 139: 5-17. 10.1104/pp.105.063743.PubMedPubMed CentralView ArticleGoogle Scholar
- Reid KE, Olsson N, Schlosser J, Peng F, Lund ST: An optimized grapevine RNA isolation procedure and statistical determination of reference genes for real-time RT-PCR during berry development. BMC Plant Biol. 2006, 6: 27-10.1186/1471-2229-6-27.PubMedPubMed CentralView ArticleGoogle Scholar
- Brunner AM, Yakovlev IA, Strauss SH: Validating internal controls for quantitative plant gene expression studies. BMC Plant Biol. 2004, 4: 14-10.1186/1471-2229-4-14.PubMedPubMed CentralView ArticleGoogle Scholar
- Paolacci AR, Tanzarella OA, Porceddu E, Ciaffi M: Identification and validation of reference genes for quantitative RT-PCR normalization in wheat. BMC Mol Biol. 2009, 10: 11-10.1186/1471-2199-10-11.PubMedPubMed CentralView ArticleGoogle Scholar
- Faccioli P, Ciceri GP, Provero P, Stanca AM, Morcia C, Terzi V: A combined strategy of "in silico" transcriptome analysis and web search engine optimization allows an agile identification of reference genes suitable for normalization in gene expression studies. Plant Mol Biol. 2007, 63: 679-688. 10.1007/s11103-006-9116-9.PubMedView ArticleGoogle Scholar
- Barsalobres-Cavallari CF, Severino FE, Maluf MP, Maia IG: Identification of suitable internal control genes for expression studies in Coffea arabica under different experimental conditions. BMC Mol Biol. 2009, 10: 1-10.1186/1471-2199-10-1.PubMedPubMed CentralView ArticleGoogle Scholar
- Exposito-Rodriguez M, Borges AA, Borges-Perez A, Perez JA: Selection of internal control genes for quantitative real-time RT-PCR studies during tomato development process. BMC Plant Biol. 2008, 8: 131-10.1186/1471-2229-8-131.PubMedPubMed CentralView ArticleGoogle Scholar
- Nicot N, Hausman JF, Hoffmann L, Evers D: Housekeeping gene selection for real-time RT-PCR normalization in potato during biotic and abiotic stress. J Exp Bot. 2005, 56: 2907-2914. 10.1093/jxb/eri285.PubMedView ArticleGoogle Scholar
- Jian B, Liu B, Bi Y, Hou W, Wu C, Han T: Validation of internal control for gene expression study in soybean by quantitative real-time PCR. BMC Mol Biol. 2008, 9: 59-10.1186/1471-2199-9-59.PubMedPubMed CentralView 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 Biol. 2002, 3: RESEARCH0034-10.1186/gb-2002-3-7-research0034.PubMedPubMed CentralView ArticleGoogle Scholar
- Andersen CL, Jensen JL, Orntoft TF: Normalization of real-time quantitative reverse transcription-PCR data: a model-based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data sets. Cancer Res. 2004, 64: 5245-5250. 10.1158/0008-5472.CAN-04-0496.PubMedView ArticleGoogle Scholar
- Li J, Li X, Guo L, Lu F, Feng X, He K, Wei L, Chen Z, Qu LJ, Gu H: A subgroup of MYB transcription factor genes undergoes highly conserved alternative splicing in Arabidopsis and rice. J Exp Bot. 2006, 57: 1263-1273. 10.1093/jxb/erj094.PubMedView ArticleGoogle Scholar
- Zhong R, Ye ZH: Regulation of cell wall biosynthesis. Curr Opin Plant Biol. 2007, 10: 564-572. 10.1016/j.pbi.2007.09.001.PubMedView ArticleGoogle Scholar
- Zhou J, Lee C, Zhong R, Ye ZH: MYB58 and MYB63 are transcriptional activators of the lignin biosynthetic pathway during secondary cell wall formation in Arabidopsis. Plant Cell. 2009, 21: 248-266. 10.1105/tpc.108.063321.PubMedPubMed CentralView ArticleGoogle Scholar
- Gorshkova TA, Sal'nikov VV, Chemikosova SB, Ageeva MV, Pavlencheva NV, van Dam J: The snap point: a transition point in Linum usitatissimum bast fiber development. Industrial Crops and Products. 2003, 18: 213-221. 10.1016/S0926-6690(03)00043-8.View ArticleGoogle Scholar
- Jain M, Khurana JP: Transcript profiling reveals diverse roles of auxin-responsive genes during reproductive development and abiotic stress in rice. FEBS J. 2009, 276: 3148-3162. 10.1111/j.1742-4658.2009.07033.x.PubMedView ArticleGoogle Scholar
- Jain M, Nijhawan A, Tyagi AK, Khurana JP: Validation of housekeeping genes as internal control for studying gene expression in rice by quantitative real-time PCR. Biochem Biophys Res Commun. 2006, 345: 646-651. 10.1016/j.bbrc.2006.04.140.PubMedView ArticleGoogle Scholar
- Hong SY, Seo PJ, Yang MS, Xiang F, Park CM: Exploring valid reference genes for gene expression studies in Brachypodium distachyon by real-time PCR. BMC Plant Biol. 2008, 8: 112-10.1186/1471-2229-8-112.PubMedPubMed CentralView ArticleGoogle Scholar
- Martin R, Hollenbeck V, Dombrowski J: Evaluation of Reference Genes for Quantitative RT-PCR in Lolium perenne. Crop Science. 2008, 48: 1881-1887. 10.2135/cropsci2007.10.0597.View ArticleGoogle Scholar
- Gallie DR, Le H, Caldwell C, Browning KS: Analysis of translation elongation factors from wheat during development and following heat shock. Biochem Biophys Res Commun. 1998, 245: 295-300. 10.1006/bbrc.1998.8427.PubMedView ArticleGoogle Scholar
- Gonzaìlez-Verdejo C, Die J, Nadal S, Jimeìnez-Mariìn A, Moreno M, Romaìn B: Selection of housekeeping genes for normalization by real-time RT-PCR: Analysis of Or-MYB1 gene expression in Orobanche ramosa development. Analytical Biochemistry. 2008, 379: 38-43.Google Scholar
- Mu RL, Cao YR, Liu YF, Lei G, Zou HF, Liao Y, Wang HW, Zhang WK, Ma B, Du JZ, et al: An R2R3-type transcription factor gene AtMYB59 regulates root growth and cell cycle progression in Arabidopsis. Cell Res. 2009, 19: 1291-1304. 10.1038/cr.2009.83.PubMedView ArticleGoogle Scholar
- Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000, 132: 365-386.PubMedGoogle Scholar
- Zuker M: Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 2003, 31: 3406-3415. 10.1093/nar/gkg595.PubMedPubMed CentralView 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.