Skip to main content

TILLING in the two-rowed barley cultivar 'Barke' reveals preferred sites of functional diversity in the gene HvHox1



The economic importance of cereals such as barley, and the demand for improved yield and quality require a better understanding of the genetic components that modulate biologically and commercially relevant traits. While Arabidopsis thaliana is the premiere model plant system, the spectrum of its traits cannot address all of the fundamental questions of crop plant development. Unlike Arabidopsis, barley is both a crop and a model system for scientific research, and it is increasingly being used for genetic and molecular investigations into the conserved biological processes of cereals. A common challenge in genetic studies in plants with large genomes arises from the very time-consuming work of associating mutant phenotypes with gene sequence information, especially if insertion mutagenesis is not routine, as in barley. Reverse genetics based on chemical mutagenesis represents the best solution to this obstacle.


In barley, we generated a new TILLING (Targeting Local Lesions IN Genomes) resource comprising 10,279 M2 mutants in the two-rowed malting cultivar 'Barke,' which has been used in the generation of other genomic resources in barley (~150,000 ESTs, DH mapping population). The value of this new resource was tested using selected candidate genes. An average frequency of approximately one mutation per 0.5 Mb was determined by screening ten fragments of six different genes. The ethyl methanesulphonate (EMS)mutagenesis efficiency was studied by recording and relating the mutagenesis-dependent effects found in the three mutant generations (M1-M3). A detailed analysis was performed for the homeodomain-leucine-zipper (HD-ZIP) gene HvHox1. Thirty-one mutations were identified by screening a 1,270-bp fragment in 7,348 M2 lines. Three of the newly identified mutants exhibited either a six-rowed or an intermedium-spike phenotype, and one mutant displayed a significantly altered spikelet morphology compared to that of the 'Barke' wild type. Our results indicate a bias in the frequency of independent functional mutations at specific base pair (bp) positions within the gene HvHox1.


A new TILLING population was developed as a resource for high-throughput gene discovery in an alternative barley germplasm. Pilot screening demonstrated a similar or even slightly higher mutation frequency when compared to previously published barley TILLING populations that should allow for the identification of diverse allelic variation. Partial phenotypic evaluation of the M2 and M3 generations has revealed the presence of a wide spectrum of morphological diversity that highlights the great potential of this resource for use in forward genetic screens. Altogether, our study shows the efficiency of screening and the applicability of the new TILLING population for genetic studies in the barley crop model system.


Barley has a long history as a model plant in mutation research and breeding [1]. One year after Muller published his pioneering work on mutations induced by X-rays in fruit flies [2], Stadler published the first results on induced mutations in barley [3]. In the same year, the Swedish barley mutation research programme was initiated to study the fundamental effects of different mutagens on the barley genome and their applicability for breeding purposes [4]. Since spontaneous mutations occur with extremely low frequency, mutation induction techniques were widely expected to rapidly increase variability in crop species. After initial attempts, more than ten years passed before the first practical mutation procedures were implemented leading to improved mildew resistance in barley [57]. In 1959, a conference at today's Leibniz Institute of Plant Genetics and Crop Plant Research (IPK) substantially promoted interest in chemical mutagenesis [8]. Today, over 300 barley varieties are listed that were either directly or indirectly derived from mutation breeding [9], reflecting the impact of mutation induction on crop improvement.

In the past, the systematic development of mutagenesis was limited by the lack of effective mutation screening techniques and the rather basic knowledge of genes underlying the designated traits, rather than by the efficiency of mutation induction [10]. The introduction of TILLING (for Targeting Local Lesions IN Genomes) [11], which combines chemical mutagenesis with high-throughput genome-wide screening for point mutations in genes of interest, effectively complements classical forward mutation screening based solely only on phenotype. Chemical mutagens like ethyl methanesulphonate (EMS) induce high-density single nucleotide changes that are randomly distributed in the genome [12]. These DNA-damaging agents induce allelic variation that includes missense mutations in coding sequences (leading to altered protein structure and function), mutations of non-coding sequences that may affect gene regulation, as well as premature stop codon mutations and splice site changes that may result in knockout/knockdown of a protein. TILLING introduces new aspects into plant breeding, such as the possibility of surveying specifically for allelic series of economically interesting genes [13, 14]. Moreover, as a highly precise and targeted approach, TILLING does not involve the generation of genetically modified plants, and hence it is a more globally accepted tool for breeding.

In the present study, we describe a new EMS-induced mutant population comprising 10,279 M2 individuals of the two-rowed spring-type barley cultivar 'Barke.' TILLING was used for ten DNA fragments derived from six different genes in order to estimate the average mutation frequency for this population. A more detailed analysis was performed for the gene HvHox1, which was previously shown to control the row-type character in barley [15], allowing for a frequency bias for functional mutation sites to be determined. Our results suggest that the new TILLING resource of the barley cultivar 'Barke' has potential for use in fundamental research as well as for applied breeding.


Generation of a TILLING population in the barley cultivar 'Barke'

The two-rowed barley cultivar 'Barke' was used to develop a mutant population. The response to treatment with ethyl methanesulphonate (EMS) had to be determined first, because the genotype-dependent sensitivity of 'Barke' to chemical mutagenesis was unknown. Mutagenic efficiency was used as the method of validation. Based on previous experience with EMS mutagenesis with TILLING in a different barley cultivar ('Optic' [16]), nine EMS concentrations between 20 (0.2%) and 60 mM (0.63%) were applied. Relatively narrow increments of 5 mM were chosen since chemical mutagens are characterised by specific threshold values, after which the survival rates may quickly decline [17]. The mutation efficiency was monitored as the relative rate of mutation frequency in comparison to the frequency of undesirable changes at a given mutagen concentration. The applied range of concentrations was expected to allow differentiation between the most efficient as well as the maximum tolerated EMS treatment for the cultivar 'Barke.'

As a first parameter, the tolerance of 'Barke' against EMS toxicity was analysed by comparing the relative fitness as a function of the applied EMS concentration; the assessment was based on the respective frequencies of germination, sterility and M2 seed yield (Table 1). Germination rates significantly decreased from 92% for 20 mM EMS to an average of 40.5% for 60 mM EMS. In the same range of EMS concentrations, M1 spike sterility was scored in three categories: completely sterile spike (100% sterility), semi-sterile spike (measured as > 50%, 50%, < 50% sterility) and fully fertile spike (0% sterility). At increasing EMS concentrations, a clearly significant negative relationship was observed with M1 spike fertility leading to nearly 100% sterile spikes at 60 mM EMS (Table 1). For practical reasons in the context of population development, treatment with 40 mM EMS was found to be the upper limit, as the level of fully sterile spikes increased from 63% at 40 mM to 84% at the next higher concentration of 45 mM. Consequently, the overall number of M1 plants carried to the next generation and that were used to build up the overall population varied substantially between treatments (Table 2). In summary, from 99,879 mutagenised barley seeds (M0), 24,600 M1 plants advanced to the M2 generation.

Table 1 Average germination and sterility rates (%) of M1 plants observed after different EMS Treatments
Table 2 M2 population structure: number of M2 lines contributing to the population from different EMS concentrations

Plant fertility was also assessed in the M2 generation since efficient seed set in M2 is a prerequisite for effective phenotyping of M3 families identified after TILLING. The M2 seed yield was negatively correlated with EMS dosage (r = -0.93 and r = -0.99). From the 20 mM-treated M2 mutants, on average more than 90 seeds were harvested regardless of whether the plants were grown under greenhouse or field conditions. Only 65% of the plants receiving 40 mM EMS produced M3 seeds, and 38% of the plants with the highest EMS concentration yielded M3 seeds that could be successfully harvested (data not shown).

