通过转录组测序技术挖掘蛋白质合成相关基因,为解析蛋白质合成机制奠定基础。以遗传背景相近但蛋白含量差异较大的大豆高蛋白品系冀HJ117(蛋白质含量52. 99%)及其回交亲本冀豆12(蛋白质含量46. 48%)为研究材料,利用转录组测序技术和生物信息学分析,以期挖掘大豆籽粒蛋白质合成相关基因。通过对转录组数据中冀HJ117和冀豆12差异表达基因的分析筛选,共得到336个差异表达基因,其中冀HJ117较冀豆12有195个上调表达基因,141个下调表达基因。通过GO功能富集分析,发现在分子功能类型中注释的差异基因主要与催化和连接等功能相关; KEGG显著富集分析发现差异表达基因主要富集在蛋白质内质网合成途径中,该途径共筛选到34个差异表达基因,其中33个基因在冀HJ117中表达量较高。从该途径中筛选出9个候选基因,采用荧光定量PCR方法检测其表达量,表达趋势与RNA-Seq检测结果基本一致,证实了RNA-Seq数据的可靠性。RNA-Seq测序结果获得了与蛋白质合成相关基因的基本信息,蛋白质内质网合成途径可能是冀HJ117与冀豆12蛋白质含量产生差异的重要通路。

为了寻找仿刺参(Apostichopus japonicus)养殖期间的"化皮病"关键调控基因,并分析这些基因所参与的信号通路,对本课题组前期已获得的仿刺参"化皮"I期(早期)、II期(中期)和III期(后期)3个阶段的病变及其同一个体正常体壁组织之间的差异表达基因(differentially expressed genes,DEGs)进行进一步分析。主成分分析(principal component analysis,PCA)结果显示,"化皮"II期病变组织与正常组织之间的差异最小,"化皮"I期与III期表达关系比较接近,"化皮"II期是一个"转折期"。KEGG富集分析结果显示,补体与凝血级联(Complement and coagulation cascades)通路和细胞外基质受体(ECM-receptor interaction)通路在"化皮"3个阶段都显著改变。通过构建"化皮"过程关键差异表达基因调控网络,发现Ig GFc-binding protein(Fc GBP)基因和Tenascin(TN)蛋白家族基因在"化皮"不同阶段参与到发生显著变化的信号通路。q RT-PCR验证结果显示,5个DEGs在仿刺参"化皮"不同阶段表达趋势与RNA-Seq结果一致,皮尔逊相关系数r值为0.7714。"化皮"过程关键调控基因的筛选将为抗逆品种选育以及"化皮病"的防控提供科学依据。

为了充分挖掘野生大豆种质资源中的高蛋白基因及其连锁标记,以来自中国、韩国和日本,涵盖第4,5,6,7,8熟期组的508份野生大豆种质资源为材料,通过全基因组关联分析挖掘与野生大豆中高蛋白基因相关的SNP。参试材料蛋白含量数据从美国农业部种质资源信息网下载,为2 a利用凯氏定氮法测定蛋白含量数据平均值,基因型数据从Soybase网站下载,利用Illumina公司大豆50K芯片(SoySNP50K BeadChip含有52 041个SNP标记)检测获得。结果表明,参试材料蛋白含量呈正态分布,介于38.1%～56.9%,平均48.1%,标准差2.71%。遗传结构分析将参试材料划分为3组,分别包含271,111,126份材料。基于混合线性模型的关联分析,共检测到与蛋白含量相关的SNP位点74个,散布在19条染色体的60个单倍型区段内。显著性SNP位点LOD平均值为3.47,SNP位点BARC_1.01_Gm_01_54656209_A_G的LOD值最大,为5.18。根据显著性SNP位点富集程度,确定第11号染色体常染色质区15 128 832～15 253 199 bp、第12号染色体异染色质区26 842 687～27 818 244 bp的单倍型区段为本研究中的2个蛋白含量显著性相关区段,命名为HAP_11_1和HAP_12_1。HAP_11_1中,SNP位点BARC_1.01_Gm_11_15167305_G_A的LOD值最大,为3.80,可解释遗传变异为2.88%。HAP_12_1中,SNP位点BARC_1.01_Gm_12_27563620_C_T的LOD值最大,为4.12,可解释遗传变异为3.23%。为野生大豆高蛋白基因育种利用提供了检测标记,为野生大豆高蛋白基因克隆提供了线索。

