<title>Background</title> <p>Legumes are second only to grasses in agricultural importance [##REF##12644639##1##]. They are a mainstay of sustainable agricultural systems because of their ability to reduce atmospheric nitrogen (N<sub>2</sub>) to ammonia via a symbiosis with bacteria called rhizobia. This provides legumes and subsequent crops with a free and renewable source of nitrogen in lieu of expensive, environmentally-unfriendly fertilizers. Development and differentiation of root nodules, the organ that accommodates nitrogen-fixing rhizobia in legumes, is orchestrated by transcription factors [##REF##10647012##2##, ####REF##10691739##3##, ##REF##12397181##4##, ##REF##12856938##5##, ##REF##15961668##6##, ##REF##15961669##7##, ##REF##17114582##8##, ##REF##17449807##9####17449807##9##]. Transcription factors are DNA-binding proteins that regulate the transcription of most, if not all genes [##REF##11118137##10##]. As a result, TFs play central roles in all aspects of plant biology, including development and differentiation of organs and adaptive responses to changes in the environment [##REF##17556517##11##]. Transcription factors as a whole are an important target of plant research because they are a key to understanding the regulation of important plant processes as well as potential tools to optimize these processes for agriculture.</p> <p>The importance of TFs in plant biology is reflected by the fact that approximately 5% of all plant genes encode such proteins [##REF##11118137##10##]. Thus, even species with relatively small genomes, such as <italic>Arabidopsis thaliana </italic>contain thousands of TF genes [##REF##11118137##10##]. This presents a real challenge for systematic approaches to decipher the function of TF genes in plants. Classical, 'forward' genetics has uncovered the roles of perhaps a hundred TF genes in Arabidopsis [##UREF##0##12##] and far fewer in other species [##REF##17556517##11##]. Reverse-genetic approaches, using T-DNA insertion mutants for instance [##REF##16755288##13##], provide a means to decipher in a systematic and relatively rapid manner the function of TF genes/proteins, although gene-redundancy often stymies this enterprise [##UREF##0##12##]. Another stumbling-block is that phenotypes associated with non-redundant TFs may be subtle in nature.</p> <p>Transcript profiling can help to uncover the functions of TF genes/proteins by revealing where and when in a plant TF genes are expressed. This information can help direct our attention to particular organs, developmental stages, or conditions under which aberrant phenotypes might become apparent in a TF mutant of interest.</p> <p><italic>Medicago truncatula </italic>is a model legume species that is currently the focus of an international genome sequencing effort [##REF##15824279##14##]. Several generations of cDNA [##REF##15466239##15##] and oligonucleotide arrays [##REF##15778460##16##] have been developed for transcriptome analysis of <italic>Medicago truncatula</italic>, including most-recently an Affymetrix GeneChip that contains 51,000 probe-sets representing a large proportion of all the genes in this species [##REF##18302802##17##]. While these tools now provide a means to measure the transcriptional output of a large proportion of genes in Medicago, inherent limitations in the sensitivity of hybridization-based technologies [##REF##12192590##18##] mean that transcripts of a substantial number of genes cannot be detected even when probes for these transcripts are present on the array/chip. Furthermore, expansion of arrays to encompass novel genes uncovered by genome sequencing is not a trivial task. An alternative to arrays that is 2–3 orders of magnitude more sensitive and more flexible in terms of expansion to encompass novel genes is quantitative reverse transcription-polymerase chain reaction (qRT-PCR). Platforms for qRT-PCR analysis of thousands of Arabidopsis and rice TF genes have been developed by us and others [##REF##15078338##19##,##REF##17559651##20##], and utilized to identify TF genes involved in Arabidopsis responses to nutrient stress and pathogen attack [##REF##15375205##21##, ####REF##16183832##22##, ##REF##17177879##23##, ##REF##17722694##24####17722694##24##]. Here we describe a bioinformatics pipeline to identify putative TF genes in <italic>Medicago truncatula </italic>and to design gene-specific oligonucleotide primers for qRT-PCR analysis of all predicted TF transcripts. Over 1000 TF primer pairs were tested and used to identify sets of organ-enhanced TF genes that may play important roles in organ development or differentiation in this species.</p>
<title>Methods</title> <title>Plant material and growth conditions</title> <p><italic>Medicago truncatula </italic>cv. Jemalong, line A17 wild type plants were vernalized for 3 days in the dark at 4°C on sterile, wet filter paper. Germinated seedlings were transferred to pots containing Turface (BWI Texarcana, Texarcana, TX). Plants were grown in growth chambers under a 16 h day and 8 h night regime, at 200 μE light intensity, 24°C and 40% relative humidity.</p> <p>Vegetative organs (leaves, stems, roots, and nodules) were harvested 28 days after planting. Leaf material did not include petioles and stems did not include buds. Roots consisted of the entire root system with laterals. Several plants grown at the same time were pooled for each of the three biological replicates. Biological replicates were planted on separate days. Nodules were harvested from plants inoculated with <italic>Sinorhizobium meliloti </italic>strain 1021 one and seven days after sowing. Reproductive organs were harvested from plants that were vernalized for two weeks to decrease the time to flowering. Flowers were harvested on the day of opening. Pods were harvested from 1 to 21 days after the appearance of the floral bud to cover a wide range of developmental stages. Harvested plant material was frozen in liquid nitrogen before storage at -80°C.</p> <title>RNA isolation and cDNA synthesis</title> <p>Total RNA was extracted using Trizol reagent [##REF##8747660##36##], following the manufacturer's instructions (Invitrogen GmbH, Karlsruhe, Germany). RNA was quantified using a Nanodrop Spectrophotometer ND-100 (NanoDrop Technologies, Wilington, DE). Sixty μg of total RNA were digested with RNase free DNase1 (Ambion Inc., Houston, TX), according to manufacturer's protocol. RNA integrity was checked using an Agilent 2100 BioAnalyser and RNA 6000 Nano Chips (Agilient Technologies, Waldbronn, Germany), and by electrophoresis on a 3% (v/w) agarose gel before and after DNase I treatment. The absence of contaminating genomic DNA after DNase I treatment was verified by PCR analysis, using primer pairs designed to amplify a 107 bp genomic fragment of the control gene, <italic>Ubiquitin </italic>(TC102473intronF, 5'-GTCCTCTAAGGTTTAATGAACCGG-3'; TC102473intronR, 5'-GAAAGACACAGCCAAGTTGCAC-3').</p> <p>First-strand complementary DNA was synthesized by priming with oligo-dT<sub>12–18 </sub>(Qiagen, Hilden, Germany), using SuperScript III reverse transcriptase (Invitrogen GmbH, Karlsruhe, Germany) following the instructions of the provider. To assess cDNA synthesis efficiency, qPCR was used to amplify segments in the 5' and 3' regions of <italic>Ubiquitin </italic>cDNA approx. 1600 and 400 bp from the 3'-end, respectively (primers: TC102473_5'F, 5'-TTGGAGACGGATTCCATTGCT-3'; TC102473_5'R, 5'-GCCAATTCCTTCCCTTCGAA-3; TC102473_3'F, 5'-GGCCCTAGAACATTTCCTGTGG-3'; and TC102473_3'R, 5'-TTGGCAACCAAAATGTTCCC-3'). If ΔCt (Ct3'-Ct5') &lt; 2, then cDNA synthesis efficiency was judged to be satisfactory, and the cDNA was considered suitable for qRT-PCR analysis.</p> <title>PCR primer design</title> <p>The primer design pipeline was implemented in object-oriented PERL modules supported by a MySQL relational database. Primers iterated through three phases before approval: design, specificity, and selection.</p> <p>The design phase interrogated TF genes with a sliding window 250 bp across that stepped 50 bp along the entire target sequence, generating primer candidates at each window. Experimental conditions, as outlined in the Results section above, were enforced by the following MIT Primer3 parameters: PRIMER_MIN_TM 58, PRIMER_OPT_TM 60, PRIMER_MAX_TM 61, PRIMER_SELF_ANY 6, PRIMER_SELF_END 2, PRIMER_MAX_POLY_X 3, and PRIMER_PRODUCT_SIZE_RANGE '100–150' [##UREF##6##37##]. The specificity phase aligned primer candidates via WU-Blast to a database of all known Medicago sequences. The selection phase sorted primer candidates by the number of possible secondary hits, self-complementarity, and poly-X characteristics. Secondary hits were defined as specificity alignments that contained at least one of the terminal ends of the primer and achieved 80% or greater identity over the length of the primer. The sequences of each primer pair are given in Supplementary Material (see Additional file ##SUPPL##0##1##).</p> <title>Real-time PCR conditions and analysis</title> <p>PCR reactions were carried out in an ABI PRISM<sup>® </sup>7900 HT Sequence Detection System (Applied Biosystems, Foster City, CA, USA). SYBR<sup>® </sup>Green was used to quantify dsDNA synthesis. Reactions (5 μl total volume) were performed in an optical 384-well plate containing 2.5 μl 2 × SYBR<sup>® </sup>Green Power Master Mix reagent (Applied Biosystems, Warringen, UK), 5 ng cDNA and 200 n<sc>M</sc> of each gene-specific primer. Primer pairs were aliquoted using a pipetting robot (Evolution P3 liquid handling system, Perkin Elmer, MA, USA) to minimize pipetting errors. cDNA was aliquoted as a master mix of cDNA and 2 × SYBR<sup>® </sup>Green reagent, using an electronic Eppendorf multipipette. Reaction plates were sealed with a transparent adhesive cover before proceeding (Applied Biosystems, Foster City, CA, USA). All templates were amplified using the following standard PCR protocol: 50°C for 2 min; 95°C for 10 min; 40 cycles of 95°C for 15 sec and 60°C for 1 min, and SYBR<sup>® </sup>Green fluorescence was measured continuously. Melting curves were generated after 40 cycles by heating the sample up to 95°C for 15 sec followed by cooling down to 60°C for 15 s and heating the samples to 95°C for 15 sec.</p> <p>Data analysis was performed with the SDS 2.2.1 software (Applied Biosystems). To determine the threshold cycle value (C<sub>T</sub>) for each PCR reaction, the threshold (Δ<italic>R</italic><sub>n</sub>) was set within the logarithmic amplification phase. All amplification plots were analyzed with an Δ<italic>R</italic><sub>n </sub>of 0.2. PCR efficiency (E) was estimated using LinReg software with data obtained from the exponential phase of each individual amplification plot and the equation (1+E) = 10<sup>slope </sup>[##REF##12618301##31##]. To compare data from different PCR runs and different cDNA samples, C<sub>T </sub>values were normalized against the geometric mean of four reference genes (<italic>Ubquitin</italic>, <italic>PPRep</italic>, <italic>PDF2</italic>, and <italic>PTB</italic>), whose transcript levels were most stable across the biological samples analyzed. The average of the geometric mean of these four genes for all 18 samples was C<sub>T </sub>21.23 ± SD1.15. For normalization, the mean reference gene C<sub>T </sub>value was substracted from the C<sub>T </sub>value of the TF gene of interest, yielding a ΔC<sub>T </sub>value. The expression ratios for the identification of organ-enhanced genes were obtained using the following formula on all 30 organ combinations: , where ΔΔC<sub>T </sub>was calculated by ΔC<sub>TA </sub>minus ΔC<sub>Tb</sub>, A and B are averages of three biological replicates of the two organs being compared, and <italic>E </italic>is the PCR efficiency. Dissociation curves were analysed using SDS 2.2.1 software (Applied Biosystems). RT-PCR products were resolved on 3% (w/v) agarose gels (LE Agarose, Biozym, Oldendorf, Germany) run at 4 V cm<sup>-1 </sup>in TAE Tris-Acetate-EDTA buffer, along with a 200-bp DNA-standard ladder (Promega GmbH). A subset of 178 RT-PCR products was sequenced at the JC Venter Institute (Rockville, MD, USA).</p>
<title>Results and Discussion</title> <title>Identification of putative transcription factors</title> <p>TF protein families are generally defined by the type(s) of DNA-binding domain they contain and putative TF genes are often identified on the basis of DNA sequences that encode known DNA-binding domains [##REF##11118137##10##,##REF##17556517##11##,##REF##15731212##25##,##REF##16769687##26##]. We utilized this approach to identify putative TFs of Medicago amongst the set of proteins predicted from genomic sequence by the International Medicago Genome Annotation Group (IMGAG). Proteins of IMGAG release 1, which contained over 40,000 predicted proteins, were screened for the presence of known or presumed DNA-binding domains (Table ##TAB##0##1##), using INTERPRO [##UREF##1##27##]. Medicago proteins containing putative DNA binding domains and other domains associated with TFs were then used as query sequences in WU-BLASTX [##UREF##2##28##] which included searches of both the non-redundant DNA database of NCBI [##UREF##3##29##] and the well-curated protein database, UniProt [##UREF##4##30##] to check annotations of related proteins in support of tentative Medicago TF assignments. This process resulted in a list of 1045 putative TF genes (see Additional file ##SUPPL##0##1##). We utilized genomic sequences rather than the large collection of partial cDNA sequences present in Expressed Sequence Tag (EST) databases for Medicago as the starting point for TF gene discovery because protein sequences derived from genomic sequence are more complete and the set of IMGAG proteins essentially contains no redundancy. Although identification of the 'complete' set of Medicago TFs from IMGAG-annotated proteins will only be possible upon completion of genome sequencing, we expect little or no redundancy in the protein set targeted by our primer set. This approach avoids wasting money on redundant primer sets and re-organization of primers when redundancy is detected, both of which would have been inevitable if we chose to use ESTs in addition to genomic sequences to identify Medicago TFs.</p> <title>PCR primer design</title> <p>To ensure maximum specificity and efficiency during PCR amplification of TF cDNA under a standard set of reaction conditions, a stringent set of criteria was used for primer design. This included predicted melting temperatures (<italic>T</italic><sub>m</sub>) of 58°C to 61°C, limited self-complementarity and poly-X, and PCR amplicon lengths of 100–150 base pairs (bp). Secondary hits were minimized by aligning primer candidates to all known Medicago sequences via WU-BLAST [##UREF##2##28##] and eliminating primer pairs with multiple potential hits.</p> <title>PCR primer testing: gene-specificity and amplification efficiency</title> <p>PCR primers were tested on Medicago cDNA free of genomic DNA contamination as follows. First, total RNA was extracted from various organs using Trizol reagent (Invitrogen GmbH, Karlsruhe, Germany), which yielded high quality RNA as judged by gel electrophoresis and by Agilent 2100 BioAnalyser using RNA 6000 Nano Chips (Agilient Technologies, Waldbronn, Germany). Typical RNA yields ranged from 0.5–1.0 μg RNA/mg fresh mass for nodules and leaves, respectively. Isolated RNA was treated with DNAse I (Ambion, product number 1907) to remove all contaminating genomic DNA, and this was always confirmed by PCR using primers to non-coding regions of the <italic>Ubiquitin </italic>gene (TC102473; AC137828-19.4). After inactivation of DNAse I, RNA was reverse transcribed using SuperScript III reverse transcriptase (Invitrogen GmbH, Karlsruhe, Germany) and oligo-dT<sub>12–18 </sub>to prime the reaction.</p> <p>Specificity of PCR primers was assessed in three ways: by melting curve analysis of PCR reaction products; by separating the products of all reactions via electrophoresis in 3% agarose gels; and by sequencing a sub-set of PCR reaction products (Figure ##FIG##0##1##). 94.5% (998/1045) of primer pairs gave unique PCR products of the expected size. Only 3.3% (34/1045) of primer pairs yielded no product and 2.2% (23/1045) gave non-specific products. (see Additional file ##SUPPL##0##1##). Sequencing was performed on 178 randomly-chosen PCR products amplified from a 1:1 mixture of leaf and root cDNA. In the vast majority of cases (92.7% or 165/178), the sequence of the PCR product was identical to that of the intended target gene. In 5.1% of cases, the amplicon sequence matched multiple related genes, including the target gene, while in only 2.2% of cases the amplicon sequence did not match the target gene sequence.</p> <p>Ideally, PCR results in an exact doubling of the amount of dsDNA after each temperature cycle. In practice, however, this is generally not the case because the reactions are less than 100% efficient. Primer sequences can affect PCR efficiency, so we determined the efficiency of each TF primer pair from amplification plots, using LinRegPCR software [##REF##12618301##31##]. First, the correlation coefficient derived from linear regression analysis of each amplification plot (e.g. see Figure ##FIG##1##2##) was used to assess the 'quality' of each reaction, and all reactions with an R<sup>2 </sup>&lt; 0.990 were excluded from further analysis (10.6% of reactions). Next, average PCR efficiencies (E) were computed for each individual primer pair across all analyzed samples. 53.4% (558 TF genes) displayed PCR efficiencies greater than 0.80, while 39.7% (415 TF genes) had efficiencies between 0.51–0.80. Only 2.6% (27 TF genes) had mean E values below 0.4; these all yielded R<sup>2 </sup>&lt; 0.99 in LinRegPCR analysis and mostly represented reactions that lacked detectable fragment amplification (C<sub>T </sub>&gt; 40) or that generated unspecific PCR products (Figure ##FIG##1##2##; see Additional file ##SUPPL##0##1##). A similar range of PCR efficiencies were determined for Arabidopsis and rice TF primers previously [##REF##15078338##19##,##REF##17559651##20##].</p> <title>Selection of reference genes</title> <p>Reference genes with stable expression/transcript levels throughout development and in the face of environmental challenge are crucial for the normalization of expression data of other genes. Potentially useful reference genes were chosen based on published data for Medicago (e.g. <italic>Msc27 </italic>[##REF##1623194##32##]) and <italic>Arabidopsis thaliana </italic>[##REF##16166256##33##]. The closest Medicago homologues of Arabidopsis genes were identified by BlastN [##UREF##3##29##]. Gene-specific primer pairs for the genes encoding elongation factor 1α (EST317575), glyceralaldehyde-3-phosphate dehydrogenase (MtC00030_GC; CT573421_3.4), β-tubulin (TC106341), Pentatricopeptide repeat protein (TC96273), actin2 (TC107326; AC137836_27.5), Ubiquitin (TC102473; AC137828-19.4), Helicase (CB892427), and the genes <italic>PDF2 </italic>(TC107161), <italic>UPL7 </italic>(TC111218), <italic>PTB </italic>(TC111751), <italic>UBC </italic>(AW686873), <italic>bHLH </italic>(CX538576), and <italic>UBC9 </italic>(TC106312) were designed using the criteria described above (Table ##TAB##1##2##). The specificity of PCR primers was tested using 18 first-strand cDNAs from six different organs of <italic>Medicago </italic>(three biological replicates each). All primer pairs produced a single PCR product of the expected size, as shown by gel electrophoresis and unique dissociation curves generated by the PCR machine after 40 cycles (Figure ##FIG##2##3##). To determine which reference genes were best suited for transcript normalisation, we used the software geNORM [##UREF##5##34##], which uses pair-wise comparison and geometric averaging across a matrix of biological samples to determine gene expression stability (M; [##REF##12184808##35##]). The genes <italic>PDF2</italic>, <italic>PPRep</italic>, <italic>Ubiquitin, and PTB h</italic>ad the lowest M (greatest transcript stability) and, therefore, were judged to be the best reference genes for this diverse set of developmental samples (Figure ##FIG##3##4##; Table ##TAB##1##2##).</p> <title>Identification of organ-enhanced TFs of Medicago</title> <p>To get an overview of TF gene expression in <italic>Medicago truncatula </italic>and to identify TFs induced in specific organs, we used the real-time RT-PCR platform described above. Transcript profiling was performed on six different organs of <italic>Medicago </italic>(leaves, stems, flowers, pods, roots, and nodules) with three independent biological replicates for each (see Additional file ##SUPPL##1##2##). The fraction of genes for which transcripts were detected within 40 cycles ranged from 77.2% in leaves to 90.8% in pods. Transcripts from nearly all putative TF genes (96.8% or 1011/1045) were detected in at least one organ. Genes were called detected if they were expressed in at least two biological replicates with a C<sub>T </sub>&lt; 40. Approximately half of all TF genes exhibited differential expression during plant development, based on significant differences (p ≤ 0.05) in transcript levels between organs. Few TF genes (1.19% or 12/1011) were expressed exclusively in vegetative organs (leaves, stems, roots or nodules), and even fewer (0.5% or 5/1011) were expressed only in reproductive organs (flowers or pods). A relatively small number of TF genes exhibited greater than ten-fold ratios in expression level in one organ compared to any other organ (Table ##TAB##2##3##). For comparison, we have included gene expression ratios derived from Affymetrix array data from the same RNA samples. While there is reasonable qualitative agreement between gene expression ratios obtained using the two methods, the lack of quantitative agreement is likely due to the limited sensitivity and low signal to noise ratio near the detection limit of Affymetrix arrays [##REF##15078338##19##]. The genes listed in Table ##TAB##2##3## may control development and/or differentiation in Medicago and are interesting targets for future research.</p>
<title>Results and Discussion</title> <title>Identification of putative transcription factors</title> <p>TF protein families are generally defined by the type(s) of DNA-binding domain they contain and putative TF genes are often identified on the basis of DNA sequences that encode known DNA-binding domains [##REF##11118137##10##,##REF##17556517##11##,##REF##15731212##25##,##REF##16769687##26##]. We utilized this approach to identify putative TFs of Medicago amongst the set of proteins predicted from genomic sequence by the International Medicago Genome Annotation Group (IMGAG). Proteins of IMGAG release 1, which contained over 40,000 predicted proteins, were screened for the presence of known or presumed DNA-binding domains (Table ##TAB##0##1##), using INTERPRO [##UREF##1##27##]. Medicago proteins containing putative DNA binding domains and other domains associated with TFs were then used as query sequences in WU-BLASTX [##UREF##2##28##] which included searches of both the non-redundant DNA database of NCBI [##UREF##3##29##] and the well-curated protein database, UniProt [##UREF##4##30##] to check annotations of related proteins in support of tentative Medicago TF assignments. This process resulted in a list of 1045 putative TF genes (see Additional file ##SUPPL##0##1##). We utilized genomic sequences rather than the large collection of partial cDNA sequences present in Expressed Sequence Tag (EST) databases for Medicago as the starting point for TF gene discovery because protein sequences derived from genomic sequence are more complete and the set of IMGAG proteins essentially contains no redundancy. Although identification of the 'complete' set of Medicago TFs from IMGAG-annotated proteins will only be possible upon completion of genome sequencing, we expect little or no redundancy in the protein set targeted by our primer set. This approach avoids wasting money on redundant primer sets and re-organization of primers when redundancy is detected, both of which would have been inevitable if we chose to use ESTs in addition to genomic sequences to identify Medicago TFs.</p> <title>PCR primer design</title> <p>To ensure maximum specificity and efficiency during PCR amplification of TF cDNA under a standard set of reaction conditions, a stringent set of criteria was used for primer design. This included predicted melting temperatures (<italic>T</italic><sub>m</sub>) of 58°C to 61°C, limited self-complementarity and poly-X, and PCR amplicon lengths of 100–150 base pairs (bp). Secondary hits were minimized by aligning primer candidates to all known Medicago sequences via WU-BLAST [##UREF##2##28##] and eliminating primer pairs with multiple potential hits.</p> <title>PCR primer testing: gene-specificity and amplification efficiency</title> <p>PCR primers were tested on Medicago cDNA free of genomic DNA contamination as follows. First, total RNA was extracted from various organs using Trizol reagent (Invitrogen GmbH, Karlsruhe, Germany), which yielded high quality RNA as judged by gel electrophoresis and by Agilent 2100 BioAnalyser using RNA 6000 Nano Chips (Agilient Technologies, Waldbronn, Germany). Typical RNA yields ranged from 0.5–1.0 μg RNA/mg fresh mass for nodules and leaves, respectively. Isolated RNA was treated with DNAse I (Ambion, product number 1907) to remove all contaminating genomic DNA, and this was always confirmed by PCR using primers to non-coding regions of the <italic>Ubiquitin </italic>gene (TC102473; AC137828-19.4). After inactivation of DNAse I, RNA was reverse transcribed using SuperScript III reverse transcriptase (Invitrogen GmbH, Karlsruhe, Germany) and oligo-dT<sub>12–18 </sub>to prime the reaction.</p> <p>Specificity of PCR primers was assessed in three ways: by melting curve analysis of PCR reaction products; by separating the products of all reactions via electrophoresis in 3% agarose gels; and by sequencing a sub-set of PCR reaction products (Figure ##FIG##0##1##). 94.5% (998/1045) of primer pairs gave unique PCR products of the expected size. Only 3.3% (34/1045) of primer pairs yielded no product and 2.2% (23/1045) gave non-specific products. (see Additional file ##SUPPL##0##1##). Sequencing was performed on 178 randomly-chosen PCR products amplified from a 1:1 mixture of leaf and root cDNA. In the vast majority of cases (92.7% or 165/178), the sequence of the PCR product was identical to that of the intended target gene. In 5.1% of cases, the amplicon sequence matched multiple related genes, including the target gene, while in only 2.2% of cases the amplicon sequence did not match the target gene sequence.</p> <p>Ideally, PCR results in an exact doubling of the amount of dsDNA after each temperature cycle. In practice, however, this is generally not the case because the reactions are less than 100% efficient. Primer sequences can affect PCR efficiency, so we determined the efficiency of each TF primer pair from amplification plots, using LinRegPCR software [##REF##12618301##31##]. First, the correlation coefficient derived from linear regression analysis of each amplification plot (e.g. see Figure ##FIG##1##2##) was used to assess the 'quality' of each reaction, and all reactions with an R<sup>2 </sup>&lt; 0.990 were excluded from further analysis (10.6% of reactions). Next, average PCR efficiencies (E) were computed for each individual primer pair across all analyzed samples. 53.4% (558 TF genes) displayed PCR efficiencies greater than 0.80, while 39.7% (415 TF genes) had efficiencies between 0.51–0.80. Only 2.6% (27 TF genes) had mean E values below 0.4; these all yielded R<sup>2 </sup>&lt; 0.99 in LinRegPCR analysis and mostly represented reactions that lacked detectable fragment amplification (C<sub>T </sub>&gt; 40) or that generated unspecific PCR products (Figure ##FIG##1##2##; see Additional file ##SUPPL##0##1##). A similar range of PCR efficiencies were determined for Arabidopsis and rice TF primers previously [##REF##15078338##19##,##REF##17559651##20##].</p> <title>Selection of reference genes</title> <p>Reference genes with stable expression/transcript levels throughout development and in the face of environmental challenge are crucial for the normalization of expression data of other genes. Potentially useful reference genes were chosen based on published data for Medicago (e.g. <italic>Msc27 </italic>[##REF##1623194##32##]) and <italic>Arabidopsis thaliana </italic>[##REF##16166256##33##]. The closest Medicago homologues of Arabidopsis genes were identified by BlastN [##UREF##3##29##]. Gene-specific primer pairs for the genes encoding elongation factor 1α (EST317575), glyceralaldehyde-3-phosphate dehydrogenase (MtC00030_GC; CT573421_3.4), β-tubulin (TC106341), Pentatricopeptide repeat protein (TC96273), actin2 (TC107326; AC137836_27.5), Ubiquitin (TC102473; AC137828-19.4), Helicase (CB892427), and the genes <italic>PDF2 </italic>(TC107161), <italic>UPL7 </italic>(TC111218), <italic>PTB </italic>(TC111751), <italic>UBC </italic>(AW686873), <italic>bHLH </italic>(CX538576), and <italic>UBC9 </italic>(TC106312) were designed using the criteria described above (Table ##TAB##1##2##). The specificity of PCR primers was tested using 18 first-strand cDNAs from six different organs of <italic>Medicago </italic>(three biological replicates each). All primer pairs produced a single PCR product of the expected size, as shown by gel electrophoresis and unique dissociation curves generated by the PCR machine after 40 cycles (Figure ##FIG##2##3##). To determine which reference genes were best suited for transcript normalisation, we used the software geNORM [##UREF##5##34##], which uses pair-wise comparison and geometric averaging across a matrix of biological samples to determine gene expression stability (M; [##REF##12184808##35##]). The genes <italic>PDF2</italic>, <italic>PPRep</italic>, <italic>Ubiquitin, and PTB h</italic>ad the lowest M (greatest transcript stability) and, therefore, were judged to be the best reference genes for this diverse set of developmental samples (Figure ##FIG##3##4##; Table ##TAB##1##2##).</p> <title>Identification of organ-enhanced TFs of Medicago</title> <p>To get an overview of TF gene expression in <italic>Medicago truncatula </italic>and to identify TFs induced in specific organs, we used the real-time RT-PCR platform described above. Transcript profiling was performed on six different organs of <italic>Medicago </italic>(leaves, stems, flowers, pods, roots, and nodules) with three independent biological replicates for each (see Additional file ##SUPPL##1##2##). The fraction of genes for which transcripts were detected within 40 cycles ranged from 77.2% in leaves to 90.8% in pods. Transcripts from nearly all putative TF genes (96.8% or 1011/1045) were detected in at least one organ. Genes were called detected if they were expressed in at least two biological replicates with a C<sub>T </sub>&lt; 40. Approximately half of all TF genes exhibited differential expression during plant development, based on significant differences (p ≤ 0.05) in transcript levels between organs. Few TF genes (1.19% or 12/1011) were expressed exclusively in vegetative organs (leaves, stems, roots or nodules), and even fewer (0.5% or 5/1011) were expressed only in reproductive organs (flowers or pods). A relatively small number of TF genes exhibited greater than ten-fold ratios in expression level in one organ compared to any other organ (Table ##TAB##2##3##). For comparison, we have included gene expression ratios derived from Affymetrix array data from the same RNA samples. While there is reasonable qualitative agreement between gene expression ratios obtained using the two methods, the lack of quantitative agreement is likely due to the limited sensitivity and low signal to noise ratio near the detection limit of Affymetrix arrays [##REF##15078338##19##]. The genes listed in Table ##TAB##2##3## may control development and/or differentiation in Medicago and are interesting targets for future research.</p>
<title>Conclusion</title> <p>We have established a flexible platform for high-throughput qRT-PCR analysis of Medicago TF gene expression that is based on gene-specific primers arrayed in 384-well plates and SYBR<sup>® </sup>Green detection of gene-specific PCR amplicons. Currently, the platform has primer pairs for 1045 TF genes and we have plans to extend this to all predicted Medicago TF genes as genome sequencing progresses. At this stage, the resource has been utilized by several groups in Europe, Australia, and the USA, and we expect it will become the 'gold-standard' for TF transcript profiling in <italic>Medicago truncatula</italic>.</p>
<p>This is an Open Access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <title>Background</title> <p><italic>Medicago truncatula </italic>is a model legume species that is currently the focus of an international genome sequencing effort. Although several different oligonucleotide and cDNA arrays have been produced for genome-wide transcript analysis of this species, intrinsic limitations in the sensitivity of hybridization-based technologies mean that transcripts of genes expressed at low-levels cannot be measured accurately with these tools. Amongst such genes are many encoding transcription factors (TFs), which are arguably the most important class of regulatory proteins. Quantitative reverse transcription-polymerase chain reaction (qRT-PCR) is the most sensitive method currently available for transcript quantification, and one that can be scaled up to analyze transcripts of thousands of genes in parallel. Thus, qRT-PCR is an ideal method to tackle the problem of TF transcript quantification in Medicago and other plants.</p> <title>Results</title> <p>We established a bioinformatics pipeline to identify putative TF genes in <italic>Medicago truncatula </italic>and to design gene-specific oligonucleotide primers for qRT-PCR analysis of TF transcripts. We validated the efficacy and gene-specificity of over 1000 TF primer pairs and utilized these to identify sets of organ-enhanced TF genes that may play important roles in organ development or differentiation in this species. This community resource will be developed further as more genome sequence becomes available, with the ultimate goal of producing validated, gene-specific primers for all Medicago TF genes.</p> <title>Conclusion</title> <p>High-throughput qRT-PCR using a 384-well plate format enables rapid, flexible, and sensitive quantification of all predicted Medicago transcription factor mRNAs. This resource has been utilized recently by several groups in Europe, Australia, and the USA, and we expect that it will become the 'gold-standard' for TF transcript profiling in <italic>Medicago truncatula</italic>.</p>
<title>Competing interests</title> <p>The authors declare that they have no competing interests.</p> <title>Authors' contributions</title> <p>KK performed the experimental work and helped draft the manuscript. W-RS, TC, and MS helped to conceive the project and provided practical advice. MW, FC, and CDT carried out bioinformatic analysis of Medicago TFs, and HW designed gene-specific primers. JCR and YX were responsible for the sequencing and analysis of the PCR products. IT carried out qRT-PCR and TG performed statistical analyses. MU designed and coordinated the project and wrote the manuscript.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>This work was supported in part by the European Commission FP6 Framework Programme Grain Legume Integrated Project (FOOD-CT-2004-506223), the Max Planck Society, the National Research Initiative (NRI) Plant Genome Program of the USDA Cooperative State Research, Education and Extension Service (CSREES), and the Samuel Roberts Noble Foundation.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p><bold>Specificity of transcription factor PCR primers</bold>. Specificity was confirmed by dissociation curves with a single peak (A) while double peaks (B) indicated off-target ampification. The derivative of fluorescence intensity is shown on the y-axis. Separation of PCR products on 3% (v/w) agarose gels following electrophoresis (C) confirmed the presence of unique amplicons of the expected size for most reactions. Few reactions yielded no products (indicated by arrow).</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p><bold>Amplification efficiency of transcription factor-specific primer pairs</bold>. Typical real-time RT-PCR amplification plots of 384 TF genes (left) and distribution of PCR efficiencies for all 1045 TF primer pairs (right).</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p><bold>Specificity, efficiency, and reproducibility of PCR primers designed to amplify reference gene transcripts</bold>. Specificity of primers was confirmed by the presence of unique amplicons of the expected size following electrophoresis on 3% (v/w) agarose gels (A) and by dissociation curves with a single peak (B to D). Typical real-time RT-PCR amplification plots of three reference gene transcripts (E to G).</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p><bold>Ranking of 8 reference genes in <italic>M. truncatula</italic></bold>. Transcript levels of all 8 genes were measured by qRT-PCR, using 18 independent cDNA preparations from six different organs with three replicate measurements of each cDNA preparation. A low value for the average expression stability M, as calculated by geNORM software, indicates more stable expression throughout the various organs.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>Classification of putative transcription factors of Medicago into families and sub-families</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center"><bold>TF Family</bold></td><td align="center"><bold>No. of </bold><break/><bold>Genes</bold></td><td align="center"><bold>Characteristic Domain </bold><break/><bold>(InterPro No.)</bold></td><td align="center"><bold>Domain</bold><break/><bold>Function</bold></td><td align="center"><bold>Domain Description</bold></td></tr></thead><tbody><tr><td align="left">MYB/HD-like</td><td align="center">76</td><td align="center">IPR001005, IPR009057</td><td align="center">D</td><td align="center">Myb, DNA-binding; homeodomain like</td></tr><tr><td align="left">MYB</td><td align="center">58</td><td align="center">IPR001005</td><td align="center">D</td><td align="center">Myb, DNA-binding</td></tr><tr><td align="left">C<sub>2</sub>H<sub>2 </sub>(Zn)</td><td align="center">64</td><td align="center">IPR007087</td><td align="center">NA</td><td align="center">Zn-finger, C2H2 type</td></tr><tr><td align="left"><bold>AP2/EREBP</bold></td><td align="center"><bold>55</bold></td><td align="center"><bold>IPR001471</bold></td><td align="center"><bold>D</bold></td><td align="center"><bold>Pathogenesis-related transcriptional factor and ethylene response factor</bold></td></tr><tr><td align="left">BHLH</td><td align="center">49</td><td align="center">IPR001092</td><td align="center">D</td><td align="center">Basic helix-loop-helix dimerisation region bHLH</td></tr><tr><td align="left">HD-like</td><td align="center">50</td><td align="center">IPR009057</td><td align="center">D</td><td align="center">Homeodomain like</td></tr><tr><td align="left">HD family</td><td/><td align="center">IPR001356</td><td align="center">D</td><td align="center">Homeobox</td></tr><tr><td align="right">HD</td><td align="center">25</td><td/><td/><td/></tr><tr><td align="right">HD-ZIP</td><td align="center">5</td><td align="center">IPR006712</td><td align="center">P</td><td align="center">HD-ZIP protein, N terminus</td></tr><tr><td align="right">HD-PHD finger</td><td align="center">2</td><td align="center">IPR001965</td><td align="center">P</td><td align="center">Zn-finger like, PHD-finger</td></tr><tr><td align="left">MADS</td><td align="center">48</td><td align="center">IPR002100</td><td align="center">D</td><td align="center">TF, MADS-box</td></tr><tr><td align="left">BZIP</td><td align="center">41</td><td align="center">IPR004827</td><td align="center">D</td><td align="center">Basic Leu zipper (bZIP) TF</td></tr><tr><td align="left">PHD</td><td align="center">34</td><td align="center">IPR001965</td><td align="center">P</td><td align="center">Zn-finger like, PHD-finger</td></tr><tr><td align="left"><bold>WRKY family</bold></td><td/><td align="center"><bold>IPR003657</bold></td><td align="center"><bold>D</bold></td><td align="center"><bold>DNA-binding WRKY</bold></td></tr><tr><td align="right"><bold>WRKY</bold></td><td align="center"><bold>29</bold></td><td/><td/><td/></tr><tr><td align="right"><bold>LLR WRKY</bold></td><td align="center"><bold>1</bold></td><td align="center"><bold>IPR001611</bold></td><td/><td align="center"><bold>Leu-rich repeat</bold></td></tr><tr><td align="left"><bold>ABI3/VP1</bold></td><td align="center"><bold>29</bold></td><td align="center"><bold>IPR003340</bold></td><td align="center"><bold>D</bold></td><td align="center"><bold>TF B3</bold></td></tr><tr><td align="left"><bold>NAC</bold></td><td align="center"><bold>29</bold></td><td align="center"><bold>IPR003441</bold></td><td align="center"><bold>D</bold></td><td align="center"><bold>No apical meristem (NAM) protein</bold></td></tr><tr><td align="left">C<sub>3</sub>H-type1(Zn)</td><td align="center">27</td><td align="center">IPR000571</td><td align="center">D</td><td align="center">Zn-finger, C-x8-C-x5-C-x3-H type</td></tr><tr><td align="left"><bold>ARF</bold></td><td align="center"><bold>23</bold></td><td align="center"><bold>IPR003340, IPR010525, IPR011525</bold></td><td align="center"><bold>D</bold></td><td/></tr><tr><td align="left">JUMONJI</td><td align="center">20</td><td align="center">IPR003347</td><td align="center">D</td><td align="center">TF jumonji, jmjC</td></tr><tr><td align="left"><bold>GRAS</bold></td><td align="center"><bold>19</bold></td><td align="center"><bold>IPR005202</bold></td><td align="center"><bold>P</bold></td><td align="center"><bold>GRAS TF</bold></td></tr><tr><td align="left">HMG</td><td align="center">15</td><td align="center">IPR000637</td><td align="center">D</td><td align="center">HMG-I and HMG-Y, DNA binding</td></tr><tr><td align="left">AS2</td><td align="center">14</td><td align="center">IPR004883</td><td align="center">P</td><td align="center">Lateral organ boundaries</td></tr><tr><td align="left">C<sub>2</sub>C<sub>2 </sub>(Zn)</td><td/><td/><td/><td/></tr><tr><td align="right"><bold>Dof</bold></td><td align="center"><bold>13</bold></td><td align="center"><bold>IPR003851</bold></td><td align="center"><bold>D</bold></td><td align="center"><bold>Zn-finger, Dof type</bold></td></tr><tr><td align="right">GATA</td><td align="center">7</td><td align="center">IPR000679</td><td align="center">D</td><td align="center">Zn-finger, GATA type</td></tr><tr><td align="right"><bold>CO-like</bold></td><td align="center"><bold>6</bold></td><td align="center"><bold>IPR000315</bold></td><td align="center"><bold>D</bold></td><td align="center"><bold>Zn-finger, B-box</bold></td></tr><tr><td align="right"><bold>YABBY</bold></td><td align="center"><bold>5</bold></td><td align="center"><bold>IPR006780</bold></td><td align="center"><bold>D</bold></td><td align="center"><bold>YABBY protein</bold></td></tr><tr><td align="left">CCAAT-HAP3 type</td><td align="center">12</td><td align="center">IPR003958</td><td align="center">D</td><td align="center">TF CBF/NF-Y/archaeal histone</td></tr><tr><td align="left"><bold>GRF</bold></td><td align="center"><bold>8</bold></td><td align="center"><bold>IPR010666</bold></td><td align="center"><bold>D</bold></td><td align="center"><bold>Zn-finger, GRF type</bold></td></tr><tr><td align="left"><bold>SBP</bold></td><td align="center"><bold>8</bold></td><td align="center"><bold>IPR004333</bold></td><td align="center"><bold>D</bold></td><td align="center"><bold>SBP</bold></td></tr><tr><td align="left"><bold>EIL</bold></td><td align="center"><bold>7</bold></td><td align="center"><bold>IPR006957</bold></td><td align="center"><bold>D</bold></td><td align="center"><bold>Ethylene insensitive 3</bold></td></tr><tr><td align="left">LIM</td><td align="center">7</td><td align="center">IPR001781</td><td align="center">P</td><td align="center">Zn-binding protein, LIM</td></tr><tr><td align="left">SNF2</td><td align="center">6</td><td align="center">IPR000330</td><td align="center">D</td><td align="center">SNF2 family N-terminal domain</td></tr><tr><td align="left">E2F/DP</td><td align="center">5</td><td align="center">IPR003316</td><td align="center">D</td><td align="center">TF E2F/dimerisation partner (TDP)</td></tr><tr><td align="left"><bold>TCP</bold></td><td align="center"><bold>5</bold></td><td align="center"><bold>IPR005333</bold></td><td align="center"><bold>D</bold></td><td align="center"><bold>TCP TF</bold></td></tr><tr><td align="left"><bold>FHA</bold></td><td align="center"><bold>5</bold></td><td align="center"><bold>IPR000253</bold></td><td align="center"><bold>D</bold></td><td align="center"><bold>Forkhead-associated</bold></td></tr><tr><td align="left">ARID</td><td align="center">4</td><td align="center">IPR001606</td><td align="center">D</td><td align="center">AT-rich interaction region</td></tr><tr><td align="left">HSF</td><td align="center">4</td><td align="center">IPR000232</td><td align="center">D</td><td align="center">Heat shock factor (HSF)-type, DNA binding</td></tr><tr><td align="left"><bold>AUX/IAA</bold></td><td align="center"><bold>3</bold></td><td align="center"><bold>IPR003311</bold></td><td align="center"><bold>D</bold></td><td align="center"><bold>AUX/IAA protein</bold></td></tr><tr><td align="left"><bold>SRS</bold></td><td align="center"><bold>3</bold></td><td align="center"><bold>IPR006510</bold></td><td align="center"><bold>D</bold></td><td align="center"><bold>Zn-finger, LRP1 type</bold></td></tr><tr><td align="left">TUB</td><td align="center">3</td><td align="center">IPR000007</td><td align="center">D</td><td align="center">Tubby</td></tr><tr><td align="left">ZIM</td><td align="center">3</td><td align="center">IPR010399</td><td align="center">D</td><td align="center">ZIM</td></tr><tr><td align="left">DDT</td><td align="center">3</td><td align="center">IPR004022</td><td align="center">D</td><td align="center">DDT</td></tr><tr><td align="left"><bold>ZF-HD</bold></td><td align="center"><bold>2</bold></td><td align="center"><bold>IPR006455</bold></td><td align="center"><bold>D</bold></td><td align="center"><bold>Homeobox domain, ZF-HD class</bold></td></tr><tr><td align="left"><bold>MBF1</bold></td><td align="center"><bold>2</bold></td><td align="center"><bold>IPR001387</bold></td><td align="center"><bold>D</bold></td><td align="center"><bold>Helix-turn-helix type 3</bold></td></tr><tr><td align="left"><bold>S1Fa-like</bold></td><td align="center"><bold>2</bold></td><td align="center"><bold>IPR006779</bold></td><td align="center"><bold>D</bold></td><td align="center"><bold>DNA binding protein S1FA</bold></td></tr><tr><td align="left"><bold>CAMTA</bold></td><td align="center"><bold>2</bold></td><td align="center"><bold>IPR005559</bold></td><td align="center"><bold>D</bold></td><td align="center"><bold>CG-1</bold></td></tr><tr><td align="left"><bold>LFY</bold></td><td align="center"><bold>1</bold></td><td align="center"><bold>IPR002910</bold></td><td align="center"><bold>D</bold></td><td align="center"><bold>Floricaula/leafy protein</bold></td></tr><tr><td align="left"><bold>Nin-like</bold></td><td align="center"><bold>1</bold></td><td align="center"><bold>IPR003035</bold></td><td align="center"><bold>D</bold></td><td align="center"><bold>Plant regulator RWP-RK</bold></td></tr><tr><td align="left"><bold>TAZ</bold></td><td align="center"><bold>1</bold></td><td align="center"><bold>IPR000197</bold></td><td align="center"><bold>P</bold></td><td align="center"><bold>Zn-finger, TAZ-type</bold></td></tr><tr><td colspan="5"><hr/></td></tr><tr><td align="left" colspan="5"><bold>Potentially novel plant TFs and transcriptional regulators</bold></td></tr><tr><td colspan="5"><hr/></td></tr><tr><td align="left">CCHC (Zn)</td><td align="center">112</td><td align="center">IPR001878</td><td align="center">NA</td><td align="center">Zn-finger, CCHC-type</td></tr><tr><td align="left">RR</td><td align="center">16</td><td align="center">IPR001789, IPR011006</td><td align="center">RD</td><td align="center">Response regulator receiver</td></tr><tr><td align="left">DHHC (Zn)</td><td align="center">14</td><td align="center">IPR001594</td><td align="center">D or P</td><td align="center">Zn-finger, DHHC-type</td></tr><tr><td align="left">HTH</td><td/><td/><td/><td/></tr><tr><td align="right">FIS</td><td align="center">11</td><td align="center">IPR002197</td><td align="center">D</td><td align="center">Helix-turn-helix, Fis-type</td></tr><tr><td align="right">AraC</td><td align="center">2</td><td align="center">IPR000005</td><td align="center">D</td><td align="center">Helix-turn-helix, AraC type</td></tr><tr><td align="left">BTB/POZ</td><td align="center">7</td><td align="center">IPR000210</td><td align="center">P</td><td align="center">BTB</td></tr><tr><td align="left">TTF-type (Zn)</td><td align="center">6</td><td align="center">IPR006580</td><td align="center">D</td><td align="center">Zn-finger, TTF-type</td></tr><tr><td align="left">BD</td><td align="center">6</td><td align="center">IPR001487</td><td align="center">P</td><td align="center">Bromodomain</td></tr><tr><td align="left">Lambda-DB</td><td align="center">3</td><td align="center">IPR010982</td><td align="center">D</td><td align="center">Lambda_DNA_bd</td></tr><tr><td align="left">TrpR</td><td align="center">3</td><td align="center">IPR010921</td><td align="center">D</td><td align="center">Trp repressor/replication initiator</td></tr><tr><td align="left">TPR</td><td align="center">3</td><td align="center">IPR001440</td><td align="center">P</td><td align="center">Tetratricopeptide TPR_1</td></tr><tr><td align="left">KRAB-box</td><td align="center">2</td><td align="center">IPR001909</td><td align="center">P</td><td align="center">KRAB box</td></tr><tr><td align="left">NRs</td><td align="center">2</td><td align="center">IPR008946</td><td align="center">LBD</td><td align="center">Steroid nuclear receptor, ligand binding</td></tr><tr><td align="left">R3H</td><td align="center">2</td><td align="center">IPR001374</td><td align="center">NA</td><td align="center">Single-stranded nucleic acid binding R3H</td></tr><tr><td align="left">YEATS</td><td align="center">2</td><td align="center">IPR005033</td><td align="center">TA</td><td align="center">YEATS</td></tr><tr><td align="left">U1-type (Zn)</td><td align="center">2</td><td align="center">IPR003604</td><td align="center">NA</td><td align="center">Zn-finger, U1-type</td></tr><tr><td align="left">A20-like</td><td align="center">2</td><td align="center">IPR002653</td><td align="center">P</td><td align="center">Zn-finger, A20-type</td></tr><tr><td align="left">Euk_TF</td><td align="center">1</td><td align="center">IPR008917</td><td align="center">D</td><td align="center">Euk_TF, DNA binding</td></tr><tr><td align="left">NGN</td><td align="center">1</td><td align="center">IPR006645</td><td align="center">D</td><td align="center">NGN</td></tr><tr><td align="left">p53-like</td><td align="center">1</td><td align="center">IPR008967</td><td align="center">D</td><td align="center">p53-like TF, DNA binding</td></tr><tr><td align="left">SSB protein</td><td align="center">1</td><td align="center">IPR011344</td><td align="center">D</td><td align="center">Single-strand binding protein</td></tr><tr><td align="left">ssDB TR</td><td align="center">1</td><td align="center">IPR009044</td><td align="center">D</td><td align="center">Single-strand DNA binding transcriptional regulator</td></tr><tr><td align="left">TCoAp15</td><td align="center">1</td><td align="center">IPR003173</td><td align="center">D</td><td align="center">Transcriptional coactivator p15</td></tr><tr><td align="left">BED-type (Zn)</td><td align="center">1</td><td align="center">IPR003656</td><td align="center">D</td><td align="center">Zn-finger, BED-type predicted</td></tr><tr><td align="left">TCoA</td><td align="center">1</td><td align="center">IPR009255</td><td align="center">TA</td><td align="center">Transcriptional coactivation</td></tr><tr><td align="left">Tc/PD</td><td align="center">1</td><td align="center">IPR001533</td><td align="center">TA</td><td align="center">Transcriptional coactivator</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Medicago reference genes and primers for qRT-PCR</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="center"><bold>Gene Name</bold></td><td align="center"><bold>TC</bold></td><td align="center"><bold>Accession</bold><break/><bold>Number</bold></td><td align="center"><bold>Forward/Reverse Primer (5'-3')</bold></td><td align="center"><bold>PCR </bold><break/><bold>Product</bold><break/><bold> Size (bp)</bold></td><td align="center"><bold>PCR </bold><break/><bold>efficiency</bold><break/><bold> (E)</bold></td><td align="center"><bold>R</bold><sup>2</sup></td></tr></thead><tbody><tr><td align="center">β Tubulin</td><td align="center">TC106341</td><td align="center">N</td><td align="center">TTTGCTCCTCTTACATCCCGTG / GCAGCACACATCATGTTTTTGG</td><td align="center">100</td><td align="center">1.08</td><td align="center">1.00</td></tr><tr><td align="center">PPRrep</td><td align="center">TC96273</td><td align="center">N</td><td align="center">GGAAAACTGGAGGATGCACGTA / ACAAGCCCTCGACACAAAACC</td><td align="center">100</td><td align="center">0.93</td><td align="center">1.00</td></tr><tr><td align="center">PDF 2</td><td align="center">TC107161</td><td align="center">N</td><td align="center">GTGTTTTGCTTCCGCCGTT / CCAAATCTTGCTCCCTCATCTG</td><td align="center">100</td><td align="center">0.99</td><td align="center">1.00</td></tr><tr><td align="center">bHLH</td><td align="center">CX538576</td><td align="center">N</td><td align="center">TAGCGAGTACCATGATGCCAGA / GCGCCTCTTTTGTTTTCAGC</td><td align="center">100</td><td align="center">0.89</td><td align="center">1.00</td></tr><tr><td align="center">UBC</td><td align="center">AW686873</td><td align="center">N</td><td align="center">CTGACAGCCCACTGAATTGTGA / TTTTGGCATTGCTGCAAGC</td><td align="center">100</td><td align="center">0.96</td><td align="center">1.00</td></tr><tr><td align="center">PTB</td><td align="center">TC111751</td><td align="center">N</td><td align="center">CGCCTTGTCAGCATTGATGTC / TGAACCAGTGCCTGGAATCCT</td><td align="center">100</td><td align="center">0.85</td><td align="center">1.00</td></tr><tr><td align="center">Ubiquitin</td><td align="center">TC102473</td><td align="center">AC137828_19.4</td><td align="center">GCAGATAGACACGCTGGGA / AACTCTTGGGCAGGCAATAA</td><td align="center">100</td><td align="center">0.95</td><td align="center">1.00</td></tr><tr><td align="center">UBC9</td><td align="center">TC106312</td><td align="center">AC137602_2.4</td><td align="center">GGTTGATTGCTCTTCTCTCCCC / AAGTGATTGCTCGTCCAACCC</td><td align="center">100</td><td align="center">1.13</td><td align="center">0.99</td></tr><tr><td align="center">Helicase</td><td align="center">CB892427</td><td align="center">N</td><td align="center">GTACGAGGTCGGTGCTCTTGAA / GCAACCGAAAATTGCACCATAC</td><td align="center">100</td><td align="center">0.91</td><td align="center">1.00</td></tr><tr><td align="center">ELF1α</td><td align="center">EST317575</td><td align="center">N</td><td align="center">GACAAGCGTGTGATCGAGAGATT / TTTCACGCTCAGCCTTAAGCT</td><td align="center">100</td><td align="center">0.68</td><td align="center">0.98</td></tr><tr><td align="center">UPL7</td><td align="center">TC111218</td><td align="center">N</td><td align="center">CCAGTTGTTCTCGTGGTCCATT / CCTCCAATTGTCGCCCAAA</td><td align="center">100</td><td align="center">0.93</td><td align="center">1.00</td></tr><tr><td align="center">GAPDH</td><td align="center">MtC00030_GC</td><td align="center">CT573421_3.4</td><td align="center">TGCCTACCGTCGATGTTTCAGT / TTGCCCTCTGATTCCTCCTTG</td><td align="center">100</td><td align="center">1.04</td><td align="center">0.99</td></tr><tr><td align="center">Actin2</td><td align="center">TC107326</td><td align="center">AC137836_27.5</td><td align="center">TCAATGTGCCTGCCATGTATGT / ACTCACACCGTCACCAGAATCC</td><td align="center">100</td><td align="center">1.12</td><td align="center">0.99</td></tr><tr><td align="center">MSC27</td><td align="center">X63872 (<italic>M. sativa</italic>)</td><td align="center">N</td><td align="center">GTTGAAGTAGACATTGGTGCTAACG / AGCTGAGTCATCAACACCCTCAT</td><td align="center">100</td><td align="center">0.76</td><td align="center">0.99</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>Organ-enhanced TF genes</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left" colspan="12"><bold>A </bold>Root-enhanced TF genes identified by real-time RT-PCR</td></tr></thead><tbody><tr><td/><td align="center" colspan="5"><bold>Real-time RT-PCR Expression Ratio</bold></td><td align="center" colspan="5"><bold>Affymetrix Expression Ratio</bold></td><td/></tr><tr><td colspan="12"><hr/></td></tr><tr><td align="center"><bold>Accession Number</bold></td><td align="center"><bold>R/L</bold></td><td align="center"><bold>R/S</bold></td><td align="center"><bold>R/F</bold></td><td align="center"><bold>R/P</bold></td><td align="center"><bold>R/N</bold></td><td align="center"><bold>R/L</bold></td><td align="center"><bold>R/S</bold></td><td align="center"><bold>R/F</bold></td><td align="center"><bold>R/P</bold></td><td align="center"><bold>R/N</bold></td><td align="center"><bold>Affymetrix Chip ID</bold></td></tr><tr><td colspan="12"><hr/></td></tr><tr><td align="center">AC140721_12.1</td><td align="center">121.5</td><td align="center">320.1</td><td align="center">1623.4</td><td align="center">2676.5</td><td align="center"><bold>64.4</bold></td><td align="center">161.2<sup>b</sup></td><td align="center">144.5<sup>b</sup></td><td align="center">219.2<sup>b</sup></td><td align="center">158.1<sup>b</sup></td><td align="center">99.9</td><td align="center">Mtr.50075.1.S1_s_at</td></tr><tr><td align="center">AC135101_25.1</td><td align="center">312.2</td><td align="center">438.0</td><td align="center">291.6</td><td align="center">787.0</td><td align="center"><bold>61.3</bold></td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td></tr><tr><td align="center">AC140721_14.1</td><td align="center">1204.9</td><td align="center">751.2</td><td align="center">956.6</td><td align="center">6296.6</td><td align="center"><bold>32.6</bold></td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td></tr><tr><td align="center">AC140031_3.1</td><td align="center">275.4</td><td align="center">1157.3</td><td align="center">923.4</td><td align="center">258916.9<sup>a</sup></td><td align="center"><bold>31.1</bold></td><td align="center">0.4</td><td align="center">0.7</td><td align="center">0.8</td><td align="center">0.9</td><td align="center">2.2</td><td align="center">Mtr.47227.1.S1_s_at</td></tr><tr><td align="center">AC140721_13.1</td><td align="center">39.9<sup>a</sup></td><td align="center">179.8</td><td align="center">614.8</td><td align="center">935.1<sup>a</sup></td><td align="center"><bold>26.5</bold></td><td align="center">35.9</td><td align="center">48.1<sup>b</sup></td><td align="center">41.2</td><td align="center">44.2<sup>b</sup></td><td align="center">39.5</td><td align="center">Mtr.50074.1.S1_at</td></tr><tr><td align="center">AC140031_7.1</td><td align="center">360.2<sup>a</sup></td><td align="center">830.6</td><td align="center">1542<sup>a</sup></td><td align="center">6868.8</td><td align="center"><bold>19.7</bold></td><td align="center">0.7<sup>b</sup></td><td align="center">1<sup>b</sup></td><td align="center">1.1<sup>b</sup></td><td align="center">1<sup>b</sup></td><td align="center">0.9<sup>b</sup></td><td align="center">Mtr.47229.1.S1_at</td></tr><tr><td align="center">AC146574_6.1</td><td align="center">1095.7</td><td align="center"><bold>11.9</bold><sup>a</sup></td><td align="center">48.8</td><td align="center">51.8</td><td align="center">88.8</td><td align="center">1.2<sup>b</sup></td><td align="center">1.1<sup>b</sup></td><td align="center">1.3<sup>b</sup></td><td align="center">1.1<sup>b</sup></td><td align="center">1.2<sup>b</sup></td><td align="center">Mtr.40781.1.S1_s_at</td></tr><tr><td align="center">AC125478_13.7</td><td align="center">1117.4</td><td align="center">7008.0</td><td align="center">3264.9<sup>a</sup></td><td align="center">7256.6</td><td align="center"><bold>11.6</bold></td><td align="center">211.5<sup>b</sup></td><td align="center">248.8<sup>b</sup></td><td align="center">296.4<sup>b</sup></td><td align="center">228.7<sup>b</sup></td><td align="center">45.3</td><td align="center">Mtr.15416.1.S1_at</td></tr><tr><td align="center">AC125478_7.2</td><td align="center">1052.3</td><td align="center">503.3</td><td align="center">938.4<sup>a</sup></td><td align="center">2754.5</td><td align="center"><bold>11.2</bold></td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td></tr><tr><td align="center">AC122726_21.111</td><td align="center">50.0</td><td align="center">48.0</td><td align="center">20943.9</td><td align="center">2658.6</td><td align="center"><bold>10.4</bold></td><td align="center">74.3<sup>b</sup></td><td align="center">76.4<sup>b</sup></td><td align="center">72.8<sup>b</sup></td><td align="center">82.8<sup>b</sup></td><td align="center">26.6</td><td align="center">Mtr.15568.1.S1_s_at</td></tr><tr><td colspan="12"><hr/></td></tr><tr><td align="left" colspan="12"><bold>B </bold>Nodule-enhanced TF genes identified by real-time RT-PCR</td></tr><tr><td colspan="12"><hr/></td></tr><tr><td/><td align="center" colspan="5"><bold>Real-time RT-PCR Expression Ratio</bold></td><td align="center" colspan="5"><bold>Affymetrix Expression Ratio</bold></td><td/></tr><tr><td colspan="12"><hr/></td></tr><tr><td align="center"><bold>Accession Number</bold></td><td align="center"><bold>N/L</bold></td><td align="center"><bold>N/S</bold></td><td align="center"><bold>N/F</bold></td><td align="center"><bold>N/P</bold></td><td align="center"><bold>N/R</bold></td><td align="center"><bold>N/L</bold></td><td align="center"><bold>N/S</bold></td><td align="center"><bold>N/F</bold></td><td align="center"><bold>N/P</bold></td><td align="center"><bold>N/R</bold></td><td align="center"><bold>Affymetrix Chip ID</bold></td></tr><tr><td colspan="12"><hr/></td></tr><tr><td align="center">AC148816_3.2</td><td align="center">9966.7</td><td align="center">1826.2</td><td align="center">6262.9</td><td align="center">21709.1</td><td align="center"><bold>662.7</bold></td><td align="center">434.6<sup>b</sup></td><td align="center">516.8<sup>b</sup></td><td align="center">455<sup>b</sup></td><td align="center">521.7<sup>b</sup></td><td align="center">551.4<sup>b</sup></td><td align="center">Mtr.14503.1.S1_at</td></tr><tr><td align="center">AC147774_3.2</td><td align="center"><bold>55.2</bold><sup>a</sup></td><td align="center">466.7</td><td align="center">726.3</td><td align="center">881.2</td><td align="center">498.5</td><td align="center">17.7<sup>b</sup></td><td align="center">16<sup>b</sup></td><td align="center">17.8<sup>b</sup></td><td align="center">19.1<sup>b</sup></td><td align="center">16.1<sup>b</sup></td><td align="center">Mtr.19554.1.S1_at</td></tr><tr><td align="center">AC138056_33.241</td><td align="center">43.9</td><td align="center">93.2</td><td align="center">18.3</td><td align="center">17.3</td><td align="center"><bold>13.1</bold></td><td align="center">0.9<sup>b</sup></td><td align="center">0.8<sup>b</sup></td><td align="center">0.8<sup>b</sup></td><td align="center">1.1<sup>b</sup></td><td align="center">1.1<sup>b</sup></td><td align="center">Mtr.17993.1.S1_at</td></tr><tr><td align="center">AC124214_39.2</td><td align="center">121.8</td><td align="center">92.1</td><td align="center">44.8</td><td align="center">64.9</td><td align="center"><bold>10.5</bold></td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td></tr><tr><td colspan="12"><hr/></td></tr><tr><td align="left" colspan="12"><bold>C </bold>Pod-enhanced TF genes identified by real-time RT-PCR</td></tr><tr><td colspan="12"><hr/></td></tr><tr><td/><td align="center" colspan="5"><bold>Real-time RT-PCR Expression Ratio</bold></td><td align="center" colspan="5"><bold>Affymetrix Expression Ratio</bold></td><td/></tr><tr><td colspan="12"><hr/></td></tr><tr><td align="center"><bold>Accession Number</bold></td><td align="center"><bold>P/L</bold></td><td align="center"><bold>P/S</bold></td><td align="center"><bold>P/F</bold></td><td align="center"><bold>P/N</bold></td><td align="center"><bold>P/R</bold></td><td align="center"><bold>P/L</bold></td><td align="center"><bold>P/S</bold></td><td align="center"><bold>P/F</bold></td><td align="center"><bold>P/N</bold></td><td align="center"><bold>P/R</bold></td><td align="center"><bold>Affymetrix Chip ID</bold></td></tr><tr><td colspan="12"><hr/></td></tr><tr><td align="center">AC143340_4.7</td><td align="center">154.1</td><td align="center"><bold>44.8</bold></td><td align="center">173.3</td><td align="center">180.3<sup>a</sup></td><td align="center">407.3</td><td align="center">1.6<sup>b</sup></td><td align="center">1.8<sup>b</sup></td><td align="center">1.7<sup>b</sup></td><td align="center">1.7<sup>b</sup></td><td align="center">1.7<sup>b</sup></td><td align="center">Mtr.17931.1.S1_at</td></tr><tr><td colspan="12"><hr/></td></tr><tr><td align="left" colspan="12"><bold>D </bold>Flower-enhanced TF genes identified by real-time RT-PCR</td></tr><tr><td colspan="12"><hr/></td></tr><tr><td/><td align="center" colspan="5"><bold>Real-time RT-PCR Expression Ratio</bold></td><td align="center" colspan="5"><bold>Affymetrix Expression Ratio</bold></td><td/></tr><tr><td colspan="12"><hr/></td></tr><tr><td align="center"><bold>Accession Number</bold></td><td align="center"><bold>F/L</bold></td><td align="center"><bold>F/S</bold></td><td align="center"><bold>F/P</bold></td><td align="center"><bold>F/N</bold></td><td align="center"><bold>F/R</bold></td><td align="center"><bold>F/L</bold></td><td align="center"><bold>F/S</bold></td><td align="center"><bold>F/P</bold></td><td align="center"><bold>F/N</bold></td><td align="center"><bold>F/R</bold></td><td align="center"><bold>Affymetrix Chip ID</bold></td></tr><tr><td colspan="12"><hr/></td></tr><tr><td align="center">AC129092_13.1</td><td align="center"><bold>40.3</bold></td><td align="center">383.3</td><td align="center">65.1</td><td align="center">1683.9</td><td align="center">318.1</td><td align="center">24.2</td><td align="center">47.1<sup>b</sup></td><td align="center">38.2</td><td align="center">45.4<sup>b</sup></td><td align="center">55.1<sup>b</sup></td><td align="center">Mtr.16432.1.S1_at</td></tr><tr><td align="center">AC148485_10.1</td><td align="center">476.2</td><td align="center">203.7</td><td align="center"><bold>27.2</bold></td><td align="center">174.6</td><td align="center">44.5</td><td align="center">36.5<sup>b</sup></td><td align="center">32.9<sup>b</sup></td><td align="center">22.0</td><td align="center">34.7<sup>b</sup></td><td align="center">29.7<sup>b</sup></td><td align="center">Mtr.20392.1.S1_at</td></tr><tr><td align="center">AC140915_20.1</td><td align="center">169.5</td><td align="center">101.1</td><td align="center"><bold>18.6</bold></td><td align="center">481.0</td><td align="center">178.7</td><td align="center">1.4<sup>b</sup></td><td align="center">1.4<sup>b</sup></td><td align="center">1.3<sup>b</sup></td><td align="center">1.3<sup>b</sup></td><td align="center">1.6<sup>b</sup></td><td align="center">Mtr.51688.1.S1_at</td></tr><tr><td align="center">AC141107_50.2</td><td align="center">38.2</td><td align="center">33.0</td><td align="center"><bold>17.9</bold></td><td align="center">114.0</td><td align="center">97.4</td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td></tr><tr><td align="center">AC144731_15.21</td><td align="center">64.4</td><td align="center">44.3</td><td align="center">33.6</td><td align="center">61.5</td><td align="center"><bold>16.3</bold></td><td align="center">0.6<sup>b</sup></td><td align="center">0.6<sup>b</sup></td><td align="center">0.8<sup>b</sup></td><td align="center">0.8<sup>b</sup></td><td align="center">0.8<sup>b</sup></td><td align="center">Mtr.19093.1.S1_at</td></tr><tr><td align="center">AC150978_12.1</td><td align="center">235.3</td><td align="center">198.5</td><td align="center">48.4</td><td align="center"><bold>14.0</bold></td><td align="center">16.9</td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td></tr><tr><td align="center">AC141107.5.61</td><td align="center">261.1</td><td align="center">127.8</td><td align="center"><bold>13.3</bold></td><td align="center">96.4</td><td align="center">71.5</td><td align="center">0.9<sup>b</sup></td><td align="center">0.9<sup>b</sup></td><td align="center">0.8<sup>b</sup></td><td align="center">0.1<sup>b</sup></td><td align="center">1.1<sup>b</sup></td><td align="center">Mtr.51651.1.S1_at</td></tr><tr><td align="center">AC157472_19.1</td><td align="center">60.2<sup>a</sup></td><td align="center">92.1</td><td align="center"><bold>13.2</bold></td><td align="center">61.5<sup>a</sup></td><td align="center">113.6</td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td></tr><tr><td align="center">AC148527_19.141</td><td align="center">188.4</td><td align="center">42.1</td><td align="center">17.3</td><td align="center">26.5</td><td align="center"><bold>13.0</bold></td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td></tr><tr><td align="center">AC144726_6.1</td><td align="center">254.7</td><td align="center">55.2</td><td align="center"><bold>10.9</bold></td><td align="center">112.7<sup>a</sup></td><td align="center">754.7</td><td align="center">59.2</td><td align="center">73.7<sup>b</sup></td><td align="center">7.6</td><td align="center">89.2</td><td align="center">97.9<sup>b</sup></td><td align="center">Mtr.19024.1.S1_at</td></tr><tr><td align="center">AC157488_16.1</td><td align="center">7053.2<sup>a</sup></td><td align="center">141.1</td><td align="center">18.7</td><td align="center">81.8</td><td align="center"><bold>10.7</bold></td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td><td align="center">n</td></tr><tr><td align="center">AC123899_15.181</td><td align="center"><bold>10.2</bold></td><td align="center">92.2</td><td align="center">19.1</td><td align="center">62.8</td><td align="center">24.9</td><td align="center">0.9<sup>b</sup></td><td align="center">0.5<sup>b</sup></td><td align="center">0.8<sup>b</sup></td><td align="center">0.1<sup>b</sup></td><td align="center">0.8<sup>b</sup></td><td align="center">Mtr.52015.1.S1_at</td></tr></tbody></table></table-wrap>
<inline-formula><mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M1" name="1746-4811-4-18-i1" overflow="scroll"><mml:semantics><mml:mrow><mml:msup><mml:mrow><mml:mo stretchy="false">(</mml:mo><mml:mn>1</mml:mn><mml:mo>+</mml:mo><mml:mtext>E</mml:mtext><mml:mo stretchy="false">)</mml:mo></mml:mrow><mml:mrow><mml:mi>Δ</mml:mi><mml:mi>Δ</mml:mi><mml:msub><mml:mtext>C</mml:mtext><mml:mtext>T</mml:mtext></mml:msub></mml:mrow></mml:msup></mml:mrow></mml:semantics></mml:math></inline-formula>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional file 1</title><p>Complete list of TF genes, primer sequences and corresponding PCR efficiencies. The TF primer platform was established based on the data for the gene models according to SA, specific amplification; NA, no amplification; NS, non-specific amplification. RT-PCR products of primer sequences indicated in bold were sequenced.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional file 2</title><p>Complete list of TF genes, gene families and experimental data. Shown are the Medicago gene Accession Numbers (TIGR) in ascending order, the TC number if available, as well as the transcription factor family and the subfamily names. The next columns show experimental results for the real-time RT-PCR reactions performed on six different organs of Medicago in three indipendent biological replicates, as explained in the manuscript. CT = not normalized CT value, ΔCT = CT value normalized against the geometric mean of 4 house keeping genes; ΔΔCT = power(PCReff;-ΔCT); log2 ΔΔCT = logarithmus of ΔCT.</p></caption></supplementary-material>
<table-wrap-foot><p>IMGAG (International Medicago Gene Annotation Group)-proteins were classified as putative TFs if they contained characteristic DNA-binding or other characteristic TF domains and if annotations of matching proteins obtained by BLAST searches were consistent with such a classification. TF families previously identified in plants are presented in the first part of the table while potentially novel plant TF families, which were identified by the presence of domains associated with TFs and other transcriptional regulators outside the plant kingdom, are presented in the latter part of the table. D = DNA binding domain; P = protein-protein interaction domain; NA = nucleic acid (DNA and RNA) binding domain; RD = receiver domain; LBD = ligand binding; transcriptional co-activator. Plant-specific TF families and sub-families are indicated in bold (according to [##UREF##0##12##])</p></table-wrap-foot> <table-wrap-foot><p>Mean PCR efficiency (E) was determined from three biological replicates of each of six organs, using LinRegPCR [##REF##12618301##31##], which also yielded mean R<sup>2</sup>. N, no corresponding GenBank accession number.</p></table-wrap-foot> <table-wrap-foot><p>A TF gene was considered organ-enhanced if transcript levels for that gene were more than 10-fold higher in one organ than in any other organ. Transcript ratios were calculated using the mean of three biological replicates for each organ. Data from qRT-PCR are compared with data for the same RNA samples obtained from Affymetrix Gene chips [##UREF##7##38##]. Affymetrix data were normalized using the Robust Multiarray Average (RMA) method, as described by [##UREF##7##38##], prior to calculation of ratios. Data in bold represent the lowest transcript ratio of the corresponding gene across all organs. n = not present on Affymetrix chip; a = Ct &gt; 40 in two or three biological replicates of denominator organ; b = transcript called 'absent' by Affymetrix software in two or three biological replicates of denominator organ. L = leaf; S = stem; P = pod; F = flower; R = root; N = nodule.</p></table-wrap-foot>
<graphic xlink:href="1746-4811-4-18-1"/> <graphic xlink:href="1746-4811-4-18-2"/> <graphic xlink:href="1746-4811-4-18-3"/> <graphic xlink:href="1746-4811-4-18-4"/>
<media xlink:href="1746-4811-4-18-S1.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media> <media xlink:href="1746-4811-4-18-S2.xls" mimetype="application" mime-subtype="vnd.ms-excel"><caption><p>Click here for file</p></caption></media>
[{"surname": ["Riechmann", "Somerville CR, Meyerowitz EM"], "given-names": ["JL"], "article-title": ["Transcriptional Regulation: a Genomic Overview"], "source": ["The Arabidopsis Book"], "year": ["2002"], "publisher-name": [" The American Society of Plant Biologists, Rockville"]}, {"article-title": ["InterPro: Database of protein families, domains and sites"]}, {"article-title": ["WU-BLAST: Washington University BLAST Archives"]}, {"article-title": ["NCBI: The National Center for Biotechnology Information"]}, {"article-title": ["UniProt: The Universal Protein Resource"]}, {"article-title": ["geNORM software "]}, {"surname": ["Rozen", "Skaletsky", "Krawetz S, Misener S"], "given-names": ["S", "H"], "article-title": ["Primer 3 on the WWW for general users and for biologist programmers"], "source": ["Bioinformatics Methods and Protocols: Methods in Molecular Biology"], "year": ["2000"], "publisher-name": ["Totowa,  NJ , Humana Press"], "fpage": ["365"], "lpage": ["386"]}, {"surname": ["Benedito", "Torres-Jerez", "Murray", "Andriankaja", "Allen", "Kakar", "Wandrey", "Verdier", "Zuber", "Ott", "Moreau", "Niebel", "Frickey", "Weiller", "He", "Dai", "Zhao", "Tang", "Udvard"], "given-names": ["VA", "I", "JD", "A", "S", "K", "M", "J", "H", "T", "S", "A", "T", "G", "J", "X", "PX", "Y", "MK"], "article-title": ["A gene expression atlas of the model legume Medicago truncatula"], "source": ["Plant Journal"], "year": ["2008"], "volume": ["(in press)"]}]
Plant Methods. 2008 Jul 8; 4:18