Quantitative RT-PCR based platform for rapid quantification of the transcripts of highly homologous multigene families and their members during grain development
BMC Plant Biology volume 12, Article number: 184 (2012)
Cereal storage proteins represent one of the most important sources of protein for food and feed and they are coded by multigene families. The expression of the storage protein genes exhibits a temporal fluctuation but also a response to environmental stimuli. Analysis of temporal gene expression combined with genetic variation in large multigene families with high homology among the alleles is very challenging.
We designed a rapid qRT-PCR system with the aim of characterising the variation in the expression of hordein genes families. All the known D-, C-, B-, and γ-hordein sequences coding full length open reading frames were collected from commonly available databases. Phylogenetic analysis was performed and the members of the different hordein families were classified into subfamilies. Primer sets were designed to discriminate the gene expression level of whole families, subfamilies or individual members. The specificity of the primer sets was validated before successfully applying them to a cDNA population derived from developing grains of field grown Hordeum vulgare cv. Barke. The results quantify the number of moles of transcript contributed to a particular gene family and its subgroups. More over the results indicate the genotypic specific gene expression.
Quantitative RT-PCR with SYBR Green labelling can be a useful technique to follow gene expression levels of large gene families with highly homologues members. We showed variation in the temporal expression of genes coding for barley storage proteins. The results imply that our rapid qRT-PCR system was sensitive enough to identify the presence of alleles and their expression profiles. It can be used to check the temporal fluctuations in hordein expressions or to find differences in their response to environmental stimuli. The method could be extended for cultivar recognition as some of the sequences from the database originated from cv. Golden Promise were not expressed in the studied barley cultivar Barke although showed primer specificity with their cloned DNA sequences.
Cereal-derived products constitute a major part of human and livestock diets. In 2009 the annual production of all cereals exceeded 2,400 million tonnes . Assuming an average protein content of 10%, over 240 million tonnes of grain protein was harvested in 2009. The average protein content of cereals indicates that relatively small differences exist within and between species and that these can be amplified by environmental factors. For example the protein content on a dry weight basis is in the range of 10–15% in wheat ; 8–15% for barley  or 9.1-13.3% for rye .
Improvement in complex traits such as cereal yield and quality through plant breeding and the associated molecular and biochemical tools is vital to keep pace with population growth and nutritional requirements. Considerable effort has been directed at developing a range of genomic DNA markers with the goal of supporting marker assisted plant breeding. The value of anonymous genetic markers such as random DNA markers (SSRs, AFLPs, RFLPs etc.) depends on the known linkage between marker and target locus alleles . In contrast, “functional markers” (FMs) are derived from polymorphic sites within genes, quantitative trait nucleotides (QTN) or quantitative trait insertion – deletion mutations (QTINDEL) and have significant potential to translate genomic technologies into improved crop varieties. Once genetic effects have been assigned to functional sequence motifs, FMs can be used for fixation of beneficial alleles [6, 7]. Syvänen  has described five phases of functional marker development (1) functionally characterised genes, (2) allele sequences from such genes, (3) identification of polymorphic, functional motifs affecting plant phenotype within these genes, (4) validation of associations between DNA polymorphisms and trait variation, and (5) conversion into technical assays using, e.g., any of the single nucleotide polymorphism (SNP) or INDEL detection technologies. To date, only a limited number of genes have been isolated for functional markers development in wheat  and to the authors knowledge there are no reports of functional marker development with respect to storage proteins in barley, the focus of our research.
The genetic control of grain protein content in barley has been reviewed by Ullrich  and the trait is clearly polygenic, with quantitative trait loci (QTLs) mapping onto all seven chromosomes. Hordeins, the main storage proteins of the barley endosperm are encoded by a single multigenic locus located on chromosome 5 and are divided into four groups: B-hordein (sulphur-rich), C-hordein (sulphur-poor), γ-hordein (sulphur-rich) and D-hordein (high molecular weight), distinguishable by their electrophoretic mobility and amino acid composition. The four groups are encoded by the genes: Hor2 (B-fraction), Hor1 (C-fraction), Hor3 (D-fraction) and Hor5 (γ-fraction), located on barley chromosome 5 (1 H). The B-hordeins account for 70-80% (mol wt 35–46 kDa); the C-hordeins 10–20% (mol wt 55–75 kDa); the D hordeins 2–4% (mol wt 100 kDa) and the γ-hordeins amount <5% (35–46 kDa) of the total hordein fraction . The Hor2 gene family encoding the B-hordeins and Hor1 loci encoding C-hordeins are thought contain between 20–30 genes per haploid genome [11, 12]. The D-and γ-hordein groups, encoded by the Hor3 and Hor5 loci respectively [3, 13], are minor components and little is known of the extent of polymorphism of the genes although their products extensively studied [14, 15].
The variation in barley storage protein allelic complement and the associated contribution to the storage protein profiles both in terms of amino acid and polypeptide composition provides the bases for the observed significant variation within and between that barley cultivars and wild relatives with respect to the number and type of storage proteins/polypeptides [16, 17].
To address the need for studying the expression patterns of hordein genes, a rapid qRT-PCR screening method was developed and verified. Our work provides a tool for identifying presence of the storage protein alleles and their expression in developing barley grain. We assume that it can be also used to study changes in a response to different environmental conditions or be assistance in researching genetic markers and more specifically, functional markers, in barley. Although the impact of cultivar-specific polymorphism could be considered as a limitation of the method, at the same time it could be used to highlight the inherited problems of the primer design when database sequences that originate from many different cultivars are used.
Hordein sequence analysis and primer design
The first step toward developing the high throughput platform was to curate all known hordeins sequences from EMBL Nucleotide Sequence Database (EMBL), DNA Data Bank of Japan (DDBJ), GenBank at the NCBI (GenBank) and HarvEST databases. Incomplete, partial and EST sequences were filtered out of the collection to ensure that only accessions that coded for the full proteins were used. All the accession numbers of the genes used in this study can be found in Table1.
Hordein genes are characterised by repeated sequence motifs and the impact and evolution of this feature has been extensively reviewed [3, 18, 19]. A full discussion of the impact of repeats on evolutionary analysis is beyond the scope of this work. However a typical phylogenetic analysis of sequence data involves five distinct steps: (a) data collection, (b) inference of homology, (c) sequence alignment, (d) alignment trimming, and (e) phylogenetic analysis . As part of the ‘alignment trimming’ steps a core sequence was chosen to improve the quality of the multiple sequence alignment but also assist the primer design. The core sequences were selected by cutting off most of the repetitive regions. The 29 different B-hordein alleles (full lengths are between 762–933 bp) were trimmed to 596 to 609 bp, while the chosen “core” sequences of the three different C-hordein alleles were 640 bp (full lengths are between 867–1017 bp), 1220 bp for the five D-hordein alleles (full length sequences are between 2124–2274 bp) and 649 bp for the different γ-hordein alleles (full lengths are between 768–918 bp) (Additional file 1). The ‘core’ sequences were used to create multiple sequence alignment and design the primers. The phylogenetic analysis was performed on the full length clones, core sequences and proteins, all producing the same arrangements within the families. The results of using the protein sequences were visualised in unrooted phylogenetic tree (Figure1).
The sequence identity of hordein alleles is very high: 84 to 99% for B-, 80 to 92% for C-, 99% for D-, from 98 to 99% for γ1 and 94% for γ3-hordeins, therefore designing specific primers for groups, subfamilies and individual members are challenging. Primers were designed to recognise whole groups and where it was appropriate subfamilies or individual members of the gene families (Figure1; Table2).
B-hordein represents 70-80% of the hordein storage protein, as such is the most significant class of hordeins in terms of amount protein found in the mature grain . The phylogenetic analysis of the B-hordein sequences divided the population into two major subfamilies: B3 and B1, and a minor subfamily B2, these subfamilies could be further subdivided in to groups according to their sequence differences (Figure1). Primer sets were designed to recognise the whole gene family (Common B primer set); the subfamilies (B1-, B2- and B3-), subgroups of subfamilies (B1a; B3a; B3b; B3c) and individual members (GQ342976, DQ826387, JQ867088, DQ148297, X87232, JQ867084) (Figure1; Table2).
In mature barley grains C-hordein accounts for approximately 10-20% of the total hordein protein . However, to date, only three full length sequences have been deposited on the databases from three different cultivars (Table1). Phylogeny analysis of the C-hordein sequences identified two subfamilies and the two groups showed MW differences as well (Figure1; Table1). A common C-hordein primer set was designed to study the expression of the whole family (Common C primer set) and in addition primers specific for the subfamilies were also created (C1, S66938) (Figure1; Table2).
D-hordein accounts for less than 5% of total grain protein in the mature grain . Trawling the databases revealed five highly similar coding sequences corresponding to D-hordein derived from four different cultivars. Given the sequence architecture of D-hordein we designed a common primer set for all the genes. The phylogenetic analysis generated two subfamilies, which we annotated D1 and D2 but were unable to design primers to distinguish the specific alleles as the differences among the sequences were in the number of repeats (Figure1; Table1; Table2). It is suggested that an alternative method such as standard RT-PCR could be used to differentiate the expression of the D-hordein alleles .
Gamma-hordeins are represented by a small group of protein, their contribution to the total grain protein has not been precisely determined . We obtained six sequences from five cultivars and these were apportioned into two subfamilies, γ1 and γ3 by phylogenetic analysis (Figure1; Table1). As the coding sequences of two subfamilies are very different we were not able to design a common primer sets for the family. Common primers were designed for γ1- and γ3-subfamilies and four specific primers were designed for individual subgroups (Figure1; Table2).
Quantitation of the DNA concentration by standard curve using qRT-PCR reaction
In order to construct standard-curves, serial dilutions of known amounts of cloned actin DNA was-amplified in the qRT-PCR reactions. Actin DNA was diluted from 1.44 ng/μl to 1.44 fg/μl (10 times dilution series) (Figure2). The standard was included in every plate and the PCR efficiencies were between 94.25% and 109.49%. From the dilution curves and calculated PCR efficiencies the range of reliable and acceptable Ct values was established between 5 and 28. Based on the slopes for hordein and reference primers we selected the optimal cDNA concentrations to calculate expression level of hordein transcripts.
Validation of primer specificity and standard curve by qRT-PCR reaction using cloned DNA targets
In order to validate the efficacy of allelic specific primers to accurately discriminate between alleles within a cDNA population it was necessary to demonstrate the primer specificity on the DNA template of the corresponding allele. We have 15 different hordein DNA clones representing 9 B-, 1 C-, 2 D-, and 3 γ-hordein alleles. Each primer set, which had been designed towards a specific allele or alleles was tested both on the DNA of the target allele but also in the presence of non-target allelic DNA. A primer set was designated as specific when only target gene or genes were amplified. Validated primers, specific to the selected hordein groups are presented in Table2 and in the Additional file 1. The phylogenetic relationship and the primers sets used are illustrated in Figure3.
We used the above mentioned 15 different hordein DNA clones to validate the qRT-PCR with actin standard curve (Figure2). The validation experiment was performed by adding same amount of DNA to the DNA pool from the individual clones and their amol/μl concentrations were calculated. For example, common B group had 9 individual members with the total 5.16 amol/μl DNA in the pooled DNA sample and the DNA concentration calculated from the Ct value was 5.71 ± 0.18 amol/μl (Table3). B1 and B3 subgroups showed high similarities between the added and calculated DNA concentrations: 2.94 vs. 2.92 ±0.11 and 2.21 vs. 2.41 ± 0.61 amol/μl, respectively. Furthermore the calculated B1 and B3 values added up to the total calculated common B values 5.33 ±0.36 vs. 5.71 ± 0.18 amol/μl (Table3). The method was working as well when individual members like the B1- (JQ867084) or γ1- (JQ867080) hordein were evaluated from the pool. Overall, we obtained similar values when the target DNA concentrations were measured versus calculated from the Ct values using the actin standard curve and these experiments proved the specificity of designed primers and efficacy of the system (Table3).
The stability of reference genes
The genes for actin, ubiquitin or protein translation factor SUI1 homolog (GOS2 protein) were used as internal references to normalise cDNA concentration between grain samples taken at different developmental stages. The fluctuations of the housekeeping genes were calculated relative to each other during grain development (Table4). Actin gene showed the most stabile expression, followed by the ubiquitin gene, while the gene for protein translation factor SUI1 homolog (GOS2 protein) showed higher expression level at 25DAP but stable expression in the earlier stages (Table4).
Differential expression of hordein alleles in developing barley (cv. Barke) grain
Quantitative-PCR experiments were performed to; 1) validate the high throughput platform using the same material described in microarray study conducted by Hansen et al. study ; 2) extend the gene expression study with our newly designed experimental set up (Figure1; Table2). Although three housekeeping genes (actin, ubiquitin and protein translation factor SUI1 homolog) were used to normalise the data acquired across the development series, we presented the results using actin as the reference gene as this gene proved to be the most stabile during in the studied period (Table4).
The number of moles of transcript corresponding to a specific class of hordein alleles expressed during grain development was calculated using a standard curve created using actin as an internal reference and standard. Our work indicated that the genes encoding hordeins, the major storage proteins of barley, were expressed at different levels. Figure4 illustrates that throughout development B-hordein transcripts are the most abundant, followed by C-, γ- and D-hordein transcripts.
Further analysis of the temporal expression data shown in Table5 reveals the percentage contribution of each family to the total amount of hordein transcripts. The total amount of hordein transcripts (100%) was added up from Common B, C, D and γ hordein transcript (amol of hordein/amol of actin) (Additional file 2). The total B-hordein gene expression contribution decreased from 82.86% at 10 DAP to 80.43% at 25 DAP, the total C-hordein gene expression contribution increased from 5.73% to 8.72% over the same period, while D-hordein gene expression levels were 1.55% of the total hordeins at 10 DAP; 2.57% at 15 DAP; 2.48% at 18 DAP and 1.89% at 25 DAP. The γ-hordein group (γ1 and γ3 together) decreased from 9.86 at 10 DAP to 8.96% at 25 DAP (Table5).
The multiple sequence alignments of the different families coding hordein sequences revealed both polymorphisms with respect to nucleotides but also overall length of the clone and therefore predicted number of amino acids (Table1). Studying the phylogenetic relationships of the B-hordein groups suggested two major subfamilies (B1 and B3), one minor subfamily (B2) and several subgroups in the major groups but the MW of the transcript and the coded protein had no correlation with the groups (Figure1; Table1).
Analysing the percentage contribution of subfamilies and subgroups of subfamilies within the B-hordein family revealed the absolute concentration of B1-group transcripts was high at the beginning of development from 35.02% (10 DAP) to 34.94% (15 DAP) and decreased to 15.84% at 25 DAP. The primer set for the B1a subgroup covers a significant proportion of the B1-subfamily and its contribution to this subfamily was 35.09%, 34.92%, 34.62% and 47.03% for 10, 15, 18 and 25 DAP, respectively. The transcripts from B3-group had a similar expression to B1-group at early grain developmental stages but in contrast to B1-hordeins their expression decreased to about 30% at 15 and 18 DAP and significantly increased and reached again over 41% of the total B-hordein fraction at 25 DAP (Table5). The B3a subgroup forms a large part of the B3-subfamily and constitutes ~20% of all B-hordeins at early developmental stages and increased to 34.51% at 25 DAP. The increase of the number of transcripts of the B3a subgroup during late development accounts for the change in transcript level for B3-group as a whole.
D1-hordeins has 4 representatives while the D2 group contains 1 member. No working primer sets were established to enable the characterisation of the individual members of the two subfamilies as the differences were in the number of the repetitive motifs (Figure1; Table5).
Presence and absence of alleles in cv. Barke
Table1 illustrates that the sequences sourced from the databases where actually derived from a wide range of cultivars. It was noted that certain primer sets designed to allelic sequences drawn from the database were not able to detect the corresponding sequence in cv. Barke. We were able to test some of these primer sets isolated from cultivar Golden Promise on their own DNA clones. In spite of the fact that the primer sets supported the production of an appropriate PCR product when Golden Promise DNA was used (Table3) we did not observe a PCR product when tested on the cultivar Barke DNA (Table5). For example, the allele of the B1-subfamily isolated from cultivar Golden Promise was not expressed in cv. Barke (Table5). Similarly, we did not obtain a PCR product from the cDNA of cv. Barke when using primer sets designed for the members of B3b subgroup cloned from cv. Golden Promise (JQ867086, JQ867087, JQ867089) (Table5). Furthermore some of the alleles from other cultivars, for examples B1-subfamily members isolated from cultivars Carlsberg II (X87232) and XQ053 (DQ148297) or γ-sequences from cultivars Carina (X13508), Riso 56 (AJ580585) Carlsberg II (X72628) did not produce a PCR product from Barke cDNA (Table5). Although DNA representing those alleles was not available for testing we suppose that similar results would be obtained, namely the expression pattern is cultivar dependent.
Genetic sequence alignment is the basis of many evolutionary and comparative studies. When creating a multiple sequence alignment of the collated hordein genes, the number and type of repeats within the genes resulted in gaps in otherwise highly homologous sequences. However when performing phylogenetic analysis of the hordein genes using the multiple sequence alignment of full sequences it became apparent that the algorithm removes the gaps when making evolutionary related comparisons [22, 23]. Phylogeny-aware gap placement software has been developed to reduce errors in sequence alignment and evolutionary analysis as Löytynoja and Goldman . However it appears that a typical phylogenetic analysis of sequence data involves five distinct steps, one of which is ‘alignment trimming’  choosing of core areas of genes without extensive repeat supported both the phylogenic analysis and primer design. To study long repetitive regions requires techniques different from those which were described in this study.
The DNA sequence identity of hordein alleles is frequently higher than 95%, therefore designing specific primers for families, subfamilies and individual members is challenging. In the light of the fact the primers set had to be designed to often discriminate single nucleotide changes within alleles we adopted the strategy, which underpins the development of SNP (Single Nucleotide Polymorphism) detection. The SNP detection is based on the ability to discriminate single point mutations and relies on DNA polymerase with proof reading activity, to extend a primer only when its 3'-end is perfectly complementary to the template . Whiley and Sloots  used Taqman probes and combination of variously modified primers and they noticed a destabilizing effect for a single base mismatch in the 3'-end while 5'-end is less likely to introduce error. The observation was considered in our experiment and was confirmed in our results: most of the specific primers have a mismatched nucleotide at the 3'-end (Figure1; Table2; Additional file 1). In contrast to one fluorescence dye (SYBR® Green I) qPCR assay, Taqman probes can be based either on regular oligonucleotides or on Locked Nucleic Acid (LNA) and detect specific reaction products only [27, 28]. However, the method is relatively expensive and a different probe has to be synthesised for each unique target sequence which is a complicated task in the case of large multigene families. In our experiment we proved that SYBR®-Green based detection is sensitive when appropriate primer sets were used and was therefore more economic (Figure1; Table2).
Validation of actin standard curve and primer specificity using cloned DNA targets by qRT-PCR reaction showed the robustness of the method. It is recommended that if possible to check the specificity of designed primers with original cloned DNA template but very often this kind of control is beyond the bounds of possibility for a researcher. Furthermore it is advised to design primer pairs were both of primers of the set can distinguish the SNP. It was observed in our experiment as well that the specificity increased when both primers of the set were unique for desired amplicons (Table2; Additional file 1).
Hansen et al. , using microarray derived data, reported variation in the temporal expression of genes coding for barley storage protein family members within the cultivar Barke. The data sets resulting from microarray were validated by using qRT-PCR and the primers were chosen to recognise most members of the same gene family . Using the same field grown material our principle objectives were two fold 1) quantify the total mRNA of specific classes of hordein during development 2) attempt to dissect out the contribution of subclasses and specific alleles within a class during development. We found substantial fluctuation in the contribution of the different families to the total in the different hordein gene fractions; furthermore we were able to distinguish different contribution of the family members to the total hordein gene pool during grain development.
In the case of B-hordein family, it has been suggested that it is a multigene family with approximately 34 members . Previous reports have identified at least three classes of B-hordein on the basis of their cyanogen bromide (CNBr) cleavage patterns and considerable variation was observed in the numbers and amounts of polypeptides of each class present in different genotypes . The work of Kreis et al.  verified at least two major subfamilies of B-hordein mRNAs, and thus of genes, associated with the Hor2 locus. We established similar major subfamilies by the phylogenic analysis of the currently available 29 sequences (Figure1; Table1). Most of the available sequences from the databases were not annotated according to groups and just described as B-hordein. Our sequences, recently cloned and submitted to the database (JQ series in Table1), have the classification presented in this paper. Although we found no correlation between the groups and the MW of the coded proteins we adhered to the classification of B-hordein as described by Shewry  who distinguished two groups: B1 (class I and II) and B3 (class III) by MW. Further to the two major groups (B1 and B3) suggested by phylogenetic analysis, one minor group (B2) and several subgroups in the major groups were obtained (Figure1). The contribution of B1- and B3-subfamilies to the total pool of B-hordein transcripts differs during grain development however it would appear that the B3-subfamily contributes the greater proportion of total B-hordein content in the mature grain (Table5). Hansen et al.  observed a similar tendency for some of B1- and B3-representatives in their microarray analysis. The total of amount C-hordein gene transcripts increased during the developmental period reaching 8.72% of the total, while D-hordein gene expression levels reached the highest level at 15DAP (2.57%). The amount of transcripts for the γ1-hordein subfamily decreased during the studied period while the percentage of γ3 –hordeins appeared to increase throughout development. It was reported by Rechinger et al.  that there is role for γ3-hordein in the transport and targeting of prolamin polypeptides so they do not serve only as a storage proteins.
The selection of housekeeping genes is critical for gene expression studies. Actin, GAPDH, tubulin, 18S rRNA and heat shock protein 70 are common reference genes for barley [33–35]. In our experiment we verified two new housekeeping genes: ubiquitin and protein translation factor SUI1 (GOS2 protein) homolog (Table1). We checked the stability of these reference genes at different time points and our results implied that the ubiquitin and protein translation factor SUI1 (GOS2 protein) homolog coding genes (Table1) can be used for normalisation when gene expression level is studied in grain growth stages and different barley cultivars.
The often neglected fact of the database sequences that they usually originate from different cultivars. It is true for the available hordein sequences as well; they were cloned from many different cultivars. Our result highlighted that presence of some of the alleles are cultivar and/or developmental stage dependent and it should be considered when gene expression studies are performed.
The method described enabled rapid characterisation of the allelic contribution to the total hordein storage protein transcript population during grain development. Using cheaper SYBR Green labelling in the qRT-PCR reactions was sufficient to distinguish expression levels of large gene families and their members even with high sequence identity. The qRT-PCR validation experiments using cloned DNA targets proved the specificity of designed primers and with the application of the actin standard curve the efficacy of the system was proven as well. We confirmed the stability of the expression of the chosen reference genes during the studied barley grain development period and found substantial fluctuation in the contribution of the different families to the total in the different hordein gene fractions; furthermore we were able to distinguish different contribution of the family members to the total hordein gene pool during grain development. Our result highlighted that presence of some of the alleles are cultivar and/or developmental stage dependent and it should be considered when gene expression studies are performed. The described primer sets could be used as functional marker to help the breeding effort for better storage protein qualities.
Hordeum vulgare L cv. Barke was grown under field conditions as described Hansen et al. . All the plant material was morphologically and chronologically staged in accordance with internationally recognised criteria of Zadoks code . Individual spikes were tagged at flowering and harvested in the morning (09.00-10.00 h) at 10, 15, 18, and 25 d after pollination (DAP). Developing grains were immediately frozen in liquid nitrogen and stored at −80°C until analysis. Two grains were sampled from the middle of a spike. Three spikes per treatment where sampled and the grains pooled before analysis. Each measurement was repeated three times.
DNA, RNA isolation and cDNA synthesis
DNA coding individual hordein alleles and actin gene from barley was prepared from plasmid clones using GenElute Plasmid Miniprep kit (Sigma-Aldrich). DNA was measured using DNA Quantitation Kit, Fluorescence Assay (Sigma-Aldrich).
Total RNA was extracted from milled material according to manufacturer’s protocol (FastRNA Pro Green Kit, Bio101 Systems, France). The isolated RNA was treated with DNase according to the manufacturer protocol (Qiagen) to ensure that all genomic DNA was removed. The RNA was re-isolated with FastRNA Pro Green Kit. RNA quality was checked using an Agilent 2001 Bioanalyzer (Agilent Technologies, Inc.). Samples with RNA Integrity Number (RIN) above 7  were used for mRNA extraction with Dynabeads-Oligo (dT)12–18 according to manufacturer’s protocol (Invitrogen, Norway). First strand cDNA was prepared using 500 ng (500 ng/μl) of Oligo (dT)12–18 primer and Superscript II reverse transcriptase according to manufacturer’s protocol (Invitrogen, USA). The resulting cDNA mixture was diluted to 200 μl by adding 180 μl of MilliQ- H2O and stored at −20°C.
Designing of specific primers for qRT-PCR expression analysis
All available full length gene sequences for hordeins- B, C, D and γ were collected from: EMBL Nucleotide Sequence Database (EMBL), DNA Data Bank of Japan (DDBJ), GenBank at the NCBI (GenBank) and HarvEST database. The accession numbers of the chosen DNA sequences are listed in Table1. The molecular sequences were aligned using online version of ClustalW2 software . The phylogenetic relationships of the different hordein family members were analysed using the software compiled at Phylogeny.fr using standard module [39, 40]. Primer pairs, specific for the whole family, subfamilies of the family and individual sequences were designed manually and their quality was checked by Oligonucleotide primer check software . We were looking for differences between alleles and prioritized primers with mismatched nucleotide at 3′ end. The selected primers are oligonucleotides with a length between 18 and 22 bases and 40-60% of GC content. The amplicon length is between 50 and 150 bp. The selection of actin, ubiquitin and protein translation factor SUI1 homolog (GOS2 protein) as reference genes was based on the report of Sreenivasulu et al.  and the genes were used for normalisation and quantification (Table1). A list of the primers designed towards the hordein gene families and the reference genes is given in Table2.
Quantitative RT-PCR conditions
Quantitative RT-PCR reactions were carried out in triplicate in 384 well microtiter plates (ABI PRISM TM; Applied Biosystems). The total reaction volume was 10 μl which comprised of 5 μl Power SYBR Green Master Mix (Applied Biosystems), 0.5 μM forward and reverse primers (Invitrogen), 1 μl appropriately diluted plasmid DNA, plasmid DNA mix or cDNA. No-template control (NTC) reactions were carried out to check the potential of primer-dimers formation. The qRT- PCR reactions were performed using a 7900HT Sequence Detection System (Applied Biosystems) programmed with the following thermal profile setup: one cycle at 50°C for 2 min; one cycle at 95°C for 10 min; 40 cycles at 95°C for 15 s and 60°C for 1 min. ‘Absolute Quantification’ assay type was used. Data analysis was performed with SDS 2.2.1 software (Applied Biosystems) followed by Microsoft Office Excel 2007 and outliers had been removed.
The DNA standard curve was prepared from a dilution series (10-1 to 10-8) using DNA isolated from the cloned actin gene (HVSMEi0002G07f). The qRT-PCR reactions with primers specific for actin gene were performed in triplicate for each concentration. The PCR efficiency was calculated from the slope of the standard curve according to the following formula: E = 10(−1/slope)-1 where an efficiency of 1 corresponds to 100% . The Ct value obtained for each hordein and reference gene was an average of three PCR reactions on the same cDNA pool. The number of attomoles of individual transcripts was calculated from the weight of specific cDNA derived from the standard curve and corresponding length of coding region (Table1). The number of moles (amol) of a specific hordein were normalised to the number of moles of housekeeping genes (actin or ubiquitin or protein translation factor SUI1 homolog (GOS2 protein). The reference genes and standard curve were present on each plate to detect and remove inter-run variation.
We portrayed the hordein content (amol) per actin according to the days after pollination (Figure4). Due to the wide range of hordein data, we used logarithmic scale for the representation in R 2.15.0 .
Our study conforms to the Minimum Information for Publication of Quantitative Real-Time PCR Experiments (MIQE).
Validation of primer specificity toward target alleles within qRT-PCR reaction
We possess the JQ series (see Table1) of hordein alleles cloned into plasmid vector. To test the specificity of each primer set which had been designed towards a specific allele or alleles qRT-PCR reaction were carried out on a) the DNA template of the target allele(s) b) non-target allelic DNA templates and c) with pooled samples i.e. all alleles. When the target allele was amplified the amplification plot was visually inspected to verify that the Ct values were in optimal range as described previously in the section of qRT-PCR conditions of the M&M. The number of ng amplified by the allele specific primers in the presence of the target template was determined by referencing a standard curve presented in Figure2. The number of amol was calculated from the ng of DNA derived from the standard curve and corresponding length of coding region (Table1).
We measured DNA concentrations of 15 hordein clones as described previously in DNA, RNA and cDNA section and carried out the qRT-PCR reactions with selected primers to check validation of the system. The differences between the number of amol estimated from the amplification product and that actually added to the reaction were used to estimate the amplification efficiency (internal control for the quality of the allele specific primer set) (Table3).
Amplified fragment length polymorphisms
Days after pollination
Locked nucleic acid
Minimum Information for Publication of Quantitative Real-Time PCR Experiments
Restriction fragment length polymorphisms
Single nucleotide polymorphism
Simple sequence repeats
Quantitative reverse transcription-polymerase chain reaction
Quantitative trait insertion – deletion
Quantitative trait locus
Quantitative trait nucleotides.
FAO, statistics. http://www.fao.org/corp/statistics/en/.
Konzak CF: Genetic control of the content, amino acid composition, and processing properties of proteins in wheat. Adv Genet. 1977, 19: 407-582.
Shewry PR: Barley seed proteins. Barley: Chemistry and Technology. Edited by: MacGregor AW, Bhatty RS. 1993, American Association of Cereal Chemists (AACC), St. Paul, Minnesota, 131-197.
Saastamoinen M, Plaami S, Kumpulainen J: Pentosan and [beta]-glucan content of finnish winter rye varieties as compared with rye of six other countries. J Cereal Sci. 1989, 10 (3): 199-207. 10.1016/S0733-5210(89)80049-9.
Lübberstedt T, Melchinger AE, Klein D, Fähr S, Dally A, Westhoff P: QTL mapping in testcrosses of flint lines of maize: III. Comparison across populations for forage traits. Crop Sci. 1998, 38 (5): 1278-1289. 10.2135/cropsci1998.0011183X003800050027x.
Andersen JR, Lübberstedt T: Functional markers in plants. Trends Plant Sci. 2003, 8 (11): 554-560. 10.1016/j.tplants.2003.09.010.
Brenner EA, Zein I, Chen Y, Andersen JR, Wenzel G, Ouzunova M, Eder J, Darnhofer B, Frei U, Barrière Y: Polymorphisms in O-methyltransferase genes are associated with stover cell wall digestibility in European maize (Zea mays L.). BMC Plant Biol. 2010, 10 (1): 27. 10.1186/1471-2229-10-27.
Syvanen A: Accessing genetic variation: genotyping single nucleotide polymorphisms. Nature Rev Genet. 2001, 2 (12): 930-942. 10.1038/35103535.
Bagge M, Xia X, Lübberstedt T: Functional markers in wheat. Curr Opin Plant Biol. 2007, 10 (2): 211-216. 10.1016/j.pbi.2007.01.009.
Ullrich SE: Genetics and breeding of barley quality attributes. Barley Science: Recent advances from molecular biology to agronomy of yield and quality. Edited by: Slafer GA, Molina-Cano JL, Savin R, Araus JL, Romagosa I. The Haworth Press, Binghamton, New York, 2002:115-142.
Bunce N, Forde B, Kreis M, Shewry P: DNA restriction fragment length polymorphism at hordein loci: application to identifying and fingerprinting barley cultivars. Seed Sci Technol. 1986, 14 (2): 419-429.
Shewry P, Bunce NAC, Kreis M, Forde B: Polymorphism at the Hor 1 locus of barley (Hordeum vulgare L.). Biochem Genet. 1985, 23 (5): 391-404. 10.1007/BF00499082.
Shewry P, Finch R, Parmar S, Franklin J, Miflin B: Chromosomal location of Hor3, a new locus governing storage proteins in barley. Heredity. 1983, 50 (2): 179-189. 10.1038/hdy.1983.19.
Gu YQ, Anderson OD, Londeorë CF, Kong X, Chibbar RN, Lazo GR: Structural organization of the barley D-hordein locus in comparison with its orthologous regions of wheat genomes. Genome. 2003, 46 (6): 1084-1097. 10.1139/g03-071.
Pistón F, Shewry PR, Barro F: D hordeins of Hordeum chilense: a novel source of variation for improvement of wheat. Theor Appl Genet. 2007, 115 (1): 77-86. 10.1007/s00122-007-0542-0.
Hansen M, Friis C, Bowra S, Holm PB, Vincze E: A pathway-specific microarray analysis highlights the complex and co-ordinated transcriptional networks of the developing grain of field-grown barley. J Exp Bot. 2009, 60 (1): 153-167.
Han Z, Wu F, Deng G, Qian G, Yu M, Jia Y: Structural and expressional analysis of the B-hordein genes in Tibetan hull-less barley. Genetica. 2010, 138 (2): 227-239. 10.1007/s10709-009-9415-6.
Shewry PR, Darlington H: The proteins of the mature barley grain and their role in determining malting performance. Barley Science: Recent Advances from Molecular Biology to Agronomy of Yield and Quality. Edited by: Slafer GA, Molina-Cano JL, Savin R, Araus JL, Romagosa I. 2002, The Haworth Press, Binghamton, New York, 503-521.
Ullrich SE: Barley: Production, Improvement, and Uses. Wiley-Blackwell, Oxford 2011.
Rokas A: Phylogenetic Analysis of Protein Sequence Data Using the Randomized Axelerated Maximum Likelihood (RAXML) Program. Curr Protoc Mol Biol. 2011, 19: 1-19. 11
Van Pelt-Verkuil E, Van Belkum A, Hays JP: Principles and technical aspects of PCR amplification. 2008, Springer Verlag, Berlin.
McGuire G, Denham MC, Balding DJ: Models of sequence evolution for DNA sequences containing gaps. Mol Biol Evol. 2001, 18 (4): 481-490. 10.1093/oxfordjournals.molbev.a003827.
Bandyopadhyay S, Saha S: GAPS: A clustering method using a new point symmetry-based distance measure. Pattern Recogn. 2007, 40 (12): 3430-3451. 10.1016/j.patcog.2007.03.026.
Löytynoja A, Goldman N: Phylogeny-aware gap placement prevents errors in sequence alignment and evolutionary analysis. Science. 2008, 320 (5883): 1632-1635. 10.1126/science.1158395.
Kwok PY, Chen X: Detection of single nucleotide polymorphisms. Curr Issues Mol Biol. 2003, 5: 43-60.
Whiley DM, Sloots TP: Sequence variation in primer targets affects the accuracy of viral quantitative PCR. J Clin Virol. 2005, 34 (2): 104-107. 10.1016/j.jcv.2005.02.010.
Holland PM, Abramson RD, Watson R, Gelfand DH: Detection of specific polymerase chain reaction product by utilizing the 5'––3'exonuclease activity of Thermus aquaticus DNA polymerase. Proc Natl Acad Sci U S A. 1991, 88 (16): 7276. 10.1073/pnas.88.16.7276.
Braasch DA, Corey DR: Locked nucleic acid (LNA): fine-tuning the recognition of DNA and RNA. Chem Biol. 2001, 8 (1): 1-7. 10.1016/S1074-5521(00)00058-2.
Faulks A, Shewry P, Miflin B: The polymorphism and structural homology of storage polypeptides (hordein) coded by the Hor-2 locus in barley (Hordeum vulgare L.). Biochem Genet. 1981, 19 (9): 841-858. 10.1007/BF00504250.
Kreis M, Rahman S, Forde B, Pywell J, Shewry P, Miflin B: Sub-families of hordein mRNA encoded at the Hor 2 locus of barley. Mol Gen Genet. 1983, 191 (2): 194-200. 10.1007/BF00334813.
Shewry PR: Barley Seed Storage Proteins: Structure, Synthesis, and Deposition. Nitrogen metabolism of plants. (Proceedings of the Phytochemical Society of Europe: 33). Edited by: Mengel K, Pilbeam DJ. 1992, Clarendon, Oxford, 201-227.
Rechinger K, Simpson D, Svendsen I, Cameron-Mills V: A role for γ3 hordein in the transport and targeting of prolamin polypeptides to the vacuole of developing barley endosperm. Plant J. 1993, 4 (5): 841-853. 10.1046/j.1365-313X.1993.04050841.x.
Close TJ, Wanamaker SI, Caldo RA, Turner SM, Ashlock DA, Dickerson JA, Wing RA, Muehlbauer GJ, Kleinhofs A, Wise RP: A new resource for cereal genomics: 22 K barley GeneChip comes of age. Plant Physiol. 2004, 134 (3): 960-968. 10.1104/pp.103.034462.
Burton RA, Shirley NJ, King BJ, Harvey AJ, Fincher GB: The CesA gene family of barley. Quantitative analysis of transcripts reveals two groups of co-expressed genes. Plant Physiol. 2004, 134 (1): 224-236. 10.1104/pp.103.032904.
Walia H, Wilson C, Wahid A, Condamine P, Cui X, Close TJ: Expression analysis of barley (Hordeum vulgare L.) during salinity stress. Funct Integr Genomics. 2006, 6 (2): 143-156. 10.1007/s10142-005-0013-0.
Zadoks JC, Chang TT, Konzak CF: A decimal code for the growth stages of cereals. Weed Res. 1974, 14 (6): 415-421. 10.1111/j.1365-3180.1974.tb01084.x.
Schroeder A, Mueller O, Stocker S, Salowsky R, Leiber M, Gassmann M, Lightfoot S, Menzel W, Granzow M, Ragg T: The RIN: an RNA integrity number for assigning integrity values to RNA measurements. BMC Mol Biol. 2006, 7 (1): 3. 10.1186/1471-2199-7-3.
Larkin M, Blackshields G, Brown N, Chenna R, McGettigan P, McWilliam H, Valentin F, Wallace I, Wilm A, Lopez R: Clustal W and Clustal X version 2.0. Bioinformatics. 2007, 23 (21): 2947-2948. 10.1093/bioinformatics/btm404.
Dereeper A, Guignon V, Blanc G, Audic S, Buffet S, Chevenet F, Dufayard JF, Guindon S, Lefort V, Lescot M: Phylogeny. fr: robust phylogenetic analysis for the non-specialist. Nucleic Acids Res. 2008, 36 (suppl 2): W465-W469.
Dereeper A, Audic S, Claverie JM, Blanc G: BLAST-EXPLORER helps you building datasets for phylogenetic analysis. BMC Evol Biol. 2010, 10 (1): 8. 10.1186/1471-2148-10-8.
The OLIGONUCLEOTIDE PRIMER CHECK software. http://depts.washington.edu/bakerpg/primertemp/primermelttemp.html.
Sreenivasulu N, Usadel B, Winter A, Radchuk V, Scholz U, Stein N, Weschke W, Strickert M, Close TJ, Stitt M: Barley grain maturation and germination: metabolic pathway and regulatory network commonalities and differences highlighted by new MapMan/PageMan profiling tools. Plant Physiol. 2008, 146 (4): 1738-1758. 10.1104/pp.107.111781.
Pfaffl MW: A new mathematical model for relative quantification in real-time RT–PCR. Nucleic Acids Res. 2001, 29 (9): e45. 10.1093/nar/29.9.e45. e45
R Development Core Team: R: A language and environment for statistical computing. 2012, R Foundation for Statistical Computing, Vienna, Austria, http://www.R-project.org/.
We would like to thank K.B. Nellerup and O.B. Hansen for their excellent technical support. The work was supported by grant from The Danish Directorate for Food, Fødevareforskningsprogrammet 2008.
1Department of Genetics and Biotechnology, Aarhus University, Research Centre Flakkebjerg, Forsøgsvej 1, Slagelse, DK-4200, Denmark. 2Verzyme (UK) Ltd., Plas Gogerddan, Aberystwyth, Wales, SY23 3EB, United Kingdom. 3MTA-ELTE-MTM Ecology Research Group, Biological Institute, Eötvös Loránd University, Pázmány Péter sétány 1C, Budapest, H-1117, Hungary.
The authors declare that they have no competing interests.
AK isolated the RNA, made the cDNAs, designed the primers, carried out the qRT-PCR experiments and evaluated the data. EZ helped with the statistical analysis. SB and EV participated in the planning the experiments, evaluating the data and was involved writing the article together with AK. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Alignment of B-, C-, D- and γ-hordein sequences found in NCBI databases. A location of qRT-PCR primers is highlighted. (PDF 200 KB)
About this article
Cite this article
Kaczmarczyk, A., Bowra, S., Elek, Z. et al. Quantitative RT-PCR based platform for rapid quantification of the transcripts of highly homologous multigene families and their members during grain development. BMC Plant Biol 12, 184 (2012). https://doi.org/10.1186/1471-2229-12-184