【目的】获取芋（Colocasia esculenta）全长转录本序列和结构信息，并挖掘淀粉生物合成相关基因和转录本的结构信息，为后续阐明芋淀粉生物合成的分子机制及深入利用芋基因资源提供科学依据。【方法】以荔浦芋1号为试材，采集其3月龄植株的不同组织（叶片、叶柄、球茎、匍匐茎和根）开展混合样三代全长转录组测序，对获得的转录本进行功能注释，通过生物信息学方法分析转录本可变剪接（AS）、可变多聚腺苷酸化（APA）、长链非编码RNA（lncRNA）、转录因子（TF）和简单重复序列（SSR）等的结构信息，通过京都基因和基因组百科全书（KEGG）注释获得芋淀粉生物合成相关的基因和转录本，并对其进行AS和APA分析。【结果】通过对芋叶片、叶柄、球茎、匍匐茎和根混合样开展三代全长转录组测序，共获得全长转录本序列38043条；通过与参考基因组进行比对，鉴定出新转录本31058条，其中28785条获得功能注释；通过对转录本结构分析，共鉴定出9360个AS事件、6436个基因存在poly(A)位点、25081个SSR位点、304个lncRNA、1911个融合转录本和1608个TF。通过KEGG注释获得芋淀粉生物合成相关基因14个和转录本60条；对获得的转录本进行AS和APA分析发现，6个基因发生了AS事件，包括内含子保留（IR）、5'端可变剪接（A5SS）、3'端可变剪接（A3SS）、外显子跳跃（ES）和互斥可变外显子（MEE）5种类型，有10个基因存在poly(A)位点，其中7个基因存在2个及2个以上poly(A)位点。【结论】通过三代全长转录组测序分析获得芋全长转录本序列和结构信息，挖掘到芋淀粉生物合成相关基因14个，转录本60条，可作为后续阐明芋淀粉生物合成分子机制及深入利用芋基因资源的科学依据。

【目的】基于全基因组重测序结果,开发与高蛋白、耐荫、抗倒伏等性状紧密相关的分子标记,同时利用开发的分子标记构建遗传连锁图谱,并对籽粒蛋白质含量进行QTL定位,为后续高蛋白、耐荫、抗倒育种研究提供参考和分子标记资源。【方法】以大面积栽培品种南豆12和地方品种十月黄为亲本,构建F2分离群体。对亲本材料进行覆盖度约为40×的全基因组重测序,用BWA、GATK及Breakdancer等软件比对,检测亲本材料在全基因组范围内的突变类型,挖掘相关变异基因。结合种子不同发育时期和荫蔽处理获得的转录组数据,结合qRT-PCR对发生突变的储藏蛋白、环境适应相关基因进行表达规律分析。同时,基于重测序数据,挖掘亲本间存在于基因编码区的SNP位点,对其进行酶切位点分析,将SNP标记转化为CAPS或dCAPS标记。此外,搜索亲本间存在的插入/缺失变异位点,在插入/缺失位点两侧高度保守的区域设计引物开发InDel标记。对开发的CAPS标记和InDel标记进行多态性筛选,选取具有多态性的CAPS分子标记和InDel标记,对F2材料进行基因分型。根据分型结果,利用JoinMap 4.0软件进行遗传连锁图谱的构建。依据构建的遗传图谱,结合近红外分析获得F2材料的籽粒蛋白质含量数据,使用Windows QTL Cartographer V2.5软件对大豆籽粒蛋白质含量进行QTL分析。【结果】测序结果显示,南豆12大量储藏蛋白、环境适应相关的重要基因或同源基因发生突变。转录组数据分析结果显示部分变异基因呈现不同的表达模式且差异显著,qRT-PCR分析进一步验证了该结果。此外,经检测开发的540个CAPS分子标记中有332个具有酶切多态性,300对InDel引物中有201对引物能扩增出多态性。基于533个多态性分子标记构建了一张包含20个连锁群的遗传连锁图谱,覆盖长度2 973.87 cM,标记间平均遗传距离5.58 cM。利用此图谱对大豆籽粒蛋白质含量进行QTL定位,共检测到QTL位点6个,可解释4.68%—18.25%的表型变异。【结论】基于亲本间的变异位点,共开发了533个多态性分子标记(包含8个基因特异性分子标记),检测到6个大豆籽粒蛋白质含量QTL位点,其中,主效QTL位点1个(qSPC-6)。

