诺禾致源-ChIP-seq 流程
ChIP-seq
染色体免疫共沉淀（ChIP）是一种用于研究蛋白质与DNA的体内相互作用的经典实验技术。采用特异性抗体将目的蛋白进行免疫沉淀，由此可以把目的蛋白所结合的基因组DNA片段也富集下来。通过与高通量测序技术的结合，对ChIP后的DNA产物进行测序分析，从全基因组范围内寻找目的蛋白的DNA结合位点，以高效率的测序手段得到高通量的数据结果。
1.1C hIP免疫沉淀实验流程
目前主要有两种不同的ChIP实验方法，大致流程如下（均以细胞样品的处理过程为例）：
Cross-liking Chromatin Immunoprecipitation (X-ChIP)
1.甲醛处理细胞，使DNA-protein的相互结合作用被交联固定。
2.裂解细胞，得到全细胞裂解液。
3.超声处理，将基因组DNA打断至100-500 bp。
4.抗体免疫沉淀：在细胞裂解液中加入一抗和beads，并进行孵育。
5.采用合适的实验条件进行洗脱，并解交联。
6.通过qPCR对ChIP结果进行验证。
7.准备好的ChIP后的DNA样品可以用于ChIP Sequencing建库。
Native Chromatin Immunoprecipitation (N-ChIP)
1.通过非变性的方式得到核裂解液。
2.微球菌核酸酶（Micrococcal nuclease）消化染色质，得到单核小体或核小体
寡聚体。
3.抗体免疫沉淀：在细胞裂解液中前后加入一抗和beads，并进行孵育。
4.DNA分离。
5.通过qPCR对ChIP结果进行验证。
6.准备好的ChIP后的DNA样品可以用于ChIP Sequencing建库。
下面步骤由我们公司做：
1.2C hIP Sequencing文库构建流程
1.DNA片段末端修复,3’端加A碱基，连接测序接头（详细步骤请参考Illumina 公司Paired-End DNA Sample Prep kit）。
2.PCR扩增及DNA产物的片段大小选择（一般为100-300 bp，包括接头序列在内）
合格的文库用于上机测序。
诺禾致源高分文章集锦-植物转录组
温带和热带莲根状茎形成过程中的转录组分析 Transcriptomic Analysis of the Regulation of Rhizome Formation in Temperate and Tropical Lotus (Nelumbo nucifera ) 研究对象：莲根状茎 期刊：Scientific Reports 影响因子：5.578 合作单位：中国科学院武汉植物园 发表时间：2015年7月 摘 要 Rhizome is the storage organ of lotus derived from modified stems. The development of rhizome is a complex process and depends on the balanced expression of the genes that is controlled by environmental and endogenous factors. However, little is known about the mechanism that regulates rhizome girth enlargement. In this study, using RNA-seq, transcriptomic analyses were performed at three rhizome developmental stages—the stolon, middle swelling and later swelling stage —in the cultivars ‘ZO’ (temperate lotus with enlarged rhizome) and ‘RL’ (tropical lotus with stolon). About 348 million high-quality reads were generated, and 88.5% of the data were mapped to the reference genome. Of 26783 genes identified, 24069 genes were previously predicted in the reference, and 2714 genes were novel transcripts. Moreover, 8821 genes were differentially expressed between the cultivars at the three stages. Functional analysis identified that these genes were significantly enriched in pathways carbohydrate metabolism and plant hormone signal transduction. Twenty-two genes involved in photoperiod pathway, starch metabolism and hormone signal transduction were candidate genes inducing rhizome girth enlargement. Comparative transcriptomic analysis detected several differentially expressed genes and potential candidate genes required for rhizome girth enlargement, which lay a foundation for future studies on molecular mechanisms underlying rhizome Formation. 关键词 根状茎；变态发育； DGE 研究背景 莲根状茎，即莲藕，作为一种变态茎，是莲 的贮藏器官。根状茎的发育是一个复杂的过 程，受到与环境及内源因素调控的基因平衡 表达的影响。关于根状茎膨大的调控机制很 少为人所知。
全基因组重测序探索刚地弓形虫致病基因
首页 科技服务 医学检测 科学与技术 市场与支持 加入我们 关于我们 安徽医科大学研究人员携手诺禾致源重测序团队，通过对2种刚地弓形虫的全基因组重测序变异检测研究， 从基因组水平解释了2种虫株产生表型差异的原因，为弓形虫病的治疗和疫苗研发提供了理论依据。 该研究成果发表于2015年10月的BMC Genomics杂志（IF: 3.986）。 研究背景 刚地弓形虫（Toxoplasma gondii Nicolle&Manceaux, 1908）寄生于人和许多种动物的有核细胞，但只能在猫科动物的肠道内繁衍，能够引起人畜共患的弓形虫病。与北美和欧洲群体遗传结构不同，Chinese 1 (ToxoDB#9)是中国的弓形虫优势基因型。在Chinese 1型弓形虫中，Wh3（强毒株）和Wh6（弱毒株）对小鼠表现出了不 同的毒力。本研究拟通过全基因组重测序技术，从基因组水平探究两种虫株表型及致病性差异的原因。 研究结果 1 SNP、indel检测及注释与参考基因组比对发现，Wh3中共有505,856个SNPs，30,004个indels；Wh6中共有505,654个SNPs，30,658个indels。进一步分析两样本特有变异，发现Wh3中特有SNP和indels分别位于2847和2452个基因中，Wh6中特有SNP和indels分别位于2868和2613个基因中（图1，图3）。 图1 SNPs、indels的对比分析及分布情况统计 注：a为SNP韦恩图；b为indels韦恩图；c为SNP突变型分布情况；d为编码区indels长度分布情况 图2 CNVs（左）和SVs（右）的分布情况统计 图3 Wh3 （左）和Wh6（右）的全基因组变异情况 图4 I、II、III型弓形虫与Chinese 1型弓形虫的ROP16和GRA15序列比对分析 图5 三种弓形虫虫株的表达模式分析 NGS项目文章 全基因组重测序探索 探索刚地弓形虫毒力相关基因 BMC Genomics 2 CNV、SV检测及注释与参考基因组比对发现，Wh3中共有2320个SVs，1942个CNVs；Wh6中共有4661个SVs，3080个CNVs。其中，Wh3含有85个片段插入（总长度：282,700bp），2995个片段缺失（总长度：4,940,000 bp）；而Wh6含有90个片段插入（总长度：328,800bp），1852个片段缺失（总长度：7,157,700 bp）（图2，图3）。 3 毒力相关因子的变异信息分析通过对与弓形虫毒力和侵染性相关的一系列关键因子（R O P s 、GRAs、MICs、RONs和SAGs）的变异信息分析，发现与其他影响因子相比，GRA3和RON3的编码基因中含有更多的SNPs和indels；其中，G R A 3编码基因含有35个SNPs和2个indels，RON3编码基因含有89个SNPs和6 个indels。同时，与I、II和III型弓形虫相比，Chinese1型弓形虫的ROP16和GRA15表现为多态性的ROP16I/III 和GRA15 II（图4）。 4 qRT-PCR分析 为探究与Wh3和Wh6表型差异相关的基因，分别对Wh3、Wh6和RH这三种虫株进行qRT-PCR分析。与强毒株Wh3相比，发现在弱毒株Wh6中，GRA3和RON3的基因表达量显著上调，而ROM4, profilin, M2AP, AMA1, RON2, RON3和RON4的基因表达量显著下调。 与参考基因组虫株RH相比，在Wh3和Wh6中的SRS9, ROP8,MIC8和RON5的基因表达量均上调，而SAG1, ROP5和ROP18的基因表达量均下调（图5）。
诺禾致源高分文章集锦-植物基因组
陆地棉基因组测序揭示四倍体棉进化与纤维发育机制Sequencing of allotetraploid cotton (Gossypium hirsutum L. acc. TM-1) provides a resource for fiber improvement 研究对象：陆地棉遗传标准系TM-1 期刊：Nature Biotechnology 影响因子：41.514 合作单位：南京农业大学 发表时间：2015年4月 摘 要 Upland cotton is a model for polyploid crop domestication and transgenic improvement. Here we sequenced the allotetraploid Gossypium hirsutum L. acc. TM-1 genome by integrating whole-genome shotgun reads, bacterial artificial chromosome (BAC)-end sequences and genotype-by-sequencing genetic maps. We assembled and annotated 32,032 A-subgenome genes and 34,402 D-subgenome genes. Structural rearrangements, gene loss, disrupted genes and sequence divergence were more common in the A subgenome than in the D subgenome, suggesting asymmetric evolution. However, no genome-wide expression dominance was found between the subgenomes. Genomic signatures of selection and domestication are associated with positively selected genes (PSGs) for fiber improvement in the A subgenome and for stress tolerance in the D subgenome. This draft genome sequence provides a resource for engineering superior cotton lines.关键词 陆地棉；de novo；四倍体 研究背景 陆地棉（Gossypium hirsutum L.）隶属锦葵目（Malvales），锦葵科（Malvaceae），棉属（Gossypium），因最早在美洲大陆种植而得名，是世界上最重要的棉花栽培品种，占全球棉花种植面积的90%以上。尽管陆地棉在棉花产业中占据核心地位，但由于其为异源四倍体，相关的全基因组测序工作一直难以开展。来自南京农业大学、北京诺禾致源、美国德克斯大学的国际团队，利用最新测序技术，成功构建了高质量的陆地棉全基因组图谱，为进一步改良棉花的农艺性状提供了基础，同时也为多倍体植物的形成和演化机制提供了新的启示。
诺禾致源真核无参转录组生物信息分析结题报告2013年8月
真核无参转录组生物信息分析结题报告 建库测序流程 Total RNA样品检测 文库构建 上机测序 F:/结题报告+老销售培训/结题报告模板修改/…/真核无参转录组_Report.html 1/38
2/38 F:/结题报告+老销售培训/结题报告模板修改/…/真核无参转录组_Report.html 北京诺禾致源生物信息科技有限公司 一、建库测序流程 从RNA 样品到最终数据获得，样品检测、建库、测序每一个环节都会对数据质量和数量产生影响，而数据质量又会直接影响后续信息分析的结果。为了从源头上保证测序数据的准确性、可靠性，诺禾致源对样品检测、建库、测序每一个生产步骤都严格把控，从 根本上确保了高质量数据的产出。实验流程图如下： 1 Total RNA 样品检测 诺禾致源对RNA 样品的检测主要包括4种方法：(1) 琼脂糖凝胶电泳分析RNA 降解程度以及是否有污染(2) Nanodrop 检测RNA 的纯度（OD260/280比值）(3) Qubit 对RNA 浓度进行精确定量(4) Agilent 2100精确检测RNA 的完整性 2 文库构建及库检 样品检测合格后，用带有Oligo （dT ）的磁珠富集真核生物mRNA （若为原核生物，则通过试剂盒去除rRNA 来富集mRNA ）。随后加入fragmentation buffer 将mRNA 打断成短片段，以mRNA 为模板，用六碱基随机引物（random hexamers ）合成一链cDNA ，然后加入缓冲液、dNTPs 、RNase H 和DNA polymerase I 合成二链cDNA ，随后利用AMPure XP beads 纯化双链cDNA 。纯化的双链cDNA 再进行末端修复、加A 尾并连接测序接头，然后用AMPure XP beads 进行片段大小选择，最后进行PCR 富集得到最终的cDNA 文库。文库构建完成后，先使用Qubit2.0进行初步定量，稀释文库至1ng/ul ，随后使用Agilent 2100对文库的insert size 进行检测，insert size 符合预期后，使用Q-PCR 方法对文库的有效浓度进行准确定量（文库有效浓度 ＞2nM ），以保证文库质量。文库构建原理图如下：
转录组测序
真核mRNA测序是基于HiSeq平台，对真核生物特定组织或细胞在某个时期转录出来的所有mRNA进行测序，既可研究已知基因，亦能发掘新基因，全 面快速地获得mRNA序列和丰度信息。真核mRNA测序方法可以分为：有参考转录组、无参考转录组以及数字基因表达谱（DGE）三大类。 技术参数 案例解析 ［案例一］ mRNA和small RNA转录组揭示新合成异源六倍体小麦杂种 优势的动态部分同源调控 诺禾致源携手中国农业科学院作物科学研究所，利用转录组测序技术，对杂交亲本、新合成异源六倍体小麦的幼苗、穗和种子进行了mRNA和smallRNA测序及信息分析，发现新合成异源六倍体小麦绝大部分基因表现为12类基因表达模式，包括加性表达，少部分的基因表现为非加性，基因的非加性表现出非常强的发育时期特异性，与生长势密切相关；miRNA的丰度随着倍性的增加逐渐下降，新合成异源六倍体小麦中非加性表达的 miRNA也同样表现出亲本显性表 达，miRNA的表达敏感性与生长势和适应性密切相关。该研究揭示了不同倍性 非对等杂种优势的分子基础。 ［案例二］ 磷酸三（2,3-二氯丙基）酯（TDCPP）对四膜虫生长繁殖的 抑制作用与核糖体相关 诺禾携手华中农业大学,利用转录组测序和信息分析技术,研究了TDCPP处理组和对照组差异基因表达,并对差异表达基因进行KEGG通路分析,发现核糖体基因通路显著富集, 同时伴随胞浆和粗面内质网上核糖体数量减少体积增大。这些探索表明四膜虫可以作为TDCPP反应的生物指标，为后续研究TDCPP作用其他生物的毒理机制提供了新视角。 ［案例三］ 转录组揭示寄主植物与宿主之间进行RNA交换的机制 参考文献 菟丝子被称作勒死草，会用被称作吸根的专用器官穿透宿主组织与其建立联系，可以吸取宿主的水份与营养物质，也能吸取RNA(mRNA)分子。本研究分别选取菟丝子和拟南芥及番茄的共生体茎上的三段组织进行转录组学的研究，发现寄生植物与寄主之间mRNA的转移量很大且是一种双向转移的模式；两种宿主相比，更多的拟南芥RNA被转移到菟丝子植物之中，而且菟丝子与拟南芥之间较自由的交换，可表明调节菟丝子吸根选择性的机制可能是宿主特异性的，从而揭示了寄主与宿主之间进行RNA转移的遗传机制。 [1] Li A, Liu D, Wu J, et al . mRNA and small RNA transcriptomes reveal insights into dynamic homoeolog regulation of allopolyploid heterosis in nascent hexaploid wheat [J]. The Plant Cell, 2014: tpc. 114.124388.[2] Jing Li, John P , Giesy, Liqin Yu, et al . Effects of Tris (1,3-dichloro-2-propyl) Phosphate (TDCPP) in Tetrahymena Thermophila: Targeting the Ribosome. Scientific Reports. 2015, 5:10562. [3] Kim G, LeBlanc M L, et al . Genomic-scale exchange of mRNA between a parasitic plant and its hosts [J]. Science, 2014, 345(6198): 808-811. 图1 非加性表达miRNA与亲本显性表达miRNA的 等级聚类分析和两者的关联 图2 显著富集的KEGG通路 图3 菟丝子与拟南芥、番茄转移RNA和非转移RNA的表达和富集分析 样品要求文库类型测序策略数据量类型 分析内容 项目周期 真核有参转录组测序 真核无参转录组测序 6 Gb、8 Gb、10 Gb、12 Gb clean data 6 M clean reads 3 Gb clean data 项目数据至少12 Gb clean data 数字基因表达谱（DGE） HiSeq PE150 HiSeq PE150 HiSeq SE50HiSeq PE125普通转录组文库； 链特异性转录组文库 40天50天30天 35天（有参）45天（无参） RNA样品总量≥1.5 μg； RNA样品浓度≥50 ng/μL 参考基因组比对 新转录本预测可变剪切分析SNP/InDel分析 基因表达水平分析RNA-seq整体质量评估 转录因子注释GO/KEGG富集分析蛋白互作网络分析基因共表达网络构建可视化结果展示 参考转录组拼接 转录本/Unigene长度统计 基因功能注释NR，NT，Swiss Prot GO，KEGG，KOG Protein Family CDS预测分析SNP/SSR分析
诺禾致源高分文章集锦-动物基因组
川金丝猴全基因组测序解析其植食性机制与进化史 Whole-genome sequencing of the snub-nosed monkey provides insights into folivory and evolutionary history 研究对象：川金丝猴期刊：Nature Genetics 影响因子：29.352 合作单位：中国科学院动物研究所发表时间：2014年11月 摘 要 Colobines are a unique group of Old World monkeys that principally eat leaves and seeds rather than fruits and insects. We report the sequencing at 146× coverage, de novo assembly and analyses of the genome of a male golden snub-nosed monkey (Rhinopithecus roxellana ) and resequencing at 30× coverage of three related species (Rhinopithecus bieti , Rhinopithecus brelichi and Rhinopithecus strykeri ). Comparative analyses showed that Asian colobines have an enhanced ability to derive energy from fatty acids and to degrade xenobiotics. We found evidence for functional evolution in the colobine RNASE1 gene, encoding a key secretory RNase that digests the high concentrations of bacterial RNA derived from symbiotic microflora. Demographic reconstructions indicated that the profile of ancient effective population sizes for R. roxellana more closely resembles that of giant panda rather than its congeners. These findings offer new insights into the dietary adaptations and evolutionary history of colobine primates. 关键词 金丝猴；重测序；植食性；进化 研究背景 金丝猴（Rhinopithecu spp.)隶属于灵长目（Primates）、疣猴亚科（Colobinae）、仰鼻猴属 （Rhinopithecus ） ，目前共有5个种，即川、滇、黔、缅甸和越南金丝猴。通过对金丝猴基因组以及肠道基因组进行全面系统的研究，解析了金丝猴的植食性分子遗传机制，为了解疣猴亚科的系统进化、功能适应性奠定了遗传基础，同时开展了仰鼻猴属的进化历史和遗传多态性分析。
群体重测序与大山雀的群体进化
诺禾致源最新“高性价比”群体进化研究成果 继2013年合作完成地山雀基因组测序之后，北京诺禾致源重测序事业部团队与中国科学院动物研究所研究人 员再次携手，通过对13个地区的32只大山雀进行全基因组重测序，解析了喜马拉雅山脉东部大山雀对随海拨 变化的气候的适应机制。研究成果发表于2015年9月的Scientific Reports杂志（IF：5.578）。其中，中国科 学院动物研究所屈延华研究员、诺禾致源田仕林为论文的共同第一作者。 群体重测序揭示 大山雀适应随海拨而变的气候的机制 NGS项目文章
研究背景 大山雀（Parus major ）隶属于雀形目（Passeriformes）山雀科（Paridae）山雀属（Parus spilonotus ）。在东亚，主要分布在低海拔地区，也有部分种群生活在喜马拉雅山脉东部的高海拔地区。有季节性迁徙习性，在海拔4000m地区繁殖，在海拔2000m地区越冬。本研究采用群体重测序技术，从基因组水平上揭示了喜马拉雅山脉东部大山雀的起源及对这种季节性的、随海拨高度变化的气候的适应性机制。 研究方法 基于Illumina HiSeq 2000 测序平台，对来自13个地区的32只大山雀进行全基因组重测序，其中，11只大山雀来自喜马拉雅山脉东部地区，11只来自中国中/东部地区，10只来自内蒙古和蒙古，测序深度5X/样。以近缘物种地山雀（Pseudopodoces humilis ）基因组作为参考基因组，对大山雀群体进行了遗传多样性、种群历史动态、选择消除等分析。研究结果 1. 大山雀的群体进化分析 群体遗传多样性分析表明，来自蒙古（MON）、东喜马拉雅山脉（EH） 和中国中/东地区（CE）的大山雀各自聚为一类。EH和CE的亲缘关系较MON 更近。大山雀和地山雀约在5.8-13.3百万年前发生了分化；0.7-2.8百万年前 大山雀中分化出了MON分支；0.4-1.9百万年前EH和CE发生了分化。 2. 大山雀种群历史动态分析 EH、CE和MON在0.3-0.4百万年前种群遭遇了瓶颈效应。CE在0.06百万 年前有效群体大小迅速扩张，EH在同一时间开始扩张，但变化较平缓。而 MON经历了一个漫长的瓶颈，一直持续到0.02百万年前，在末次盛冰期有效 群体大小稍有增加。以上结果表明，EH种群动态受冰期气候影响小，这是由 于东喜马拉雅山脉局部环境相对稳定，长期生活在该地区的大山雀发生了高 海拔适应性进化。 3. 大山雀对高海拔的适应机制 通过选择消除分析，在东喜马拉雅山脉大山雀检测到183个基因受到强选 择，主要参与能量代谢过程和低氧反应。 （1）能量代谢基因进化：东喜马拉雅大山雀糖类代谢相关的基因发生了 快速进化，包括氨基糖和核苷酸糖代谢（5个基因）和胰岛素信号途径（11个 基因）。大山雀这种适应机制适合其在中、高海拔季节性迁徙的生活习性。 （2）低氧适应：低氧反应基因富集在MAPK信号通路，调节东喜马拉雅 大山雀体内相关基因的表达来适应低氧环境。 （3）形态进化：东喜马拉雅大山雀骨骼发育相关的基因发生了快速进 化，与分布在中国中/东低地区域的大山雀相比，体型更大，有利于保存热 量。此外，体型大倾向于有较高的氧亲和力，这也是东喜马拉雅大山雀能够 更好地调节体温适应高海拔的寒冷气候的一个原因。参考文献 Qu YH, Tian SL, et al. Genetic responses to seasonal variation in altitudinal stress: whole-genome resequencing of great tit in eastern Himalayas. Scientific Reports, 2015. 图1 大山雀遗传多样性及分化时间 图2 调控低氧反应的受选择基因
De novo 转录组揭示真菌与微生物的共生机制
然而，利用以上研究数据，另有文章发表： NGS项目文章 猪苓，非褶菌目多孔菌科树花属，是一种药用真菌； 蜜环菌，伞菌目蜜环菌属，夏秋季常寄生于树丛根部，是一种根腐菌。 猪苓的生长需要蜜环菌提供营养，二者之间存在共生关系。 在共生关系建立的过程中，蜜环菌的菌索会吸附并侵染猪苓的菌核， 而猪苓则会激发自身防御机制，如菌核细胞壁不规则地增厚以抵抗蜜环菌进一步侵染。 这个抵抗侵染过程的分子机制是怎样的？ 2015年10月，猪苓对蜜环菌的防御机制研究成果发表在Scientific Reports上， 该研究由中国医学科学院药用植物研究所郭顺星教授课题组负责， 其中De novo 转录组测序工作由诺禾致源完成。 研究内容与结果 首页 科技服务 医学检测 科学与技术 市场与支持 加入我们 关于我们 De novo 转录组 揭示猪苓对蜜环菌的防御机制 3. 差异基因功能分析 利用DESeq软件进行CT和CK的差异基因筛选，得到10,933个DEG，包括 8,780个上调基因和2,213个下调基因（图2）。 表达上调的差异基因GO富集分析发现：在molecular function部分，富集 最明显的是苏氨酸型肽链内切酶活性、肽酶活性以及特异性DNA结合等 term；在cellular component部分，富集显著的是蛋白酶复合体、膜组成 等term；而在biological process部分，蛋白质折叠、三羧酸循环以及染色 质组装等term下显著富集。而下调的DEG被富集到对有机物反应、细胞组 分合成等功能term，这预示着猪苓受到侵染时会限制自身的生长与繁殖。 对DEG进行KEGG富集分析发现：精氨酸和脯氨酸代谢通路富集显著；免疫 相关基因精氨酸酶（ARG）显著上调，推测其与猪苓抵抗蜜环菌侵染相关。 而在核糖体合成途径通路中，含有11个下调基因，这与GO富集中发现下调 基因在生物学过程相关term富集的结果是一致的，从而推测猪苓受到侵染 会减弱自身的生物学过程以保存能量。 2. 转录本功能注释 对拼接得到的转录本进行功 能注释，采用Nr, Nt, Swiss -Prot, KEGG, GO, KOG 和 Pfam 七大数据库。并对 unigene进行GO和KOG功 能分类（图1和图2）。 1. 转录本拼接 四个文库测序下机的原始 数据进行质控，将各样本 的clean reads进行混合拼 接，得到38444条 unigene，作为参考序列用 于后续分析。 4. qPCR验证差异表达基因 选取13个DEG进行实时荧光定量PCR，检测结果与RNA-seq一致。 结论与讨论 本研究首次通过RNA-seq研究蜜环菌侵染和未侵染的猪苓样本，筛选出差异表达基因10933个； 通过GO富集和KEGG富集分析，鉴定到了可能与抵御反应相关的若干差异基因。 其中，且热休克蛋白（Hsps）基因、抗氧化防御相关基因、凝集素基因、致病相关蛋白（TLP）基因、 次级代谢相关基因、细胞壁水解与融合相关基因、PDR家族基因、WD40蛋白基因等 在蜜环菌侵染的猪苓样本中呈现表达上调，均可能与抵御反应相关。 图3 转录数据的SSR分析 图2 对差异基因进行筛选 图1 unigene的GO功能分类 郭顺星教授课题组继续利用本次RNA-seq的数据结 果进行SSR分析，丰富了药用猪苓的因地域而异的 遗传多态性，并为猪苓群体遗传和行为生态学研究 提供了依据。 阅读原文>> 阅读原文>> Scientific Reports
基于 GBS 技术的玉米高密度遗传图谱构建及 QTL 定位
诺禾致源重测序 基于 GBS 技术的玉米高密度遗传图谱构建及 QTL 定位 BMC Genomics 首页 科技服务 医学检测 科学与技术 市场与支持 加入我们 关于我们 中国农业科学院作物研究所研究人员携手诺禾致源重测序团队， 采用 GBS 技术，对314株高世代群体（RILs）进行低深度测序，检测 SNP，开发 Bin 标记，构建高密度遗传图谱， 并对株型相关性状进行了定位，筛选出候选基因。 研究成果发表于2016年3月的 BMC Genomics 杂志（IF：3.986）。 玉米株型相关性状（株高、穗高等）与玉米产量、抗倒伏等密切相关，研究其遗传特点具有重要意义。 现存的低密度遗传图谱限制了 QTL 作图的高效性和准确性。基于二代测序的 GBS 技术已成为一种构建高密度遗传图谱的有力工具。 本文采用 GBS 技术与高世代作图群体（RILs）相结合，极大的提升了 QTL 作图的有效性以及对复杂农艺性状的研究。 研究背景 NGS项目文章 实验材料和方法 研究结果 取 材 测序技术测序平台 Ye478（母本）和 Qi319（父本）RILs (F11代)，314株 亲本全基因组重测序，测序深度分别为29.5X 和33.7X; 子代 GBS 测序，平均测序深度约0.07X 亲本（Illumina HiSeq 2000 PE125测序），子代（Illumina HiSeq 2500 PE125测序） 1 亲本重测序和子代 GBS 测序 Ye478（母本）和 Qi319（父本）进行全基因组重测序，测序深度分别为29.5x 和33.7x，与参考基因组（B73RefGen_V3）进行比对后，分别获得678,819,425和803,698,828条reads，亲本间纯合且有差异的 SNP 共有3,549,088个（图1）。 子代 RILs 群体进行 GBS 测序，共获得137,699,000条 reads，平均每个个体357,376条 reads，相当于玉米基因组大小的0.07x。通过标记筛选，剩余88,268个 SNPs 可用于 bin 标记开发。 2 基于bin 标记的遗传图谱构建 以100Kb大小为间隔，若相邻的100Kb间隔在整个 RILs 群体中都没有发生重组事件，则相邻的100Kb间隔被连接成 bins，这些 bins 被当做遗传标记（markers），共开发出4,183个 bin 标记。基于这些 bin 标记构建玉米高密度遗传图谱，遗传距离总长为1,545.65cM，标记间平均遗传距离为0.37cM，平均物理间隔0.51Mb（图2，表1）。 3 图谱质量评估 将遗传图谱与玉米参考基因组进行共线性分析，获得了良好的结果。为评估该图谱用于性状定位的效果，使用该图谱对玉米棒颜色性状（cob color）进行了 QTL 分析。在1号染色体上检测到QTL qC1，该 QTL 所在区域包含已克隆的基因 pericarp color 1(P1)，LOD 值为80.78。 4 株型性状相关的Bin 标记 株型相关的三个性状：株高（PH），穗高（EH）和节间数量（IN）都具有较高的遗传力，以 PH 遗传力最高，为94.87%，其次为 EH 和 IN 分别是85.24%和82.33%。基于 bin map，对三个环境下的三个性状进行定位，共鉴定出35个 QTL，每个 QTL 解释表型变异2.6％～15.68%。 在三个环境中都检测到多效 QTL（pQTL10），位于标记mk4012和 mk4037间，物理距离约14.6 Mb，表明其与三个性状都相关。 5 候选基因预测 根据玉米基因注释结果，qPH10包含45个蛋白编码的基因，但是只有7个具有功能注释。其中，2个为 MYB 转录因子家族的基因：GRMZM2G325907和GRMZM2G108892，可能调节株高性状，是qPH10的候选基因。 表1 玉米高密度遗传图谱结果统计 Chr.a No.markers b Physical distance (Mb) Genetic distance (cM) Avg.distance between markers (cM) <5cM Gap Max.gap (cM) 18.483723.084.93234.103837139.463354.064.15198.732733275.547443.033.36132.232674315.654473.022.36130.242744425.268453.068.07139.712784511.544353.083.02183.961643649.249363.042.34118.671593769.27534.066.24153.571853851.1112383.050.22120.751323911.547274.079.82136.941672015 1.119 6147 3.05 6.54517 .95023 814l a t o T a Chr.,indicates chromosome b No.markers,the number of markers on chromosome 图4 三个环境下的PH性状相关QTL在染色体上的分布 图3 玉米棒颜色性状相关QTL在染色体上的分布 图1 亲本间多态性SNP在全基因组及外显子区域的分布 图2 玉米 bin map (横轴表示染色体编号，纵轴表示样本数；红色表示与 亲本Qi319基因型相同，蓝色表示与亲本Ye478相同；黄色：杂合基因型)
诺禾致源高分文章集锦-人类篇
RNA 测序技术揭秘继发胶质母细胞瘤中新的重复出现的PTPRZ1-MET 融合基因 RNA-seq of 272 gliomas revealed a novel, recurrent PTPRZ1-MET fusion transcript in secondary glioblastomas 研究对象：人脑胶质瘤 期刊：Genome Research 影响因子：14.63 合作单位：北京大学 发表时间：2014年8月 摘 要 Studies of gene rearrangements and the consequent oncogenic fusion proteins have laid the foundation for targeted cancer therapy. To identify oncogenic fusions associated with glioma progression, we catalogued fusion transcripts by RNA-seq of 272 gliomas. Fusion transcripts were more frequently found in high-grade gliomas, in the classical subtype of gliomas, and in gliomas treated with radiation/temozolomide. Sixty-seven in-frame fusion transcripts were identified, including three recurrent fusion transcripts: FGFR3-TACC3, RNF213-SLC26A11, and PTPRZ1-MET (ZM ). Interestingly, the ZM fusion was found only in grade III astrocytomas (1/13; 7.7%) or secondary GBMs (sGBMs, 3/20; 15.0%). In an independent cohort of sGBMs, the ZM fusion was found in three of 20 (15%) specimens. Genomic analysis revealed that the fusion arose from translocation events involving introns 3 or 8 of PTPRZ and intron 1 of MET . ZM fusion transcripts were found in GBMs irrespective of isocitrate dehydrogenase 1 (IDH1) mutation status. sGBMs harboring ZM fusion showed higher expression of genes required for PIK3CA signaling and lowered expression of genes that suppressed RB1 or TP53 function. Expression of the ZM fusion was mutually exclusive with EGFR overexpression in sGBMs. Exogenous expression of the ZM fusion in the U87MG glioblastoma line enhanced cell migration and invasion. Clinically, patients afflicted with ZM fusion harboring glioblastomas survived poorly relative to those afflicted with non-ZM -harboring sGBMs (P < 0.001). Our study profiles the shifting RNA landscape of gliomas during progression and reveled ZM as a novel, recurrent fusion transcript in sGBMs. 关键词 胶质瘤；融合基因；PTPRZ1-MET 研究背景 脑胶质瘤由神经上皮细胞恶化演进形成，是成人中 枢神经系统中最常见的恶性肿瘤，致死率极高，恶性脑 胶质瘤5年存活率仅5%，具有浸润生长、易复发、预 后差等特点。根据WHO 形态学分类标准，脑胶质瘤分 为四个级别和多种类型，其恶性程度不同，常规放化疗 的预后差别较大，不同分型以及原发和复发胶质瘤发生 发展的分子机制也有所不同，因此亟需发展针对分子病 理的新的分型标准，指导临床治疗和个体化医疗，提高 病人的生存期和生存质量。
诺禾致源-ChIP-seq 流程
ChIP-seq 染色体免疫共沉淀（ChIP）是一种用于研究蛋白质与DNA的体内相互作用的经典实验技术。采用特异性抗体将目的蛋白进行免疫沉淀，由此可以把目的蛋白所结合的基因组DNA片段也富集下来。通过与高通量测序技术的结合，对ChIP后的DNA产物进行测序分析，从全基因组范围内寻找目的蛋白的DNA结合位点，以高效率的测序手段得到高通量的数据结果。 1.1C hIP免疫沉淀实验流程 目前主要有两种不同的ChIP实验方法，大致流程如下（均以细胞样品的处理过程为例）： Cross-liking Chromatin Immunoprecipitation (X-ChIP) 1.甲醛处理细胞，使DNA-protein的相互结合作用被交联固定。 2.裂解细胞，得到全细胞裂解液。 3.超声处理，将基因组DNA打断至100-500 bp。 4.抗体免疫沉淀：在细胞裂解液中加入一抗和beads，并进行孵育。 5.采用合适的实验条件进行洗脱，并解交联。 6.通过qPCR对ChIP结果进行验证。 7.准备好的ChIP后的DNA样品可以用于ChIP Sequencing建库。 Native Chromatin Immunoprecipitation (N-ChIP) 1.通过非变性的方式得到核裂解液。 2.微球菌核酸酶（Micrococcal nuclease）消化染色质，得到单核小体或核小体 寡聚体。 3.抗体免疫沉淀：在细胞裂解液中前后加入一抗和beads，并进行孵育。 4.DNA分离。 5.通过qPCR对ChIP结果进行验证。 6.准备好的ChIP后的DNA样品可以用于ChIP Sequencing建库。 下面步骤由我们公司做：
诺禾致源有参转录组结题报告
NHXXXXXX_species转录组生物信息分析结题报告建库测序流程 Total RNA样品检测 文库构建 库检 上机测序 生物信息分析流程 结果展示及说明 原始序列数据 测序数据质量评估 参考序列比对分析 可变剪切分析 新转录本预测 SNP和InDel分析 基因表达水平分析 RNA-seq整体质量评估 基因差异表达分析 差异基因GO富集分析 差异基因KEGG富集分析 差异基因蛋白互作网络分析 参考文献 附录 文件目录列表 软件列表 Methods英文版 备注
一、建库测序流程 从RNA样品到最终数据获得，样品检测、建库、测序每一个环节都会对数据质量和数量产生影响，而数据质量又会直接影响后续信息分析的结果。为了从源头上保证测序数据的准确性、可靠性，诺禾致源对样品检测、建库、测序每一个生产步骤都严格把控，从根本上确保了高质量数据的产出。流程图如下：
1 Total RNA样品检测 诺禾致源对RNA样品的检测主要包括4种方法： (1) 琼脂糖凝胶电泳分析RNA降解程度以及是否有污染 (2) Nanodrop检测RNA的纯度（OD260/280比值） (3) Qubit对RNA浓度进行精确定量 (4) Agilent 2100精确检测RNA的完整性 2 文库构建 样品检测合格后，用带有Oligo（dT）的磁珠富集真核生物mRNA（若为原核生物，则通过试剂盒去除rRNA来富集mRNA）。随后加入fragmentation buffer将mRNA打断成短片段，以mRNA为模板，用六碱基随机引物（random hexamers）合成一链cDNA，然后加入缓冲液、dNTPs和DNA polymerase I合成二链cDNA，随后利用AMPure XP beads纯化双链cDNA。纯化的双链cDNA再进行末端修复、加A尾并连接测序接头，然后用AMPure XP beads进行片段大小选择，最后进行PCR富集得到最终的cDNA文库。构建原理图如下： 3 库检 文库构建完成后，先使用Qubit2.0进行初步定量，稀释文库至1ng/ul，随后使用Agilent 2100对文库的insert size进行检测，insert size符合预期后，使用Q-PCR方法对文库的有效浓度进行准确定量（文库有效浓度＞2nM），以保证文库质量。 4 上机测序 库检合格后，把不同文库按照有效浓度及目标下机数据量的需求pooling后进行HiSeq/MiSeq测序。