2020-04ml 目录

目的幼苗幼树空间分布是影响森林群落动态变化的重要因素,对云冷杉不同苗高进行空间分布格局分析,并分别研究幼苗幼树与周围小径木和大径木的关联性,以探究云冷杉幼苗空间变化规律和空间结构对幼苗幼树的影响,提出促进云冷杉天然更新措施,为云冷杉林可持续经营提供理论依据。方法本研究于2017年8月在金沟岭林场设置长白山云冷杉针阔混交林标准地（0.36 hm2）,应用SADIE空间分析方法,分析了林分不同苗高级云冷杉的分布规律。结果云冷杉林分总体直径分布呈反“J”型;色木槭与冷杉直径分布呈增长型,云杉和红松为多峰波动曲线。林分内小径木聚集指数I_a > 1,且随机化检验概率P_a < 0.025,总体呈聚集分布;大径木P_a值为0.025 ~ 0.975,呈随机分布;随林木生长发育,林分空间分布由聚集逐渐趋于随机。冷杉作为主要树种之一,其小径木与大径木空间分布与林分总体空间分布相似;云杉呈随机分布;冷杉幼苗与冷杉大径木具有明显的空间关联性,并与云杉呈正相关,其幼苗具有较好的耐荫性;云杉幼苗总体呈空间分离或空间不相关,仅与冷杉大径木呈空间正相关性,云杉与冷杉为良好的伴生树种。苗高Ⅰ、Ⅱ、Ⅲ级冷杉幼苗空间关联性相似,苗高Ⅳ、Ⅴ级的幼树空间关联性相似;苗高Ⅰ、Ⅱ级云杉幼苗空间关联性相似。结论林分内,冷杉幼苗株数较多,与大多林木呈空间正相关,有利于冷杉的林下更新。云杉木材质量较好,作为主要目的树种,其大径木及幼苗株数较少,且幼苗幼树与小径木及大径木均呈空间分离,林下云杉更新状况较差,因此,在云杉幼苗时期采取适当的遮蔽措施,幼树时期进行采伐等经营措施,促进云杉幼苗生长,并通过对冷杉株数的合理控制,调整林分空间结构,为其释放生长空间,以增加云杉母树株数,提高产种能力及产量,促进云杉更新。      

2020-02ml 目录

目的采用非线性分位数回归方法构建樟子松树干削度方程,并对比分析9个分位数（τ = 0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9）模型和传统的非线性回归削度方程的预测精度。方法以七台河市林业局金沙林场154株人工樟子松干形数据为研究对象,选取简单削度方程、分段削度方程和可变指数削度方程,利用非线性回归和非线性分位数回归方法构建樟子松树干削度方程。采用确定系数（R2）、平均误差（MAB）、相对误差（MPB）、均方根误差（RMSE）为统计指标对构建的削度方程进行对比分析。结果（1）在9个分位点（τ = 0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9）下的各削度方程都可以收敛,分位数回归方法可以灵活预测各分位点树干曲线的变化。（2）与非线性回归相比,基于中位数（τ = 0.5）时的各削度方程在拟合过程中表现最好,其中以可变指数削度方程表现最优。（3）检验结果也表明:相对于非线性回归的各削度方程,基于中位数（τ = 0.5）的简单削度模型的MAB和MPB均下降26.7%,RMSE下降19.9%;基于中位数（τ = 0.5）的分段削度方程和可变指数方程预测能力较强。（4）中位数回归的各削度方程在树干大部分的预测能力都优于相应的非线性削度方程。结论分位数回归方法是一种稳健的建模方式,基于中位数（τ = 0.5）的可变指数削度方程的预测精度最高,适合该区域樟子松树干干形的预测。      

2020-01ml 目录