Based on prior experience with reduced fertility in advanced mutant generations, ten plants per M2 family were initially cultivated from 15,346 M1 plants to construct the final TILLING population. After tillering, the M2 families were evaluated for viability, and one M2 individual from each family was selected for further cultivation and subsequent development of M2 DNA and M3 seed stocks. This selection was done in order to avoid redundancy of mutation events in the final population. For barley, it was assumed that only one primordial germ line cell of the embryonic main shoot apex was present in M1 [4]. Therefore, M1 seeds of a single spike could share the same mutational event, and a single mutation could appear as a cluster in the segregating spike progenies [18]. The main selection criterion at this stage of population development was general vitality (reduced risk of growing sterile plants), since high yields of M3 seeds were targeted as an essential prerequisite for the development of a useful and sustainable TILLING resource. In some cases (1.7% of all M2 families), when no clear decision could be made for the selection of a single M2, two plants of the same M2 family were kept until maturity. Such siblings were further handled as independent mutant line accessions. Finally, a set of 10,279 M2 plants derived in different proportions from various EMS concentrations (Table 2) was selected to construct the overall population.

For the efficient development of a mutant population, it is desirable to get the earliest possible assessment of the achieved mutation density. Since results of previous barley mutagenesis experiments revealed the correlation between EMS dosage and M2 chlorophyll seedling mutant frequencies [19], M2 families were evaluated for such early phenotypic markers. The analysis focused on the five most frequent phenotypes: albina, xantha, viridis (all three belonging to sub-class unicolour [20]), viridoalbina and striata (both belonging to sub-class bicolour [20]) (Additional file 1). The frequency of such phenotypes was determined in progeny obtained at EMS concentrations ranging from 20 to 50 mM (Table 3). The relative frequencies were calculated based on the number of mutant individuals per 100 germinated M2 plants. Although this method does not consider cluster mutations within the M2 families, compared to spike progeny-based methods [20] it has the advantage of being insensitive to variations in the respective spike progeny sizes, which are affected in M2 families especially at higher EMS treatments (data not shown). The albina phenotype was most frequently observed (Table 3), and a significant relationship between EMS treatment and frequency of chlorophyll-deficient mutants was exclusively seen for this phenotype. The other four investigated phenotypes seemed to be distributed randomly, suggesting that those mutations might be conditional and vary depending on the environment, similarly to the results seen for several chlorophyll seedling mutants in previous studies [21, 22].

Table 3 Average M2 frequencies of chlorophyll seedling mutants (%) after different EMS treatments

Screening for morphological mutants (forward genetics) in the TILLING population

In total, 12,703 M2 individuals (EMS dosages between 20 and 35 mM) were scored for their phenotype during field (11,661 M2) or greenhouse (1,042 M2) cultivation. All plants were observed at two-week intervals from the seedling stage to full maturity with specific emphases on altered spike morphology, heading date (and plant development in general), plant colour, plant height, and other obvious variations form the 'Barke' wild-type. About 20% of the surveyed M2 plants displayed noticeable phenotypes, which were classified into six morphological categories (Additional file 2A). The most common category, plant height (Additional file 2A), contained 45% of all visible phenotypes and mainly comprised all different kinds of dwarfism. Among the 238 M2 plants classified as plant colour mutants, 21% (50 of 238) were impaired in surface wax structure, including such types as glossy sheath and eceriferum. Another 18% (43 of 238) were assigned as viable Chlorina seedlings, showing light- or yellow-green leaves and a stunted growth habit, and these comprised the second most frequent type of the plant colour mutant category.

In addition to phenotypic screening in M2, 1,200 M3 families (derived from 20 to 40 mM EMS treatment) were screened as batches of 1,000 and 200 families, respectively, under field conditions from the subsequent years 2006 and 2007. Screening and classification were carried out as described for the M2 population (Additional file 2B). Approximately 37% of all analysed M3 families showed a mutant phenotype; of these, 5.6% were homozygous for the observed trait, and another 5.6% of the families were segregating two independent phenotypes (data not shown). In 2.8% of the families, plants were observed that displayed more than one aberrant phenotype that could be ascribed to known monogenic mutations (data not shown). In contrast to the phenotypic spectrum observed in the M2 generation, alterations of tillering and plant colour were observed more frequently, while the frequency of dwarfism was notably lower (Additional file 2B).

The frequency of visible morphological phenotypes correlated positively with EMS dosage. In the M2 generation, this frequency ranged from 10% among 20 mM-treated plants to 23% among plants receiving 35 mM EMS. By contrast, in the M3 population, this ranged from 18% among 20 mM-treated plants to 63% among plants receiving 35 mM EMS.

Determination of molecular mutation frequency--TILLING of candidate genes (reverse genetics)

In order to survey the mutation frequency of the composite population at the molecular level, ten fragments from six genes were screened. This work was performed in parallel to the late stages of population development, so only a limited set of 7,389 M2 lines could be used for TILLING (Table 4) (except for the gene mlo for which only 1,920 M2 were surveyed). The mutation frequency per gene fragment was estimated by dividing the total base pairs screened by the total number of mutations identified. The density of mutations was calculated in the established barley population as an average frequency of all screened fragments. Scoring of mutations located either in the proximal or distal 50 bp of each amplicon was compromised because of priming and systematic artefacts on the LiCOR gels. Therefore, in order to calculate the effective number of base pairs screened, each amplicon length was reduced by 100 bp. Gene fragments ranging from 321 bp to 1,400 bp in length were designed to cover either whole open reading frames (ORFs) or selected coding regions of the respective target genes.

Table 4 List of TILLING targets, size of amplicons, and the number and distribution of mutation types

The screen revealed a total of 81 independent mutations (Table 4), including 64 (79%) located in exons. Twenty-nine (45%) of the exon-located mutations induced a change of amino acid (AA) sequence, and three (5%) affected splice sites or inserted a premature stop codon, which most likely led to truncated reading frames. Based on all of the screened fragments, the average mutation frequency for the population was found to be one mutation per 0.5 Mb. Therefore, the entire population of 10,279 M2 lines translates to an average of 20 mutations for any 1 kb DNA fragment. The probability of having at least one truncation mutation is about 60% [1-(1-0.04)20; based on the overall truncation frequency of 4% (3/81)].

Since the TILLING population was derived from seed batches treated with different EMS concentrations, the results could be correlated to EMS dosage [mutation frequency [EMS] = ∑ SNPs per gene/∑ screened Mb per gene]. The single base pair substitution rate was positively correlated (rs = 0.94) with increasing EMS concentration (Figure 1) and ranged (based on TILLING of 8 PCR-derived DNA fragments from 5 genes) from 1 mutation/0.8 Mb at 20 mM EMS to 1 mutation/0.20 Mb and 1 mutation/0.14 Mb after treatment with 40 and 50 mM EMS, respectively. Likewise, an increase in the respective levels of missense mutations was observed between 20 mM- and 40 mM-plants, while there were not any non-synonymous mutations observed with the highest treatment (Figure 1).

Figure 1

Mutation frequency in relation to EMS concentration. The relative SNP rates for EMS treatments between 20 and 50 mM were calculated as particular sub-population specific frequencies. Data are based on the mutation screening of eight gene fragments in 7,389 M2 lines. The white and black bars show the accumulated synonymous and non-synonymous (AA changes and truncation) mutations detected for each sub-population. The relationship between SNP frequencies and EMS dosages was calculated based on a non-parametric Spearman rank correlation (coefficient of rs = 0.94 with a significance at p ≤ 0.01).

New alleles and location bias of functional mutations in the gene HvHox1

A detailed analysis was performed for the HvHox1 gene (homeodomain-leucine-zipper (HD-ZIP) containing transcription factor), which has been shown to encode the row-type controlling locus Vrs1 in barley [15]. Wild-type barley has a two-rowed ear morphology where only the central of three sessile spikelets is fertile and produces seeds. The recessive allele vrs1 affects the fertility of the lateral spikelets, leading to a six-rowed ear. The cultivar 'Barke' possesses a two-rowed ear morphology, so TILLING has the potential to survey for de novo functional alleles with six-rowed phenotypes.

A 1,270-bp region covering the entire open reading frame (ORF) of HvHox1 was screened in 7,348 M2 lines by analysing two overlapping fragments, giving a total of 9.4 Mb of gene sequence analysed. Thirty-one alleles, including ten missense mutations and one splice junction mutation, were identified (Additional file 3). The majority of these mutations (25 of 31) were GC to AT transitions, and 68% (21 of 31) of all mutations in HvHox1 were found in the homozygous state. The average mutation frequency was 1 mutation/0.25 Mb (Additional file 3; Figure 2A). A comparable departure from the expected 1:3 homozygous mutant versus heterozygous and homozygous wild-type ratio for M2 plants was reported by TILLING in the barley cultivar 'Morex,' and it was interpreted as a possible sensitivity limitation for the detection of heterozygotes in 8-fold DNA-pools [23]. Seven missense mutations induced changes to AAs with altered physicochemical properties in regards to polarity and/or hydrophobicity. Four of the missense mutations were detected in the conserved and functionally relevant homeodomain-leucine zipper motif (HD-ZIP).

