Genetic architecture of limit dextrinase inhibitor (LDI) activity in Tibetan wild barley
© Huang et al.; licensee BioMed Central Ltd. 2014
Received: 27 January 2014
Accepted: 28 April 2014
Published: 1 May 2014
Limit dextrinase inhibitor (LDI) inhibits starch degradation in barley grains during malting because it binds with limit dextrinase (LD). There is a wide genetic variation in LDI synthesis and inactivation during barley grain development and germination. However, the genetic control of LDI activity remains little understood.
In this study, association analysis was performed on 162 Tibetan wild accessions by using LDI activity, 835 Diversity Arrays Technology (DArT) markers and single nucleotide polymorphisms (SNPs) of the gene HvLDI encoding LDI. Two DArT markers, bpb-8347, bpb-0068, and 31 SNPs of HvLDI were significantly associated with LDI activity, explaining 10.0%, 6.6% and 13.4% of phenotypic variation, respectively. Bpb-8347 is located on chromosome 6H, near the locus of HvLDI, and bpb-0068 is located on 3H.
The current results confirmed the locus of the gene controlling LDI activity and identified a new DArT markers associated with LDI activity. The SNPs associated with LDI activity may provide a new insight into the genetic variation of LDI activity in barley grains.
Malting is the first step in brewing or beer production of malt barley, and the properties of malting barley are generally evaluated on the basis of malting quality. High malting quality is characterized by high malt extract (ME) and diastatic power (DP), as well as high wort fermentation by yeast, i.e. complete degradation of starch and polysaccharide in malt . The major chemical compound in barley grains is starch, which consists of 30% amylose and 70% amylopectin .
Starch is degraded by coordinative action of α-amylase, β-amylase, limit dextrinase (LD) and α-glucosidase in the endosperm of the germinating grains . LD is the only enzyme of cleaving α-1-6 linkages in branched dextrins molecules . Without the action of LD, the branched dextrins could not be fermented by yeast . Hence, LD plays an important role in the brewing and distilling industries for producing fermentable sugar.
LD in barley grains exists in three different forms: insoluble bound, soluble inactive and active free, and only the active free form participates in starch mobilization . However, the free LD accounts for a small portion of total LD. The low activity of available LD in grains is due to its combination with endogenous inhibitors, limit dextrinase inhibitor (LDI) . LD and LDI interact with each other in a 1:1 molar ratio complex . The LD bound to LDI is thought to be a limiting factor for complete degradation of starch .
LDI is synthesized in the developing grains, and degraded gradually during the malting process . Consequently, LD activity in malt is not only dependent on the amount of the enzyme itself in mature grains, and also is closely related to the amount and type of LDI. HvLDI was identified to encode the LDI protein, and located at chromosome 6H [10, 11]. However, genetic variation and controlling of LDI synthesis and degradation in barley grains and malt is not totally understood. In view of a wider variation of LD activity in Tibetan wild barley than in cultivated barley (unpublished data), it is imperative to determine the genetic variation of LDI content in Tibetan wild barley, which has been recently proved to be a progenitor of cultivated barley and is rich in genetic diversity [12, 13].
Association analysis has been widely used to identify genes or loci of important traits, such as abiotic stress tolerance [14, 12, 15], agronomic traits [16, 17] and also barley quality [18, 19]. However, there has been no report on the discovery of novel loci or elite alleles attributed to LDI activity using association analysis.
The objectives of this study were to (1) determine the genetic variation of LDI activity in Tibetan wild barley; (2) identify the DArT markers and SNPs of HvLDI associated with LDI activity in Tibetan wild barley.
A total of 162 Tibetan wild barley accessions, kindly provided by professor Sun of Huazhong Agricultural University, China, were used for LDI activity and association analysis. All accessions were planted in early November 2011 and 2012 in Zijingang campus, Zhejiang University (Hangzhou, China). Every genotype was cultivated in accordance with local agronomic practices with three replications. The 162 barley accessions were planted in a block with each accession consisting of three lines (2 m length and line distance is 0.25 m).
Assay of LDI activity
Grain samples were micro-malted in a Phoenix System Micro-malting Apparatus (Adelaide, Australia) in the order of steeping, germination and drying. LD was partially purified from barley malt according to Kristensen et al. . The fractions were applied to gel filtration chromatography after the ion exchange chromatography step, and partially purified LD was collected and used for the measurement of LDI activity.
LDI activity was determined according to MacGregor et al.  with some modification. One ml of 0.1 M sodium acetate (pH 5.5) containing 10 mM 1,10-phenanthroline was added to 0.1 g barley powder, incubated at 4°C for 30 min. The extract was heated at 70°C for 40 min, centrifuged and the supernatant was collected. The protein content of the extracts was measured using Bradford assay kits (Sangon Biotech). Twenty micrograms protein of LDI extract were mixed with 10 mU partially purified LD and the volume was made up to 0.5 ml in 0.1 M maleic acid containing 0.02% Na azide (pH 5.5). The remaining LD was determined using the Limit-Dextrinase assay kit (Megazyme). The LDI activity was calculated as the reduced LD activity.
Population structure and kinship analysis
Totally, 771 DArT markers, with minor allele frequency (MAF) higher than 0.03 (Additional file 1: Table S1), were used for population structure analysis using STRUCTURE software (v2.3.3) , in which the number of clusters (k) was set from 1 to 12 and ten iterations were performed in an admixture model with 10,000 burning period and 100,000 MCMC (Markov Chain Monte Carlo). The DArT markers used were derived from Diversity Arrays Technology Pty Ltd, Australia and distributed over the whole genome [23–25]. The most probable number of clusters (k) was estimated according to the value of Δk. When Δk had the highest value, the value of k was the number of clusters . Kinship was estimated using SPAGeDi software . Genetic distance and neighbor-joining tree were developed with NTSYSpc (version 2.10e).
Variance component and heritability estimation
Variance analysis was conducted using SPSS software. The used model was: y = mu + GENOTYPE + ERROR (fixed in low case, random in capitals). Heritability was estimated according to H2 = Vg/(Vg + Ve).
Genome-wide association study (GWAS) of LDI activity
Genome-wide association study between LDI activity (mean value of 2011 and 2012) and DArT markers was conducted using TASSEL software (v3.0), where Q, K and Q + K methods were applied . The structure matrix was included as covariate to correct population structure in Q model: y = Qν + Xβ + e, where X is the DArT marker matrix, Q is the structure matrix, β and ν are coefficient vectors for DArT marker vector and population structure vector, and e is the random error vector. In K model: y = Xβ + Zu + e. The Q + K model can be written in a matrix form as: y = Qν + Xβ + Zu + e, where X, β, Q and ν are the same as those mentioned above, Z is the kinship matrix, and u is a vector of random genetic effects u ~ N (0, 2 K), where K is the kinship matrix. Manhattan plots were displayed using R software (v2.14.2).
PCR amplification and sequencing
Genomic DNA was extracted from leaves of 150 barley seedlings using CTAB method . The primers used for DNA amplification were designed using the primer design tool of NCBI (http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/tools/primer-blast/). Primers for HvLDI were: forward, TTTTCGCATGTCACCAAAAATGT; reverse, TCCGCTTCATTACCTTGGCG. The amplified DNA fragment was about 900 bp.
The PCR reactions were completed as follows: 2.5 μl of 10 × TransTaq HiFi buffer I (200 mM Tris–HCl (pH 8.4), 100 mM (NH4)2SO4, 20 mM MgCl2, 200 mM KCl, 2 μl of 2.5 mM dNTPs, 1 μl of 10 μM forward primers, 1 μl of 10 μM reverse primer, 0.5 μl of 5 units μl-1 of TransTaq polymerase High Fidelity (Beijing TransGen Biotech Co., Ltd.), and 1 μl of 50 ng of genomic DNA. The PCR amplification program started at 5 min at 95°C, followed by 32 cycles of 30 s at 95°C, 30 s at 60°C, and 1 min at 72°C, and then 10 min at 72°C for a final extension. The PCR products were analyzed by 1% agarose gel electrophoresis. They were sequenced using 3730 × l DNA Analyzer (Applied Biosystems Inc., USA).
Polymorphism of HvLDI and haplotype analysis
Availability of supporting data
The sequences of the 6 haplotypes were deposited in GenBank, and the accession numbers of these haplotypes, H1-H6 are KJ710426, KJ710427, KJ710427, KJ710429, KJ710430 and KJ710431, respectively.
LDI activities of 162 Tibetan wild accessions
Variance component and heritability (H 2 ) estimates for LDI activity of 2011 and 2012
Genome-wide association study
Bpb-8347 and bpb-0068, located at 73.44 cM on chromosome 6H and 66.50 cM on 3H, showed significant associations with LDI activity, explaining 6.8% and 5.3% of phenotypic variation in 2011 and 10.4% and 5.5% in 2012 under Q + k method. Similarly, the two DArT markers could explain 10.0% and 6.6% of phenotypic variation for the mean values of 2011 and 2012.
Association between HvLDI and LDI activity and haplotype analysis
In this study, we investigated the genetic variation of LDI activity and performed GWAS analysis of HvLDI, in order to identify the loci controlling LDI activity in Tibetan wild barley. A wide variation in LDI activity was found in the 162 Tibetan wild barley accessions, and some of which displayed a minus value. It should be noted that some dextrins released by the de-branching enzymes during germination could also be detected by the Limit-Dextrizyme method, and these dextrins may account for the overestimation of LD activity [11, 32, 33]. Thus, the detected LD activities in extraction solutions of some accessions were higher than those of added LD, resulting in minus values of LDI activity, which was calculated as the reduced LD activity. However, the minus values of LDI in this experiment only occurred in very few accessions, and were very small (close to zero). Therefore the accessions with minus values of LDI are quite low in LDI activity. Actually, the impact of these dextrins produced during short-time extraction on the final LDI values could be negligible, as reported by .
Two DArT markers, i.e. bpb-8347 and bpb-0068 were found to be associated with LDI activities in Tibetan wild barley. Bpb-8347, a marker with highest GWA peak, is located near the locus of HvLDI gene obtained from IPK database, where the whole barley genome sequences lay . This result confirms that HvLDI is a major gene controlling LDI activity in barley grain. Bpb-0068 might be a potential locus contributing to LDI activity, and it has not been reported up to date. Hence, it may be assumed that some minor effect genes affecting LDI activity may exist in Tibetan wild barley, and need to be confirmed by QTL analysis in segregating population in future work.
A total of 31 SNPs of HvLDI were significantly associated with LDI activity. It is also found that 11 of these SNPs are significantly associated with the malt extract in 56 cultivated genotypes in a previous study , in which association analysis between HvLDI gene and LDI activity was not conducted. These results indicate that the LDI activity has a close relationship with malt extract, because the two traits are affected by the same 11 SNPs. Indeed, it was previously reported that the HvLDI antisense lines showed unpredicted pleiotropic effects on numerous enzyme activities, including α- and β-amylases, starch synthases, as well as starch granule types, all of which affect malt extract . In addition, the LD activity-associated SNP in the previous study  was not detected in the Tibetan wild population, suggesting that LD activity in malt may be independent on polymorphism of HvLDI gene in Tibetan population.
Interestingly, HA and HB showed a significant difference of the LDI activity only in Q4 subpopulation (Figure 6b,c). Q4 contains all kinds of haplotypes, while every other subpopulation mainly consists of one haplotype (Figure 6a). In other word, Q4 has a wider genetic diversity than every other subpopulation. It seems that these associated SNPs are false positive due to highly structured LD (Additional file 5: Figure S2) in this chromosomal region, and similar phenomenon were also observed in association analysis of frost tolerance [36, 37]. It may be hypothesized that the variation of LDI activity in Q4 may be caused by other variations, which share linkage disequilibrium with the detected SNPs, in or around this region [34, 38]. These variations should only belong to the accessions of HA in Q4 subpopulation, and may exist in the upstream or promoter region of HvLDI, or the regulated factors, which co-segregated with the detected SNPs. However, the assumption needs to be confirmed in the further work.
Four SNPs (13, 24, 72, and 352) within a coding region generate the alteration of the amino acids sequence. But these two kinds of haplotypes showed no significant difference in LDI activity among all accessions. It can be concluded that the alteration of LDI activity may be not due to the change of amino acid composition.
This study confirmed the locus of LDI gene and detected a DArT marker, bpb-0068, associated with LDI activity. The results provide useful information for identifying minor effect genes affecting LDI activity, and also prove that Tibetan wild barley is an elite germplasm, which may act as an abundant gene pool for barley breeding. Furthermore, the detected SNPs in this study should be helpful for better understanding the genetic control of LDI activity in barley grain.
This research was supported by National Natural Science Foundation of China (Nos. 31330055, 31271656, 31129005), Zhejiang Provincial Natural Science Foundation of China (No. Z3110054) and China Agriculture Research System. We thank Mr Xianyin Zhang and Ms. Mei Li, the technicians of 985-Institute of Agrobiology for providing convenience in using the experimental equipments.
- Fincher GB: Molecular and cellular biology associated with endosperm mobilisation in germinating cereal grains. Annu Rev Plant Physiol Plant Mol Biol. 1989, 40: 305-346. 10.1146/annurev.pp.40.060189.001513.View ArticleGoogle Scholar
- Vasanthan T, Hoover R: Barley starch: production, properties, modification and uses. Starch: Chemistry and Technology. 3rd edition. Edited by: BeMiller J, Whistler R. New York: Academic press, 2009:601-628.View ArticleGoogle Scholar
- Manners DJ: Starch degradation during malting mashing. Brewers Digest. 1974, 49: 56-62.Google Scholar
- Manners DJ, Marshall JJ, Yellowlees D: The specificity of cereal limit dextrinase. Biochem J. 1970, 116: 539-541.PubMed CentralView ArticlePubMedGoogle Scholar
- Panchal CJ, Stewart GG: Utilisation of wort carbohydrates. Brewers Digest. 1979, 54: 36-43.Google Scholar
- Sissons MJ, Lance RCM, Wallace W: Bound and free forms of barley limit dextrinase. Cereal Chem. 1994, 71: 520-521.Google Scholar
- Macri LJ, MacGregor AW, Shroeder SW, Bazin SL: Detection of a limit dextrinase inhibitor in barley. J Cereal Sci. 1993, 18: 103-106. 10.1006/jcrs.1993.1038.View ArticleGoogle Scholar
- MacGregor AW, Donald LJ, MacGregor EA, Duckworth HW: Stoichiometry of the complex formed by barley limit dextrinase with its endogenous inhibitor. Determination by electrospray time-of-flight mass spectrometry. J Cereal Sci. 2003, 37: 357-362. 10.1006/jcrs.2002.0500.View ArticleGoogle Scholar
- Sissons M, Taylor M, Proudlove M: Barley malt limit dextrinase: its extraction, heat stability and activity during malting and mashing. J Am Soc Brew Chem. 1995, 53: 104-110.Google Scholar
- Karakousis A, Barr AR, Kretschmer JM, Manning S, Logue SJ, Roumeliotis S, Collins HM, Chalmers KJ, Li CD, Lance RCM, Langridge P: Mapping and QTL analysis of the barley population Galleon × Haruna Nijo. Aust J Agr Res. 2003, 54 (12): 1131-1135. 10.1071/AR02179.View ArticleGoogle Scholar
- Stahl Y, Alexander RD, Coates S, Bryce JH, Jenkinson HR, Morris PC: The barley limit dextrinase inhibitor: Gene expression, protein location and interaction with 14-3-3 protein. Plant Sci. 2007, 172: 452-461. 10.1016/j.plantsci.2006.10.008.View ArticleGoogle Scholar
- Wu DZ, Qiu L, Xu LL, Ye LZ, Chen MX, Sun DF, Chen ZH, Zhang HT, Jin XL, Dai F, Zhang GP: Genetic Variation of HvCBF Genes and Their Association with Salinity Tolerance in Tibetan Annual Wild Barley. PLoS One. 2011, 6 (7): e22938-10.1371/journal.pone.0022938.PubMed CentralView ArticlePubMedGoogle Scholar
- Dai F, Nevo E, Wu D, Comadran J, Zhou M, Qiu L, Chen ZH: Avigdor Beiles, Chen GX, Zhang GP: Tibet is one of the centers of domestication of cultivated barley. Proc Natl Acad Sci U S A. 2012, 109 (42): 16969-16973. 10.1073/pnas.1215265109.PubMed CentralView ArticlePubMedGoogle Scholar
- Nguyen Viet L, Oene D, Marcos M, Benjamin K, Andreas G, Richard GF, Visser C, Gerard Van Der L: Association mapping of salt tolerance in barley (Hordeum vulgare L.). Theor Appl Genet. 2013, doi:10.1007/s00122-013-2139-0Google Scholar
- Long NV, Dolstra O, Malosetti M, Kilian B, Graner A, Visser RG, Van Der Linden CG: Association mapping of salt tolerance in barley (Hordeum vulgare L.). Theor Appl Genet. 2013, doi:10.1007/s00122-013-2139-0Google Scholar
- Neumann K, Kobiljski B, Dencic S, Varshney RK, Borner A: Genome-wide association mapping: a case study in bread wheat (Triticum aestivum L.). Mol Breeding. 2011, 27: 37-58. 10.1007/s11032-010-9411-7.View ArticleGoogle Scholar
- Pasam RK, Sharma R, Malosetti M, Van Eeuwijk FA, Haseneyer G, Kilian B, Graner A: Genome-wide association studies for agronomical traits in a world wide spring barley collection. BMC Plant Biol. 2012, 12: 16-10.1186/1471-2229-12-16.PubMed CentralView ArticlePubMedGoogle Scholar
- Cai SG, Yu G, Chen XH, Huang YC, Jiang XG, Zhang GP, Jin XL: Grain protein content variation and its association analysis in barley. BMC Plant Biol. 2013, 13: 35-10.1186/1471-2229-13-35.PubMed CentralView ArticlePubMedGoogle Scholar
- Jin XL, Cai SG, Ye LZ, Chen ZH, Zhou MX, Zhang GP: Association of HvLDI with limit dextrinase activity and malt quality in barley. Biotechnol Lett. 2013, 35: 639-645. 10.1007/s10529-012-1106-1.View ArticlePubMedGoogle Scholar
- Kristensen M, Planchot V, Abe J, Svensson B: Large-scale purification and characterization of barley limit dextrinase, a member of the alpha amylase structural family. Cereal Chem. 1998, 75: 473-479. 10.1094/CCHEM.1922.214.171.1243.View ArticleGoogle Scholar
- MacGregor AW, Macri LJ, Schroeder SW, Bazin SL: Purification and characterization of limit dextrinase inhibitors from barley. J Cereal Sci. 1994, 20: 33-41. 10.1006/jcrs.1994.1042.View ArticleGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.PubMed CentralPubMedGoogle Scholar
- Wenzl P, Li H, Carling J, Zhou M, Raman H, Paul E, Hearnden P, Maier C, Xia L, Caig V, Jaroslava O, Mehmet C, David P, Junping W, Rosy R, Smith KP, Muehlbauer GJ, Chalmers KJ, Andris K, Eric H, Andrzej K: A high-density consensus map of barley linking DArT markers to SSR, RFLP and STS loci and agricultural traits. BMC Genomics. 2006, 7 (1): 206-10.1186/1471-2164-7-206.PubMed CentralView ArticlePubMedGoogle Scholar
- Wenzl P, Carling J, Kudrna D, Jaccoud D, Huttner E, Kleinhofs A, Kilian A: Diversity Arrays Technology (DArT) for whole-genome profiling of barley. Proc Natl Acad Sci U S A. 2004, 15 (26): 9915-9920. doi:10.1073/pnas.0401076101View ArticleGoogle Scholar
- Alsop BP, Farre A, Wenzl P, Wang JM, Zhou MX, Romagosa I, Kilian A, Steffenson BJ: Development of wild barley-derived DArT markers and their integration into a barley consensus map. Mol Breeding. 2011, 27: 77-92. 10.1007/s11032-010-9415-3. doi:10.1007/s11032-010-9415-3View ArticleGoogle Scholar
- Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005, 14: 2611-2620. 10.1111/j.1365-294X.2005.02553.x.View ArticlePubMedGoogle Scholar
- Hardy OJ, Vekemans X: Spagedi: a versatile computer program to analyze spatial genetic structure at the individual or population levels. Mol Ecol Notes. 2002, 2: 618-620. 10.1046/j.1471-8286.2002.00305.x.View ArticleGoogle Scholar
- Bradbury PJ, Zhang ZW, Kroon DE, Casstevens TM, Ramdoss Y, Buckler ES: TASSEL: software for association mapping of complex traits in diverse samples. Bioinformatics. 2007, 23 (19): 2633-2635. 10.1093/bioinformatics/btm308.View ArticlePubMedGoogle Scholar
- Huang JC, Ge XJ, Sun M: Modified CTAB protocol using a silica matrix for isolation of plant genomic DNA. Biotechniques. 2000, 28: 432-434.PubMedGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.PubMed CentralView ArticlePubMedGoogle Scholar
- Librado P, Rozas J: DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinform Appl Note. 2009, 25: 1451-1452. 10.1093/bioinformatics/btp187.View ArticleGoogle Scholar
- MacGregor AW, Bazin SL, Schroedor SW: Effect of starch hydrolysis products on the determination of limit dextrinase and limit dextrinase inhibitors in barley and malt. J Cereal Sci. 2002, 35: 17-28. 10.1006/jcrs.2001.0408.View ArticleGoogle Scholar
- McDougall GJ, Ross HA, Swanston JS, Davies HV: Limit dextrinase from germinating barley has endotransglycosylase activity, which explains its activation by maltodextrins. Planta. 2004, 218: 542-551. 10.1007/s00425-003-1141-1.View ArticlePubMedGoogle Scholar
- Mayer KF, Waugh R, Brown JW, Schulman A, Langridge P, Platzer M, Fincher GB, Muehlbauer GJ, Sato K, Close TJ, Wise RP, Stein N, International Barley Genome Sequencing Consortium: A physical, genetic and functional sequence assembly of the barley genome. Nature. 2012, 491: 711-717.PubMedGoogle Scholar
- Stahl Y, Coates S, Bryce JH, Morris PC, Morris PC: Antisense downregulation of the barley limit dextrinase inhibitor modulates starch granule size distribution, starch composition and amylopectin structure. Plant J. 2004, 39: 599-611. 10.1111/j.1365-313X.2004.02159.x.View ArticlePubMedGoogle Scholar
- Fricano A, Rizza F, Faccioli P, Pagani D, Pavan P, Stella A, Rossini L, Piffanelli P, Cattivelli L: Genetic variants of HvCbf14 are statistically associated with frost tolerance in a European germplasm collection of Hordeum vulgare. Theor Appl Genet. 2009, 119: 1335-1348. 10.1007/s00122-009-1138-7.PubMed CentralView ArticlePubMedGoogle Scholar
- Visioni A, Tondelli A, Francia E, Pswarayi A, Malosetti M, Russell J, Thomas W, Waugh R, Pecchioni N, Romagosa I, Jordi C: Genome-wide association mapping of frost tolerance in barley (Hordeum vulgare L.). BMC Genomics. 2013, 14: 424-10.1186/1471-2164-14-424.PubMed CentralView ArticlePubMedGoogle Scholar
- Myles S, Davison D, Barrett J, Stoneking M, Timpson N: World wide population differentiation at disease associated SNPs. BMC Genomics. 2008, 1: 22-10.1186/1755-8794-1-22.View ArticleGoogle Scholar
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 (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.