甘蓝型优质油菜籽粒含油量和千粒质量相关性分析

检测了2004—2005年国家油菜区域试验8个区的甘蓝型油菜品种(系)的籽粒含油量与千粒质量,并进行相关性分析,结果表明,长江上游区(1,2组)、长江中游区(A,B,R组)、黄淮区和长江下游区(C,D组)共计8个区的甘蓝型油菜品种(系)的籽粒含油量与千粒质量的相关系数分别为-0.011 3,0.357 9,-0.268 4,0.194 6,0.150 0,0.523 0,-0.326 6和0.205 3,即除了黄淮区的油菜籽粒含油量与千粒质量呈中度相关外,其它7个区都呈弱度相关或基本无相关性。     

油菜耐盐相关性状的全基因组关联分析及其候选基因预测

【目的】通过对甘蓝型油菜耐盐相关性状进行全基因组关联分析,寻找可能与油菜耐盐性相关的SNP位点,发掘与油菜耐盐性有关的候选基因。【方法】以1.2%NaCl溶液作为培养液,去离子水为对照,对307个不同品系的甘蓝型油菜进行发芽试验。播种后7 d测定幼苗根长、鲜重及发芽率,计算盐胁迫下各性状相对值,并作为评价耐盐的指标。结合油菜60K SNP芯片,利用SPAGeDi v1.4软件对该群体307份甘蓝型油菜进行亲缘关系分析,并计算亲缘关系值的矩阵。利用软件STRUCTURE v2.3.4对关联群体进行了群体结构分析。为有效排除假关联的影响,将群体结构和材料间的亲缘关系考虑进关联分析中,同时进行了PCA+K模型、Q+K模型以及K模型3种混合线性模型分析和比较,根据所有SNP的–lg(P)观察值和期望值,确定每个性状GWAS分析的最优模型。采用TASSEL 5.0软件,在最优模型下对307份材料耐盐各性状的相对值分别与SNP标记进行全基因组关联分析。利用油菜基因组数据库,在显著SNP位点侧翼序列200 kb范围内提取基因。根据拟南芥中已经明确功能的耐盐相关基因,筛选出目标基因组区段内与耐盐相关的油菜同源基因。【结果】全基因组关联分析共检测到164个与根长显著关联的SNP位点,23个与鲜重显著关联的SNP位点,38个与发芽率显著关联的SNP位点。其中与根长、鲜重、发芽率最显著关联的SNP位点分别位于染色体A08、A02和A06,贡献率分别为23.84%、18.59%和31.81%。在这些显著SNP位点侧翼序列200 kb范围内发掘出可能与油菜耐盐性有关的50个候选基因。这些候选基因主要包括转录因子MYB、WRKY、ABI1、b ZIP、ERF1、CZF1、XERICO等以及一些下游受转录因子调控的不同功能基因NHX1、PTR3、CAT1、HKT、CAX1、ACER、STH、STO等。在根长和发芽率2个不同耐盐性状的分析结果中均筛选出位于A03染色体上的耐盐基因BnaA03g14410D。另外,这些耐盐候选基因中包含两组串联重复基因,分别是位于A03染色体上的BnaA03g18900D和BnaA03g18910D,位于C09染色体上的BnaC09g19080D、BnaC09g19090D和BnaC09g19100D。除此之外,耐盐候选基因中还包含2个距离非常近(中间只间隔2个基因)的重复基因BnaC02g39600D和BnaC02g39630D。【结论】共检测到225个与耐盐性状显著关联的SNP位点,筛选出50个可能与油菜耐盐性有关的候选基因。     

甘蓝型油菜品质相关性状全基因组关联分析

采用近红外方法,分析104份甘蓝型油菜种子中的芥酸、油酸与硫代葡萄糖苷含量;运用简化基因组技术(ddRADseq)对全部材料进行基因分型;结合性状与基因型的数据,利用Tassel混合线性模型进行全基因组关联分析。结果表明:芥酸含量与硫代葡萄糖苷呈极显著正相关,油酸含量与芥酸、硫代葡萄糖苷呈极显著负相关;群体结构分析将104份材料划分为4个亚群;通过全基因组关联分析发现,与芥酸、油酸含量显著关联的SNP标记主要分布在A08与C03染色体上,与硫代葡萄糖苷含量显著关联的标记分布在A05、A02与C09染色体上。     