目的传统的森林生物多样性保育价值缺乏统一的评估标准,不便于结果比较。目前,物种保育价值模型主要通过林分丰富度核算森林当前保育价值,但未考虑土地利用类型对树种生长的潜在影响,不能体现森林的长期保育价值。近年来,逐步完善的生物多样性价值评价新理论中,尚未出现适用于森林生物多样性保育价值评估的方法。方法通过分析福建土地利用情况计算生境质量并划分等级,首次提出生境质量调整系数。利用生境质量和森林分布之间的地理对应关系,构建出基于生境质量的森林生物多样性保育价值模型（HQ-BPV模型）,由此将土地利用的影响融入福建省森林生物多样性保育价值评估中。结果（1）福建省生境质量分布整体呈现内北高南低、内陆高沿海低和自然保护区高非自然保护区低的态势。（2）福建主要植被类型的平均生境质量按大小排序为针阔混交林 > 竹林 > 阔叶混交林 > 阔叶林 > 灌丛 > 针叶林 > 经济林。（3）调整后栎类和栲类的物种单位面积保育价值最高,为3.58万元/hm2,调整后各优势树种（组）的物种保育总值中阔叶混交林价值最高,为658.42亿元,构建模型核算的福建省森林生物多样性保育价值为2 056.39亿元。结论基于能够量化土地利用影响的生境质量指数,构建出HQ-BPV模型用于评价福建省森林生物多样性保育价值,为森林生物多样性的保育价值核算提供新思路,有利于森林资源监测和结构优化。      

2020-05ml 目录

    目的   理解林木蒸腾日内变化及其对主要环境因子的响应规律,进一步阐释短时间尺度下环境条件调控森林蒸腾的机理。    方法   以宁夏六盘山香水河小流域华北落叶松人工林为研究对象,在2018年生长季（5—10月）连续监测样树的树干液流变化,并同步观测气象条件和土壤湿度,分析小时尺度的林分蒸腾对环境因子的响应,并建立多因素影响的蒸腾模型。    结果   （1）在小时尺度上,林分蒸腾量（T）对太阳辐射（R_s）和饱和水汽压差（VPD）的响应均呈二次多项式函数关系;随R_s和VPD的增加,T均先增加,当R_s和VPD分别达到666.7 W/m2和1.86 kPa后达到峰值,然后逐渐减小。（2）T对土壤可利用水分（REW）的响应符合趋于饱和的指数关系,T随REW的增加表现为先增加,当REW > 0.3后,T逐渐趋于稳定。（3）在确定T响应R_s、VPD和REW的类型并耦合形成蒸腾模型后,利用生长季内的奇数天小时观测值进行拟合参数,并用偶数天小时观测值进行验证,得到T响应多因素变化的耦合模型:T = (? 6.347 0 × 10? 5\begin{document}${R_{\rm{s}}^2}$\end{document}? 0.637 0R_s ? 208.734 8) × (? 0.003 2VPD2 + 0.013 8VPD + 0.001 7) × (? 0.008 1 ? 0.004 6(1 ? exp(? 12.469 6REW))),该模型在校准阶段（R2 = 0.74,纳什效率系数（NSE）= 0.82）和验证阶段（R2 = 0.77,NSE = 0.84）均表现出较好的模拟效果。    结论   在小时尺度上,林分蒸腾量可以由耦合了太阳辐射、饱和水汽压差和土壤可利用水分影响的耦合模型进行较好的预测。本研究结果可为精确预测变化环境下的华北落叶松林分日内蒸腾提供理论基础,同时模型的构建方法可为其他区域和其他树种的林分蒸腾模型的建立提供参考依据。      

2020-01ml 目录

