芝麻产量相关性状的多位点全基因组关联分析及候选基因预测

【目的】通过对芝麻产量相关性状的全基因组关联分析,挖掘与产量性状关联的SNP位点及预测候选基因,为通过分子标记辅助选择育种等方式提高芝麻产量提供技术基础。【方法】以363份不同遗传背景和地理来源的芝麻种质资源构成的自然群体为研究对象,调查2年2点4环境下8个产量相关性状（单株产量、单株蒴数、蒴粒数、千粒重、株高、主茎果轴长、始蒴高度和表观收获指数）的表型值,借助覆盖全基因组的42 781个SNP标记,利用多位点SNP随机效应混合线性模型（multi-locus random-SNP-effect mixed linear model,mrMLM）对8个产量相关性状进行全基因组关联分析,检测与产量相关性状显著关联的SNP位点,并预测候选基因。【结果】在4个不同环境下,8个产量相关性状表现出广泛的表型变异,变异系数为6.51%—33.57%;相关性分析表明单株产量与单株蒴数、株高、主茎果轴长、表观收获指数呈极显著正相关;方差分析表明产量相关性状的基因型效应、环境效应、基因型与环境互作效应均达到了极显著水平。通过多位点全基因组关联分析共检测到210个与产量相关性状显著关联的SNP,在2018年南阳环境下检测到47个SNP,解释表型变异的1.63%—17.29%;在2019年南阳环境下检测到35个SNP,解释表型变异的1.94%—11.90%;在2018年平舆环境下检测到35个SNP,解释表型变异的2.15%—15.90%;在2019年平舆环境下检测到53个SNP,解释表型变异的1.25%—11.13%;在4个环境的综合BLUP条件下检测到75个SNP,解释表型变异的1.44%—13.58%。上述210个SNP涉及到175个位点,其中10个位点在3个及以上环境中被重复检测到。在这10个位点基因组区域内,共鉴定到214个候选基因,其中156个候选基因具有功能注释,主要涉及植物代谢、生物调控、生长发育等生物学过程。根据功能注释筛选出4个可能与芝麻产量相关的候选基因,其中SIN_1006338编码1-氨基环丙烷-1-羧酸合酶3（1-aminocyclopropane-1-carboxylate synthase 3-like）,参与乙烯的生物合成;SIN_1024330编码碱性螺旋-环-螺旋（basic helix-loop-helix）转录因子,负向调控植物细胞和器官的伸长;SIN_1014512编码吲哚-3-乙酸-酰胺合成酶GH3.6（indole-3-acetic acid-amido synthetase GH3.6）,参与调控茎和下胚轴细胞的伸长生长;SIN_1011473编码泛素受体蛋白DA1（protein DA1-like）,参与调节植物细胞增殖周期。【结论】通过多位点SNP随机效应混合线性模型的全基因组关联分析,检测到175个与产量相关性状显著关联的位点,筛选出4个可能与产量相关的重要候选基因。     

甜瓜叶绿素含量全基因组关联分析及候选基因预测

【目的】挖掘与叶绿素含量显著关联的单核苷酸多态性(single nucleotide polymorphism, SNP)位点和候选基因，为甜瓜叶绿素含量改良提供分子靶点和基因资源。【方法】以118份具有广泛变异的甜瓜种质为自然群体，采用2 531 449个高质量SNP标记对叶片叶绿素含量进行全基因组关联分析，挖掘优异等位变异，并预测候选基因。【结果】甜瓜自然群体叶绿素含量趋向正态分布，包含5个明显的亚群。利用Q模型对2次试验的叶绿素含量及其最佳线性无偏预测(best linear unbiased prediction, BLUP)值进行关联分析，共检测到15个显著位点，分布在甜瓜第1、2、4、8、11、12号染色体上，表型贡献解释率为5.62%～6.69%。其中，8个位点的不同基因型之间存在显著表型差异。结合关联位点的候选区域和转录组数据，共鉴定到28个差异表达基因，其中MELO3C018513.2和MELO3C003666.2是改良甜瓜叶绿素含量的潜在候选基因。【结论】通过高质量SNP标记Q模型的全基因组关联分析，共检测到15个与叶绿素含量显著关联的位点，筛选出2个可能与叶绿素含...     