Figure 2

Tilling of the gene HvHox1 reveals new alleles of the two-rowed/six-rowed locus Vrs1. (A) Model of the gene HvHox1 (Vrs1) designed from output files of the software PARSESNP [58]. Gene structure and distribution as well the type of mutation discovered in a 1,270-bp region covering the entire open reading frame (ORF) of the gene are illustrated: exons = red open boxes, introns = red lines, conserved homeodomain-leucine zipper I-class homeobox domain (obtained from the BLOCKS database, [67]) = green boxes above the gene structure model. Mutations are indicated as small arrowheads below the gene structure model: purple = non-functional mutations, i.e., intron-located or synonymous mutations; black = non-synonymous mutations; and red = a splice junction mutation. Numbered arrowheads indicate mutations leading to phenotypic changes in M2 individuals 11910-1 (1), 3930-1 (2), 8408-1 (3), and 11657-1 (4), respectively. Detailed views of total inflorescences (B-E) and lateral spikelet morphology (F-I) are given for wild-type (two-rowed) 'Barke' (B, F), a six-rowed (hex-v) mutant (8408-1) (C, G), an Intermedium spike-d (Int-d) mutant (11657-1) (D, H), and mutant 3930-1 (E, I). Mutant 8408-1 possessed fully fertile and awned lateral spikelets (G), whereas the lateral spikelets of mutant 11657-1 exhibited shorter awns and only partial fertility, mostly in the basal lateral spikelets (H). The lateral spikelets of 3930-1 were significantly enlarged and tip-pointed compared to 'Barke,' but were sterile and not awned (I).

For all of the M2 plants that exhibited a non-synonymous mutation, 16 to 20 M3 individuals were grown for phenotyping of the row-type character. The progeny of each line were subsequently genotyped. In three of the four missense mutations affecting the HD-ZIP domain, a change in spike appearance, especially in the lateral spikelets, was observed (Figure 2). Mutant line 3930-1, which carried a G to A transition at the 1,039-bp position, exhibited significantly enlarged [ANOVA F-values of 334*** (spikelet length) and 614*** (spikelet width) at a significance value of p ≤ 0.001] and tip-pointed but sterile lateral spikelets without awns in the M3 progeny (Figure 2E, I). The mutant line 8408-1 (T1079A) showed a hexastichon (hex-v) six-rowed phenotype (Figure 2C, G), and the mutant 11657-1 (G1115A) displayed an Intermedium-d spike (Int-d) phenotype (Figure 2D, H). Mutant 11657-1 was the only case where a conservative substitution (Arg to His) led to a visible impact. Finally, a second hexastichon (hex-v) six-rowed phenotype was obtained from mutant line 11910-1 carrying a splice junction mutation at the 3' end of intron 1.

All phenotypically relevant missense mutations were found in the functionally important homeodomain (HD) region of the HvHOX1 protein, while mutants with AA changes outside the HD-ZIP domain displayed normal two-rowed spikes. The four HvHox1 mutants newly identified in the 'Barke' TILLING population that displayed a mutant phenotype were compared to a set of 18 missense alleles previously induced and reported for the cultivars 'Bonus,' 'Foma,' and 'Kristina' [15]. Among all 22 mutant lines, 9 of the 21 highly conserved AAs of the homeodomain domain were affected (Figure 3).

Figure 3

Amino acid conservation in the HvHOX1 homeodomain. AA conservation in the HvHOX1 homeodomain was assessed by a sequence comparison of 100 different homeodomains (homeodomain smart accession No. 00389) provided by NCBI's Conserved Domain Database (CDD) [61, 62]. Highly conserved AAs are highlighted in red, whereas less conserved AAs are highlighted in blue. Positions of the AAs in each protein are indicated with numbers flanking the corresponding peptide, and AAs are given in single letter code. Asterisks above the HvHOX1 sequence indicate induced and non-synonymous substituted AAs detected in 21 barley mutants. The corresponding mutant lines together with the respective AA substitutions are shown in the boxes above. Mutations obtained from the 'Barke' mutant population are depicted in bold letters. The remaining mutations originate from alleles of the cultivars 'Bonus,' 'Foma' and 'Kristina,' as previously described by [15]. The homeodomains of the proteins Antennapedia (Antp) (Schistocerca americana, American grasshopper) and Engrailed (En) (Drosophila virilis) were used as references for the localisation of the functionally relevant AAs in HvHOX1, as previously described [2932]. The main structural features are indicated by α, core AAs involved in the determination of the three-dimensional protein structure; and β, AAs implicated in DNA binding. A model of the secondary structure of the HD peptide is indicated below the alignment, beginning with the N-terminus.

Some of the new 'Barke' TILLING mutations coincided with AA substitutions in previously induced vrs1 alleles (Figure 3). A guanine to adenine transition was found at bp position 1115 in mutant 11657-1, leading to a change from Arg to His. The same nucleotide was affected by a G to T transversion (Arg substituted by Leu) in the mutant lines hex-v.42 and 43 induced in the cultivar 'Kristina' by fast-neutron and EMS treatments, respectively [15].

'Barke' mutant 8408-1 shared an A to T transversion at bp position 1,079 (L95Q) with previously identified mutants hex-v.28 (ethylene imine (EI), cultivar 'Foma') and Int-d.68 (EMS, cultivar 'Kristina') [15]. All three of these mutant lines carried an identical AA substitution, raising the question of whether the newly observed 'Barke' allele was induced independently or if it originated from pollen/seed contamination by a common allele of six-rowed barley cultivars. To test this, we analysed the entire DNA sequence of the 'Barke' Vrs1 allele. Instead of having one of the three haplotypes vrs1.a1, vrs1.a2 and vrs1.a3 common for modern six-rowed cultivars, the HvHox1 gene of the mutant 8408-1 exhibited the haplotype of the two-rowed genotype 'Barke' (except for the transversion at bp position 1,079). Furthermore, we never observed any six-rowed plant among the > 90,000 M1 'Barke' plants, supporting the hypothesis of a newly derived allele. Furthermore, a 'Barke' genetic background was confirmed for 8408-1 by genotyping 25 EST-based SSR markers (data not shown). Altogether, we were able to rule out that the allele of 8408-1 was due to any kind of contamination and to confirm that the 8408-1 allele was derived from 'Barke' by EMS mutagenesis

The alleles of the mutants 3930-1 and 11910-1, which lead to enlarged pointed lateral spikelets or a true hex-v phenotype, respectively, also turned out to be novel vrs1 alleles, as they involved previously unaffected AAs and nucleotide positions. Interestingly, the splice-junction mutation of 11910-1 (G856A) was separated only by 1 bp from a previously described mutation at bp position 855 that was induced in the genotype 'Kristina' by isopropyl methanesulphonate, leading to mutant hex-v.38.


Genetic mutation is a powerful tool that establishes a direct link between the biochemical function of a gene product and its role in vivo. For decades, genes have been identified and functionally characterised by the isolation and study of mutants that are defective in a specific process of interest. With the advent of molecular biology, several reverse genetic techniques have been developed for the functional characterisation of cloned genes. In plants, the most commonly used reverse genetic approaches are post-transcriptional gene silencing (PTGS) [24] and insertional mutagenesis [25]. However, PTGS is labour-intensive and can give ambiguous results. For insertional mutagenesis, the frequency of mutations per genome is typically low, and the mutant alleles are likely to result in a complete loss-of-function of the gene, which might limit the analysis if the effect is lethal or detrimental. By contrast, the TILLING technology induces diverse allelic variation, including missense mutations that can lead to altered protein function, as well as premature stop codons and splice site changes that can result in knockout or knockdown of a gene [26].