为了评价抗虫转基因甘蔗优良株系Bt2、Bt17对土壤生态环境造成的生态安全风险,采集抗虫转基因甘蔗Bt2、Bt17号株系及其受体非转基因品种ROC22甘蔗根际附近土壤,研究抗虫转基因甘蔗对其根际土壤蔗糖酶、酸性磷酸酶、中性磷酸酶及碱性磷酸酶的影响。结果表明,抗虫转基因甘蔗对其根际土壤酶活性的影响会因甘蔗生长周期、抗虫转基因甘蔗株系以及酶的种类而大有不同。Bt2的土壤酶活性在甘蔗生长的各个时期均未与对照存在差异,而Bt17对土壤酶活性的影响则较为复杂。与受体非转基因甘蔗品种ROC22相比,Bt17根际的土壤蔗糖酶活性在甘蔗生长的任何时期都无差异;而土壤酸性磷酸酶活性在甘蔗生长的各个时期均显著高于对照;中性磷酸酶活性则在苗期、分蘖期、成熟期显著高出对照,而在生长期与对照无差异;土壤碱性磷酸酶活性在苗期、生长期时与对照无差异,但分蘖期和成熟期显著低于对照,说明Bt2号株系对土壤酶活性并未产生影响,Bt17号株系对土壤酶活性可能产生较小的影响。     

第二代测序技术在植物遗传研究中的应用

第二代测序技术,又称"新一代测序技术",其较Sanger测序具有通量高、测序成本低和测序时间短等特点。介绍了以454、Illumina/Solexa和SOLi D为代表的第二代高通量测序技术的基本特性及其在植物全基因组、转录组和miRNA测序等领域的研究进展,并对新一代高通量测序技术的发展趋势进行了展望。     

基于高通量测序分析肺炎支原体感染姜曲海猪肺组织miRNA表达谱

为探明姜曲海猪感染肺炎支原体后肺组织miRNA(microRNA)表达谱和分子机制,选取50日龄姜曲海猪为实验猪,随机分成感染组和对照组,人工感染肺炎支原体28 d后,采集肺部组织,进行高通量miRNA测序,采用生物信息学软件进行miRNA鉴定和靶基因预测。结果显示:感染组和对照组的肺组织分别筛选到14 265 786条和14 000 588条小RNA纯净序列(clean reads)。与对照组相比,感染组中筛选到73个显著差异表达的miRNAs,其中39个表达上调,34个表达下调,从中随机选取4个miRNAs进行定量PCR(quantitative real time PCR,qRT-PCR)验证,感染组和对照组的表达水平与测序结果基本一致。73个差异表达miRNAs预测到1 685个靶基因和4 220个靶位点,靶基因预测筛选到8个与免疫调控相关的miRNAs。靶基因GO(gene ontology)分析显示,miRNA广泛参与生物过程、细胞组成和分子功能的调控。靶基因KEGG(Kyoto encyclopedia of genes and genomes)分析显示,miRNA参与调控细胞凋亡、ECM受体相互作用、粘附斑通路等信号通路。本研究通过高通量测序获得姜曲海猪感染肺炎支原体后肺组织miRNA表达谱,通过生物信息学分析筛选到与免疫调控相关的miRNAs和信号通路,为进一步阐明姜曲海猪的肺炎支原体感染机制奠定了基础。     

干旱胁迫下花生转录组与miRNA测序及相关基因的表达

[目的]探明干旱胁迫下花生转录组与miRNA测序及相关基因的表达,为深入分析花生抗旱分子机制提供理论依据.[方法]采用转录组测序和miRNA测序技术,对桂花37花生在15％ PEG6000模拟干旱胁迫下基因和miRNA的表达模式进行分析.[结果]转录组测序发现有13个显著持续上调的差异基因,其中7个为HSP分子伴侣家族...     

2019-3ml目录

