基于转录组测序的槜李EST-SSR标记开发

为开发槜李EST-SSR标记,本研究利用MISA软件筛选了槜李花转录组测序获得的35584条Unigenes,对其SSR信息进行分析后,利用Primer Premier 3.0软件设计EST-SSR引物,并随机选取40对SSR引物对12个李品种进行EST-SSR引物筛选及多态性分析。结果发现,在槜李花转录组中共搜索到个10791个SSR位点,分布于8433条Unigenes,SSR发生频率为23.70%,平均每3.71 kb含有1个SSR;SSR重复基元中二核苷酸重复出现频率最高,占总SSR数量的52.98%,其次为三核苷酸重复(占24.00%)和单核苷酸重复(占20.95%);二核苷酸重复基元以AG/CT为主(85.95%),三核苷酸重复基元以AAG/CTT为主(31.24%)。利用Primer Premier 3.0软件共设计出9870对候选引物,随机选择40对引物对12个李品种进行SSR引物筛选及多态性分析。40对引物均能扩增出预期大小的条带,有效扩增效率为100%,40对引物中有5对引物在12个李品种中表现出多态性。本研究开发的EST-SSR标记可为李属植物遗传多样性分析提供丰富的候选标记,同时可为槜李发育相关功能基因定位、遗传图谱构建、及分子标记辅助育种等研究提供帮助。     

云南金花茶转录组序列分析及功能注释

【目的】深入研究云南金花茶优良性状、基因功能和基因分布特征,可为云南金花茶乃至山茶属植物功能基因的探索以及该物种的遗传保护提供基础生物信息学理论参考。【方法】采用Illumina Hi Seq 2000技术,对云南金花茶叶片进行转录组测序及相关数据分析,在对测序数据整理、de novo组装后,获得Unigene序列,继而利用相关生物信息学数据库进行序列比对。【结果】组装后,共获得95 979条Unigenes,其中N50(拼接转录本不小于总长50%的长度)为1 660 nt,平均长度为1 124 nt,Q20和Q30序列(处理后质量高于20和30的碱基)分别占96.39%和91.28%。经比对,在Nr、Nt、Swiss-Prot中能得到注释的Unigene分别为58 830、43 623、44 315条。在GO中,有41 905条(43.66%)Unigenes注释224 129个GO功能,将其分为3个大类和56亚类,所占比例最多的为生物过程这一类功能。在KOG中,有23 499条(24.48%)Unigenes注释26 430个KOG功能信息,将其分为26个基因功能大类,其中表达量较高的分别是一般功能基因和翻译后修饰、蛋白质转化、伴侣功能的相关基因。另外,共有23 214条(24.18%)Unigenes在KEGG中得到注释,根据其涉及的相关通路可将其归为19个亚类,其中以代谢通路较为丰富。此外,利用相关数据库和ESTScan软件对云南金花茶转录组Unigene进行CDS比对和预测,共预测到26 428条CDS,其长度集中在100～500 nt,占总CDS的82.10%。【结论】云南金花茶所含基因信息丰富,利用分析得到的所有注释信息可以更深层次地探索其基因组信息和基因分布情况。     

不同水稻种质中渗透胁迫抗性基因DREB2A的遗传多样性分析

本研究旨在分析转录因子DREB2A基因在不同水稻种质中的遗传多样性,以期为水稻耐渗透胁迫遗传改良提供分子工具。利用单倍型分析、系统进化树、遗传距离和密码子偏好性分析,对85份不同类型水稻种质中DREB2A基因的功能性核苷酸序列变异及遗传多样性进行了研究。共鉴定出55个核苷酸变异位点,其中12个位于编码区,43个位于非编码区;鉴定出12个DREB2A等位基因型,其中来自非洲栽培稻(Oryza glaberrima)的3个等位基因型表现大片段变异;根据DREB2A基因的序列变异鉴定出39个单倍型,其中33个为新单倍型;将所有单倍型分为三组（Group I、II和III）,其中非洲栽培稻(Oryza glaberrima)中的10个单倍型单独分为一组(Group III);系统进化树、遗传距离和密码子偏好性分析均表明Group III与其他两组具有较大差异。对85份不同类型水稻种质中DREB2A基因的序列分析表明,非洲栽培稻(Oryza glaberrima)中的DREB2A等位基因在序列变异、系统进化关系和密码子偏好性方面均明显不同于其他种质材料中的等位基因。     

植物冷驯化转录调控的影响因素

植物在冷驯化过程中通过转录调控发生大量生理生化上的变化,遗传物质、光照、温度等是影响转录调控的关键因素。本文从遗传物质及环境因子方面对植物冷驯化转录调控的影响进行了概述,希望人们在设计实验时能对其重要性予以重视。     

花生AhbHLH63基因的克隆与表达分析