棉花产量构成因素性状的全基因组关联分析

【目的】 对铃重、衣分、单株铃数和籽指等棉花产量构成因素性状进行全基因组关联分析（genome-wide association study,GWAS）,发掘与其关联的标记位点、优异等位变异及候选基因,为棉花产量的分子育种提供理论依据。【方法】 以408份陆地棉品种（系）资源为材料,利用Cotton SNP 80K芯片,对6个环境的铃重、衣分、单株铃数和籽指4个产量构成因素性状进行基于混合线性模型（mixed linear model,MLM）的全基因组关联分析,检测与产量构成因素性状显著关联的位点、优异等位变异;进一步依据转录组数据的基因表达量,在显著关联的位点侧翼序列1 Mb区间挖掘可能的候选基因。【结果】 4个产量构成因素性状在不同环境下均表现出广泛的表型变异,其中,单株铃数变异系数最大为16.67%—22.66%,各性状的遗传率为48.4%—92.2%;除铃重与衣分间相关性不显著外,其他性状间均呈显著或极显著相关性;基于6个环境各性状表型数据的最佳线性无偏预测值（best linear unbiased prediction,BLUP）,GWAS共检测到分布于基因组的7个区间内23个与目标性状关联的SNP位点,其中,与铃重关联的位点5个,与衣分关联的位点1个,与单株铃数关联的位点9个,与籽指关联的位点8个,有3个位点（TM21094、TM21102和TM57382）同时与多个目标性状关联;鉴定到7个最优SNP位点的优异等位变异,分别为TM21099（TT）、TM57382（GG）、TM78920（CC）、TM53448（TT）、TM59015（AA）、TM43412（GG）和TM69770（AA）;利用转录组数据分析,在基因组的7个区间筛选到158个与产量形成可能的候选基因,GO富集分析和KEGG代谢途径分析发现,候选基因功能类别多样并参与了多种代谢途径。【结论】 在陆地棉品种（系）群体中共鉴定到23个与产量构成因素性状关联的SNP位点,筛选到158个可能与产量性状相关的候选基因。     

50份马铃薯种质资源表型鉴定和遗传多样性分析

为丰富福建地区马铃薯种质资源和遗传基因改良奠定材料基础,参照《马铃薯种质资源描述规范和数据标准》,对50份种质资源的25个表型性状进行多样性、相关性和聚类分析。结果表明,50份马铃薯种质资源质量性状多样性指数差异很大,变幅为0.52～2.01,其中肉色、开花繁茂性、薯形、叶色的遗传多样性指数分别达到2.01、1.81、1.61、1.52、0.52;数量性状遗传多样性指数除了主茎分枝数、蛋白质含量外,其他9个性状均达到2.10以上,单薯重、单株产量、单株结薯数、主茎分枝数变异系数较大,均达到27%以上。相关性分析方面,生育期与单株结薯数、株高、干物质含量呈显著正相关;商品薯率与单株结薯数呈极显著负相关,与单薯重呈极显著正相关,而单薯重与单株结薯数呈极显著负相关;单株产量与商品率、单薯重呈极显著正相关;株高与单株结薯数、茎粗呈极显著正相关,与商品率呈显著负相关,与单薯重呈极显著负相关;干物质含量与淀粉含量呈极显著正相关。通过主成分分析得到6个主成分,累计贡献率达90.691 0%,能较好反映该11个数量性状的主要信息。聚类分析将50份马铃薯种质资源划分为4个类群,第Ⅰ类群为中早熟低产型材料...     

甘蓝型油菜角果长度全基因组关联分析

