随着高通量测序技术的出现,测序数据量呈现指数级增长,各类公共数据规模日益庞大,公共数据整合分析已逐渐成为很多研究中的一项重要分析手段。利用公共数据一方面可以降低研究成本,另一方面可以补充一些受限于研究者研究背景与技术水平而难以自主检测的关键数据。截至目前,高通量测序数据库SRA数据库已收录超过10Pbase的NGS数据,但是由于公共数据整合分析过程中下载、存储和分析各个环节都存在较高的技术门槛,导致研究者对公共数据的利用不充分。对于这些问题百迈客已经提供了非常成熟的解决方案,一键式导入即可完成数据的下载和存储,辅以平台上多达25个高度集成化的分析流程,可视化界面助您轻松搞定分析难题。中国农业科学院马有志老师课题组的徐兆师研究员等研究人员在百迈客云平台结合公共数据完成拟南芥抗病性研究,成果于2018年1月22日发表在《Frontiers in Plant Science》上,影响因子4.298。 摘要 Bax抑制剂-1(BI-1)是真核生物中进化保守的内质网(ER),定位细胞死亡抑制因子。 BI-1抑制生物和非生物胁迫的能力在拟南芥中得到了充分研究,而小麦中BI-1的功能在很大程度上是未知的。在这项研究中,将禾谷镰刀菌(Fg)处理的小麦进行RNA-seq分析,然后分离小麦BI-1基因TaBI-1.1。通过水杨酸(SA)处理诱导TaBI-1.1表达,并通过脱落酸(ABA)处理诱导TaBI-1.1的下调。基于β-葡糖醛酸糖苷酶(GUS)染色,TaBI-1.1在成熟叶和根中表达,但不在下胚轴或幼叶中表达。 TaBI-1.1在拟南芥中的组成型表达增强了其对丁香假单胞菌病毒(Pst)DC3000感染的抗性并诱导SA相关基因表达。此外,TaBI-1.1转基因拟南芥显示出由高浓度SA引起的损害的减轻,并降低了对ABA的敏感性。与表型一致,35S :: TaBI-1.1和Col-0植物的RNA-seq分析显示TaBI-1.1参与生物胁迫。这些结果表明,TaBI-1.1正向调节SA信号并在对生物胁迫的响应中起重要作用。此外,TaBI-1.1与水通道蛋白TaPIP1相互作用,并且它们都定位于ER膜(内质网膜)。此外,研究人员证明了SA处理后TaPIP1上调,并且TaPIP1转基因拟南芥增强了对Pst DC3000感染的抗性。由此表明,在ER膜上TaBI-1.1和TaPIP1之间的相互作用可能发生在对SA信号和防御反应的响应中。 使用拟南芥Columbia-0(Col-0)作为TaBI-1.1过表达的背景, 突变体atbi1-2从拟南芥生物资源中心(ABRC)获得。 从栽培的小麦品种小白麦扩增TaBI-1.1和TaPIP1的cDNA, 将PCR产物克隆到pLB载体中,使用DNAMAN 6.0软件进行氨基酸序列同一性比较。 RNA-Seq Analysis 收集来自4周龄Col-0和转基因系35S :: TaBI-1.1的叶片用于RNA-seq分析,使用HiSeq 2500平台测序。测序数据与TAIR 10拟南芥参考基因组比对。使用TopHat和Cufflink软件包进行mRNAseq数据分析以及DEG的鉴定和分析。通过GOseq软件进行差异表达基因的GO富集分析。根据KEGG数据库通路进行KEGG富集分析,寻找差异表达基因的富集通路。在百迈客云平台(BMKCloud)完成小麦Fg处理的RNA-seq数据下载和分析, SRA数据库(登录号:PRJNA289545)。 其他实验:1.qRT-PCR分析;2.转基因拟南芥植株在胁迫处理下的应答及性状评价;3.β-葡糖醛酸糖苷酶(GUS)活性测定;4.细菌接种和细菌生长的监测;5.酵母双杂交系统;6.Pull-Down assay(免疫沉淀法);7.亚细胞定位测定;8.ELISA Assay(酶联免疫吸附试验) 分析结果 1.通过qRT-PCR和β-葡糖醛酸糖苷酶(GUS)染色确定TaBI-1.1的表达模式 研究人员分析了来自Fg处理的小麦RNA-seq数据,以研究植物防御信号通路。从RNA-seq分析中分离出前10个差异表达基因。在这10个基因中鉴定了两个小麦BI-1基因:TRIAE_CS42_U_TGACv1_644608_AA2140670和TRIAE_CS42_6BS_TGACv1_515717_AA1671980。 TRIAE_CS42_6BS_TGACv1_515717_AA1671980在之前的研究中被命名为TaBI-1。在小麦Ensembl数据库中仅筛选出三种BI-1基因:TaBI-1,TRIAE_CS42_U_TGACv1_644608_AA2140670和TRIAE_CS42_6AS_TGACv1_488014_AA1573990。在两个差异表达的BI-1基因中, TaBI-1.1的差异最大。根据氨基酸比对,TaBI-1.1与TaBI-1同一性达到99.19%。与对照相比,Fg处理对TaBI-1.1的表达水平上调24倍。许多关于AtBI-1的研究已经证实AtBI-1在植物中对生物和非生物胁迫的抗性中起关键作用。 TaBI-1.1蛋白的序列与AtBI-1有69.88%的同一性,表明该序列在BI-1家族中基本保守。使用qRT-PCR监测TaBI-1.1的表达模式以确定TaBI-1.1在生物和非生物胁迫中可能的作用。TaBI-1.1的表达在响应SA处理时显着上调,在响应ABA处理时下调。 SA处理48小时后TaBI-1.1表达达到8倍高峰(图1A)。 响应ABA处理的下调水平在8小时达到初始水平的约1/5(图1C)。 响应NaCl处理,表达水平在4小时后下降,在2小时稍微增加并且在8小时后回到其初始水平(图1B)。 研究者创建了含有1.7kb TaBI-1.1启动子并生成转基因拟南芥的PBI :: GUS融合构建体,以研究TaBI-1.1的空间表达模式。使用GUS染色测定组织GUS活性。在成熟叶和根中有TaBI-1.1表达,但在下胚轴和幼叶中观察不到(图1D)。还检测了各种小麦组织中的TaBI-1.1表达水平,包括根,茎,叶和小穗,结果表明TaBI-1.1的表达在这些小麦组织中普遍存在,而且在叶中最高,在小穗中最低(图1I)。在SA和NaCl处理之后,成熟叶中的表达与对照相比增加,并且SA处理的表达更高。当SA和NaCl处理时,在下胚轴中也检测到表达(图1E,F)。在ABA处理的植物的叶子中检测到较弱的表达(图1G)。通过qRT-PCR进一步证实GUS表达水平(图1H)。以上结果表明,TaBI-1.1在各种胁迫下均有表达。 图1. 通过qRT-PCR和GUS染色评估TaBI-1.1的相应表达模式 2.TaBI-1.1的表达增强了拟南芥对Pst DC3000感染的抗性 经过Fg和SA处理后表达上调,假设TaBI-1.1可能参与生物胁迫反应。研究人员在拟南芥中异位表达TaBI-1.1,在花椰菜花叶病毒(CaMV)35S启动子的控制下验证这一假设。选择两个TaBI-1.1表达水平相对较高的纯合株系3和7(35S :: TaBI-1.1#1和35S :: TaBI-1.1#2)用于进一步分析(图2A)。将atbi1-2,Col-0和两个转基因品系的四周龄叶子进行Pst DC3000感染或10mM MgCl 2处理(模拟物)。在模拟实验中,在用10mM MgCl 2处理3天后,在atbi1-2,Col-0和两个转基因品系的叶中没有观察到明显差异(图2B)。在用600nm(OD600)0.002的光密度接种Pst DC3000,3天后在这些叶中检测到疾病症状。 atbi1-2突变体表现出的症状较严重,因为几乎所有的叶片都表现出严重的萎黄和坏死,而两种转基因品系表现出比atbi1-2和Col-0更轻微的疾病症状。转基因植物的一小部分叶子是绿色的,没有表现出萎黄和坏死。 Col-0叶中疾病症状的程度介于atbi1-2和两个转基因品系之间(图2C)。在接种Pst DC3000或10mM...

百迈客微生物宏基因组云平台正式上线啦!专业的研发团队精心打造出操作简单、运行快速、分析丰富的微生物宏基因组云平台。充分满足您的个性化需求,让您秒变生信高级分析工程师! 微云特点: 省时省力:极短分析时间,节省中间沟通环节,余留更多的精力用于数据分析。 操作简单:分析平台以简单式操作为原则,保证快速上手进行个性化信息挖掘。 专业可靠:专业的百迈客云+微生物信息开发团队精心研发,做到分析更专业。 分析全面:20+基本分析内容,全面上线,多角度多参数分析。   微生物宏基因组云平台功能简介: 基本分析支持用户自定义参数进行分析。可自定义宏基因组基本分析的主要参数,包括“综合选项”、“原始数据导入”、“数据质控过滤”,“基本参数设置”、“环境因子设置”、“样品分组信息”等六个参数模块 ,填写并确认参数信息后点击“提交”即可运行该项目基本分析,可在“总览/我的项目”或 “总览/我的任务”中查看。 微云平台界面展示: 结题报告图例模板: [gallery size="large" ids="3610,3611,3612,3613,3614,3615,3616,3617,3618,3619,3620"] 欢迎点击下方按钮,预约申请体验!       加入百迈客微生物云交流平台(QQ群:439540040)参加免费培训,12h在线答疑。    ...

棉花是世界上最重要的经济作物之一,在2018年5月8日中国农业科学院棉花研究所所长李付广研究员、武汉大学朱玉贤院士、中国农业科学院棉花研究所杜雄明研究员、中国农业科学院农业基因组研究所所长黄三文研究员、林涛博士与北京百迈客生物科技有限公司关于亚洲棉的合作成果发表在Nature Genetics上,论文题目为“Resequencing of 243 diploid cotton accessions based on an updated A genome identifies the genetic basis of key agronomic traits”。该研究以全新亚洲棉基因组为基础,增加遗传进化和GWAS研究,对棉的品质、产量、抗病等重要农艺性状进行研究。其中朱玉贤院士与李付广研究员为通讯作者,杜雄明研究员为第一作者。 以下为文献详细解读: 英文题目:Sequencing of 243 diploid cotton accessions based on an updated A genome identifies the genetic basis of key agronomic traits. 中文题目:以更新的亚洲棉A基因组为基础的243份二倍体棉花的重要农艺性状的研究 1. 摘要: 亚洲棉(Gossypium arboreum)和草棉(Gossypium herbaceum)的祖先是现代栽培异源四倍体棉花A亚基因组的供体。本研究中通过整合了不同的技术,提升了亚洲棉的基因组组装水平;同时对243株二倍体棉花(亚洲棉和草棉)进行全基因组重测序分析,绘制基因组变异图谱,并发现亚洲棉和草棉(A)与雷蒙德氏棉(D)同时进行了分化;单独的对亚洲棉分析表明亚洲棉起源于中国南部,随后被引入长江和黄河地区,大多数具有驯化相关特性的种质都经历了地理隔离;通过亚洲棉的全基因组关联分析(GWAS),鉴定了亚洲棉11个重要农艺性状的98个显著关联位点,GaKASIII的非同义替换(半胱氨酸/精氨酸替换)使得棉籽中的脂肪酸组成(C16:0和C16:1)发生了变化;棉花枯萎病抗性与GaGSTF9基因的表达激活相关。本研究对理解棉花A亚基因组的进化具有重要的意义。 2. 研究背景: 棉花是世界上最重要的商业作物之一,同时也是研究植物多倍化的有价值的资源。亚洲棉最可能在马达加斯加或印度河流域文明(巴基斯坦摩亨佐达罗)开始驯化,随后分散到非洲和亚洲一些地区。亚洲棉最初在1000多年前作为观赏植物引入中国。当地方的农业生态环境的适应和人类选择影响的过程中,中国的Gossypium arboreum形成了独特的地理种群,称之为“sinense...

2017年12月,吉林大学袁宝老师在Scientific Reports发表了一篇关于大鼠垂体的研究,鉴定了不同发育阶段的垂体中的lncRNAs,并构建了促卵泡生长激素FSH的调控网络。 Identification of long non-coding RNAs in the immature and mature rat anterior pituitary 成熟与未成熟的大鼠垂体前叶中的lncRNA鉴定 期刊: Scientific Reports    影响因子:4.259 1.研究背景 垂体是哺乳动物最重要的内分泌腺,分前叶和后叶两部分。它分泌多种激素,如生长激素、促甲状腺激素、促肾上腺皮质激素、促卵泡生长激素(FSH)等。这些激素对代谢、生长、发育和生殖等有重要作用。在未成熟的垂体中,激素的含量应该在适当的范围内,以适应身体的发育和生长。在成熟的脑垂体中,激素主要影响生殖能力。在大鼠的垂体研究中,更多的是关注基因和miRNA,对lncRNAs的研究比较少。最近的研究报道lncRNAs参与了垂体的肿瘤发生,说明lncRNA很可能通过调控垂体相关生物学过程来调节生长发育与生殖。 2.材料方法 样品:15日龄和120日龄大鼠(Rattus norvegicus)各取3只,同日龄的3只大鼠的垂体前叶混合成1个样品,2个样品测序数据量分别是14G和15G。 平台:百迈客Illumina HiSeq X ten平台 分析平台:百迈客云平台(BMKCloud),利用云平台做了个性化分析 分析内容: 图1  分析内容 3.研究结果 1.鉴定未成熟和成熟垂体前叶形态学差异 我们检测了两个发育阶段垂体前叶形态学上的差异,成熟的垂体前叶有更清晰的边界和更多的毛细血管,因此我们推测有某些因子导致了大鼠垂体的形态和功能差异,根据近几年的研究,这种因子可能是lncRNA。 图2 未成熟和成熟的垂体组织的形态学鉴定 (A.15日龄垂体前叶放大100倍;B. 15日龄垂体前叶放大500倍;C. 120日龄垂体前叶放大100倍;D. 120日龄垂体前叶放大500倍) 2.lncRNA的鉴定和差异分析 对两个日龄的垂体前叶分别测序,比对效率分别为95.06% 和95.15% ,本研究中我们鉴定共到7039个lncRNA,lncRNA的平均长度为2052bp,比mRNA略长。而76.5%的lncRNA含有2个外显子,mRNA平均含有4.1个外显子。 图3  lncRNA和mRNA的长度(A)、外显子个数(B)的比较 使用fold change≥2且FDR< 0.05,筛选到1181个差异表达的lncRNA,其中有607个是上调的。基于位置关系进行靶基因预测,发现靶基因主要参与了抗原处理和呈递途径,自身免疫甲状腺疾病途径和I型糖尿病途径。 图4  lncRNAs 和mRNA的表达量分布情况 促卵泡生长激素(FSH) 是垂体前叶分泌的一种糖蛋白激素,其主要功能是卵巢的卵泡发育和成熟。我们在反式靶基因中找到了调控FSHb基因的lncRNA,并用其中3个新鉴定到的lncRNA构建调控网络图。 图5 qPCR验证10个高表达的lncRNAs和3种与FSHb基因相互作用的新lncRNA (A. 10个高表达的lncRNA的定量结果;B. lncRNAs和FSHb基因的相互作用网络;C-E.三种新lncRNA的定量结果) 3.qPCR验证 我们挑选调控FSHb的3个新lncRNA 和十倍差异的lncRNA 做qPCR验证,定量结果和测序结果吻合。前人的研究表明FSHb基因的表达量在性成熟阶段明显增加,这三个新lncRNA的表达趋势和FSHb mRNA一致。 结论 本项目系统地研究了lncRNA在15日龄未成熟小鼠和120日龄成熟大鼠垂体前叶的情况,共鉴定到7039个lncRNA,1181个差异表达的lncRNA,并预测靶基因。顺式靶基因主要富集到抗原处理和呈递途径,自身免疫甲状腺疾病途径和I型糖尿病途径;根据反式靶基因构建了调控网络,发现有3个新lncRNA能够靶向调控促卵泡生长激素FSHb基因,为后续的深入研究奠定了基础。 研究亮点 垂体分泌的几种激素,在发育阶段起着关键作用,尤其是促卵泡生长激素FSH,能够促进卵泡的发育和成熟。lncRNA在生物学过程中有重要作用,如表观修饰,细胞分化,转录调控,生长发育,抗病等。最近的研究表明,lncRNAs参与垂体肿瘤的发生,这项研究鉴定了在大鼠垂体前叶各发育阶段的lncRNA的情况,并找到了调控FSHb的lncRNAs。 参考文献 BaoY, Zhang J B et...

今天小编带大家看一看百迈客云物种数据库的核心功能模块。 百迈客云物种数据库采用模块化设计,标准功能模块主要有三大类   1.存储&查询模块         存储&查询模块设置了种质资源、突变体信息、性状座位信息、分子标记、基因功能信息、基因组数据、转录组数据的存储和查询模块 01物种种质资源存储&查询模块 种质资源是一种遗传资源。种质系指农作物亲代传递给子代的遗传物质,它往往存在于特定品种之中。如古老的地方品种、新培育的推广品种、重要的遗传材料以及野生近缘植物,都属于种质资源的范围。用户可在该模块中按照物种的名称、编号、父母本、省份、地理位置、表型查询相关材料品种概况。   02物种突变体信息存储&查询模块 用户在该模块通过突变体名称或表型检索相关突变体,搜索结果可展示出突变体名称、编号、突变基因信息等。其中,点击突变基因可链接至“基因查询模块”查询基因功能、序列、基因组位置等详细信息。   03物种性状座位信息存储&查询模块 在该模块用户通过输入表型数据查询该性状控制座位在基因组上的定位结果,该结果整合自不同研究不同定位方法定位结果,并可视化展示相应性状座位对于于最新版本基因组上的位置及相关基因。其中,相关基因也可链接至“基因查询模块”查询基因功能、序列、基因组位置等详细信息。   04分子标记存储&查询模块 用户可按照精确编号、模糊编号(正则匹配)、或连锁性状搜索分子标记,搜索结果含有标记序列信息、标记扩增引物序列、标记上下游基因(可链接至“基因查询模块”)、标记相关性状等。   05基因功能信息存储&查询模块 在该模块中,用户可以基因名称、基因ID、功能关键词、物种基因组起始位置、转录本序列、是否定位、是否克隆为条件检索目标基因,检索结果包括基因转录本序列(核苷酸or氨基酸序列)、染色体、起止位置、各功能数据库注释结果、相关文献、相关性状(可链接至“物种性状座位信息存储&查询模块”)。   06基因组数据存储&查询模块 1)使用基因组版本号检索相关基因组信息,检索结果包括该基因组序列组装、结构注释、功能注释、测序数据等,所有结果均支持本地下载; 2)可按照物种名检索相关泛基因组信息,检索结果包括各品种基因组序列组装、结构注释、功能注释、变异位点、测序数据等,所有结果均支持本地下载; 3)按品种名检索该品种全基因组重测序、简化基因组重测序测序数据及变异位点检测结果。   07物种转录组数据存储&查询模块 1)用户可通过基因ID与基因name检索相关基因下转录本序列、项目中的转录水平、差异表达、共表达、互作信息等,检索结果支持本地下载。 2)按项目名称检索相关转录组学项目测序数据及完整分析结果,检索结果支持本地下载。   2.分析工具模块 分析工具均可从网页调用,无需用户掌握任何编程基础。包含数据库自带分析工具,比如 blast、blat、引物设计、转录因子、metaQTL等;当然也可调用部署在百迈客公有云计算平台(www.biocloud.net)的分析流程(APP)与生信软件(Tools。APP包括有参/无参转录组分析APP、lncRNA分析APP、microRNA分析APP、circleRNA分析APP、重测序分析APP、GWAS个性化分析APP、BSA个性化分析APP、遗传图谱个性化分析APP、微生多样性分析APP等;Tools模块包含100余款各类生信分析软件,涵盖比对、注释、绘图、统计等;  BMKCloud个性化分析示例 3.公共数据库模块 包含7大模块: 01高通量测序数据库 • 该公共数据库模块同步NCBI维护的SRA、GEO两大高通量测序数据库,目前该模块12PBase,1100个属的高通量测序数据,数据类型涵盖RNA-seq、全基因组重测序、WES、基因组denovo组装、甲基化、ChIP-seq等。 • 用户可按物种名称、数据编号、研究方向检索相关高通量测序公共数据,检索结果无需下载,一键导入数据库自带分析工具或百迈客公有云计算平台(www.biocloud.net)分析工具即可进行数据分析。   02功能基因数据库 该公共数据库模块整合包括GO、KEGG、miRbase、lncRNAdb、MiRTarBase、UniProt、genbank等在内的30多个基因功能&序列类数据库,目前收录了约1700万条基因相关信息,用户输入基因名称、ID、功能关键词后,一键即可获取基因的功能、序列、同源基因、表达、相关ncRNA、变异、相关文献等信息。   03功能定位数据库 该数据库收录了所有公共的性状相关数据,目前已收录11个物种347个性状相关的QTL位点数据,如QTL、eQTL、候选基因、关联数据等信息,QTL性状位置可视化展示以及物种间的比较,为各个物种间定位信息的比较提供重要的数据。   04文献数据库 该公共数据库模块同步NCBI维护的PubMed文献索引数据库,目前收录2600万+文献信息,该模块支持按照发表日期、影响因子、文章类型、信息完整度(有无测序数据、有无全文等)进行精确检索,检索结果包括文献题目、摘要、全文链接、文献所涉及测序数据等。   05其他公共数据库模块 • 参考基因组数据库:包含有826个物种的924套参考基因组信息; • 变异数据库:整合了dbSNP库中的数据信息,目前收录了8亿+变异数据,支持通过rs、ss号进行信息的检索,并且支持筛选特定物种的、特定染色体、特定区域的变异信息; • 互作网络数据库:整合了PPI、StarBase、LncRNADisease三个数据库中的4类互作网络信息,目前共收录了306万互作网络关系,包含了lncRNA、miRNA、mRNA、circRNA四类数据。 ...

发表杂志:Plant Growth Regulation 影响因子:2.646 摘要 荔枝(Litchi chinensis Sonn.)是一种重要的水果作物,起源于中国,并在世界热带和亚热带地区商业化种植。荔枝红色果皮的颜色是荔枝果实市场接受度的重要品质,果皮的粉红色/红色是花青素生物合成和积累的结果,花青素-3-葡萄糖苷和花青素-3-芸香糖苷是荔枝红果皮中的主要花青素。荔枝花青素的数量和组成在不同的品种间差异很大,也很大程度上受各种环境因素的影响。因此,深入了解荔枝中花青素生物合成的调控具有重要的科学意义和经济意义。 许多研究表明,花青素生物合成很大程度上受植物激素的影响。外源脱落酸(ABA)处理促进荔枝花青素生物合成,而外源性N-(2-氯吡啶-4-基)-N'-苯基脲(CPPU)抑制荔枝花青素生物合成。但是,ABA或CPPU调控花青素生物合成的机制仍不清楚。为了进一步了解外源ABA和CPPU调控荔枝果实花青素生物合成的分子基础,华南农业大学园艺学院赵杰堂老师课题组对ABA或CPPU处理的荔枝果皮进行转录组测序,并进行了全面的分析。 材料和方法 选取生长十五年的L. chinensis cv,妃子笑。 三种处理方式:ABA处理,CPPU处理和对照处理(自来水),每组处理约15个果实簇。处理后分别取0,10和20天果皮圆盘,共取7组样本分别是:Cont-0,Cont-10,Cont-20,ABA-10,ABA-20,CPPU-10,CPPU-20,分别提取总RNA,用于RNA测序。 测序平台:北京百迈客生物科技有限公司Illumina Hiseq平台 150PE 分析平台:百迈客云平台(BMKCloud) 分析内容:1.原始数据质控;2.参考基因组比对(未发表);3.差异表达基因分析(韦恩图绘制,KEGG富集分析);4.基因功能注释(使用Blast2Go 基于Nr,Nt,COG,KO,GO数据库的最高相似性进行基因注释);5.统计分析 花青素和叶绿素含量的测定以及内源ABA分析 研究人员根据Wrolstad(1982)等人的描述进行相应的修改对果皮中总花青素含量进行定量,根据Arnon(1949)描述的方法计算总叶绿素含量,根据Jia(2011)等人所述,进行了一些修改,使用气相色谱 - 质谱(GC-MS)测定ABA含量,每个过程重复三次。 结果 1.经外源ABA和CPPU处理后,荔枝果皮的表型、色素含量和内源ABA水平的变化 处理10天后,ABA处理和对照处理的果实开始着色(图1a),尽管花青素含量没有显着差异(图1b)。然而,CPPU处理的果实只有一点点红色(图1a),并且花青素含量远远低于ABA处理和对照的果实(图1b)。处理后20天,ABA处理的果实变成均匀的红色,并且CPPU处理的果实仍然是绿色的,带有一点红色(图1a)。同时,对照果实呈现不均匀的红色(图1a),这是cv妃子笑的典型特征。如图1b所示,ABA处理20天后明显促进荔枝果皮花青素积累,而CPPU处理显着抑制荔枝果皮花青素积累。与花青素含量变化相反,果实成熟过程中叶绿素含量降低。显然,ABA处理加速叶绿素降解和CPPU处理延缓了这一过程(图1c)。ABA处理的果实与对照之间的内源ABA水平没有显着差异,而CPPU处理的果相比而言具有较低的ABA水平(图1d)。 图1.各处理组果皮颜色、花青素、叶绿素及ABA含量变化 2.RNA测序以及参考基因组比对 共构建7个文库,测序过滤后每个样本数据均不低于6.6Gb,平均GC含量为45.36%(表1)。与荔枝参考基因组比对(未发表),比对率是74.7%(表1)。 表1. RNA-Seq结果展示 3.外源ABA和CPPU处理后荔枝果皮基因表达谱的综合分析 三种处理后0天,10天和20天的荔枝果皮中发现了许多差异表达转录本(图2)。荔枝果皮色素沉着过程中,在对照组鉴定了2263个差异表达基因(DEGs)(图2a)。除此之外, 经过ABA(图2b)和CPPU(图2c)处理后的三个阶段分别鉴定了1986和2862个DEGs。 图2.三组处理中差异表达基因韦恩图 为了鉴定ABA和CPPU调控的花青素生物合成的DEGs,分别选择ABA / CPPU处理后10天和20天的DEGs与对照组进行比较。 与对照相比,在ABA和CPPU处理的果皮中分别有579和827个DEGs。 表达模式分析表明,在ABA处理10天和20天后,上调的转录本数量均比下调的转录本数量多。而 CPPU处理后20天上调的转录本数量比ABA处理的低很多。 4.外源ABA处理后的荔枝果皮表达分析 外源ABA处理组,对显著差异表达的基因进行GO分析,发现了与细胞过程,代谢过程,细胞,细胞部分,催化活性和转运蛋白活性有关的GO terms的富集。也做了DEGs的KEGG通路富集分析。 579个DEGs中的351个被富集到101个通路中,前5个通路组分别为:植物 - 病原体相互作用通路(49DEGs),植物激素信号转导(32DEGs),类黄酮生物合成(27DEGs),内质网中的蛋白质加工(27DEGs) 碳代谢(20 DEGs),淀粉和蔗糖代谢(20 DEGs)。 外源ABA处理的DEGs中,大部分与花青素合成有关的DEGs上调,如PAL,C4H,CHS,CHI,DFR,LDOX和GTs。此外, FLS和LAR等黄酮醇和原花青素合成基因上调, C3'H,COMT,F5H,CCR,CAD等木质素合成基因也上调了。 ABA处理后另一类显着差异表达的unigene参与植物激素代谢和信号传导。两个编码ABA生物合成途径中的关键酶的NCEDs unigenes显著上调,而ABA信号转导的基因在外源ABA处理后没有明显的变化。外源ABA处理也影响了与生长素信号途径有关的基因表达。例如,属于Aux / IAAs家族的两个unigenes上调。另外,编码DELLA蛋白,负调控GA信号转导的unigenes经ABA处理后上调。 5.外源CPPU处理的荔枝果皮表达分析 将CPPU处理组与对照组之间的DEG进行GO分析。 与外源ABA处理相比,差异不显着。KEGG富集分析, 897个DEGs中的493DEGs个被富集到117个通路,前六个通路组是碳代谢(54DEGs),植物 - 病原体相互作用(49DEGs),光合作用(42DEGs),氨基酸生物合成(38DEGs),苯丙素生物合成 (33DEGs),植物激素信号转导(32DEGs)。 与ABA处理不同的是,外源CPPU处理后许多上调基因参与碳代谢,氨基酸生物合成和光合作用的基因,特别是在CPPU处理后10天较明显。例如,unigene Litchi__GLEAN_10019646在外源CPPU处理10天后显着上调6.1倍,它编码捕光复合体II叶绿素a / b结合蛋白1。另外一类在外源性CPPU处理中显著差异表达的unigenes涉及叶绿素生物合成代谢。 CPPU处理后,大部分参与类黄酮和花青素生物合成的DEGs均下调,如PAL,C4H,CHS和LDOX。此外,黄酮醇和原花青素合成基因,如FLS和LAR也下调。有29个DEGs比对到植物激素信号转导途径。其中,ABA信号中的两个ABA受体PYR / PYL均被CPPU下调。另外,与生长素,GA和乙烯信号有关的基因大部分也下调。 6.ABA和CPPU处理的荔枝果皮的表达分析 受外源ABA和CPPU共同响应的DEGs共199个,其中大部分unigenes在两种处理组中显示出类似的改变。 值得注意的是,有10个unigenes在转录水平上表现出相反的模式。 其中,编码GST4蛋白的unigene(Litchi_GLEAN_10051861)被报道与荔枝中花青素的液泡转移有关。 ABA处理和CPPU处理之间的显着差异之一是叶绿素代谢。 鉴定了与叶绿素生物合成和降解有关的24个候选基因,其响应外源ABA和CPPU的表达模式如图3所示。总体而言,ABA处理对叶绿素生物合成和降解的基因表达没有显著影响。 与对照相比,CPPU处理显着提高大多数叶绿素合成基因,并下调TF...

弹指一挥间,2017年转眼就到了要跟大家告别的时刻,在上一期的年终大盘点中跟大家分享了百迈客云平台用户的傲人成绩,在累累硕果的背后,一直有着一群兢兢业业的生信人,这一年里,他们是根据广大用户的体验反馈,不断改进,不断更新,使咱们的云平台来了一场华丽丽的大变身。细数这一年,云平台总共更新53次,更新项目达200多个,包括新主流程的上线,页面版块的调整,分析功能的完善,数据库的改进,平台速度的提升等。每一次改进都是为了使云平台的用户有一个更好的分析体验,每一次更新都希望用户可以有更高效的科研产出。今天,咱们一起来瞅瞅百迈客云平台17年到底有哪些大变样。   2017全新上线篇 “让基因分析更简单”是百迈客云平台前进的不竭动力,相信大家最关心的应该是云平台这一年有哪些新的分析流程上线。 2017 4/14  云平台上线lncRNA v3.0主流程; 4/21  全基因组关联(GWAS)分析平台(Beta版,云服务个性化)正式上线; 6/9   BSA 分析平台 v1.0(Beta版,云服务个性化)正式上线; 9/22  全基因组重测序分析平台3.0正式上线。做lncRNA,GWAS,BSA,全基因组重测序分析的用户有福啦! 另外,早在2017年3月3日,咱们云平台的云组第一期已正式上线,实现了组内成员计算、存储、工具的资源共享等功能,一个课题组的成员们可以同时享用同一个分析平台,再也不用为互传数据和结果而烦忧; 2017年6月30日 ,移动端微信公众号全新改版,在百迈客云微信公众号中添加文献检索、文献解读、文献收藏、云账号绑定、项目进展查看和评价等功能,同期,在线支付功能上线,套餐用户支持使用支付宝、微信在线购买分析平台、小工具,收货地址与个人信息关联,用户可以在线填写发票信息。   APP主流程更新篇 这么多新的APP成功上线了,显然是一件值得庆贺的事情,无独有偶,咱们已经在线上的APP也从来没有停止前进的步伐。 2017 1/20  无参转录组分析APP进行了全面更新还实现了图片和基因之间的交互,新增“WGCNA”个性化工具,新增基因结构挖掘功能等功能; 3/3  微生物多样性分析APP更新至2.0版,新增6项分析内容;并且美化了流程中各步骤生成的图片,确保所有图形结果有png图和矢量图等功能; 3/17  微生物多样性分析APP再次更新:增加5个个性化分析内容;同期,小RNA测序分析APP也进行了20多项更新:增加小RNA所对应的染色体、起始点、终止点信息和miRNA对应的前体位置信息,支持更加精细化的检索,添加一键绘图等功能; 3/26  微生物多样性分析APP更新了引导教程图片; 6/16  微生物多样性分析APP进行了功能更新:基本分析和个性化中RDA/CCA环境因子支持文件上传,增加保存历史参数功能等; 7/7  ChIP-Seq 分析APP的基本分析,增加了Input/IP设置,且支持有生物学重复的差异分析等功能; 7/21  真核生物有参/无参分析APP个性化分析进行了更新:基因挖掘中相应添加eggNOG分类图、KEGG分类图、KEGG富集图、GO富集图、MA图和火山图等; 7/28  真核生物有参/无参分析APP个性化添加差异分组时增加Pvalue选项; 8/25  lncRNA分析APP进行了20多项的更新:新增图片交互功能,增加火山图、MA图,增加矢量图, Readme.pdf文件,WGCNA分析等; 9/1  微生物多样性分析APP个性化更新,增加21个个性化且每个个性化都有多个参数可以调整分析; 12/1  微生物多样性APP第四期更新:对OTU个数分布图、物种分布柱状图、PCA分析等15项进行了优化。   极速体验更新篇 为了方便用户能够在短时间内体验到云平台的标准分析和个性化分析内容,研发人员特意准备了极速体验项目,标准分析在两分钟内得到结题报告,随后即可进入个性化分析流程。 2017 3/31  有参、无参、表达谱,小RNA,四个分析平台更新了极速体验项目,增加示例数据; 6/30  群体三个分析平台(BSA、全基因组关联、基于第二代测序标记的遗传图谱)增加极速体验功能。目前已有极速体验的分析平台包括:真核生物有参转录组分析平台、真核生物无参转录组分析平台、小RNA测序分析平台、LncRNA测序分析平台、微生物多样性分析平台、数字基因表达谱分析平台、医学有参转录组分析平台、外显子组测序分析平台。 是不是要为咱们的研发人员打个call? 云课堂更新篇 百迈客云课堂,绝对是您科研路上的好帮手。2017年百迈客云课堂先后共进行了5次更新,现在的云课堂集185+线上培训课程于一体,囊括RNA系列、DNA系列、生物信息软件使用和大家讲坛等课程,针对分析平台每款APP都有详细说明及步骤提示,无论有没有生物信息学基础,都能帮助用户在短期内玩转各类分析。   数据库更新篇 整合利用公共数据发表论文是近年来的研究热点,2017年百迈客云平台数据库先后更新9次,包含8大数据库、12PB公共数据、373万样本,共收录826个参考基因组、2000+种医学数据、306万互作网络关系、1700万基因、8亿变异数据。   小工具更新篇 2017年云平台小工具先后更新了4次,包含PCA分析、WGCNA分析、转录因子预测、遗传共线性分析、GO富集分析、聚类热图、柱状图、点线图、miRNA基因家族分类等104款小工具和3大工具集,覆盖12大功能类别:基因分析、遗传进化、ncRNA、质控、组装、比对、数据提取、突变、统计、表格处理、绘图、序列分析等。   文献功能更新篇 2017年云平台对文献模块先后进行了3次更新,如今查找文献和数据导入都非常方便,共整合了2700+万篇文献,同步NCBI数据,设置多项过滤栏,方便快速查找,还可轻松获得全文解读,文章数据可一键导入分析。   云平台其他功能更新篇 2017 3/10 云平台建立项目进展查询机制,用户通过云平台可以及时了解项目进度,增加项目线上评价功能,用户可以在线对项目各个维度的服务进行评价,我们可以通过评价提升项目服务质量,以及对页面内容进行整理更新,使项目的管理更清晰; 4/10  云平台结题报告进行全面优化,包括加载速度和流畅度,右侧导航可以智能展开和收起,图片可调大,添加专有名词解释; 5/19  Venn图支持多个区域选择,WGCNA添加分类模块功能富集结果图的展示,转录因子支持动物物种的预测,有参支持按照指定基因名称替换; 5/26  云科技官网上线文献检索、功能基因数据库和学术课堂功能,不用登陆即可体验超强功能; 7/14  真核有参、小RNA测序、微生物多样性、ChIP-Seq分析平台增加pdf格式结题报告下载功能; 7/21  云平台登录支持第三方微信便捷登录; 9/22  为了用户在体验过程中的流畅度,在平台登录页面添加了浏览器使用建议。贴心无处不在~   告别成绩斐然的2017,迎来了充满希望的2018,我们坚信新的一年百迈客云平台(BMKCloud)一定会再创辉煌,为广大用户提供一个更加完美的数据分析平台。作为一个生物大数据分析的开放云平台,拥有庞大的专业用户群和开发者用户群(15000+),为用户提供完整的生物信息分析以及整合利用公共数据的解决方案,我们一直在努力!...

百迈客云作为针对生命科学家的“Turn-key BioCloud”,2015年初出茅庐,2016年开始全国推广,一路走来得到了像中科院陈润生院士等生信领域权威专家的大力支持和帮助。百迈客云日臻完善,集成式的计算资源和工具、交互性的操作界面、开放性的公共数据库、实用的培训讲堂、专业的开发支持团队、一站式体验模式,成为广大科研工作者生信分析智囊团,极大地加速了科研创新孵化进程,仅2017年一年内基于BMKCloud分析发表的SCI可谓是层出不穷,捷报频传。 研究热点揭秘 截止到11月3日,已发SCI数达到19篇,研究对象涵盖了粮食作物、经济作物、人、家禽家畜、水产、林木、花卉、瓜果等类型。植物类占比52.63%,动物类占比42.11%,医学类占比5.26%。发表的杂志有Archives of Toxicology(IF=5.905)、Frontiers in Plant Science(IF=4.298)、Scientific Reports (IF=4.259)等。累计影响因子为65.3。 云平台对这些文章都做了详尽梳理解析,点击即可查看: OsbZIP46、SAPK6基因共同过表达提升水稻抗旱能力Frontiers in Plant Science IF:4.298 通过转录组分析揭示马杜霉素在原代鸡心肌细胞中的毒性Archives of Toxicology  IF:5.905  小麦与禾谷胞囊线虫在侵染初始阶段的转录组响应 Scientific Reports IF:4.259 芜菁花芽转录组比较分析解析多倍体异常减数分裂过程Frontiers in Plant Science IF:4.298 单侧眼组织缺损和视网膜裂症的家族遗传原因Sci Rep .IF:4.259 图1 2017年基于BMKCloud分析发表的部分SCI文章 产出单位分布 从通讯作者所属地域来看的话,这些云平台的用户主要分布在全国12个省、直辖市,其中产出最多的单位是南京农业大学,17年基于BMKCloud发表了3篇SCI(累计IF 13.23),短期连发多篇文章的老师也不少,南海水产所喻达辉老师4个月内连发两篇(累计IF 6.298),南京农业大学的张艳丽老师更是在短短2个月内连发两篇(累计IF 7.329) 图2 部分科研单位的地理分布 分析平台热度排行 百迈客BMKCloud已推出25款分析平台,包括真核生物有参转录组、真核生物无参转录组、小RNA、长链非编码RNA、微生物多样性、重测序、ChIP-Seq、CircRNA、全转录组联合分析等多个平台,覆盖转录调控、微生物多样性、个体重测序、群体个性化分析等多种组学分析项目,可满足所有分析高通量测序数据整合分析及深度挖掘的需求。 从2017年已发表的SCI文章来看,使用频率最高的是有参转录组分析平台,占比63%,lncRNA分析平台和小RNA分析平台占比11%,另外医学外显子分析平台、微生物多样性分析平台、甲基化分析平台各占比5%。 图3 百迈客分析平台使用热度一览 云课程 让生物信息更简单 BMKCloud部署了185+课时的线上培训,囊括RNA系列、DNA系列、生物信息软件使用和大家讲坛等课程,每款APP(分析平台)都有详细说明及步骤提示,无论有没有生物信息学基础,用户都能够在短期内玩转各类分析。这些作者中有42.86%的老师是完全不懂生信的,约有50%的老师略知一二,仅7.14%对生信分析较了解。 BMKCloud热门个性化分析 除了主流程APP之外呢,BMKCloud上也部署了104个分析工具、3款工具集。涵盖12大功能类别(基因分析、遗传进化、ncRNA、质控、组装、比对、数据提取、突变、统计、表格处理、绘图、序列分析等)。基本上可满足所有转录组学分析的需求。图4展示了19篇文章中热门的个性化分析 图4 BMKCloud个性化分析示例 仅从这19篇SCI内相关分析内容来看的话,使用次数最高的分析依次是GO分析、热图、KEGG分析、韦恩图和COG分析。其中热图和GO分析使用率最高,均达到52.63%(图5);此外共表达趋势分析,蛋白互作分析等分析也应用较多。 图5 BMKCloud个性化分析热度一览   BMKCloud的便捷度 BMKCloud使用便捷,助您多角度进行分析挖掘,随时随地抓住分析灵感。 这些研究者中有35.71%的老师选择在8:00-19:00进行分析,选择晚上进行分析的老师占比21.43%,凌晨进行分析的占比7.14%。 这些研究者中所有拥有BMKCloud年账号的用户会对同一批样本进行2-4次分析和挖掘,而实现这些仅需要修改参数、阈值或更换参考基因组等后重新投递任务即可。避免了传统科技服务模式中多环节沟通,多次收费,出错率高的问题。 BMKCloud的速度 几乎所有的研究者从拿到数据到分析完成所需周期都在3个月内。35.7%的研究者1个月内就完成了文章内的全部分析内容,其他的研究者拿到数据到发表文章的周期都控制在6个月内。其中,中国农科院油料研究所胡琼老师项目组使用有参转录组分析平台投递了200G数据,7天就完成了整个主流程及个性化分析。 BMKCloud的口碑 云平台分析和文章转化的高效赢得了15000+用户以及广泛的认可度,所有研究者都对我们的演示、售后及便捷度做出了很高的评价,并表示愿意推荐给身边的科研同事及圈内朋友。2017年借助BMKCloud取得累累硕果的研究者们对云平台评价如下: (1)简单方便(流程标准化,图形化,操作简单,使用快捷,节省时间); (2)灵活自由(参数设置可随意更改,可随时随地分析以及查看进度); (3)可重复(可多次分析,分析工具种类多样,可以帮助挖掘出新的分析点); (4)云课堂以及文献资源很实用; (5)数据储存空间大; (6)公共数据库是大亮点。 广大科研用户的肯定,科研探索的累累硕果是对我们最好的肯定与鼓励。百迈客将一如既往地秉承“成就客户”的服务理念,持续迭代开发“云模式”,继续引领科技服务2.0时代,更好地服务于广大科研工作者。也希望越来越多的课题组能够享受到云平台的便捷与高效,发表更多更优质的文章!   年末福利 年底我们对老客户发文有感恩活动,如果您基于平台分析发表了SCI文章并引用了BMKCloud或www.biocloud.net,即可获赠1个月现有账号使用权限!   BMKCloud有别于传统的“服务器+命令行”的生信分析模式,是一种人人都可以分析的新型分析模式。百迈客云平台代表了科技服务的新趋势,我们已步入云端可视化自主分析的基因科技服务2.0时代了! 为感恩回馈各位新老用户对百迈客云的支持与关注,百迈客云现推出72款生物信息分析小工具供免费使用,希望在科研路上助您一臂之力!注册百迈客云平台即可免费使用。...

Unraveling the genetic cause of a consanguineous family with unilateral coloboma and retinoschisis: expanding the phenotypic variability of RAX mutations  云平台助力解析单侧眼组织缺损和视网膜裂症的家族遗传原因   杂志: Scientific Reports  影响因子:4.259 PMID:28831107   研究背景: 眼组织缺损(OC)是一种常见的发育结构性缺陷,因视神经不完全闭合发展而成。通常在虹膜、脉络膜视网膜或视盘组织部分发生眼组织缺损。这一疾病进程由多种基因突变引起,具有遗传异质性和复杂性。理解OC发病分子机制及基因型与表型相关性,对于OC的分子诊断具有重要意义。 近期云平台助力温州医科大学黄秀峰老师在Scientific Reports发表了一篇关于眼组织缺损致病性突变的研究,通过对OC患者及其家族成员进行了全外显子组测序(WES),纯合子定位分析、全面变异分析,找到了血缘家族中单侧眼组织缺损与视网膜裂症的遗传原因,阐述了OC基因型表型相关性。跟随小编去看看吧。   材料和方法: 材料:采集患者及其家族成员血液,提取DNA   测序平台:Illumina HiSeq 2000   分析平台:BMKCloud,具体分析如下: 1.全基因组测序结果与人参考基因组进行比对(SOAPaligner),筛选SNV和Indel(SOAPaligner、GATK ) 2.SNP阵列分析纯合子定位 3.变异型分析 根据ExAC、NHLBI ESP、1000 Genome、dbSNP137数据库检索,对致病变异型进行初始筛选,并进行Sanger测序验证。评估候选变异型作用。分析错义突变体(SIFT、Polyphen-2、MutationTaster)。预测每个候选变异型多肽的拓扑模型(SMART、RaptorX、PyMol)另外进行序列比对(Clustal Omega),确定核苷酸保守性程度(PhyloP)   研究结果: 1.表型检测 患者是一名14岁男性,携带常染色体隐性遗传病,生于近缘家族,全面眼科检查显示其最近左眼视力逐渐丧失,双眼出现轻微白内障症状。眼前段裂隙灯检查显示右眼正常(图1A),但左眼6点钟方向虹膜缺损(图1D)。眼底检查右眼视神经盘缺损(图1B),左眼脉络膜缺损(图1E)。在视觉范围上,右眼生理盲区扩大(图1C),双眼视觉对比敏感度低(图1C,F)。其他家庭成员没有OC或眼部畸形。除OC表型外,眼底荧光血管造影(FFA)和光学相干断层扫描(OCT)结果提示患双侧视网膜血管炎和继发性视网膜劈裂症(图1G,H,K和L)。接受2个月激素治疗后,双眼黄斑水肿逐渐消失(图1I和J)。 图1 RAX基因突变患者临床特征 2.WES及纯合子定位分析表明RAX发生了突变 WES平均测序深度30X,平均覆盖率>95%. 变异型分析的详细过程(图2A)。变体数据库中排除次要等位基因频率(MAF)>0.005的变异型。纯合子分析显示有12个大于2 Mb的纯合子位点(图2B)。在使用逐步筛选策略后,发现在纯合子区域有3个候选变异型(RAX、TRPM5、PCDH17)。在家族成员中扩大检测后筛选得到一个RAX错义突变(c.113 T> C,p.I38T)。且在dbSNP137、1000 Genome、ESP6500、内部数据库、ExAC数据库中没有此错义突变(图3A-B)。 纯合子分析表明RAX基因位于18号染色体上较大的纯合子区域(18.19 Mb)(图2B)。 c.113 T> C变异使疏水性异亮氨酸转变为亲水性苏氨酸。该家族中基因型与表型分离(图3A,B):患者父母是杂合子携带者,其健康的姐妹无此突变。在38位的异亮氨酸高度保守(图3C),核苷酸位点c.113T> C的PhyloP得分为2.986。表明RAX基因的纯合子突变 c.113 T> C是OC的病因。 图2 变异分析流程及纯合子分析图谱 图3 RAX突变鉴定(A.谱系遗传;B.Sanger验证;C.RAX基因 Ile38 氨基酸残基进化保守性) 3.RAX突变比较分析和结构分析 患者在RAX基因的两个等位基因的exon1内有突变(图4A)。RAX c.113 T > C突变导致38位的氨基酸错义突变,氨基酸I38位于八肽模体中,生信分析预测显示I38T错义突变具有破坏性(SIFT得分:0; MutationTaster得分:1),可能对蛋白质,视网膜和前神经折叠同源序列造成损害(PolyPhen2得分:0.873)。子域识别结果表明这个突变不在同源序列、配对序列中(SMART)。对RAX进行潜在功能分析(RaptorX,PyMol),构建结构模型(图4B),结合遗传分析及结构分析结果,我们推测序列改变增加了I38残基的螺旋构象,从而使RAX蛋白空间结构改变。 图4 RAX基因突变结构分析 结果: 在本研究中,我们对UC和视网膜裂症患者及其家庭成员进行了分子遗传学检测。利用全外显子组测序和纯合性定位组合方法,筛选后得到唯一候选致病基因RAX致病突变(c.113 T> C,p.I38T)。结果得到临床、功能模型数据的支持。本研究扩大了对RAX突变表型变异性的理解,有利于进一步理解OC发病分子机制及基因型表型相关性。   创新点: 本研究不仅证明了WES和纯合子绘图是解析遗传异质性疾病患者致病性突变的有用工具,也有助于进一步阐明OC疾病发病分子机制及基因型表型相关性。   参考文献:Huang X F, Huang Z Q, Lin D, et al. Unraveling the genetic cause of a consanguineous family with unilateral coloboma and retinoschisis: expanding the phenotypic variability of RAX mutations.[J]. Sci Rep, 2017, 7(1)....

ranscriptome Analysis of Floral Buds Deciphered an Irregular Course of Meiosis in Polyploid Brassica rapa  云平台助力芜菁花芽转录组比较分析解析多倍体异常减数分裂过程 杂志: Frontiers in Plant Science  影响因子:4.298 PMID:28553302 云平台助力郑州大学的老师们发表了一篇关于四倍体与二倍体芜菁转录组比较分析的文章,这是第一个同时在细胞和转录组水平证明染色体多倍性对减数分裂过程有不利影响的研究,有助于全面了解二倍体和多倍体芜菁减数分裂时花芽转录组的一致性和差异性。   研究背景: 多倍体化对于动植物进化及物种形成具有重要意义。染色体倍性增加在大多数植物中十分常见,约30-80%的被子植物在进化过程中历经这一过程,种内基因组复制可产生同源多倍体,种间杂交可产生异源多倍体。理解多倍体对植物繁殖的影响对于多倍体育种项目具有重要意义。 本文在细胞和分子水平对同源四倍体、二倍体芜菁的生殖组织(未成熟的花芽)进行了转录组比较分析。先进行细胞学分析,发现结果显著异常后,进行RNA测序分析,进一步阐明这种细胞学差异的分子机制。   材料和方法: 材料:温室中培育芜菁同源四倍体和二倍体,采集1-1.5mm的花芽。 细胞学分析:花芽制片观察染色体行为。 转录组分析:二倍体花芽样本T1,T2,T3,同源四倍体样本T4,T5,T6进行测序,筛选DEGs,进行功能注释及相关验证。   测序平台:Illumina HiSeq 2000 分析平台:BMKCloud,具体分析如下: 测序原始数据去接头、低质量序列后与参考基因组比对(TopHat2),对比对上的序列进行转录本重构,评估表达水平(Cufflinks)。鉴定DEGs(logFC≥2,FDR<0.01)并与Nr、UniProtKB\Swissprot、KOG、COG、GO、KEGG等数据库进行比对注释。   研究结果: 1.同源四倍体芜菁异常的减数分裂过程 减数分裂过程中染色体行为分析结果显示:二倍体中,同源染色体中期配对为二倍体,后期分离(图1A-D)。四倍体中,在减数分裂中期I形成多倍体和单倍体,在后期I、II发生不均等分离(图1E-L)。 图1 四倍体芜菁异常染色体行为 统计结果表明与二倍体相比,芜菁多倍体化后整个减数分裂过程发生了不规则变化。 图2 四倍体芜菁减数分裂中异常染色体行为分析 2.测序结果 每个样本文库得到29.07 GB clean read,Q30≥88.55% 。同源四倍体、二倍体比对到参考基因组效率分别为74.34%,75.88%。样本组的表达模式分析共得到40927个基因,其中有1001个新基因。 3.DEGs分析 40927个基因中,4601个基因是差异表达的(图3A),2343个上调,2259个下调(图3B),两组间DEGs占比较少((11.24%)。随机挑选DEGs进行层次聚类,分析表达模式(图3C),发现四倍体中K1、K3、K6、K9类别下调,K2、K4、K5、K7,K8类别显著上调。 图3 DEGs表达分析 与COG、GO、KOG、Nr、KEGG、Swiss-Prot比对后,97%的DEGs至少能与一个数据库比对上(表1) 表1 差异表达Unigene注释 COG分类中,2615个DEGs可归到25个COG类别(图4),其中包含最多DEGs的类别为:复制重组与修复(222个,8.49%),转录(261个,9.98%),一般功能(515个,19.69%),与KEGG结果一致。 GO分类中大多数DEGs聚类到细胞过程(61.97%),繁殖过程(15.49%),结合过程(42.96%)(图5)。 KEGG功能注释结果表明,4,601个DEGs中有1,453个可归到50个生物途径。高度富集的通路为:氨基酸合成(48个,3.3%),植物激素信号转导(61个,4.2%)和蛋白质加工内质网途径(47个,3.2%)(图6)。 图4 COG分类 图5 新DEGs的GO分类 图6 KEGG富集通路 4.减数分裂相关DEGs分析 为阐述四倍体芜菁减数分裂过程改变相关的遗传信息,选择以下明显与减数分裂相关的COG类别进行细分:(1)复制,重组和修复;(2)染色质结构和动力学;(3)细胞周期控制,细胞分裂,染色体分区;(4)细胞骨架。在4,601个位点中,共鉴定出288个与减数分裂相关的基因(图7A)。 应用拟南芥数据库(TAIR)的BLASTN搜索,进一步鉴定芜菁中减数分裂直系同源基因,11个两组间差异表达的已知减数分裂基因。如图7所示,显著富集的L组(复制,重组和修复)包含121个上调、102个下调基因,这组聚类结果(图8B)显示已知减数分裂基因(包括RAD54,DMC1和RPA)在四倍体芜菁中表达显著下调。D组(细胞周期控制,细胞分裂和染色体分配)中,10个下调,13个上调,已知的减数分裂基因DIF1/SYN1和CyclinA1-2/TAM上调(图7C)。在B组(染色质结构和动态变化)(图7D),有7个基因上调,包括已知调节开花时间和植物发育的NF-YB8基因,14个基因下调,包括组蛋白相关基因(H2AV和H4)。Z组(细胞骨架)中,8个基因上调,7个基因下调(图7E),且没有已知减数分裂相关基因。已知减数分裂基因MMD1和HOP2在一般功能组(R)中富集(图4)。HOP2,XRI1,ZYP1a和ZYP1b未归到任何COG组。 KEGG比对后将288个中的108个DEGs分配到28个途径,包括减数分裂相关同源重组(2.8%),DNA修复和重组(6.5%),DNA复制(3.7%),染色体和相关蛋白(4.6%)。 图7 COG分类中减数分裂相关基因分布及表达分析 同源重组对于真核生物减数分裂中SPO11蛋白质引起的DSBs的精确修复十分重要,作为调节关键介质的DMC1,RAD54和RPA在这个途径中显著下调。Cyclin-A1-2 / TAM是DNA复制途径中唯一与细胞周期蛋白有关的基因,表达上调。这些结果表明染色体集在同源四倍体中增加了一倍,与减数分裂有关的大多数基因可能通过上调或至少正常表达以满足减数分裂期间基因组含量增加的需求。否则,减数分裂关键基因下调可能会干扰四倍体减数分裂中染色体行为。 qRT-PCR 对24个减数分裂相关基因进行验证,20多个基因表达与测序结果一致。对拟南芥中直系同源基因分析结果显示两个物种减数分裂相关基因的总体表达模式相近。   结果: 本研究细胞学分析发现,合成多倍体整体减数分裂过程是显著不规则的。为从分子水平阐明这一过程的遗传基础,本文进行了同源四倍体和二倍体芜菁花芽间的遗传调节差异的转录组比较分析。在40927个表达基因中,同源四倍体芜菁花芽中鉴定出4601个DEGs,其中288个与减数分裂相关。DMC1(已知减数分裂特异性基因,与DSBs同源染色体依赖性修复相关)在同源四倍体芜菁中表达显著下调,可能与减数分裂I期的异常有关。某些RNA解旋酶,细胞周期、体细胞DNA修复相关的DEG在基因组复制后表达上调,减数分裂DSB修复有关基因表达显著下调。芜菁和拟南芥中减数分裂相关基因的总体表达模式相近。   创新点: 这是第一个在细胞和转录组水平证明染色体多倍性对减数分裂过程有不利影响的研究,有助于全面了解二倍体和多倍体芜菁减数分裂时花芽转录组的一致性和差异性。   参考文献:Braynen J, Yang Y, Wei F, et al. Transcriptome Analysis of Floral Buds Deciphered an Irregular Course of Meiosis in Polyploid Brassica rapa[J]. Frontiers in Plant Science, 2017, 8:768....