目的地表温度是区域与全球尺度地表过程分析与模拟的重要参数,在地表与大气能量交换的过程中扮演着重要的角色。本文使用3种算法对北京地区Landsat-8影像进行地表温度反演,并使用MODIS 地表温度产品对反演结果进行交叉验证,评估Landsat-8用于地表温度反演的精度与适用性,为后续使用Landsat-8反演地表温度的研究提供参考。方法对反演地表温度所需的3个重要参数（大气平均作用温度、地表比辐射率、水汽含量）进行获取,得到3个算法的反演结果后,对各个算法的结果进行敏感性和差值分析,并将结果与同期MOD11A1地表温度产品进行对比分析。主要分析手段包括北京市不同行政区之间算法结果与温度产品的平均温度结果比较、不同地类之间算法结果与温度产品的平均温度结果比较,以及选取尺度效应较低、温度随时间变化较小的密云水库中心区域比较两者的温差。结果3个算法的反演结果总体平均温差不超过1 K,其中Cristóbal等提出的改进后单通道算法与其他两个算法的温差最小,Wang等提出的改进后单窗算法与其他两个算法的温差最大。Landsat-8地表温度反演结果普遍高于MODIS温度产品。通过选取密云水库中心区域对本研究反演结果和温度产品进行对比可以得出,3种算法结果与MODIS温度产品的总平均差值为1.373 K。结论反演结果总体上具有较理想的反演精度。Jiménez-Mu?oz提出的劈窗算法具有最好的敏感性分析结果,且与MODIS温度产品的结果最接近。Landsat-8反演结果与MODIS温度产品在总体地温分布规律上保持一致,但Landsat-8因具有更高的分辨率而能更好地分辨小地块不同地类的地温差异,在精确反演地表温度领域拥有更大的优势。      

玉米花粉高温胁迫相关miRNA的筛选及其靶基因分析

以高耐热玉米品种郑单958、低耐热玉米品种先玉335为材料,以正常生长条件为对照,在花期进行高温胁迫,通过miRNA高通量测序筛选玉米花粉中的差异表达miRNA,然后预测其靶基因,并对靶基因的本体特征和代谢通路进行富集分析。结果表明,共筛选到818个miRNA前体序列。在郑单958高温胁迫花粉与对照花粉对比组（HT958 vs CK958）中共筛选到19个显著差异表达miRNA序列,其中15个miRNA序列上调表达,4个下调表达,3个miRNA序列达到极显著水平（P<0.01）。对这19个显著差异表达miRNA的靶基因进行预测,共获得了503个基因转录本,其富集较多的GO生物学过程条目分别为转录调控DNA-模板、微管生物学过程、磷酸化作用、RNA聚合酶Ⅱ正向调控转录过程、甲基化作用等,KEGG富集较显著的代谢通路分别是谷胱甘肽代谢、碳代谢、花生四烯酸代谢、糖酵解/糖异生、叶酸生物合成等。在先玉335高温胁迫花粉与对照花粉对比组（HT335 vs CK335）中共筛选到15个显著差异表达miRNA序列,其中7个miRNA序列上调表达,8个下调表达,1个miRNA序列达到了极显著水平...     

胡杨异形叶抗氧化能力的比较

