中国青藏高原青稞抗穗发芽资源评价与α-淀粉酶基因的克隆及表达

穗发芽（PHS，preharvest sprouting）是影响禾本科作物生产的重要的灾害之一。收获时期如遇潮湿天气容易导致穗发芽发生。发生穗发芽的种子内部水解酶（主要是α-淀粉酶）活性急剧升高，胚乳贮藏物质开始降解，造成作物产量和品质严重降低。因此，选育低穗发芽风险的品种是当前作物育种工作中面临的重要任务。 青稞(Hordeum vulgare ssp. vulgare)主要分布于青藏高原，自古以来就是青藏高原人民的主要粮食。近年来，由于青稞丰富的营养成分和特有的保健品质、在燃料工业中的潜力以及在啤酒酿造工业中的利用前景，在发达国家日趋受到重视，掀起综合研究利用的热潮。我国拥有占全世界2/3 以上的青稞资源，具有发展青稞产业的得天独厚的条件。然而，由于青稞收获期间恰逢青藏高原雨季来临，常有穗发芽灾害发生，使青稞生产损失巨大。目前对青稞穗发芽研究很少，适用于育种的穗发芽抗性材料相对缺乏，不能很好的满足青稞穗发芽抗性育种的需要。本研究以青藏高原青稞为材料，对其穗发芽抗性的评价指标和体系进行构建，同时筛选青稞抗穗发芽品种并对其抗性进行评价，还利用分子生物学手段对青稞穗发芽抗性的分子机理进行了初步探讨。主要研究结果如下： 1. 本试验以来自于我国青藏高原地区的青稞为材料，对休眠性测定的温度范围进行探讨，并对各种穗发芽抗性测定方法的对青稞的适用性进行评测。通过探讨温度对13 个不同基因型的青稞籽粒发芽和休眠性表达的影响，对筛选青稞抗穗发芽资源的温度条件进行探索，并初步分析了其休眠性表达的机理。在10，15，20，25，30℃的黑暗条件下，选用新收获的13 个青稞品种为材料进行籽粒发芽实验，以发芽指数（GI）评价其休眠性。结果发现，不同品种对温度敏感性不同，其中温度不敏感品种，在各温度条件下均表现很低的休眠性；而温度敏感品种，其休眠性表达受低温抑制，受高温诱导。15℃至25℃是进行青稞休眠性鉴定的较适宜的温度范围。通过对供试材料发芽后的α-淀粉酶活性，发现温度对青稞种子的休眠性表达的影响至少在一定程度上表现在对α-淀粉酶活性的调控上。随后，对分别在马尔康和成都进行种植的34 份青稞穗发芽指数（SI），穗发芽率（SR），籽粒发芽指数（GI）和α-淀粉酶活性（AA）进行了测定和分析，发现它们均受基因型×栽培地点的极显著影响，且四个参数之间具有一定相关性。GI 参数由于其变异系数较低，在不同栽培地点稳定性好，且操作简便，是较可靠和理想的穗发芽评价参数。SI 参数可作为辅助，区别籽粒休眠性相似的材料（基因型）或全面评价材料（基因型）的穗发芽抗性特征。AA 参数稳定性较差，并且检测方法复杂，因此不建议在育种及大量材料筛选和评价时使用。此外，青稞穗发芽抗性受环境影响较大，评价时应考虑到尽可能多的抗性影响因素及其在不同栽培条件下的变异。 2. 对来自青藏高原的青稞穗发芽抗性特征及其与其它农艺性状间的关系进行研究。通过测定穗发芽指数（SI）、籽粒发芽指数（GI）和α-淀粉酶活性（AA），表明113 份青稞材料的穗发芽抗性具有显著差异。SI、GI 和AA 参数的变幅分别为1.00~8.86、0.01~0.97 和0.00~2.76，其均值分别为4.72、0.63 和1.22。根据SI 参数，六个基因型，包括‘XQ9-5’，‘XQ33-9’，‘XQ37-5’，‘XQ42-9’，‘XQ45-7’和‘JCL’被鉴定为抗性品种。综合SI、GI 和AA 参数，可以发现青稞的穗发芽抗性机制包含颖壳等穗部结构的抗性和种子自身的抗性（即种子休眠性），且供试材料中未发现较强的胚休眠品种，除‘XQ45-7’外，所有品种在发芽第四天均能检测出α-淀粉酶活性。穗部结构和种子休眠的抗性机制因基因型不同而不同，在穗发芽抗性中可单独作用或共同作用。农家品种和西藏群体分别比栽培品种和四川群体的穗发芽抗性强，而在不同籽粒颜色的青稞中未发现明显差异。相关性检验发现，青稞的穗发芽抗性，主要是种子休眠性，与百粒重、开花期、成熟期、穗长、芒长和剑叶长呈显著负相关关系，与株高相关性不显著。农艺性状可以作为穗发芽抗性材料选育中的辅助指标。本试验为青稞穗发芽抗性育种研究提供了必要的理论基础和可供使用的亲本材料。 3. α-淀粉酶是由多基因家族编码的蛋白质，在植物种子萌发时高度表达，与植物种子的萌发能力密切相关。在大麦种子发芽时，高等电点α-淀粉酶的活性远大于低等电点的α-淀粉酶。为了研究不同穗发芽抗性青稞品种中编码高等电点α-淀粉酶Amy1 基因结构与抗性间的关系，我们以筛选得到的抗性品种‘XQ32-5’（TR1）、‘XQ37-5’（TR2）、‘XQ45-7’（TR3），易感品种‘97-15’（TS1）、‘9657’（TS2）以及强休眠大麦品种‘SAMSON’（SAM）为材料，对其Amy1 基因的编码区序列进行克隆和结构分析，并对它们推导的氨基酸序列进行比较。结果显示，青稞Amy1 基因具有三个外显子、两个内含子，编码区中有13 个核苷酸变异位点，均位于2、3 号外显子，2 个变异位点位于2 号外显子。SAM 和TS1 分别在2 号外显子相应位置有5 个相同的碱基(GAACT)的插入片段。相应α-淀粉酶氨基酸序列推导发现，所有核苷酸变异中有8 个导致相应氨基酸残基的改变，其余位点为同义突变。青稞Amy1 基因编码区序列品种间相似度高达99%以上，部分序列变异可能与其穗发芽抗性有关。随后，我们又通过SYBR Green 荧光定量技术对该基因在不同发芽时间（1d~7d）的相对表达水平进行了差异性检测。结果发现，7 天内不能检测到SAM 的Amy1 基因表达，5 个青稞品种间的Amy1 基因的相对表达量均随着发芽时间延长而上升，但上升方式有所不同。弱抗品种该基因表达更早，转录本增加速率更大，且在4~5 天可达到平台期。发芽7 天中，抗性品种总转录水平明显低于易感品种。本研究结果表明，青稞Amy1 基因的转录水平是与其穗发芽抗性高度相关。 我国青藏高原青稞，尤其是农家品种的穗发芽抗性具有丰富的变异，蕴藏着穗发芽抗性育种的宝贵资源。本研究为青稞穗发芽抗性育种建立了合理抗性评价体系，筛选出可供育种使用的特殊材料，阐明了农艺性状可辅助穗发芽抗性育种，同时还对穗发芽抗性与α-淀粉酶基因的结构和表达关系进行分析，为青稞穗发芽抗性资源筛选奠定了基础。 Preharvest sprouting (PHS) is a serious problem in crop production. It often takes place when encountering damp, cold conditions at harvest time and results in the decrease of grain quality and great loss of yield by triggering the synthesis of endosperm degrading enzymes (mostly the α-amylase). Therefore, PHS is regarded as an important criterion for crop breeding. In order to minimize the risk of PHS, resistant genotypes are highly required. Hulless barley (Hordeum vulgare ssp. vulgare) is the staple food crop in Qinghai-Tibetan Plateau from of old, where is one of the origin and genetic diversity centers of hulless barley. Recently, interest in hulless barley has been sparked throughout the world due to the demonstrations of its great potential in health food industry and fuel alcohol production. Indeed, hulless barley can also be utilized to produce good quality malt if the appropriate malting conditions are used. In China, overcast and rainy conditions often occur at maturity of hulless barley and cause an adverse on its production and application. PHS resistant genotypes, therefore, are highly required for the hulless barley breeding programs. However, few investigations have been made so far on this issue. The objectives of this study were: 1) to assessment of methods used in testing preharvest sprouting resistance in hulless barley; 2) to evaluate the variability and characteristics of PHS resistance of hulless barley from Qinghai-Tibet Plateau in China; 3) to select potential parents for PHS resistance breeding; 4) to primarily study on the molecular mechanism of PHS resistance of hulless barley. Our results are as followed: 1. We investigated the temperature effects on seed germination and seed dormancy expression of hulless barley, discussed appropriate temperature range for screening of PHS resistant varieties, and analyzed the mechanism of seed dormancy expression of hulless barley. The dormancy level of 13 hulless barley were evaluated by GI (germination index) values calculating by seed germination tests at temperature of 10，15，20，25，30℃ in darkness. There were great differences in temperature sensitivity among these accessions. The insensitive accessions showed low dormancy at any temperature while the dormancy expression of sensitive accessions could be restrained by low temperature and induced by high temperature. The temperature range of 15℃ to 25℃ was workable for estimating of dormancy level of hulless barley according to our data. Analysis of α-amylase activity showed that the temperature effects on seed germination and the expression of seed dormancy be achieved probable via regulating of α-amylase activity. Furthermore, we evaluated the differences in sprouting index (SI), sprouting rate (SR), germination index (GI) and α-amylase activity (AA) between Maerkang and Chengdu among 34 accessions of hulless barley from Qinghai-Tibetan Plateau in China. These PHS sprouting parameters were significantly affected by accession×location, and they had correlation between each other. GI was the most reliable parameter because of its low CV value, good repeatability and simple operation. SI could assist in differentiating between accessions of similar dormancy or overall evaluation of the resistance. AA was bad in repeatability and had relatively complex testing method, therefore, not appropriate for breeding and evaluation and screening of PHS resistant materials. Besides, since PHS resistance of hulless barley was greatly influenced by its growth environment, possibly much influencing factors and variations between cultivated conditions should be considered. 2. In this study, large variation was found among 113 genotypes of hulless barley (Hordeum vulgare ssp.vulgare) from Qinghai-Tibetan Plateau in China, based on the sprouting index (SI), germination index (GI) and α-amylase activity (AA) which derived from sprouting test of intact spikes, germination test of threshed seeds and determination of α-amylase activity, respectively. The range of SI, GI and AA was 1.00~8.86, 0.01~0.97 and 0.00~2.76，the mean was 4.72, 0.63 and 1.22 espectively. Six resistant genotypes, including ‘XQ9-5’, ‘XQ33-9’, ‘XQ37-5’, ‘XQ42-9’, ‘XQ45-7’ and ‘JCL’, were identified based on SI. Integrating the three parameters, it was clear that both hulls and seeds involved in PHS resistance in intact spikes of hulless barley and there was no long-existent embryo dormancy found among the test genotypes. All the genotypes, except ‘XQ45-7’, had detectable α-amylase activity on the 4th day after germination. There was PHS resistance imposed by the hull and seed per se and the two factors can act together or independent of each other. Besides, landraces or Tibet hulless barley had a wider variation and relatively more PHS resistance when compared with cultivars or Sichuan hulless barley. No significant difference was found among hulless barley of different seed colors. The correlation analysis showed PHS resistance was negatively related to hundred grain weight, days to flowering, days to maturity, spike length, awn length and flag length but not related to plant height. This study provides essential information and several donor parents for breeding of resistance to PHS. 3. Alpha-amylase isozymes are encoded by a family of multigenes. They highly express in germinating seeds and is closely related to seed germination ability. In barley germinating seeds, the activity of high pI α-amylase is much higher than low pI α-amylase. The aim of this study was to determine the relationship between preharvest sprouting resistance of hulless barley and the gene structure of Amy1 gene which encodes high pI α-amylase. The coding region and cDNA of Amy1 gene of three resistant accessions, including ‘XQ32-5’ (TR1), ‘XQ37-5’ (TR2), ‘XQ45-7’ (TR3), two susceptible accessions ‘97-15’ (TS1), ‘9657’ (TS2) and one highly dormant barley accession ‘SAMSON’ (SAM) was cloned. Analysis of their DNA sequences revealed there were three exons and two introns in Amy1 gene. Thirteen variable sites were in exon2 and exon3, 2 variable sites were in intron2. SAM and TS1 had a GAACT insert segment in the same site in intron2. Only 8 variable sites caused the change of amino acid residues. There were 99% of similarity between the tested hulless barley and some of the variable sites might be related with preharvest sprouting resistance. Then, we investigated the expression level of Amy1 gene in the 7-day germination test. Results of quantitative real-time PCR indicated that the relative expression trends of Amy1 gene were the same but had significant differences in the increase fashion between hulless barleys and no detectable expression was found in SAM. Susceptible accessions had earlier expression and faster increase and reached the maximum on day 4 ~ day 5. Besides, total transcripts level was found lower in resistant accessions than susceptible accessions. This study indicated that α-amylase activity was highly related to the transcription level of Amy1 gene which not correlated to missense mutation sites. In conclusion, hulless barley, especially the landraces from Qinghai-Tibetan Plateau in China possesses high degree of variation in PHS performance, which indicates the potential of Tibetan hulless barley as a good source for breeding of resistance to PHS. This study provides several donor parents for breeding of resistance to PHS. Our results also demonstrate that agronomic traits may be used as assistants for PHS resistance selection in hulless barley. Besides, analysis of high pI α-amylase coding gene Amy1 revealed the relative high expression of was Amy1 one of the mainly reason of different PHS resistance level in hulless barley.

