渤海湾7种虾虎鱼科鱼类COI基因序列特征及分子分类研究.pdf

天津农业科学TianjinAgricultural Sciences畜牧兽医与水产养殖渤海湾7 种虾虎鱼科鱼类COI基因序列特征及分子分类研究2023,29(9):32-37房恩军，郭彪，郑德斌，张雪，王硕，王宇（天津市水产研究所天津市海洋牧场技术工程中心，天津30 0 2 2 1）摘要：为获得渤海湾天津海域常见虾虎鱼线粒体条形码序列及系统分类信息，采集渤海湾7 种虾虎鱼样本，通过形态初步鉴定物种，并提取DNA进行线粒体COI基因部分序列扩增测序，用系统进化分析及分子鉴定软件MEGA对渤海湾常见12 种虾虎鱼COI基因部分片断进行序列特征分析及分子分类分析。扩增得到7 种虾虎鱼COI基因片断，分析得到其平均GC含量低于AT含量,碱基在密码子位置分布中具有偏向性,且GC含量最高为第2 位碱基,与多种鱼类COI基因碱基特点相似。基因密码子碱基变异率为19.4%，转换颠换比为1.3。遗传距离分析及系统进化树分析理清了渤海湾天津海域12 种主要虾虎鱼的系统分类及进化关系。遗传距离显示，矛尾刺虾虎鱼(Amblychaeturichthys hexanema)和斑尾复虾虎鱼(Synechogobius ommaturus)为同种鱼。虾虎鱼COI编码基因密码子高变异率可能与其良好的环境适应性有关。关键词：渤海湾虾虎鱼；线粒体细胞色素c氧化酶亚基(COI)；系统进化中图分类号：Q917.4Sequence and Molecular Classification of COI Gene of Seven Species of Goby from the Bohai BayFANG Enjun,GUO Biao,ZHENG Debin,ZHANG Xue,WANG Shuo,WANG Yu(Tianjin Fisheries Research Institute,Center for Marine Ranching Engineering Science Research of Tianjin,Tianjin 300221,China)Abstract:In order to clarify the mitochondrial barcode sequence and taxonomic information of common goby fish in Tianjin area ofBohai Bay,seven species of goby fish were collected from Bohai Bay,identified by morphology,and DNA was extracted for amplifica-tion and sequencing of mitochondrial COI gene fragments.Phylogenetic analysis and molecular identification software MEGA wereused to analyze the sequence characteristics and molecular classification of COl gene fragments of 12 common goby species in BohaiBay.Seven COI gene fragments of Goby were amplified and sequenced,the average GC content was lower than AT content,the basehad a bias in the codon position distribution,and the highest GC content was the second base,which was similar to the COI gene baseof many fishes.The codon base variation rate was 19.4%,and the ratio of transitions to transversions was 1.3.Phylogenetic distanceanalysis and phylogenetic tree analysis clarified the phylotaxonomic and evolutionary relationships of 12 major goby species in Tianjinarea of Bohai Bay.Genetic distance indicated that Amblychaeturichthys hexanema and Synechogobius ommaturus were the samespecies.The high codon variation rate of codon encoding gene in Goby may be related to its good environmental adaptability.Key words:Goby of the Bohai Bay;cytochrome oxidase subunit I(CO)gene;phylogenetic analysis文献标识码：ADOI编码：10.39 6 9/j.issn.1006-6500.2023.09.006渤海湾是一个浅水淤泥质半封闭海湾，历史上饵料丰富,基础条件优越,是渤海鱼类的产卵场。近年来，由于围填海工程、水域污染和过度捕捞等因素的影响，渔业资源的生存环境恶化,影响了其生态系统结构和功能，渔业资源衰退，小型化、低值化问题突出，种类和数量也在减少。年度调查表明，天津海域产卵场逐渐退化，与历史调查数据相比，渤海湾鱼类浮游生物的丰度和多样性显著减少，经济鱼类渔获和产卵量也逐渐减少-2。而虾虎鱼由于环境适应力强，繁殖力强和生命周期短等特点3,其资源量稳定增加，已逐渐代替其他经济鱼类，成为渤海湾天津海域春季休渔期前鱼类群落及仔稚鱼的优势种,在海洋食物网中发挥着重要的作用。但是虾虎鱼种类繁多,且有关仔稚鱼的形态相近，相关资料较少，鉴定难度大,准确性难以保证,虾虎鱼仔稚鱼鉴定也是世界难点和是鱼类早期资源研究的重要内容。20世纪9 0 年代以来，国内外科学家发现线粒体一些短的基因（如细胞色素氧化酶第一亚基收稿日期：2 0 2 3-0 4-17基金项目：国家重点研发计划（2 0 19 YFE01223000；2 0 19 Y FD 0 9 0 2 10 2）；天津市海水养殖现代农业产业技术体系创新团队项目(ITTMRS2021005)作者简介：房恩军（19 6 8 一），男，天津人，高级工程师，主要从事渔业资源修复及评估技术研究。通讯作者简介：王宇（19 8 6 一），女，湖北鄂州人，工程师，硕士，主要从事渔业资源修复及评估技术研究。第9 期房恩军等：渤海湾7 种虾虎鱼科鱼类COI基因序列特征及分子分类研究33(COI)基因、12 s rRNA、细胞色素B(CtyB)等)进化速率快且相对保守，其中COI基因5端6 50 bp片断可有效区分种间物种,可作为物种分类的DNA条形码(DNABarcoding）,随后与之配套的数据库也日趋完善，广泛应用于生命科学、法医学、食品质量控制等领域，也为鱼类浮游生物种类的分子鉴定奠定了良好的基础4-9 。虾虎鱼亚目（Gobioidei)有两千余种虾虎鱼，腹鳍合成一个吸盘,是鲈形目中最大的一个类群。虾虎鱼多为小型底栖鱼类，主要以小型底栖甲壳动物为食，在渤海湾生态系统中能量流动和物质循环中有重要作用，是底层食物链的关键物种。根据渤海鱼类生物学记载，渤海湾天津海域历史上能采集到的虾虎鱼有12 种左右,其中拉氏狼牙虾虎鱼（Odontamblyopus rubicundus）、斑尾腹虾虎鱼（Sy n e c h o g o b i u s o mma t u r u s）具有一定的经济价值。除了小头栉孔虾虎鱼。常见的虾虎鱼的线粒体(COI)基因片断均能在数据库中获取。但是虾虎鱼形态多变、各地名称不统一等原因，需要通过扩增进一步明确序列信息。本研究选择渤海湾常见7 种虾虎鱼类作为研究对象,通过获取线粒体(COI)基因序列分析比较其序列变异情况，进一步进行系统发育及遗传分析,从而为鱼类DNA条形码库及鱼类浮游生物的鉴定积累基础数据，为渤海海洋食物网的研究和稳定，为渤海湾虾虎鱼资源保护和合理开发利用提供科学依据。1材料与方法1.1样品采集样品来源于2 0 18 年春季和夏季渤海湾水域地笼或拖网采样，参照渤海鱼类生物学共鉴定出矛尾刺虾虎鱼（Amblychaeturichthys hexanema)、裸颈镐虾虎鱼（TridentigerNudicervicus）、矛尾虾虎鱼(Chaeturichthysstigmatias)红狼牙虾虎鱼(Odontamblyopusrubicundus）弹涂鱼（Periophthalmus modestus）钟道虾虎鱼（Triaenopogonbarbatus）、小头栉孔虾虎鱼(Ctenotrypauchen microcephalus)7种7 个虾虎鱼。样本用9 5%的乙醇混合保存后取鱼鳍。1.2DNA 提取采用酚-氯仿法提取DNA,取鱼尾鳍0.5g放人1.5mL的离心管中,将鱼鳍剪碎后加入裂解液和蛋白酶涡旋混匀，离心后于56 消化3h，酚-氯仿抽提后2 倍无水乙醇沉淀，7 5%乙醇洗涤吹干后以ddH20溶解，4冰箱保存。1.3引物设计及PCR扩增在上海生工合成Ward2文献中鱼类COI通用引物:FishF1(5-TCAACCAACCACAAAGACATTGGCAC-3),FishR1(5-TAGACTTCTGGGTGGCCAAAGAATCA-3)。组建一个2 5LPCR反应体系：模板DNA2.0 ddH,06.5 lPremix Taq12.5 L引物各2.0 LPCR扩增的程序为：952 min95 变性 40 s92退火1 min72 延伸 1 min72 延伸5 minPCR扩增产物通过琼脂糖凝胶(1.5%)电泳1520min，凝胶电泳成像仪检测后送至上海欧易生物技术有限公司进行双向测序，测序引物为扩增引物。1.4数据分析从BOLD数据库中获取鲈形目的石首鱼科黄姑鱼（Ni b e a a l b i f l o r a）、马科的四指马（Eleutheronematetradactylum）鲽形目的舌鳎科短吻红舌鳎(Cynoglossusjoyneri）、半滑舌鳎（Cynoglossussemilaevis）共4种鱼类的同源序列作为外群，同时从NCBI和BOLD数据库下载天津可能存在的虾虎鱼序列，结合测序序列进行系统进化分析(表1）。对所获得的序列利用MEGA进行多序列比对，人工核查后，用“Statistics”程序计算序列碱基组成、遗传颠换值、变异位点、密码子偏向性等。不同个体间的遗传距离用 Distance程序下的Kimura双参数模型计算。用Modeltest检测出最适合的核苷酸序列分析模型后采用邻位连接法(Neighbor-Joining,NJ)构建进化树Bootstrap置信值重复抽样10 0 0 次。2结果与分析2.1渤海湾7 种虾虎鱼的 COI基因序列特征将扩增测序得到的COI序列进行blast比对,相似度大于9 8%的即为同一种，比对结果支持形态鉴定结果。多序列比对得到6 52 bp共有序列，证实COI可以用作物种鉴定条形码8。碱基特征分析得到7 种虾虎鱼的COI基因部35个循环34物种名称四指马黄姑鱼短吻红舌鳎半滑舌鳎矛尾刺虾虎鱼矛尾虾虎鱼钝尖尾虾虎鱼晴尾蝌蚪虾虎鱼红狼牙虾虎鱼弹涂鱼斑尾复虾虎鱼钟道虾虎鱼裸颈虾虎鱼暗虾虎鱼纹虾虎鱼小头栉孔虾虎鱼分序列符合硬骨鱼类G+C含量（44.6%）小于A+T含量(55.4%)的特点4。密码子的碱基使用频率存在一定的偏向性，平均GC含量最高为第2 位碱基（55.50%），4种碱基在第1、2、3位密码子分布含量序号1234567密码子第1位碱基第2 位碱基第3位碱基总计2.2虾虎鱼科12 种鱼类的种间遗传距离12种虾虎鱼的种间遗传距离平均值为0.2 1。矛尾刺虾虎鱼和斑尾复虾虎鱼的种间遗传距离最小（0.0 14),小于0.0 2 的临界值31,应为同种虾虎鱼，斑尾复虾虎鱼和纹虾虎鱼种间遗传距离最大(0.435)。天津农业科学表1进化树中16 种鱼类COI序列信息拉丁名缩写EleutheronematetradactylumETNibeaalbifloraNACynoglossus joyneriCJCynoglossus semilaevisCSAcanthogobius hastaAHChaeturichthys stigmatiasCSAmblychaeturichthys hexanemaAHLophiogobius ocellicaudaLOOdontamblyopus rubicundusORPeriophthalmusmodestusPMSynechogobius ommaturusSOTriaenopogonbarbatusTBTridentigerNudicervicusTNTridentigerobscurusTOTridentigertrigonocephalusTTCtenotrypauchenmicrocephalusCM都不平均。本研究中COI基因密码子变异情况表明，共有不变位点（52 2 个)多于变异位点（12 6 个），转换位点(7 1个)多于颠换位点（55个）,转换与颠换比为1.3,其中第1位碱基变异位点最多,第3位最少。表2 潮渤海湾7 种虾虎鱼COI基因部分序列碱基GC含量平均GC含量/%物种平均值矛尾刺虾虎鱼45.95矛尾虾虎鱼44.82弹涂鱼42.59红狼牙虾虎鱼44.68钟道虾虎鱼47.07裸颈虎鱼44.51小头栉孔虾虎鱼42.67平均44.61表3渤海湾7 种虾虎鱼COI基因序列的密码子碱基变异情况不变位点/个(i)颠换位点/个(sv)10956201112123522702.3虾虎鱼科鱼类的分子系统进化树系统树显示，7 种虾虎鱼中，红狼牙虾虎鱼与矛尾刺虾虎鱼距离较远，小头栉孔虾虎鱼和弹涂鱼聚为一支，矛尾刺虾虎鱼和斑尾复虾虎鱼聚为一支，纹虾虎鱼和暗虾虎鱼遗传距离较近。第2 9 卷分类地位序列来源鲈形目马鲮科EU595106鲈形目石首鱼科HM180717鲽形目舌鳎科JQ738456鲽形目舌鳎科JN242743刺虾虎鱼属本研究、HQ536241矛尾虾虎鱼属本研究、HQ711866钝尾虾虎鱼属JQ738533.JQ738527蝌蚪虾虎鱼属NC_020783狼牙虾虎鱼属本研究、JX679044弹涂鱼属本研究、KP638476复虾虎鱼属NC_020767编虾虎鱼属本研究、KF642960镐虾虎鱼属本研究JX679063虾虎鱼属HQ536532、H Q 536 533编虾虎鱼属KX254504、H M 18 0 9 2 8栉孔虾虎鱼属本研究第1位碱基第2 位碱基39.3556.7437.9654.8831.8054.4234.1058.3343.7854.4236.1155.3531.9754.3736.4455.50转换位点/个(si)513155第3位碱基41.6741.6741.6741.6743.0642.1341.6741.93R(si/sv)1.13.73.01.3第9 期房恩军等：渤海湾7 种虾虎鱼科鱼类COI基因序列特征及分子分类研究35表4虾虎鱼科12 种鱼类两两比对的种间遗传距离物种ORORCSAH0.247TB0.222OR0.177TT0.369TO0.236LO0.235PM0.205SO0.256TN0.240CM0.203382912233537895985860.0503讨论与结论3.17种渤海湾虾虎鱼类的COI基因序列特点COI基因部分序列第2 位碱基的GC含量最CS0.2370.2090.2460.2470.4030.2560.2600.2210.2260.2690.21187979178AH0.2310.2340.4350.2350.1910.2270.0140.2240.273100 Tridentiger.Nudicervicus-1裸颈编虾虎鱼-11Tridentiger.nudicervicus-2裸颈虾虎鱼-2Periophthalmus.modestus-2弹涂鱼-2-Tridentiger.barbatus-1钟道虾虎鱼-1Tridentiger.barbatus-2钟道虾虎鱼-2100Tridentiger.obscurus-1暗虾虎鱼-1LTridentiger.obscurus-2暗编虾虎鱼-2Tridentiger.trigonocephalus-1纹镐虾虎鱼-110olTridentiger.trigonocephalus-2纹编虾虎鱼-2Periophthalmus.modestus-1弹涂鱼-1,Odontamblyopus.rubicundus-1红狼牙虾虎鱼-110olOdontamblyopus.rubicundus-2红狼牙虾虎鱼-2H图1NJ法构建的虾虎鱼科鱼类分子系统树高,与石首鱼科COI基因碱基特征一致,与于亚男长江口虾虎鱼科鱼类COI基因碱基特征不一致等15。第3位碱基的GC含量比较恒定，而第1位碱基则因种类不同而呈现含量差异，第1位碱基能否作为不同种类区分的特征也需要进一步扩大样本量做深人TB0.2330.3970.1740.2330.2020.2350.2080.24071 Acanthogobius.hasta-2矛尾刺虾虎鱼-2100Synechogobius.ommaturus斑尾腹虾虎鱼LAcanthogobius.hasta-1矛尾刺虾虎鱼-1Lophiogobius.ocellicauda晴尾蝌蚪虾虎鱼-Eleutheronema tetradactylum 四指马鲮Cynoglossusjoyneri短吻红舌鳎-Cynoglossussemilaevis半滑舌鳎Nibeaalbiflora黄姑鱼100jAcanthogobius.hexanema-1钝尖尾虾虎鱼-11Acanthogobius.hexanema-2钝尖尾虾虎鱼-2Chaeturichthys.stigmatias-1矛尾虾虎鱼-1100lChaeturichthys.stigmatias-2矛尾虾虎鱼-2OR0.3990.2340.2340.2190.2360.2440.244TT0.4210.4000.3900.4350.3810.227Ctenotrypauchen.microcephalus小头栉孔虾虎鱼TO0.2260.2020.2230.2120.238LO0.2230.1950.2380.256PM0.2190.1940.015SO0.2270.260TN0.213CM36的研究。本研究中7 种虾虎鱼的密码子碱基使用存在一定偏向性,但与李明晖等16 研究结果不一致。由于其与基因长度、碱基组分、翻译起始效应、基因的表达水平等多种因素有关16-17 ,需要系统地对虾虎鱼亚目不同地理群体的该基因密码子碱基使用情况进行分析再下结论，或者精确到不同科属的具体情况来进行探讨。7 种虾虎鱼的COI基因颠换大于转换，转换/颠换值(R)值为1.3小于2,且第1位碱基承担着主要选择压力，表明虾虎鱼可能通过高度基因突变实现快速进化，从而增强其抗逆性和适应性:9,18 。3.2COI基因在渤海湾天津海域虾虎鱼科鱼类分类中的应用虾虎鱼的COI基因的种间和种内遗传距离相差2 0 倍以上，可作为天津海域虾虎鱼科鱼类分子分类的条形码和仔稚鱼分子鉴定的依据13。在系统进化树上，12 种虾虎鱼节点支持率较高，除了刺虾虎鱼属、晴尾蝌蚪虾虎鱼与其他虾虎鱼遗传距离较远,其他虾虎鱼均聚为一大类。于亚男等4研究成果也出现了有的虾虎鱼先与外群聚在一起,之后再与其他虾虎鱼相聚的现象。该基因条形码可以很好地区分不同种，但对于不同科或更大的分类界元的效果不太明显。由于有的虾虎鱼比较稀有、难以采集，本研究采集样品有限，需要加大样本量进一步优化进化树。形态学上，孔虾虎鱼和拉氏狼牙虾虎鱼具有相似的体型、体色、眼晴高度退化等特点，有些科学家将它们归为近盲虾虎鱼亚科，这与进化树分类结果比较一致。但是它们遗传距离大于0.2，确定是不同的种。本研究中,背眼虾虎鱼科的弹涂鱼与近盲虾虎鱼亚科的小头栉孔虾虎鱼聚为一支,相关研究也证实他们亲缘关系较近，具体还需加大样本进行确认。本研究得出，矛尾刺虾虎鱼和斑尾复虾虎鱼为同一物种，种间遗传距离仅为0.0 14,小于Hebert等13提出的临界值0.0 2 0；进化树上，节点支持率达到100%聚为一支。虾虎鱼是虾虎鱼科的一个属，渤海常见的有钟虾虎鱼（虾虎鱼）、纹虾虎鱼、裸颈虾虎鱼等。本研究中所有的虾虎鱼遗传距离较近，与已知分类结果一致。晴尾蝌蚪虾虎鱼由于头部形态特征与其他虾虎鱼相差较大,从形态上很好区分，与其他虾虎鱼科属遗传距离较远。天津农业科学由于虾虎鱼种类众多，变异较大，其鉴定除了分子鉴定还需结合形态学和生态地理分布等因素，也可以结合其他条形码基因综合分析其进化、地理分布及系统分类等。参考文献：1李宝华，张素青，王秀芹，等.2 0 0 4-2 0 0 9 年渤海湾（天津区域)鱼虾贝产卵场水质环境综合评价.水产科学，2010,29(11):661-665.2徐海龙，谷德贤，乔秀亭，等.渤海湾主要渔业资源长度与体质量关系分析J.南方水产科学，2 0 14，10(1)：57-63.3刘静，陈咏霞，马琳.黄渤海鱼类图志M.北京：科学出版社,2 0 15.4 SACCONE C,DE GIORGI C,GISSI C,et al.Evolutionarygenomics in Metazoa:the mitochondrial DNA as a modelsystemJ.Gene,1999,238(1):195-209.5柳淑芳，陈亮亮，戴芳群，等.基于线粒体CO1基因的DNA条形码在石首鱼科(Sciaenidae)鱼类系统分类中的应用J.海洋与湖沼，2 0 10,41(2)：2 2 3-2 32.6高天翔，毕潇潇，赵林林，等.基于线粒体Cytb基因全序列的松江鲈群体遗传结构分析.水生生物学报，2 0 13，37(2):199-207.7彭居俐，王绪祯，王丁，等.基于线粒体CO1基因序列的DNA条形码在鲤科属鱼类物种鉴定中的应用J.水生生物学报，2 0 0 9,33(2)：2 7 1-2 7 6.8彭士明，施兆鸿，侯俊利，等.银3个野生群体线粒体COI基因的序列差异分析.上海海洋大学学报,2 0 0 9,18(4):398-402.9毕潇潇，高天翔，肖永双，等.4种鳕鱼线粒体16 SrRNA、C O I 和Cytb基因片段序列的比较研究J.南方水产,2 0 0 9,5(3):46-52.10杨纪明.渤海鱼类的食性和营养级研究J.现代渔业信息，2 0 0 1,16(10)：10-19.11李明德，张洪杰，等.渤海鱼类生物学M.北京：中国科学技术出版社，19 9 1：6 2-6 8.12 WARD R D,ZEMLAK T S,INNES B H,et al.DNAbarcoding Australias fish species J.PhilosophicalTransactions of the Royal Society B:Biological Sciences,2005,360(1462):1847-1857.13 HEBERT P D N,STOECKLE M Y,ZEMLAK T S,et al.Identification of birds through DNA barcodes J.PLoSBiology,2004,2(10):e312.14 AQUINO L M G,TANGO J M,CANOY R J C,et al.DNA barcoding of fishes of Laguna de Bay,PhilippinesJ.Mitochondrial DNA,2011,22(4):143-153.第2 9 卷第9 期15于亚男，宋超，侯俊利，等.基于线粒体COI基因部分序列的长江口虾虎鱼科鱼类系统分类J.淡水渔业，2014,44(5):3-8.16李明晖，安长廷，李昂，等.黄渤海虾虎鱼DNA条形码分类体系构建J.中国水产科学，2 0 2 2，2 9(8):1179-1188.房恩军等：渤海湾七种虾虎鱼科鱼类COI基因序列特征及分子分类研究3717 FEDOROV A,SAXONOV S,GILBERT W.Regularitiesof context-dependent codon bias in eukaryotic genes J.Nucleic Acids Research,2002,30(5):1192-1197.18郭新红，刘少军，刘巧，等.鱼类线粒体DNA研究新进展J.遗传学报，2 0 0 4,31(9)：9 8 3-10 0 0.(上接第31页）参考文献：1张树庚，冯宜梅，陈丽萍.武威市粮改饲产业发展现状及推进对策J.甘肃科技，2 0 2 0,36(14)：1-3.2周波，韩小花，李小红，等.河南省夏播青贮玉米品种筛选与综合评价J.草业科学，2 0 2 1，38(2)：316-32 6.3侯湃，陈彩锦，张静妮，等.宁夏固原地区不同青贮玉米品种生产性能及营养品质研究J.草地学报，2 0 2 1，29(10):2346-2354.4高文辉，赵利妮，张玉飞.12 个青贮玉米品种在榆阳区种植对比试验J.养殖与饲料，2 0 2 1,2 0(12)：11-14.5席杰军，梁子栋，张钰靖，等.陕西关中地区31个青贮玉米品种比较试验J.草地学报，2 0 18,2 6(6)：136 3-136 7.6段震宇，张小伟，王婷，等.青贮玉米生物产量与主要农艺性状的相关性J.大麦与谷类科学，2 0 17，34(1):24-27.7雷志刚，王业建，梁晓玲，等.16 个青贮玉米品种农艺性状与产量的通径分析J.草食家畜，2 0 16(4)：43-49.8舒秋萍，包兴国，杨文玉，等.武威井灌区粮饲兼用玉米品种比较试验结果J.甘肃农业科技，2 0 0 2(5)：14-15.9周德录，李博文，董禄信，等.中部旱作区粮饲兼用玉米新品种比较试验初报J.农业科技与信息，2 0 18(18)：5-7,12.10李斌，张国龙.15个粮饲兼用玉米品种在武威市?州区沙漠化区域的品比试验J.农业科技与信息，2 0 19(15):19-21.11李博文，孙多鑫，李城德，等.河西灌区不同青贮玉米品种产量表现试验研究初报J.农业科技与信息，2021(6):37-39.12李圆，李城德，李世成，等.甘肃省玉米粮改饲发展情况与思考J.中国农技推广，2 0 2 1,37(10)：5-7.13刘卓，邵怀峰，温万，等.宁南地区34个青贮玉米品种农艺性状及营养品质评价研究J.饲料研究，2 0 2 1,44(11):98-104.14娄芬，李小冬，尚以顺，等.毕节地区适宜青贮玉米品种（系)筛选及营养价值评价J.草业学报，2 0 2 0，2 9(6):214-224.15张书兴，王云平，王筱，等.内蒙古半干旱区青贮玉米的适应性J.草业科学，2 0 2 2,39(6)：12 17-12 2 4.16中华人民共和国国家质量监督检验检疫总局，中国国家标准化管理委员会.青玉米品质分级：GB/T25882-2010S.北京：中国标准出版社，2 0 11.17娜日苏，项错锋，梁庆伟，等.8 个青贮玉米品种生产性能和饲用品质比较J.饲料研究，2 0 2 2，45(13)：105-110.18雷志刚，王业建，浩江，等.不同青贮玉米品种品质性状与农艺性状的相关分析J.新疆农业科学，2 0 17,54(4):694-699.19杨秀芳，陈玲玲，乌艳红，等.应用灰色关联度综合评价2 6 个青贮玉米的生产性能J.草业科学，2 0 12，2 9(1):105-111.