【目的】挖掘与油菜角果长度性状显著相关的SNP位点及候选基因,为揭示油菜角果长度性状的遗传基础和分子机制提供理论依据,为油菜产量分子标记辅助选择育种奠定基础。【方法】在江西农业大学试验地和江西省红壤研究所试验地2个环境下考察300份甘蓝型油菜自交系的角果长度性状,利用简化基因组测序技术(specific locus amplified fragment sequencing,SLAF-seq)对300份甘蓝型油菜自交系基因组DNA进行测序并分析,利用获得的均匀分布于甘蓝型油菜基因组上的201 817个群体SNP(single nucleotide polymorphism,SNP)对角果长度性状进行全基因组关联分析(genome-wide association study,GWAS),探测与油菜角果长度显著相关的SNP位点,并基于群体连锁不平衡分析结果搜寻显著SNP位点两侧100 kb范围内的基因,通过BLAST获得关联区域内基因的注释信息,根据注释信息找出与性状相关的候选基因。【结果】农大试验地角果长度表型变异幅度为46.35—107.07 mm;红壤所试验地角果长度表型变异幅度为39.41—101.35 mm,两性状在2个环境下均表现出广泛表型变异。通过一般线性模型(general linear model,GLM)关联分析,农大环境下共检测到121个角果长度显著关联的SNP位点,分布在A04、A06、A08、A09、C02、C03、C06和C09等8条染色体上,其中,A09染色体上分布最多(83个SNP),红壤所环境下检测到22个角果长度显著关联的SNP位点,其中,1个在C09染色体上,其余21个均分布于A09染色体,在两地探测到20个一致性SNP位点;通过混合线性模型(mixed linear model,MLM)分析,农大环境下共检测到5个角果长度显著关联的SNP位点,其中,3个SNP位点与红壤所环境下检测到3个SNP位点一致,所有位点均位于A09染色体上。对MLM关联分析得到的显著SNP位点两侧100 kb区域内基因进行搜寻并进行功能注释,发现多个候选基因参与调节碳水化合物的运输与合成、花器官和种子的发育、信号转导等,它们可能通过上述功能影响油菜角果的生长,导致角果长度的差异。【结论】通过GLM和MLM两种分析方法探测到多个与油菜角果长度性状显著关联的基因位点,并在显著性位点附近搜寻到相关候选基因。     

催芽天数、GA和BR互作对秋作紫马铃薯萌芽、出苗和产量的影响

【目的】研究催芽天数(D)、赤霉素(GA)、油菜素内酯(BR)互作对秋作紫马铃薯种薯萌芽、田间出苗和产量的影响。【方法】以‘黑美人’品种的原种为试验材料,采用三因素裂裂区试验设计,对秋播马铃薯种薯进行催芽处理。【结果】催芽温天数、赤霉素和油菜素内酯及三因素的交互效应对马铃薯种薯的发芽率、单薯发芽数、芽粗、芽长、田间出苗率、单株主茎数、公顷主茎数和产量均存在显著的影响。播种前25℃催芽天数越长,单个种薯的发芽数和单株主茎数越大,产量越高。使用赤霉素催芽的种薯发芽数、单株主茎数和产量显著高于不使用赤霉素催芽处理。BR_(100)水平的种薯发芽数、单株主茎数和产量显著高于BR_(500)水平。不同催芽处理组合中,D2+GA_(60)+BR_(100)处理的单薯发芽数最高,单株主茎数和鲜薯产量处于较高水平。【结论】对产量贡献率较高的因子依次为出苗率、单株主茎数、株高和单株结薯数。产量与单株主茎数、公顷主茎数、单株结薯数都呈一元二次函数关系,理论最适单株主茎数为2.35个。     

基于SLAF-seq的华山松球果数量性状全基因组关联分析

【目的】挖掘与华山松球果数量相关的基因。【方法】以华山松球果数量性状表型数据筛选华山松高、低球果数量单株，利用前期已开发的单核苷酸多态性(single-nucleotide polymorphism,SNP)标记结合表型数据进行全基因组关联分析(genome-wide association study,GWAS)，构建华山松样品的系统发育树，发掘与华山松球果数量性状相关的SNP标记与候选基因，并进一步分析筛选出的基因在华山松针叶、木质部、韧皮部、树皮、幼嫩球果和树根的qPCR基因组织特异性表达。【结果】系统发育树显示：同一种源的单株基本处于相近的位置，而来自巍山的单株广泛分布于各个群体中，与主成分分析结果基本一致。全基因组关联分析共检测到12个与球果数量显著关联的SNP位点，其中，Marker193307的序列与纤维素合成相关基因Korrigan (KOR1，编码跨膜内切β-1,4-D葡聚糖酶)的序列相似度较高(74.74%)。qPCR基因组织特异性表达分析结果显示：该基因在华山松枝条韧皮部的表达量最高，幼嫩球果中次之，在根系的表达量最少，初步判断该基因与球果数量相关。【结论】华山松...     