青藏高原青稞淀粉资源的评价

青稞（Hordeum vulgare L.var．nudum Hook．f.），即裸大麦，是兼食用、饲用和酿造于一体的作物，有着重要的利用价值。淀粉是青稞籽粒中含量最多、最重要的碳水化合物，淀粉含量、直支淀粉比将会直接影响淀粉的功能特性，进而影响淀粉的应用领域。我国青藏高原青稞的栽培和食用历史悠久，特色青稞资源极其丰富。目前关于青藏高原青稞淀粉特性的报道还不多见，筛选和培育特色淀粉青稞利于拓展青稞的应用领域, 从而提高其经济价值。 本研究以114份青藏高原青稞品种(系)为实验材料，通过SDS-PAGE对材料的胚乳淀粉颗粒结合蛋白（SGAPs）进行分离，确定各蛋白的分子量大小、组合类型和多态性等。然后按照国标法测试材料的籽粒总淀粉含量和直链淀粉含量，通过微型糊化粘度仪分析相应的淀粉糊化特性，最后使用显微镜观察比较了青稞的淀粉颗粒形态特征。主要结果如下： 1、114种青稞中共分离出20种不同的SGAP条带，条带分子量为35.00～112.39 KDa，分布频率为12.28～97.37%。材料含有的SGAPs条带数从10到14不等，超过一半的材料含11种SGAP条带。20种条带形成16种组合类型，其中西藏地区青稞包含所有16个组合类型，四川地区青稞包含其中12个组合类型。青藏高原青稞籽粒淀粉颗粒结合蛋白的差异很大，遗传多样性丰富。 2、114份青稞的总淀粉含量、直链淀粉含量、直支淀粉比、峰值粘度、糊化温度和峰值温度的变幅分别为51.26～66.70%、14.64～29.74%、0.17～0.42、194～1135BU、58.8～65.2℃和81.4～92.4℃，相应的平均值分别为59.82%、23.60%、0.31、722.30BU、62.1℃和88.8℃。群体在总淀粉含量、直链淀粉含量、直支淀粉比、峰值粘度、糊化温度和峰值温度上的分布具有明显的正态性；所有胚乳淀粉体的淀粉粒都呈复粒结构。对西藏和四川的材料进行了分组比较， 两地区的青稞在直链淀粉含量和直支淀粉比上的差异达到显著水平。 3、筛选出18份具有特殊淀粉特性的青稞品种，其中5份材料的总淀粉含量超过65%，包括NB63-1、NB67、甘孜白六棱、98221-1和NB63；3份材料的直链淀粉含量大于29%，包括藏青85、藏青3号和喜马拉6号；8份材料的直支淀粉比小于0.25，包括99033-6、春青稞、阿坝330、Jan-03、米麦114、396、NB63-1和92013；7份材料的糊化温度低于60℃，同时材料的峰值粘度大于1000BU,并且峰值温度低于90℃，包括足捉春、Jan-03、阿坝330、米麦114、春青稞、20003和阿青5号。 4、各淀粉特性间存在高度相关性。直链淀粉含量和直支淀粉比与糊化温度成极显著正相关，与峰值粘度成极显著负相关，与A型淀粉粒数量和大小呈负相关。不同SGAPs组合的品种之间，淀粉含量和淀粉糊化特性间差异均达显著水平。SGAP2、SGAP5、SGAP6和SGAP7可能对籽粒直链淀粉含量、直支淀粉比和糊化温度有正向效应；SGAP3、SGAP9∼SGAP20可能对峰值粘度有正向效应。 本研究对青藏高原青稞淀粉资源进行了较为全面的评价，对该区青稞淀粉特性有了系统的认识。研究筛选出的特殊青稞品种可作为青稞育种和青稞淀粉工业应用的潜在资源，淀粉特性差异巨大的众多青稞品种也为拓宽青稞应用领域提供了丰富的资源保障。本研究对部分SGAPs在性质上的鉴定和功能上的初步推断为青稞材料的筛选提供了指导，也为品质育种提供了理论参考。 Hulless barley (naked barley, Hordeum vulgare L.) is a short- season, early maturing crop with a wide range of adaptation. It has been attracting more and more attention due to its superior nutrition and extensive industrial applications. Starch is the main ingredient in hulless barley seeds which makes up 65 percent of hulless barley���s dry weight. The ratio of the amylose/amylopectin and the size, shape, distribution of starch granules can affect the physico-chemical and functional properties of starch, which may turn affect its utilizations. The Qinghai-Tibet Plateau, which is located in southwestern China, is a typical area of vertical agricultural ecosystem and one of the barley origin centers with abundant hulless barley resources. There are little reports about hulless barley in Qinghai-Tibet Plateau at present. To screen and cultivate some characteristic hulless barley can improve its value. An improved SDS-PAGE was used to identify SGAPs combination of 114 hulless barley varieties. Starch content (total starch and amylose starch) was determined according to the standard methods GB5006-85 and GB/T 15683 using PerkinElmer M341 Precision Automatic Polarimeter and UV spectrophotometer 755B respectively. The pasting properties were measured by BRABENDER Micrio Visco-Amylo- Graph 803201. The morphology of starch granules were observed and compared with Axioplan 2 Imaging light microscopy. The following were the results obtained: 1. There were 20 major SGAPs presented in 114 varieties, with the molecular weight ranged from 35.00 to 112.39 KDa, and the frequencies ranged from 12.28% to 97.37%. The number of SGAP bands in each accession varied from 10 to 14, more than half of the population had 11 bands. There were 16 distinct SGAP patterns in the 114 varieties, the Tibet hulless barley had all of the 16 types and the Sichuan hulless barley had 12 types. The results indicated the Qinghai-Tibet Plateau hulless barley had a polymorphism of the SGAPs. 2. The ranges of the total starch content, amylose content, Am/Ap, peak viscosity, pasting temperature and peak temperature of the 114 hulless barley were 51.26～66.70%,14.64～29.74%,0.17～0.42,194～1135BU,58.8～65.2 and 81.4℃～92.4, with an average of ℃59.82%, 23.60%, 0.31, 722.30BU, 62.1 and 88.8,℃℃ respectively. The distributions of the total starch content, amylose content, Am/Ap, peak viscosity, pasting temperature and peak temperature were visibly normal school. All of the amyloplasts in endosperm of varieties showed bimodal size distributions.The main starch properties of hulless barley from Tibet and Sichuan were separated and compared, the differences on amylose content and Am/Ap were obvious. 3. Eighteen accessions which had special starch properties were screened out. Five accessions with total starch content beyond 65%, including NB63-1, NB67, Ganzibailiuleng, 98221-1 and NB63; three accessions, Zangqing85, Zangqing3 and Ximala6, with the highest amylose content (>29%); five accessions with Am/Ap less than 0.25, including 99033-6, Chun Qingke, A Ba 330, Jan-03, Mi Mai114, 396, NB63-1 and 92013; seven accessions had a pasting temperature under 60, ℃meanwhile their peak viscosity beyond 1000BU and their peak temperature under 90℃，including Zu Cuochun, Jan-03, A Ba 330, Mi Mai 114, Chun Qingke, 20003 and A Qing 5. 4. There were high correlations between starch properties. Amylose content and Am/Ap were positively correlated to pasting temperature, negatively correlated to peak viscosity, negatively correlated to the number and granule size of A-type granule. Different SGAP combinations caused significant diversities in starch content and pasting properties. SGAP2, SGAP5, SGAP6 and SGAP7 may have positive effect on amylose content, Am/Ap and pasting temperature; SGAP3, SGAP9∼SGAP20 may have positive effect on peak viscosity. Our research made a comprehensive evaluation on the hulless barley starch from the Qinghai-Tibet Plateau, we can get a systemic understanding. Some special accessions were screened out can be used on the hulless barley breeding lines and industries utilization.The combination of the SGAPs may become a criterion to evaluate the hulless barley endosperm starch quality. Consequently, the results will be good information for further studies on the hulless barley.