甘蓝型油菜一次分枝角度主效基因染色体定位和候选基因挖掘

一次分枝角度是构建甘蓝型油菜理想株型的重要性状之一.本研究以分枝角度大的纯系黄花亲本62与分枝角度小的纯系粉花亲本77为亲本构建114个双单倍体(Double Haploid, DH)群体株系,通过基于测序的基因分型(Genotyping by Sequence, GBS)技术简化测序获得单核苷酸多态性(Single Nucleotide Polymorphism, SNP)数据,考察DH群体一次分枝角度表型,利用一般线性模型(General Linear Model, GLM)进行全基因组关联分析.结果表明：该群体GLM模型检测到8个显著SNP位点都分布在A07染色体,解释了25.21%～35.05%的表型变异.定位区间内分析和注释了49个基因,通过京都基因与基因组百科全书(KEGG)通路富集分析推测甘蓝型油菜一次分枝角度主效基因为编码GH3家族蛋白的BnaA07g01910D基因,此基因负责细胞繁殖扩增可能影响分枝角度表型.该研究结果有助于解析甘蓝型油菜分枝角度的遗传基础和调控机制,为指导分枝角度遗传改良奠定了基础.     

盐胁迫下甘蓝型油菜发芽期下胚轴和根长的全基因组关联分析

【目的】解析甘蓝型油菜发芽期根和下胚轴发育及耐盐性的调控位点,筛选油菜耐盐性相关的候选基因,可为油菜耐盐性改良提供依据。【方法】以317份具有代表性的甘蓝型油菜自交系为材料,在正常生长和盐胁迫条件下进行沙培鉴定,利用芸薹属60K SNP芯片和全基因组关联分析鉴定正常生长与盐胁迫下甘蓝型油菜发芽期根和下胚轴长度显著关联的SNP,并确定其连锁不平衡区间。通过区间内基因的功能注释及盐胁迫下油菜幼苗根和叶片转录组差异表达基因筛选连锁不平衡区间内的重要候选基因,并以实时荧光定量PCR分析候选基因的组织特异性和盐胁迫诱导表达模式,提高候选基因筛选的准确性。【结果】正常生长和盐胁迫下甘蓝型油菜发芽期下胚轴和根长在不同材料间变异较大,频次分布表明目标性状均为数量性状,受多基因调控。全基因组关联分析模型比较表明,MLM+P+K模型为最优模型。以此模型对目标性状进行全基因组关联分析,检测到45个显著关联SNP,其中40个与下胚轴长度显著关联,5个与根长显著关联,单个SNP解释的表型变异分别为9.12%—14.46%和7.67%—8.93%。重复检测的显著相关SNP中,值得注意的是C04染色体的rs8970,同时与4个性状显著关联,表型贡献率为7.67%—12.35%,是唯一在下胚轴长和根长间重复检测到的显著关联SNP。11个重要关联SNP中有6个位于10—442 kb的连锁不平衡区块中。转录组分析表明,11个连锁不平衡区间共包含447个基因,其中15个受盐胁迫诱导表达。转录组和基因功能注释综合分析表明,BnaSRO1、BnaPAGR2、BnaNPH3、BnaMYB124、BnaSAM-Mtase、BnaBIN2、BnaUMAMIT11、BnaEXPA7、BnaRPT3、BnaEF-hand和BnaF3H很可能为各自区间的候选基因。实时荧光定量PCR结果证实除BnaNPH3外,其他基因均在根或下胚轴中受盐胁迫诱导上调表达。组织特异性分析还发现BnaUMAMIT11、BnaPAGR2和BnaEXPA7主要在萌发的根和下胚轴中特异表达,BnaRPT3、BnaBIN2和BnaMYB124虽然呈组成型表达,但在萌发阶段的下胚轴中表达量最高,证实这些基因很可能参与油菜发芽期根和下胚轴生长发育及耐盐性的调节。【结论】全基因组关联分析共鉴定出45个控制油菜发芽期根和下胚轴发育及耐盐性的显著关联SNP。连锁不平衡、转录组和基因功能注释综合分析初步鉴定出11个重要候选基因。     