目的胡杨是我国西部荒漠地区抗逆性较强的树种,具有异叶性。前人的研究发现,胡杨叶片从披针形叶到锯齿卵圆形叶其结构趋于旱生性,且抗性逐渐增强。本文拟从活性氧清除机制着手,深入探讨胡杨异形叶的抗氧化能力,以期阐释其抗逆性的生化机制。方法以成年胡杨披针形叶、卵圆形叶和锯齿卵圆形叶3种典型形态叶为试验材料,对其·O₂?产生速率、MDA含量、抗氧化酶（SOD、POD、CAT、APX、GR）活性和抗氧化物质（黄酮、类胡萝卜素、AsA、GSH）含量随季节（5、7和9月）变化的规律进行了比较研究。结果（1）·O₂?产生速率随季节而变化,·O₂?产生速率在锯齿卵圆形与卵圆形叶中递增,在披针形叶中先增后减;MDA含量在锯齿卵圆形叶中递减,在卵圆形叶中先减后增,而在披针形叶中持续递增。值得注意的是,9月份披针形叶中·O₂?产生速率最低而MDA含量最高。（2）锯齿卵圆形叶与卵圆形叶的SOD活性先降后升而CAT活性先升后降,POD活性持续升高而GR活性递减。APX活性在锯齿卵圆形叶中先升后降,在卵圆形叶中一直升高。而披针形叶的SOD、POD、CAT、APX和GR活性均呈先升后降的变化趋势。整体而言,锯齿卵圆形叶的抗氧化酶活性 > 卵圆形叶 > 披针形叶。其中,POD活性差异最为明显。（3）黄酮和类胡萝卜素在7月份大量积累,9月份AsA含量最高而GSH含量最低。黄酮和类胡萝卜素含量在锯齿卵圆形与卵圆形叶中先增后减,在披针形叶中递增;AsA含量在锯齿卵圆形与卵圆形叶中递增,在披针形叶中先减后增;3种形态叶的GSH含量均递减。（4）锯齿卵圆形叶中·O₂?产生速率与POD活性及AsA含量呈极显著正相关（P < 0.01）,卵圆形叶中·O₂?产生速率与POD和APX活性及类胡萝卜素和AsA含量存在极显著正相关（P < 0.01）,披针形叶中·O₂?产生速率则与POD、CAT和APX活性及黄酮和类胡萝卜素含量都存在极显著正相关（P < 0.01）。结论胡杨3种形态叶的抗氧化能力不同。锯齿卵圆形叶抗氧化能力最强,卵圆形叶次之,披针形叶最弱。3种形态叶随季节变化的抗氧化应答机制也存在差异。夏季锯齿卵圆形叶与卵圆形叶中以POD、CAT和APX作为主要清除ROS的抗氧化酶,秋季则以高活性的SOD和POD作为主要的抗氧化酶。同时,秋季卵圆形叶中APX活性也维持在较高水平。而披针形叶中以SOD、POD、CAT、APX和GR协同清除ROS。此外夏季胡杨叶片大量积累黄酮和类胡萝卜素消除ROS;秋季则通过AsA-GSH循环清除更多的ROS。      

1月龄五指山猪与长白猪骨骼肌miRNA转录组比较

为研究microRNA(miRNA)在五指山猪和长白猪肌肉生长发育中的差异和探究差异miRNA对骨骼肌发育转录后调控的影响,以1月龄的五指山猪和长白猪背最长肌为材料,通过转录组测序和生物信息学分析筛选与骨骼肌发育相关的差异miRNA。结果显示,从2个文库中共鉴定获得318种己知的猪miRNA,五指山猪和长白猪分别鉴定出312个和306个已知miRNA,同时分别发现了83个和88个新的miRNA。获得17个差异显著的miRNA,其中16个miRNA表达上调,1个miRNA表达下调。靶基因预测和生物学功能预测发现17个差异miRNA,共预测到535个靶基因,靶向作用于mTOR信号通路和肌动蛋白细胞骨架调控的ssc-miR-486,ssc-miR-204,ssc-miR-338和ssc-miR-885-3p 4个miRNA可能参与骨骼肌生长发育过程。     

2019-12ml 目录

林业遥感经历了航片判读、卫片目视解译、蓄积量定量估测阶段后,已经进入参数定量反演阶段。在林业对遥感业务化监测和精度提升的强烈需求下,定量遥感逐步与林业遥感交叉融合,林业遥感定量化研究的人才队伍、理论模型、数据源和应用方法逐渐成熟。本文提出了林业定量遥感的概念和框架,指出了其中关键的科学问题:（1）如何使遥感解译、建模和反演适应复杂的森林状况;（2）如何提高参数反演的准确度;（3）如何丰富林业遥感数据源;（3）如何发展更为智能化的遥感数据自动化算法。在介绍适合林业的定量遥感模型和通用反演方法的基础上,阐述了高光谱、热红外、激光雷达和微波遥感数据源的林业应用方法。未来林业定量遥感将在全波段数据统一建模和信息融合机制、机理模型反演、大数据融合等方面进行突破。      

玉米抗感粗缩病重要miRNA及其靶基因的诠释