【目的】南方大豆皱叶症（southern soybean crinkle leaf disease,SSCLD）严重时可导致大豆减产40%左右，利用高世代分离群体转录组测序技术（high-generation segregating populations RNA-seq,HGRNA-seq）挖掘控制SSCLD的候选基因，为揭示SSCLD发生的分子机制提供数据支撑。【方法】以2个南方皱叶症症级差异材料及其衍生的F2:7株系为研究材料，对双亲进行重测序，对12个F2:7单株（7个皱叶，5个正常叶）进行单个样本转录组测序，利用转录组及亲本的SNP/InDel数据进行混池法定位分析，利用转录组表达量差异数据进行GO和KEGG功能注释和富集分析，并在定位区间附近开发7个KASP分子标记，利用230个F2群体构建局部连锁图谱，对转录组数据定位结果进行验证。联合基因定位和转录组分析结果，筛选控制SSCLD的候选基因。【结果】利用混池法把控制皱叶症的基因位点CL12定位在大豆第12染色体末端39 231 651—40 705 115 bp的1 473 464 bp区间内，利用F2群体将控制皱叶症的基因定位于39 743 275—40 948 295 bp的1 205 020 bp区间内，与混池法定位结果基本一致。GO注释结果显示，代谢过程包括免疫系统过程，以及对刺激的反应等，细胞组分主要与膜等相关，KEGG注释结果显示，生物系统通路中主要包括植物-病原菌互作和环境适应等通路。GO富集的表达差异基因主要同跨膜受体蛋白活性、蛋白磷酸化及信号受体活性等方面相关，KEGG富集到的最多差异表达基因（differently expressed genes,DEGs）主要在植物-病原菌互作和植物MAPK信号通路上。结合南方大豆皱叶症诱因特点，选择定位候选区间内与抗病等相关且在外显子上存在非同义突变或表达量有差异的基因作为候选基因，结合qRT-PCR验证，最终确定GLYMA＿12G223100、GLYMA＿12G223900、GLYMA＿12G224100、GLYMA＿12G231800和GLYMA＿12G233000等5个基因为控制南方大豆皱叶症的候选基因。【结论】HGRNA-seq实现了RNA-seq和BSA-seq相结合，成功挖掘了控制SSCLD的候选基因。

为挖掘我国名优鱼类长吻（鱼危）(Leiocassis longirostris)生长相关基因，本研究运用Illumina高通量测序技术比较分析了快速生长组[平均体质量为(534.02±53.68)g]和缓慢生长组[平均体质量为(108.41±4.96) g]各9尾长吻（鱼危）的脑组织基因表达谱。测序共获得267 404 674个高质量测序片段(clean reads)，通过2种不同生长速率长吻（鱼危）脑组织转录组比较筛选出518个差异表达基因，其中，412个基因表达量上调，106个基因表达量下调。对12个差异表达基因进行实时荧光定量PCR验证的结果与转录组测序结果一致。GO功能分类显示，大量差异表达基因富集到生长(growth)、生长因子活性(growth factor activity)和激素介导的信号通路(hormone-mediated signaling pathway) GO条目中。KEGG富集分析显示，一些差异表达基因在MAPK信号通路(MAPK signaling pathway)、转化生长因子β信号通路(TGF-beta signaling pathway)、钙离子信号通路(calcium signaling pathway)和神经活性配体–受体相互作用(neuroactiveligand-receptorinteraction)等途径中富集。根据GO功能注释和KEGG富集分析，筛选出gnrh、thr、egr1、fgf18、sst、gipr、cart和crf等基因是调控长吻（鱼危）生长发育的关键候选基因。本研究结果为后续深入研究长吻（鱼危）生长调控机制提供了重要的参考资料。