甘蓝型油菜籽粒产量和品质性状对氮肥用量的响应

为探索甘蓝型油菜籽粒产量和品质性状对氮肥用量的响应,以宁杂1818和秦优10号为试验材料,设0 kg/hm~2(N0)、120 kg/hm~2(N120)、240 kg/hm~2(N240)和360 kg/hm~2(N360)4个氮肥水平处理,于2014-2016年度在南京市、扬州市和泰州市进行大田试验。结果表明,油菜籽粒产量、产油量、含油率和蛋白质含量对氮肥用量的响应均达到极显著水平,油酸和芥酸含量的响应较小,多不饱和脂肪酸和硫苷含量的响应因品种、试点和年度的不同均有明显差异,部分试点达到显著水平。籽粒产量和产油量对施氮量的响应符合一元二次方程。研究结果还表明籽粒产量、产油量和蛋白质含量与施氮量呈极显著正相关,多不饱和脂肪酸含量与施氮量呈显著正相关,含油率与施氮量呈极显著负相关,硫苷含量与施氮量呈显著负相关。试验中发现每增施120 kg/hm~2氮肥, 3个施氮处理的籽粒增产效率依次为81.11%、34.12%和-2.24%,籽粒增产效率随施氮量增加逐渐下降,品种间表现一致。在南京、扬州和泰州试点,每1 kg氮肥供给平均可增收菜籽11.94 kg、11.35 kg和12.36 kg,最适施氮量宁杂1818分别为300～310 kg/hm~2、300～320 kg/hm~2和270～300 kg/hm~2,秦优10号分别为315～335 kg/hm~2、295～315 kg/hm~2,300～310 kg/hm~2,最适施氮量下不仅产量高而且品质优。     

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

【目的】挖掘与油菜角果长度性状显著相关的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两种分析方法探测到多个与油菜角果长度性状显著关联的基因位点,并在显著性位点附近搜寻到相关候选基因。     

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