借助高通量测序和生物学分析手段,筛选不同抗性品种间的差异表达miRNA,预测相关靶基因.对玉米抗、感粗缩病不同种质材料高通量测序鉴定出的差异miRNA的靶基因,在Genome、Gene Ontology等生物信息数据库中进行比对,同时做GO(Gene Ontology,简称GO)、京都基因和基因组百科全书(Kyoto Encyclopedia of Genes and Genomes,简称KEGG)富集性分析.从某些相关差异miRNA的靶基因及其作用位点的有关GO ID及其生物学功能注释中,发现一些可能与玉米粗缩病抗性有较大关系的miRNA,如zma-miR156a、zma-miR156k、zma-miR1432等;通过进一步的miRNA靶基因KEGG富集性分析,获得代谢通路ID,根据其生物功能注释,推测也有一定数量的miRNA与玉米种质抗(感)病相关.结合所涉差异miRNA在玉米抗、感种质接毒前后上、下调的表现,对这些重要miRNA靶基因可能涉及粗缩病的致病机制或玉米抗性相关的GO ID、KEGG通路ID所蕴含的功能进行了诠释.研究旨在探明抗病相关miRNA通过靶向作用其靶基因,在抗病过程中可能发挥的重要生物学功能和作用途径,研究结果可为探讨玉米抗粗缩病的抗病机制研究奠定理论基础.     

胡杨叶与根中特异性表达基因的表达谱分析

以2年生胡杨实生苗为研究材料,通过NaCl溶液处理24 h,分别收集NaCl处理后的细嫩白根及叶片,采用CTAB法抽提胡杨盐胁迫下的叶与根的RNA,用基因芯片技术研究胡杨盐胁迫后叶与根中基因的差异表达。经过比对,差异10倍以上的基因共有175个,其中在叶中表达上调的有120个,在根中表达上调的有55个。     

基于高通量测序的低温胁迫下赤桉miRNAs的挖掘与分析

目的 预测、挖掘和分析涉及赤桉Eucalyptus camaldulensis低温胁迫应答的miRNA,为研究其调控赤桉低温胁迫应答的分子网络奠定基础。方法 采用高通量测序对低温处理组和对照(CK)组的赤桉组培苗茎尖进行小RNA测序。以miRBase21.0、Rfam14.1和巨桉E. grandis基因组为参考数据库,利用Bowtie、miREAP和miRDeep2等软件进行miRNA预测,使用RNAfold对预测到的miRNA前体进行二级结构的折叠;采用psRNATarget预测靶基因,通过DEGSeq包分析差异表达的miRNA,并对它们进行GO注释和KEGG富集分析。结果 在赤桉中,共预测到隶属于54个家族的392个已知miRNA和97个新miRNA;其中,CK组共预测到282个已知miRNA,65个新miRNA;低温处理组共预测到329个已知miRNA,51个新miRNA。挖掘到80个在低温处理下显著差异表达的miRNA,包括55个上调和25个下调。GO基因功能注释和KEGG富集分析的结果表明,这些差异表达miRNA可能通过参与代谢通路、次生代谢物的生物合成、细胞膜的改变、信号转导和生物调节等响应低温胁迫。此外,还挖掘到25个可能与ICE1-CBFs-COR通路有关的miRNA。结论 借助高通量测序和生物信息学软件初步得到了低温胁迫下差异表达的赤桉miRNA,为进一步分析这些miRNA在赤桉低温胁迫中的分子功能提供一些参考。     

胡杨异形叶CYT和EST同工酶特征的初步研究

以胡杨不同发育阶段植株(五种树冠类型)的不同冠层异形叶为材料,调查同一单株上不同冠层同一节位叶片的叶形指数,测定其细胞色素氧化酶同工酶酶谱(CYT)和酯酶(EST)同工酶酶谱。探讨不同发育阶段胡杨植株叶形的变化和CYT、EST同工酶酶谱间的关系。结果表明:(1)不同发育阶段的胡杨植株都表现出叶形指数随冠层的升高而减小;(2)同一株不同冠层叶片的CYT同工酶酶谱带型基本一致,但各冠层叶片间酶谱带的表达量上存在细小差异。不同发育阶段植株CYT同工酶在质(酶带数目和迁移率)和量(酶的活性)的表达上各不相同,表现为随发育进程的推进,一些条带的表达或不表达;(3)EST同工酶的质(酶带数目和迁移率)和量(酶的活性)的表达在同一株不同冠层及不同发育阶段的胡杨叶片中均有差异。     