TILLING of the homeodomain-leucine-zipper (HD-ZIP) gene HvHox1 [15] in our newly developed 'Barke' cultivar TILLING population demonstrated the feasibility of generating an allelic series of a gene in barley with subtle phenotypic variation. In our study, one splice-junction mutation and three missense mutations were identified that provided three different types of mutant spike morphology. Almost all phenotypically relevant missense mutations of this and previous studies [15] were found to be situated in the functionally relevant homeodomain (HD) region. Those mutants displaying a full or intermediate six-rowed phenotype affected AAs that are universally conserved within members of the HOX protein family [27]. The different mutant phenotypes might be explained not only by the distinct positions of the substitutions in the protein sequence, but also by the character of the amino acid replacement. HOX proteins are DNA-binding transcriptional regulators that contain a highly conserved homeodomain, which is characterised by three α-helical regions (helices 1 - 3) that fold into a tight globular bundle to form a so-called helix-turn-helix (HTH) motif [28]. The AAs of the HvHOX1 homeodomain with putative functional relevance can be predicted by referring to the ANTENNAPEDIA (Antp; AAB03236) and ENGRAILED (En; XP_002050130) protein sequences, a method applied in previous studies [2932]. The arginine at position 107 (Arg107), which is predicted to be involved in direct DNA binding by water-mediated hydrogen bonds between the homeodomain and DNA [31, 32], was affected in two previously identified mutants as well as in one of the 'Barke' mutants. In all three cases, the mutation led to phenotypic changes, but of varying severity. The hex-v.42 and hex-v.43 mutant alleles (originating from the 'Kristina' cultivar) were altered in hydrophobicity (Arg to Leu) and displayed a fully six-rowed phenotype (Takao Komatsuda, unpublished data). A conservative substitution (Arg to His) was identified in the 'Barke' mutant 11657-1, and this led to an Intermedium spike-d (Int-d) phenotype with less developed, short-awned and fully sterile lateral spikelets.

Nine of the 21 (43%) strictly conserved AAs of the homeodomain of HvHox1 were affected by substitutions in previously analysed phenotypic mutants [15]. From the 'Barke' population, we obtained a total of nine mutants that were affected in this region; two mutants, 11657-1 and 8408-1, contained mutations at AA positions identical to induced alleles previously published. This coincidence was surprising but not necessarily unexpected. The probability of detecting a mutation at the same AA of at least one of the previously found substitutions is 99% [1-(1-0.43)9].

The study of mutants originating from different germplasms can provide added value in the functional analysis of a trait or gene. However, it needs to be kept in mind that phenotypic differences of mutants at the same position of the gene may also be the result of modifying loci. Row-type in barley is controlled by the major gene vrs1 (HvHox1), but interaction with the unlinked intermedium (int)-c locus has been demonstrated. Plants heterozygous for vrs1 exhibited different degrees of fertility and seed size in lateral spikelets within or between different six-rowed plants in the presence of different int-c alleles [3335]. Therefore, even mutations in the same AA of the HvHOX1 protein can lead to different phenotypes depending on the int-c allele. As a consequence, the availability of TILLING resources in different germplasms can be important while studying specific traits. Our TILLING population was generated in the two-rowed central European malting variety 'Barke,' whereas previously published TILLING populations were established either in the two-rowed Scottish malting cultivar 'Optic' or the six-rowed North-American malting variety 'Morex' [16, 23]. TILLING of HvHox1 in the latter population would not have been feasible.

Access to TILLING resources in different genetic backgrounds could also be of importance if mutant traits will be utilised in plant breeding. Every individual of a TILLING population typically accumulated a high density of mutations. Given an average mutation frequency of 1 mutation/0.5 Mb, a haploid genome size of 5,000 Mb for barley [36], and an overall number of 40,000 barley genes with an average ORF length of 2.5 kb [37, 38] at an average G/C content of approximately 50% [39, 40], the presence of approximately 100 induced SNPs that affect ORFs can be assumed for every plant of the TILLING population [Total SNP in ORFs = 40,000 ORFs × 1.25 kb~G/Ccontent per ORF × 1 SNP /500 kb]. Of these, about 40% will be non-synonymous mutations, and perhaps one fourth of these will be homozygous. Therefore, each M2 plant will theoretically contain at least ten profoundly affected genes per genome. This background load of mutations is often seen as a major obstacle for utilising TILLING-derived mutant alleles in plant breeding. The background mutational load needs to be removed by recombination in multiple generations of backcrossing or cross-breeding [41]. After backcrossing for six generations, theoretically only 1.5% residual heterozygosity will remain in a mutant, so the true efficiency of removing linked, undesired mutant variation may differ significantly depending on the region of the genome. Linkage disequilibrium (LD), the effect that closely linked genes will be inherited as haplotype blocks at higher than statistically expected probability, may extend to modern barley germplasms over distances of more than 50 cM [42, 43]. Notably, LD is typically high in the centromeric regions of barley chromosomes where recombination frequency is low [44]. However, the impact of LD is not only relevant for predicting the efficiency of removing background mutational load. The presence of natural diversity in a detrimental modifier locus in LD to a desired mutant locus could have a similar effect. In such cases, the principle availability of induced allelic diversity in different genetic backgrounds may facilitate the use of TILLING mutants in a breeding program.

A high mutation frequency is of paramount importance for the development of a TILLING population because it is the key factor influencing the necessary effective population size and thus also the labour investment required for screening. A moderate mutation frequency of approximately one mutation per 0.5 Mb, which is characteristic of the present 'Barke' population as well as of previously published barley TILLING populations [16, 23], requires a relatively large number of lines (~10,000) to achieve genome saturation. Therefore, access to barley TILLING populations featuring high mutation frequencies comparable to those reported for Arabidopsis [12, 45] is an alluring prospect. Interestingly, during our efficiency analysis, we were able to demonstrate such high mutation frequencies for barley because we encountered levels of approximately one mutation per 0.2 and 0.14 Mb among mutants obtained from 40 and 50 mM EMS treatments. Although promising, the reproducibility of such high mutation frequencies needs to be confirmed in future experiments since the subpopulations obtained from 40 mM and 50 mM EMS concentrations were only 4% and 1.5% of the overall population, respectively.

Indeed, high mutagen dosages can impose practical problems, and improved screening efficiency may be compromised due to undesirable effects [46] since high mutagen dosage is typically associated with reduced fitness and fertility in advanced generations. In EMS-derived M1 plants, these effects are mainly caused by the toxicity of EMS metabolites, which leads to reduced M1 seedling germination, decreased growth and survival rates, and partially reduced M1 spike fertility. Over 80% fully sterile M1 'Barke' plants were obtained using EMS concentrations above 40 mM. Space constraints for raising a barley TILLING population are significantly higher compared to Arabidopsis; thus increased M1 and M2 sterility at high EMS concentrations cannot easily be compensated for by increasing the M1 population size. Therefore, building the final bulk of the 'Barke' TILLING population from treatments in the range of 25 - 40 mM EMS was driven by a compromise between effectiveness and practicality.


We developed a TILLING population of 10,279 M2 plants as a resource for functional analysis of genes and for high-throughput gene discovery in an alternate germplasm of the crop model system barley. The screening of six target genes demonstrated a moderate mutation frequency of approximately one mutation per 0.5 Mb. This frequency should allow for the identification of multiple randomly distributed mutations in chosen gene regions. Phenotypic evaluation of portions of the M2 and M3 generations revealed the presence of a wide spectrum of morphological diversity, which is a further indication of the forward genetic potential of this resource.

A detailed TILLING analysis was performed for the homeodomain-leucine-zipper (HD-ZIP) gene HvHox1, which is the major factor controlling the row-type morphology of the barley spike. Multiple alleles causing phenotypic changes of the two-rowed spike morphology of the 'Barke' cultivar were obtained. A comparison of newly identified and previously characterised HvHox1 mutants indicated a bias for preferred nucleotide/AA positions that lead to altered row-type character if affected by mutation. These results contribute to a better understanding of the functionally relevant sites of the HvHOX1 protein. Furthermore, this exemplary case illustrates the specific potential of TILLING to provide multiple independent alleles of a single gene with varying levels of phenotypic expressivity.