Isolation and characterization of a novel variant of HMW glutenin subunit gene from the St genome of Pseudoroegneria stipifolia

The x- and y-type high molecular weight (HMW) glutenin subunits are conserved seed storage proteins in wheat and related species. Here we describe investigations on the HMW glutenin subunits from several Pseudoroegneria accessions. The electrophoretic mobilities of the HMW glutenin subunits from Pd. stipifolia, Pd tauri and Pd strigosa were much faster than those of orthologous wheat subunits, indicating that their protein size may be smaller than that of wheat subunits. The coding sequence of the Glu-1St1 subunit (encoded by the Pseudoroegneria stipifolia accession PI325181) was isolated, and found to represent the native open reading frame (ORF) by in vitro expression. The deduced amino acid sequence of Glu-1St1 matched with that determined from the native subunit by mass spectrometric analysis. The domain organization in Glu-1St1 showed high similarity with that of typical HMW glutenin subunits. However, Glu-1St1 exhibited several distinct characteristics. First, the length of its repetitive domain was substantially smaller than that of conventional subunits, which explains its much faster electrophoretic mobility in SDS-PAGE. Second, although the N-terminal domain of Glu-1St1 resembled that of y-type subunit, its C-terminal domain was more similar to that of x-type subunit. Third, the N- and C-terminat domains of Glu-1St1 shared conserved features with those of barley D-hordein, but the repeat motifs and the organization of its repetitive domain were more similar to those of HMW glutenin subunits than to D-hordein. We conclude that Glu-1St1 is a novel variant of HMW glutenin subunits. The analysis of Glu-1St1 may provide new insight into the evolution of HMW glutenin subunits in Triticeae species. (C) 2007 Elsevier Ltd. All rights reserved.