旨在通过对布莱凯特黑牛大小睾丸转录组数据进行分析，挖掘影响睾丸生殖功能的关键候选基因，探究睾丸大小对牛生殖功能维持的分子机制。采集12月龄布莱凯特黑牛睾丸，根据质量、长径和短径将其分为2组(大小睾丸组各3头)进行高通量转录组测序；以|log_2(Fold Change)|≥1且P-value≤0.01为阈值筛选差异基因，对差异基因进行GO和KEGG功能富集分析，筛选出影响睾丸生殖功能的关键基因；对其进行CDS区克隆测序，并运用qRT-PCR、Western Blot、免疫组化等方法进行进一步研究。结果显示，差异表达基因共353个，其中114个基因在小睾丸中表达上调，239个基因表达下调。GO功能注释表明，差异基因参与代谢、发育、生殖等过程；KEGG富集分析显示，差异基因主要富集在Wnt信号传导、PI3K-Akt信号通路等与睾丸生殖功能相关的途径。最终筛选出睾丸生殖功能基因CCN4进行进一步研究。qRT-PCR和Western Blot结果显示，CCN4基因和蛋白在大睾丸中表达量显著低于小睾丸。免疫组化结果显示，CCN4蛋白主要在精细胞外的区域(特别是支持细胞)中表达，且大睾丸中的相对表达量显著低于小睾丸，推测其高表达会抑制支持细胞发育，从而影响睾丸生殖功能。综上，研究结果可为深入开展牛睾丸大小对生殖功能影响的分子机制研究提供基础资料。

【目的】通过对大豆受豆卷叶螟幼虫胁迫下的转录组和蛋白质组结果进行联合分析,筛选出一些与大豆抗豆卷叶螟相关的候选基因,为深入认识大豆抗豆卷叶螟的分子调控机制奠定基础。【方法】以高抗材料赶泰-2-2(HR)和高感材料皖82-178(HS)为研究对象,运用RNA-Seq技术和iTRAQ技术鉴定出豆卷叶螟幼虫取食诱导0和48 h时样品间的差异表达基因(DEGs)和差异表达蛋白(DEPs),将在蛋白水平和转录水平关联到的所有可定量数据进行关联分析,计算蛋白水平和转录水平间的相关系数。【结果】蛋白质鉴定结果表明,HR

利用差示扫描量热法 (DSC) ,对豆浆中的 2种主要大豆蛋白质β伴大豆球蛋白 (7S球蛋白 )和大豆豆球蛋白 (11S球蛋白 )的变性温度进行了试验研究。结果表明 :豆浆中 7S球蛋白的变性温度为 (70± 2 )℃ ,11S球蛋白的变性温度为 (90± 2 )℃ ;豆浆中固形物的质量分数 (5 .76～ 10 .37 )与 7S球蛋白和 11S球蛋白的变性温度呈显著的正相关关系。当加热到 70℃时 ,7S球蛋白变性 ;当加热到 95℃时 ,7S球蛋白和 11S球蛋白都发生了完全变性