【目的】利用全基因组关联分析(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年总趋势基本一致,并且均具有较大的表型分布。通过全基因组关联分析,在P<5.46×10-6阈值条件下,从第8、9和10染色体上共鉴定3个与粳稻分蘖数相关的QTL(qTiller8、qTiller9和qTiller10),其中,在...     

甘蓝型油菜籽粒RNA提取方法的探讨

用CTAB法、SDS法、上海华瞬植物叶RNA抽提试剂盒分别提取甘蓝型油菜花后40 d的籽粒RNA.用电泳法和A260/A280与A260/A280的比值比较提取RNA的质量.结果发现:用SDS法、上海华瞬植物叶RNA抽提试剂盒都没有成功提取RNA,只有CTAB法成功地提取出了完整的甘蓝型油菜籽粒RNA.用CTAB法提取的总RNA的A260/A230与A260/A280的值分别为2.06和1.94,说明采用该方法提取的RNA污染小,纯度较高,获得的RNA产量也较高.通过逆转录实验和基因片段的克隆也证明获得的RNA质量比较高,可以直接满足一般的分子实验.     

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

【目的】挖掘与叶绿素含量显著关联的单核苷酸多态性(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个可能与叶绿素含...     

黄瓜幼苗下胚轴长度GWAS分析及候选基因挖掘

【目的】挖掘与黄瓜幼苗下胚轴长度显著相关的SNP位点及候选基因,为揭示下胚轴长度的遗传基础和分子机制提供理论依据,为短下胚轴分子标记辅助选择育种奠定基础。【方法】以95份黄瓜核心种质为试验材料,分别于2016年春季、2017年春季、2017年秋季和2018年春季在中国农业科学院南口试验基地塑料大棚进行种植,在两叶一心期调查黄瓜幼苗的下胚轴长度;利用Structure 2.3.4软件分析群体结构,Haploview软件分析连锁不平衡的衰减;基于最优模型对下胚轴长度进行全基因组关联分析(GWAS),依据关联SNP位点的LD区间序列,预测与下胚轴长度相关的重要关联候选基因,并利用荧光定量PCR对预测基因进行表达模式分析。【结果】共检测到8个显著关联的位点(Hl1.1、Hl1.2、Hl2.1、Hl3.1、Hl3.2、Hl4.1、Hl5.1、Hl6.1),分别位于1、2、3、4、5、6号染色体,其中,Hl2.1、Hl3.1、Hl3.2、Hl5.1、Hl6.1等5个位点被重复检测到两次以上。通过分析关联SNP位点的LD区间序列,获得Csa1G074930、Csa1G475980、Csa2G3816...     

甘蓝型油菜polCMS育性恢复位点的全基因组关联分析

【目的】甘蓝型油菜波里马细胞质雄性不育(pol CMS)在中国已被广泛应用于杂交种育种,其育性恢复程度表现出受1对主效基因的控制,并受微效修饰基因的影响。通过全基因组关联分析方法挖掘育性恢复位点,并对候选基因进行比较分析。【方法】通过芸薹属60K SNP芯片对308份甘蓝型油菜自然群体进行基因型分型,并用pol CMS系301A作母本,与上述材料分别进行杂交得到308份F1,每份F-_1分别于2013年和2014年进行种植,每年2次重复,于始花期根据花粉育性和花蕊发育情况调查F1植株的育性等级,同时对测交父本自然群体进行群体结构分析和亲缘关系评估,并结合测交父本的基因型分型结果和F_1的育性等级进行全基因组关联分析(GWAS)。从GWAS分析中显著的SNP左右100 kb区间或与显著SNP处于同一单体型块(R~20.5)的区间内预测候选基因,并对候选基因进行QTL比较分析和单体型或等位基因的效应分析。【结果】方差分析结果显示,两年F_1的育性等级存在显著差异(P0.01),但相关分析发现,两年的育性等级存在显著的正相关(r=0.52,P0.001)。群体结构分析显示,所有测交父本被分为3个亚群(冬性、春性和半冬性),亲缘关系分析发现,任何2个材料之间平均亲缘关系值为0.072,73%的任意材料间亲缘关系值小于0.1,其中,约53%的材料亲缘关系值为0。GWAS分析共检测到13个与育性恢复程度显著关联的SNP,构成了6个候选区间,分别位于A01、A09、C03、C06和C08 5条染色体上,单个SNP解释的表型变异介于2.53%—9.96%。从中共预测到6个与育性恢复位点相关的候选基因,其中4个编码的蛋白含有恢复基因特有的PPR保守基序。共线性分析发现,4个候选基因中的2个(Bna A09g46700D和Bna C08g40710D)位于A09和C08染色体部分同源区间,且与已克隆的pol CMS育性恢复位点ORF2同源。另外2个新鉴定到的候选基因(Bna C03g45840D和Bna C06g13000D)连锁的SNP等位基因或单体型变化都与育性等级显著相关(P0.001)。【结论】通过GWAS分析鉴定到多个与油菜育性恢复有关的候选基因,开发基于与这些基因连锁位点或SNP的功能标记将有助于对该不育系统进行恢复系和保持系的筛选。     

甜玉米籽粒体积和粒重的全基因组关联分析

为挖掘控制甜玉米籽粒体积和粒重的相关基因,2017—2019年,以209份甜玉米自交系构成的关联群体为材料,测定干籽粒的体积和粒重,结合全基因组重测序(De novo sequencing)获得的980万个单核苷酸多态性(Single-nucleotide polymorphism,SNP)标记,采用GLM模型进行全基因组关联分析。结果表明,甜玉米籽粒体积变异范围为0.06~0.15 mL,变异倍数为2.5倍,单粒重变异范围为0.08~0.17 g,变异倍数为2.1倍,均符合正态分布。共检测到15个SNP位点与籽粒体积或粒重显著关联,分别解释籽粒体积55.7%和粒重52.1%的表型变异,其中4个位点和普通玉米报道的QTL重叠,其余11个是新鉴定的位点,通过基因注释发现大部分候选基因为转录因子等调控基因。此外,首次检测到Br2(矮化基因2)与籽粒体积、粒重两性状均显著关联,并且籽粒体积和粒重均随着该基因表达量的升高而增加。     

小麦籽粒镉元素含量全基因组关联分析及候选基因预测

以国内外207份小麦种质为材料,利用660K SNP芯片对其进行基因型检测,并结合不同环境下表型数据和最佳线性无偏预测值 (BLUP,Best linear unbiased prediction) 对小麦籽粒镉元素含量进行全基因组关联分析。结果表明：与小麦籽粒镉元素含量显著关联的SNP 310个,这些SNP分布于除3D和4D外的19条染色体上,单个SNP解释变异率为10.95%～14.66%。不同环境下检测到的关联SNP结果存在差异,其中在原阳地区检测到186个SNP,开封地区检测到71个SNP。基于BLUP值分析获得53个SNP。基于SNP物理位置,将距离较近的SNP进行整合,共获得有效QTL位点52个。同时发现了7个在多环境下表现稳定的SNP,并对其进行单标记效应分析。最后对基于获得的关联SNP进行了候选基因预测,共获得7个与小麦籽粒镉元素含量相关的候选基因,其中 TraesCS1B01G321700和TraesCS1B01G320200可能与镉元素调控相关基因转录有关,而TraesCS7B01G459000和TraesCS7B01G456900可能与镉元素的吸收和转运等代谢过程有关。还筛选出了对镉具有良好避性的部分小麦优异种质,如‘云麦51’‘郑麦379’‘白穗白’‘云麦53’‘双丰收’。     

高温胁迫下油菜籽粒成熟期基因差异表达分析

利用数字化基因表达谱技术,对花后不同时期及不同环境温度条件下差异表达基因进行统计分析.将3份 甘蓝型油菜于开花后20d移入到40 ℃高温和低温2个环境,分别处理12d和25d后取角果种子提取RNA,对 花后不同时期及不同环境温度条件下差异表达基因进行统计分析.结果再以低温条件为对照,在高温下花后32d 共检测出968个差异表达基因,其中上调表达基因为501个,下调表达基因为467个;在花后45d检测到1064 个差异表达基因,其中上调表达基因为274个,下调表达基因为790个.结果表明,随着角果成熟,绝大多数基 因的表达量减少,少数表达量显著上升的基因值得关注,它们可能与角果储藏物质的合成和角果成熟有关;在籽 粒充实的高峰期(花后32d),高温胁迫下上调表达的基因较多,其中可能有一部分基因与环境应答有关.本研究 结果与前人的结果在一定程度相似,表明本研究所采用的研究方法是可行的,获得的与环境胁迫相关基因为后 续研究奠定了基础.     

加权关联共表达网络在筛选致病候选基因上的优势分析

面对多基因疾病的致病基因筛选的难题,加权关联共表达网络模型作为一种新的基于图论和统计学的生物信息挖掘方法,较其他传统方法以其较好的预测精确度,完备的生物学理论以及实现的简洁可行,有着较好应用前景。重点将加权关联共表达网络与基于网络距离和基于疾病表性相似性这2大主流方法进行实例对比分析,以体现加权关联共表达网络在致病候选基因筛选上的优势性。     

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

【目的】通过对芝麻产量相关性状的全基因组关联分析,挖掘与产量性状关联的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个可能与产量相关的重要候选基因。     

母猪初情期全基因组关联分析和位置候选基因研究

通信作者陈从英,Tel/Fax：0791-83813080;E-mail：chcy75@  hotmail.com【目的】分离影响母猪初情期的主效基因或分子标记。【方法】以白色杜洛克×二花脸资源群体F2代母猪为研究材料,利用Illumina猪60K SNP芯片分别对F0和F1及316头有初情期表型记录的F2代母猪进行基因型判定,通过单标记全基因组关联分析（GWAS）和连锁-连锁不平衡（LDLA）分析检测与母猪初情期显著关联的SNP位点或单倍型。采用标准关联分析、标记辅助关联分析和F值下降检验分析lin-28同源B基因（LIN28B）和跨膜蛋白38B基因（TMEM38B）3个SNP位点与母猪初情期的关联性。【结果】①与母猪初情期关联性最强的SNP为ASGA0032316,位于7号染色体（SSC7）33.07 Mb处RAB23基因内含子中,同时在1、6、12、15和17号染色体也检测到与母猪初情期显著关联的SNP;②达基因组显著水平的单倍型均位于SSC7,其中关联性最强的单倍型位于SSC7的38.39 -38.47 Mb处F1RVR7和ZFAND3基因之间的基因间隔区;③LIN28B和TMEM38B基因的3个SNP与母猪初情期均未达到显著相关（P＞0.05)。【结论】在白色杜洛克×二花脸资源群体中,与母猪初情期关联性最强的SNP位点和单倍型均定位于7号染色体, 1、6、12、15和17号染色体也定位到与母猪初情期显著关联的SNP,LIN28B和TMEM38B基因不是1号染色体QTL区间内的因果基因或需要搜寻更多的SNP进行分析验证。     

整合数字基因表达谱与全基因组关联分析鉴定猪血液性状候选基因

【目的】整合数字基因表达谱与全基因组关联分析鉴定白色杜洛克×二花脸F₂资源群体的血液性状候选基因。【方法】白色杜洛克×二花脸F₂资源群体在(240±3)d屠宰,收集血液于抗凝管中进行血常规检测。利用Illumina 60K SNP芯片对1 020头F₂资源群体进行基因分型。剔除基因型检出率< 90%和孟德尔错误检出率> 5%的个体。检出率< 95%、次等位基因频率< 5%、哈代-温伯格检验(HWE) P < 5×10^-6、与性染色体连锁疑似常染色体的SNP被筛除。利用Illumina GA II 测序仪测序对502头F₂资源群体的肝脏进行数字基因表达谱测序。测序得到的原始数据经过滤获得清洁标签后与参考标签数据库比对,将能唯一比对到参考基因序列的清洁标签数量进行标准化处理以获得标准化的基因表达量。每个转录本的表达水平进一步转化为lg2值。在少于20%的个体中表达的转录本被滤去。表型性状和基因表达性状使用R程序包中GenABEL内polygentic功能进行性别、批次和亲缘关系的校正。其残差使用R程序包中斯皮尔曼系数评估基因表达水平与表型数据的关联性,设定保守阈值P < 5×10^-4时调整多重检验。将检测到的表达数量性状位点(eQTL)及其对应基因根据其位置相对照的关系进行绘图。搜寻前期GWAS最高点5.0 Mb区域内eQTL结合GWAS结果进行综合分析。Gene Ontology & KEGG pathway富集分析使用在线工具DAVID。基因共表达网络使用在线工具GeneMANIA进行构建。【结果】白色杜洛克×二花脸F₂资源群体中502个个体的20 108个肝脏转录本通过了质检。当P < 5×10^-4时鉴别到与血红蛋白(HGB)、红细胞数目(RBC)、红细胞压积(HCT)、平均红细胞体积(MCV)、平均红细胞血红蛋白含量(MCH)和白细胞数目(WBC)关联的转录本共259个。有34个转录本与一个以上表型关联。用上述血细胞性状关联的转录本进行eQTL定位,当阈值P < 10^-5时得到304个eQTL,每个转录本映射到1-6个eQTL,其中有35个顺式eQTL,120个反式eQTL。MCH和MCV的顺式eQTL位置重叠位于8号染色体。7号染色体存有数量最多的eQTL,其中多数为反式eQTL。通过eQTL定位确定了KIT为候选基因。Gene Ontology & KEGG pathway富集分析鉴别到与红细胞性状相关的基因KIT、PSEN2和TFRC,与白细胞性状相关的基因THBS1和CYR61。通过整合前期GWAS数据及eQTL定位并建立基因共表达网络,鉴定到与白细胞性状相关的基因RPS10。【结论】利用基于白色杜洛克×二花脸F₂资源群体的eQTL定位并结合前期GWAS数据,鉴别到与红细胞性状相关的基因KIT、PSEN2和TFRC,与白细胞性状相关的基因THBS1、CYR61和RPS10。