Seed material

Seeds of Hordeum vulgare L. cultivar 'Barke' (obtained from: Saatzucht Josef Breun GdbR, Herzogenaurach, Germany) were used for chemical mutagenesis. 'Barke' is a spring-type, two-rowed, European malting variety. 'Barke' was generated by cross-breeding between the two-rowed spring barley varieties 'Alexis' (BSA 1102) and 'Libelle' (1256) in 1996. The pedigrees of both 'Barke' parents are: 'Alexis' ('Br1747' × 'Rupee') × 'Br1622'from 1990, and 'Libelle' ('Br 1622' × 'Trumpf') from 1968 (Saatzucht Josef Breun GdbR, personal communication).

Chemical mutagenesis

Barley seeds were treated at different concentrations with the chemical mutagen ethyl methanesulfonate (EMS) according to David G. Caldwell (personal communication). Batches of ~1,600 seeds were used to fill 2000-ml glass flasks and were pre-soaked in 500 ml of deionised water (ddH2O) for 4 h at room temperature (20-25°C). Then, the water was replaced by 350 ml of an EMS solution [20-60 mM in ddH2O], and seeds were gently shaken (125 rpm, tabletop shaker) for 16 h at room temperature. Subsequently, the EMS solution was collected for decontamination. The seeds were washed two times with 250 ml of 200 mM sodium thiosulphate (30 min for each step) and subsequently two times with 1 l of ddH20 (first for 30 min, second for 1 h). After removing the supernatant, the seeds were transferred to trays covered with Whatman paper and air-dried at 4°C (16 h) prior to sowing.

Creation of the TILLING population

Cultivation of M1

M1 plants were greenhouse-grown until maturity in 299-well seed trays placed on tables that were flooded daily for 12 minutes. After heading, the flooding was reduced to a two-day cycle. Twenty-three and twenty-nine days after sowing, seedlings were sprayed with a fertiliser solution (0.2% Wuxal, nitrogen-phosphate-potassium 8-8-6, AGLUKON GmbH & Co. KG, Düsseldorf, Germany). Subsequently, plants were watered by flooding the tables every second or third day in the summer and winter, respectively, with water-soluble Hakaphos blue (nitrogen-phosphate-potassium 15-10-15 & nutrient salts, COMPO GmbH & Co.KG, Münster, Germany). Every M1 generation was evaluated for germination and sterility rates. Germination was measured as the percentage of developed seedlings, and sterility was measured as seed set per harvested main spike. Germination data were collected 8 days after sowing for EMS treatments between 20 and 40 mM and 12 days after sowing for concentrations between 45 and 60 mM (due to the delayed development of those seedlings). Plants were allowed to self-pollinate, and the main spike of each plant was harvested as a source of M2 seeds.

Cultivation of M2

From each M1, a maximum of ten M2 individuals was cultivated as a family to the tillering stage (BBCH 30, [47]). Subsequently, plants were evaluated, and one or two individuals of each family were selected based on viability to ensure a sufficient seed set for further cultivation. In cases of two selected individuals per family, the siblings were integrated into the population as individual mutant lines.

The M2 plants were cultivated as batches either under greenhouse (1,000 M2 families) or field (3,780 M2 families) conditions. The M2 seedlings were regularly monitored for the presence of chlorophyll defects and other visible mutant phenotypes. For the analysis of chlorophyll mutant frequencies, we calculated the relative frequencies as mutations per 100 germinated M2 plants [48]. All phenotypes were scored in reference to the parent cultivar 'Barke.' The M3 seeds from individual M2 mutants were collected, catalogued, vacuum-packed and stored at 4°C until use as a resource for phenotyping.

Cultivation of M3

The M3 progeny of identified mutants were grown for genotyping and for the determination of visible mutant phenotypes. Forward genetic mutant screening was performed on subsets of M3 families (16 seeds per single M2 individual) derived from EMS concentrations between 20 and 40 mM. These were cultivated under field conditions. Plants were scored for visible phenotypes during the seedling stage, tillering, awn emergence, time of heading and final maturity. All phenotypes were scored in reference to the parent cultivar 'Barke.'

Genomic DNA isolation

Genomic DNA was prepared from young leaves after lyophilisation of tissue. Twenty to thirty milligrams of lyophilised leaf tissue (the remainder of each sample was stored as a backup) was milled in 2 ml tubes, and the DNA was extracted essentially according to a previously described procedure [49, 50]. The DNA was subsequently transferred to 96-well plates and concentrations were measured with an FLX 800 Microplate Fluorescence Reader (384 square well plates; BIO-TEC Instruments, Inc.) using Hoechst 33258 blue dye (Sigma, Deisenhofen, Germany) according to the manufacturer's instructions (Version 04/2004). Aliquots of the DNA were diluted to a final concentration of 20 ng/μl for PCR and arranged in 2570 two-dimensional (2D) eight-fold pools [51] for mutation screening.

Cel I-based mutation screening

Primers, PCR

Oligonucleotides for PCR were designed either directly using the Primer3 software [52] or via the program CODDLE (Codons Optimised to Discover Deleterious Lesions) [53]. Gene-specific primer design was either based on genomic or EST sequence information for the candidate genes. Unlabelled and identical primers labelled at the 5'end with the fluorescent dye IRD700 (forward) or IRD800 (reverse), respectively, were mixed and used for PCR amplification as follows: 3:2 (labelled : unlabeled) ratio for the 100 μM IRD700-labelled forward primer and 3:1 (labelled : unlabeled) ratio for the 100 μM IRD800-labelled reverse primer. Primers were designed to have melting temperatures between 60°C and 70°C. PCR amplification was carried out in a 30-μl volume containing 20 ng of individual pooled DNA, 1× buffer (Qiagen, Hilden, Germany), 0.2 mM dNTPs, 0.3 μM primers, and 0.03 U of Taq polymerase (Qiagen, Hilden, Germany). The PCR reactions were conducted using a thermal cycler (Applied Biosystems 9800 Fast Thermal Cycler, Foster City, USA) as follows: heat denaturation at 95°C for 3 min, followed by 8 cycles of touchdown PCR (94°C for 30 s, annealing at primer Tm +3°C to Tm+4°C for 30 s, decreasing -1°C per cycle, and extension at 72°C for 1 to 1.30 min [for 300- to 1.500-bp products]); 35 cycles of: heat denaturation at 94°C for 30 s, annealing at primer Tm -4°C to Tm -5°C for 30 s, decreasing -1°C per cycle, and extension at 72°C for 1 to 1.30 min; and final extension at 72°C for 10 min. The amplification step was followed by heteroduplex formation: inactivation and denaturation at 99°C for 10 min; and a re-annealing process of 70 cycles for 20 s at 70°C to 49°C, decreasing 0.3°C per cycle. The primers used for PCR amplification of DNA pools and sequencing are listed in Additional file 4.

CEL1 nuclease mismatch cleavage assay

After PCR amplification, samples were incubated for 45 min at 45°C in a 10-μl volume of 50-100 ng of DNA with 0.06 U of Surveyor CEL1 enzyme (Transgenomics, Omaha, USA) in a 10× buffer (10 mM HEPES (pH 7.0), 10 mM KCl, 10 mM MgCl4*7H20, 0.002% Triton X-100 and 10 μg/ml BSA). CEL1 digestion was stopped by adding 5 μl of 75 mM EDTA (pH 8.0) followed by freezing (-20°C) of the probes for 30 min.