【目的】皂苷是苦瓜果实苦味的主要来源，具有降血糖、抗癌等多种药用价值。鉴定和挖掘参与调控苦瓜皂苷生物合成的代谢通路及基因，为进一步深入解析苦瓜皂苷形成的分子机制提供参考。【方法】以苦瓜高代自交系GK24为试验材料，分别在子房期（T1）、幼果期（T2）、商品果期（T3）和完全成熟期（T4）取果实组织，通过香草醛-冰醋酸法测定不同时期的苦瓜皂苷含量，利用转录组测序的方法鉴定差异表达基因。【结果】从T1到T4，苦瓜内源皂苷含量随着果实发育而呈现显著下降的趋势。转录组测序共鉴定到17 504个基因，在T1-vs-T2、T2-vs-T3和T3-vs-T4三组比较中，下调表达基因数量均高于上调表达基因数量。GO富集分析表明含磷复合物代谢过程、驱动蛋白复合体和2-琥珀酰-6-羟基-环己二烯-1-羧酸合成酶活性分别是生物过程、细胞组分及分子功能模块最显著富集的条目。KEGG分析显示全局与概述图谱、碳水化合物代谢和氨基酸代谢是3组比较中共同的最显著富集的代谢通路。根据表达模式可将基因划分为20个模块，其中profile 0中基因的表达量从T1到T4逐渐降低，与苦瓜果实发育过程中皂苷含量变化规律一致。从profile 0中鉴定出与苦瓜皂苷生物合成相关的基因14个，包括萜类骨架生物合成通路上的基因AAT1、HMG1、MVK、PMK、MVD2和FPS1，倍半萜与三萜生物合成途径中的基因SS12、SQE1和CPQ及后修饰酶基因CYP97A3、CYP71AN24、UGT94E5及UGT73C6。实时荧光定量PCR(qRT-PCR)结果与RNA-seq结果一致。【结论】苦瓜果实中皂苷含量随着果实不断成熟而逐渐降低。苦瓜皂苷生物合成首先通过上游的萜类骨架生物合成（ko00900）代谢途径完成萜类皂苷骨架的积累，随后通过倍半萜与三萜生物合成（ko00909）代谢途径以及下游的氧化还原酶、糖基转移酶的修饰最终生成皂苷，本研究通过转录组测序共鉴定到14个与苦瓜皂苷生物合成相关的基因。

【目的】研究巨型稻幼苗期响应盐胁迫差异基因的表达特异性，解析巨型稻耐盐机理及分子机制，为水稻耐盐性遗传改良提供理论依据。【方法】以湘巨1号为材料，设置对照和1个盐浓度水平(0.5%NaCl)，处理72 h，随后应用转录组测序技术分析盐胁迫下响应基因的动态变化与主要富集途径。【结果】盐处理后湘巨1号幼苗的苗高和鲜重均受到低于30%的胁迫损害。每个样本共获得的干净片段为58.93 G Clean bases，平均GC含量为53.51%，平均Q20和Q30分别为97.84%和93.89%。基因差异表达分析筛选出2327个差异基因，其中1363个上调表达基因，964个下调表达基因。对差异基因进行功能注释分析发现，GO富集分析结果显示差异基因显著富集在166个GO条目中，其中最显著富集的条目分别是生物过程中的氧化还原过程、核代谢过程、转录调控等；细胞组分中的膜组分；分子功能中的氧化还原酶活性、单氧酶活性、转录因子活性等。KEGG分析结果表明，差异基因显著富集在次生代谢物合成、植物激素信号转导、苯丙烷生物合成、谷胱甘肽代谢等途径。【结论】湘巨1号具有较强耐盐性，GO分析表明湘巨1号通过调控抗氧化酶活性和转录因子表达来响应盐胁迫。KEGG分析表明植物激素信号转导、苯丙烷生物合成、谷胱甘肽代谢等途径在响应盐胁迫中起主导作用。本研究利用RNA-Seq测序技术筛选了水稻响应盐胁迫的候选基因，为深入研究巨型稻耐盐调控机制提供了基础，也为后续新品种研发培育提供了重要的理论依据。