玉米花期根系结构的表型变异与全基因组关联分析

【目的】根系作为植株吸收水分和养分的重要器官,对玉米生长及产量的形成至关重要。研究玉米根系结构的遗传机制指导玉米高产育种实践。【方法】以111份玉米优异自交系为材料,于2017年在北京、陕西永寿、山西定襄和河南原阳4个环境下对玉米地下节根层数（RLN）、地下节根总条数（TRN）、地下节根角度（RA）、地下节根面积（RS）、地下节根体积（RV）和地下节根干重（RDW）等6个玉米根系相关性状进行调查。取4个环境的平均值作为6个根系相关性状的表型数据,对6个相关性状进行统计分析和相关性分析,对不同年代、不同类群自交系的地下节根相关性状进行差异分析。基于该群体全基因组152 352个高质量SNP标记,利用FarmCPU模型进行全基因组关联分析获得显著关联SNP位点,并在LD衰减距离范围内查找候选基因,对候选基因的功能进行富集分析。【结果】表型分析表明,6个地下节根性状均呈现正态分布,且均显示出较高的遗传力;相关性分析结果表明,地下节根层数和总条数均与地下节根角度和面积呈负相关,地下节根的角度、面积、体积和干重等4个性状之间相互呈现显著正相关关系;不同年代的玉米地下根系结构存在差异,地下节根层数和总条数在年代的更替间表现出下降的趋势,地下节根角度和面积在年代更替间表现出上升的趋势,根干重和根体积在各年代间无显著差异;玉米地下根系结构在类群间也存在差异,旅大红骨类群的6个地下节根性状值均高于其余类群。全基因组关联分析共检测到26个SNP位点与地下节根层数、总条数、体积和干重性状显著关联（P<0.00001）,其中11个显著关联位点定位于前人报道的根系QTL区间内,2个显著关联SNP在地下节根层数和总条数中均被检测到。基于显著关联SNP位点共挖掘到177个候选基因,其中135个具有功能注释,Zm00001d037368可能为控制地下节根层数和总条数的一因多效候选基因。候选基因功能的富集分析结果显示,候选基因的功能主要涉及植物体内的代谢调节、应激反应、运输活性、催化活性、结合蛋白及细胞成分等。【结论】玉米自交系的根系结构在不同年代间和不同类群间存在不同程度的差异,采用全基因组关联分析策略挖掘控制玉米根系结构的相关遗传位点及候选基因,共检测到26个显著关联的SNP位点。     

粳稻分蘖数全基因组关联分析及候选基因的挖掘