The effect of supplementary feeds on the bodyweight of yaks in cold season

The present study was conducted to determine the effects of supplementary feeds, oat hay (OH), highland barley straw (HBS) and multi-nutrient blocks supplementation (UMMB) on reducing liveweight losses of both yak cows and calves grazed on low quality pastures during cold season. The trials of OH and HBS supplementation were conducted by using completely random design on 104 yak cows between 6 and 12 years of age as the following treatments: pure grazing (41 animals, body weight 230 67 kg) as control (CK); grazing+1.5 kg DM of OH per head daily (30 animals, body weight 216 28 kg); gazing. 1.5 kg DM of HBS per head daily (33 animals, body weight 221 34 kg). The trial of UMMB was conducted on three types of yaks, 1-year calves (8-12 months old, body weight 61.1 6.9 kg), 2-year calves (18-24 months old, 98.0 11.3 kg) and yak cows (164.5 27.1 (S.D.) kg) with 20 animals in control group (CK) and 20 animals in supplement group for each type by using completely random design as the following treatments: pure grazing for CK group; grazing+ 150, 250 and 500 g UMMB per day averagely for 1-year calf, 2-year calf and cow at night. The results indicate that the animals supplemented with oat hay received body weight gain (32 20.7 g day(-1)), while those supplemented with highland barley straw still suffered from body weight loss (-56.7 39.3 a day(-1)); UMMB supplementation can decrease the body weight loss by 109.7%, 86.6% and 63.4% for the 1-year calves, 2-year calves and yak cows, respectively, as compared with pure grazing. Around US$1.60 output can be achieved on the basis of US$1 input for UMMB supplementation in the farming systems of the 1-year calves, 2-year calves and yak cows, while US$1 input can produce US$1.55 and 1.14 output for OH and FIBS supplementations, respectively, in yak cows' farming system. It can be preliminary concluded that UMMB supplementation was the most economic way to alleviate body weight loss of grazing yaks over cold season, and the higher productive returns were obtained from OH supplementation for grazing yak cows during winter/spring months. © 2004 Elsevier B.V All rights reserved.