【背景】绿壳蛋深受消费者喜爱，蛋壳绿色的深浅在市场上是影响绿壳蛋定价销售的重要参考指标。绿壳蛋的形成受多基因共同调控，蛋壳绿色深浅不一，其分子机理尚不清楚。通过对赤水乌骨鸡蛋壳腺组织进行转录组测序，挖掘其调控绿壳蛋蛋壳颜色深浅的候选基因以及关键信号通路，探究影响蛋壳颜色遗传性，以期发展绿壳蛋的选种选育和提高经济效益。【目的】探讨赤水乌骨鸡的遗传基础，并通过SLCO1B3基因分型对其进行鉴定和筛选，以期通过分子标记在青壳鸡育种规划中提供新的见解，并在后期的选择策略中帮助控制和提高赤水乌骨鸡蛋壳品质的同质性。【方法】以纯系的280日龄赤水乌骨鸡为研究对象，屠宰产浅绿色蛋（QL）和深绿色蛋（SL）的母鸡各3只，采集蛋壳腺以RNA-seq技术检测分析，筛选与蛋壳颜色密切相关的差异表达基因（differentially expressed genes, DEGs），并对这些DEGs进行GO和KEGG富集分析。利用qRT-PCR技术检测与蛋壳颜色相关的6个候选基因的转录水平变化以验证转录组数据可靠性。【结果】在深绿组与浅绿组中共筛选到93个显著DEGs，有59个基因在深绿组中显著上调，34个显著下调。DEGs注释到GO数据库进行比对，主要显著性富集的是钠离子转运、负离子结合、肌质网膜等。KEGG分析富集程度最显著的是醛固酮调节钠的重吸收，以及富集的矿物质吸收、亚油酸的新陈代谢等信号通路。通过qRT-PCR验证，结果显示这些基因的表达趋势与转录组的测序结果一致。【结论】经功能分析，TF基因、SCNN1家族基因、CYP450家族基因、SLC家族基因和FAM家族基因，以及醛固酮调节钠的重吸收信号通路参与蛋壳色素合成、转运和沉积。上述基因和信号通路可能是影响赤水乌骨鸡蛋壳绿色深浅不一的候选基因和关键信号通路。

【目的】定位大豆蛋白质和油分含量QTL及互作分析,为大豆品质性状QTL精细定位和分子辅助育种提供基础。【方法】以Charleston和东农594为亲本,构建了含147个株系的重组自交系,以F2:19—F2:20代重组自交系为试验材料,利用Windows QTL Cartographer V.2.5软件的复合区间作图法和多重区间作图法,对该群体的蛋白质和油分含量进行QTL定位分析,并利用QTL Network 2.1软件分析QTL间的上位性效应及环境互作效应。【结果】采用CIM和MIM 2种算法在2011和2012年哈尔滨、红兴隆、佳木斯和牡丹江每年3个地点共6个种植环境下共定位了9个蛋白质和1...

采用ＩｌｌｕｍＩｎａ高通量测序平台对黑壳与金壳葡萄牙牡蛎（Ｃｒａｓｓｏｓｔｒｅａ ａｎｇｕｌａｔａ）外套膜组织开展转录组测序研究，其中黑壳组总计产出５ ８５７ ７１８ ０９２ ｎｔ数据，筛选后的干净阅读子５８ ００４ ５６０个；金壳组产出６ ７２７ ４４７ ３２６ ｎｔ数据，干净阅读子６６ ６１４ ８５４个。将发掘的新基因与Ｃｏｇ、ｇｏ、Ｋｅｇｇ、ｓｗＩｓｓ－Ｐｒｏｔ、ｎｒ数据库进行序列比对，最终得到各数据库注释的新基因数量分别为１４７、３９７、２０３、５６７、１ ３４７个。通过对差异表达基因的序列进行生物学功能注释和聚类分析以及利用已知贝壳基质蛋白序列进行搜索，鉴定出１４条可能与壳色表达有关．．．