Precipitation of 30-μl samples in 96-well PCR plates was carried out by adding 5 μl of 3 M Na-Acetate (pH 5.2) and 75 μl of EtOH (99.8%), subsequent shaking for 15 - 20 min, and centrifugation for 30 min at 4470.6 × g. The supernatant was discarded by turning the plates facedown on top of a stack of filter paper and centrifuging for 1 min at 497 × g. The precipitated samples were washed in 100 μl of ethanol (75%) at room temperature (20 - 25°C) while centrifuging for 30 min at 4,470.6 × g. Removal of the supernatant was performed as described above. Finally, samples were dried for 20 min at room temperature, resuspended in 8 μl of formamide loading buffer (33% (w/v) deionised formamide (Roth, Karlsruhe, Germany), 25 mM Tris (pH 7.5), 25 mM EDTA and ~0.02% (w/v) bromphenol blue) with constant shaking (300 rpm) for 5-10 min, denatured for 5 min at 95°C and subsequently placed on ice until electrophoresis. Electrophoresis of the denatured DNA samples was performed (6.0% Long Ranger polyacrylamide slab gels, FMC Corporation, composition according to manufacturer's data; 1× TBE buffer: 89 mM Tris Base, 89 mM boric acid, 2 mM Na2EDTA*2H20) at 1,500 V, 40 mA and 40 W settings on a Li-COR 4300 DNA Analyser (LI-COR Biosciences, Lincoln, NE). Gel images were visualised using the GelBuddy software [54, 55] and inspected visually for the presence of cleavage products.

TILLING screens were conducted by using a CEL1-based heteroduplex analysis for mutation detection. The CEL1-based assay relies on using the LI-COR two laser/PA-slabgel system. The optimal mutation detection procedure and pool size for 'Barke' DNA was determined by performing test screenings using known mutations in the Hv-eIF4E gene [56] of the barley genotypes 'IGRI' and 'Franka.' Thereby, we found that for pooled DNA, a 1:10 dilution of CEL1 in a specific buffer solution combined with a 45 min digestion at 45°C provides an optimal ratio between the background and signal strength. The 8-fold arrangement was found to be most reliable, although mutations were still detectable in the tested 16-fold pools. To improve the mutation screening, the population was screened in a two-dimensional format. In each 96-well pool plate, the DNA of a total of 384 M2 individuals was arranged in row and column 8-fold pools. As each sample is present in two different pools, any mutation should be detected twice.

DNA sequencing and sequence data processing

For confirmation of presumed mutated loci, amplicons of the respective target gene were generated from putative mutants by utilising the same PCR conditions as established for CEL I analysis. Amplicons were purified by ultrafiltration using a QIAvac 96 vacuum (Qiagen, Hilden, Germany), processed with a NucleoFast-96 PCR Plate (MACHEREY-NAGEL GmbH & Co. KG, Düren, Germany) and directly cycle-sequenced using the ABI Big Dye terminator v3.1 sequencing standard kit according to the manufacturer's protocol (Applied Biosystems, Foster City, USA). Sequencing reactions were resolved on a 96-capillary sequencing device (3730xl DNA Analyzer, Applied Biosystems, Foster City, USA). Sequences were aligned against the "wild-type" reference of the 'Barke' cultivar with Sequencher 4.6 software (Gene Codes, Ann Arbor, MI).

Genes were routinely analysed using the program CODDLE to obtain a gene model and to identify the region that would have the highest likelihood of being functionally affected by EMS mutagenesis. The software PARSESNP (Project Aligned Related Sequences and Evaluate SNPs) [57, 58] was used to predict the severity of each identified mutation. If necessary, the Web version of the FGENESH program was used for gene structure prediction, and parameters were set for monocot plants [59]. The Conserved Domain Database (CDD) [60] accessible on the National Center for Biotechnology Information (NCBI) was used for detailed analysis of identified mutations, as it provides multiple sequence alignments and derived database search models that suggest protein domains conserved in molecular evolution [61, 62].

Statistical Analysis

Statistical analysis was conducted to investigate possible EMS dosage effects on genomic mutation frequencies, M1 germination and sterility, M2 chlorophyll mutant occurrence, and M3 seed yield. The STATISTICA data analysis software version 6 (StatSoft, Inc. (2001), was used to conduct correlation analyses [(Pearson (r) and Spearman rank (rs)], tests of linear relationships [63] and one-factorial ANOVAs with the respective software protocols.

How to access the TILLING population

The 'Barke' TILLING population was developed in the framework of the German GABI (G enome A nalysis in the Bi ological system plant) program. Screening of the resource is not established as a service yet but can be performed on the basis of collaboration. For further information please contact Nils Stein, IPK Gatersleben.


  1. 1.

    Hagberg A, Persson G: Practical use of mutations in genetics, taxonomy and breeding. Barley genetics I. Proceedings First International Barley Genetics Symposium 26-31 August 1963. Edited by: Broekhuizen S. 1964, Wageningen, the Netherlands, 55-67.

    Google Scholar 

  2. 2.

    Muller HJ: Artificial transmutation of the gene. Science. 1927, 66: 84-87. 10.1126/science.66.1699.84.

    Article  CAS  PubMed  Google Scholar 

  3. 3.

    Stadler LJ: Mutations in barley induced by X-rays and radium. Science. 1928, 68: 186-187. 10.1126/science.68.1756.186.

    Article  CAS  PubMed  Google Scholar 

  4. 4.

    Harten van AM: Mutation breeding: theory and practical applications. 1998, Cambridge: Cambridge University Press, 101-

    Google Scholar 

  5. 5.

    Freisleben R, Lein A: Über die Auffindung einer mehltauresistenten Mutante nach Röntgenbestrahlung einer anfälligen reinen Linie von Sommergerste. Naturwiss. 1942, 30: 608-10.1007/BF01488231.

    Article  Google Scholar 

  6. 6.

    Freisleben R, Lein A: Vorarbeiten zur züchterischen Auswertung röntgeninduzierter Mutationen I. Z Pflanzenz. 1943, 25: 235-254.

    Google Scholar 

  7. 7.

    Freisleben R, Lein A: Vorarbeiten zur züchterischen Auswertung röntgeninduzierter Mutationen II. Z Pflanzenz. 1943, 25: 255-283.

    Google Scholar 

  8. 8.

    Stubbe H: Chemische Mutagenese. Proceedings Erwin-Baur-Gedächtnisvorlesungen I; Gatersleben. Edited by: Stubbe H. 1960, Berlin: Akademie-Verlag

    Google Scholar 

  9. 9.

    FAO/IAEA Mutant Variety Database. []

  10. 10.

    Micke A, Donini B, Maluszynski M: Induced Mutations for Crop Improvement. Mutat Breed Rev. 1990, 7: 1-41.

    Google Scholar 

  11. 11.

    McCallum CM, Comai L, Greene EA, Henikoff S: Targeted screening for induced mutations. Nature Biotechnol. 2000, 18: 455-457. 10.1038/74542.

    Article  CAS  Google Scholar 

  12. 12.

    Greene EA, Codomo CA, Taylor NE, Henikoff JG, Till BJ, Reynolds SH, Enns LC, Burtner C, Johnson JE, Odden AR, Comai L, Henikoff S: Spectrum of chemically induced mutations from a large-scale reverse-genetic screen in Arabidopsis. Genetics. 2003, 164: 731-740.

    PubMed Central  CAS  PubMed  Google Scholar 

  13. 13.

    Slade AJ, Fuerstenberg SI, Loeffler D, Steine MN, Facciotti D: A reverse genetic, nontransgenic approach to wheat crop imporvement by TILLING. Nat Biotech. 2005, 23: 75-81. 10.1038/nbt1043.

    Article  CAS  Google Scholar 

  14. 14.

    Slade AJ, Knauf VC: TILLING without a plough: a new method with applications for reverse genetics. Transgenic Res. 2005, 14: 109-15. 10.1007/s11248-005-2770-x.

    Article  CAS  PubMed  Google Scholar 

  15. 15.

    Komatsuda T, Pourkheirandish M, He C, Azhaguvel O, Kanamori H, Perovic D, Stein N, Graner A, Wicker T, Tagiri A, Lundqvist U, Fujimura T, Matsuoka M, Matsumoto T, Yano M: Six-rowed barley originated from a mutation in a homeodomain-leucine zipper I-class homeobox gene. Proc Natl Acad Sci USA. 2007, 104: 1424-1429. 10.1073/pnas.0608580104.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  16. 16.

    Caldwell DG, McCallum N, Shaw P, Muehlbauer GJ, Marshall DF, Waugh R: A structured mutant population for forward and reverse genetics in barley (Hordeum vulgare L.). Plant J. 2004, 40: 143-50. 10.1111/j.1365-313X.2004.02190.x.

    Article  CAS  PubMed  Google Scholar 

  17. 17.

    Zimmer KG: Studies on quantitative radiation biology. 1961, Oliver & Boyd: Edinburgh-London

    Google Scholar 

  18. 18.

    Sharma B, Bansal HC: The multiple-cluster Mutation complex in mutagenesis with higher plants. Theor Appl Genet. 1972, 42: 25-31. 10.1007/BF00306074.

    Article  CAS  PubMed  Google Scholar 

  19. 19.

    Sheeba A, Ibrahim SM, Yogameenakshi P: Induced chlorophyll mutation studies in sesame. Sesame and Safflower Newsletter. 2003, 18: 33-38.

    Google Scholar 

  20. 20.

    Gustafsson Å: The mutation system of the chlorophyll apparatus. Lunds Universitets Årsskrift NF Avd 2. 1940, 36: 1-40.

    Google Scholar 

  21. 21.

    von Wettstein D: Multiple allelism in induced chlorophyll mutants of barley. Hereditas. 1960, 47: 767-770.

    Google Scholar 

  22. 22.

    Nybom N: The pigment characteristics of chlorophyll mutations in barley. Hereditas. 1955, 41: 483-498.

    Article  Google Scholar 

  23. 23.

    Talamè V, Bovina R, Sanguineti MC, Tuberosa R, Lundqvist U, Salvi S: TILLMore, a resource for the discovery of chemically induced mutants in barley. Plant Biotechnol J. 2008, 6: 477-485. 10.1111/j.1467-7652.2008.00341.x.

    Article  PubMed  Google Scholar 

  24. 24.

    Chuang C-F, Meyerowitz EM: Specific and heritable genetic interference by double-stranded RNA in Arabidopsis thaliana. Proc Natl Acad Sci USA. 2000, 97: 4985-4990. 10.1073/pnas.060034297.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  25. 25.

    Alonso JM, Stepanova AN, Leisse TJ, Kim CJ, Chen H, Shinn P, Stevenson DK, Zimmerman J, Barajas P, Cheuk R, et al: Genome-wide insertional mutagenesis of Arabidopsis thaliana. Science. 2003, 301: 653-657. 10.1126/science.1086391.

    Article  PubMed  Google Scholar 

  26. 26.

    Gilchrist EJ, Haughn GW: TILLING without a plough: a new method with applications for reverse genetics. Curr Opin Plant Biol. 2005, 8: 211-215. 10.1016/j.pbi.2005.01.004.

    Article  CAS  PubMed  Google Scholar 

  27. 27.

    Henriksson E, Olsson ASB, Johannesson H, Johansson H, Hanson J, Engström P, Söderman E: Homeodomain leucine zipper class I genes in Arabidopsis. Expression patterns and phylogenetic relationships. Plant Physiol. 2005, 139: 509-518. 10.1104/pp.105.063461.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  28. 28.

    Affolter M, Percival-Smith A, Muller M, Leupint W, Gehring WJ: DNA binding properties of the purified Antennapedia homeodomain (helix-turn-helix). Proc Natl Acad Sci USA. 1990, 87: 4093-4097. 10.1073/pnas.87.11.4093.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  29. 29.

    Vershon AK: Protein interactions of homeodomain proteins. Curr Opin Biotechnol. 1996, 7: 392-396. 10.1016/S0958-1669(96)80113-3.

    Article  CAS  PubMed  Google Scholar 

  30. 30.

    Vainius D: Investigation of interactions between homeodomain proteins and DNA. PhD thesis. 2004, Georg-August-Universität Göttingen, Fakultät für Physik

    Google Scholar 

  31. 31.

    Wolberger C: Homeodomain interactions. Curr Opin Struct Biol. 1996, 6: 62-68. 10.1016/S0959-440X(96)80096-0.

    Article  CAS  PubMed  Google Scholar 

  32. 32.

    Gehring WJ, Affolter M, Büirglin T: Homeodomain proteins. Annu Rev Biochem. 1994, 63: 487-526. 10.1146/

    Article  CAS  PubMed  Google Scholar 

  33. 33.

    Ayoub M, Symons S, Edney M, Mather D: QTLs affecting kernel size and shape in a two-rowed by six-rowed barley cross. Theor Appl Genet. 2002, 105: 237-247. 10.1007/s00122-002-0941-1.

    Article  CAS  PubMed  Google Scholar 

  34. 34.

    Komatsuda T, Mano Y: Molecular mapping of the intermedium spike-c (int-c) and non-brittle rachis 1 (btr1) loci in barley (Hordeum vulgare L.). Theor Appl Genet. 2002, 105: 85-90. 10.1007/s00122-001-0858-0.

    Article  CAS  PubMed  Google Scholar 

  35. 35.

    Lundqvist U, Lundqvist A: The co-operation between intermedium genes and the six-row gene hex-v in a six-row variety of barley. Hereditas. 1989, 110: 227-233. 10.1111/j.1601-5223.1989.tb00785.x.

    Article  Google Scholar 

  36. 36.

    Arumuganathan K, Earle ED: Estimation of nuclear DNA content of plants by flow cytometry. Plant Mol Biol Rep. 1991, 9: 229-233. 10.1007/BF02672073.

    Article  CAS  Google Scholar 

  37. 37.

    Stein N: Triticeae genomics: advances in sequence analysis of large genome cereal crops. Chromosome Research. 2007, 15: 21-31. 10.1007/s10577-006-1107-9.

    Article  CAS  PubMed  Google Scholar 

  38. 38.

    Miklos GLG, Rubin GM: The role of the genome project in determining gene function: Insights from model organisms. Cell. 1996, 86: 521-529. 10.1016/S0092-8674(00)80126-9.

    Article  CAS  PubMed  Google Scholar 

  39. 39.

    Wong GK-S, Wang J, Tao L, Tan J, Zhang JG, Passey DA, Yu J: Compositional Gradients in Gramineae Genes. Genome Res. 2002, 12: 851-856. 10.1101/gr.189102.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  40. 40.

    Barakat A, Matassi G, Bernardi G: Distribution of genes in the genome of Arabidopsis thaliana and its implications for the genome organization of plants. Proc Natl Acad Sci USA. 1998, 95: 10044-10049. 10.1073/pnas.95.17.10044.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  41. 41.

    Hänsel H: Model for a theoretical estimate of optimal mutation rates per M1-nucleus with a view to selecting beneficial mutations in different M-generation. Induced mutations and their utilization. 1967, Berlin, Akademic-Verlag, 79-87.

    Google Scholar 

  42. 42.

    Malysheva-Otto L, Ganal M, Roder M: Analysis of molecular diversity, population structure and linkage disequilibrium in a worldwide survey of cultivated barley germplasm (Hordeum vulgare L.). BMC Genet. 2006, 7: 6-10.1186/1471-2156-7-6.

    PubMed Central  Article  PubMed  Google Scholar 

  43. 43.

    Rostoks N, Ramsay L, MacKenzie K, Cardle L, Bhat PR, Roose ML, Svensson JT, Stein N, Varshney RK, Marshall DF, Graner A, Close TJ, Waugh R: Recent history of artificial outcrossing facilitates whole-genome association mapping in elite inbred crop varieties. Proc Natl Acad Sci. 2006, 103: 18656-18661. 10.1073/pnas.0606133103.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  44. 44.

    Kunzel G, Korzun L, Meister A: Cytologically integrated physical restriction fragment length polymorphism maps for the barley genome based on translocation breakpoints. Genetics. 2000, 154: 397-412.

    PubMed Central  CAS  PubMed  Google Scholar 

  45. 45.

    Colbert T, Till BJ, Tompa R, Reynolds S, Steine MN, Yeung AT, McCallum CM, Comai L, Henikoff S: High-throughput screening for induced point mutations. Plant Physiol. 2001, 126: 480-484. 10.1104/pp.126.2.480.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  46. 46.

    Konzak CF, Nilan RA, Wagner J, Foster J: Efficient chemical mutagenesis. The use of induced mutations in plant breeding. Report of FAO/IAEA Technical Meeting; Rome, Italy. 1965, New York: Pergamon Press, 49-70.

    Google Scholar 

  47. 47.

    Meier U, Bleiholder H: BBCH Skala. Phänologische Entwicklungsstadien wichtiger landwirtschaftlicher Kulturen, einschließblich Blattgemüse und Unkräuter. 2006, Bergen/Dumme: Agrimedia GmbH, 70-

    Google Scholar 

  48. 48.

    Gaul H: Critical analysis of the methods for determining the mutation frequency after seed treatment with mutagens. Genetica Agraria. 1960, 12: 297-318.

    Google Scholar 

  49. 49.

    Doyle JJ, Doyle JL: A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochem Bull. 1987, 19: 11-15.

    Google Scholar 

  50. 50.

    Doyle JJ, Doyle JL: Isolation of plant DNA from fresh tissue. Focus. 1990, 12: 13-15.

    Google Scholar 

  51. 51.

    Till BJ, Zerr T, Comai L, Henikoff S: A protocol for TILLING and Ecotilling in plants and animals. Nature Protocols. 2006, 1: 2465-2477. 10.1038/nprot.2006.329.

    Article  CAS  PubMed  Google Scholar 

  52. 52.

    Primer3 software. []

  53. 53.

    CODDLE (Codons Optimized to Discover Deleterious Lesions). []

  54. 54.

    GelBuddy software. []

  55. 55.

    Zerr T, Henikoff S: Automated band mapping in electrophoretic gel images using background information. Nucleic Acids Res. 2005, 33: 2806-2812. 10.1093/nar/gki580.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  56. 56.

    Stein N, Perovic D, Kumlehn J, Pellio B, Stracke S, Streng S, Ordon F, Graner A: The eukaryotic translation initiation factor 4E confers multiallelic recessive Bymovirus resistance in Hordeum vulgare (L.). Plant J. 2005, 42: 912-922. 10.1111/j.1365-313X.2005.02424.x.

    Article  CAS  PubMed  Google Scholar 

  57. 57.

    PARSESNP (Project Aligned Related Sequences and Evaluate SNPs). []

  58. 58.

    Taylor NE, Greene EA: PARSESNP: A tool for the analysis of nucleotide polymorphisms. Nucleic Acids Res. 2003, 31: 3808-3811. 10.1093/nar/gkg574.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  59. 59.

    Salamov A, Solovyev V: Ab initio gene finding in Drosophila genomic DNA. Genome Res. 2000, 10: 516-522. 10.1101/gr.10.4.516.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  60. 60.

    NCBI's Conserved Domain Database (CDD). []

  61. 61.

    Marchler-Bauer A, Anderson JB, Chitsaz F, Derbyshire MK, Deweese-Scott C, Fong JH, Geer LY, Geer RC, Gonzales NR, Gwadz M, He S, Hurwitz DI, Jackson JD, Ke Z, Lanczycki CJ, Liebert CA, Liu C, Lu F, Lu S, Marchler GH, Mullokandov M, Song JS, Tasneem A, Thanki N, Yamashita RA, Zhang D, Zhang N, Bryant SH: CDD: specific functional annotation with the Conserved Domain Database. Nucleic Acids Res. 2009, 37: D205-D210. 10.1093/nar/gkn845.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  62. 62.

    Marchler-Bauer A, Anderson JB, Derbyshire MK, DeWeese-Scott C, Gonzales NR, Gwadz M, Hao L, He S, Hurwitz DI, Jackson JD, Ke Z, Krylov D, Lanczycki CJ, Liebert CA, Liu C, Lu F, Lu S, Marchler GH, Mullokandov M, Song JS, Thanki N, Yamashita RA, Yin JJ, Zhang D, Bryant SH: CDD: a conserved domain database for interactive domain family analysis. Nucleic Acids Res. 2007, 35: 237-40. 10.1093/nar/gkl951.

    Article  Google Scholar 

  63. 63.

    Linder A, Berchthold W: Statistische Auswertung von Prozentzahlen. 2007, Basel: UTB Birkhäuser Verlag, 35: 237-40.

    Google Scholar 

  64. 64.

    Griffiths S, Dunford RP, Coupland G, Laurie DA: The evolution of CONSTANS-like gene families in barley, rice, and Arabidopsis. Plant Physiol. 2003, 131: 1855-1867. 10.1104/pp.102.016188.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  65. 65.

    Bueschges R, Hollricher K, Panstruga R, Simons G, Wolter M, Frijters A, van Daelen R, Lee van der T, Diergaarde P, Groenendijk J, Toepsch S, Vos P, Salamini F, SchulzeLefert P: The barley Mlo gene: a novel control element of plant pathogen resistance. Cell. 1997, 88: 695-705. 10.1016/S0092-8674(00)81912-1.

    Article  Google Scholar 

  66. 66.

    Panstruga R, Buschges R, Piffanelli P, Schulze-Lefert P: A contiguous 60 kb genomic stretch from barley reveals molecular evidence for gene islands in a monocot genome. Nucleic Acids Res. 1998, 26: 1056-1062. 10.1093/nar/26.4.1056.

    PubMed Central  Article  CAS  PubMed  Google Scholar 

  67. 67.

    Henikoff JG, Pietrokovski S, McCallum CM, Henikoff S: Blocks-based methods for detecting protein homology. Electrophoresis. 2000, 21: 1700-1706. 10.1002/(SICI)1522-2683(20000501)21:9<1700::AID-ELPS1700>3.0.CO;2-V.

    Article  CAS  PubMed  Google Scholar 

Download references


The authors thank Mary Ziems, Bettina Brückner and Mechthild Pürschel for their excellent technical assistance, and Peter Schreiber, Kathrin Thiebe and Katrin Menzel for substantial help and support in field and greenhouse work. Financial support of the project was provided by the Federal Ministry of Education and Research (BMBF) within the GABI program (FKZ: 0313123 C). Thanks to KWS Lochow AG for helping to establish high-density cultivation and to Saatzucht Josef Breun GdbR for providing stocks of 'Barke' founder seeds. We thank David Caldwell, David Leader, and Robbie Waugh for helpful instruction during the establishment phase.

Author information



Corresponding author

Correspondence to Nils Stein.

Additional information

Authors' contributions

SG established the barley TILLING platform. PB and NS conceived the study. NS contributed to the design of the study. SG and NS wrote the manuscript. TK provided information on the HvHox1 gene, including the sequence and phenotypic data of additional barley mutants. UL was involved in the forward genetic screening of the mutant populations. All authors read and approved the manuscript.

Electronic supplementary material

Additional figure 1 - Examples of chlorophyll-deficient seedling mutants in M

Additional file 1: 2 . Supplemental figure showing chlorophyll-deficient seedling mutants. Mutants belonging to the unicolour sub-class: (A) albina, (B) xantha, and (C) viridis. Mutants belonging to the bicolour sub-class: (D) viridoalbina and (E) striata. The classification corresponds to [20]. (PDF 902 KB)

Additional figure 2 - Frequency of mutant phenotype categories observed in M

Additional file 2: 2 and M 3 screenings. This supplemental figure is an overview of the occurrence of different mutant phenotypes observed in the mutant populations. Mutant phenotypes observed in M2 and M3 derived from mutagenesis at 20 - 35 mM EMS were classified into six morphological categories. (A) In the M2 generation, approximately 20% of the 12,703 individuals showed a mutant phenotype vs. the 'Barke' wild-type. (B) In the M3 generation, 1,200 M3 families (16 plants per M3 line) were grown, and approximately 37% of the total M3 families displayed a mutant phenotype vs. the 'Barke' wild-type. (PDF 72 KB)

Additional table 1 - A series of mutations in the

Additional file 3: HvHox1 gene identified by TILLING. The additional table shows the position and mutation types identified in the HvHox1 gene. (DOC 62 KB)

Additional file 4: Additional table 2 - Sequences of primers used for TILLING of candidate genes. This additional table specifies primers used for PCR amplification of DNA pools and sequencing. (DOC 25 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and Permissions

About this article

Cite this article

Gottwald, S., Bauer, P., Komatsuda, T. et al. TILLING in the two-rowed barley cultivar 'Barke' reveals preferred sites of functional diversity in the gene HvHox1. BMC Res Notes 2, 258 (2009).

Download citation


  • Barley Cultivar
  • Ethyl Methanesulphonate
  • High Mutation Frequency
  • Tilling Population
  • Lateral Spikelet