【目的】利用全基因组关联分析(genome-wide association study,GWAS)检测与粳稻分蘖数显著相关的SNP位点,筛选影响分蘖数的候选基因,为分子辅助育种提高产量提供理论基础。【方法】利用295份来自世界各地的粳稻品种,于2018和2019年分蘖盛期调查水稻分蘖数,结合高通量重测序获得的788 396个高质量多态性SNP,利用TASSEL 5.0软件的MLM模型进行全基因组关联分析,对GEC软件计算的有效独立SNP数目进行阈值确定,判定SNP标记与目标性状关联的显著性。根据2年检测到的峰值SNP和水稻每条染色体LD衰减距离,确定2年共定位分蘖数主效QTL,并进一步提取QTL区间内所有基因外显子区非同义突变SNP和启动子区SNP进行单倍型分析,结合基因注释筛选影响粳稻分蘖数候选基因。【结果】295份粳稻品种的分蘖数在2018和2019年总趋势基本一致,并且均具有较大的表型分布。通过全基因组关联分析,在P5.46×10-6阈值条件下,从第8、9和10染色体上共鉴定3个与粳稻分蘖数相关的QTL(qTiller8、qTiller9和qTiller10),其中,在2年中均检测到qTiller9,在2018和2019年的表型贡献率分别为11.86%和10.61%,而qTiller8和qTiller10仅在2018年被检测到,表型贡献率分别为9.36%和9.10%。对qTiller9区间内的全部15个基因进行单倍型分析,结果表明,在qTiller9区间内共有6个基因(LOC_Os09g25090、LOC_Os09g25100、LOC_Os09g25150、LOC_Os09g25190、LOC_Os09g25200和LOC_Os09g25220)的不同单倍型分蘖数之间存在极显著差异,LOC_Os09g25090被启动子区SNP分为2种单倍型,hap2(TAA)分蘖数极显著大于hap1(AGG);LOC_Os09g25100被非同义突变SNP分为2种单倍型,hap2(GAGA)分蘖数极显著大于hap1(AGCC);LOC_Os09g25150被非同义突变SNP分为2种单倍型,hap2(ATG)分蘖数极显著大于hap1(GCC);LOC_Os09g25190被启动子区SNP分为2种单倍型,hap2(GCATCGCATCGACGCCGA)分蘖数极显著大于hap1(ATGCTGATGAAGTCATCC);LOC_Os09g25200被非同义突变SNP分为2种单倍型,hap2(TAG)分蘖数极显著大于hap1(AGA);LOC_Os09g25220被非同义突变SNP分为2种单倍型,hap1(GG)分蘖数极其显著大于hap2(AA)。结合基因注释发现,LOC_Os09g25090和LOC_Os09g25100均预测编码钙调蛋白依赖性蛋白激酶,是脱落酸(ABA)表达所必需Ca2+的传感器,推测LOC_Os09g25090和LOC_Os09g25100为影响粳稻分蘖数的候选基因。【结论】筛选出LOC_Os09g25090和LOC_Os09g25100为影响粳稻分蘖数的候选基因。     

柚果肉颜色遗传变异分析及候选基因挖掘

【目的】果肉颜色是柚类品种重要的外观性状以及品质性状,挖掘与柚果肉颜色显著相关的变异位点及相关基因,为进一步理解柚类品种果肉呈色机理以及分子标记辅助育种奠定基础。【方法】通过表型调查及色差仪测量,对100份柚类种质的果肉颜色进行评价和分类,利用GBS(genotyping-by-sequencing)测序技术对所有材料进行简化基因组测序。将测序得到的基因型数据通过GCTA软件计算特征值及特征向量分析群体结构,利用plink 2.0软件计算不同果肉颜色群体的遗传分化系数（Fst）,选用GEMMA软件中GLM模型进行全基因组关联分析,选取与颜色表型显著关联的变异位点进行等位变异分析,根据柑橘LD大小对变异位点邻近位置进行基因注释,筛选出与果肉颜色形成相关的候选基因,随机选择4份白肉柚和4份红肉柚材料在不同发育时期对候选基因表达进行初步验证。【结果】根据果肉颜色将100份柚类种质分为白肉柚和红肉柚两大类,其中包含58份白肉柚和42份红肉柚。通过Fst遗传分化指数分析和GWAS全基因组关联分析共筛选到Fst系数大于0.4且-log10(P)>9的位点6个,对6个变异位点进行基因分型,根据...     

浅谈饲料中有效能的测定技术

饲料中有效能是供动物生长发育的基础.不同动物所用的有效能体系不同,目前大多数动物采用消化能、代谢能体系,但随着研究的发展与深入,发现最能反映饲料有效能的是净能体系.无论哪种体系,采用合理的测定技术准确测定饲料中的有效能值显得尤其重要,通过对饲料有效能值的准确测定可以实现动物所需能量的精确供给,减少养殖成本,使经济效益最大化.文章综述了几种有效能评价体系的测定技术.     

高校贫困生认定工作面面观

国家贫困生资助政策实施以来,对贫困生帮助很大,同时在实际运行中还存在着一些问题。本文提出贫困生认定工作仍需要进一步采取各种相关配套措施,以推动和保障贫困生资助工作更好地开展。     

《河北农业大学学报》创刊年代考