藏绵羊不同毛色皮肤组织miRNA表达谱及靶基因分析

【目的】藏绵羊又称藏系绵羊,是中国三大粗毛羊品种之一,主要分布在西藏及其毗邻的高寒牧区,如青海、甘肃、四川、云南等地。藏绵羊由于其被毛纤维粗长,毛质细腻且保暖性好,是制造藏式地毯的上好原料。而毛色作为一种重要的经济性状,同时制约羊毛的经济价值。然而,目前对于绵羊毛色的分子调控机制尚不明确,以藏绵羊为研究对象,对其不同颜色皮肤组织（黑色和白色）进行转录组测序,旨在探讨miRNA在藏绵羊不同毛色皮肤组织转录后水平中发挥的作用及其可能参与的调控通路。【方法】选取2只具有黑白花毛色的健康藏绵羊,减去体侧被毛使皮肤完全裸露并用75%酒精消毒处理,屠宰后快速切取皮肤组织并保存在组织样品保存液中防止RNA降解,提取RNA后通过miRNA测序和生物信息学分析,获得藏绵羊不同毛色（黑色和白色）皮肤组织miRNA表达谱,筛选组织差异表达miRNA并预测相关靶基因。通过靶基因分析,获得与调控藏绵羊毛色性状可能相关的信号通路。【结果】黑色和白色皮肤组织共获得85.76兆条原始读长,经过滤后得到85.08兆条clean reads;共鉴定出334个已知的miRNA和59个新的miRNA;从中获得23个差异表达miRNA,其中14个差异表达miRNA在白色皮肤组织中表达显著上调,9个表达下调。miR-2284b和miR-744仅在藏绵羊白色皮肤中表达,而miR-23b、miR-411a-5p、miR-30c、miR-423-3p和miR-324-5p则仅在藏绵羊黑色皮肤中表达,但表达量相对较低。表达量最高的miRNA为miR-10a,而倍数差异最大的为miR-23b。其中,miR-10a可参与调节多种信号通路,其作用涉及增殖、分化、凋亡、细胞粘附等各个方面,目前并无调控绵羊毛色方面的相关报道,但其功能及作用机制的研究仍在不断扩展。miR-23b则与Wnt、Notch等信号通路有关,这些信号通路中的部分基因也和黑色素的生成有密切关系。结合倍数差异和miRNA表达量分析,miR-411a-5p、miR103和miR-200b可能也和毛色调控密切相关。本研究共预测出981个靶基因可能参与毛色调控,多数基因和WNT信号通路、MAPK通路、EDNRB信号通路及黑素原生成通路（melanogenesis pathway）相关,其中黑素原生成通路显示和WNT信号通路、KIT信号通路及EDNRB信号通路相互关联。黑素原生成通路显示藏绵羊不同毛色可能最终通过MITF基因调控下游的TYR基因（酪氨酸酶基因）影响黑色素的合成。同时,筛选7个差异表达miRNA进行定量验证,结果表明,5个miRNA在黑色皮肤组织中上调表达,2个miRNA在白色组织中上调表达,定量结果与测序结果一致。【结论】获得了藏绵羊皮肤组织miRNA表达谱,并通过生物信息学分析得到可能和调控毛色性状相关的miRNA和信号通路。这些结果表明,藏绵羊毛色性状可能受到多种miRNA的调控并涉及多个信号通路,这有助于更进一步了解毛色转录后水平的调控过程,并对进一步的功能验证奠定基础。