Diversity in boron toxicity tolerance of Australian barley (Hordeum vulgare L.) genotypes
© Hayes et al. 2015
Received: 18 June 2015
Accepted: 6 September 2015
Published: 26 September 2015
Boron (B) is an important micronutrient for plant growth, but is toxic when levels are too high. This commonly occurs in environments with alkaline soils and relatively low rainfall, including many of the cereal growing regions of southern Australia. Four major genetic loci controlling tolerance to high soil B have been identified in the landrace barley, Sahara 3771. Genes underlying two of the loci encode the B transporters HvBot1 and HvNIP2;1.
We investigated sequence and expression level diversity in HvBot1 and HvNIP2;1 across barley germplasm, and identified five novel coding sequence alleles for HvBot1. Lines were identified containing either single or multiple copies of the Sahara HvBot1 allele. We established that only the tandemly duplicated Sahara allele conferred B tolerance, and this duplicated allele was found only in a set of nine lines accessioned in Australian collections as Sahara 3763–3771. HvNIP2;1 coding sequences were highly conserved across barley germplasm. We identified the likely causative SNP in the 5’UTR of Sahara HvNIP2;1, and propose that the creation of a small upstream open reading frame interferes with HvNIP2;1 translation in Sahara 3771. Similar to HvBot1, the tolerant HvNIP2;1 allele was unique to the Sahara barley accessions. We identified a new source of the 2H B tolerance allele controlling leaf symptom development, in the landrace Ethiopia 756.
Ethiopia 756, as well as the cultivar Sloop Vic which carries both the 2H and HvBot1 B tolerance alleles derived from Sahara 3771, may be valuable as alternative parents in breeding programs targeted to high soil B environments. There is significant diversity in B toxicity tolerance among contemporary Australian barley varieties but this is not related to variation at any of the four known B tolerance loci, indicating that novel, as yet undiscovered, sources of tolerance exist.
High soil boron (B) can affect yields of barley (Hordeum vulgare L.) across southern Australia by up to 17 % , depending on a multitude of site, seasonal and genetic factors [2–5]. Along with disease ratings, standard information for new barley varieties released in South Australia, Victoria and Western Australia often includes a boron tolerance rating [6, 7], allowing farmers to select varieties tolerant to high soil B. Genetic variability for high B tolerance has long been known [8, 9]. The most tolerant barley identified amongst breeding material in Australia is the unadapted six-row North African landrace, Sahara 3771. This genotype was accessioned in Australian collections in the early 1900s , one of a set of nine barley lines listed as Sahara 3763 – Sahara 3771. It has been considered an important source of B tolerance for barley breeding programs over many years. Four major QTL for B tolerance were identified in Sahara 3771, in a genetic study using a doubled haploid (DH) population derived from a cross between the South Australian malting variety, Clipper, and Sahara 3771 . Subsequent research to fine-map two of the regions revealed the identity of the tolerance genes HvBot1 (chromosome 4H)  and HvNIP2;1 (chromosome 6H) . They encode two types of transporter that function to minimise the amount of B in barley roots. These genes have been partially characterised, but the prevalence of the tolerant alleles across Australian germplasm was not known. It was also not known if there is significant diversity in HvBot1 and HvNIP2;1 contributing to B tolerance, such as has been found in wheat for TaBot-B5 . Therefore, the aims of this study were to: 1) determine the prevalence of known B tolerance alleles in Australian barley germplasm; 2) develop an improved set of markers for tracking the introgression of B tolerance from Sahara 3771; and 3) identify alternative sources of B tolerance in barley. This was a broad study and, although the set of germplasm assessed was not exhaustive, our data suggest that the tolerance alleles found in Sahara 3771 are rare. The significance of a QTL on chromosome 2H controlling leaf symptom expression is highlighted as a target for future breeding and selection for B tolerance in barley.
Tolerance to high soil B in the barley landrace Sahara 3771 has been attributed to four major QTL, on chromosomes 2H, 3H, 4H and 6H. We screened a set of 65 diverse barley genotypes (Additional file 1: Table S1) for variation at these loci using genomic Southern analysis, which also enabled us to assess gene copy number variation. Coding sequence for the genes encoding B transport proteins HvBot1 and HvNIP2;1, and which lie beneath the 4H and 6H tolerance loci, respectively, was also amplified and sequenced. In sourcing diverse germplasm to screen, we obtained seed for nine barleys accessioned in the Australian Grains Genebank as Sahara 3763 to Sahara 3771. Our analyses suggest that the Sahara accessions possess a unique set of B tolerance alleles.
Genetic variation at the 4H locus (HvBot1)
Southern analysis of the 65 barley genotypes revealed that the tandem HvBot1 gene duplication found in Sahara 3771 is rare, although other genotypes (eg. California Mariout and derivative cultivars CM67 and CM72, and the Japanese cultivars Haruna Nijo and Amagi Nijo) showed a Sahara-like restriction pattern without gene duplication (panel A, Additional file 2: Figure S1). All nine Sahara genotypes from the Australian Grains Genebank possessed the HvBot1 gene duplication, and displayed similarly high B-tolerant phenotypes in hydroponic experiments (Additional file 2: Figure S1).
We compared HvBot1 expression levels in roots of seedlings of barley genotypes carrying the different alleles (Fig. 1b). HvBot1 expression was similar for all variants except for genotypes carrying the duplicated Sahara allele (Bot1(Dp).b), which in this experiment showed more than 200-fold higher levels of expression (Fig. 1b). Multiple sequencing efforts failed to uncover any sequence differences, in either the HvBot1 gene or its promoter, between copies in lines carrying the HvBot1 duplication, or between the duplicated allele and lines containing a single copy of the Sahara HvBot1 allele.
We also assessed the function of each of the variants by heterologous expression in yeast (Saccharomyces cereviseae). With the exception of the Haruna Nijo allele (Bot1.d), all variants conferred a similar level of tolerance to media containing high concentrations of B when expressed in yeast (Panel A in Additional file 3: Figure S2). Clones containing the Haruna Nijo variant were no more tolerant to high B than empty vector control yeast, indicating that one or both of the residue substitutions in the Haruna Nijo HvBot1 allele disrupted protein function. We mutated each of the residues separately in the Sahara HvBot1 ORF and tested the two mutations in a second yeast expression experiment, which demonstrated that Leucine at position 234 in HvBot1 is critical for function (Fig. 1c). By contrast, the Thr541Met substitution had no effect on transporter function. The cultivars Haruna Nijo and Amagi Nijo which carry the Haruna Nijo HvBot1 allele are intolerant to high B when grown in hydroponics. Root B concentrations were similar to Clipper, and suggested that no net efflux of B from the roots was occurring in either cultivar (Panel B in Additional file 3: Figure S2).
We developed a KASP™ marker to track a SNP in exon 11 unique to Sahara (wri57; Additional file 4: Table S2, with marker location shown in Fig. 1a). The previously reported xBot1 marker  (location shown in Fig. 1a) was designed around a SNP that is not specific to the Sahara HvBot1 allele, but is also common to the Alexis, Morex, WI4304 and Haruna Nijo alleles (Fig. 1a).
Genetic variation at the 6H locus (HvNIP2;1)
A second SNP, unique to the Sahara accessions and located 44 bp upstream of the ATG start codon in the 5’ UTR, was identified in HvNIP2;1. At this position the thymidine base in Sahara creates a short, upstream ORF (uORF) encoding a translatable peptide of thirteen residues immediately upstream of, and in frame with the HvNIP2;1 ORF (Fig. 3a). We hypothesised that this uORF in Sahara may interfere with HvNIP2;1 translation, and be the principal mechanism for reduced function of HvNIP2;1 leading to reduced B permeability and a greater tolerance to high B. The 133 bp 5’UTR sequences from Clipper and Sahara HvNIP2;1 were cloned into the SP6 vector, immediately upstream of the firefly (Photinus pyralis) luciferase gene driven by the SP6 RNA polymerase promoter. In vitro transcription/translation reactions with each construct demonstrated that the Sahara 5’ UTR sequence inhibited translation of the luciferase reporter gene by 43 % relative to the Clipper 5’ UTR sequence (Fig. 3c).
Interestingly, the uORF SNP was found only in Sahara 3771 and the other eight Sahara accessions, suggesting that tolerance to high soil B attributable to this locus is very rare. We developed a KASP™ marker to track the Sahara uORF SNP in HvNIP2;1 (wri59; Additional file 4: Table S2).
Genetic variation at the 3H and 2H loci
The genetic determinants at these loci are unknown. However, a member of the Bot B transporter gene family, HvBot2, is located within the region of the 3H QTL controlling relative root length at high B and is considered to be a likely candidate gene  ([GenBank:KR817581] lists the coding sequence for HvBot2 from Clipper). The HvBot2 sequence in Sahara 3771 is disrupted by a large deletion (approx. 5.1 kB) beginning at the 3’ end of the CDS, which may result in a loss of function. We designed two sets of primers, to amplify within and across the deletion, respectively (Additional file 5: Table S3). Discriminating PCR determined that with the exception of the nine Sahara accessions, the HvBot2 CDS was intact in all tested barley genotypes. A KASP™ marker was developed to track the Sahara 3771 deletion (wri58; Additional file 4: Table S2).
The 2H QTL for B tolerance controls leaf symptom expression in the Clipper X Sahara 3771 population and is located close to the centromere . Fine mapping in this region is difficult due to low recombination frequency. Furthermore, within the Clipper X Sahara 3771 DH population the QTL was poorly defined, largely because leaf symptoms could not be accurately assessed in lines carrying Sahara alleles at the other B tolerance loci. Additional background genetic effects might also contribute to an unclear phenotype in the DH lines. We identified seven DH lines that were recombinant within the 2H interval and backcrossed them to Clipper. BC1F2 plants retaining the recombinant Sahara haplotype across the 2H region but with Clipper alleles at the 3H, 4H and 6H loci were selected by genotyping. Evaluation of leaf symptoms of B toxicity in BC1F3 progeny clearly defined the interval for the 2H QTL to approximately 10.3 cM, between the markers GBMS0160 and Bmag381 (Additional file 6: Fig. S3, Panel A). We tested the markers contig888 and G9-138C in an attempt to characterise the barley genotypes at the 2H QTL region. G9-138C appeared to associate with a B tolerance phenotype. For this marker, the Sahara accessions and the line Ethiopia 756 shared a similar restriction pattern (Additional File 6: Fig. S3, Panel B).
Prevalence and impact of B tolerance alleles in Australian barley varieties
A set of 80 Australian barley varieties included in National Variety Trial testing between 2008 and 2012 were screened for allele type at each of the four known B tolerance loci, leaf symptoms of B toxicity, and leaf B concentrations when grown in a glasshouse with an elevated supply of B. The data for each variety included in the screen are provided in Additional file 7: Table S4. A small number of varieties were heterogeneous for either the HvNIP2;1 or HvBot1 alleles, and were separated into sub-lines. The most common HvBot1 allele amongst the Australian barley varieties was the Clipper allele (67 genotypes). The Morex (10 genotypes) and Alexis (7 genotypes) alleles were also relatively common. None of the genotypes contained the Sahara HvNIP2;1 allele on chromosome 6H or the Sahara variant of HvBot2 on chromosome 3H.
In our varietal screen, several other barley cultivars also showed both low leaf symptoms of B toxicity and low leaf B concentrations, including Fleet, Capstan and Buloke (Fig. 4; Additional file 7: Table S4). We phenotyped a barley DH population (233 lines) derived from reciprocal crosses of Fleet with Commander. QTL analysis failed to reveal an association between leaf symptoms of B toxicity and any of the known B tolerance loci (data not shown). Thus, the chromosomal regions responsible for B tolerance observed in Parent 19, Fleet and other cultivars with lower levels of tolerance remain unknown, and do not appear to correspond to any of the B tolerance QTL found in Sahara 3771.
The identities of two genes conferring tolerance to high soil B in barley have been determined [12, 13, 16]. One encodes HvNIP2;1, a transport protein belonging to the NIP family of aquaporins, and is found on chromosome 6H. The other is on chromosome 4H and encodes a different kind of transporter, HvBot1. The tolerant alleles for both genes are present in a set of barley accessions from North Africa named Sahara, including the well-studied genotype Sahara 3771. A genetic study using a population derived from an intolerant barley variety, Clipper, and Sahara 3771, identified four chromosome locations associated with B tolerance, with all of the favourable alleles coming from Sahara 3771 . In this study, we investigated the prevalence of these four alleles in barley germplasm, and the allelic diversity of HvBot1 and HvNIP2;1. Amongst the studied germplasm, we found that the Sahara lines contain a rare set of tolerance alleles, including a tandem duplication of HvBot1 (4H), a critical 5’UTR SNP in the coding sequences of HvNIP2;1 (6H), a deletion in the gene HvBot2 that is located under the 3H tolerance locus, and a rare haplotype across the 2H QTL region.
Tolerance due to HvBot1 (chromosome 4H) was originally attributed to either greater functionality of the transporter as a result of sequence differences between the Sahara and Clipper HvBot1 coding sequences, or higher levels of expression . In addition to promoter sequence differences, tandem duplication of HvBot1 in the genome of Sahara has created an estimated four copies. Our study indicates that this duplication occurred relatively recently in barley evolution, and is rare. Southern analysis showed that the duplication is only present in Sahara 3771, eight other genotypes also accessioned with the Sahara name, and two Australian cultivars that contain Sahara 3771 in their pedigree. Sequence analysis of many genotypes in this study also revealed that the original published HvBot1 sequence for Clipper contained a sequencing error. There are in fact no residue differences between the alleles from Clipper and Sahara 3771, and both alleles confer a similar level of B tolerance to yeast when heterologously expressed (Fig. 1). All HvBot1 sequences lodged with GenBank have been amended and updated accordingly. Thus, differences between the two alleles are entirely attributable to expression level differences.
We identified seven HvBot1 CDS alleles in total. Only one of these, found in two Japanese cultivars Haruna Nijo and Amagi Nijo, showed reduced function when expressed in yeast. This is due to a residue substitution at position 234 in the non-functional Haruna Nijo allele, of histidine for leucine (Fig. 1c). All single-copy HvBot1 alleles were expressed at similarly low levels, including the single-copy Sahara HvBot1 allele (Fig. 1b). Genomic sequencing of a number of single-copy Sahara HvBot1 genotypes failed to identify any sequence differences between the single-copy and multi-copy Sahara alleles, or between copies present in Sahara. Moreover, a line derived from the Clipper X Sahara 3771 DH population and with a recombination within the HvBot1 gene cluster, showed a dependency of HvBot1 gene expression level on the number of Sahara HvBot1 copies present (Fig. 2). We conclude that the increased expression of HvBot1 found in Sahara is due to the duplication of this gene. A number of other specific examples of gene copy number expansion occurring through evolution have been described in plant genomes, which have also resulted in increased gene expression [17–19]. Recent whole genome analyses have revealed that copy number variation is common in plants, and appears to be biased towards genes involved in abiotic and biotic stress responses [20, 21].
Further, we propose that duplication of the Sahara HvBot1 allele is responsible for B tolerance at the 4H locus. Segregation analysis in two independent F2 barley populations (Parent 19 × Clipper; Ethiopia 756 × Clipper) failed to demonstrate linkage between the presence of a single copy of the Sahara HvBot1 allele and B tolerance as determined by leaf symptom expression (Fig. 6). By contrast, in the original genetic study the duplicated HvBot1 allele from Sahara was associated with leaf symptom development, as well as a number of other traits associated with B tolerance .
We also examined diversity in coding sequence and gene expression level for HvNIP2;1, the gene underlying a major QTL in Sahara 3771 on chromosome 6H , revealed that the HvNIP2;1 sequences were highly conserved amongst the studied germplasm. We identified only two SNPs in coding sequence of this gene; one synonymous SNP that was present in around half of the genotypes examined and could not be related to a B tolerance phenotype, and a second SNP in the 5’UTR that was unique to the set of nine Sahara accessions (Fig. 3a). We propose that the 5’UTR SNP is causative of B tolerance in Sahara. An in vitro transcription/translation assay showed reduced luciferase activity in the presence of the 5’UTR derived from Sahara HvNIP2;1 (Fig. 3c), suggesting that a small uORF created by the Sahara 5’UTR SNP interferes with translation of HvNIP2;1. In our previous study we detected a modest reduction in HvNIP2;1 transcript levels in root tissues of B-treated Sahara 3771 compared to Clipper . However, in broadening the analysis of expression to a wider range of genotypes, we were not able to reproduce these differences (Fig. 3b). Future experiments to quantify active protein levels in roots of Sahara and other, intolerant barley genotypes will help to confirm if the uORF in Sahara inhibits HvNIP2;1 translation. Upstream ORFs are common elements of plant transcript sequences and have been shown to mediate translation for a number of genes (reviewed in ), including transcription factors [23, 24] and enzymes involved in polyamine biosynthesis .
It has been suggested that when using Sahara as a breeding parent, there may be linkage drag of unfavourable traits together with the 4H allele, or pleiotropic effects of the 4H tolerance allele in different backgrounds [26, 27]. However, we identified two released barley varieties carrying the 4H allele from Sahara, Navigator and Sloop Vic (Additional file 7: Table S4), demonstrating there are unlikely to be severe penalties in either yield or quality traits associated with the introgression of the HvBot1 gene duplication. We investigated by genotype analysis the extent of the Sahara-derived 4H segment in Navigator (WI4262; released 2009), and estimate it to cover a region of at least 14 cM. While developed for cultivation in higher rainfall, long-season environments not typically associated with B toxicity, Navigator consistently shows low levels of leaf symptom expression when challenged with high B (this study and ).
The cultivar Sloop Vic was specifically selected for B tolerance . We identified genetic heterogeneity for the 4H locus in this variety (Fig. 5), but also the presence of a chromosomal segment from Sahara that includes the 2H B tolerance allele. The two other tolerance alleles from Sahara (3H and 6H) were not retained in Sloop Vic. Our findings suggest that the 2H allele is highly useful for imparting B tolerance. We also identified another unadapted six-row barley, Ethiopia 756, as an alternative source of the 2H tolerance allele. Ethiopia 756 shares the same haplotype for a number of markers across the 2H interval with Sahara, including a rare G9-138C allele. Segregation analysis in a population derived from Ethiopia 756 crossed with Clipper revealed an association between the critical 2H region in Ethiopia 756 and leaf symptoms of B toxicity (Fig. 6). In an alternative population derived from Parent 19 and Clipper the same region was not linked with B tolerance.
This study of allelic diversity for the four known boron tolerance loci revealed that Sahara 3771 has a unique set of boron tolerance genes, rare amongst germplasm. It allowed us to identify the causative features of two of the genes responsible for B tolerance in Sahara. These are a tandem duplication of HvBot1 and a short, upstream open reading frame in the coding sequences of HvNIP2;1. Our findings facilitate the development of markers for tracking B tolerance in barley with greater specificity, and we designed three KASP™ marker assays for tracking tolerance derived from Sahara. Our detection of the 2H locus in Ethiopia 756 is the first reported validation of the Sahara 3771 B tolerance QTL in another population.
The genetic locations controlling B tolerance in Parent19, as well as lower levels of tolerance in cultivars such as Fleet, Buloke and Capstan that were observed in our study and in the field by others [6, 7], remain unknown; selection for B tolerance derived from these sources can only be made on the basis of phenotype. However, the genotypes Ethiopia 756 and Sloop Vic may be used in breeding programs as sources of B tolerance attributable to a significant QTL for leaf symptom expression on chromosome 2H, and may be valuable alternatives to Sahara.
Seed for most genotypes of barley (Hordeum vulgare L.) was obtained from either the Australian Grains Genebank or from our own collections. Seed of Parent 19 and a number of ICARDA genotypes showing B tolerance in the field was provided by Drs Jason Eglinton and Stewart Coventry, The University of Adelaide’s Barley Breeding Program. The varieties Compass, Flinders, Grange, Henley, LaTrobe, SY Rattler and Westminster were obtained from Dr Hugh Wallwork (South Australian Research and Development Institute), and Macumba was obtained from Amanda Box (The University of Adelaide).
Allele diversity analyses
Seeds were germinated on moist filter paper in the laboratory. After 4–7 days, whole roots and shoots from five uniform seedlings for each genotype were pooled and snap-frozen in liquid nitrogen, for RNA and DNA extractions, respectively. Pooling from a number of individual seedlings facilitated the detection of genetic heterogeneity within lines. The 65 genotypes included for allele diversity analysis are listed in Additional file 1: Table S1.
RNA was extracted using TRIzol (Invitrogen). For cDNA synthesis, we used Superscript III Reverse Transcriptase (Invitrogen); samples for quantitative real-time reverse transcription polymerase chain reaction (q-RT-PCR) were DNase-treated prior to cDNA synthesis using DNA-Free (Ambion, USA). Coding sequence for HvBot1 and HvNIP2;1 was amplified from cDNA obtained from 68 genotypes, and Sanger sequenced at the Australian Genome Research Facility (AGRF; Adelaide). q-RT-PCR analysis of expression of HvBot1 and HvNIP2;1 was performed as previously described , using gene-specific primers (Additional file 5: Table S3). For some experiments, semi-quantitative reverse transcription PCR (semi-q-RT-PCR) was performed on cDNAs, where we used HvGAP amplification (Additional file 5: Table S3) to indicate between-sample variation in template concentration.
Genomic DNA was isolated using phenol-chloroform extraction, digested with the restriction enzyme Dra I and the products separated by gel electrophoresis before transfer to a nylon membrane for Southern hybridisation using standard methods. Membranes were probed with radio-labelled nucleic acid fragments of the B tolerance genes HvBot1 and HvNIP2;1, and the putative tolerance gene HvBot2 (co-locating with a tolerance QTL on chromosome 3H). Several probes located within the region of the 2H B tolerance QTL were also utilised. Probe details are included in Additional file 5: Table S3. Membranes were stripped with a 0.1 % sodium dodecyl sulphate (SDS), 2 mM ethylenediaminetetraacetic acid (EDTA) solution heated to boiling point between each sequential hybridisation.
Heterologous expression and in vitro transcription/translation experiments
Open reading frames for each of the identified HvBot1 alleles were PCR-amplified from cDNA prepared from roots of the barley genotypes Sahara 3771, Clipper, Haruna Nijo, WI4304, Alexis and Tadmor using primers listed in Additional file 5: Table S3, and products were cloned in the Gateway entry vector pCR8 (Invitrogen, Carlsbad, CA, USA). Sequences of the inserts were verified by Sanger sequencing before recombining into the destination vector pYES3.DEST (Invitrogen) for yeast expression. The Sahara HvBot1 expression vector was also used as a template for in vitro site-directed mutagenesis (Quikchange II Site-Directed Mutagenesis kit, Stratagene, La Jolla CA, USA) using specific primer pairs detailed in Additional file 5: Table S3. Two variants of HvBot1 were created, each with a single residue substitution of either Leu234His or Thr541Met. Both substitutions are present in the non-functional Haruna Nijo HvBot1 allele. Experiments for functional assessment of the HvBot1 alleles and variants in yeast (Saccharomyces cerevisiae) were performed as previously described .
5’UTR sequences from HvNIP2;1 genes of Clipper and Sahara were amplified from cDNA using primers listed in Additional file 5: Table S3. The PCR products and SP6 control luciferase vector DNA (Promega, USA) were digested with restriction enzymes Not I and Bam HI. Digested products were then cloned into the linearised SP6 vector, between the SP6 RNA polymerase promoter and the firefly (Photinus pyralis) luciferase gene. Insert sequences were confirmed by Sanger sequencing. We then used a TnT Coupled Wheat Germ Extract System (Promega, USA) to perform in vitro transcription/translation reactions using SP6 control DNA and the modified vector. We made six reactions for each construct. Translation of luciferase protein for each reaction was measured as luciferase activity, by recording luminescence using a POLARstar Optima plate reader (BMG LabTech) programmed to perform multiple reads on each sample at 1 s intervals for a total of approximately 2 min. All plotted data were linear with slopes close to zero, and were extrapolated to time 0 to determine initial rates of activity. Four luminescence assays were recorded for each transcription/translation reaction.
Screening of current varieties
Seeds were sown into 50:50 UC mix: coco peat potting medium, in 25 cm diameter pots without replication. Following germination, each pot was thinned to three uniform seedlings. Plants were grown through to maturity during spring, in a temperature-controlled greenhouse (day/night average temperatures of 23 °C/19 °C) with natural lighting. Pots were moved during the growing period to avoid positional effects on light interception or transpiration, which might influence the development of symptoms of B toxicity. Plants were watered regularly with tap water, and fertilised twice with a multi-element, slow-release fertiliser. Tap water delivered to the University of Adelaide’s Waite Campus between August and December 2013 contained B concentrations of 0.1–0.5 mg B L−1, which was sufficient to induce B toxicity-attributable leaf necrosis in intolerant barley varieties in this experiment. We have consistently observed B toxicity symptoms in barley plants grown across the Waite Campus since 2012, co-incident with elevated B concentrations in tap water (data not supplied). In this experiment, symptoms of B toxicity were assessed visually three times during the growth of plants to full maturity and an average score determined (0 = no necrosis; 6 = severe necrosis). The penultimate leaves from five tillers were sampled at mid-grain fill, oven-dried and chopped finely with scissors. A sub-sample of dried leaf tissue was acid-digested for analysis of B concentration by Radial View Inductively Coupled Plasma-Optical Emission Spectrometry (ICP-OES) . Selected varieties were grown hydroponically with high B using a nutrient solution described in . B concentrations in root tissues of these plants were determined using an azomethine-H colorimetric assay .
Development of KASP™ markers to detect the Sahara 3771 allele at the HvBot1, HvBot2 and HvNIP2;1 loci
KASP™ markers were designed based on polymorphisms between Sahara 3771 and B intolerant cultivars at the B tolerance loci HvBot1, HvBot2 and HvNIP2;1 (Additional file 4: Table S2). The markers wri57 and wri59 (for HvBot1 and HvNIP2;1, respectively) assay SNP polymorphisms between the Sahara accessions and B intolerant lines, while the wri58 marker (for HvBot2) assays a deletion present only in Sahara accessions. All primer sets were designed using Kraken (LGC Limited, London, UK) with default parameters. Assays were validated in the parental lines and a set of varieties for which alleles were previously determined using an automated SNPLine system (LGC Limited, London, UK).
F2 segregant analysis of B toxicity symptoms
Crosses were made using two barley lines with low levels of leaf B toxicity symptoms as female parents (Parent 19 and Ethiopia 756), and the intolerant variety Clipper as the male parent. F1 hybrids from these crosses were confirmed by genotyping using the xBot1 CAPS marker (Additional file 5: Table S3). F2 progeny were germinated on moist filter paper, and transferred to aerated nutrient solution containing 3 mM B. Seedlings were grown for 15 d in a growth chamber (18 °C/15 °C and 50 %/70 % humidity day/night settings, with a 12 h photoperiod and light intensity of approx. 200 mol m−2 s−1 at plant height). Individual seedlings were then scored for percentage necrosis on the first leaf, and sampled for DNA isolation using a freeze-dried method of extraction . Material was genotyped at the 4H and 2H B tolerance loci using the KASP™ assays wri57 (Additional file 4: Table S2) and ABC02403 (Additional file 5: Table S3), respectively. The 2H KASP™ marker was designed and validated as discriminating between Clipper and the B tolerant parents, and verified by genotyping of selected recombinant lines to reside within the defined QTL interval. Association between genotype and phenotype for each population and each of the 4H and 2H regions was analysed by separate one-way ANOVA.
Availability of supporting data
The data sets supporting the results of this article are included within the article and its additional files. Nucleic acid sequences described in this research have been lodged with GenBank: [GenBank:EF660435, GenBank:EF660437, GenBank:KR605456, GenBank:KR605457, GenBank:KR605458, GenBank:KR605459, GenBank:KR605460, GenBank:KR605461, GenBank:KR817581].
We thank Alison Hay for technical assistance, Ute Baumann, Carolyn Schultz and Nadim Shadiac for their advice regarding the HvNIP2;1 in vitro transcription/translation experiment, and Penny Tricker for critically reading the manuscript. Tefera Angessa (University of Western Australia) suggested inclusion of genotype Ethiopia 756 in screening for B tolerance, and Stewart Coventry (The University of Adelaide’s Barley Breeding Program) provided seed of Parent 19 and a selection of ICARDA lines observed to have B tolerance in field screens. Waite Analytical Services performed the acid extraction and ICP-OES analysis of B concentrations in leaf material, and Yuan Li and Stephen Fletcher performed the q-RT-PCR experiments. We also acknowledge Tim March for running a QTL analysis on genotyping-by-sequencing marker data for a barley DH population derived from reciprocal crosses of Fleet and Commander. This research was jointly funded by the Australian Research Council, the Grains Research and Development Corporation, The University of Adelaide and the South Australian Government.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Cartwright B, Zarcinas BA, Mayfield AH. Toxic concentrations of boron in a red-brown earth at Gladstone. South Australia Aust J Soil Res. 1984;22:261–72. doi:10.1071/SR9840261.View ArticleGoogle Scholar
- McDonald GK, Eglinton JK, Barr AR. Assessment of the agronomic value of QTL on chromosomes 2H and 4H linked to tolerance to boron toxicity in barley (Hordeum vulgare L.). Plant Soil. 2010;326:275–90. doi:10.1007/s11104-009-0006-1.View ArticleGoogle Scholar
- Reid R, Fitzpatrick K. Influence of leaf tolerance mechanisms and rain on boron toxicity in barley and wheat. Plant Physiol. 2009;151:413–420. doi: http://dx. doi. org/ 10.1104/ pp.109.141069Google Scholar
- Nuttall JG, Hobson KB, Materne M, Moody DB, Munns R, Armstrong RD. Use of genetic tolerance in grain crops to overcome subsoil constraints in alkaline cropping soils. Soil Res. 2010;48:188–199. doi:10.1071/SR09081
- Mertens J, Van Laer L, Salaets P, Smolders E. Phytotoxic doses of boron in contrasting soils depend on soil water content. Plant Soil. 2011;342:73–82. doi:10.1007/s11104-010-0666-×.View ArticleGoogle Scholar
- Paynter B, Hills A, Gupta S, Collins S, Dhammu H, Malik R, Trainor G. 2015 Barley Variety Sowing Guide for Western Australia. Bulletin 4860, WA Department of Agriculture and Food; 2014. ISSN: 1833 7366.Google Scholar
- Wheeler R. Barley variety sowing guide 2015. In: South Australian Sowing Guide 2015. South Australian Research and Development Institute; 2014. pp. 16–20Google Scholar
- Nable RO. Resistance to boron toxicity amongst several barley and wheat cultivars: A preliminary examination of the resistance mechanism. Plant Soil. 1988;112:45–52. doi:10.1007/BF02181751.View ArticleGoogle Scholar
- Paull JG, Cartwright B, Rathjen AJ. Responses of wheat and barley genotypes to toxic concentrations of soil boron. Euphytica. 1988;39:137–44. doi:10.1007/BF00039866.View ArticleGoogle Scholar
- Waterhouse WL. Studies in the inheritance of resistance to leaf rust, Puccinia anomala Rostr., in crosses of barley. I. Reprinted from the J Proc Royal Soc New South Wales, 1927; Vol. LXI:218-247.Google Scholar
- Jefferies S, Barr AR, Karakousis A, Kretshmer JM, Manning S, Chalmers KJ, et al. Mapping of chromosome regions conferring boron toxicity tolerance in barley (Hordeum vulgare L.). Theor Appl Genet. 1999;98:1293–303. doi:10.1007/s001220051195.View ArticleGoogle Scholar
- Sutton T, Baumann U, Hayes J, Collins NC, Shi B-J, Schnurbusch T, et al. Boron toxicity tolerance in barley arising from efflux transporter amplification. Science. 2007;318:1446–9. doi:10.1126/science.1146853.View ArticlePubMedGoogle Scholar
- Schnurbusch T, Hayes JE, Hrmova M, Baumann U, Ramesh SA, Tyerman SD, et al. Boron toxicity tolerance in barley through reduced expression of the multifunctional aquaporin HvNIP2;1. Plant Physiol. 2010;153:1706–15. doi:10.1104/pp.110.158832.
- Pallotta M, Schnurbusch T, Hayes J, Hay A, Baumann U, Paull J, et al. Molecular basis of adaptation to high soil boron in wheat landraces and elite cultivars. Nature. 2014;514:88–91. doi:10.1038/nature13538.View ArticlePubMedGoogle Scholar
- Öz MT. Microarray based expression profiling of barley under boron stress and cloning of 3H boron tolerance gene. PhD Thesis, Middle East Technical University. 2012. Available: https://etd.lib.metu.edu.tr/upload/12614203/index.pdf
- Reid R. Identification of boron transporter genes likely to be responsible for tolerance to boron toxicity in wheat and barley. Plant Cell Physiol. 2007;48:1673–8. doi:10.1093/pcp/pcm159.View ArticlePubMedGoogle Scholar
- Hanikenne M, Talke IN, Haydon MJ, Lanz C, Nolte A, Motte P, et al. Evolution of metal hyperaccumulation required cis-regulatory changes and triplication of HMA4. Nature. 2008;453:391–5. doi:10.1038/nature06877.View ArticlePubMedGoogle Scholar
- Cook DE, Lee TG, Guo X, Melito S, Wang K, Bayless AM, et al. Copy number variation of multiple genes at Rhg1 mediates nematode resistance in soybean. Science. 2012;338:1206–9. doi:10.1126/science.1228746.View ArticlePubMedGoogle Scholar
- Maron LG, Guimarães CT, Kirst M, Albert PS, Birchler JA, Bradbury PJ, et al. Aluminum tolerance in maize is associated with higher MATE1 gene copy number. Proc Nat Acad Sci. 2013;110:5241–6. doi:10.1073/pnas.1220766110.PubMed CentralView ArticlePubMedGoogle Scholar
- McHale LK, Haun WJ, Xu WW, Bhaskar PB, Anderson JE, Hyten DL et al. Structural variants in the soybean genome localize to clusters of biotic stress-response genes. Plant Physiol. 2012;159:1295–1308. doi:10.1104/pp.112.194605
- Muñoz-Amatriaín M, Eichten SR, Wicker T, Richmond TA, Mascher M, Steuernagel B, et al. Distribution, functional impact, and origin mechanisms of copy number variation in the barley genome. Genome Biol. 2013;14:R58. doi:10.1186/gb-2013-14-6-r58.PubMed CentralView ArticlePubMedGoogle Scholar
- Von Arnim AG, Jia Q, Vaughn JN. Regulation of plant translation by upstream open reading frames. Plant Sci. 2014;214:1–12. doi:10.1016/j.plantsci.2013.09.006.View ArticleGoogle Scholar
- Lohmer S, Maddaloni M, Motto M, Salamini F, Thompson RD. Translation of the mRNA of the maize transcriptional activator Opaque-2 is inhibited by upstream open reading frames present in the leader sequence. Plant Cell. 1992;5:65–73. doi:10.1105/tpc.5.1.65.View ArticleGoogle Scholar
- Locatelli F, Magnani E, Vighi C, Lanzanova C, Coraggio I. Inhibitory effect of myb7 uORF on downstream gene expression in homologous (rice) and heterologous (tobacco) systems. Plant Mol Biol. 2002;48:309–18. doi:10.1023/A:1013340004348.View ArticlePubMedGoogle Scholar
- Hanfrey C, Francescetti M, Mayer MJ, Illingworth C, Michael AJ. Abrogation of upstream open reading frame-mediated translational control of a plant S-adenosylmethionine decarboxylase results in polyamine disruption and growth perturbations. J Biol Chem. 2002;277:44131–9. doi:10.1074/jbc.M206161200.View ArticlePubMedGoogle Scholar
- Emebiri L, Michael P, Moody D. Enhanced tolerance to boron toxicity in two-rowed barley by marker-assisted introgression of favourable alleles derived from Sahara 3771. Plant Soil. 2009;314:77–85. doi:10.1007/s11104-008-9707-0.View ArticleGoogle Scholar
- Schnurbusch T, Hayes J, Sutton T. Boron toxicity tolerance in wheat and barley: Australian perspectives. Breeding Sci. 2010;60:297–304. doi:10.1270/jsbbs.60.297.View ArticleGoogle Scholar
- Coventry SJ, Smith DS, McKay A, Eckermann PJ, Eglinton JK. Investigating root architecture in barley and responses to salinity and high boron. In: Ceccarelli S, Grando S, editors. Proc 10th Int Barley Genet Symp. Alexandria, Egypt: ICARDA; 2008. pp. 180–189.Google Scholar
- Moody D. ‘SLOOP VIC’ Application No. 2002/066. Plant Varieties J. 2003;16:288–92.Google Scholar
- Burton RA, Shirley NJ, King BJ, Harvey AJ, Fincher GB. The CesA gene family of barley. Quantitative analysis of transcripts reveals two groups of co-expressed genes. Plant Physiol. 2004;134:224–36. doi:10.1104/pp. 103.032904.PubMed CentralView ArticlePubMedGoogle Scholar
- Hayes JE, Pallotta M, Baumann U, Berger B, Langridge P, Sutton T. Germanium as a tool to dissect boron toxicity effects in barley and wheat. Functional Plant Biol. 2013;40:618–27. doi:10.1071/FP12329.View ArticleGoogle Scholar
- Wheal MS, Fowles TO, Palmer LT. A cost-effective acid digestion method using closed polypropylene tubes for inductively coupled plasma optical emission spectrometry (ICP-OES) analysis of plant essential elements. Anal Methods. 2011;3:2854–63. doi:10.1039/C1AY05430A.View ArticleGoogle Scholar
- Wolf B. Improvements in the azomethine-H method for the determination of boron. Comm Soil Sci Plant Anal. 1974;5:39–44. doi:10.1080/00103627409366478.View ArticleGoogle Scholar
- Pallotta MA, Warner P, Fox RL, Kuchel H, Jefferies SJ, Langridge P. Marker assisted wheat breeding in the southern region of Australia. In: Pogna NE, editor. Proc 10th Int Wheat Genet Symp. Paestum, Italy: Istituto Sperimentale per la Cerealicoltura; 2003. p. 789–91.Google Scholar
- Franckowiak JD, Lundqvist U. Coordinator’s report: Rules for nomenclature and gene symbolization in barley. Barley Genet Newsletter. 2010;40:178–82.Google Scholar