清光绪二十八年(1902)河北农业大学前身—直隶农务学堂诞生,经几易其名,于1958年更名为河北农业大学至今。清光绪三十一年(1905)直隶高等农业学堂时期创办了《北直农话报》,清光绪三十四年(1908)更名为《直隶农务官报》,中华民国七年(1918)改出《农学月刊》,中华民国十七年(1928)易名为《河大农刊》,中华民国二十三年(1934)更名为《河北通俗农刊》,中华民国二十四年(1935)易名为《河北农林学刊》,1948年更名为《河北农学院研究专刊》,1959年更名为《河北农业大学学报》至今。《河北农业大学学报》前身诸刊都与现时的《河北农业大学学报》有着一脉相承的历史渊源,各刊之间联系紧密,连续性、继承性强。因此,《河北农业大学学报》的创刊时间应追溯至1905年创办的《北直农话报》。     

生态旅游潜在客源市场特征的调查研究--以湖南永州金洞森林旅游开发为例

为探明客源市场生态旅游消费的潜在特征,采用问卷调查的形式,就长沙市居民对湖南金洞生态旅游开发的意向等问题进行抽样调查．结果显示,生态旅游符合人们“回归自然”的旅游新时尚,有着极大的开发空间,指出生态旅游的开发要注重环境保护和可持续发展．开发的产品要以休闲度假类的大众产品为主,开发生态旅游都市客源市场还要多种渠道并用,尤其是要注重媒体的宣传．     

探究板栗的科学贮藏方法

板栗属壳斗科栗属(Castanea mollisima Blume),其种子属于顽拗型种子,不耐贮藏。基于近年来板栗贮藏保鲜技术研究成果,从合理采收、贮前处理、贮藏方法等方面进行论述。     

动物疫苗接种异常反应的处理

动物疫病是畜牧业生产的大敌,要发展畜牧业生产就要防治动物疫病。但在疫苗接种时,经常出现正常反应外的其他不利于机体的反应,如废食、皮疹、休克、死亡等,正确处理或避开这些问题,对推行动物疫病的计划免疫,实施强制免疫,保护人畜安全具有十分重要的现实意义和作用。     

厚胶合板弹性模量预测模型

为了预测厚胶合板弹性模量,通过简化层合板理论,该文建立了胶合板弹性模量预测模型,并以单板条、经过涂胶热压处理的单板条和相同工艺条件下的单板层积材的弹性模量,采用4种不同铺层方式的19层桉树胶合板对模型进行了验证。结果表明:3种预测值与实测值的趋势一致,相关系数R2顺纹在0.86以上,横纹在0.88以上,但是精度不同。采用单板条弹性模量预测的胶合板弹性模量比实测值偏低;采用经过涂胶热压处理后的单板条弹性模量预测的胶合板弹性模量比实测值偏高;采用相同工艺条件下的单板层积材弹性模量预测的胶合板弹性模量与实测值偏差较小,顺纹平均误差为4.64%,横纹平均误差10.94%。因此,采用相同工艺条件下的单板层积材的弹性模量来预测胶合板的弹性模量是可行的。      

增强我国水果业出口竞争力的思考

分析了入世4年来,我国水果出口在量上实现了突破,但并没有实现质变的主要原因,指出我国水果业存在的问题。论述了引进外资对发展我国水果业的意义,并提出了具体可行的对策与建议。     

山茱萸多糖提取工艺优化及马钱苷含量测定

采用L（934）正交设计试验,对山茱萸浸提液中山茱萸多糖的酶水解法提取工艺进行了优化研究,并对浸提液的中有效成分马钱苷含量进行了HPLC法分析。结果表明,山茱萸多糖浸提的最佳工艺为：液料比1∶5,浸提时间4 h,浸提温度80℃,果胶酶添加量0.55 g/L。用HPLC法测定出的山茱萸浸提液中马钱苷平均含量为0.512 ...     

甘肃省日光温室生产发展对策的探讨

甘肃省设施农业发展历史悠久,古代创造了麦草覆盖生产韭 黄及泥碗护苗等传统设施农业技术,至今仍然受到农民欢迎,在 甘肃省中部应用面积约5万多hm2。建国后,甘肃设施农业获得 了新生,大致经历了三个阶段;第一阶段为引进应用北京改良式 温室阶段,时间在20世纪50-60年代;第二个阶段为塑料拱棚 与地膜覆盖栽培阶段,时间在20世纪70-80年代:第三阶段为