Evolution of the Grain Dispersal System in Barley.

About 12,000 years ago in the Near East, humans began the transition from hunter-gathering to agriculture-based societies. Barley was a founder crop in this process, and the most important steps in its domestication were mutations in two adjacent, dominant, and complementary genes, through which grains were retained on the inflorescence at maturity, enabling effective harvesting. Independent recessive mutations in each of these genes caused cell wall thickening in a highly specific grain "disarticulation zone," converting the brittle floral axis (the rachis) of the wild-type into a tough, non-brittle form that promoted grain retention. By tracing the evolutionary history of allelic variation in both genes, we conclude that spatially and temporally independent selections of germplasm with a non-brittle rachis were made during the domestication of barley by farmers in the southern and northern regions of the Levant, actions that made a major contribution to the emergence of early agrarian societies.

Wheat Ym2 originated from Aegilops sharonensis and confers resistance to soil-borne Wheat yellow mosaic virus infection to the roots.

Wheat yellow mosaic virus (WYMV) is a pathogen transmitted into its host's roots by the soil-borne vector Polymyxa graminis. Ym1 and Ym2 genes protect the host from the significant yield losses caused by the virus, but the mechanistic basis of these resistance genes remains poorly understood. Here, it has been shown that Ym1 and Ym2 act within the root either by hindering the initial movement of WYMV from the vector into the root and/or by suppressing viral multiplication. A mechanical inoculation experiment on the leaf revealed that the presence of Ym1 reduced viral infection incidence, rather than viral titer, while that of Ym2 was ineffective in the leaf. To understand the basis of the root specificity of the Ym2 product, the gene was isolated from bread wheat using a positional cloning approach. The candidate gene encodes a CC-NBS-LRR protein and it correlated allelic variation with respect to its sequence with the host's disease response. Ym2 (B37500) and its paralog (B35800) are found in the near-relatives, respectively, Aegilops sharonensis and Aegilops speltoides (a close relative of the donor of bread wheat's B genome), while both sequences, in a concatenated state, are present in several accessions of the latter species. Structural diversity in Ym2 has been generated via translocation and recombination between the two genes and enhanced by the formation of a chimeric gene resulting from an intralocus recombination event. The analysis has revealed how the Ym2 region has evolved during the polyploidization events leading to the creation of cultivated wheat.

Heterosis in crop improvement.

Heterosis, also known as hybrid vigor, is the phenomenon wherein a progeny exhibits superior traits relative to one or both parents. In terms of crop breeding, this usually refers to the yield advantage of F1 hybrids over both inbred parents. The development of high-yielding hybrid cultivars across a wider range of crops is key to meeting future food demands. However, conventional hybrid breeding strategies are proving to be exceptionally challenging to apply commercially in many self-pollinating crops, particularly wheat and barley. Currently in these crops, the relative performance advantage of hybrids over inbred line cultivars does not outweigh the cost of hybrid seed production. Here, we review the genetic basis of heterosis, discuss the challenges in hybrid breeding, and propose a strategy to recruit multiple heterosis-associated genes to develop lines with improved agronomic characteristics. This strategy leverages modern genetic engineering tools to synthesize supergenes by fusing multiple heterotic alleles across multiple heterosis-associated loci. We outline a plan to assess the feasibility of this approach to improve line performance using barley (Hordeum vulgare) as the model self-pollinating crop species, and a few heterosis-associated genes. The proposed method can be applied to all crops for which heterotic gene combinations can be identified.

Back to the future for drought tolerance.

Global agriculture faces increasing pressure to produce more food with fewer resources. Drought, exacerbated by climate change, is a major agricultural constraint costing the industry an estimated US$80 billion per year in lost production. Wild relatives of domesticated crops, including wheat (Triticum spp.) and barley (Hordeum vulgare L.), are an underutilized source of drought tolerance genes. However, managing their undesirable characteristics, assessing drought responses, and selecting lines with heritable traits remains a significant challenge. Here, we propose a novel strategy of using multi-trait selection criteria based on high-throughput spectral images to facilitate the assessment and selection challenge. The importance of measuring plant capacity for sustained carbon fixation under drought stress is explored, and an image-based transpiration efficiency (iTE) index obtained via a combination of hyperspectral and thermal imaging, is proposed. Incorporating iTE along with other drought-related variables in selection criteria will allow the identification of accessions with diverse tolerance mechanisms. A comprehensive approach that merges high-throughput phenotyping and de novo domestication is proposed for developing drought-tolerant prebreeding material and providing breeders with access to gene pools containing unexplored drought tolerance mechanisms.

Guard Cell Transcriptome Reveals Membrane Transport, Stomatal Development and Cell Wall Modifications as Key Traits Involved in Salinity Tolerance in Halophytic Chenopodium quinoa.

A comparative investigation was conducted to evaluate transcriptional changes in guard cells (GCs) of closely related halophytic (Chenopodium quinoa) and glycophytic (Spinacia oleracea) species. Plants were exposed to 3 weeks of 250 mM sodium chloride treatment, and GC-enriched epidermal fragments were mechanically prepared. In both species, salt-responsive genes were mainly related to categories of protein metabolism, secondary metabolites, signal transduction and transport systems. Genes related to abscisic acid (ABA) signaling and ABA biosynthesis were strongly induced in quinoa but not in spinach GCs. Also, expression of the genes encoding transporters of amino acids, proline, sugars, sucrose and potassium increased in quinoa GCs under salinity stress. Analysis of cell-wall-related genes suggests that genes involved in lignin synthesis (e.g. lignin biosynthesis LACCASE 4) were highly upregulated by salt in spinach GCs. In contrast, transcripts related to cell wall plasticity Pectin methylesterase3 (PME3) were highly induced in quinoa. Faster stomatal response to light and dark measured by observing kinetics of changes in stomatal conductance in quinoa might be associated with higher plasticity of the cell wall regulated by PME3 Furthermore, genes involved in the inhibition of stomatal development and differentiation were highly expressed by salt in quinoa, but not in spinach. These changes correlated with reduced stomatal density and index in quinoa, thus improving its water use efficiency. The fine modulation of transporters, cell wall modification and controlling stomatal development in GCs of quinoa may have resulted in high K+/Na+ ratio, lower stomatal conductance and higher stomatal speed for better adaptation to salinity stress in quinoa.

Grain Disarticulation in Wild Wheat and Barley.

Our industrial-scale crop monocultures, which are necessary to provide grain for large-scale food and feed production, are highly vulnerable to biotic and abiotic stresses. Crop wild relatives have adapted to harsh environmental conditions over millennia; thus, they are an important source of genetic variation and crop diversification. Despite several examples where significant yield increases have been achieved through the introgression of genomic regions from wild relatives, more detailed understanding of the differences between wild and cultivated species for favorable and unfavorable traits is still required to harness these valuable resources. Recently, as an alternative to the introgression of beneficial alleles from the wild into domesticated species, a radical suggestion is to domesticate wild relatives to generate new crops. A first and critical step for the domestication of cereal wild relatives would be to prevent grain disarticulation from the inflorescence at maturity. Discovering the molecular mechanisms and understanding the network of interactions behind grain retention/disarticulation would enable the implementation of approaches to select for this character in targeted species. Brittle rachis 1 and Brittle rachis 2 are major genes responsible for grain disarticulation in the wild progenitors of wheat and barley that were the target of mutations during domestication. These two genes are only found in the Triticeae tribe and are hypothesized to have evolved by a duplication followed by neo-functionalization. Current knowledge gaps include the molecular mechanisms controlling grain retention in cereals and the genomic consequences of strong selection for this essential character.

Grain dispersal mechanism in cereals arose from a genome duplication followed by changes in spatial expression of genes involved in pollen development.

KEY MESSAGE: Grain disarticulation in wild progenitor of wheat and barley evolved through a local duplication event followed by neo-functionalization resulting from changes in location of gene expression. One of the most critical events in the process of cereal domestication was the loss of the natural mode of grain dispersal. Grain dispersal in barley is controlled by two major genes, Btr1 and Btr2, which affect the thickness of cell walls around the disarticulation zone. The barley genome also encodes Btr1-like and Btr2-like genes, which have been shown to be the ancestral copies. While Btr and Btr-like genes are non-redundant, the biological function of Btr-like genes is unknown. We explored the potential biological role of the Btr-like genes by surveying their expression profile across 212 publicly available transcriptome datasets representing diverse organs, developmental stages and stress conditions. We found that Btr1-like and Btr2-like are expressed exclusively in immature anther samples throughout Prophase I of meiosis within the meiocyte. The similar and restricted expression profile of these two genes suggests they are involved in a common biological function. Further analysis revealed 141 genes co-expressed with Btr1-like and 122 genes co-expressed with Btr2-like, with 105 genes in common, supporting Btr-like genes involvement in a shared molecular pathway. We hypothesize that the Btr-like genes play a crucial role in pollen development by facilitating the formation of the callose wall around the meiocyte or in the secretion of callase by the tapetum. Our data suggest that Btr genes retained an ancestral function in cell wall modification and gained a new role in grain dispersal due to changes in their spatial expression becoming spike specific after gene duplication.

Molecular genetics of leaf rust resistance in wheat and barley.

The demand for cereal grains as a main source of energy continues to increase due to the rapid increase in world population. The leaf rust diseases of cereals cause significant yield losses, posing challenges for global food security. The deployment of resistance genes has long been considered as the most effective and sustainable way to control cereal leaf rust diseases. While genetic resistance has reduced the impact of these diseases in agriculture, losses still occur due to the ability of the respective rust pathogens to change and render resistance genes ineffective plus the slow pace at which resistance genes are discovered and characterized. This article highlights novel recently developed strategies based on advances in genome sequencing that have accelerated gene isolation by overcoming the complexity of cereal genomes. The leaf rust resistance genes cloned so far from wheat and barley belong to various protein families, including nucleotide binding site/leucine-rich repeat receptors and transporters. We review recent studies that are beginning to reveal the defense mechanisms conferred by the leaf rust resistance genes identified to date in cereals and their roles in either pattern-triggered immunity or effector-triggered immunity.

Elucidation of the origin of 'agriocrithon' based on domestication genes questions the hypothesis that Tibet is one of the centers of barley domestication.

Wild barley forms a two-rowed spike with a brittle rachis whereas domesticated barley has two- or six-rowed spikes with a tough rachis. Like domesticated barley, 'agriocrithon' forms a six-rowed spike; however, the spike is brittle as in wild barley, which makes the origin of agriocrithon obscure. Haplotype analysis of the Six-rowed spike 1 (vrs1) and Non-brittle rachis 1 (btr1) and 2 (btr2) genes was conducted to infer the origin of agriocrithon barley. Some agriocrithon barley accessions (eu-agriocrithon) carried Btr1 and Btr2 haplotypes that are not found in any cultivars, implying that they are directly derived from wild barley through a mutation at the vrs1 locus. Other agriocrithon barley accessions (pseudo-agriocrithon) carried Btr1 or Btr2 from cultivated barley, thus implying that they originated from hybridization between six-rowed landraces carrying btr1Btr2 and Btr1btr2 genotypes followed by recombination to produce Btr1Btr2. All materials we collected from Tibet belong to pseudo-agriocrithon and thus do not support the Tibetan Plateau as being a center of barley domestication. Tracing the evolutionary history of these allelic variants revealed that eu-agriocrithon represents six-rowed barley lineages that were selected by early farmers, once in south-eastern Turkmenistan (vrs1.a1) and again in the eastern part of Uzbekistan (vrs1.a4).

miR172 downregulates the translation of cleistogamy 1 in barley.

Background and Aims: Floret opening in barley is induced by the swelling of the lodicule, a trait under the control of the cleistogamy1 (cly1) gene. The product of cly1 is a member of the APETALA2 (AP2) transcription factor family, which inhibits lodicule development. A sequence polymorphism at the miR172 target site within cly1 has been associated with variation in lodicule development and hence with the cleistogamous phenotype. It was unclear whether miR172 actually functions in cly1 regulation and, if it does, which miR172 gene contributes to cleistogamy. It was also interesting to explore whether miR172-mediated cly1 regulation occurs at transcriptional level or at translational level. Methods: Deep sequencing of small RNA identified the miR172 sequences expressed in barley immature spikes. miR172 genes were confirmed by computational and expression analysis. miR172 and cly1 expression profiles were determined by in situ hybridization and quantitative expression analysis. Immunoblot analysis provided the CLY1 protein quantifications. Definitive evidence of the role of miR172 in cleistogamy was provided by a transposon Ds-induced mutant of Hv-miR172a. Key Results: A small RNA analysis of the immature barley spike revealed three isomers, miR172a, b and c, of which miR172a was the most abundant. In situ hybridization analysis showed that miR172 and cly1 co-localize in the lodicule primordium, suggesting that these two molecules potentially interact with one another. Immunoblot analysis showed that the sequence polymorphism at the miR172 target site within cly1 reduced the abundance of the CLY1 protein, but not that of its transcript. In a Ds-induced mutant of Hv-miR172a, which generates no mature miR172a, the lodicules fail to grow, resulting in a very small lodicule. Conclusions: Direct evidence is presented to show that miR172a acts to reduce the abundance of the CLY1 protein, which enables open flowering in barley.

Six-rowed spike4 (Vrs4) controls spikelet determinacy and row-type in barley.

Inflorescence architecture of barley (Hordeum vulgare L.) is common among the Triticeae species, which bear one to three single-flowered spikelets at each rachis internode. Triple spikelet meristem is one of the unique features of barley spikes, in which three spikelets (one central and two lateral spikelets) are produced at each rachis internode. Fertility of the lateral spikelets at triple spikelet meristem gives row-type identity to barley spikes. Six-rowed spikes show fertile lateral spikelets and produce increased grain yield per spike, compared with two-rowed spikes with sterile lateral spikelets. Thus, far, two loci governing the row-type phenotype were isolated in barley that include Six-rowed spike1 (Vrs1) and Intermedium-C. In the present study, we isolated Six-rowed spike4 (Vrs4), a barley ortholog of the maize (Zea mays L.) inflorescence architecture gene RAMOSA2 (RA2). Eighteen coding mutations in barley RA2 (HvRA2) were specifically associated with lateral spikelet fertility and loss of spikelet determinacy. Expression analyses through mRNA in situ hybridization and microarray showed that Vrs4 (HvRA2) controls the row-type pathway through Vrs1 (HvHox1), a negative regulator of lateral spikelet fertility in barley. Moreover, Vrs4 may also regulate transcripts of barley SISTER OF RAMOSA3 (HvSRA), a putative trehalose-6-phosphate phosphatase involved in trehalose-6-phosphate homeostasis implicated to control spikelet determinacy. Our expression data illustrated that, although RA2 is conserved among different grass species, its down-stream target genes appear to be modified in barley and possibly other species of tribe Triticeae.

High-resolution genetic mapping and physical map construction for the fertility restorer Rfm1 locus in barley.

KEY MESSAGE: High-resolution genetic linkage mapping and BAC physical mapping narrowed the fertility restorer locus Rfm1 in barley to a sub-centimorgan genetic interval and a 208-kb physical interval. Rfm1 restores the fertility of msm1 and msm2 male-sterile cytoplasms in barley. The fertility restoration gene is located on the short arm of chromosome 6H (6HS), and we pursued a positional cloning of this gene. Starting from a previous result that has delimited Rfm1 within a 10.8 cM region on 6HS, we developed novel CAPS and SSR markers tightly linked to the gene in barley using the sequence information from the syntenic region of rice and barley genome assemblies. Next, we performed fine mapping of the Rfm1 locus. To isolate recombinants, we surveyed 3,638 F2 plants derived from a cross between the CMS strain and the Rf strain with adjacent markers (NAS2090 and NAS1080). This analysis identified 175 recombinant plants from the F2 population to build a high-resolution map with nine markers tightly linked to the Rfm1 locus. Rfm1 was located within the 0.14 cM region delimited by two markers (NAS9113 and NAS9200). Using these flanking markers as well as marker cosegregating with Rfm1 (NAS9133), we screened the BAC libraries of the cultivar Morex, an rfm1 carrier. We isolated 11 BAC clones and constructed a BAC physical map using their fingerprints. Finally, we delimited the Rfm1 locus encompassing the rfm1 allele on a 208-kb contig composed of three minimally overlapping BAC clones. This precise localization of the Rfm1 locus in the barley genome is expected to greatly accelerate the future map-based cloning of the Rfm1 gene by sequence analysis and its genetic transformation for the complementation of cytoplasmic male-sterile plants.

An ATP-binding cassette subfamily G full transporter is essential for the retention of leaf water in both wild barley and rice.

Land plants have developed a cuticle preventing uncontrolled water loss. Here we report that an ATP-binding cassette (ABC) subfamily G (ABCG) full transporter is required for leaf water conservation in both wild barley and rice. A spontaneous mutation, eibi1.b, in wild barley has a low capacity to retain leaf water, a phenotype associated with reduced cutin deposition and a thin cuticle. Map-based cloning revealed that Eibi1 encodes an HvABCG31 full transporter. The gene was highly expressed in the elongation zone of a growing leaf (the site of cutin synthesis), and its gene product also was localized in developing, but not in mature tissue. A de novo wild barley mutant named "eibi1.c," along with two transposon insertion lines of rice mutated in the ortholog of HvABCG31 also were unable to restrict water loss from detached leaves. HvABCG31 is hypothesized to function as a transporter involved in cutin formation. Homologs of HvABCG31 were found in green algae, moss, and lycopods, indicating that this full transporter is highly conserved in the evolution of land plants.

Divergence of expression pattern contributed to neofunctionalization of duplicated HD-Zip I transcription factor in barley.

Barley (Hordeum vulgare) spikes are developmentally switched from two-rowed to six-rowed by a single recessive gene, six-rowed spike 1 (vrs1), which encodes a homeodomain-leucine zipper I class transcription factor. Vrs1 is a paralog of HvHox2 and both were generated by duplication of an ancestral gene. HvHox2 is conserved among cereals, whereas Vrs1 acquired its current function during the evolution of barley. It was unclear whether divergence of expression pattern or protein function accounted for the functionalization of Vrs1. Here, we conducted a comparative analysis of protein functions and gene expression between HvHox2 and Vrs1 to clarify the functionalization mechanism. We revealed that the transcriptional activation activity of HvHOX2 and VRS1 was conserved. In situ hybridization analysis showed that HvHox2 is localized in vascular bundles in developing spikes, whereas Vrs1 is expressed exclusively in the pistil, lemma, palea and lodicule of lateral spikelets. The transcript abundance of Vrs1 was > 10-fold greater than that of HvHox2 during the pistil developmental stage, suggesting that the essential function of Vrs1 is to inhibit gynoecial development. We demonstrated the quantitative function of Vrs1 using RNAi transgenic plants and Vrs1 expression variants. Expression analysis of six-rowed spike mutants that are nonallelic to vrs1 showed that Vrs1 expression was up-regulated by Vrs4, whereas HvHox2 expression was not. These data demonstrate that the divergence of gene expression pattern contributed to the neofunctionalization of Vrs1.

An alternative mechanism for cleistogamy in barley.

Cleistogamy in barley is genetically determined by the presence of the recessive allele cly1, but the dominant allele at the linked locus Cly2 is epistatic over cly1. Although the molecular basis for cly1 action is well understood, that of Cly2 is not. Here we show that anther non-extrusion can occur not just when the lodicules fail to expand adequately (a trait which is fully determined by the allelic state at the cly1 locus), but by the premature timing of anthesis before the spike has emerged from the boot. The transcription of HvAP2 at cly1 is unaffected by the timing of anthesis. Where this occurs prematurely, by the time that the spike has emerged from the boot, the lodicules have already become shrunken and have lost the capacity to push the lemma and palea apart. Premature anthesis appears to be governed by a dominant gene, probably Cly2. Of the three phases of development of a non-cleistogamous barley floret (spike emergence from the boot, floret gaping induced by lodicule expansion and anther extrusion), genetic variation is available regarding at least the former two.

Structure, transcription and post-transcriptional regulation of the bread wheat orthologs of the barley cleistogamy gene Cly1.

The majority of genes present in the hexaploid bread wheat genome are present as three homoeologs. Here, we describe the three homoeologous orthologs of the barley cleistogamy gene Cly1, a member of the AP2 gene family. As in barley, the wheat genes (designated TaAP2-A, -B and -D) map to the sub-telomeric region of the long arms of the group 2 chromosomes. The structure and pattern of transcription of the TaAP2 homoeologs were similar to those of Cly1. Transcript abundance was high in the florets, and particularly in the lodicule. The TaAP2 message was cleaved at its miR172 target sites. The set of homoeolog-specific PCR assays developed will be informative for identifying either naturally occurring or induced cleistogamous alleles at each of the three wheat homoeologs. By combining such alleles via conventional crossing, it should be possible to generate a cleistogamous form of bread wheat, which would be advantageous both with respect to improving the level of the crop's resistance against the causative pathogen of fusarium head blight, and for controlling pollen-mediated gene flow to and from genetically modified cultivars.

An eceriferum locus, cer-zv, is associated with a defect in cutin responsible for water retention in barley (Hordeum vulgare) leaves.

Drought limits plant growth and threatens crop productivity. A barley (Hordeum vulgare) ethylene imine-induced monogenic recessive mutant cer-zv, which is sensitive to drought, was characterized and genetically mapped in the present study. Detached leaves of cer-zv lost 34.2 % of their initial weight after 1 h of dehydration. The transpiration was much higher in cer-zv leaves than in wild-type leaves under both light and dark conditions. The stomata of cer-zv leaves functioned normally, but the cuticle of cer-zv leaves showed increased permeability to ethanol and toluidine blue dye. There was a 50-90 % reduction in four major cutin monomers, but no reduction in wax loads was found in the cer-zv mutant as compared with the wild type. Two F(2) mapping populations were established by the crosses of 23-19 × cer-zv and cer-zv × OUH602. More polymorphisms were found in EST sequences between cer-zv and OUH602 than between cer-zv and 23-19. cer-zv was located in a pericentromeric region on chromosome 4H in a 10.8 cM interval in the 23-19 × cer-zv map based on 186 gametes tested and a 1.7 cM interval in the cer-zv × OUH602 map based on 176 gametes tested. It co-segregated with EST marker AK251484 in both maps. The results indicated that the cer-zv mutant is defective in cutin, which might be responsible for the increased transpiration rate and drought sensitivity, and that the F(2) of cer-zv × OUH602 might better facilitate high resolution mapping of cer-zv.

Cleistogamous flowering in barley arises from the suppression of microRNA-guided HvAP2 mRNA cleavage.

The cleistogamous flower sheds its pollen before opening, forcing plants with this habit to be almost entirely autogamous. Cleistogamy also provides a means of escape from cereal head blight infection and minimizes pollen-mediated gene flow. The lodicule in cleistogamous barley is atrophied. We have isolated cleistogamy 1 (Cly1) by positional cloning and show that it encodes a transcription factor containing two AP2 domains and a putative microRNA miR172 targeting site, which is an ortholog of Arabidopsis thaliana AP2. The expression of Cly1 was concentrated within the lodicule primordia. We established a perfect association between a synonymous nucleotide substitution at the miR172 targeting site and cleistogamy. Cleavage of mRNA directed by miR172 was detectable only in a noncleistogamous background. We conclude that the miR172-derived down-regulation of Cly1 promotes the development of the lodicules, thereby ensuring noncleistogamy, although the single nucleotide change at the miR172 targeting site results in the failure of the lodicules to develop properly, producing the cleistogamous phenotype.

Variation in the wheat AP2 homoeologs, the genes underlying lodicule development.

The bread wheat genome harbors three homoeologs of the barley gene HvAP2, which determines the cleistogamous/non-cleistogamous flowering. The three homoeologs, TaAP2-A, TaAP2-B and TaAP2-D, are derived from the A, B and D genomes. The importance of lodicule swelling in assuring non-cleistogamous flowering in a range of wild and domesticated wheat accessions of varying ploidy level was established. Re-sequencing of wheat AP2 homoeologous genes was carried out to identify natural variation at both the nucleotide and polypeptide level. The sequences of wheat AP2 homoeologs are highly conserved even across different ploidy levels and no functional variants at the key miR172 targeting site were detected. These results indicate that engineering of cleistogamous wheat will require the presence of a functional TaAP2 modification at each of the three homoeologs.

Population-genetic analysis of HvABCG31 promoter sequence in wild barley (Hordeum vulgare ssp. spontaneum).

BACKGROUND: The cuticle is an important adaptive structure whose origin played a crucial role in the transition of plants from aqueous to terrestrial conditions. HvABCG31/Eibi1 is an ABCG transporter gene, involved in cuticle formation that was recently identified in wild barley (Hordeum vulgare ssp. spontaneum). To study the genetic variation of HvABCG31 in different habitats, its 2 kb promoter region was sequenced from 112 wild barley accessions collected from five natural populations from southern and northern Israel. The sites included three mesic and two xeric habitats, and differed in annual rainfall, soil type, and soil water capacity. RESULTS: Phylogenetic analysis of the aligned HvABCG31 promoter sequences clustered the majority of accessions (69 out of 71) from the three northern mesic populations into one cluster, while all 21 accessions from the Dead Sea area, a xeric southern population, and two isolated accessions (one from a xeric population at Mitzpe Ramon and one from the xeric 'African Slope' of "Evolution Canyon") formed the second cluster. The southern arid populations included six haplotypes, but they differed from the consensus sequence at a large number of positions, while the northern mesic populations included 15 haplotypes that were, on average, more similar to the consensus sequence. Most of the haplotypes (20 of 22) were unique to a population. Interestingly, higher genetic variation occurred within populations (54.2%) than among populations (45.8%). Analysis of the promoter region detected a large number of transcription factor binding sites: 121-128 and 121-134 sites in the two southern arid populations, and 123-128,125-128, and 123-125 sites in the three northern mesic populations. Three types of TFBSs were significantly enriched: those related to GA (gibberellin), Dof (DNA binding with one finger), and light. CONCLUSIONS: Drought stress and adaptive natural selection may have been important determinants in the observed sequence variation of HvABCG31 promoter. Abiotic stresses may be involved in the HvABCG31 gene transcription regulations, generating more protective cuticles in plants under stresses.