bHLH转录因子在调节植物生长发育中起着重要作用。本研究利用转录组测序结果从栽培种花生丰花1号中克隆得到AhbHLH63基因。序列分析显示,AhbHLH63基因有两个剪接体AhbHLH63-1.1和AhbHLH63-1.2,其ORF分别长1 188 bp和1 152 bp,分别编码395个氨基酸和383个氨基酸组成的蛋白。亲缘关系分析表明,AhbHLH63蛋白含有一个预测的bHLH结构域,与其他植物的bHLH63具有较高的同源性,与苜蓿等的同源蛋白亲缘关系较近。实时荧光定量PCR分析表明,AhbHLH63-1.1和AhbHLH63-1.2表达模式基本一致,均为组成型表达,在花中表达量最高,茎和叶中次之,在根中表达量最低;在花生种子的不同发育时期,发育中期表达量较高。     

野生Paxillus ammoniavirescens菌的分子鉴定及生物学特性研究

[目的]对野生食用菌进行分子鉴定及生物学特性研究,为其进一步液体深层发酵培养和开发利用奠定基础。[方法]采用rDNA ITS序列分析方法对一野生菌进行分子鉴定,通过单因素和正交试验确定了该菌株的生物学特性。[结果]由新鲜子实体分离得到的纯菌株为Paxillus ammoniavirescens,其菌丝体生长的最适碳源和氮源分别为可溶性淀粉和蛋白胨,最适培养温度为25℃,最适pH为4.5。[结论]卷边网褶菌菌丝体液体培养的最佳配方为:可溶性淀粉20 g,蛋白胨2 g,KH_2PO_4 3.0 g,MgSO_4 4.5 g,水1 L。按此配方培养,可获得平均干质量为80 mg的菌丝体,明显高于其他配方。     

水杉长链非编码RNA分析

水杉(Metasequoia glyptostroboides)是中国著名的一级保护植物,素有"活化石"之称。为了更好地保护和利用这一重要种质资源,本研究以水杉为实验材料,在全长转录组测序的基础上,利用四个软件对61057个Unigenes进行了蛋白编码潜能预测,最终获得3385个长链非编码RNA(long non-coding RNAs,lncRNAs)分子。进一步对这些lncRNAs的长度和数量作统计分析,结果显示,最短为200 nt,最长为11134 nt,平均长度为1956 nt;大部分lncRNAs的长度介于200~2000 nt,占总量的59.5%;在2001~4000 nt的长度范围内,也发现较多的lncRNAs,达到总量的32.2%;在长度为4001~6000 nt的范围内,只有7.6%;长度超过6000 nt的lncRNAs很少,总共只有25个(占总量的0.7%)。本研究结果为水杉资源的保护和利用研究提供了一些分子依据,为其它植物在相关领域的研究提供了一些参考。     

杜仲雌花芽2个发育时期转录组分析

杜仲(Eucommia ulmoides)的果实是提取杜仲胶和杜仲籽油的重要原材料,在工业、医疗、食品等领域开发前景广阔。为了揭示雌蕊原基发育相关基因的表达情况,本研究以杜仲果用良种‘华仲6号’(Huazhong No.6)为材料,利用lllumina Hiseq X-10高通量测序平台,分别对花序原基分化期和雌蕊原基分化期的花芽进行转录组测序,通过生物信息学对2个发育时期的雌花芽转录组进行比较分析,筛选出与杜仲雌蕊原基发育相关的差异基因。结果显示,转录组测序共获得68.11 Gb数据,各样品的Clean reads与杜仲基因组进行序列比对,比对效率为91.59%~92.05%,同时2个发育时期共筛选出485个差异表达基因,其中219个差异基因在雌蕊原基分化期表达上调,266个表达下调。差异基因GO富集结果表明,花的发育、光周期现象、激素生物合成以及蛋白结合转录因子活性等相关的生物途径被富集。KEGG富集结果表明淀粉和蔗糖代谢、碳代谢、昼夜节律和植物激素信号转导等代谢通路被富集。结果表明光周期途径是诱导杜仲成花的重要途径,且雌蕊原基发育受碳水化合物、植物激素和其他代谢物质调控。本研究为杜仲花器官发育调控基因的挖掘提供了基础数据,也为果用杜仲的分子育种提供了参考。     

团头鲂tf和tfr1a基因启动子克隆及转录调控分析

为探索鱼类转铁蛋白基因tf和转铁蛋白受体基因tfr1a的转录调控机制,本实验以团头鲂为研究对象,在其全基因组数据库中获取tf和tfr1a基因序列,对2个基因候选启动子区转录因子结合位点及CpG岛进行预测,通过PCR方法克隆得到tf和tfr1a基因近端启动子区不同长度片段,连接至pGL3-Basic/pEGFP-1载体,瞬时转染入Hela细胞,并采用双荧光素酶报告基因检测系统进行检测。结果发现,团头鲂tf基因启动子区无CpG岛位点,而tfr1a基因启动子区有2个CpG岛位点。成功构建9个tf和10个tfr1a不同长度启动子片段的重组质粒,经双荧光素酶报告基因系统检测发现,tf启动子核心区域为-268^+56 bp,且-1 308^-1 102 bp片段可能存在正调控该基因表达的转录因子结合位点;tfr1a启动子核心区域为-224^+48 bp,且+48^+92 bp可能存在抑制该基因转录的负调控元件,而-1 229^-1 219 bp区域可能存在促进tfr1a基因表达的正调控转录因子结合位点。